Quantitative multi-energy micro-CT: A simulation and phantom study for simultaneous imaging of four different contrast materials using an energy integrating detector

Emerging from the development of single-energy Computed Tomography (CT) and Dual-Energy Computed Tomography, Multi-Energy Computed Tomography (MECT) is a promising tool allowing advanced material and tissue decomposition and thereby enabling the use of multiple contrast materials in preclinical research.The scope of this work was to evaluate whether a usual preclinical micro-CT system is applicable for the decomposition of different materials using MECT together with a matrix-inversion method and how different changes of the measurement-environment affect the results.A matrix-inversion based algorithm to differentiate up to five materials (iodine, iron, barium, gadolinium, residual material) by applying four different acceleration voltages/energy levels was established. We carried out simulations using different ratios and concentrations (given in fractions of volume units, VU) of the four different materials (plus residual material) at different noise-levels for 30 keV, 40 keV, 50 keV, 60 keV, 80 keV and 100 keV (monochromatic). Our simulation results were then confirmed by using region of interest-based measurements in a phantom-study at corresponding acceleration voltages. Therefore, different mixtures of contrast materials were scanned using a micro-CT. Voxel wise evaluation of the phantom imaging data was conducted to confirm its usability for future imaging applications and to estimate the influence of varying noise-levels, scattering, artifacts and concentrations.The analysis of our simulations showed the smallest deviation of 0.01 (0.003-0.15) VU between given and calculated concentrations of the different contrast materials when using an energy-combination of 30 keV, 40 keV, 50 keV and 100 keV for MECT. Subsequent MECT phantom measurements, however, revealed a combination of acceleration voltages of 30 kV, 40 kV, 60 kV and 100 kV as most effective for performing material decomposition with a deviation of 0.28 (0-1.07) mg/ml. The feasibility of our voxelwise analyses using the proposed algorithm was then confirmed by the generation of phantom parameter-maps that matched the known contrast material concentrations. The results were mostly influenced by the noise-level and the concentrations used in the phantoms.