Lucy-Richardson

Lucy Richardson algorithm for image restoration

clc;
close all;
clear all

f = im2double（imread（‘D:\vide\Code\Weiner\cameraman.tif‘））;
subplot（321）
imshow（f）
title（‘Orignal Image‘）;

PSF = fspecial（‘gaussian‘720）;
f = edgetaper（fPSF）;
SD = 0.001;
g = imnoise（imfilter（fPSF）‘gaussian‘0SD^2）;
DAMPAR = SD;
LIM = ceil（size（PSF1）/2）;
WEIGHT = zeros（size（g））;
WEIGHT（LIM+1:end-1LIM+1:end-1）=1;


subplot（322）
imshow（g）
title（‘Image blurrd using Gaussian blurr+Gaussian noise of mean 0 and var .0001 ‘）;
NUMIT=5;
fr = deconvlucy（gPSFNUMITDAMPARWEIGHT）;
subplot（323）
imshow（fr）
title（‘Image restored using 5 iterations‘）;
NUMIT=10;
fr = deconvlucy（gPSFNUMITDAMPARWEIGHT）;
subplot（324）
imshow（fr）
title（‘Image restored using 10 iterations‘）;
NUMIT=15;
fr = deconvlucy（gPSFNUMITDAMPARWEIGHT）;
subplot（325）
imshow（fr）
title（‘Image restored using 15 iterations‘）;
NUMIT=20;
fr = deconvlucy（gPSFNUMITDAMPARWEIGHT）;
subplot（326）
imshow（fr）
title（‘Image restored using 20 iterations‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      36336  2012-08-06 12:41  Lucy-Richardson\gaussian_deblur.jpg

     文件      97245  2012-08-06 09:52  Lucy-Richardson\image_0.jpg

     文件      96527  2012-08-06 09:53  Lucy-Richardson\image_1.jpg

     文件      95800  2012-08-06 09:53  Lucy-Richardson\image_2.jpg

     文件      95515  2012-08-06 09:54  Lucy-Richardson\image_3.jpg

     文件      96153  2012-08-06 09:54  Lucy-Richardson\image_4.jpg

     文件       1007  2012-08-14 12:53  Lucy-Richardson\L_R.asv

     文件       1003  2012-08-16 17:03  Lucy-Richardson\L_R.m

     文件       1131  2012-09-28 10:24  Lucy-Richardson\L_R_motion.asv

     文件       1131  2012-09-28 10:26  Lucy-Richardson\L_R_motion.m

     文件        246  2012-09-28 08:57  Lucy-Richardson\PSNR_calc.m

    ..A.SH.      8192  2012-10-11 13:38  Lucy-Richardson\Thumbs.db

     目录          0  2012-10-11 13:38  Lucy-Richardson

----------- ---------  ---------- -----  ----

               530286                    13