第06章 图像复原2.pptVIP

第6章 图像复原(2)image restoration ——运动模糊复原 匀速运动引起的图像模糊 由匀速运动引起的模糊 由匀速运动引起的模糊 由匀速运动引起的模糊 edgetaper(I,PSF)函数 clear all; close all I = im2double(imread(cameraman.tif)); imshow(I); title(Original Image (courtesy of MIT)); % Simulate a motion blur. LEN = 21;THETA = 11; PSF = fspecial(motion, LEN, THETA); blurred = imfilter(I, PSF, conv, circular); % Simulate additive noise. noise_mean = 0; noise_var = 0.0001; blurred_noisy = imnoise(blurred, gaussian, noise_mean, noise_var); figure, imshow(blurred_noisy),title(Simulate Blur and Noise) % Try restoration assuming no noise. estimated_nsr = 0; wnr2 = deconvwnr(blurred_noisy, PSF, estimated_nsr); figure, imshow(wnr2) title(Restoration of Blurred, Noisy Image Using NSR = 0) %Try restoration using a better estimate of the noise-to-signal-power ratio. estimated_nsr = noise_var / var(I(:)); wnr3 = deconvwnr(blurred_noisy, PSF, estimated_nsr); figure, imshow(wnr3) title(Restoration of Blurred, Noisy Image Using Estimated NSR); 维纳滤波复原MATLAB编程 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(:)/numel(noise)); NCORR=fftshift(real(ifftn(NP))); IP=abs(fftn(I)).^2; 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)); 例6.1 维纳滤波复原的具体实现方法 【例6.1】原始图像如图6.2（a），使用函数DECONVWNR对图6.2（b）所示的无噪声模糊图像进行复原重建，观察所得结果，并将不同PSF产生的复原效果进行比较。 图6.3 不同PSF产生的复原效果比较 Lucy-Richardson复原程序举例 %Lucy-Richardson复原程序I=checkerboard(8);PSF=fspecial(gaussian,7,10);V=.0001;BlurredNoisy=imnoise(imfilter(I,PSF),gaussian,0,V);WT=zeros(siz