Energy-Threshold Bias Calculator: A Physics-Model-Based Adaptive Correction Scheme for Photon-Counting CT

Photon-counting detector based computed tomography (PCCT) has greatly advanced in recent years. However, spectral inconsistency, referring to inter-pixel variations in detected counts per energy bin, can easily leads to ring or band artifacts and inaccuracies in CT reconstructed images. This work proposes a novel physics-model based method to correct for spectral inconsistency by modeling it through two terms: 1) a fixed spectral skew term (energy threshold-independent filtration function) determined at a given energy threshold, and 2) a variable energy-threshold bias term that can be directly calculated by using our spectral model as the threshold changes. After the two terms being computed out in the calibration stage, they will be incorporated into our spectral model to adaptively generate the spectral correction vectors as well as the material decomposition vectors if needed, pixel-by-pixel for PCCT projection data. Using a minimum set of parameters with explicit physics meaning, such an energy-threshold bias calculator (ETB-Cal) has advantages of computational efficiency, robustness in implementation, and convenience with no need of X-ray fluorescence materials in calibration. To validate our method, both numerical simulations and physical experiments using multiple phantoms were carried out on a tabletop PCCT system, with preliminary results showing a significant reduction in non-uniformity, from 29.3 to 5.8 HU for Gammex multi-energy phantom versus no correction (comparatively, 8.3 HU was achieved by a polynomial-involving model-based approach with no explicit modeling and calculating of energy threshold bias but more calibration data required), and from 27.9 to 3.2 HU for the Kyoto head phantom.