1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
|
function [recon_id, recon] = astra_create_reconstruction(rec_type, proj_id, sinogram, iterations, use_mask, mask, use_minc, minc, use_maxc, maxc)
%--------------------------------------------------------------------------
% [recon_id, recon] = astra_create_reconstruction(rec_type, proj_id, sinogram, iterations, use_mask, mask, use_minc, minc, use_maxc, maxc)
%
% Create a CPU based iterative reconstruction.
%
% rec_type: reconstruction type, 'ART', 'SART' 'SIRT' or 'CGLS', not all options are adjustable
% proj_id: identifier of the projector as it is stored in the astra-library
% sinogram: sinogram data OR sinogram identifier
% iterations: number of iterations to perform
% use_mask: use a reconstrucionmask? 'yes' or 'no'
% mask: mask data OR mask identifier.
% use_minc: use a minimum constraint? 'yes' or 'no'
% minc: minimum constraint value
% use_maxc: use a maximum constraint? 'yes' or 'no'
% maxc: maximum constraint value
% recon_id: identifier of the reconstruction data object as it is now stored in the astra-library
% recon: MATLAB data version of the reconstruction
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
% This file is part of the ASTRA Toolbox
%
% Copyright: 2010-2016, iMinds-Vision Lab, University of Antwerp
% 2014-2016, CWI, Amsterdam
% License: Open Source under GPLv3
% Contact: astra@uantwerpen.be
% Website: http://www.astra-toolbox.com/
%--------------------------------------------------------------------------
if nargin <= 4
use_mask = 'no';
mask = [];
use_minc = 'no';
minc = 0;
use_maxc = 'no';
maxc = 255;
end
if nargin <= 6
use_minc = 'no';
minc = 0;
use_maxc = 'no';
maxc = 255;
end
if numel(sinogram) == 1
sinogram_id = sinogram;
else
proj_geom = astra_mex_projector('projection_geometry', proj_id);
sinogram_id = astra_mex_data2d('create', '-sino', proj_geom, sinogram);
end
% create reconstruction object
vol_geom = astra_mex_projector('volume_geometry', proj_id);
recon_id = astra_mex_data2d('create', '-vol', vol_geom, 0);
% configure
cfg = astra_struct(rec_type);
cfg.ProjectorId = proj_id;
cfg.ProjectionDataId = sinogram_id;
cfg.ReconstructionDataId = recon_id;
if strcmp(use_mask,'yes')
if numel(mask) == 1
mask_id = mask;
else
mask_id = astra_mex_data2d('create', '-vol', vol_geom, mask);
end
cfg.options.ReconstructionMaskId = mask_id;
end
cfg.options.UseMinConstraint = use_minc;
cfg.options.MinConstraintValue = minc;
cfg.options.UseMaxConstraint = use_maxc;
cfg.options.MaxConstraintValue = maxc;
cfg.options.ProjectionOrder = 'random';
alg_id = astra_mex_algorithm('create', cfg);
% iterate
astra_mex_algorithm('iterate', alg_id, iterations);
% return object
recon = astra_mex_data2d('get', recon_id);
% garbage collection
astra_mex_algorithm('delete', alg_id);
if numel(sinogram) ~= 1
astra_mex_data2d('delete', sinogram_id);
end
if strcmp(use_mask,'yes')
if numel(mask) ~= 1
astra_mex_data2d('delete', mask_id);
end
end
|
