Statistical image reconstruction for polyenergetic X-ray computed tomography

This paper describes a statistical image reconstruction method for X-ray computed tomography (CT) that is based on a physical model that accounts for the polyenergetic X-ray source spectrum and the measurement nonlinearities caused by energy-dependent attenuation. We assume that the object consists of a given number of nonoverlapping materials, such as soft tissue and bone. The attenuation coefficient of each voxel is the product of its unknown density and a known energy-dependent mass attenuation coefficient. We formulate a penalized-likelihood function for this polyenergetic model and develop an ordered-subsets iterative algorithm for estimating the unknown densities in each voxel. The algorithm monotonically decreases the cost function at each iteration when one subset is used. Applying this method to simulated X-ray CT measurements of objects containing both bone and soft tissue yields images with significantly reduced beam hardening artifacts.

[1]  Jeffrey A. Fessler,et al.  Fast Monotonic Algorithms for Transmission Tomography , 1999, IEEE Trans. Medical Imaging.

[2]  P. Joseph,et al.  A Method for Correcting Bone Induced Artifacts in Computed Tomography Scanners , 1978, Journal of computer assisted tomography.

[3]  Wang,et al.  Iterative X-ray Cone-Beam Tomography for Metal Artifact Reduction and Local Region Reconstruction , 1999, Microscopy and Microanalysis.

[4]  O Nalcioglu,et al.  Post-reconstruction method for beam hardening in computerised tomography. , 1979, Physics in medicine and biology.

[5]  J. H. Hubbell,et al.  Tables of X-Ray Mass Attenuation Coefficients and Mass Energy-Absorption Coefficients 1 keV to 20 MeV for Elements Z = 1 to 92 and 48 Additional Substances of Dosimetric Interest , 1995 .

[6]  Harold L. Kundel,et al.  Handbook of Medical Imaging, Volume 1. Physics and Psychophysics , 2000 .

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

[8]  Neal H. Clinthorne,et al.  Design of an experimental system for dual energy X-ray CT , 1999, 1999 IEEE Nuclear Science Symposium. Conference Record. 1999 Nuclear Science Symposium and Medical Imaging Conference (Cat. No.99CH37019).

[9]  Hakan Erdogan,et al.  Ordered subsets algorithms for transmission tomography. , 1999, Physics in medicine and biology.

[10]  P M Joseph,et al.  A method for simultaneous correction of spectrum hardening artifacts in CT images containing both bone and iodine. , 1997, Medical physics.

[11]  Patrick Dupont,et al.  Reduction of metal streak artifacts in X-ray computed tomography using a transmission maximum a posteriori algorithm , 1999 .

[12]  Patrick Dupont,et al.  Reduction of metal artifacts in X-ray computed tomography , 2000 .

[13]  M. Abidi,et al.  Reconstruction of multi-energy X-ray computed tomography images of laboratory mice , 1999 .

[14]  幹男 山本.  “IEEE Trans. on Medical Imaging”誌と本誌 , 1983 .

[15]  Yoram Bresler,et al.  Globally convergent edge-preserving regularized reconstruction: an application to limited-angle tomography , 1998, IEEE Trans. Image Process..

[16]  Jeffrey A. Fessler,et al.  Ieee Transactions on Image Processing: to Appear Hybrid Poisson/polynomial Objective Functions for Tomographic Image Reconstruction from Transmission Scans , 2022 .

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

[18]  M. Defrise,et al.  Iterative reconstruction for helical CT: a simulation study. , 1998, Physics in medicine and biology.

[19]  A. Macovski,et al.  Energy-selective reconstructions in X-ray computerised tomography , 1976, Physics in medicine and biology.

[20]  Ken D. Sauer,et al.  A unified approach to statistical tomography using coordinate descent optimization , 1996, IEEE Trans. Image Process..

[21]  S Napel,et al.  Modeling of polychromatic attenuation using computed tomography reconstructed images. , 1999, Medical physics.

[22]  Isaac N. Bankman,et al.  Handbook of medical imaging , 2000 .

[23]  Jeffrey A. Fessler,et al.  Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction , 1997, IEEE Transactions on Medical Imaging.

[24]  H. Skinner,et al.  CT image correction for beam hardening using simulated projection data , 1990 .

[25]  K. Lange,et al.  EM reconstruction algorithms for emission and transmission tomography. , 1984, Journal of computer assisted tomography.

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

[27]  Ronald J. Jaszczak,et al.  Bayesian reconstruction and use of anatomical a priori information for emission tomography , 1996, IEEE Trans. Medical Imaging.

[28]  Chye Hwang Yan,et al.  Reconstruction algorithm for polychromatic CT imaging: application to beam hardening correction , 2000, IEEE Transactions on Medical Imaging.

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

[30]  R. Brooks,et al.  Beam hardening in X-ray reconstructive tomography , 1976 .

[31]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[32]  S. Manglos,et al.  Transmission maximum-likelihood reconstruction with ordered subsets for cone beam CT. , 1995, Physics in medicine and biology.

[33]  Alvaro R. De Pierro,et al.  A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography , 1995, IEEE Trans. Medical Imaging.

[34]  Jeffrey A. Fessler,et al.  Experimental evaluation for joint estimation approach , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[35]  K. Lange,et al.  A Theoretical Study of Some Maximum Likelihood Algorithms for Emission and Transmission Tomography , 1987, IEEE Transactions on Medical Imaging.

[36]  John M. Boone,et al.  X-Ray Production, Interaction, and Detection in Diagnostic Imaging , 2000 .

[37]  Hakan Erdogan,et al.  Monotonic algorithms for transmission tomography , 1999, IEEE Transactions on Medical Imaging.

[38]  Jolyon A. Browne,et al.  Developments with maximum likelihood X-ray computed tomography , 1992, IEEE Trans. Medical Imaging.

[39]  Predrag Sukovic,et al.  Penalized weighted least-squares image reconstruction for dual energy X-ray transmission tomography , 2000, IEEE Transactions on Medical Imaging.

[40]  Ken D. Sauer,et al.  A local update strategy for iterative reconstruction from projections , 1993, IEEE Trans. Signal Process..

[41]  J. Hsieh,et al.  An iterative approach to the beam hardening correction in cone beam CT. , 2000, Medical physics.

[42]  Patrick Dupont,et al.  An iterative maximum-likelihood polychromatic algorithm for CT , 2001, IEEE Transactions on Medical Imaging.

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

[44]  A. R. De Pierro,et al.  On the relation between the ISRA and the EM algorithm for positron emission tomography , 1993, IEEE Trans. Medical Imaging.

[45]  C. Helstrom,et al.  Compensation for readout noise in CCD images , 1995 .

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