Improving model-data mismatch for photon-counting detector model using global and local model parameters

BACKGROUND

An energy-discriminating capability of a photon counting detector (PCD) can provide many clinical advantages, but several factors, such as charge sharing (CS) and pulse pileup (PP), degrade the capability by distorting the measured x-ray spectrum. To fully exploit the merits of PCDs, it is important to characterize the output of PCDs. Previously proposed PCD output models showed decent agreement with physical PCDs; however, there were still scopes to be improved: a global model-data mismatch and pixel-to-pixel variations.

PURPOSES

In this study, we improve a PCD model by using count-rate-dependent model parameters to address the issues and evaluate agreement against physical PCDs.

METHODS

The proposed model is based on the cascaded model, and we made model parameters condition-dependent and pixel-specific to deal with the global model-data mismatch and the pixel-to-pixel variation. The parameters are determined by a procedure for model parameter estimation with data acquired from different thicknesses of water or aluminum at different x-ray tube currents. To analyze the effects of having proposed model parameters, we compared three setups of our model: a model with default parameters, a model with global parameters, and a model with global-and-local parameters. For experimental validation, we used CdZnTe-based PCDs, and assessed the performance of the models by calculating the mean absolute percentage errors (MAPEs) between the model outputs and the actual measurements from low count-rates to high count-rates, which have deadtime losses of up to 24%.

RESULTS

The outputs of the proposed model visually matched well with the PCD measurements for all test data. For the test data, the MAPEs averaged over all the bins were 49.2-51.1% for a model with default parameters, 8.0-9.8% for a model with the global parameters, and 1.2-2.7% for a model with the global-and-local parameters.

CONCLUSION

The proposed model can estimate the outputs of physical PCDs with high accuracy from low to high count-rates. We expect that our model will be actively utilized in applications where the pixel-by-pixel accuracy of a PCD model is important.

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)² × 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.

A data-driven maximum likelihood classification for nanoparticle agent identification in photon-counting CT

The nanoparticle agent, combined with a targeting factor reacting with lesions, enables specific CT imaging. Thus, the identification of the nanoparticle agents has the potential to improve clinical diagnosis. Thanks to the energy sensitivity of the photon-counting detector (PCD), it can exploit the K-edge of the nanoparticle agents in the clinical x-ray energy range to identify the agents. In this paper, we propose a novel data-driven approach for nanoparticle agent identification using the PCD. We generate two sets of training data consisting of PCD measurements from calibration phantoms, one in the presence of nanoparticle agent and the other in the absence of the agent. For a given sinogram of PCD counts, the proposed method calculates the normalized log-likelihood sinogram for each class (class 1: with the agent, class 2: without the agent) using theKnearest neighbors (KNN) estimator, backproject the sinograms, and compare the backprojection images to identify the agent. We also proved that the proposed algorithm is equivalent to the maximum likelihood-based classification. We studied the robustness of dose reduction with gold nanoparticles as the K-edge contrast media and demonstrated that the proposed method identifies targets with different concentrations of the agents without background noise.

Advanced Monte Carlo simulations of emission tomography imaging systems with GATE

Built on top of the Geant4 toolkit, GATE is collaboratively developed for more than 15 years to design Monte Carlo simulations of nuclear-based imaging systems. It is, in particular, used by researchers and industrials to design, optimize, understand and create innovative emission tomography systems. In this paper, we reviewed the recent developments that have been proposed to simulate modern detectors and provide a comprehensive report on imaging systems that have been simulated and evaluated in GATE. Additionally, some methodological developments that are not specific for imaging but that can improve detector modeling and provide computation time gains, such as Variance Reduction Techniques and Artificial Intelligence integration, are described and discussed.