实验三图像滤波

. . .

. .

实验三 图像滤波实验

一、实验目的

1.熟练掌握空域滤波中常用滤波器的原理；

2.理解邻域平均、中值滤波降噪法等的图像增强原理；

3.加深对图像增强的认识，巩固所学图像增强理论知识和相关算法； 4.利用MATLAB程进行图像滤波，观察增强效果。 二、实验原理

以图像平滑为例：图像平滑方法是一种实用的图像处理技术，能减弱或消除图像中的高频率分量，但不影响低频率分量。因为高频率分量主要对应图像中的区域边缘等灰度值较大变化较快的部分，平滑滤波将这些分量滤去可减少局部灰度起伏，使图像变得比较平滑。实际应用中，平滑滤波还可用于消除突发噪声，或者在提取较大目标前去除过小的细节或将目标内的小间断连接起来。 1. 邻域平均法

用窗口像素的平均值取代中心像素原来的灰度值。邻域的选取通常有4-邻域和8-邻域。

g(i,j)?1?f(x,y)L(x,y)?A

2. 中值滤波

与加权平均方式的平滑滤波不同，中值滤波是抑制噪声的非线性处理方法。对于给定的n个数值{a1,a2,…,an}，将它们按大小有序排列。当n为奇数时，位于中间位置的那个数值称为这n各数值的中值。当n为偶数时，位于中间位置的两个数值的平均值称为这n个数值的中值。中值滤波，就是图像中滤波后的某像素输出值等于该像素邻域中各像素灰度的中值。 三、实验内容与步骤：按下面要求编写程序并运行

对peppers.bmp图像添加零均值“高斯”噪声，用3×3或9×9滑动平均模板进行降噪处理，观察平均降噪图像的效果。 1.对peppers.bmp图像添加零均值“高斯”噪声 程序代码如下：

I=imread('d:\\peppers.bmp');

.

word . .

. . .

. .

J=imnoise(I,'gaussian',0,0.02);

subplot(1,2,1);imshow(I);title('原图');

subplot(1,2,2);imshow(J);title('加噪处理后的图');

运行结果如下：

2.用3×3滑动平均模板进行降噪处理，观察平均降噪图像的效果 程序代码如下：

I=imread('d:\\peppers.bmp'); J=imnoise(I,'gaussian',0,0.02); [H,W]=size(J); M=double(J); M1=M; for j=2:H-1 for k=2:W-1

M1(j,k)=[M(j-1,k-1)+M(j-1,k)+M(j-1,k+1)+M(j,k-1)+M(j,k)+M(j,k+1)+M(j+1,k-1)+M(j+1,k)+M(j+1,k+1)]/9; end; end;

subplot(1,3,1);imshow(I);title('原图');

.

word . .

. . .

. .

subplot(1,3,2);imshow(J);title('加噪处理');

subplot(1,3,3);imshow(uint8(M1),[]);title('3*3邻域平均降噪处理'); 运行结果如下：

3.对peppers.bmp图像添加“椒盐”噪声，分别采用5×5的矩形、圆形及十字形窗口进行中值滤波处理，观察处理结果。（滤波过程采用迭代操作，对输入图像重复进行同样的中值滤波，直到输出不再有变化为止。） 3.1.对peppers.bmp图像添加“椒盐”噪声 程序代码如下：

I=imread('d:\\peppers.bmp');

J = imnoise(I,'salt & pepper',0.2); subplot(1,2,1);imshow(I);title('原图'); subplot(1,2,2);imshow(J);title('椒盐加噪处理'); 运行结果如下：

.

word . .

. . .

. .

3.2采用5×5的矩形中值滤波处理 程序代码如下：

I= imread('d:\\peppers.bmp'); J=imnoise(I,'salt & pepper',0.2); P=double(J); [H,W]=size(J); flag=1;

while( flag==1) %若输出跟前一个图有变化则进行中值滤波处理 flag=0; for i=3:(H-2) for j=3:(W-2)

A=P(i-2:i+2,j-2:j+2); %5*5矩阵窗口模板 mid=median(A(:)); %求矩阵模板中值 if(P(i,j)~=mid) P(i,j)=mid; flag=1; end

.

word . .

. . .

. .

end end

if flag==0 %若输出不再变化停止滤波操作 break; end; end;

subplot(1,3,1),imshow(I);title('原图');

subplot(1,3,2),imshow(J);title('椒盐加噪处理 ');

subplot(1,3,3),imshow(uint8(P));title('5*5矩形中值滤波'); 运行结果如下：

3.3采用圆形窗口中值滤波处理 程序代码如下：

I= imread('d:\\peppers.bmp'); J=imnoise(I,'salt & pepper',0.2); P=double(J); [H,W]=size(J); flag=1; while( flag==1) flag=0;

.

word . .