Quadratic regularization design for fan beam transmission tomography

Statistical methods for tomographic image reconstruction have shown considerable potential for improving image quality in X-ray CT. Penalized-likelihood (PL) image reconstruction methods require maximizing an objective function that is based on the log-likelihood of the sinogram measurements and on a roughness penalty function to control noise. In transmission tomography, PL methods (and MAP methods) based on conventional quadratic regularization functions lead to nonuniform and anisotropic spatial resolution, even for idealized shift-invariant imaging systems. We have previously addressed this problem for parallel-beam emission tomography by designing data-dependent, shift-variant regularizers that improve resolution uniformity. This paper extends those methods to the fan-beam geometry used in X-ray CT imaging. Simulation results demonstrate that the new method for regularization design requires very modest computation and leads to nearly uniform and isotropic spatial resolution in the fan-beam geometry when using quadratic regularization.

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

[2]  Jeffrey A. Fessler,et al.  Regularization for uniform spatial resolution properties in penalized-likelihood image reconstruction , 2000, IEEE Transactions on Medical Imaging.

[3]  Carl R. Crawford,et al.  Reconstruction Algorithm for Fan Beam with a Displaced Center-of-Rotation , 1986, IEEE Transactions on Medical Imaging.

[4]  Richard M. Leahy,et al.  A theoretical study of the contrast recovery and variance of MAP reconstructions with applications to the selection of smoothing parameters , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[5]  J.A. Fessler.  Analytical approach to regularization design for isotropic spatial resolution , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[6]  Jeffrey A. Fessler,et al.  Conjugate-gradient preconditioning methods for shift-variant PET image reconstruction , 1999, IEEE Trans. Image Process..

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

[8]  Jeffrey A. Fessler.  Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, IEEE Trans. Medical Imaging.

[9]  Guy M. Besson,et al.  CT fan-beam parametrizations leading to shift-invariant filtering , 1996 .