资源描述:
《数字图像的加噪方法与模糊方法.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
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
2、=imadd(I,im2uint8(noise));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(
3、salt));subplot(2,2,4),imshow(Noisy3),title('自定义随机噪声图像');含噪图像的去噪方法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),t
4、itle('原图像');subplot(2,2,2),imshow(noisy),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);
5、PSF=fspecial('motion',31,11);J=imfilter(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=i
6、mfilter(I,PSF,'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(blurr
7、ed),title('运动模糊图像')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);blurr
8、ed=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);wnr