实验5复原完成.docx

下载文档 降价啦

5
0
约7.71千字
约 13页
2017-07-08 发布于重庆
举报
版权申诉
保障服务

实验5复原完成.docx

1、本文档共13页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

实验5复原完成

（一）用各种方法去除不同强度的高斯噪I=imread(rice.png);subplot(4,5,1);imshow(I);title(原图);subplot(4,5,11);imshow(I);title(原图);A=imnoise(I,gaussian,0.01);subplot(4,5,2);imshow(A);title(加高斯0.01噪声污染);O=double(A);B=imnoise(I,gaussian,0.09);subplot(4,5,12);imshow(B);title(加高斯0.09噪声污染);L=double(B);A1=imfilter(O,fspecial(average,3));subplot(4,5,3);imshow(A1,[]);title( 3*3算术均值);B1=imfilter(L,fspecial(average,3));subplot(4,5,13);imshow(B1,[]);title( 3*3算术均值);A2=exp(imfilter(log(O),fspecial(average,3)));subplot(4,5,4);imshow(A2,[]);title( 3*3几何均值);B2=exp(imfilter(log(L),fspecial(average,3)));subplot(4,5,14);imshow(B2,[]);title( 3*3几何均值);Q=-1.5;A3=imfilter(O.^(Q+1),fspecial(average,3))./imfilter(O.^Q,fspecial(average,3));subplot(4,5,5);imshow(A3,[]);title( Q=-1.5逆谐波);B3=imfilter(L.^(Q+1),fspecial(average,3))./imfilter(L.^Q,fspecial(average,3));subplot(4,5,15);imshow(B3,[]);title( Q=-1.5逆谐波);Q=1.5;A4=imfilter(O.^(Q+1),fspecial(average,3))./imfilter(O.^Q,fspecial(average,3));subplot(4,5,6);imshow(A4,[]);title( Q=1.5逆谐波);B4=imfilter(L.^(Q+1),fspecial(average,3))./imfilter(L.^Q,fspecial(average,3));subplot(4,5,16);imshow(B4,[]);title(Q=1.5逆谐波);A5=medfilt2(O);subplot(4,5,7);imshow(A4,[]);title(中值滤波);B5=medfilt2(L);subplot(4,5,17);imshow(B4,[]);title(中值滤波);[r,c]=size(O);for i=2:r-1 for j=2:c-1 a=[I(i-1,j-1),I(i-1,j),I(i-1,j+1),I(i,j-1),I(i,j),I(i,j+1),I(i+1,j-1),I(i+1,j),I(i+1,j+1)]; b=[I(i-1,j-1),I(i-1,j),I(i-1,j+1),I(i,j-1),I(i,j),I(i,j+1),I(i+1,j-1),I(i+1,j),I(i+1,j+1)]; m=max(a);n=min(b); A6(i,j)=(m+n)/2; end;end;subplot(4,5,8);imshow(A6,[]);title(中点滤波);[r,c]=size(L);for i=2:r-1 for j=2:c-1 a=[I(i-1,j-1),I(i-1,j),I(i-1,j+1),I(i,j-1),I(i,j),I(i,j+1),I(i+1,j-1),I(i+1,j),I(i+1,j+1)]; b=[I(i-1,j-1),I(i-1,j),I(i-1,j+1),I(i,j-1),I(i,j),I(i,j+1),I(i+1,j-1),I(i+1,j),I(i+1,j+1)]; m=max(a);n=min(b); B6(i,j)=(m+n)/2; end;end;subplot(4,5,18);imshow(B6,[]);title( 中点滤波);[m,n]=size(O);A8=I;%d=2d=2;for i=2:m-1 for j=2:n-1 img = O(i-1