Segmentation-free x-ray energy spectrum estimation for computed tomography

X-ray energy spectrum plays an essential role in imaging and related tasks. Due to the high photon flux of clinical CT scanners, most of the spectrum estimation methods are indirect and are usually suffered from various limitations. The recently proposed indirect transmission measurement-based method requires at least the segmentation of one material, which is insufficient for CT images of highly noisy and with artifacts. To combat for the bottleneck of spectrum estimation using segmented CT images, in this study, we develop a segmentation-free indirect transmission measurement based energy spectrum estimation method using dual-energy material decomposition. The general principle of the method is to compare polychromatic forward projection with raw projection to calibrate a set of unknown weights which are used to express the unknown spectrum together with a set of model spectra. After applying dual-energy material decomposition using high-and low-energy raw projection data, polychromatic forward projection is conducted on material-specific images. The unknown weights are then iteratively updated to minimize the difference between the raw projection and estimated projection. Both numerical simulations and experimental head phantom are used to evaluate the proposed method. The results indicate that the method provides accurate estimate of the spectrum and it may be attractive for dose calculations, artifacts correction and other clinical applications.

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

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

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

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

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

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

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

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

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

[10]  Kai Niu,et al.  A patient-specific scatter artifacts correction method , 2014, Medical Imaging.

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

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

[13]  Wei Zhao,et al.  Patient-specific scatter correction for flat-panel detector-based cone-beam CT imaging , 2015, Physics in medicine and biology.

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

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