Taly Gilat Schmidt1, Barbara A Sammut1, Rina Foygel Barber2, Xiaochuan Pan3, Emil Y Sidky3. 1. Department of Biomedical Engineering, Marquette University and Medical College of Wisconsin, Milwaukee, Wisconsin, USA. 2. Department of Statistics, University of Chicago, Chicago, Illinois, USA. 3. Department of Radiology, University of Chicago, Chicago, Illinois, USA.
Abstract
PURPOSE: The constrained one-step spectral CT image reconstruction (cOSSCIR) algorithm with a nonconvex alternating direction method of multipliers optimizer is proposed for addressing computed tomography (CT) metal artifacts caused by beam hardening, noise, and photon starvation. The quantitative performance of cOSSCIR is investigated through a series of photon-counting CT simulations. METHODS: cOSSCIR directly estimates basis material maps from photon-counting data using a physics-based forward model that accounts for beam hardening. The cOSSCIR optimization framework places constraints on the basis maps, which we hypothesize will stabilize the decomposition and reduce streaks caused by noise and photon starvation. Another advantage of cOSSCIR is that the spectral data need not be registered, so that a ray can be used even if some energy window measurements are unavailable. Photon-counting CT acquisitions of a virtual pelvic phantom with low-contrast soft tissue texture and bilateral hip prostheses were simulated. Bone and water basis maps were estimated using the cOSSCIR algorithm and combined to form a virtual monoenergetic image for the evaluation of metal artifacts. The cOSSCIR images were compared to a "two-step" decomposition approach that first estimated basis sinograms using a maximum likelihood algorithm and then reconstructed basis maps using an iterative total variation constrained least-squares optimization (MLE+TV min $_{\text{min}}$ ). Images were also compared to a nonspectral TV min $_{\text{min}}$ reconstruction of the total number of counts detected for each ray with and without normalized metal artifact reduction (NMAR) applied. The simulated metal density was increased to investigate the effects of increasing photon starvation. The quantitative error and standard deviation in regions of the phantom were compared across the investigated algorithms. The ability of cOSSCIR to reproduce the soft-tissue texture, while reducing metal artifacts, was quantitatively evaluated. RESULTS: Noiseless simulations demonstrated the convergence of the cOSSCIR and MLE+TV min $_{\text{min}}$ algorithms to the correct basis maps in the presence of beam-hardening effects. When noise was simulated, cOSSCIR demonstrated a quantitative error of -1 HU, compared to 2 HU error for the MLE+TV min $_{\text{min}}$ algorithm and -154 HU error for the nonspectral TV min $_{\text{min}}$ +NMAR algorithm. For the cOSSCIR algorithm, the standard deviation in the central iodine region of interest was 20 HU, compared to 299 HU for the MLE+TV min $_{\text{min}}$ algorithm, 41 HU for the MLE+TV min $_{\text{min}}$ +Mask algorithm that excluded rays through metal, and 55 HU for the nonspectral TV min $_{\text{min}}$ +NMAR algorithm. Increasing levels of photon starvation did not impact the bias or standard deviation of the cOSSCIR images. cOSSCIR was able to reproduce the soft-tissue texture when an appropriate regularization constraint value was selected. CONCLUSIONS: By directly inverting photon-counting CT data into basis maps using an accurate physics-based forward model and a constrained optimization algorithm, cOSSCIR avoids metal artifacts due to beam hardening, noise, and photon starvation. The cOSSCIR algorithm demonstrated improved stability and accuracy compared to a two-step method of decomposition followed by reconstruction.
PURPOSE: The constrained one-step spectral CT image reconstruction (cOSSCIR) algorithm with a nonconvex alternating direction method of multipliers optimizer is proposed for addressing computed tomography (CT) metal artifacts caused by beam hardening, noise, and photon starvation. The quantitative performance of cOSSCIR is investigated through a series of photon-counting CT simulations. METHODS: cOSSCIR directly estimates basis material maps from photon-counting data using a physics-based forward model that accounts for beam hardening. The cOSSCIR optimization framework places constraints on the basis maps, which we hypothesize will stabilize the decomposition and reduce streaks caused by noise and photon starvation. Another advantage of cOSSCIR is that the spectral data need not be registered, so that a ray can be used even if some energy window measurements are unavailable. Photon-counting CT acquisitions of a virtual pelvic phantom with low-contrast soft tissue texture and bilateral hip prostheses were simulated. Bone and water basis maps were estimated using the cOSSCIR algorithm and combined to form a virtual monoenergetic image for the evaluation of metal artifacts. The cOSSCIR images were compared to a "two-step" decomposition approach that first estimated basis sinograms using a maximum likelihood algorithm and then reconstructed basis maps using an iterative total variation constrained least-squares optimization (MLE+TV min $_{\text{min}}$ ). Images were also compared to a nonspectral TV min $_{\text{min}}$ reconstruction of the total number of counts detected for each ray with and without normalized metal artifact reduction (NMAR) applied. The simulated metal density was increased to investigate the effects of increasing photon starvation. The quantitative error and standard deviation in regions of the phantom were compared across the investigated algorithms. The ability of cOSSCIR to reproduce the soft-tissue texture, while reducing metal artifacts, was quantitatively evaluated. RESULTS: Noiseless simulations demonstrated the convergence of the cOSSCIR and MLE+TV min $_{\text{min}}$ algorithms to the correct basis maps in the presence of beam-hardening effects. When noise was simulated, cOSSCIR demonstrated a quantitative error of -1 HU, compared to 2 HU error for the MLE+TV min $_{\text{min}}$ algorithm and -154 HU error for the nonspectral TV min $_{\text{min}}$ +NMAR algorithm. For the cOSSCIR algorithm, the standard deviation in the central iodine region of interest was 20 HU, compared to 299 HU for the MLE+TV min $_{\text{min}}$ algorithm, 41 HU for the MLE+TV min $_{\text{min}}$ +Mask algorithm that excluded rays through metal, and 55 HU for the nonspectral TV min $_{\text{min}}$ +NMAR algorithm. Increasing levels of photon starvation did not impact the bias or standard deviation of the cOSSCIR images. cOSSCIR was able to reproduce the soft-tissue texture when an appropriate regularization constraint value was selected. CONCLUSIONS: By directly inverting photon-counting CT data into basis maps using an accurate physics-based forward model and a constrained optimization algorithm, cOSSCIR avoids metal artifacts due to beam hardening, noise, and photon starvation. The cOSSCIR algorithm demonstrated improved stability and accuracy compared to a two-step method of decomposition followed by reconstruction.
Authors: Ruoqiao Zhang; Jean-Baptiste Thibault; Charles A Bouman; Ken D Sauer; Jiang Hsieh Journal: IEEE Trans Med Imaging Date: 2013-09-17 Impact factor: 10.048
Authors: Jeffrey F Williamson; Bruce R Whiting; Jasenka Benac; Ryan J Murphy; G James Blaine; Joseph A O'Sullivan; David G Politte; Donald L Snyder Journal: Med Phys Date: 2002-10 Impact factor: 4.071
Authors: Seung Chol Han; Yong Eun Chung; Young Han Lee; Kwan Kyu Park; Myeong Jin Kim; Ki Whang Kim Journal: AJR Am J Roentgenol Date: 2014-10 Impact factor: 3.959
Authors: Nicolas Ducros; Juan Felipe Perez-Juste Abascal; Bruno Sixou; Simon Rit; Françoise Peyrin Journal: Med Phys Date: 2017-09 Impact factor: 4.071