Maximum-Likelihood Dual-Energy TomographicImage Reconstruction

Dual-energy (DE) X-ray computed tomography (CT) has shown promise for material characterization and for providing quantitatively accurate CT values in a variety of applications. However, DE-CT has not been used routinely in medicine to date, primarily due to dose considerations. Most methods for DE-CT have used the filtered backprojection method for image reconstruction, leading to suboptimal noise/dose properties. This paper describes a statistical (maximum-likelihood) method for dual-energy X-ray CT that accommodates a wide variety of potential system configurations and measurement noise models. Regularized methods (such as penalized-likelihood or Bayesian estimation) are straightforward extensions. One version of the algorithm monotonically decreases the negative log-likelihood cost function each iteration. An ordered-subsets variation of the algorithm provides a fast and practical version.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85934/1/Fessler172.pd

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.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85895/1/Fessler74.pd

Fast kVp-Switching Dual Energy CT for PET Attenuation Correction

X-ray CT images are used routinely for attenuation correction in PET/CT systems. However, conventional CT-based attenuation correction (CTAC) can be inaccurate in regions containing iodine contrast agent. Dual-energy (DE) CT has the potential to improve the accuracy of attenuation correction in PET, but conventional DECT can suffer from motion artifacts. Recent X-ray CT systems can collect DE sinograms by rapidly switching the X-ray tube voltage between two levels for alternate projection views, reducing motion artifacts. The goal of this work is to study statistical methods for image reconstruction from both fast kVp-switching DE scans and from conventional dual-rotate DE scans in the context of CTAC for PET.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/86003/1/Fessler244.pd

Statistical X-Ray-Computed Tomography Image Reconstruction with Beam- Hardening Correction

This paper describes two statistical iterative reconstruction methods for X-ray CT. The rst method assumes a mono-energetic model for X-ray attenuation. We approximate the transmission Poisson likelihood by a quadratic cost function and exploit its convexity to derive a separable quadratic surrogate function that is easily minimized using parallelizable algorithms. Ordered subsets are used to accelerate convergence. We apply this mono-energetic algorithm (with edge-preserving regularization) to simulated thorax X-ray CT scans. A few iterations produce reconstructed images with lower noise than conventional FBP images at equivalent resolutions. The second method generalizes the physical model and accounts for the poly-energetic X-ray source spectrum and the measurement nonlinearities caused by energy-dependent attenuation. We assume the object consists of a given number of nonoverlapping tissue types. The attenuation coeÆcient of each tissue is the product of its unknown density and a known energy-dependent mass attenuation coeÆcient. We formulate a penalized-likelihood function for this polyenergetic model and develop an iterative algorithm for estimating the unknown densities in each voxel. Applying this method to simulated X-ray CT measurements of a phantom containing both bone and soft tissue yields images with signi cantly reduced beam hardening artifacts.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85939/1/Fessler165.pd