数字图像的加噪方法与模糊方法

实验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(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(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),title(原图像);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);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=imfilter(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(blurred),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);blurred=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);%计算噪声图像的能量,可以用含噪的模糊图像代替原图像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