An Algorithm for Constrained One-Step Inversion of Spectral CT Data

Authors

Rina Foygel Barber, University of Chicago
Emil Y. Sidky, University of Chicago
Taly Gilat Schmidt, Marquette UniversityFollow
Xiaochuan Pan, University of ChicagoFollow

Document Type

Article

Language

Eng

Publication Date

4-19-2016

Publisher

Institute of Physics

Source Publication

Physics in Medicine and Biology

Source ISSN

1361-6560

Abstract

We develop a primal-dual algorithm that allows for one-step inversion of spectral CT transmission photon counts data to a basis map decomposition. The algorithm allows for image constraints to be enforced on the basis maps during the inversion. The derivation of the algorithm makes use of a local upper bounding quadratic approximation to generate descent steps for non-convex spectral CT data discrepancy terms, combined with a new convex-concave optimization algorithm. Convergence of the algorithm is demonstrated on simulated spectral CT data. Simulations with noise and anthropomorphic phantoms show examples of how to employ the constrained one-step algorithm for spectral CT data.

Comments

Published version. Physics in Medicine and Biology, Vol. 61, No. 10 (2016): 3784-3818. DOI. © 2016 Institute of Physics and Engineering in Medicine. Used with permission.

Creative Commons License

This work is licensed under a Creative Commons Attribution 3.0 License.

Recommended Citation

Barber, Rina Foygel; Sidky, Emil Y.; Schmidt, Taly Gilat; and Pan, Xiaochuan, "An Algorithm for Constrained One-Step Inversion of Spectral CT Data" (2016). Biomedical Engineering Faculty Research and Publications. 526.
https://epublications.marquette.edu/bioengin_fac/526