diff options
| author | Daniil Kazantsev <dkazanc@hotmail.com> | 2017-08-24 12:46:12 +0100 |
|---|---|---|
| committer | Daniil Kazantsev <dkazanc@hotmail.com> | 2017-08-24 12:46:12 +0100 |
| commit | 69c2afc4670966c7e12420be7afd55ea69dfc8e8 (patch) | |
| tree | 2823622348828d99f4220d04b7ac11f4a84ebe86 /main_func | |
| parent | d7281b95f70dd3707f82640972a52f4155c2fde5 (diff) | |
| download | regularization-69c2afc4670966c7e12420be7afd55ea69dfc8e8.tar.gz regularization-69c2afc4670966c7e12420be7afd55ea69dfc8e8.tar.bz2 regularization-69c2afc4670966c7e12420be7afd55ea69dfc8e8.tar.xz regularization-69c2afc4670966c7e12420be7afd55ea69dfc8e8.zip | |
a minor big with TGV being constatly swtiched one removed
Diffstat (limited to 'main_func')
| -rw-r--r-- | main_func/FISTA_REC.m | 7 |
1 files changed, 4 insertions, 3 deletions
diff --git a/main_func/FISTA_REC.m b/main_func/FISTA_REC.m index 381fa34..a6e0ae5 100644 --- a/main_func/FISTA_REC.m +++ b/main_func/FISTA_REC.m @@ -98,7 +98,7 @@ else for i = 1:niter [id,x1] = astra_create_backprojection3d_cuda(sqweight.*y, proj_geomT, vol_geomT); s = norm(x1(:)); - x1 = x1/s; + x1 = x1./s; [sino_id, y] = astra_create_sino3d_cuda(x1, proj_geomT, vol_geomT); y = sqweight.*y; astra_mex_data3d('delete', sino_id); @@ -196,7 +196,7 @@ end if (isfield(params,'Regul_LambdaTGV')) LambdaTGV = params.Regul_LambdaTGV; else - LambdaTGV = 0.01; + LambdaTGV = 0; end if (isfield(params,'Ring_LambdaR_L1')) lambdaR_L1 = params.Ring_LambdaR_L1; @@ -369,6 +369,7 @@ if (subsets == 0) X2 = LLT_model(single(X), lambdaHO, tauHO, iterHO, 3.0e-05, 0); end X = 0.5.*(X + X2); % averaged combination of two solutions + end if (lambdaPB > 0) % Patch-Based regularization (can be slow on CPU) @@ -384,7 +385,7 @@ if (subsets == 0) if (LambdaTGV > 0) % Total Generalized variation (currently only 2D) lamTGV1 = 1.1; % smoothing trade-off parameters, see Pock's paper - lamTGV2 = 0.5; % second-order term + lamTGV2 = 0.8; % second-order term for kkk = 1:SlicesZ X(:,:,kkk) = TGV_PD(single(X(:,:,kkk)), LambdaTGV, lamTGV1, lamTGV2, IterationsRegul); end |