Major reorganization, updated routines

1 files changed, 0 insertions, 99 deletions

diff --git a/demo/DemoRD1.m b/demo/DemoRD1.m
deleted file mode 100644
index 9a43cb5..0000000
--- a/demo/DemoRD1.m
+++ /dev/null
@@ -1,99 +0,0 @@
-% Demonstration of tomographic reconstruction from neutron tomography
-% dataset (basalt sample) using Student t data fidelity 
-clear all
-close all 
-
-% adding paths
-addpath('data/'); 
-addpath('main_func/'); 
-addpath('supp/'); 
-
-load('sino_basalt.mat') % load real neutron data
-
-size_det = size(sino_basalt, 1); % detector size
-angSize = size(sino_basalt,2); % angles dim
-recon_size = 650; % reconstruction size
-
-FBP = iradon(sino_basalt, rad2deg(angles),recon_size);
-figure; imshow(FBP , [0, 0.45]); title ('FBP reconstruction');
-
-%%
-fprintf('%s\n', 'Reconstruction using FISTA-LS without regularization...');
-clear params
-params.sino = sino_basalt';
-params.N  = recon_size;
-params.angles = angles;
-params.iterFISTA  = 50;
-params.show = 0;
-params.maxvalplot = 0.6; params.slice = 1;
-
-tic; [X_fista] = FISTA_REC(params); toc;
-figure; imshow(X_fista , [0, 0.45]); title ('FISTA-LS reconstruction');
-%%
-fprintf('%s\n', 'Reconstruction using FISTA-LS-TV...');
-clear params
-params.sino = sino_basalt';
-params.N  = recon_size;
-params.angles = angles;
-params.iterFISTA  = 150;
-params.lambdaTV = 0.0003; % TV regularization parameter
-params.tol = 1.0e-04;
-params.iterTV = 20;
-params.show = 1;
-params.maxvalplot = 0.6; params.slice = 1;
-
-tic; [X_fista_TV] = FISTA_REC(params); toc;
-figure; imshow(X_fista_TV , [0, 0.45]); title ('FISTA-LS-TV reconstruction');
-%%
-%%
-fprintf('%s\n', 'Reconstruction using FISTA-GH-TV...');
-clear params
-params.sino = sino_basalt';
-params.N  = recon_size;
-params.angles = angles;
-params.iterFISTA  = 350;
-params.lambdaTV = 0.0003; % TV regularization parameter
-params.tol = 1.0e-04;
-params.iterTV = 20;
-params.lambdaR_L1 = 0.001; % Soft-Thresh L1 ring variable parameter
-params.show = 1;
-params.maxvalplot = 0.6; params.slice = 1;
-
-tic; [X_fista_GH_TV] = FISTA_REC(params); toc;
-figure; imshow(X_fista_GH_TV , [0, 0.45]); title ('FISTA-GH-TV reconstruction');
-%%
-%%
-fprintf('%s\n', 'Reconstruction using FISTA-Student-TV...');
-clear params
-params.sino = sino_basalt';
-params.N  = recon_size;
-params.angles = angles;
-params.iterFISTA  = 350;
-params.L_const = 7000; % Lipshitz constant
-params.lambdaTV = 0.0003; % TV regularization parameter
-params.tol = 1.0e-04;
-params.iterTV = 20;
-params.fidelity = 'student'; % choosing Student t penalty
-params.show = 1;
-params.maxvalplot = 0.6; params.slice = 1;
-
-tic; [X_fistaStudentTV] = FISTA_REC(params); toc;
-figure; imshow(X_fistaStudentTV , [0, 0.45]); title ('FISTA-Student-TV reconstruction');
-%%
-
-fprintf('%s\n', 'Segmentation using OTSU method ...');
-level = graythresh(X_fista);
-Segm_FISTA = im2bw(X_fista,level);
-figure; imshow(Segm_FISTA, []); title ('Segmented FISTA-LS reconstruction');
-
-level = graythresh(X_fista_TV);
-Segm_FISTA_TV = im2bw(X_fista_TV,level);
-figure; imshow(Segm_FISTA_TV, []); title ('Segmented FISTA-LS-TV reconstruction');
-
-level = graythresh(X_fista_GH_TV);
-BW_FISTA_GH_TV = im2bw(X_fista_GH_TV,level);
-figure; imshow(BW_FISTA_GH_TV, []); title ('Segmented FISTA-GH-TV reconstruction');
-
-level = graythresh(X_fistaStudentTV);
-BW_FISTA_Student_TV = im2bw(X_fistaStudentTV,level);
-figure; imshow(BW_FISTA_Student_TV, []); title ('Segmented FISTA-Student-LS reconstruction');
\ No newline at end of file