图像复原方法综述(3)

(4) 盲去卷积

clc

clear all close all

I=imread('C:\\Documents and Settings\\Administrator\\桌面\\原始图.jpg'); subplot(231); imshow(I); I=rgb2gray(I); subplot(232); imshow(I);

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

BlurredNoisy=imnoise(imfilter(I,PSF),'gaussian',0,V); BlurredNoisy=double(BlurredNoisy); WT=zeros(size(I));

WT(5:end-4,5:end-4)=1; INITPSF=ones(size(PSF));

FUN=inline('PSF+P1','PSF','P1');

[J P]=deconvblind(BlurredNoisy,INITPSF,5,10*sqrt(V), WT,FUN,0);

[K P]=deconvblind(BlurredNoisy,INITPSF,10,10*sqrt(V), WT,FUN,0);

[L P]=deconvblind(BlurredNoisy,INITPSF,20,10*sqrt(V), WT,FUN,0); subplot(233);imshow(mat2gray(BlurredNoisy)); title('A=Blurred and Noisy');

subplot(234);imshow(mat2gray(J)); title('True PSF');

subplot(235);imshow(mat2gray(K)); title('Deblured Image');

subplot(236);imshow(mat2gray(L)); title('Recovered PSF');