Abstract
Spectral computed tomography (CT) exploits the measurements obtained by a photon counting detector to reconstruct the chemical composition of an object. In particular, spectral CT has shown a very good ability to image K-edge contrast agent. Spectral CT is an inverse problem that can be addressed solving two subproblems, namely the basis material decomposition (BMD) problem and the tomographic reconstruction problem. In this work, we focus on the BMD problem, which is ill-posed and nonlinear. The BDM problem is classically either linearized, which enables reconstruction based on compressed sensing methods, or nonlinearly solved with no explicit regularization scheme. In a previous communication, we proposed a nonlinear regularized Gauss-Newton (GN) algorithm.1 However, this algorithm can only be applied to convex regularization functionals. In particular, the ep (p < 1) norm or the `0 quasi-norm, which are known to provider sparse solutions, cannot be considered. In order to better promote the sparsity of contrast agent images, we propose a nonlinear reconstruction framework that can handle nonconvex regularization terms. In particular, the e1/e2 norm ratio is considered.2 The problem is solved iteratively using the block variable metric forward-backward (BVMF-B) algorithm,3 which can also enforce the positivity of the material images. The proposed method is validated on numerical data simulated in a thorax phantom made of soft tissue, bone and gadolinium, which is scanned with a 90-kV x-ray tube and a 3-bin photon counting detector.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call