A CT-based analytical dose calculation method for HDR 192Ir brachytherapy.

PURPOSE This article presents an analytical dose calculation method for high-dose-rate 192Ir brachytherapy, taking into account the effects of inhomogeneities and reduced photon backscatter near the skin. The adequacy of the Task Group 43 (TG-43) two-dimensional formalism for treatment planning is also assessed. METHODS The proposed method uses material composition and density data derived from computed tomography images. The primary and scatter dose distributions for each dwell position are calculated first as if the patient is an infinite water phantom. This is done using either TG-43 or a database of Monte Carlo (MC) dose distributions. The latter can be used to account for the effects of shielding in water. Subsequently, corrections for photon attenuation, scatter, and spectral variations along medium- or low-Z inhomogeneities are made according to the radiological paths determined by ray tracing. The scatter dose is then scaled by a correction factor that depends on the distances between the point of interest, the body contour, and the source position. Dose calculations are done for phantoms with tissue and lead inserts, as well as patient plans for head-and-neck, esophagus, and MammoSite balloon breast brachytherapy treatments. Gamma indices are evaluated using a dose-difference criterion of 3% and a distance-to-agreement criterion of 2 mm. PTRAN_CT MC calculations are used as the reference dose distributions. RESULTS For the phantom with tissue and lead inserts, the percentages of the voxels of interest passing the gamma criteria (Pgamma > or = 1) are 100% for the analytical calculation and 91% for TG-43. For the breast patient plan, TG-43 overestimates the target volume receiving the prescribed dose by 4% and the dose to the hottest 0.1 cm3 of the skin by 9%, whereas the analytical and MC results agree within 0.4%. Pgamma > or = 1 are 100% and 48% for the analytical and TG-43 calculations, respectively. For the head-and-neck and esophagus patient plans, Pgamma > or = 1 are > or = 99% for both calculation methods. CONCLUSIONS A correction-based dose calculation method has been validated for HDR 192Ir brachytherapy. Its high calculation efficiency makes it feasible for use in treatment planning. Because tissue inhomogeneity effects are small and primary dose predominates in the near-source region, TG-43 is adequate for target dose estimation provided shielding and contrast solution are not used.

[1]  J F Williamson,et al.  On the validity of the superposition principle in dose calculations for intracavitary implants with shielded vaginal colpostats. , 2001, Medical physics.

[2]  J F Williamson,et al.  Monte Carlo evaluation of kerma at a point for photon transport problems. , 1987, Medical physics.

[3]  J F Williamson,et al.  Comparison of measured and calculated dose rates in water near I-125 and Ir-192 seeds. , 1991, Medical physics.

[4]  R. Sloboda,et al.  Brachytherapy scatter dose calculation in heterogeneous media: II. Empirical formulation for the multiple-scatter contribution , 2007, Physics in medicine and biology.

[5]  Brachytherapy scatter dose calculation in heterogeneous media: I. A microbeam ray-tracing method for the single-scatter contribution. , 2007, Physics in medicine and biology.

[6]  X. Li,et al.  Dose characterization in the near-source region for two high dose rate brachytherapy sources. , 2002, Medical Physics (Lancaster).

[7]  J. Williamson,et al.  Monte Carlo-aided dosimetry of a new high dose-rate brachytherapy source. , 1998, Medical physics.

[8]  D Baltas,et al.  Evaluation of a TG-43 compliant analytical dosimetry model in clinical 192Ir HDR brachytherapy treatment planning and assessment of the significance of source position and catheter reconstruction uncertainties. , 2004, Physics in medicine and biology.

[9]  Dose perturbation induced by radiographic contrast inside brachytherapy balloon applicators. , 2004, Medical physics.

[10]  Anders Ahnesjö,et al.  Accounting for high Z shields in brachytherapy using collapsed cone superposition for scatter dose calculation. , 2003, Medical physics.

[11]  Frank Verhaegen,et al.  Patient-specific Monte Carlo dose calculations for high-dose-rate endorectal brachytherapy with shielded intracavitary applicator. , 2008, International journal of radiation oncology, biology, physics.

[12]  Measurement and calculation of heterogeneity correction factors for an Ir-192 high dose-rate brachytherapy source behind tungsten alloy and steel shields. , 1996, Medical physics.

[13]  305C Physics and quality assurance for brachytherapy — part II: Low dose rate and pulsed dose rate , 1996 .

[14]  F Verhaegen,et al.  BrachyGUI: an adjunct to an accelerated Monte Carlo photon transport code for patient-specific brachytherapy dose calculations and analysis , 2008 .

[15]  J. Williamson,et al.  Update of AAPM Task Group No. 43 Report: A revised AAPM protocol for brachytherapy dose calculations. , 2004 .

[16]  J F Williamson,et al.  Dose calculations about shielded gynecological colpostats. , 1990, International journal of radiation oncology, biology, physics.

[17]  D. Baltas,et al.  An analytical dosimetry model as a step towards accounting for inhomogeneities and bounded geometries in 192Ir brachytherapy treatment planning. , 2003, Physics in medicine and biology.

[18]  H. Meertens,et al.  Screens in ovoids of a Selectron cervix applicator. , 1985, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[19]  Anders Ahnesjö,et al.  Optimization of the computational efficiency of a 3D, collapsed cone dose calculation algorithm for brachytherapy. , 2008, Medical physics.

[20]  E. Poon,et al.  Development of a scatter correction technique and its application to HDR I192r multicatheter breast brachytherapy. , 2009, Medical physics.

[21]  R. Nath,et al.  Tissue inhomogeneity correction for brachytherapy sources in a heterogeneous phantom with cylindrical symmetry. , 1992, Medical physics.

[22]  Luc Beaulieu,et al.  The evolution of brachytherapy treatment planning. , 2009, Medical physics.

[23]  R. Taylor,et al.  Benchmarking brachydose: Voxel based EGSnrc Monte Carlo calculations of TG-43 dosimetry parameters. , 2007, Medical physics.

[24]  Jan-Jakob Sonke,et al.  A fast algorithm for gamma evaluation in 3D. , 2007, Medical physics.

[25]  J F Williamson,et al.  Dosimetric modeling of the microselectron high-dose rate 192Ir source by the multigroup discrete ordinates method. , 2000, Medical physics.

[26]  Mark J Rivard,et al.  Approaches to calculating AAPM TG-43 brachytherapy dosimetry parameters for 137Cs, 125I, 192Ir, 103Pd, and 169Yb sources. , 2006, Medical physics.

[27]  D Baltas,et al.  A monte carlo dosimetry study of vaginal 192Ir brachytherapy applications with a shielded cylindrical applicator set. , 2004, Medical physics.

[28]  Firas Mourtada,et al.  Comparison of a finite-element multigroup discrete-ordinates code with Monte Carlo for radiotherapy calculations , 2006, Physics in medicine and biology.

[29]  F. Vicini,et al.  Dosimetric characteristics of the MammoSite RTS, a new breast brachytherapy applicator. , 2002, International journal of radiation oncology, biology, physics.

[30]  Anders Ahnesjö,et al.  Brachytherapy source characterization for improved dose calculations using primary and scatter dose separation. , 2005, Medical physics.

[31]  J. Williamson.  The Sievert integral revisited: evaluation and extension to 125I, 169Yb, and 192Ir brachytherapy sources. , 1996, International journal of radiation oncology, biology, physics.

[32]  Frank Verhaegen,et al.  Tissue segmentation in Monte Carlo treatment planning: a simulation study using dual-energy CT images. , 2008, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[33]  J. Williamson,et al.  Two-dimensional scatter integration method for brachytherapy dose calculations in 3D geometry. , 1997, Physics in medicine and biology.

[34]  K R Russell,et al.  Dose calculation in brachytherapy for a source using a primary and scatter dose separation technique , 1996, Physics in medicine and biology.

[35]  M. Altschuler,et al.  Rapid three-dimensional treatment planning: I. Ray-tracing approach to primary component dose calculations. , 1987, Physics in medicine and biology.

[36]  Siyong Kim,et al.  Dose errors due to inhomogeneities in balloon catheter brachytherapy for breast cancer. , 2004, International journal of radiation oncology, biology, physics.

[37]  D. Baltas,et al.  On the use of high dose rate Ir192 and Yb169 sources with the MammoSite® radiation therapy system. , 2007, Medical physics.

[38]  J. Williamson,et al.  Analytical approach to heterogeneity correction factor calculation for brachytherapy. , 1998, Medical physics.

[39]  Proceedings of the 38th Annual ASTRO Meeting305C Physics and quality assurance for brachytherapy — Part II: Low dose rate and pulsed dose rate , 1996 .

[40]  J. H. Hubbell,et al.  EPDL97: the evaluated photo data library `97 version , 1997 .

[41]  R. Siddon.  Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[42]  Omar Chibani,et al.  MCPI©: A sub-minute Monte Carlo dose calculation engine for prostate implants. , 2005, Medical physics.

[43]  A Ahnesjö,et al.  The collapsed cone superposition algorithm applied to scatter dose calculations in brachytherapy. , 2000, Medical physics.

[44]  D Baltas,et al.  The effect of patient inhomogeneities in oesophageal 192Ir HDR brachytherapy: a Monte Carlo and analytical dosimetry study. , 2004, Physics in medicine and biology.

[45]  Luc Beaulieu,et al.  Correction of CT artifacts and its influence on Monte Carlo dose calculations. , 2007, Medical physics.

[46]  J. Siebers,et al.  MO‐E‐T‐618‐01: An Integrated CT‐Based Monte Carlo Dose‐Evaluation System for Brachytherapy and Its Application to Permanent Prostate Implant Postprocedure Dosimetric Analysis , 2005 .

[47]  Michael J Price,et al.  Optimization of deterministic transport parameters for the calculation of the dose distribution around a high dose-rate 192Ir brachytherapy source. , 2008, Medical physics.