Iterative clustering material decomposition aided by empirical spectral correction for photon counting detectors in micro-CT

Purpose Photon counting detectors offer promising advancements in computed tomography (CT) imaging by enabling the quantification and three-dimensional imaging of contrast agents and tissue types through simultaneous multi-energy projections from broad X-ray spectra. However, the accuracy of these decomposition methods hinges on precise composite spectral attenuation values that one must reconstruct from spectral micro-CT. Errors in such estimations could be due to effects such as beam hardening, object scatter, or detector sensor-related spectral distortions such as fluorescence. Even if accurate spectral correction is done, multi-material separation within a volume remains a challenge. Increasing the number of energy bins in material decomposition problems often comes with a significant noise penalty but with minimal decomposition benefits. Approach We begin with an empirical spectral correction method executed in the tomographic domain that accounts for distortions in estimated spectral attenuation for each voxel. This is followed by our proposed iterative clustering material decomposition (ICMD) where clustering of voxels is used to reduce the number of basis materials to be resolved for each cluster. Using a larger number of energy bins for the clustering step shows distinct advantages in excellent classification to a larger number of clusters with accurate cluster centers when compared with the National Institute of Standards and Technology attenuation values. The decomposition step is applied to each cluster separately where each cluster has fewer basis materials compared with the entire volume. This is shown to reduce the need for the number of energy bins required in each decomposition step for the clusters. This approach significantly increases the total number of materials that can be decomposed within the volume with high accuracy and with excellent noise properties. Results Utilizing a (cadmium telluride 1-mm-thick sensor) Medipix detector with a 55-mu m pitch, we demonstrate the quantitatively accurate decomposition of several materials in a phantom study, where the sample includes mixtures of soft materials such as water and poly-methyl methacrylate along with contrast-enhancing materials. We show improved accuracy and lower noise when all five energy bins were used to yield effective classification of voxels into multiple accurate fundamental clusters which was followed by the decomposition step applied to each cluster using just two energy bins. We also show an example of biological sample imaging and separating three distinct types of tissue in mice: muscle, fat, and bone. Our experimental results show that the combination of effective and practical spectral correction and high-dimensional data clustering enhances decomposition accuracy and reduces noise in micro-CT. Conclusions This ICMD allows for quantitative separation of multiple materials including mixtures and also effectively separates multi-contrast agents.