I = checkerboard(8);
% 創建一個棋盤圖像
noise = 0.1*randn(size(I));
% 產生隨機噪聲
PSF = fspecial('motion',21,11);
Blurred = imfilter(I,PSF,'circular');
% 輸入圖像濾波,使圖像變模糊
BlurredNoisy = im2uint8(Blurred + noise);
% 添加噪聲
NP = abs(fftn(noise)).^2;
% 噪聲功率
NPOW = sum(NP(:))/prod(size(noise));
NCORR = fftshift(real(ifftn(NP)));
% 噪聲的自相關
IP = abs(fftn(I)).^2;
% 原始圖像功率
IPOW = sum(IP(:))/prod(size(I));
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('deconvwnr(A,PSF,NSR)');
subplot(223);imshow(deconvwnr(BlurredNoisy,PSF,NCORR,ICORR),[]);
title('deconvwnr(A,PSF,NCORR,ICORR)');
subplot(224);imshow(deconvwnr(BlurredNoisy,PSF,NPOW,ICORR1),[]);
title('deconvwnr(A,PSF,NPOW,ICORR_1_D)');