MATLAB在图像复原中的应用研究(11)

IPOW=sum(IP(:)/numel(noise)); ICORR=fftshift(real(ifftn(IP))); ICORR1=ICORR(:,ceil(size(I,1)/2)); NSR=NPOW/IPOW;

subplot(221);imshow(BlurredNoisy,[]);

title('模糊和噪声图像');

subplot(222);imshow(deconvwnr(BlurredNoisy,PSF,NSR),[]); title('deconbwnr(A,PSF,NSR)');

subplot(223);imshow(deconvwnr(BlurredNoisy,PSF,NCORR,ICORR),[]); title('deconbwnr(A,PSF,NCORR,ICORR)');

subplot(224);imshow(deconvwnr(BlurredNoisy,PSF,NPOW,ICORR1),[]); title('deconbwnr(A,PSF,NPOW,ICORR_1_D)'); （2）.规则化滤波复原程序源代码：

I=checkerboard(8);

PSF=fspecial('gaussian',7,10); V=.01;

BlurredNoisy=imnoise(imfilter(I,PSF),'gaussian',0,V); NOISEPOWER=V*numel(I);

[J LAGRA]=deconvreg(BlurredNoisy,PSF,NOISEPOWER); subplot(221);imshow(BlurredNoisy); title('A=Blurred and Noisy'); subplot(222);imshow(J);

title('[J LAGRA]=deconvreg(A,PSF,NP)');

subplot(223);imshow(deconvreg(BlurredNoisy,PSF,[],LAGRA/10)); title('deconvreg(A,PSF,[],0.1*LAGRA)');

subplot(224);imshow(deconvreg(BlurredNoisy,PSF,[],LAGRA/10)) title('deconvreg(A,PSF,[],10*LAGRA'); （3）.Lucy-Richardson复原滤波源代码：

I=checkerboard(8);

PSF=fspecial('gaussian',7,10);

V=.0001;

BlurredNoisy=imnoise(imfilter(I,PSF),'gaussian',0,V); WT=zeros(size(I)); WT(5:end-4,5:end-4)=1;

J1=deconvlucy(BlurredNoisy,PSF);

J2=deconvlucy(BlurredNoisy,PSF,20,sqrt(V));

J3=deconvlucy(BlurredNoisy,PSF,20,sqrt(V),[],WT); subplot(221);imshow(BlurredNoisy); title('A=Blurred and Noisy');