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