Noise and resolution in images reconstructed with FBP and OSC algorithms for CT.

This paper presents a comparison between an analytical and a statistical iterative reconstruction algorithm for computed transmission tomography concerning their noise and resolution performance. The reconstruction of two-dimensional images from simulated fan-beam transmission data is performed with a filtered back-projection (FBP) type reconstruction and an iterative ordered subsets convex (OSC) maximum-likelihood method. A special software phantom, which allows measuring the resolution and noise in a nonambiguous way, is used to simulate transmission tomography scans with different signal-to-noise ratios (SNR). The noise and modulation transfer function is calculated for FBP and OSC reconstruction at several positions, distributed over the field-of-view (FOV). The reconstruction with OSC using different numbers of subsets shows an inverse linear relation to the number of iterations that are necessary to reach a certain resolution and SNR, i.e., increasing the number of subsets by a factor x reduces the number of required iterations by the same factor. The OSC algorithm is able to achieve a nearly homogeneous high resolution over the whole FOV, which is not achieved with FBP. The OSC method achieves a lower level of noise compared with FBP at the same resolution. The reconstruction with OSC can save a factor of up to nine of x-ray dose compared with FBP in the investigated range of noise levels.

[1]  R. Proksa,et al.  A quasiexact reconstruction algorithm for helical CT using a 3-Pi acquisition. , 2003, Medical physics.

[2]  Charles L. Byrne,et al.  Iterative image reconstruction algorithms based on cross-entropy minimization , 1993, IEEE Trans. Image Process..

[3]  Freek J. Beekman,et al.  Accelerated iterative transmission CT reconstruction using an ordered subsets convex algorithm , 1998, IEEE Transactions on Medical Imaging.

[4]  F. Noo,et al.  Cone-beam reconstruction using 1D filtering along the projection of M-lines , 2005 .

[5]  R. Lewitt.  Alternatives to voxels for image representation in iterative reconstruction algorithms , 1992, Physics in medicine and biology.

[6]  Thomas Köhler,et al.  Efficient projection and backprojection scheme for spherically symmetric basis functions in divergent beam geometry. , 2006, Medical physics.

[7]  Michael A. King,et al.  An interior point iterative maximum-likelihood reconstruction algorithm incorporating upper and lower bounds with application to SPECT transmission imaging , 2001, IEEE Transactions on Medical Imaging.

[8]  B. De Man,et al.  Distance-driven projection and backprojection in three dimensions. , 2004, Physics in medicine and biology.

[9]  Robert M. Lewitt,et al.  Practical considerations for 3-D image reconstruction using spherically symmetric volume elements , 1996, IEEE Trans. Medical Imaging.

[10]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[11]  Thomas Köhler,et al.  EnPiT: filtered back-projection algorithm for helical CT using an n-Pi acquisition , 2005, IEEE Transactions on Medical Imaging.

[12]  Jeffrey A. Fessler,et al.  Ieee Transactions on Image Processing: to Appear Globally Convergent Algorithms for Maximum a Posteriori Transmission Tomography , 2022 .

[13]  R. Proksa,et al.  Systematic investigation of the reconstruction of images from transmission tomography using a filtered backprojection and an iterative OSML reconstruction algorithm , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

[14]  R M Lewitt,et al.  Multidimensional digital image representations using generalized Kaiser-Bessel window functions. , 1990, Journal of the Optical Society of America. A, Optics and image science.

[15]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[16]  Samuel Matej,et al.  2.5-D simultaneous multislice reconstruction by series expansion methods from Fourier-rebinned PET data , 2000, IEEE Transactions on Medical Imaging.

[17]  Alexander Katsevich,et al.  An improved exact filtered backprojection algorithm for spiral computed tomography , 2004, Adv. Appl. Math..

[18]  Lucas C. Parra,et al.  List-Mode Likelihood EM Algorithm and Noise Estimation Demonstrated on 2D-PET , 1998, IEEE Trans. Medical Imaging.

[19]  T. Peters,et al.  Computed Tomography with Fan Beam Geometry , 1977, Journal of computer assisted tomography.

[20]  J. Fessler.  Statistical Image Reconstruction Methods for Transmission Tomography , 2000 .

[21]  Jeffrey A. Fessler,et al.  Compensation for nonuniform resolution using penalized-likelihood reconstruction in space-variant imaging systems , 2004, IEEE Transactions on Medical Imaging.

[22]  Alvaro R. De Pierro,et al.  Fast EM-like methods for maximum "a posteriori" estimates in emission tomography , 2001, IEEE Transactions on Medical Imaging.