Wavelet Regularized Alternating Minimization Algorithm for Low Dose X-ray CT

Jingwei Lu, Joseph A. O’Sullivan, David G. Politte

DOI:10.12059/Fully3D.2017-11-3201030

Published in:Fully3D 2017 Proceedings

Pages:726-732

Keywords:
computed tomography, low dose, wavelet, alternating minimization, dual domain optimization
X-ray computed tomography reconstruction has evolved over 40 years for medical, security, and industrial applications. Compared to traditional analytic reconstruction techniques such as filtered back projection (FBP), statistical reconstruction algorithms like alternating minimization (AM) provide improved image quality and can incorporate prior information. Increasing patient safety through reduced radiation dose results in fewer measured photons. Penalized AM is a powerful tool for maintaining image quality with less data, but the weight of penalty must be chosen carefully. If the penalty weight is too low, noise may not be suppressed and artifacts may be exhibited, such as those due to sharp discontinuities in attenuation at edges of dense material. If the penalty weight is higher, noise and artifacts may be reduced, but at the expense of introducing bias into the reconstruction. These contradicting requirements for the weight of the penalty limit our ability to improve reconstructed image quality in a low dose scenario. In this paper, we develop a new algorithm called wavelet regularized alternating minimization (wav-AM) by introducing a second penalty term on wavelet coefficients. By solving this dual domain optimization problem, we are able to perform 3D reconstruction of scanned baggage with low X-ray photon intensity. A medical imaging application of the wav-AM algorithm will be provided to illustrate the performance in image quality improvement. Evaluation of these real data reconstructions show reduced noise and artifacts without biasing the estimated attenuation of objects of known attenuation. The wav-AM algorithm features guaranteed convergence and increases the computational burden compared to the usual penalized AM algorithm only negligibly, even though we are solving a dual domain optimization problem.
Jingwei Lu
Washington university in St. Louis, USA
Joseph A. O’Sullivan
Washington university in St. Louis, USA
David G. Politte
Washington university School of Medicine, USA
  1. Q. Xu, D. Yang, J. Tan, A. Sawatzky, and M. A. Anastasio, “Accelerated fast iterative shrinkage thresholding algorithms for sparsity-regularized cone-beam CT image reconstruction,” Medical Physics, vol. 43, p. 1849, Apr. 2016.
  2. J. A. O’Sullivan, D. L. Snyder, and B. R. Whiting, “Alternating minimization algorithms for transmission tomography using energy detectors,” in Proc. Conf. Record of the Thirty-Sixth Asilomar Conf. Signals, Systems and Computers, vol. 1, pp. 144–147, vol.1, Nov. 2002.
  3. J. A. Fessler, E. P. Ficaro, N. H. Clinthorne, and K. Lange, “Groupedcoordinate ascent algorithms for penalized-likelihood transmission image reconstruction,” IEEE Transactions on Medical Imaging, vol. 16, no. 2, pp. 166–175, 1997.
  4. L. I. Rudin and S. Osher, “Total variation based image restoration with free local constraints,” in Proc. 1st Int. Conf. Image Processing, vol. 1, pp. 31–35, vol.1, Nov. 1994.
  5. K. Lange, R. Carson, et al., “EM reconstruction algorithms for emission and transmission tomography,” J Comput Assist Tomogr, vol. 8, no. 2, pp. 306–16, 1984.
  6. “American national standard for evaluating the image quality of x-ray computed tomography (CT) security-screening systems,” ANSI N42.45-2011, pp. 1–58, May 2011.
  7. S. Degirmenci, D. G. Politte, C. Bosch, N. Tricha, and J. A. O’Sullivan, “Acceleration of iterative image reconstruction for x-ray imaging for security applications,” in SPIE/IS&T Electronic Imaging, pp. 94010C–94010C, International Society for Optics and Photonics, 2015.
  8. S. J. Wright, R. D. Nowak, and M. A. T. Figueiredo, “Sparse reconstruction by separable approximation,” IEEE Transactions on Signal Processing, vol. 57, pp. 2479–2493, July 2009.
  9. A. Beck and M. Teboulle, “A fast iterative shrinkage-thresholding algorithm for linear inverse problems,” SIAM Journal on Imaging Sciences, vol. 2, no. 1, pp. 183–202, 2009.
  10. H. Erdogan and J. A. Fessler, “Ordered subsets algorithms for transmission tomography.,” Physics in Medicine and Biology, vol. 44, pp. 2835–2851, Nov. 1999.