实验报告七

南京信息工程大学数字图像处理实验报告实验课程数字图像处理实验名称图像频域增强实验日期 2015-6-05 指导老师徐旦华专业信息与计算科学年级 2012 姓名杨鑫学号 20121314059 得分

一、实验目的

1.掌握噪声模拟和图像滤波函数的使用方法

2.了解图像复原的Matlab实现方法

二、实验内容

1.模糊与噪声

A)对图像‘ysp.jpg’分别采用运动PSF和均值滤波PSF进行模糊。

B)在上述模糊图像上再添加噪声

I=imread('ysp.jpg');

subplot(2,3,1);imshow(I);

title('原图像');

H=fspecial('motion',45,50);

motionblur=imfilter(I,H);

subplot(2,3,2);

imshow(motionblur);

title('运动PSF');

H=fspecial('disk',10);

averageblur=imfilter(I,H);

subplot(2,3,3);imshow(averageblur)

title('均值滤波PSF')

J=imnoise(motionblur,'salt & pepper');

subplot(2,3,5);imshow(J)

title('对运动PSF加噪声')

W=imnoise(averageblur,'salt & pepper');

subplot(2,3,6);imshow(W)

title('对均值PSF加噪声 ')

figure,freqz2(psf);

title('PSF三维图');

原图

像运动

PSF 均值滤波

PSF

对运动PSF 加噪

声对均值PSF 加噪声

F x PSF 三维图

F y M a g n i t u d e

2.维纳滤波复原

len=30; theta=10; PSF=fspecial('motion',len,theta);

blurred=imfilter(I,PSF,'circular','conv');

subplot(2,3,1),imshow(blurred)

title('真实PSF')

wnr2=deconvwnr(blurred,fspecial('motion',2*len,theta));

subplot(2,3,2),imshow(wnr2)

title('2倍len 的非真实PSF')

wnr3=deconvwnr(blurred,fspecial('motion',len,2*theta));

subplot(2,3,3),imshow(wnr3)

title('2倍theta的非真实PSF')

noise=0.1*randn(size(I));

BlurredNoisy=imadd(blurred,im2uint8(noise));

subplot(2,3,4),imshow(BlurredNoisy)

title('添加模糊噪声图像');

NP=abs(fftn(noise)).^2;

NCORR=fftshift(real(ifftn(NP)));

IP=abs(fftn(im2double(I))).^2;

ICORR=fftshift(real(ifftn(IP)));

wnr6=deconvwnr(BlurredNoisy,psf,NCORR,ICORR);

subplot(2,3,5),imshow(wnr6)

title('最终复原图');

真实PSF2倍len的非真实PSF2倍theta的非真实PSF

添加模糊噪声图像最终复原图

三、思考题

1.常见的图像退化模型有哪些种类？

答:运动模糊退化模型，离焦模糊退化模型

2.不同的PSF对复原效果有什么影响？

答：不同的PSF参数值对复原效果影响较大，对不同的PSF参数值产生的复原效果不同，真实的PSF通常是未知的，需要根据一定的只是进行估计，再将估计值作为参数进行图像复原。