Spline Driven: High Accuracy Projectors for Tomographic Reconstruction From Few Projections

• Published in: IEEE transactions on image processing Vol. 24; no. 12; pp. 4715 - 4725
• Main Authors: Momey, Fabien, Denis, Loic, Burnier, Catherine, Thiebaut, Eric, Becker, Jean-Marie, Desbat, Laurent
• Format: Journal Article
• Language: English
• Published: United States IEEE 01.12.2015; The Institute of Electrical and Electronics Engineers, Inc. (IEEE); Institute of Electrical and Electronics Engineers
• Subjects: Accuracy; Algorithms; Approximation methods; B-splines; Computational modeling; Computer Science; Finite element analysis; Image Processing, Computer-Assisted - methods; Image reconstruction; Inverse Problems; Kernel; Medical Imaging; Phantoms, Imaging; Projectors; Reconstruction; Signal and Image Processing; Signal processing; Signal Processing, Computer-Assisted; Splines (mathematics); Tomography; Tomography - methods; Tomography; reconstruction; signal processing; inverse problems; B-splines; Reconstruction; Signal processing; Inverse Problems
• ISSN: 1057-7149; 1941-0042
• DOI: 10.1109/TIP.2015.2466083

Tomographic iterative reconstruction methods need a very thorough modeling of data. This point becomes critical when the number of available projections is limited. At the core of this issue is the projector design, i.e., the numerical model relating the representation of the object of interest to the projections on the detector. Voxel driven and ray driven projection models are widely used for their short execution time in spite of their coarse approximations. Distance driven model has an improved accuracy but makes strong approximations to project voxel basis functions. Cubic voxel basis functions are anisotropic, accurately modeling their projection is, therefore, computationally expensive. Both smoother and more isotropic basis functions better represent the continuous functions and provide simpler projectors. These considerations have led to the development of spherically symmetric volume elements, called blobs. Set apart their isotropy, blobs are often considered too computationally expensive in practice. In this paper, we consider using separable B-splines as basis functions to represent the object, and we propose to approximate the projection of these basis functions by a 2D separable model. When the degree of the B-splines increases, their isotropy improves and projections can be computed regardless of their orientation. The degree and the sampling of the B-splines can be chosen according to a tradeoff between approximation quality and computational complexity. We quantitatively measure the good accuracy of our model and compare it with other projectors, such as the distance-driven and the model proposed by Long et al. From the numerical experiments, we demonstrate that our projector with an improved accuracy better preserves the quality of the reconstruction as the number of projections decreases. Our projector with cubic B-splines requires about twice as many operations as a model based on voxel basis functions. Higher accuracy projectors can be used to improve the resolution of the existing systems, or to reduce the number of projections required to reach a given resolution, potentially reducing the dose absorbed by the patient.