Quadratic regularization design for 3D axial CT: Towards isotropic noise

Potential advantages of statistical image reconstruction (SIR) methods over conventional filtered back-projection (FBP) method include reduced patient dose, improved noise and spatial resolution properties. However, the use of statistical weightings and sophisticated modeling of the system that are responsible for these improvements can lead to anisotropic and nonuniform noise characteristics and spatial resolutions. As an extension to our previous work that aimed for more isotropic and uniform spatial resolution, we propose a quadratic regularization design method for 3D axial X-ray computed tomography (CT) to achieve isotropy and uniformity of noise characteristics in the reconstructed image. In practice, noise properties may affect diagnostic quality of the reconstructed image as much as the spatial resolution. Simulations and a phantom experiment demonstrate that the proposed method leads to more isotropic and uniform noise characteristics in 3D axial CT with modest computational cost.

[1]  Richard M. Leahy,et al.  Resolution and noise properties of MAP reconstruction for fully 3-D PET , 2000, IEEE Transactions on Medical Imaging.

[2]  Jeffrey A. Fessler,et al.  Fast variance computation for quadratically penalized iterative reconstruction of 3D axial CT images , 2012, 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference Record (NSS/MIC).

[3]  Jeffrey A. Fessler,et al.  Spatial resolution properties of penalized-likelihood image reconstruction: space-invariant tomographs , 1996, IEEE Trans. Image Process..

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

[5]  Jeffrey A. Fessler Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): applications to tomography , 1996, IEEE Trans. Image Process..

[6]  Olaf Dössel,et al.  An iterative method for the reconstruction of the coronary arteries from rotational x-ray angiography , 2007, SPIE Medical Imaging.

[7]  Jean-Baptiste Thibault,et al.  A three-dimensional statistical approach to improved image quality for multislice helical CT. , 2007, Medical physics.

[8]  Jeffrey A. Fessler,et al.  Quadratic Regularization Design for 2-D CT , 2009, IEEE Transactions on Medical Imaging.

[9]  Charles L. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[10]  Jeffrey A. Fessler,et al.  Quadratic Regularization Design for 3 D Axial CT , 2013 .