A Bayesian-MAP Method Based on TV for CT Image Reconstruction from Sparse and Limited Data

Computed tomography (CT) plays an important role in the field of modern medical imaging. Reducing radiation exposure dose without significantly decreasing image’s quality is always a crucial issue. Inspired by the outstanding performance of total variation (TV) technique in CT image reconstruction, a TV regularization based Bayesian-MAP (MAP-TV) is proposed to reconstruct the case of sparse view projection and limited angle range imaging. This method can suppress the streak artifacts and geometrical deformation while preserving image edges. We use ordered subset (OS) technique to accelerate the reconstruction speed. Numerical results showed that MAP-TV is able to reconstruct a phantom with better visual performance and quantitative evaluation than classical FBP, MLEM and quadrate prior MAP algorithms. The proposed algorithm can be generalized to conebeam CT image reconstruction.

[1]  E. Sidky,et al.  Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization , 2008, Physics in medicine and biology.

[2]  P.H. Pretorius,et al.  Evaluation of the ordered-subset transmission (OSTR) algorithm for transmission imaging on SPECT systems using axially overlapping cone-beams , 2006, IEEE Transactions on Nuclear Science.

[3]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[4]  Yang-Ming Zhu Extrapolation method to obtain a full projection from a limited one for limited-view tomography problems , 2006, Medical and Biological Engineering and Computing.

[5]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[6]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .

[7]  E. Sidky,et al.  Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT , 2009, 0904.4495.

[8]  P. Green Bayesian reconstructions from emission tomography data using a modified EM algorithm. , 1990, IEEE transactions on medical imaging.

[9]  R. Siddon Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[10]  H. Stark,et al.  A comparative study of three reconstruction methods for a limited-view computer tomography problem. , 1989, IEEE transactions on medical imaging.

[11]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

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

[13]  Steve B. Jiang,et al.  GPU-based Fast Low Dose Cone Beam CT Reconstruction via Total Variation , 2010, 1001.0599.