1 function vol = astra_iradon_3d(sino, theta)
2 % this function calculates 3D iradon transform of the volume using ASTRA tools
6 % initialize ASTRA projector:
7 proj_geom = astra_create_proj_geom('parallel3d', 1, 1, sz(3), sz(1), theta);
8 vol_geom = astra_create_vol_geom(sz(1), sz(1), sz(3));
9 %proj_id = astra_create_projector('linear3d', proj_geom, vol_geom);
12 volume_id = astra_mex_data3d('create','-vol', vol_geom, 0);
13 sino_id = astra_mex_data3d('create','-sino', proj_geom, sino);
15 % initialize algorithm:
16 cfg = astra_struct('BP3D_CUDA');
17 %cfg.ProjectorId = proj_id;
18 cfg.ProjectionDataId = sino_id;
19 cfg.ReconstructionDataId = volume_id;
20 %cfg.FilterType = 'none';
23 alg_id = astra_mex_algorithm('create', cfg);
24 astra_mex_algorithm('run', alg_id);
25 vol = astra_mex_data3d('get', volume_id) * 2 * size(proj_geom.ProjectionAngles,2) / pi;
28 astra_mex_data3d('delete', sino_id);
29 astra_mex_data3d('delete', volume_id);
30 astra_mex_algorithm('delete', alg_id);
31 %astra_mex_projector('delete', proj_id);