《数字图像的加噪方法与模糊方法.docx》由会员分享，可在线阅读，更多相关《数字图像的加噪方法与模糊方法.docx（11页珍藏版）》请在三一文库上搜索。
1、实验5图像复原专业：信计1401班 姓名：黄露露推荐精选学号：1408060102一、实验目的（1）掌握数字图像的加噪方法与模糊方法（2）掌握含噪图像的去噪方法（3）掌握数字图像的复原方法二、实验内容（1）掌握数字图像的加噪方法（2）掌握含噪图像的去噪方法（3）掌握数字图像的模糊方法（4）应用复原函数进行图像的复原三、核心代码数字图像的加噪方法 给图像添加噪声I=imread('rice.png');Noisy1=imnoise(I,'speckle',0.1);noise=0.1*randn(size(I);Noisy2=imadd(I,im2uint8(no
2、ise);figure,subplot(2,2,1),imshow(I),title('原图像');推荐精选subplot(2,2,2),imshow(Noisy1),title('speckle随机噪声图像');subplot(2,2,3),imshow(Noisy2),title('自定义随机噪声图像');%生成盐噪声矩阵，即有些随机像素的值为255，其余为0salt=255*(fix(0.01+rand(size(I);Noisy3=imadd(I,im2uint8(salt);subplot(2,2,4),imshow(Noisy3),t
3、itle('自定义随机噪声图像');含噪图像的去噪方法I=imread('rice.png');%噪声生成器noisy=imnoise(I,'gaussian',0.1);推荐精选%滤波器模板生成ave_h=fspecial('disk',3);%图像滤波函数denoisy1=imfilter(noisy,ave_h);denoisy2=medfilt2(noisy,3,3);figure,subplot(2,2,1),imshow(I),title('原图像');subplot(2,2,2),imshow(noi
4、sy),title('含高斯噪声图像');subplot(2,2,3),imshow(denoisy1),title('均值滤波图像');subplot(2,2,4),imshow(denoisy2),title('中值滤波图像');对rice.png创建运动模糊、高斯模糊、均值模糊的滤波器，并对图像进行模糊处理，对比不同方法的模糊效果推荐精选I=imread('rice.png');I=I(1:120,1:100); %I=rgb2gray(I);PSF=fspecial('motion',31,11);J=im
5、filter(I,PSF,'conv');figure,subplot(2,2,1),imshow(I);title('源图像')subplot(2,2,2),imshow(J),title('运动模糊图像')PSF=fspecial('disk',5);J=imfilter(I,PSF,'conv');subplot(2,2,3),imshow(J),title('均值模糊图像')PSF=fspecial('gaussian',31,3);J=imfilter(I,PSF,'
6、;conv');subplot(2,2,4),imshow(J),title('高斯模糊图像')推荐精选对一幅含噪的模糊图像进行魏纳滤波复原I=imread('rice.png');psf=fspecial('motion',31,11);blurred=imfilter(I,psf,'full');wnr=deconvwnr(blurred,psf);figure,subplot(2,2,1),imshow(I),title('原图像')subplot(2,2,2),imshow(blurred),tit
7、le('运动模糊图像')subplot(2,2,3),imshow(wnr),title('魏纳滤波复原图像')psf2=fspecial('motion',20,11);wnr2=deconvwnr(blurred,psf2);subplot(2,2,4),imshow(wnr2),title('使用较短的PSF复原图像')推荐精选对一幅图像分别加入噪声并模糊，用不同大小的PSF进行复原，对比复原I=imread('rice.png');psf=fspecial('motion',21,11);b
8、lurred=imfilter(I,psf,'full');noise=0.1*randn(size(blurred);%生成随机噪声矩阵bn=imadd(blurred,im2uint8(noise);%bn为含噪的模糊图像nsr=sum(noise(:).2)/sum(im2double(I(:).2);%计算信噪比wnr=deconvwnr(bn,psf);wnr1=deconvwnr(bn,psf,nsr);NP=abs(fftn(noise).2);%计算噪声图像的能量IP=abs(fftn(im2double(I).2);%计算噪声图像的能量,可以用含噪的模糊图像代
9、替原图像推荐精选NCORR=fftshift(real(ifftn(NP);%计算噪声图像的相关系数ICORR=fftshift(real(ifftn(IP);%计算原图像的相关系数wnr2=deconvwnr(bn,psf,NCORR,ICORR);figure,subplot(2,3,1),imshow(I),title('原图像')subplot(2,3,2),imshow(blurred),title('运动模糊图像')subplot(2,3,3),imshow(bn),title('含噪的运动模糊图像')subplot(2,3,4),imshow(wnr),title('不考虑噪声的复原图像')subplot(2,3,5),imshow(wnr1),title('已知信噪比的复原图像')subplot(2,3,6),imshow(wnr2),title('已知NCORR和ICORR复原图像')推荐精选 （注：可编辑下载，若有不当之处，请指正，谢谢!） 推荐精选