盲去卷积

工具箱通过函数deconvblind来执行盲去卷积，它有如下语法：

[f,PSFe]=deconvblind(g,INITPSF)

其中，g代表退化函数，INITPSF是点扩散函数的出事估计。PSFe是这个函数最终计算到的估计值，fr是利用估计的PSF复原的图像。 用来去的复原图像的算法是L-R迭代复原算法[13]。PSF估计受其初始推测尺寸的巨大影响，而很少受其值的影响。

盲去卷积程序：

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;

INITPSF=ones(size(PSF));

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

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

subplot(221);

imshow(BlurredNoisy);

title('A=Blurred and Noisy');

subplot(222);

imshow(PSF,[]);

title('True PSF');

subplot(223);

imshow(J);

title('Deblured Image');

subplot(224);

imshow(P,[]);

模拟实验结果如下：

该算法优点是,同时恢复了图像和点扩张函数,在对失真情况毫无先验知识的情况下, 仍能实现对模糊图像的恢复操作。利用 MATLAB实现的图像恢复, 并对恢复图像的失真情况做了改善。在进行图像恢复时,重建 PSF,对图像进行重建, 得到恢复的图像。