对图像进行模糊处理-然后加白高斯噪声-用逆滤波和维纳滤波恢复图像

作业六用式（5.6-3)对图像进行模糊处理，然后加白高斯噪声，用逆滤波和维纳滤波恢复图像，对结果进行分析。程序代码如下：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=[];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)));figure;imshow(uint8(I0));5/622(,)kuvHuve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=[];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));figure;imshow(uint8(abs(I3)));title('维纳滤波复原图');得到的图像如下