第6章-图像复原.pptVIP

fr1=deconvwnr(C,PSF); sn=abs(fft2(noise).^2); % noise power spectrum nA=sum(sn(:))/prod(size(noise)); % noise average power sf=abs(fft2(A)).^2 % image power spectrum fA=sum(sf(:))/prod(size(A)); % image average power R=nA/fA; fr2=deconvwnr(C,PSF,R); NCORR=fftshift(real(ifft2(sn))); ICORR=fftshift(real(ifft2(sf))); fr3=deconvwnr(C,PSF,NCORR,ICORR); figure(2); subplot(2,2,1);imshow(uint8(C),[]);title(模糊噪声图像); subplot(2,2,2);imshow(uint8(fr1),[]);title(直接逆滤波); subplot(2,2,3);imshow(uint8(fr2),[]);title(常数比率维纳滤波); subplot(2,2,4);imshow(uint8(fr3),[]);title(使用自相关函数的维纳滤波); 2）三种滤波方式复原图像 五 经典复原方法