Segmentation-free x-ray energy spectrum estimation for computed tomography using dual-energy material decomposition

Abstract. An x-ray energy spectrum plays an essential role in computed tomography (CT) imaging and related tasks. Because of the high photon flux of clinical CT scanners, most of the spectrum estimation methods are indirect and usually suffer from various limitations. In this study, we aim to provide a segmentation-free, indirect transmission measurement–based energy spectrum estimation method using dual-energy material decomposition. The general principle of this method is to minimize the quadratic error between the polychromatic forward projection and the raw projection to calibrate a set of unknown weights, which are used to express the unknown spectrum together with a set of model spectra. The polychromatic forward projection is performed using material-specific images, which are obtained using dual-energy material decomposition. The algorithm was evaluated using numerical simulations, experimental phantom data, and realistic patient data. The results show that the estimated spectrum matches the reference spectrum quite well and the method is robust. Extensive studies suggest that the method provides an accurate estimate of the CT spectrum without dedicated physical phantom and prolonged workflow. This paper may be attractive for CT dose calculation, artifacts reduction, polychromatic image reconstruction, and other spectrum-involved CT applications.

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

[2]  Li Zhang,et al.  A weighted polynomial based material decomposition method for spectral x-ray CT imaging. , 2016, Physics in medicine and biology.

[3]  G. Hernández,et al.  A model of tungsten anode x-ray spectra. , 2016, Medical physics.

[4]  F Verhaegen,et al.  SpekCalc: a program to calculate photon spectra from tungsten anode x-ray tubes , 2009, Physics in medicine and biology.

[5]  E. Wilkinson,et al.  A comparison of mammography spectral measurements with spectra produced using several different mathematical models. , 2001, Physics in medicine and biology.

[6]  J. Fessler,et al.  Modelling the physics in the iterative reconstruction for transmission computed tomography , 2013, Physics in medicine and biology.

[7]  Wei Zhao,et al.  Segmentation-free x-ray energy spectrum estimation for computed tomography , 2016, SPIE Medical Imaging.

[8]  Francesc Salvat,et al.  Monte Carlo simulation of x-ray spectra generated by kilo-electron-volt electrons , 2003 .

[9]  Comparison of simulated and measured spectra of an industrial 450 kV X-ray tube , 2007 .

[10]  Monte carlo simulation of the compton scattering technique applied to characterize diagnostic x-ray spectra. , 2004, Medical physics.

[11]  J J DeMarco,et al.  A Monte Carlo based method to estimate radiation dose from multidetector CT (MDCT): cylindrical and anthropomorphic phantoms. , 2005, Physics in medicine and biology.

[12]  Katsuyuki Taguchi,et al.  Enabling photon counting clinical X-ray CT , 2009, 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC).

[13]  Li Zhang,et al.  X-ray spectrum estimation from transmission measurements using the expectation maximization method , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[14]  L K Wagner,et al.  A Laplace transform pair model for spectral reconstruction. , 1982, Medical physics.

[15]  C. McCollough,et al.  CT scanner x-ray spectrum estimation from transmission measurements. , 2011, Medical physics.

[16]  J H Siewerdsen,et al.  Spektr: a computational tool for x-ray spectral analysis and imaging system optimization. , 2004, Medical physics.

[17]  Jun Zhao,et al.  United Iterative Reconstruction for Spectral Computed Tomography , 2015, IEEE Transactions on Medical Imaging.

[18]  Xuanqin Mou,et al.  A statistical iterative reconstruction framework for dual energy computed tomography without knowing tube spectrum , 2016, Optical Engineering + Applications.

[19]  L. Garnero,et al.  X-ray dual-energy calibration based on estimated spectral properties of the experimental system , 1998 .

[20]  Frank Verhaegen,et al.  Monte Carlo simulation of a computed tomography x-ray tube , 2007, Physics in medicine and biology.

[21]  Emil Y. Sidky,et al.  X-ray spectrum estimation from transmission measurements by an exponential of a polynomial model , 2016, SPIE Medical Imaging.

[22]  Jeffrey A. Fessler,et al.  Statistical image reconstruction for polyenergetic X-ray computed tomography , 2002, IEEE Transactions on Medical Imaging.

[23]  Raphael Thierry,et al.  Monte Carlo simulations of a high-resolution X-ray CT system for industrial applications , 2007 .

[24]  Koji Maeda,et al.  Compton-scattering measurement of diagnostic x-ray spectrum using high-resolution Schottky CdTe detector. , 2005, Medical physics.

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

[26]  E R van der Graaf,et al.  Spectra of clinical CT scanners using a portable Compton spectrometer. , 2015, Medical physics.

[27]  Lei Zhu,et al.  Scatter correction for cone-beam CT in radiation therapy. , 2009 .

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

[29]  Jeffrey A. Fessler,et al.  Multi-Material Decomposition Using Statistical Image Reconstruction for Spectral CT , 2014, IEEE Transactions on Medical Imaging.

[30]  Jyh-Cheng Chen,et al.  A Single Scatter Model for X-ray CT Energy Spectrum Estimation and Polychromatic Reconstruction , 2015, IEEE Transactions on Medical Imaging.

[31]  L. Xing,et al.  Using edge-preserving algorithm with non-local mean for significantly improved image-domain material decomposition in dual-energy CT , 2016, Physics in medicine and biology.

[32]  Timo Berkus,et al.  Empirical dual energy calibration (EDEC) for cone-beam computed tomography. , 2007 .

[33]  Wei Zhao,et al.  An indirect transmission measurement-based spectrum estimation method for computed tomography , 2015, Physics in medicine and biology.

[34]  A. Sisniega,et al.  Technical Note: spektr 3.0—A computational tool for x-ray spectrum modeling and analysis , 2016, Medical physics.

[35]  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.

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

[37]  Geant4 simulation of a filtered X-ray source for radiation damage studies , 2012, 1205.1973.

[38]  Andrew M Hernandez,et al.  Tungsten anode spectral model using interpolating cubic splines: unfiltered x-ray spectra from 20 kV to 640 kV. , 2014, Medical physics.

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

[40]  S. Incerti,et al.  Geant4 developments and applications , 2006, IEEE Transactions on Nuclear Science.

[41]  Energy spectrum extraction and optimal imaging via dual-energy material decomposition , 2015, 2015 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC).

[42]  Paulo R. S. Mendonça,et al.  A Flexible Method for Multi-Material Decomposition of Dual-Energy CT Images , 2014, IEEE Transactions on Medical Imaging.

[43]  Uwe Oelfke,et al.  Imaging properties of small-pixel spectroscopic x-ray detectors based on cadmium telluride sensors , 2012, Physics in medicine and biology.

[44]  Ehsan Samei,et al.  An angle-dependent estimation of CT x-ray spectrum from rotational transmission measurements. , 2014, Medical physics.

[45]  H Zaidi,et al.  Monte carlo simulation of x-ray spectra in diagnostic radiology and mammography using MCNP4C. , 2004, Physics in medicine and biology.

[46]  Iwan Kawrakow,et al.  Efficient x-ray tube simulations. , 2006, Medical physics.

[47]  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 .

[48]  K. Taguchi,et al.  Vision 20/20: Single photon counting x-ray detectors in medical imaging. , 2013, Medical physics.

[49]  Lei Zhu,et al.  Iterative image-domain decomposition for dual-energy CT. , 2014, Medical physics.

[50]  Peng Zhang,et al.  An Extended Algebraic Reconstruction Technique (E-ART) for Dual Spectral CT , 2015, IEEE Transactions on Medical Imaging.

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

[52]  A microCT X-ray head model for spectra generation with Monte Carlo simulations , 2006 .

[53]  Marc Kachelrieß,et al.  Image-based dual energy CT using optimized precorrection functions: A practical new approach of material decomposition in image domain. , 2009, Medical physics.