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.

💡 Research Summary

The paper addresses a critical source of image degradation in photon‑counting detector computed tomography (PCCT): spectral inconsistency, which manifests as pixel‑wise variations in counts across energy bins and leads to ring or band artifacts and quantitative errors. Existing correction strategies—X‑ray fluorescence (XRF) calibration, mono‑energy phantom mapping, or polynomial‑based model fitting—either require cumbersome hardware, are limited to specific phantoms, or lack explicit physical modeling of the underlying mechanisms, especially the energy‑threshold bias (ETB).

The authors propose a physics‑based two‑term factorized model that separates spectral inconsistency into (1) a fixed spectral‑skew term (g_i(E)) that captures pixel‑dependent filtration effects (e.g., variations in filter thickness or detector material), and (2) a variable ETB term (\Delta E_{k i}) that accounts for pixel‑wise deviations of the actual energy threshold from its nominal value. The photon count measured at pixel (i) above threshold (k) is expressed as

\