CT scanner x-ray spectrum estimation from transmission measurements.

PURPOSE 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. METHODS 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. RESULTS 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. CONCLUSIONS 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.

[1]  M Yaffe,et al.  Spectroscopy of diagnostic x rays by a Compton-scatter method. , 1976, Medical physics.

[2]  R Birch,et al.  Computation of bremsstrahlung X-ray spectra and comparison with spectra measured with a Ge(Li) detector. , 1979, Physics in medicine and biology.

[3]  G. Matscheko,et al.  A Compton scattering spectrometer for determining X-ray photon energy spectra. , 1987, Physics in medicine and biology.

[4]  Compton spectroscopy in the diagnostic x-ray energy range. I. Spectrometer design. , 1989, Physics in medicine and biology.

[5]  G A Carlsson,et al.  Compton spectroscopy in the diagnostic x-ray energy range. II. Effects of scattering material and energy resolution. , 1989, Physics in medicine and biology.

[6]  G. Barnes,et al.  Semiempirical model for generating tungsten target x-ray spectra. , 1991, Medical physics.

[7]  L T Hudson,et al.  A curved crystal spectrometer for energy calibration and spectral characterization of mammographic x-ray sources. , 1996, Medical physics.

[8]  L T Hudson,et al.  Flat and curved crystal spectrography for mammographic X-ray sources. , 1996, The British journal of radiology.

[9]  J. Boone,et al.  An accurate method for computer-generating tungsten anode x-ray spectra from 30 to 140 kV. , 1997, Medical physics.

[10]  C. Ruth,et al.  Estimation of a photon energy spectrum for a computed tomography scanner. , 1997, Medical physics.

[11]  Diagnostic x-ray spectra: a comparison of spectra generated by different computational methods with a measured spectrum. , 1998, Medical physics.

[12]  R. Waggener,et al.  X-ray spectra estimation using attenuation measurements from 25 kVp to 18 MV. , 1999, Medical physics.

[13]  W D McDavid,et al.  Modification and benchmarking of MCNP for low-energy tungsten spectra. , 2000, Medical physics.

[14]  M. Stampanoni,et al.  Computer algebra for x-ray spectral reconstruction between 6 and 25 MV. , 2001, Medical physics.

[15]  A Constantinesco,et al.  Evaluation of the use of six diagnostic X-ray spectra computer codes. , 2004, The British journal of radiology.

[16]  Xiaochuan Pan,et al.  A robust method of x-ray source spectrum estimation from transmission measurements: Demonstrated on computer simulated, scatter-free transmission data , 2005 .

[17]  Characteristic features of a high-energy x-ray spectra estimation method based on the Waggener iterative perturbation principle. , 2006, Medical physics.

[18]  S. Suzuki,et al.  Calculation of 10 MV x-ray spectra emitted by a medical linear accelerator using the BFGS quasi-Newton method , 2007, Physics in medicine and biology.