均值和中值滤波_图文

1、1、均值滤波器邻域平均法是一种局部空间域处理的算法。设一幅图像f(x,y为N ×N 的阵列,处理后 的图像为g(x,y,它的每个像素的灰度级由包含(x,y领域的几个像素的灰度级的平均值所 决定,即用下式得到处理后的图像:=S j i j i f M y x g ,(,(1,(式中x,y=0,l,2,N-1;S 是以(x,y为中心的邻域的集合,M 是S 内坐标点的总数。图像邻域平均法的处理效果与所用的邻域半径有关。半径愈大,则图像的模糊程度也度大。外,图像邻域平均算法简单,计算速度快,但它的主要缺点是在降低噪声的同时使图像产生模糊,特别在边沿和细节处,邻域越大,模糊越厉害。一个均值滤波

2、的例子:I=imread('cameraman.tif'%读入图像J=imnoise(I,'salt & pepper',0.02;%给图像添加椒盐噪声K=imnoise(I,'gaussian',0,0.005;%给图像添加均值为0,方差为0.005的高斯噪声subplot(231,imshow(Ititle('原图像'subplot(232,imshow(Jtitle('添加椒盐噪声图像'subplot(233,imshow(Ktitle('添加高斯噪声图像'subplot(234,im

3、show(Ititle('原图像'K1=filter2(fspecial('average',3,J/255;%使用3×3模板均值滤波subplot(235,imshow(K1title('3*3椒盐噪声均值滤波'K2=filter2(fspecial('average',3,K/255;%使用3×3模板均值滤波subplot(236,imshow(K2title('3*3高斯噪声均值滤波' 2、中值滤波器它是一种基于排序统计理论的可有效抑制噪声的非线性平滑滤波器,其基本原理是把数字图像或数字序

4、列中一点的值用该点的一个领域中各点值的中值代换。其主要功能是让周围象素灰度值的差比较大的像素改成与周围的像素值接近的值,从而可以消除孤立的噪声点,所以中值滤波对于滤除图像的椒盐噪声非常有效。中值滤波去除噪声的效果除了与噪声的类型有关外,还与邻域的空问范围和中值计算中涉及的象素数有关。一般来说,小于滤波器面积一半的亮或暗的物体基本上会被滤除,而较大的物体几乎原封不动地保存下来,因此,中值滤波器的空间尺寸必须根据现有的问题来进行调整。较简单的模板是N×N的方形(N通常为奇数。中值滤波器可以做到既去除噪声又能保护图像的边缘,从而获得较满意的复原效果,而且,在实际运算过程中不需要图像的统计特

5、性,这也带来不少方便,但对一些细节多,特别是点、线、尖项细节较多的图像不宜采用中值滤波的方法。一个中值滤波的例子:I=imread('cameraman.tif'%读入图像J=imnoise(I,'salt & pepper',0.02;%给图像添加椒盐噪声K=imnoise(I,'gaussian',0,0.005;%给图像添加均值为0,方差为0.005的高斯噪声subplot(231,imshow(Ititle('原图像'subplot(232,imshow(Jtitle('添加椒盐噪声图像'subpl

6、ot(233,imshow(Ktitle('添加高斯噪声图像'subplot(234,imshow(Ititle('原图像'K1=medfilt2(J,3,3;%使用3×3模板中值滤波subplot(235,imshow(K1title('3*3椒盐噪声中值滤波'K2=medfilt2(K,3,3;%使用3×3模板中值滤波subplot(236,imshow(K2title('3*3高斯噪声中值滤波' 3、均值滤波与中值滤波的对比I=imread('cameraman.tif'%读入图像J=im

7、noise(I,'salt & pepper',0.02;%给图像添加椒盐噪声K=imnoise(I,'gaussian',0,0.005;%给图像添加均值为0,方差为0.005的高斯噪声subplot(331,imshow(Ititle('原图像'subplot(332,imshow(Jtitle('添加椒盐噪声图像'subplot(333,imshow(Ktitle('添加高斯噪声图像'subplot(334,imshow(Ititle('原图像'K1=filter2(fspecial('average',3,J/255;%使用3×3模板均值滤波subplot(335,imshow(K1title('3*3椒盐噪声均值滤波'K2=filter2(fspecial('average',3,K/255;%使用3×3模板均值滤波subplot(336,imshow(K2title('3*3高斯噪声均值滤波'subplot(337,imshow(Ititle('原图像'K3=medfilt2(J,3,3;%使用3×3模板中值滤波subplo