基于matlab的退化图像复原（四）------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复原‘);