summaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authoralgol <dkazanc@hotmail.com>2017-10-20 16:16:32 +0100
committeralgol <dkazanc@hotmail.com>2017-10-20 16:16:32 +0100
commitc91436873e48d531b9313f9c10fa5f89bcb90ab6 (patch)
treeb36d41e8508de08e13e468774cc0df569b03d162
parent952322eb1394f18cedfbb51488e76ed52728d98d (diff)
downloadregularization-c91436873e48d531b9313f9c10fa5f89bcb90ab6.tar.gz
regularization-c91436873e48d531b9313f9c10fa5f89bcb90ab6.tar.bz2
regularization-c91436873e48d531b9313f9c10fa5f89bcb90ab6.tar.xz
regularization-c91436873e48d531b9313f9c10fa5f89bcb90ab6.zip
demo_parallel update
-rw-r--r--demos/Demo_Phantom3D_Parallel.m12
1 files changed, 12 insertions, 0 deletions
diff --git a/demos/Demo_Phantom3D_Parallel.m b/demos/Demo_Phantom3D_Parallel.m
index 9ef1f3c..402bdd2 100644
--- a/demos/Demo_Phantom3D_Parallel.m
+++ b/demos/Demo_Phantom3D_Parallel.m
@@ -32,6 +32,18 @@ sino3D_log = log(dose./max(dataRaw,1))*multifactor; %log corrected data -> sinog
clear dataExp sino_tomophan3D
%
%%
+%-------------Astra toolbox------------%
+% one can generate data using ASTRA toolbox
+proj_geom = astra_create_proj_geom('parallel', 1, det_size, angles_rad);
+vol_geom = astra_create_vol_geom(N,N);
+sino_ASTRA3D = zeros(det_size, length(angles), N, 'single');
+for i = 1:N
+[sino_id, sinoT] = astra_create_sino_cuda(TomoPhantom(:,:,i), proj_geom, vol_geom);
+sino_ASTRA3D(:,:,i) = sinoT';
+astra_mex_data2d('delete', sino_id);
+end
+%--------------------------------------%
+%%
% using ASTRA-toolbox to set the projection geometry (parallel beam)
proj_geom = astra_create_proj_geom('parallel', 1, det_size, angles_rad);
vol_geom = astra_create_vol_geom(N,N);