An algorithm for constrained one-step inversion of spectral CT data
Medical Physics
2016-05-04 v1
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.
Cite
@article{arxiv.1511.03384,
title = {An algorithm for constrained one-step inversion of spectral CT data},
author = {Rina Foygel Barber and Emil Y. Sidky and Taly Gilat Schmidt and Xiaochuan Pan},
journal= {arXiv preprint arXiv:1511.03384},
year = {2016}
}
Comments
Submitted to Physics in Medicine and Biology