diff options
| author | Daniil Kazantsev <dkazanc@hotmail.com> | 2019-03-16 23:39:29 +0000 |
|---|---|---|
| committer | Daniil Kazantsev <dkazanc@hotmail.com> | 2019-03-16 23:39:29 +0000 |
| commit | 9633c5d701c164c3ff8f4d870624f87744d186bd (patch) | |
| tree | 23fd9a110436318aa55cbb46a6f86824940b5478 /demos/demoMatlab_3Ddenoise.m | |
| parent | 1740fb91c1842dc174492857fad88c03420f74eb (diff) | |
| download | regularization-9633c5d701c164c3ff8f4d870624f87744d186bd.tar.gz regularization-9633c5d701c164c3ff8f4d870624f87744d186bd.tar.bz2 regularization-9633c5d701c164c3ff8f4d870624f87744d186bd.tar.xz regularization-9633c5d701c164c3ff8f4d870624f87744d186bd.zip | |
matlab demos updated
Diffstat (limited to 'demos/demoMatlab_3Ddenoise.m')
| -rw-r--r-- | demos/demoMatlab_3Ddenoise.m | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/demos/demoMatlab_3Ddenoise.m b/demos/demoMatlab_3Ddenoise.m index 6b21e86..3942eea 100644 --- a/demos/demoMatlab_3Ddenoise.m +++ b/demos/demoMatlab_3Ddenoise.m @@ -145,9 +145,10 @@ fprintf('Denoise using the TGV model (CPU) \n'); lambda_TGV = 0.03; % regularisation parameter alpha1 = 1.0; % parameter to control the first-order term alpha0 = 2.0; % parameter to control the second-order term +L2 = 12.0; % convergence parameter iter_TGV = 500; % number of Primal-Dual iterations for TGV epsil_tol = 0.0; % tolerance -tic; u_tgv = TGV(single(vol3D), lambda_TGV, alpha1, alpha0, iter_TGV, epsil_tol); toc; +tic; u_tgv = TGV(single(vol3D), lambda_TGV, alpha1, alpha0, iter_TGV, L2, epsil_tol); toc; rmseTGV = RMSE(Ideal3D(:),u_tgv(:)); fprintf('%s %f \n', 'RMSE error for TGV is:', rmseTGV); figure; imshow(u_tgv(:,:,3), [0 1]); title('TGV denoised volume (CPU)'); @@ -157,7 +158,7 @@ figure; imshow(u_tgv(:,:,3), [0 1]); title('TGV denoised volume (CPU)'); % alpha1 = 1.0; % parameter to control the first-order term % alpha0 = 2.0; % parameter to control the second-order term % iter_TGV = 500; % number of Primal-Dual iterations for TGV -% tic; u_tgv_gpu = TGV_GPU(single(vol3D), lambda_TGV, alpha1, alpha0, iter_TGV, epsil_tol); toc; +% tic; u_tgv_gpu = TGV_GPU(single(vol3D), lambda_TGV, alpha1, alpha0, iter_TGV, L2, epsil_tol); toc; % rmseTGV = RMSE(Ideal3D(:),u_tgv_gpu(:)); % fprintf('%s %f \n', 'RMSE error for TGV is:', rmseTGV); % figure; imshow(u_tgv_gpu(:,:,3), [0 1]); title('TGV denoised volume (GPU)'); |