图像恢复技术之Matlab中的维纳滤波代码复现

k=0.001;
blur=g;
bf = fftshift(fft2(blur));
[M,N]=size(blur);
H=zeros(M,N);
rad= 5;
length = 2*rad+1;
psf=zeros(length,length);
for i=1:length
for j=1:length

if((  (i-rad-1)*(i-rad-1) + (j-rad-1)*(j-rad-1) ) < rad*rad )
H(i,j) =1.0/(rad*rad*3.14159);
psf(i,j) =1.0/(rad*rad*3.14159);
end

end
end

H= fftshift(fft2(H));
Hg = conj(H);
spec = Hg.*H;
spectt = spec +k;
can = spec./spectt;

jk = can./H;
Gg=bf.*jk;

Gg=ifftshift(Gg);
pp=im2uint8(mat2gray(real(ifft2(Gg))));%退化后的图像

subplot(2,2,1);
imshow(blur);
subplot(2,2,2);
imshow(pp);