An empirical model of diagnostic x-ray attenuation under narrow-beam geometry.

PURPOSE The purpose of this study was to develop and validate a mathematical model to describe narrow-beam attenuation of kilovoltage x-ray beams for the intended applications of half-value layer (HVL) and quarter-value layer (QVL) estimations, patient organ shielding, and computer modeling. METHODS An empirical model, which uses the Lambert W function and represents a generalized Lambert-Beer law, was developed. To validate this model, transmission of diagnostic energy x-ray beams was measured over a wide range of attenuator thicknesses [0.49-33.03 mm Al on a computed tomography (CT) scanner, 0.09-1.93 mm Al on two mammography systems, and 0.1-0.45 mm Cu and 0.49-14.87 mm Al using general radiography]. Exposure measurements were acquired under narrow-beam geometry using standard methods, including the appropriate ionization chamber, for each radiographic system. Nonlinear regression was used to find the best-fit curve of the proposed Lambert W model to each measured transmission versus attenuator thickness data set. In addition to validating the Lambert W model, we also assessed the performance of two-point Lambert W interpolation compared to traditional methods for estimating the HVL and QVL [i.e., semi-logarithmic (exponential) and linear interpolation]. RESULTS The Lambert W model was validated for modeling attenuation versus attenuator thickness with respect to the data collected in this study (R2 > 0.99). Furthermore, Lambert W interpolation was more accurate and less sensitive to the choice of interpolation points used to estimate the HVL and/or QVL than the traditional methods of semilogarithmic and linear interpolation. CONCLUSIONS The proposed Lambert W model accurately describes attenuation of both monoenergetic radiation and (kilovoltage) polyenergetic beams (under narrow-beam geometry).

[1]  R. Sloboda,et al.  Linear accelerator photon beam quality at off-axis points. , 1997, Medical physics.

[2]  B. Archer,et al.  On the measurement of half-value layer in film-screen mammography. , 1990, Medical physics.

[3]  Cynthia H McCollough,et al.  A method to generate equivalent energy spectra and filtration models based on measurement for multidetector CT Monte Carlo dosimetry simulations. , 2009, Medical physics.

[4]  H. Akaike.  A new look at the statistical model identification , 1974 .

[5]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[6]  Cynthia H McCollough,et al.  Measurement of half-value layer in x-ray CT: a comparison of two noninvasive techniques. , 2000 .

[7]  S C Bushong,et al.  Diagnostic X-ray shielding design based on an empirical model of photon attenuation. , 1983, Health physics.

[8]  M. McKetty,et al.  The AAPM/RSNA physics tutorial for residents. X-ray attenuation. , 1998, Radiographics : a review publication of the Radiological Society of North America, Inc.

[9]  B. Bjärngard,et al.  Attenuation in high-energy x-ray beams. , 1994, Medical physics.

[10]  Gaston H. Gonnet,et al.  On the LambertW function , 1996, Adv. Comput. Math..

[11]  Robert M Corless,et al.  Some applications of the Lambert W  function to physics , 2000, Canadian Journal of Physics.

[12]  C Kleinschmidt,et al.  Analytical considerations of beam hardening in medical accelerator photon spectra. , 1999, Medical physics.

[13]  M. McNitt-Gray,et al.  Precision of dosimetry-related measurements obtained on current multidetector computed tomography scanners. , 2010, Medical physics.