(完整版)对图像进行模糊处理,然后加白高斯噪声,用逆滤波和维纳滤波恢复图像.pdf

作业六 5/6 2 2 k u v 用式（ 5.6-3) H (u,v) e 对图像进行模糊处理，然后加白高斯噪声，用逆滤波和维纳 滤波恢复图像，对结果进行分析。 程序代码如下： clear; A=imread( picture4.jpg ); figure; imshow(A); title( 原图像 ); I = rgb2gray(A); figure; imshow(I); title( 二维灰度图像 ); %将彩色图像转换为二维灰度图像 [m,n]=size(I); F=fftshift(fft2(I)); k=0.0025; H=[]; for u=1:m for v=1:n q=((u-m/2)^2+(v-n/2)^2)^(5/6); H(u,v)=exp((-k)*q); end end G=F.*H; I0=abs(ifft2(fftshift(G))); figure; imshow(uint8(I0)); title( 退化的图像 ); I1=0.1*randn(size(I0)); I1=I0+I1; %生成高斯白噪声图像 figure; imshow(uint8(I1)); title( 退化并且添加高斯白噪声图像 ); F0=fftshift(fft2(I1)); F1=F0./H; I2=ifft2(fftshift(F1)); figure; imshow(uint8(I2)); title( 全逆滤波复原图 ); K=0.1; H=[]; H0=[]; H1=[]; for u=1:m for v=1:n q=((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)); end end F2=H1.*F0; I3=ifft2(fftshift(F2)); figure; imshow(uint8(abs(I3))); title( 维纳滤波复原图 ); 得到的图像如下