Jun 18
An improved parallel projection image reconstructor for Matlab
This post describes a replacement for iradon for Matlab and Octave. The new function, CTrecon.m,
1. Fixes a bug in iradon that results in offsets from the true values in the reconstructed image.
2. Uses a mex function from iradon_speedy for backprojection. According to a comment by Jeff Orchard, the author of iradon_speedy, on the Mathworks File Exchange entry for iradon_speedy, Matlab R2011a now uses a mex function. I do not have that version so I cannot confirm this.
3. Handles complex projections to return a complex reconstructed image – see zBackproject.c
4. Regularizes the interface to (parameter-name’,parameter-value) syntax. This implies that CTrecon cannot just be dropped in for iradon. The advantage is more readable code.
Comments Off