A general CT reconstruction algorithm for model-based material decomposition

Material decomposition in CT has the potential to reduce artifacts and improve quantitative accuracy by utilizing spectral models and multi-energy scans. In this work we present a novel Model-Based Material Decomposition (MBMD) method based on an existing iterative reconstruction algorithm derived from a general non-linear forward model. A digital water phantom with inserts containing different concentrations of calcium was scanned on a kV switching system. We used the presented method to simultaneously reconstruct water and calcium material density images, and compared the results to an image domain and a projection domain decomposition method. When switching voltage every other frame, MBMD resulted in more accurate water and calcium concentration values than the image domain decomposition method, and was just as accurate as the projection domain decomposition method. In a second, slower, kV switching scheme (changing voltage every ten frames) which precluded the use of traditional projection domain based methods, MBMD continued to produce quantitatively accurate reconstructions. Finally, we present a preliminary study applying MBMD to a water phantom containing vials of different concentrations of K₂HPO₄ which was scanned on a cone-beam CT test bench. Both the fast and slow (emulated) kV switching schemes resulted in similar reconstructions, indicating MBMD's robustness to challenging acquisition schemes. Additionally, the K2HPO4 concentration ratios between the vials were accurately represented in the reconstructed K₂HPO₄ density image.