summaryrefslogtreecommitdiffstats
path: root/demos
diff options
context:
space:
mode:
authoralgol <dkazanc@hotmail.com>2017-10-26 15:45:36 +0100
committeralgol <dkazanc@hotmail.com>2017-10-26 15:45:36 +0100
commit09f9bf9828c39bcdd870cfefbcb52e61451802eb (patch)
tree614e9badac151366da548683aa8565b674f91ce5 /demos
parente6349ed946772ee25a2f3ea1eba03e4c77e4f75e (diff)
downloadregularization-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.m2
-rw-r--r--demos/Demo_RealData3D_Parallel.m4
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;