图像复原处理技术

1、实验五图像复原处理技术实验目的1 了解图像降质退化的原因，并建立降质模型。2 理解反向滤波图像复原的原理3 理解维纳滤波图像复原的原理实验原理图像复原处理一定是建立在图像退化的数学 模型基础上的，这个退化数学模型应该能够反映图像退化的原因。图像降质过程的模型如图51所示，其表达式为g(x,y)= h (x,y)*f (x,y) +n( xy)(5.1 )图51图像降质模型1、滤波图像复原逆滤波法是最简单的图像恢复方法。对5.1式两边作二维傅立叶变换，得到G (u, v) =H (u,v) F (u,v) + N (u,v)H (u,v)为成像系统的转移函数。估算得到的恢复图像的傅立叶变换F (

2、u,v)为F? u,v 二匕F u,vH u,vH u,v(5.2)若知道转移函数H u,v ,5.2式经反变换即可得到恢复图像，其退化和恢复的全过 程用图5-2表示。图5 2频域图像降质及恢复过程逆滤波恢复法会出现病态性，若H u,v，而噪声N(u,v)旳，则N u,v；'H u,v比F(x,y)大很多，使恢复出来 ? x,y与f x,y相差很大，甚至面目全非。一种改进的 方法是在H (u, v) =0的频谱点及其附近，人为仔细设置H J u,v的值，使得在这些频谱点附近，N u,v H u,v不会对F? u,v产生太大影响。二种方法是考虑到降质系统 的转移函数H u,v的带宽比噪声

3、要窄的多，其频率特性也具有低通性质，因此可令逆滤 波的转移函数u,v为#H(u,v)(u2+v2兰 D。hi u,v,0( u2 + v2 $D0(2)维纳滤波复原逆滤波简单，但可能带来噪声的放大，而维纳滤波对逆滤波的噪声放大有抑制作用。维纳滤波是寻找一个滤波器，使得复原后图像? x, y与原始图像f x, y的方差最小，_2min=E | f x, y - ？ x,y f如果图像f x, y和噪声n x, y不相关，且h x, y有零均值，则可导出维纳滤波器的传递函数为Hw u,v =1H u,v *H u,v 2H u,v 2Pn u,VPf u,v式中Pn(u,v )和Pf(u,v )分

4、别为噪声和原始图像的功率谱。实际上Pn(u,v )和Pf(u,v)n往往是未知的，这时常用常数K来近似鹽【实验】产生一模糊图像，采用维纳滤波图像复原的方法对图像进行处理。clear; %清除变量 d=15 %设定长度h=zeros(2*d+1,2*d+1);h(d+1,1:2*d+1)=1/(2*d); %设置函数 hf=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+

5、1)=h; % 扩增 hF=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.A2)./(H.A2+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,

6、1:n),min(f_est(:) max(f_est(:);%运行结果:原图模糊后作业:恢复图象1改变维纳滤波的常数 K,比较不同K值的复原图像的效果。K=0.5K=1K=10可见，K 越大，对消除模糊的效果越差，但是可以看清轮廓。K小的时候，虽然没有模糊的感觉，但是图像内容很不清晰，难以辨别。2 对图像文件 cameraman.tif 进行运动模糊处理，然后分别采用逆滤波和维纳滤波的方法 进行复原处理。clear; %清除变量c=0.1;T=1;f=imread('cameraman.tif');%f=rgb2gray(f);m,n=size(f);H=zeros(m,n)

7、;for j=1:mH(j,:)=T/(pi*j*c)*sin(pi*j*c)*exp(-i*(pi*j*c);endF=fft2(double(f);NF=F.*H;% 模糊图像的频谱newimg=real(ifft2(NF);subplot(2,2,1);imshow(f);title 惊图');% 显示原始图像subplot(2,2,2);imshow(ui nt8( newimg);显示动态模糊后的图像title(动态模糊后的图');% 维纳滤波K=0; %设定K值F_est=(H.A2)./(H.A2+K).*NF./H;f_est=real(ifft2(F_est); %恢复后的图像subplot(2,2,3);imshow(ui nt8(f_est);%显示恢复后的图像 title('维纳滤波图');%逆滤波FN=NF./H;su