Assessment of different computational models for generation of x-ray spectra in diagnostic radiology and mammography.

Different computational methods based on empirical or semi-empirical models and sophisticated Monte Carlo calculations have been proposed for prediction of x-ray spectra both in diagnostic radiology and mammography. In this work, the x-ray spectra predicted by various computational models used in the diagnostic radiology and mammography energy range have been assessed by comparison with measured spectra and their effect on the calculation of absorbed dose and effective dose (ED) imparted to the adult ORNL hermaphroditic phantom quantified. This includes empirical models (TASMIP and MASMIP), semi-empirical models (X-rayb&m, X-raytbc, XCOMP, IPEM, Tucker et al., and Blough et al.), and Monte Carlo modeling (EGS4, ITS3.0, and MCNP4C). As part of the comparative assessment, the K x-ray yield, transmission curves, and half value layers (HVLs) have been calculated for the spectra generated with all computational models at different tube voltages. The measured x-ray spectra agreed well with the generated spectra when using X-raytbc and IPEM in diagnostic radiology and mammography energy ranges, respectively. Despite the systematic differences between the simulated and reference spectra for some models, the student's t-test statistical analysis showed there is no statistically significant difference between measured and generated spectra for all computational models investigated in this study. The MCNP4C-based Monte Carlo calculations showed there is no discernable discrepancy in the calculation of absorbed dose and ED in the adult ORNL hermaphroditic phantom when using different computational models for generating the x-ray spectra. Nevertheless, given the limited flexibility of the empirical and semi-empirical models, the spectra obtained through Monte Carlo modeling offer several advantages by providing detailed information about the interactions in the target and filters, which is relevant for the design of new target and filter combinations and optimization of radiological imaging protocols.

[1]  G. D. Valdez,et al.  ITS: the integrated TIGER series of electron/photon transport codes-Version 3.0 , 1991, Conference Record of the 1991 IEEE Nuclear Science Symposium and Medical Imaging Conference.

[2]  H. Kramers,et al.  XCIII. On the theory of X-ray absorption and of the continuous X-ray spectrum , 1923 .

[3]  J. H. Hubbell,et al.  XCOM: Photon cross sections on a personal computer , 1987 .

[4]  K. F. Eckerman,et al.  Specific absorbed fractions of energy at various ages from internal photon sources: 6, Newborn , 1987 .

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

[6]  V. Cosslett,et al.  Measurements of K, L and M shell X-ray production efficiencies , 1968 .

[7]  B. Archer,et al.  Determination of diagnostic x-ray spectra with characteristic radiation using attenuation analysis. , 1988, Medical physics.

[8]  W. H. Payne,et al.  Calculated mammographic spectra confirmed with attenuation curves for molybdenum, rhodium, and tungsten targets. , 1998, Medical physics.

[9]  J. Boone.  X-ray spectral reconstruction from attenuation data using neural networks. , 1990, Medical physics.

[10]  G T Barnes,et al.  Molybdenum target x-ray spectra: a semiempirical model. , 1991, Medical physics.

[11]  L J Schreiner,et al.  A comparison of semiempirical models for generating tungsten target x-ray spectra. , 1992, Medical physics.

[12]  J M Boone,et al.  The three parameter equivalent spectra as an index of beam quality. , 1988, Medical physics.

[13]  J. F. Briesmeister.  MCNP-A General Monte Carlo N-Particle Transport Code , 1993 .

[14]  P. Huang,et al.  Reconstruction of diagnostic x-ray spectra by numerical analysis of transmission data. , 1986, Medical physics.

[15]  C S Kwok,et al.  Monte Carlo simulation of x-ray spectra in mammography. , 2000, Physics in medicine and biology.

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

[17]  T. R. Fewell,et al.  Molybdenum, rhodium, and tungsten anode spectral models using interpolating polynomials with application to mammography. , 1997, Medical physics.

[18]  F Verhaegen,et al.  Monte Carlo modelling of radiotherapy kV x-ray units. , 1999, Physics in medicine and biology.

[19]  T. Fewell,et al.  Laplace reconstruction of experimental diagnostic x-ray spectra. , 1988, Medical physics.

[20]  P. François,et al.  Simulation of x-ray spectral reconstruction from transmission data by direct resolution of the numeric system AF = T. , 1993, Medical physics.

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

[22]  J. Pattison,et al.  The effect on dose to computed tomography phantoms of varying the theoretical x-ray spectrum: a comparison of four diagnostic x-ray spectrum calculating codes. , 1998, Medical physics.

[23]  Ludwik Silberstein,et al.  Determination of the Spectral Composition of X-ray Radiation from Filtration Data* , 1932 .

[24]  Design considerations for transmission x-ray tubes operated at diagnostic energies , 2000 .

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

[26]  Andreas Höfer,et al.  Ein Programm für die Berechnung von diagnostischen Röntgenspektren , 1985 .

[27]  B W Soole,et al.  A method of X-ray attenuation analysis for approximating the intensity distribution at its point of origin of bremsstrahlung excited in a thick target by incident electrons of constant medium energy. , 1976, Physics in medicine and biology.

[28]  D. Rogers,et al.  EGS4 code system , 1985 .

[29]  B. Wall,et al.  A national survey of doses to patients undergoing a selection of routine x-ray examinations in English hospitals. NRPB-R200 , 1986 .

[30]  S. Evans.  Catalogue of Diagnostic X-Ray Spectra and Other Data , 1998 .

[31]  R N Kulkarni,et al.  Monte Carlo calculations of mammographic x-ray spectra. , 1984, Physics in medicine and biology.

[32]  T R Fewell,et al.  Photon energy distribution of some typical diagnostic x-ray beams. , 1977, Medical physics.

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

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

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

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

[37]  D R Dance,et al.  Influence of anode/filter material and tube potential on contrast, signal-to-noise ratio and average absorbed dose in mammography: a Monte Carlo study. , 2000, The British journal of radiology.