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

作业六

用式（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=[];

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));

()5/6

22(,)k u v H u v e -+=

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('维纳滤波复原图');

得到的图像如下