MATLAB图像复原处理技术

MATLAB图像复原处理技术

实验目的

1. 了解图像降质退化的原因，并建立降质模型。
2. 理解反向滤波图像复原的原理
3. 理解维纳滤波图像复原的原理实验原理图像复原处理一定是建立在图像退化的数学模型基础上的，这个退化数学模型应该能够

反映图像退化的原因。图像降质过程的模型如图5－1所示，其表达式为

g(x,y)=h (x,y)*f (x,y) +n(xy) （5.1）

图5－1图像降质模型

1. 滤波图像复原

逆滤波法是最简单的图像恢复方法。对5.1式两边作二维傅立叶变换，得到

G (u, v) =H (u,v) F (u,v) + N (u,v)

H (u,v) 为成像系统的转移函数。估算得到的恢复图像的傅立叶变换Fˆ (u,v) 为

(5.2)

若知道转移函数H

,5.2式经反变换即可得到恢复图像，其退化和恢复的全过程用图5-2表示。

图5－2频域图像降质及恢复过程

（2）维纳滤波复原

【实验】产生一模糊图像，采用维纳滤波图像复原的方法对图像进行处理。

clear; %清除变量

d=15 %设定长度

h=zeros(2*d+1,2*d+1);

h(d+1,1:2*d+1)=1/(2*d); %设置函数h

f=imread('lena.bmp'); %读取图像

[m,n]=size(f); %求出图像大小

fe=zeros(m+2*d,n+2*d); %扩增f

fe(1:m,1:n)=f;

he=zeros(m+2*d,n+2*d);

he(1:2*d+1,1:2*d+1)=h; %扩增h

F=fft2(fe);

H=fft2(he);

ns=5*rand(m+2*d,n+2*d); %产生噪声

g=ifft2(F.*H)+ns; %产生模糊且加载噪声的图像

G=fft2(g);

K=0; %设定K 值

F_est=((H.^2)./(H.^2+K)).*G./H; % 维纳滤波

f_est=real(ifft2(F_est)); %恢复后的图像

imshow(f); %显示原始图像

figure;

imshow(g(d+1:m+d,d+1:n+d),[min(g(:)) max(g(:))]); % 显示模糊后加噪声的图像

figure;

imshow(f_est(1:m,1:n),[min(f_est(:)) max(f_est(:))]);% 显示恢复后的图像

运行结果：

原图

模糊后