Dataset Information

CT scanner x-ray spectrum estimation from transmission measurements.

ABSTRACT: In diagnostic CT imaging, multiple important applications depend on the knowledge of the x-ray spectrum, including Monte Carlo dose calculations and dual-energy material decomposition analysis. Due to the high photon flux involved, it is difficult to directly measure spectra from the x-ray tube of a CT scanner. One potential method for indirect measurement involves estimating the spectrum from transmission measurements. The expectation maximization (EM) method is an accurate and robust method to solve this problem. In this article, this method was evaluated in a commercial CT scanner.Two step-wedges (polycarbonate and aluminum) were used to produce different attenuation levels. Transmission measurements were performed on the scanner and the measured data from the scanner were exported to an external computer to calculate the spectra. The EM method was applied to solve the equations that represent the attenuation processes of polychromatic x-ray photons. Estimated spectra were compared to the spectra simulated using a software provided by the manufacturer of the scanner. To test the accuracy of the spectra, a verification experiment was performed using a phantom containing different depths of water. The measured transmission data were compared to the transmission values calculated using the estimated spectra.Spectra of 80, 100, 120, and 140 kVp from a dual-source CT scanner were estimated. The estimated and simulated spectra were well matched. The differences of mean energies were less than 1 keV. In the verification experiment, the measured and calculated transmission values were in excellent agreement.Spectrum estimation using transmission data and the EM method is a quantitatively accurate and robust technique to estimate the spectrum of a CT system. This method could benefit studies relying on accurate knowledge of the x-ray spectra from CT scanner.

SUBMITTER: Duan X

PROVIDER: S-EPMC3041810 | biostudies-other | 2011 Feb

REPOSITORIES: biostudies-other

ACCESS DATA

Similar Datasets

Project description:Virtual imaging trials (VITs), defined as the process of conducting clinical imaging trials using computer simulations, offer a time- and cost-effective alternative to traditional imaging trials for CT. The clinical potential of VITs hinges on the realism of simulations modeling the image acquisition process, where the accurate scanner-specific simulation of scatter in a time-feasible manner poses a particular challenge. To meet this need, this study proposes, develops, and validates a rapid scatter estimation framework, based on GPU-accelerated Monte Carlo (MC) simulations and denoising methods, for estimating scatter in single source, dual-source, and photon-counting CT. A CT simulator was developed to incorporate parametric models for an anti-scatter grid and a curved energy integrating detector with an energy-dependent response. The scatter estimates from the simulator were validated using physical measurements acquired on a clinical CT system using the standard single-blocker method. The MC simulator was further extended to incorporate a pre-validated model for a PCD and an additional source-detector pair to model cross scatter in dual-source configurations. To estimate scatter with desirable levels of statistical noise using a manageable computational load, two denoising methods using a (1) convolutional neural network and an (2) optimized Gaussian filter were further deployed. The viability of this framework for clinical VITs was assessed by integrating it with a scanner-specific ray-tracer program to simulate images for an image quality (Mercury) and an anthropomorphic phantom (XCAT). The simulated scatter-to-primary ratios agreed with physical measurements within 4.4% ± 10.8% across all projection angles and kVs. The differences of ∼121 HU between images with and without scatter, signifying the importance of scatter for simulating clinical images. The denoising methods preserved the magnitudes and trends observed in the reference scatter distributions, with an averaged rRMSE value of 0.91 and 0.97 for the two methods, respectively. The execution time of ∼30 s for simulating scatter in a single projection with a desirable level of statistical noise indicates a major improvement in performance, making our tool an eligible candidate for conducting extensive VITs spanning multiple patients and scan protocols.

Project description:Computed Tomography (CT) imaging of the lung, reported in Hounsfield Units (HU), can be parameterized as a quantitative image biomarker for the diagnosis and monitoring of lung density changes due to emphysema, a type of chronic obstructive pulmonary disease (COPD). CT lung density metrics are global measurements based on lung CT number histograms, and are typically a quantity specifying either the percentage of voxels with CT numbers below a threshold, or a single CT number below which a fixed relative lung volume, nth percentile, falls. To reduce variability in the density metrics specified by CT attenuation, the Quantitative Imaging Biomarkers Alliance (QIBA) Lung Density Committee has organized efforts to conduct phantom studies in a variety of scanner models to establish a baseline for assessing the variations in patient studies that can be attributed to scanner calibration and measurement uncertainty.Data were obtained from a phantom study on CT scanners from four manufacturers with several protocols at various tube potential voltage (kVp) and exposure settings. Free from biological variation, these phantom studies provide an assessment of the accuracy and precision of the density metrics across platforms solely due to machine calibration and uncertainty of the reference materials. The phantom used in this study has three foam density references in the lung density region, which, after calibration against a suite of Standard Reference Materials (SRM) foams with certified physical density, establishes a HU-electron density relationship for each machine-protocol. We devised a 5-step calibration procedure combined with a simplified physical model that enabled the standardization of the CT numbers reported across a total of 22 scanner-protocol settings to a single energy (chosen at 80 keV). A standard deviation was calculated for overall CT numbers for each density, as well as by scanner and other variables, as a measure of the variability, before and after the standardization. In addition, a linear mixed-effects model was used to assess the heterogeneity across scanners, and the 95% confidence interval of the mean CT number was evaluated before and after the standardization.We show that after applying the standardization procedures to the phantom data, the instrumental reproducibility of the CT density measurement of the reference foams improved by more than 65%, as measured by the standard deviation of the overall mean CT number. Using the lung foam that did not participate in the calibration as a test case, a mixed effects model analysis shows that the 95% confidence intervals are [-862.0 HU, -851.3 HU] before standardization, and [-859.0 HU, -853.7 HU] after standardization to 80 keV. This is in general agreement with the expected CT number value at 80 keV of -855.9 HU with 95% CI of [-857.4 HU, -854.5 HU] based on the calibration and the uncertainty in the SRM certified density.This study provides a quantitative assessment of the variations expected in CT lung density measures attributed to non-biological sources such as scanner calibration and scanner x-ray spectrum and filtration. By removing scanner-protocol dependence from the measured CT numbers, higher accuracy and reproducibility of quantitative CT measures were attainable. The standardization procedures developed in study may be explored for possible application in CT lung density clinical data.

Project description:PURPOSE:This paper lays the groundwork for linking Hounsfield unit measurements to the International System of Units (SI), ultimately enabling traceable measurements across X-ray CT (XCT) machines. We do this by characterizing a material basis that may be used in XCT reconstruction giving linear combinations of concentrations of chemical elements (in the SI units of mol/m3) which may be observed at each voxel. By implication, linear combinations not in the set are not observable. METHODS AND MATERIALS:We formulated a model for our material basis with a set of measurements of elemental powders at four tube voltages, 80 kV, 100 kV, 120 kV, and 140 kV, on a medical XCT. The samples included 30 small plastic bottles of powders containing various compounds spanning the atomic numbers up to 20, and a bottle of water and one of air. Using the chemical formulas and measured masses, we formed a matrix giving the number of Hounsfield units per (mole per cubic meter) at each tube voltage for each of 13 chemical elements. We defined a corresponding matrix in units we call molar Hounsfield unit (HU) potency, the difference in HU values that an added mole per cubic meter in a given voxel would add to the measured HU value. We built a matrix of molar potencies for each chemical element and tube voltage and performed a singular value decomposition (SVD) on these to formulate our material basis. We determined that the dimension of this basis is two. We then compared measurements in this material space with theoretical measurements, combining XCOM cross section data with the tungsten anode spectral model using interpolating cubic splines (TASMICS), a one-parameter filter, and a simple detector model, creating a matrix similar to our experimental matrix for the first 20 chemical elements. Finally, we compared the model predictions to Hounsfield unit measurements on three XCT calibration phantoms taken from the literature. RESULTS:We predict the experimental HU potency values derived from our scans of chemical elements with our theoretical model built from XCOM data. The singular values and singular vectors of the model and powder measurements are in substantial agreement. Application of the Bayesian Information Criterion (BIC) shows that exactly two singular values and singular vectors describe the results over four tube voltages. We give a good account of the HU values from the literature, measured for the calibration phantoms at several tube voltages for several commercial instruments, compared with our theoretical model without introducing additional parameters. CONCLUSIONS:We have developed a two-dimensional material basis that specifies the degree to which individual elements in compounds effect the HU values in XCT images of samples with elements up to atomic number Z = 20. We show that two dimensions is sufficient given the contrast and noise in our experiment. The linear combination of concentrations of elements that can be observed using a medical XCT have been characterized, providing a material basis for use in dual-energy reconstruction. This approach provides groundwork for improved reconstruction and for the link of Hounsfield units to the SI.

Dataset Information

CT scanner x-ray spectrum estimation from transmission measurements.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets