diff options
author | algol <dkazanc@hotmail.com> | 2017-10-26 15:45:36 +0100 |
---|---|---|
committer | algol <dkazanc@hotmail.com> | 2017-10-26 15:45:36 +0100 |
commit | 09f9bf9828c39bcdd870cfefbcb52e61451802eb (patch) | |
tree | 614e9badac151366da548683aa8565b674f91ce5 /demos | |
parent | e6349ed946772ee25a2f3ea1eba03e4c77e4f75e (diff) | |
download | regularization-09f9bf9828c39bcdd870cfefbcb52e61451802eb.tar.gz regularization-09f9bf9828c39bcdd870cfefbcb52e61451802eb.tar.bz2 regularization-09f9bf9828c39bcdd870cfefbcb52e61451802eb.tar.xz regularization-09f9bf9828c39bcdd870cfefbcb52e61451802eb.zip |
Lipshitz constant attached to the regularization parameter
Diffstat (limited to 'demos')
-rw-r--r-- | demos/Demo_Phantom3D_Parallel.m | 2 | ||||
-rw-r--r-- | demos/Demo_RealData3D_Parallel.m | 4 |
2 files changed, 3 insertions, 3 deletions
diff --git a/demos/Demo_Phantom3D_Parallel.m b/demos/Demo_Phantom3D_Parallel.m index b0737f1..4219bd1 100644 --- a/demos/Demo_Phantom3D_Parallel.m +++ b/demos/Demo_Phantom3D_Parallel.m @@ -71,7 +71,7 @@ clear params params.proj_geom = proj_geom; % pass geometry to the function
params.vol_geom = vol_geom;
params.sino = single(sino3D_log); % sinogram
-params.iterFISTA = 12; %max number of outer iterations
+params.iterFISTA = 15; %max number of outer iterations
params.X_ideal = TomoPhantom; % ideal phantom
params.weights = dataRaw./max(dataRaw(:)); % statistical weight for PWLS
params.subsets = 12; % the number of subsets
diff --git a/demos/Demo_RealData3D_Parallel.m b/demos/Demo_RealData3D_Parallel.m index d680aca..f82e0b0 100644 --- a/demos/Demo_RealData3D_Parallel.m +++ b/demos/Demo_RealData3D_Parallel.m @@ -72,12 +72,12 @@ fprintf('%s\n', 'Reconstruction using FISTA-OS-GH-TV...'); clear params params.proj_geom = proj_geom; % pass geometry to the function params.vol_geom = vol_geom; -params.sino = Sino3D; +params.sino = Sino3D(:,:,10); params.iterFISTA = 18; params.Regul_Lambda_FGPTV = 5.0000e+6; % TV regularization parameter for FGP-TV params.Ring_LambdaR_L1 = 0.002; % Soft-Thresh L1 ring variable parameter params.Ring_Alpha = 21; % to boost ring removal procedure -params.weights = Weights3D; +params.weights = Weights3D(:,:,10); params.subsets = 8; % the number of ordered subsets params.show = 1; params.maxvalplot = 2.5; params.slice = 1; |