Convolutional neural network for material decomposition in spectral CT scans

Spectral computed tomography acquires energy-resolved data that allows recovery of densities of constituents of an object. This can be achieved by decomposing the measured spectral projection into material projections, and passing these decomposed projections through a tomographic reconstruction algorithm, to get the volumetric mass density of each material. Material decomposition is a nonlinear inverse problem that has been traditionally solved using model-based material decomposition algorithms. However, the forward model is difficult to estimate in real prototypes. Moreover, the traditional regularizers used to stabilized inversions are not fully relevant in the projection domain. In this study, we propose a deep-learning method for material decomposition in the projection domain. We validate our methodology with numerical phantoms of human knees that are created from synchrotron CT scans. We consider four different scans for training, and one for validation. The measurements are corrupted by Poisson noise, assuming that at most 10⁵ photons hit the detector. Compared to a regularized Gauss-Newton algorithm, the proposed deep-learning approach provides a compromise between noise and resolution, which reduces the computation time by a factor of 100.