An Algorithm for Constrained One-Step Inversion of Spectral CT Data
Institute of Physics
Physics in Medicine and Biology
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.
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.