Dual-energy CT material decomposition in the presence of nonuniqueness

Purpose: We develop a pathlength recovery method for material decomposition to address the challenge of nonunique mappings from fluence transmission to basis material pathlengths in dual-energy computed tomography (DECT) and implement it into a two-step algorithm to reconstruct material basis maps and virtual monochromatic images (VMI). Methods: By formulating the DECT imaging model as a system of implicit equations to form level curves, we present a projection-based decomposition method that explicitly accounts for cases of fluence transmissions mapping to multiple material pathlength combinations. Using numerical simulations, we investigate the method's capability to recover and analyze multiple pathlength solutions, resulting in improved material basis map reconstruction and VMI synthesis. Results: Examples of nonunique fluence mappings are demonstrated using ploychromatic spectra (70 and 120 kV), demonstrating the potential for degraded material basis maps and VMI accuracy. We apply the pathlength recovery method to simulated DECT measurements of a digital phantom and demonstrate artifacts introduced by nonunique mappings if not properly resolved. Significance: This work develops a method for resolving nonunique fluence transmission to material pathlength mappings in DECT material decomposition while furthing the theoretical understanding of DECT system behaviors. Our method address efforts to enhance reliability and accuracy of material decomposition that has implications relevant to general spectral CT imaging.