clear?all
close?all

%?set?up?cameram?man?example
I?=?double（imread（‘cameraman.tif‘））;
H?=?fspecial（‘a(chǎn)verage‘9）;
%H=fspecial（‘Gaussian‘212）;

[MN]=size（I）;

%%?Add?blur?and?noise
B?=?imfilter（IH‘circular‘）;
sigma=3;
Bn?=?B+randn（MN）*sigma;

%%?Pre-compute?the?weight?by?tikronov?optimal?reconstruction
[Im0?sigma1?lambdaerr]?=?tikronov_optimal_lambda（Bn?H?sigma）;
%?TV?reconstruction
opts.bTol=0.99*sigma;?%?stop?when?||Au-g||\leq?btol*sigma?for?noise?case
opts.Init=Im0;
opts.I=I;
opts.nWeightUpdate=0;


%%
time?=?cputime;
opts.type=1;
opts.mu=5;?%
?[uTV_bosenergyrelmsepsnr_n]=Decon_TV（BnHopts）;
?time=cputime-time;

%uTV=deconvrof（BnH11000.01）;
?
%%
opts.h0=2*sigma1;
opts.nDenoising=20;
opts.type=1;
opts.mu=10;?%
time_nltv_bos=cputime;
?[uNLTV_bosenergyrelmsepsnr_n_nltv_bos]=Decon_NLTV（BnHopts）;
?time_nltv_bos=cputime-time_nltv_bos
%%
opts.type=2;
opts.mu=20;
opts.epsilon=0.1;
time_nltv_pbos=cputime;
?[uNLTV_pbosenergyrelmsepsnr_n_nltv_pbos]=Decon_NLTV（BnHopts）;
?time_nltv_pbos=cputime-time_nltv_pbos
%%?
?time?=?cputime;
opts.type=1;
opts.mu=5;
?[uNLH1_bosenergyrelmsepsnr_n]=Decon_NLH1（BnHopts）;
?time=cputime-time;

?%%
opts.type=2;
opts.mu=10;
opts.epsilon=0.1;
opts.nWeightUpdate=0;
time?=?cputime;
?[uNLH1_pbosenergyrelmsepsnr_n]=Decon_NLH1（BnHopts）;
?time=cputime-time;
?
%%?Other?comparison?
?u_GPSR=GPSR_decon（BnH[]）;
%%?gradient?descent?for?nltv?
time_gd=cputime;
?w=nlmeans_weight_sym（Im0sigma1210）;
?[uNLTV_gdkenergyNLTVtime]?=?nonlocalTV_decon（BnHw15500Im0）;
?time_gd=cputime-time_gd
%%
?figure?
subplot（331）imshow（I[0?255?]）title（‘Original‘）;
subplot（332）imshow（Bn[0?255?]）title（‘Blurry‘）;
%subplot（333）imshow（Im0[0?255?]）title（[‘Initial?guess?by?setting?unknown?0?PSNR=‘num2str（PSNR（IIm0））]）;
subplot（333）imshow（uNLTV_gd[0?255?]）title（[‘NLTV?with?Gradient?descent?PSNR=‘num2str（PSNR（IuNLTV_gd））]）;
subplot（334）imshow（uTV_bos[0?255?]）title（[‘TV?PSNR=‘num2str（PSNR（IuTV_bos））]）;

subplot（335）imshow（uNLH1_bos[0?255??]）title（[‘NLH1?BOS?PSNR=‘num2str（PSNR（IuNLH1_bos））]）;
subplot（336）imshow（uNLH1_pbos[0?255?]）title（[‘NLH1?PBOS?PSNR=‘num2str（PSNR（IuNLH1_pbos））]）;
subplot（337）imshow（uNLTV_bos[0?255?]）title（[‘NLTV?BOS?PSNR=‘num2str（PSNR（IuNLTV_bos））]）;
subplot（338）imshow（uNLTV_pbos[0?255?]）title（[‘NLTV?PBOS?PSNR=‘num2str（PSNR（IuNLTV_pbos））]）;
subplot（339）imshow（u_GPSR[0?255?]）title（[‘GPSR??PSNR=‘num2str（PSNR（Iu_GPSR））]）;

%%
if?0
name=‘cameraman_deblur‘;
dir=‘\results‘;
gmin=0;
gmax=255;
image_name=[name‘_original.eps‘];??
image2eps（Iimage_namedirgmingmax）;


image_name=[name‘_blur.eps‘];??
image2eps（Bnimage_namedirgmingmax）;

image_name=[name‘_init.eps‘];
image2eps（Im0image_namedirgmingmax）;


image_name=[name‘_tv.eps‘];
image2eps（uTV_bosimage_namedirgmingmax）;

image_name=[name‘_nltv_bos.eps‘];
image2eps（uNLTV_bosimage_namedirgmingmax）;

image_name=[name‘_nltv_pbos.eps‘];
image2eps（uNLTV_pbosimage_namedirgmingmax）;

image_na

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件?????432086??2008-11-20?18:20??BOS_NLTV_v1\barbar_test.mat

?????文件???????3410??2009-07-16?17:56??BOS_NLTV_v1\figure_1.m

?????文件???????5276??2009-07-16?14:05??BOS_NLTV_v1\figure_2.m

?????文件???????5381??2009-07-16?17:54??BOS_NLTV_v1\figure_3.m

?????文件???????6148??2007-12-04?11:46??BOS_NLTV_v1\GPSR_5.0\.DS_Store

?????文件??????65126??2007-01-25?13:54??BOS_NLTV_v1\GPSR_5.0\Camera.tif

?????文件???????7046??2007-09-23?20:49??BOS_NLTV_v1\GPSR_5.0\demo_continuation.m

?????文件???????5857??2008-11-21?22:51??BOS_NLTV_v1\GPSR_5.0\demo_image_deblur.asv

?????文件???????5857??2008-11-21?22:54??BOS_NLTV_v1\GPSR_5.0\demo_image_deblur.m

?????文件???????8489??2007-10-19?16:40??BOS_NLTV_v1\GPSR_5.0\figures_1_2_3.m

?????文件???????5225??2007-09-12?18:21??BOS_NLTV_v1\GPSR_5.0\figure_4.m

?????文件???????5276??2007-12-04?11:46??BOS_NLTV_v1\GPSR_5.0\figure_5.m

?????文件???????2377??2007-07-12?18:41??BOS_NLTV_v1\GPSR_5.0\figure_6.m

?????文件??????21851??2007-12-04?11:26??BOS_NLTV_v1\GPSR_5.0\GPSR_Basic.m

?????文件??????21647??2007-12-04?11:08??BOS_NLTV_v1\GPSR_5.0\GPSR_BB.m

?????文件??????12630??2007-09-06?15:10??BOS_NLTV_v1\GPSR_5.0\IST.m

?????文件???????8569??2007-06-16?15:41??BOS_NLTV_v1\GPSR_5.0\l1_ls.m

?????文件???????3178??2007-12-05?19:06??BOS_NLTV_v1\GPSR_5.0\README.txt

?????文件??????43807??2009-02-28?15:36??BOS_NLTV_v1\images\barbara256.png

?????文件??????65126??2007-01-25?13:54??BOS_NLTV_v1\images\Camera.tif

?????文件??????76763??2008-08-15?14:26??BOS_NLTV_v1\images\patches.gif

????..A.SH.?????30720??2009-06-09?13:58??BOS_NLTV_v1\images\Thumbs.db

?????文件??????17457??2009-07-16?15:00??BOS_NLTV_v1\mex\compute_nl_weights_mex.c

?????文件??????12800??2009-07-16?14:53??BOS_NLTV_v1\mex\compute_nl_weights_mex.mexw32

?????文件???????5160??2009-07-02?10:25??BOS_NLTV_v1\mex\NLH1_mex.c

?????文件???????7168??2009-07-16?14:53??BOS_NLTV_v1\mex\NLH1_mex.mexw32

?????文件???????4485??2008-09-02?15:38??BOS_NLTV_v1\mex\nlmeans_weight_sym.c

?????文件???????7680??2009-07-16?14:53??BOS_NLTV_v1\mex\nlmeans_weight_sym.mexw32

?????文件???????8713??2008-11-20?16:58??BOS_NLTV_v1\mex\SBNLTV_mex.c

?????文件???????8704??2009-07-16?14:53??BOS_NLTV_v1\mex\SBNLTV_mex.mexw32

............此處省略77個(gè)文件信息