北航数图实验报告三 图像复原实验.pdfVIP

北航数图实验报告三图像复原实验--第1页

京航空航天大学

数字图像处理实验报告

院

专业方向

班级

学号

学生姓名

指导教师

北航数图实验报告三图像复原实验--第1页

北航数图实验报告三图像复原实验--第2页

图像复原实验

1．实验目的

比较不同参数选择对复原

结果的影响。

2．实验内容

（1）利用反向滤波方法进行图像复原；

（2）利用维纳滤波方法进行图像复原。

3.实验要求

（1）输入图像采用实验1所获取的图像，对输入图像采用运动降质模

型，如下式所示

错误!不能通过编辑域代码创建对象。

与降值图像相关的参数是：错误!不能通过编辑域代码创建对象。；

（2）对每一种方法通过计算复原出来的图像的峰值信噪比，进行最优参数

的选择，包括反向滤波方法中进行复原的区域半径错误!不能通过编辑域代

码创建对象。、维纳方法中的噪声对信号的频谱密度比值K；

（3）将降质图像和利用最优参数恢复后的图像同时显示出来，以便比较。

4.程序代码

1.运动降质的源代码

saber.bmp);

I=rgb2gray(I);

F=fft2(I);

F=fftshift(F);

form=-((size(F,1)+1)/2):(size(F,1)+1)/2

forn=-size(F,2)/2:size(F,2)/2

H(m+((size(F,1)+1)/2)+1,n+size(F,2)/2+1)=5*sin(pi*(m+n))*exp(-sqrt(-1

)*pi*(m+n))/(pi*(m+n));

end

end

form=1:size(F,1)

forn=1:size(F,2)

北航数图实验报告三图像复原实验--第2页

北航数图实验报告三图像复原实验--第3页

if(isnan(real(H(m,n))))

G(m,n)=F(m,n);

H(m,n)=1;

else

G(m,n)=F(m,n)*H(m,n);

end

end

end

g=ifft2(G);

t1=abs(g);

mm=max(max(t1));mn=min(min(t1));

t1=255+255/(mm-mn)*(t1-mm);

t=uint8(t1);

axes(handles.axes2);

imshow(t)