diff options
| author | Daniil Kazantsev <dkazanc@hotmail.com> | 2017-07-03 22:57:55 +0100 |
|---|---|---|
| committer | Daniil Kazantsev <dkazanc@hotmail.com> | 2017-07-03 22:57:55 +0100 |
| commit | 75e7ebc79e86710c4cfb05b8af8e429146288c85 (patch) | |
| tree | 3d6abaccb81c44df111605474c117fb2bc6570a9 /demos | |
| parent | 329a104d4cb5ba50a59fb80e58de0453ba49f075 (diff) | |
| download | regularization-75e7ebc79e86710c4cfb05b8af8e429146288c85.tar.gz regularization-75e7ebc79e86710c4cfb05b8af8e429146288c85.tar.bz2 regularization-75e7ebc79e86710c4cfb05b8af8e429146288c85.tar.xz regularization-75e7ebc79e86710c4cfb05b8af8e429146288c85.zip | |
FISTA-GH bug fixed
Diffstat (limited to 'demos')
| -rw-r--r-- | demos/DemoRD1.m | 2 | ||||
| -rw-r--r-- | demos/DemoRD2.m | 12 |
2 files changed, 7 insertions, 7 deletions
diff --git a/demos/DemoRD1.m b/demos/DemoRD1.m index c25bb3e..e190d34 100644 --- a/demos/DemoRD1.m +++ b/demos/DemoRD1.m @@ -75,7 +75,7 @@ params.L_const = 3500; % Lipshitz constant params.Regul_LambdaTV = 0.0003; % TV regularization parameter params.fidelity = 'student'; % choosing Student t penalty params.show = 1; -params.initilize = 1; % warm start with SIRT +params.initialize = 1; % warm start with SIRT params.maxvalplot = 0.6; params.slice = 1; tic; [X_fistaStudentTV] = FISTA_REC(params); toc; diff --git a/demos/DemoRD2.m b/demos/DemoRD2.m index ab4da96..518e24b 100644 --- a/demos/DemoRD2.m +++ b/demos/DemoRD2.m @@ -75,6 +75,7 @@ params.proj_geom = proj_geom; % pass geometry to the function params.vol_geom = vol_geom; params.sino = Sino3D; params.iterFISTA = 40; +params.L_const = 7.6789e+08; params.Regul_LambdaTV = 0.005; % TV regularization parameter for FISTA-TV params.Ring_LambdaR_L1 = 0.002; % Soft-Thresh L1 ring variable parameter params.Ring_Alpha = 21; % to boost ring removal procedure @@ -108,17 +109,16 @@ figure; imshow(X_fista_GH_TVLLT(:,:,1) , [0, 2.5]); title ('FISTA-GH-TV-LLT reco %% % fprintf('%s\n', 'Reconstruction using FISTA-Student-TV...'); % clear params +% params.proj_geom = proj_geom; % pass geometry to the function +% params.vol_geom = vol_geom; % params.sino = Sino3D(:,:,10); -% params.N = recon_size; -% params.angles = angles_rad; -% params.iterFISTA = 100; +% params.iterFISTA = 50; % params.L_const = 0.01; % Lipshitz constant -% params.lambdaTV = 0.006; % TV regularization parameter for FISTA-TV -% params.tol = 1.0e-04; -% params.iterTV = 20; +% params.Regul_LambdaTV = 0.008; % TV regularization parameter for FISTA-TV % params.fidelity = 'student'; % choosing Student t penalty % params.weights = Weights3D(:,:,10); % params.show = 0; +% params.initialize = 1; % params.maxvalplot = 2.5; params.slice = 1; % % tic; [X_fistaStudentTV] = FISTA_REC(params); toc; |