Dataset Information

A cascaded model of spectral distortions due to spectral response effects and pulse pileup effects in a photon-counting x-ray detector for CT.

ABSTRACT: Energy discriminating, photon-counting detectors (PCDs) are an emerging technology for computed tomography (CT) with various potential benefits for clinical CT. The photon energies measured by PCDs can be distorted due to the interactions of a photon with the detector and the interaction of multiple coincident photons. These effects result in distorted recorded x-ray spectra which may lead to artifacts in reconstructed CT images and inaccuracies in tissue identification. Model-based compensation techniques have the potential to account for the distortion effects. This approach requires only a small number of parameters and is applicable to a wide range of spectra and count rates, but it needs an accurate model of the spectral distortions occurring in PCDs. The purpose of this study was to develop a model of those spectral distortions and to evaluate the model using a PCD (model DXMCT-1; DxRay, Inc., Northridge, CA) and various x-ray spectra in a wide range of count rates.The authors hypothesize that the complex phenomena of spectral distortions can be modeled by: (1) separating them into count-rate independent factors that we call the spectral response effects (SRE), and count-rate dependent factors that we call the pulse pileup effects (PPE), (2) developing separate models for SRE and PPE, and (3) cascading the SRE and PPE models into a combined SRE+PPE model that describes PCD distortions at both low and high count rates. The SRE model describes the probability distribution of the recorded spectrum, with a photo peak and a continuum tail, given the incident photon energy. Model parameters were obtained from calibration measurements with three radioisotopes and then interpolated linearly for other energies. The PPE model used was developed in the authors' previous work [K. Taguchi et al., "Modeling the performance of a photon counting x-ray detector for CT: Energy response and pulse pileup effects," Med. Phys. 38(2), 1089-1102 (2011)]. The agreement between the x-ray spectra calculated by the cascaded SRE+PPE model and the measured spectra was evaluated for various levels of deadtime loss ratios (DLR) and incident spectral shapes, realized using different attenuators, in terms of the weighted coefficient of variation (COVW), i.e., the root mean square difference weighted by the statistical errors of the data and divided by the mean.At low count rates, when DLR < 10%, the distorted spectra measured by the DXMCT-1 were in agreement with those calculated by SRE only, with COVW's less than 4%. At higher count rates, the measured spectra were also in agreement with the ones calculated by the cascaded SRE+PPE model; with PMMA as attenuator, COVW was 5.6% at a DLR of 22% and as small as 6.7% for a DLR as high as 55%.The x-ray spectra calculated by the proposed model agreed with the measured spectra over a wide range of count rates and spectral shapes. The SRE model predicted the distorted, recorded spectra with low count rates over various types and thicknesses of attenuators. The study also validated the hypothesis that the complex spectral distortions in a PCD can be adequately modeled by cascading the count-rate independent SRE and the count-rate dependent PPE.

SUBMITTER: Cammin J

PROVIDER: S-EPMC3979165 | biostudies-other | 2014 Apr

REPOSITORIES: biostudies-other

ACCESS DATA

Similar Datasets

Project description:Recently, photon counting x-ray detectors (PCXDs) with energy discrimination capabilities have been developed for potential use in clinical computed tomography (CT) scanners. These PCXDs have great potential to improve the quality of CT images due to the absence of electronic noise and weights applied to the counts and the additional spectral information. With high count rates encountered in clinical CT, however, coincident photons are recorded as one event with a higher or lower energy due to the finite speed of the PCXD. This phenomenon is called a "pulse pileup event" and results in both a loss of counts (called "deadtime losses") and distortion of the recorded energy spectrum. Even though the performance of PCXDs is being improved, it is essential to develop algorithmic methods based on accurate models of the properties of detectors to compensate for these effects. To date, only one PCXD (model DXMCT-1, DxRay, Inc., Northridge, CA) has been used for clinical CT studies. The aim of that study was to evaluate the agreement between data measured by DXMCT-1 and those predicted by analytical models for the energy response, the deadtime losses, and the distorted recorded spectrum caused by pulse pileup effects.An energy calibration was performed using 99mTc (140 keV), 57Co (122 keV), and an x-ray beam obtained with four x-ray tube voltages (35, 50, 65, and 80 kVp). The DXMCT-1 was placed 150 mm from the x-ray focal spot; the count rates and the spectra were recorded at various tube current values from 10 to 500 microA for a tube voltage of 80 kVp. Using these measurements, for each pulse height comparator we estimated three parameters describing the photon energy-pulse height curve, the detector deadtime tau, a coefficient k that relates the x-ray tube current I to an incident count rate a by a = k x I, and the incident spectrum. The mean pulse shape of all comparators was acquired in a separate study and was used in the model to estimate the distorted recorded spectrum. The agreement between data measured by the DXMCT-1 and those predicted by the models was quantified by the coefficient of variation (COV), i.e., the root mean square difference divided by the mean of the measurement.Photon energy versus pulse height curves calculated with an analytical model and those measured using the DXMCT-1 were in agreement within 0.2% in terms of the COV. The COV between the output count rates measured and those predicted by analytical models was 2.5% for deadtime losses of up to 60%. The COVs between spectra measured and those predicted by the detector model were within 3.7%-7.2% with deadtime losses of 19%-46%.It has been demonstrated that the performance of the DXMCT-1 agreed exceptionally well with the analytical models regarding the energy response, the count rate, and the recorded spectrum with pulse pileup effects. These models will be useful in developing methods to compensate for these effects in PCXD-based clinical CT systems.

Project description:To investigate the feasibility of using an image-based method to correct for distortions induced by various artifacts in the x-ray spectrum recorded with photon-counting detectors for their application in breast computed tomography (CT).The polyenergetic incident spectrum was simulated with the tungsten anode spectral model using the interpolating polynomials (TASMIP) code and carefully calibrated to match the x-ray tube in this study. Experiments were performed on a Cadmium-Zinc-Telluride (CZT) photon-counting detector with five energy thresholds. Energy bins were adjusted to evenly distribute the recorded counts above the noise floor. BR12 phantoms of various thicknesses were used for calibration. A nonlinear function was selected to fit the count correlation between the simulated and the measured spectra in the calibration process. To evaluate the proposed spectral distortion correction method, an empirical fitting derived from the calibration process was applied on the raw images recorded for polymethyl methacrylate (PMMA) phantoms of 8.7, 48.8, and 100.0 mm. Both the corrected counts and the effective attenuation coefficient were compared to the simulated values for each of the five energy bins. The feasibility of applying the proposed method to quantitative material decomposition was tested using a dual-energy imaging technique with a three-material phantom that consisted of water, lipid, and protein. The performance of the spectral distortion correction method was quantified using the relative root-mean-square (RMS) error with respect to the expected values from simulations or areal analysis of the decomposition phantom.The implementation of the proposed method reduced the relative RMS error of the output counts in the five energy bins with respect to the simulated incident counts from 23.0%, 33.0%, and 54.0% to 1.2%, 1.8%, and 7.7% for 8.7, 48.8, and 100.0 mm PMMA phantoms, respectively. The accuracy of the effective attenuation coefficient of PMMA estimate was also improved with the proposed spectral distortion correction. Finally, the relative RMS error of water, lipid, and protein decompositions in dual-energy imaging was significantly reduced from 53.4% to 6.8% after correction was applied.The study demonstrated that dramatic distortions in the recorded raw image yielded from a photon-counting detector could be expected, which presents great challenges for applying the quantitative material decomposition method in spectral CT. The proposed semi-empirical correction method can effectively reduce these errors caused by various artifacts, including pulse pileup and charge sharing effects. Furthermore, rather than detector-specific simulation packages, the method requires a relatively simple calibration process and knowledge about the incident spectrum. Therefore, it may be used as a generalized procedure for the spectral distortion correction of different photon-counting detectors in clinical breast CT systems.

Project description:This paper describes a noniterative estimator for the energy dependent information from photon counting detectors with multibin pulse height analysis (PHA). The estimator uses the two function decomposition of the attenuation coefficient [R. E. Alvarez and A. Macovski, Phys. Med. Biol. 21, 733-744 (1976)] and its output is the line integrals of the basis set coefficients. The output noise variance and bias is compared to other noniterative estimators and to the Cramèr-Rao lower bound (CRLB).The estimator first computes an initial estimate from a linearized maximum likelihood estimator. The errors in the initial estimates are determined at a set of points from measurements on a calibration phantom. The errors at these known points are interpolated to create two-dimensional look up tables of corrections to the initial estimates. During image acquisition, the linearized maximum likelihood estimate for each data point is used as an input to the correction look up tables, and the final output is the sum of the estimate and the correction. The performance of the estimator is compared to generalizations of the polynomial and rational polynomial estimators for multibin data. The estimators are compared by the mean square error (MSE) and its components, the bias, and the variance of the output. The variance is also compared to the CRLB. The performance is simulated with two to five bins PHA data. The CRLB at a fixed object thickness is also computed as a function of the number of bins.For two bin data, all the estimators' variances are equal to the CRLB. With three or more bins, only the proposed estimator achieves the CRLB while the others, which were not optimized for noise performance, have much larger output variance. The bias of the proposed estimator is equal to the polynomial estimator for calibration phantoms with 40 or more steps, that is, 1600 combinations of basis materials, but is larger than the rational polynomial bias. In all cases at the photon counts tested, the MSE is essentially equal to the variance, indicating that the bias errors are negligible compared to the variance.The estimator provides a noniterative method to compute the energy dependent information from multibin PHA data that achieves the CRLB over a wide range of operating conditions and has low output bias. The estimator can be calibrated based on the measurements of a calibration phantom; so, it does not require measurements of the x-ray energy spectrum or the detector response functions.

Project description:PurposePhoton counting detectors (PCDs) are being introduced in advanced x-ray computed tomography (CT) scanners. From a single PCD-CT acquisition, multiple images can be reconstructed, each based on only a part of the original x-ray spectrum. In this study, we investigated whether PCD-CT can be used to estimate stopping power ratios (SPRs) for proton therapy treatment planning, both by comparing to other SPR methods proposed for single energy CT (SECT) and dual energy CT (DECT) as well as to experimental measurements.MethodsA previously developed DECT-based SPR estimation method was adapted to PCD-CT data, by adjusting the estimation equations to allow for more energy spectra. The method was calibrated directly on noisy data to increase the robustness toward image noise. The new PCD SPR estimation method was tested in theoretical calculations as well as in an experimental setup, using both four and two energy bin PCD-CT images, and through comparison to two other SPR methods proposed for SECT and DECT. These two methods were also evaluated on PCD-CT images, full spectrum (one-bin) or two-bin images, respectively. In a theoretical framework, we evaluated the effect of patient-specific tissue variations (density and elemental composition) and image noise on the SPR accuracy; the latter effect was assessed by applying three different noise levels (low, medium, and high noise). SPR estimates derived using real PCD-CT images were compared to experimentally measured SPRs in nine organic tissue samples, including fat, muscle, and bone tissues.ResultsFor the theoretical calculations, the root-mean-square error (RMSE) of the SPR estimation was 0.1% for the new PCD method using both two and four energy bins, compared to 0.2% and 0.7% for the DECT- and SECT-based method, respectively. The PCD method was found to be very robust toward CT image noise, with a RMSE of 2.7% when high noise was added to the CT numbers. Introducing tissue variations, the RMSE only increased to 0.5%; even when adding high image noise to the changed tissues, the RMSE stayed within 3.1%. In the experimental measurements, the RMSE over the nine tissue samples was 0.8% when using two energy bins, and 1.0% for the four-bin images.ConclusionsIn all tested cases, the new PCD method produced similar or better results than the SECT- and DECT-based methods, showing an overall improvement of the SPR accuracy. This study thus demonstrated that PCD-CT scans will be a qualified candidate for SPR estimations.

Dataset Information

A cascaded model of spectral distortions due to spectral response effects and pulse pileup effects in a photon-counting x-ray detector for CT.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets