Lucy-Richarson滤波复原
I0=imread(‘src.png‘);
I=rgb2gray(I0);
subplot(341);
imshow(I);
title(‘src‘);
subplot(345);
imshow(I);
title(‘src‘);
subplot(349);
imshow(I);
title(‘src‘);
H=fspecial(‘motion‘,20,30); %运动卷积的偏移矩阵
MotionBlur=imfilter(I,H); %卷积
subplot(342);
imshow(MotionBlur);
title(‘MotionBlur‘)
G2=imnoise(MotionBlur,‘gaussian‘,0,0.005);
subplot(346);
imshow(G2);
title(‘Motion+noise‘);
LR1=deconvlucy(G2,H,5)
subplot(3,4,10);
imshow(LR1);
title(‘Motion+noise复原‘);
H1=fspecial(‘disk‘,10); %圆盘状偏移矩阵
disk=imfilter(I,H1); %卷积
subplot(343);
imshow(disk);
title(‘disk‘)
G3=imnoise(disk,‘gaussian‘,0,0.001);
subplot(347);
imshow(G3);
title(‘disk+noise‘);
LR2=deconvlucy(G3,H1,10)
subplot(3,4,11);
imshow(LR2);
title(‘disk+noise复原‘);
H2=fspecial(‘unsharp‘); %钝化模糊的偏移矩阵
unsharp=imfilter(I,H2); %卷积
subplot(344);
imshow(unsharp);
title(‘unsharp‘)
G4=imnoise(unsharp,‘gaussian‘,0,0.005);
subplot(348);
imshow(G4);
title(‘unsharp+noise‘);
LR3=deconvlucy(G4,H2,5)
subplot(3,4,12);
imshow(LR3);
title(‘unsharp+noise复原‘);
感觉只是对光学聚焦模糊有比较好的效果吧。。
原文:http://www.cnblogs.com/MnsterLu/p/5638683.html