Accurate Multi-Material Decomposition in Dual-Energy CT: A Phantom Study

• Published in: IEEE transactions on computational imaging Vol. 5; no. 4; pp. 515 - 529
• Main Authors: Xue, Yi, Sun, Xiaonan, Hu, Xiuhua, Sheng, Ke, Niu, Tianye, Jiang, Yangkang, Yang, Chunlin, Lyu, Qihui, Wang, Jing, Luo, Chen, Zhang, Luhan, Desrosiers, Catherine, Feng, Kun
• Format: Journal Article
• Language: English
• Published: Piscataway IEEE 01.12.2019; The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
• Subjects: Algorithms; Attenuation coefficients; block matching; Computed tomography; Correlation analysis; Decomposition; Dual-energy CT; Energy spectra; Euclidean distance; Euclidean geometry; Image reconstruction; Imaging phantoms; material triplets; Medical imaging; multi-material decomposition; Noise; Pixels; X-ray imaging
• ISSN: 2573-0436; 2333-9403
• DOI: 10.1109/TCI.2019.2909192

DUAL-energy computed tomography (DECT) differentiates materials by exploiting the varying material linear attenuation coefficients (LACs) for different x-ray energy spectra. Multi-material decomposition (MMD) is a particularly attractive DECT clinical application to distinguish the complicated material components within the human body. One prior material assisted (PMA) image domain MMD method was implemented, but has suffered from inaccurate decomposition, magnified noise, and expensive computation. To suppress the noise, we implemented a statistical MMD (SMMD) algorithm, which applied the statistical weight to account for the noise variance in the DECT images. Its decomposition accuracy heavily relies on the initial value. In this paper, we propose a novel method to overcome these challenges. Based on the piecewise constant property of CT images with energy-dependent LAC, we assume that the pixels with high similarity have the same material composition. We cluster pixel patches into groups using the block-matching technique. The material composition in each group is preselected according to the shortest Euclidean distance in the energy map between the center of mass of the similar patch groups and the LAC of the object with known material composition pre-assigned by the clinician. MMD is performed on the central pixel of each patch using the preselected material composition. In a preliminary study, the proposed method is evaluated using the digital and water phantoms. The proposed method increases the volume fraction by 25.2% and decreases the standard deviation by 66.2% compared with the PMA method and increases the volume fraction by 19.6% compared with the SMMD method. The proposed method achieves an overall improvement of the normalized cross-correlation matrix diagonality by 34.8% and 69.4% compared with the PMA and SMMD methods. The phantom results indicate that the proposed method has the potential to be applied to clinical practice due to its increased decomposition accuracy, and suppressed noise and cross contamination.