基于matlab的退化图像复原（三）------维纳滤波复原

还是照着书上推了一遍，

反正写一遍之后很好理解，

也许是只是这点算法比较简单基础。。。

讲到底还是统计的知识，二维的最小二乘法，

让均方差最小

 

放代码：

I=imread('src.png');
I0=rgb2gray(I);
subplot(151);
imshow(I0);
title('src');
H=fspecial('motion',30,45); %运动卷积的偏移矩阵
MotionBlur=imfilter(I0,H,0,'conv'); %卷积
MF=imfilter(I0,H,'circular','conv');%图像大小通过将图像看成是一个二维周期函数的一个周期来扩展
subplot(152);
imshow(MotionBlur);
title('MotionBlur')
subplot(153);
imshow(MF);
title('MF')
wnr=deconvwnr(MotionBlur,H);
subplot(154);
imshow(wnr)
title('MotionBlur复原')
wnr1=deconvwnr(MF,H);
subplot(155);
imshow(wnr1)
title('MF复原')

　　

有个问题就是imfilter的补充矩阵，里面的circular是按周期来补矩阵，这个复原效果比较好。

直接补0，是MotionBlur这个效果。

 参考《MATLAB图像处理超级学习手册》