Model-based pulse pileup and charge sharing compensation for photon counting detectors: A simulation study.

PURPOSE: We aim at developing a model-based algorithm that compensates for the effect of both pulse pileup (PP) and charge sharing (CS) and evaluates the performance using computer simulations. METHODS: The proposed PCP algorithm for PP and CS compensation uses cascaded models for CS and PP we previously developed, maximizes Poisson log-likelihood, and uses an efficient three-step exhaustive search. For comparison, we also developed an LCP algorithm that combines models for a loss of counts (LCs) and CS. Two types of computer simulations, slab- and computed tomography (CT)-based, were performed to assess the performance of both PCP and LCP with 200 and 800 mA, (300 µm)2  × 1.6-mm cadmium telluride detector, and a dead-time of 23 ns. A slab-based assessment used a pair of adipose and iodine with different thicknesses, attenuated X-rays, and assessed the bias and noise of the outputs from one detector pixel; a CT-based assessment simulated a chest/cardiac scan and a head-and-neck scan using 3D phantom and noisy cone-beam projections. RESULTS: With the slab simulation, the PCP had little or no biases when the expected counts were sufficiently large, even though a probability of count loss (PCL) due to dead-time loss or PP was as high as 0.8. In contrast, the LCP had significant biases (>±2 cm of adipose) when the PCL was higher than 0.15. Biases were present with both PCP and LCP when the expected counts were less than 10-120 per datum, which was attributed to the fact that the maximum likelihood did not approach the asymptote. The noise of PCP was within 8% from the Cramér-Rao lower bounds for most cases when no significant bias was present. The two CT studies essentially agreed with the slab simulation study. PCP had little or no biases in the estimated basis line integrals, reconstructed basis density maps, and synthesized monoenergetic CT images. But the LCP had significant biases in basis line integrals when X-ray beams passed through lungs and near the body and neck contours, where the PCLs were above 0.15. As a consequence, basis density maps and monoenergetic CT images obtained by LCP had biases throughout the imaged space. CONCLUSION: We have developed the PCP algorithm that uses the PP-CS model. When the expected counts are more than 10-120 per datum, the PCP algorithm is statistically efficient and successfully compensates for the effect of the spectral distortion due to both PP and CS providing little or no biases in basis line integrals, basis density maps, and monoenergetic CT images regardless of count-rates. In contrast, the LCP algorithm, which models an LC due to pileup, produces severe biases when incident count-rates are high and the PCL is 0.15 or higher.

Spatio-energetic cross-talk in photon counting detectors: N × N binning and sub-pixel masking.

PURPOSE: Smaller pixel sizes of x-ray photon counting detectors (PCDs) benefit count rate capabilities but increase cross-talk and "double-counting" between neighboring PCD pixels. When an x-ray photon produces multiple (n) counts at neighboring (sub-)pixels and they are added during post-acquisition N × N binning process, the variance of the final PCD output-pixel will be larger than its mean. In the meantime, anti-scatter grids are placed at the pixel boundaries in most of x-ray CT systems and will decrease cross-talk between sub-pixels because the grids mask sub-pixels underneath them, block the primary x-rays, and increase the separation distance between active sub-pixels. The aim of this paper was, first, to study the PCD statistics with various N × N binning schemes and three different masking methods in the presence of cross-talks, and second, to assess one of the most fundamental performances of x-ray CT: soft tissue contrast visibility. METHODS: We used a PCD cross-talk model (Photon counting toolkit, PcTK) and produced cross-talk data between 3 × 3 neighboring sub-pixels and calculated the mean, variance, and covariance of output-pixels with each of N × N binning scheme [4 × 4 binning, 2 × 2 binning, and 1 × 1 binning (i.e., no binning)] and three different sub-pixel masking methods (no mask, 1-D mask, and 2-D mask). We then set up simulation to evaluate the soft tissue contrast visibility. X-rays of 120 kVp were attenuated by 10-40 cm-thick water, with the right side of PCDs having 0.5 cm thicker water than the left side. A pair of output-pixels across the left-right boundary were used to assess the sensitivity index (SI or d'), which typically ranges 0-1 and is a generalized signal-to-noise ratio and a statistics used in signal detection theory. RESULTS: Binning a larger number of sub-pixels resulted in larger mean counts and larger variance-to-mean ratio when the lower threshold of the energy window was lower than the half of the incident energy. Mean counts are in the order of no mask (the largest), 1-D mask, and 2-D mask but the difference in variance-to-mean ratio was small. For a given sub-pixel size and masking method, binning more sub-pixels degraded the normalized SI values but the difference between 4 × 4 binning and 1 × 1 binning was typically less than 0.06. 1-D mask provided better normalized SI values than no mask and 2-D mask for side-by-side case and the improvements were larger with fewer binnings, although the difference was less than 0.10. 2-D mask was the best for embedded case. The normalized SI values of combined binning, sub-pixel size, and masking were in the order of 1 × 1 (900 µm)2 binning, 2 × 2 (450 µm)2 binning, and 4 × 4 (225 µm)2 binning for a given masking method but the difference between each of them were typically 0.02-0.05. CONCLUSION: We have evaluated the effect of double-counting between PCD sub-pixels with various binning and masking methods. SI values were better with fewer number of binning and larger sub-pixels. The difference among various binning and masking methods, however, was typically less than 0.06, which might result in a dose penalty of 13% if the CT system were linear.

Spatio-energetic cross-talk in photon counting detectors: Numerical detector model (PcTK) and workflow for CT image quality assessment.

PURPOSE: The interpixel cross-talk of energy-sensitive photon counting x-ray detectors (PCDs) has been studied and an analytical model (version 2.1) has been developed for double-counting between neighboring pixels due to charge sharing and K-shell fluorescence x-ray emission followed by its reabsorption (Taguchi K, et al., Medical Physics 2016;43(12):6386-6404). While the model version 2.1 simulated the spectral degradation well, it had the following problems that has been found to be significant recently: (1) The spectrum is inaccurate with smaller pixel sizes; (2) the charge cloud size must be smaller than the pixel size; (3) the model underestimates the spectrum/counts for 10-40 keV; and (4) the model version 2.1 cannot handlen-tuple-counting withn > 2 (i.e., triple-counting or higher). These problems are inherent to the design of the model version 2.1; therefore, we developed a new model and addressed these problems in this study. METHODS: We propose a new PCD cross-talk model (version 3.2; Pc TK for "photon counting toolkit") that is based on a completely different design concept from the previous version. It uses a numerical approach and starts with a 2-D model of charge sharing (as opposed to an analytical approach and a 1-D model with version 2.1) and addresses all of the four problems. The model takes the following factors into account: (1) shift-variant electron density of the charge cloud (Gaussian-distributed), (2) detection efficiency, (3) interactions between photons and PCDs via photoelectric effect, and (4) electronic noise. Correlated noisy PCD data can be generated using either a multivariate normal random number generator or a Poisson random number generator. The effect of the two parameters, the effective charge cloud diameter (d0 ) and pixel size (dpix ), was studied and results were compared with Monte Carlo simulations and the previous model version 2.1. Finally, a script for the workflow for CT image quality assessment has been developed, which started with a few material density images, generated material-specific sinogram (line integrals) data, noisy PCD data with spectral distortion using the model version 3.2, and reconstructed PCD- CT images for four energy windows. RESULTS: The model version 3.2 addressed all of the four problems listed above. The spectra withdpix  = 56-113 µm agreed with that of Medipix3 detector withdpix  = 55-110 µm without charge summing mode qualitatively. The counts for 10-40 keV were larger than the previous model (version 2.1) and agreed with MC simulations very well (root-mean-square difference values with model version 3.2 were decreased to 16%-67% of the values with version 2.1). There were many non-zero off-diagonal elements withn-tuple-counting withn > 2 in the normalized covariance matrix of 3 × 3 neighboring pixels. Reconstructed images showed biases and artifacts attributed to the spectral distortion due to the charge sharing and fluorescence x rays. CONCLUSION: We have developed a new PCD model for spatio-energetic cross-talk and correlation between PCD pixels. The workflow demonstrated the utility of the model for general or task-specific image quality assessments for the PCD- CT.Note: The program (Pc TK) and the workflow scripts have been made available to academic researchers. Interested readers should visit the website (pctk.jhu.edu) or contact the corresponding author.

Estimation of Basis Line-Integrals in a Spectral Distortion-Modeled Photon Counting Detector Using Low-Rank Approximation-Based X-Ray Transmittance Modeling: K-Edge Imaging Application.

Photon counting detectors (PCDs) provide multiple energy-dependent measurements for estimating basis line-integrals. However, the measured spectrum is distorted from the spectral response effect (SRE) via charge sharing, K-fluorescence emission, and so on. Thus, in order to avoid bias and artifacts in images, the SRE needs to be compensated. For this purpose, we recently developed a computationally efficient three-step algorithm for PCD-CT without contrast agents by approximating smooth X-ray transmittance using low-order polynomial bases. It compensated the SRE by incorporating the SRE model in a linearized estimation process and achieved nearly the minimum variance and unbiased (MVU) estimator. In this paper, we extend the three-step algorithm to K-edge imaging applications by designing optimal bases using a low-rank approximation to model X-ray transmittances with arbitrary shapes (i.e., smooth without the K-edge or discontinuous with the K-edge). The bases can be used to approximate the X-ray transmittance and to linearize the PCD measurement modeling and then the three-step estimator can be derived as in the previous approach: estimating the x-ray transmittance in the first step, estimating basis line-integrals including that of the contrast agent in the second step, and correcting for a bias in the third step. We demonstrate that the proposed method is more accurate and stable than the low-order polynomial-based approaches with extensive simulation studies using gadolinium for the K-edge imaging application. We also demonstrate that the proposed method achieves nearly MVU estimator, and is more stable than the conventional maximum likelihood estimator in high attenuation cases with fewer photon counts.

Estimation of Basis Line-Integrals in a Spectral Distortion-Modeled Photon Counting Detector Using Low-Order Polynomial Approximation of X-ray Transmittance.

Photon counting detector (PCD)-based computed tomography exploits spectral information from a transmitted x-ray spectrum to estimate basis line-integrals. The recorded spectrum, however, is distorted and deviates from the transmitted spectrum due to spectral response effect (SRE). Therefore, the SRE needs to be compensated for when estimating basis line-integrals. One approach is to incorporate the SRE model with an incident spectrum into the PCD measurement model and the other approach is to perform a calibration process that inherently includes both the SRE and the incident spectrum. A maximum likelihood estimator can be used to the former approach, which guarantees asymptotic optimality; however, a heavy computational burden is a concern. Calibration-based estimators are a form of the latter approach. They can be very efficient; however, a heuristic calibration process needs to be addressed. In this paper, we propose a computationally efficient three-step estimator for the former approach using a low-order polynomial approximation of x-ray transmittance. The low-order polynomial approximation can change the original non-linear estimation method to a two-step linearized approach followed by an iterative bias correction step. We show that the calibration process is required only for the bias correction step and prove that it converges to the unbiased solution under practical assumptions. Extensive simulation studies validate the proposed method and show that the estimation results are comparable to those of the ML estimator while the computational time is reduced substantially.

Spatio-energetic cross talk in photon counting detectors: Detector model and correlated Poisson data generator.

PURPOSE: An x-ray photon interacts with photon counting detectors (PCDs) and generates an electron charge cloud or multiple clouds. The clouds (thus, the photon energy) may be split between two adjacent PCD pixels when the interaction occurs near pixel boundaries, producing a count at both of the pixels. This is called double-counting with charge sharing. (A photoelectric effect with K-shell fluorescence x-ray emission would result in double-counting as well). As a result, PCD data are spatially and energetically correlated, although the output of individual PCD pixels is Poisson distributed. Major problems include the lack of a detector noise model for the spatio-energetic cross talk and lack of a computationally efficient simulation tool for generating correlated Poisson data. A Monte Carlo (MC) simulation can accurately simulate these phenomena and produce noisy data; however, it is not computationally efficient. METHODS: In this study, the authors developed a new detector model and implemented it in an efficient software simulator that uses a Poisson random number generator to produce correlated noisy integer counts. The detector model takes the following effects into account: (1) detection efficiency; (2) incomplete charge collection and ballistic effect; (3) interaction with PCDs via photoelectric effect (with or without K-shell fluorescence x-ray emission, which may escape from the PCDs or be reabsorbed); and (4) electronic noise. The correlation was modeled by using these two simplifying assumptions: energy conservation and mutual exclusiveness. The mutual exclusiveness is that no more than two pixels measure energy from one photon. The effect of model parameters has been studied and results were compared with MC simulations. The agreement, with respect to the spectrum, was evaluated using the reduced χ2 statistics or a weighted sum of squared errors, χred2(≥1), where χred2=1 indicates a perfect fit. RESULTS: The model produced spectra with flat field irradiation that qualitatively agree with previous studies. The spectra generated with different model and geometry parameters allowed for understanding the effect of the parameters on the spectrum and the correlation of data. The agreement between the model and MC data was very strong. The mean spectra with 90 keV and 140 kVp agreed exceptionally well: χred2 values were 1.049 with 90 keV data and 1.007 with 140 kVp data. The degrees of cross talk (in terms of the relative increase from single pixel irradiation to flat field irradiation) were 22% with 90 keV and 19% with 140 kVp for MC simulations, while they were 21% and 17%, respectively, for the model. The covariance was in strong agreement qualitatively, although it was overestimated. The noisy data generation was very efficient, taking less than a CPU minute as opposed to CPU hours for MC simulators. CONCLUSIONS: The authors have developed a novel, computationally efficient PCD model that takes into account double-counting and resulting spatio-energetic correlation between PCD pixels. The MC simulation validated the accuracy.