图像处理翻译

       é’ˆå¯¹ç»™å‡ºçš„图像（moon.tif）或者自行选择的灰度图像：

        1）：给图像分别添加高斯噪声和椒盐噪声。

        2）：对加噪图像的中心区域（*）进行空间滤波，尽最大可能消除噪声。

        3）：对加噪图像的中心区域（*）进行频域滤波，尽最大可能消除噪声。

       æŠ€æœ¯æè¿°ï¼š

       å¯¹å›¾åƒè¿›è¡ŒåŠ é«˜æ–¯å™ªå£°å’Œæ¤’盐噪声处理；对包含高斯噪声和椒盐噪声的图片进行处理，使处理后的图像比原图像清晰。

       æ‰€éœ€åº”用到的技术，包括：

       a>对图片加噪声

       b>选取中心区域

       c>对选取的区域进行降噪处理

       d>重新生成图像。

       e>构造高斯低通滤波器时用到了高斯公式：exp(-(u^2+v^2)/(2*(D0^2)))

       ç»“果讨论：

       ä»¥ä¸‹æ˜¯å¯¹ä¸åŒçš„滤波器针对不同噪点处理的测试结果。

       å‚考下面的试验结果，进行讨论：

       A()是使用fspecial('gaussian’)平滑空域滤波处理效果，不过效果不是最好，由于最大程度降噪，导致图像模糊；

       A()是频域滤波处理后的结果，因为使用了高斯低通滤波器，所以会有条黑线，处理一般；

       A()是使用medfilt2（）空域中值滤波器效果，降噪效果很不错，图像也很清晰；

       A() 是频域滤波处理后的结果，同A()，因为使用了高斯低通滤波器，所以会有条黑线，效果一般。

       è¯•éªŒç»“æžœ:

       é«˜æ–¯åŠ å™ªå’Œæ¤’盐加噪处理图分别如下:

        如图:

        图(A):原图

        图(A):高斯加噪

        图(A):对图(A)进行中心*空域滤波

        图(A):对图(A)进行中心*频域滤波

       ï¼ˆA0） （A）

       ï¼ˆA） (A)

        如图:

        图(A):原图

        图(A):椒盐加噪

        图(A):对图(A)进行中心*空域滤波

        图(A):对图(A)进行中心*频域滤波

       (A) (A)

       (A) (A)

       é™„录：

       æºä»£ç 1 ：对高斯噪声的处理

       f=imread('moon.tif');

       J=imnoise(f,源码最牛K线图源码'gaussian',0.,0.);%添加高斯噪声

       %空域滤波

       r=[ ];

       c=[ ];

       BW=roipoly(J,c,r);

       h=fspecial('gaussian',[5 5]);

       A=roifilt2(J,h,BW);

       figure,imshow(A);

       %频域滤波

       f1=imcrop(fn,[ ]);

       %截取*大小的窗口图片

       f2=[ ];

       %建立一个新的图像

       f2=uint8(f2);

       f2=padarray(f2,[ ],);

       %将新建图像拓展到*的黑色图片

       f2=padarray(f2,[ ],0);

       %在新建图片周围添加白色使之大小为moon图片的大小

       fn=fn-f2;

       %得到中心*区域内为黑色的moon图片

       PQ=paddedsize(size(f1));

       [u,v]=dftuv(PQ(1),PQ(2));

       D0=0.2*PQ(2);

       hh=exp(-(u.^2+v.^2)/(2*(D0^2)));

       %构造高斯低通滤波器

       h1=dftfilt(f1,hh);

       A=padarray(h1,[ ],0);

       %将h1拓展到moon图片大小

       A=uint8(A)+fn;

       %得到中心*区域处理后的moon图片

       figure,imshow(A);

       æºä»£ç 2 ：对椒盐噪声的处理

       f=imread('moon.tif');

       fn=imnoise(f,'salt & pepper',0.);

       %添加椒盐噪声

       %空域滤波

       f1=imcrop(fn,[ ]);

       %截取*大小的窗口图片

       f2=[ ];

       %建立一个新的图像

       f2=uint8(f2);

       f2=padarray(f2,[ ],);

       %将新建图像拓展到*的黑色图片

       f2=padarray(f2,[ ],0);

       %在新建图片周围添加白色使之大小为moon图片的大小

       fn=fn-f2;

       %得到中心*区域内为黑色的moon图片

       h=medfilt2(f1,'symmetric');

       %对f1进行中值处理

       A=padarray(h,[ ],0);

       %将h拓展到moon图片大小

       A=A+fn;

       %得到中心*区域处理后的moon图片

       figure,imshow(A);

       %频域滤波

       f1=imcrop(fn,[ ]);

       %截取*大小的窗口图片

       f2=[ ];

       %建立一个新的图像

       f2=uint8(f2);

       f2=padarray(f2,[ ],);

       %将新建图像拓展到*的黑色图片

       f2=padarray(f2,[ ],0);

       %在新建图片周围添加白色使之大小为moon图片的大小

       fn=fn-f2;

       %得到中心*区域内为黑色的moon图片

       PQ=paddedsize(size(f1));

       [u,v]=dftuv(PQ(1),PQ(2));

       D0=0.1*PQ(2);

       hh=exp(-(u.^2+v.^2)/(2*(D0^2)));%构造高斯低通滤波器

       h1=dftfilt(f1,hh);

       A=padarray(h1,[ ],0);

       %将h1拓展到moon图片大小

       A=uint8(A)+fn;

       %得到中心*区域处理后的moon图片

       figure,imshow(A);