逆滤波与维纳滤波

数字图像处理学院：电信学院专业：信息工程班级：2008级3班姓名：朱志俊学号：200830101124实验作业5用式（5.6-3)对图像进行模糊处理，然后加白高斯噪声，用逆滤波和维纳滤波恢复图像，对结果进行分析。（1）程序：clcI=imread('E:\大三下\数字图像处理\MATLAB图像处理\水果2.jpg');figureimshow(I);title('原图像');[m,n]=size(I);F=fftshift(fft2(I));k=0.0025;H=[];foru=1:mforv=1:nq=((u-m/2)^2+(v-n/2)^2)^(5/6);H(u,v)=exp((-k)*q);endendG=F.*H;I0=abs(ifft2(fftshift(G)));figureimshow(uint8(I0));title('退化的图像');I1=imnoise(uint8(I0),'gaussian',0,0.01)figureimshow(uint8(I1));title('退化并且添加高斯噪声的图像');F0=fftshift(fft2(I1));5/622(,)kuvHuveF1=F0./H;I2=ifft2(fftshift(F1));figureimshow(uint8(I2));title('全逆滤波复原图');K=0.1;H=[];H0=[];H1=[];foru=1:mforv=1:nq=((u-m/2)^2+(v-n/2)^2)^(5/6);H(u,v)=exp((-k)*q);H0(u,v)=(abs(H(u,v)))^2;H1(u,v)=H0(u,v)/(H(u,v)*(H0(u,v)+K));endendF2=H1.*F0;I3=ifft2(fftshift(F2));figureimshow(uint8(abs(I3)));title('维纳滤波复原图');（2）结果：（3）结果分析：当退化函数：中的K=0.0025时，维纳滤波效果比逆滤波效果要好。逐渐减小K值亦可改善维纳滤波5/622(,)kuvHuve效果，如当K=0.00025时，逆滤波结果如下：维纳滤波结果如下：逆滤波效果虽有改善，但远不及威娜滤波效果。在实验中遇到这样的问题：进行傅里叶反变换后对于显示图像得取变换结果的绝对值，否则会出现下面的问题