数字图像处理实验报告3

实验三数字图像的空间域滤波

院（系）：电气学院姓名：刘丽霞学号：200902050068

实验时间：2012/11/20 实验地点：A07-602 指导老师：李静

一、实验目的

1）掌握图像滤波的基本定义及目的。

2）理解空间域滤波的基本原理及方法。

3）掌握进行图像的空域滤波的方法。

二、实验环境

1）在MATLAB的环境下进行实验；

2）复制若干图形文件至MATLAB目录下work文件夹中。

三、实验方法

1.平滑空间滤波：

1)读出eight.tif这幅图像，给这幅图像分别加入椒盐噪声和高

斯噪声后并与前一张图显示在同一图像窗口中。

I1=imnoise(I,'salt & pepper',0.02); %加入椒盐噪声

I2=imnoise(I,'gaussian',0,0.005); %加入高斯噪声

2)对加入噪声图像选用不同的平滑（低通）模板做运算，对比不

同模板所形成的效果，要求在同一窗口中显示。

Hs = fspecial('sobel');

Sobel = imfilter(I,Hs,'replicate'); %选用不同的平滑模板

3)使用函数imfilter时，分别采用不同的填充方法（或边界选

项，如零填充、’replicate’、’symmetric’、’circular’）进行低通滤波，显示处理后的图像。

4)运用for循环，将加有椒盐噪声的图像进行10次，20次均值

滤波，查看其特点,显示均值处理后的图像。

II1=rgb2gray(I1);

a=1/9.*[1 1 1

1 1 1

1 1 1];

b=conv2(a,II1); %均值处理

5)对加入椒盐噪声的图像分别采用均值滤波法，和中值滤波法对

有噪声的图像做处理，要求在同一窗口中显示结果。

2.锐化空间滤波

读出moon.tif这幅图像，采用3×3的拉普拉斯算子w = [ 1, 1, 1; 1 – 8 1; 1, 1,

1]对其进行滤波。

w =[1,1,1;

1,-8,1;

1,1,1];

K=conv2(T,w,'same');%3×3的拉普拉斯算子滤波

四、实验结果

不同尺寸梯度算子的锐化滤波1 不同尺寸梯度算子的锐化滤波2

不同尺寸梯度算子的锐化滤波3 不同尺寸梯度算子的锐化滤波4

5

不同尺寸梯度算子的锐化滤波5

I=imread('moon.tif');

T=double(I);

subplot(1,2,1),imshow(T,[]);title('设计原图');

domain=[8 8 0 8 8;

8 8 0 8 8;

0 0 0 0 0;

8 8 0 8 8;

8 8 0 8 8];

K1= ordfilt2(T,5,domain);

subplot(1,2,2),imshow(K1,[]);title('自行设计图');

'设计原图自行设计图

五、实验结果分析和结论

六、参考文献

[1] R.C. Gonzalez & R.E Woods. Digital Image Processing(2E). Prentice Hall PTR，2002 中译本：阮秋琦等译. 数字图像处理(第二版). 北京：电子工业出版社，2005.3

[2] 冈萨雷斯. 数字图像处理(MA TLAB版). 北京：电子工业出版社，2005.9

[3] K.R. Castleman. 数字图像处理. 朱志刚等译. 北京：电子工业出版社，1998.9

[4] 章毓晋. 图象工程(上册)—图象处理和分析. 北京：清华大学出版社，2003.7

[5] 何斌. Visual C++数字图像处理(第二版). 北京：人民邮电出版社，2002.1