Matlab数字图像处理尝试之：[7]邻域中值

1、计算去噪前后的峰值信噪比，做出主观和客观评价。1），邻域平均化法：clear all;>>I=imread('C:\Users\zjy\Desktop\test1.tif');subplot(2,3,1),imshow(I);>>xlabel('a)原始图像');>>J=imnoise(I,'gaussian',0,0.025);subplot(2,3,2),imshow(J);>>xlabel('b)添加白噪声');>>K=imnoise(I,'salt & pepper',0.025);subplot(2,3,3),imshow(K);xlabel('b)添加椒盐噪声');>> subplot(2,3,4)>> subplot(2,3,4),imshow(I);xlabel('a)原始图像');>> h=[1 1 1;1 0 1;1 1 1];>> h=h/8;>> J1=conv2(J,h);>> subplot(2,3,5),imshow(J1,[]);>> xlabel('邻域平均除白噪声');>>K1=conv2(K,h);>> subplot(2,3,6),imshow(K1,[]);>> xlabel('邻域平均除椒盐噪声');以上程序段实现了用邻域平均化去除高斯白噪声和椒盐噪声，效果如下图所示。

2、从上图可以看出，邻域平均法对高斯白噪声的消除效果比较好，对椒盐噪声的消除效果不理想。

3、2），中值滤波法：clear all;I=imread('C:\Users\zjy\Desktop\test1.tif');subplot(2,3,1),imshow(I);xlabel('a)原始图像');J=imnoise(I,'gaussian',0,0.025);subplot(2,3,2),imshow(J);xlabel('b)添加白噪声');K=imnoise(I,'salt & pepper',0.025);subplot(2,3,3),imshow(K);xlabel('b)添加椒盐噪声');subplot(2,3,4)subplot(2,3,4),imshow(I);xlabel('a)原始图像');J1=medfilt2(J);subplot(2,3,5),imshow(J1,[]);xlabel('中值滤波除白噪声');J2=medfilt2(K);subplot(2,3,6),imshow(J2,[]);xlabel('中值滤波除椒盐噪声');以上程序段实现了用邻域平均化去除高斯白噪声和椒盐噪声，效果如下图所示。

4、 从上图可以看出，中值滤波法对高斯白噪声如果用5*5窗口滤波效果好于3*3窗口滤波窗口，但图像模糊会加重。中值滤波法对椒盐噪声的消除效果很理想。