数字图像处理第五章

I_1=imread('lena.jpg');
PSF=fspecial('motion',21,11);
I3=imfilter(I_1,PSF,'conv','circular');
I3=imnoise(I3,'gaussian',0,0.0001);
I4=deconvwnr(I3,PSF);
I_1=im2double(I_1);
estimated_nsr=0.0001/var(I_1(:));
I5=deconvwnr(I3,PSF,estimated_nsr);
figure;
subplot(2,2,1);
imshow(I_1);
title('原图');
subplot(2,2,2);
imshow(I3);
title('模糊+噪声');
subplot(2,2,3);
imshow(I4);
title('逆滤波');
subplot(2,2,4);
imshow(I5);
title('winner滤波');