Dataset Information

Estimator for photon counting energy selective x-ray imaging with multibin pulse height analysis.

ABSTRACT: 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.

SUBMITTER: Alvarez RE

PROVIDER: S-EPMC3098892 | biostudies-other | 2011 May

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:The objective of the study was to demonstrate that, in x-ray computed tomography (CT), more than two types of materials can be effectively separated with the use of an energy resolved photon-counting detector and classification methodology. Specifically, this applies to the case when contrast agents that contain K-absorption edges in the energy range of interest are present in the object. This separation is enabled via the use of recently developed energy resolved photon-counting detectors with multiple thresholds, which allow simultaneous measurements of the x-ray attenuation at multiple energies.To demonstrate this capability, we performed simulations and physical experiments using a six-threshold energy resolved photon-counting detector. We imaged mouse-sized cylindrical phantoms filled with several soft-tissue-like and bone-like materials and with iodine-based and gadolinium-based contrast agents. The linear attenuation coefficients were reconstructed for each material in each energy window and were visualized as scatter plots between pairs of energy windows. For comparison, a dual-kVp CT was also simulated using the same phantom materials. In this case, the linear attenuation coefficients at the lower kVp were plotted against those at the higher kVp.In both the simulations and the physical experiments, the contrast agents were easily separable from other soft-tissue-like and bone-like materials, thanks to the availability of the attenuation coefficient measurements at more than two energies provided by the energy resolved photon-counting detector. In the simulations, the amount of separation was observed to be proportional to the concentration of the contrast agents; however, this was not observed in the physical experiments due to limitations of the real detector system. We used the angle between pairs of attenuation coefficient vectors in either the 5-D space (for non-contrast-agent materials using energy resolved photon-counting acquisition) or a 2-D space (for contrast agents using energy resolved photon-counting acquisition and all materials using dual-kVp acquisition) as a measure of the degree of separation. Compared to dual-kVp techniques, an energy resolved detector provided a larger separation and the ability to separate different target materials using measurements acquired in different energy window pairs with a single x-ray exposure.We concluded that x-ray CT with an energy resolved photon-counting detector with more than two energy windows allows the separation of more than two types of materials, e.g., soft-tissue-like, bone-like, and one or more materials with K-edges in the energy range of interest. Separating material types using energy resolved photon-counting detectors has a number of advantages over dual-kVp CT in terms of the degree of separation and the number of materials that can be separated simultaneously.

Project description:To develop and test a method to quantify the effect of dimensionality on the noise in energy selective x-ray imaging.The Cramèr-Rao lower bound (CRLB), a universal lower limit of the covariance of any unbiased estimator, is used to quantify the noise. It is shown that increasing dimensionality always increases, or at best leaves the same, the variance. An analytic formula for the increase in variance in an energy selective x-ray system is derived. The formula is used to gain insight into the dependence of the increase in variance on the properties of the additional basis functions, the measurement noise covariance, and the source spectrum. The formula is also used with computer simulations to quantify the dependence of the additional variance on these factors. Simulated images of an object with three materials are used to demonstrate the trade-off of increased information with dimensionality and noise. The images are computed from energy selective data with a maximum likelihood estimator.The increase in variance depends most importantly on the dimension and on the properties of the additional basis functions. With the attenuation coefficients of cortical bone, soft tissue, and adipose tissue as the basis functions, the increase in variance of the bone component from two to three dimensions is 1.4 × 10(3). With the soft tissue component, it is 2.7 × 10(4). If the attenuation coefficient of a high atomic number contrast agent is used as the third basis function, there is only a slight increase in the variance from two to three basis functions, 1.03 and 7.4 for the bone and soft tissue components, respectively. The changes in spectrum shape with beam hardening also have a substantial effect. They increase the variance by a factor of approximately 200 for the bone component and 220 for the soft tissue component as the soft tissue object thickness increases from 1 to 30 cm. Decreasing the energy resolution of the detectors increases the variance of the bone component markedly with three dimension processing, approximately a factor of 25 as the resolution decreases from 100 to 3 bins. The increase with two dimension processing for adipose tissue is a factor of two and with the contrast agent as the third material for two or three dimensions is also a factor of two for both components. The simulated images show that a maximum likelihood estimator can be used to process energy selective x-ray data to produce images with noise close to the CRLB.The method presented can be used to compute the effects of the object attenuation coefficients and the x-ray system properties on the relationship of dimensionality and noise in energy selective x-ray imaging 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.

Dataset Information

Estimator for photon counting energy selective x-ray imaging with multibin pulse height analysis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets