diff options
Diffstat (limited to 'test_geometry.m')
-rw-r--r-- | test_geometry.m | 22 |
1 files changed, 15 insertions, 7 deletions
diff --git a/test_geometry.m b/test_geometry.m index f846d52..8952660 100644 --- a/test_geometry.m +++ b/test_geometry.m @@ -15,14 +15,23 @@ proj_geom_parallel = astra_create_proj_geom('parallel', 1, 64, 1:180); proj_geom_fanflat = astra_create_proj_geom('fanflat', 1, 64, 1:180, 0, 2000); proj_geom_fanflat_vec = astra_geom_2vec(proj_geom_fanflat); +proj_geom_parallel3d = astra_create_proj_geom('parallel3d', 1, 1, 64, 64, 1:180); +proj_geom_cone = astra_create_proj_geom('cone', 1, 1, 64, 64, 1:180, 0, 2000); +proj_geom_parallel3d_vec = astra_geom_2vec(proj_geom_parallel3d); +proj_geom_cone_vec = astra_geom_2vec(proj_geom_cone); + %% create data objects vol2d_id = astra_mex_data2d('create', '-vol', vol_geom2d, 0); vol3d_id = astra_mex_data3d('create', '-vol', vol_geom3d, 0); proj_parallel_id = astra_mex_data2d('create', '-sino', proj_geom_parallel, 0); proj_fanflat_id = astra_mex_data2d('create', '-sino', proj_geom_fanflat, 0); -proj_fanflatvec_id = astra_mex_data2d('create', '-sino', proj_geom_fanflat_vec, 0); +proj_fanflat_vec_id = astra_mex_data2d('create', '-sino', proj_geom_fanflat_vec, 0); +proj_parallel3d_id = astra_mex_data3d('create', '-sino', proj_geom_parallel3d, 0); +proj_cone_id = astra_mex_data3d('create', '-sino', proj_geom_cone, 0); +proj_parallel3d_vec_id = astra_mex_data3d('create', '-sino', proj_geom_parallel3d_vec, 0); +proj_cone_vec_id = astra_mex_data3d('create', '-sino', proj_geom_cone_vec, 0); %% get geometries vol_geom2d_new = astra_mex_data2d('get_geometry', vol2d_id); @@ -30,13 +39,12 @@ vol_geom3d_new = astra_mex_data3d('get_geometry', vol3d_id); proj_geom_parallel_new = astra_mex_data2d('get_geometry', proj_parallel_id); proj_geom_fanflat_new = astra_mex_data2d('get_geometry', proj_fanflat_id); -proj_geom_fanflat_vec_new = astra_mex_data2d('get_geometry', proj_fanflatvec_id); - -proj_geom_fanflat_vec -proj_geom_fanflat_vec_new +proj_geom_fanflat_vec_new = astra_mex_data2d('get_geometry', proj_fanflat_vec_id); -proj_geom_fanflat_vec.Vectors(110,:) -proj_geom_fanflat_vec_new.Vectors(110,:) +proj_geom_parallel3d_new = astra_mex_data3d('get_geometry', proj_parallel3d_id); +proj_geom_cone_new = astra_mex_data3d('get_geometry', proj_cone_id); +proj_geom_parallel3d_vec_new = astra_mex_data3d('get_geometry', proj_parallel3d_vec_id); +proj_geom_cone_vec_new = astra_mex_data3d('get_geometry', proj_cone_vec_id); %% astra_clear;
\ No newline at end of file |