

ORIGINAL ARTICLE 

Year : 2008  Volume
: 4
 Issue : 4  Page : 173177 

A dose verification method for highdoserate brachytherapy treatment plans
Rajesh Kumar^{1}, SD Sharma^{1}, C Vijaykumar^{2}, Sudesh Deshpande^{3}, PK Sharma^{4}, S Vandana^{1}, A Philomena^{1}, Ravi H Chilkulwar^{1}
^{1} Radiological Physics and Advisory Division, Bhabha Atomic Research Centre, CT and CRS, Anushaktinagar, Mumbai400 094, India ^{2} Department of Radiology, Armed Force Medical College, Pune, India ^{3} Department of Radiology, PD Hinduja National Hospital and Research Centre, Mahim, Mumbai, India ^{4} Department of Radiology, Tata Memorial Hospital, Parel, Mumbai, India
Date of Web Publication  1Dec2008 
Correspondence Address: Rajesh Kumar Radiological Physics and Advisory Division, Bhabha Atomic Research Centre, CT and CRS, Anushaktinagar, Mumbai400 094 India
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09731482.44288
Aim: To evolve a fast dose verification method for highdoserate (HDR) brachytherapy treatment plans and to demonstrate its applicability in different clinical cases. Materials and Methods: We developed a software tool in VC++ for the Varisource HDR unit for HDR dosimetry plan verification using TG43 parameters. HDR treatment dosimetry of a number clinical cases using Varisource was verified by comparison with the treatment planning system (TPS). Results: A number of different types of clinical cases treated by Varisource were evaluated. TPS calculated dose values and verification code calculated dose values were found to agree to within 3% for most of the dose calculation points. Conclusions: We have validated with clinical cases a fast and independent dose verification method of the dosimetry at selected points for HDR brachytherapy treatments plan using TG43 parameters. This can be used for the verification of the TPS calculated dose at various points. The code is written to work with Varisource, but it can conceivably be modified for other sources also by using the fitted constant of the respective source. Keywords: Brachytherapy, dose verification, dosimetry, TG43 formalism, treatment planning
How to cite this article: Kumar R, Sharma S D, Vijaykumar C, Deshpande S, Sharma P K, Vandana S, Philomena A, Chilkulwar RH. A dose verification method for highdoserate brachytherapy treatment plans. J Can Res Ther 2008;4:1737 
How to cite this URL: Kumar R, Sharma S D, Vijaykumar C, Deshpande S, Sharma P K, Vandana S, Philomena A, Chilkulwar RH. A dose verification method for highdoserate brachytherapy treatment plans. J Can Res Ther [serial online] 2008 [cited 2020 Jan 29];4:1737. Available from: http://www.cancerjournal.net/text.asp?2008/4/4/173/44288 
> Introduction   
Highdoserate (HDR) brachytherapy has proven to be an effective treatment in the definitive management of different types of cancers and is a common treatment modality in most radiotherapy clinics. In HDR brachytherapy treatment planning, the identification of applicators is done using simulator films or 3D images such as CT, and the dwell times along these catheters are then optimized in order to deliver prescribed doses at one or more anatomical points while simultaneously satisfying various constraints. The importance of independent verification of the dosimetry prior to the treatment delivery has been recognized worldwide and is also specified in the AERB safety code as well as in the guidelines of other international regulatory agencies.^{ [1]} The planning system provides many different possibilities for dose optimization; for example, it allows alteration of the optimized plan by either pulling the isodose lines with the mouseor by directly modifying the dwell times. The complex nature of the treatment planning process means that many approaches can be used to arrive at a final plan, making it possible that a software bug that escaped discovery during commissioning may introduce error in the dosimetry if the planning steps are completed in a different order than usual. This makes a second check of the dose especially useful. Treatment planning systems using a stepping source, which can be interfaced with multimodality images (CT/MRI/ultrasound) and sophisticated dose optimization software, enable the planner to maximize the dose uniformity, while minimizing the implant volume needed to cover the target volume adequately and, at the same time, reduces the dose to the organs at risk. Such flexibility creates a challenge for the verification of the optimized calculations with practical manual calculation techniques that take only a few minutes but at the same time give a high probability of detecting significant errors. Verification of the dose calculation on a pretreatment basis, however, ensures that the correct source is being used; that the source data has not been modified and correct activity, treatment date, and decay are used; and that any bugs in the planning software did not affect the dose calculation. HDR treatment plans involve a number of dwell positions, which makes manual calculations of dose distribution, and hence of the treatment time, impractical. Although TPS facilitates the determination of the dose optimization and treatment time calculation, verification of the TPS calculated dose can pose a challenge. Most of the dose verification methods previously published are specific to particular types of clinical treatments. In this paper we describe an independent method for use in verifying HDR dosimetry of different types of clinical cases based on the AAPM TG43 formalism^{ [2] } and Cartesian coordinates obtained from the TPS.
> Materials and Methods   
The verification software code was written in VC++ for the Varisource (new design) HDR brachytherapy unit. In order to verify errors which may occur between the planning and delivery stages, the verification code uses geometrical information such as dwell time (in seconds), Cartesian coordinates (in centimeters) of dwell position, and calculation point coordinates (in centimeters) directly from the TPS, as well as air kerma strength (AKS; cGy/cm^{ 2} /h) at treatment date, which has to be entered independently by the user. The following AAPM TG43 formalism^{ [2]} was used for dose calculation in the verification code:
where S_{k} is the AKS of the source, is the dose rate constant, G ( r ,θ) is the geometry factor, g ( r ) is the radial dose function, F ( r ,θ) is the anisotropy function, r is the distance between the dose point and the centre of the source, t is dwell time in second, θ is the angle subtended by the central axis of the source and line connecting the centre of the source and the dose point, and r_{} 0 and θ_{0} are reference parameters taken to be 1 cm and 90°, respectively [Figure 1]. Analytically fitted formulae based on the Monte Carlo data of Angelopoulos et al.^{} [3,5] for the source were used to find out g ( r ) and F ( r ,θ) for the desired point of dose calculation.
For the radial dose function, the general functional form used in this study was^{ [3]} as follows:
where h, i, j, and k are fitted constants from Monte Carlo data of the respective HDR source.
For the geometry function the following expression was used in this study^{ [4]} :
where L is the active length of source.
The following expression was used to calculate anisotropy function in this study:^{ [3]}
where
are fitted constants from Monte Carlo data of the respective HDR source. Values of fitted constants for Varisource are given in [Table 1] and [Table 2].
For Cartesian dwell coordinates (Xi,Yi,Zi) and corresponding dwell time, the dose at any point (Xr, Yr, Zr) can be calculated using the following formula:
where D(r,θ) can be calculated from Equation 1, while r and θ can be calculated as:
and
where r_{ 1} and r_{ 2} are vectors, as shown in [Figure 2].
Dose rate constant reported by Angelopoulos et al .^{ [5]} 1.101 for new Varisource was used in this code during the evaluation. A screenshot of the software can be seen in [Figure 3]. Input data required by the software code can be taken from the plan report generated by the TPS (such as dwell position and calculation point coordinate, dwell time at each dwell position, AKS on the date of treatment, and TPS calculated dose at the dose calculation point). It subsequently calculates dose and associated percentage error for the points of interest, which is used as a formal verification of the treatment dosimetry. We evaluated this software on a number of different types of clinical cases treated with Varisource HDR brachytherapy. Calculation points were selected randomly from the list of the dose points entered during the treatment planning and compared with the TPS calculated dose at that point. It is notable that no correction was applied for attenuation of radiation by the applicator in this code. [Figure 4] shows the visual representation of the geometrical position of the source as well as the dose calculation point of the different types of clinical cases for which the code was evaluated.
> Results   
The code was tested for a number of clinical cases at various positions against the TPS calculated dose values (BrachyVision, Varian Medical System, USA) for the Varisource. [Table 3a] shows the TPS and verification code calculated dose values and their percentage deviation for intracavitary brachytherapy treatment plans. The data in this table indicate that the percentage variation between the TPS calculated dose and verification code calculated dose is of the order of 2%. [Table 3b] shows the TPS and verification code calculated dose values and their percentage deviation for central vaginal source (CVS) brachytherapy treatment plans. Percentage deviations for CVS quoted in this table are within 1% for most of the points, except for one point where the deviation is about 2.4%. [Table 3c] shows TPS and verification code calculated dose values and their percentage deviation for intraluminal brachytherapy treatment plans. Percentage deviations for intraluminal brachytherapy quoted in this table are within 1.5% for most of the calculation points, except at few points where the deviation is up to 2.3%. [Table 3ac] indicate that the variation in the TPS calculated dose and verification code calculated dose values is less than 1.5% at most of the evaluation points. Based on these observations, one can conclude that the TPS calculated dose values and verification code calculated dose values for the different clinical cases considered here are in good agreement (to within 3%) at most of the points of clinical importance.
> Discussion   
Pretreatment dose verification is an important aspect of patientspecific quality assurance. A simple, fast, and accurate method of dose verification is required to fulfil this requirement. The AAPM TG43 formalismbased verification code described here is a quick and simple method to check the dose calculation by the TPS. It can be used for verifying the accuracy of dose calculations during the commissioning of the TPS and, subsequently, for periodical quality control checks. Checking the dose calculation for each and every patient ensures the accuracy of dose delivery and the safety of the patients during brachytherapy treatments. Though the method is tested for the Varisource HDR brachytherapy source, it is a general verification method which can be applied for verifying the treatment plans of other HDR brachytherapy sources. Even though the accuracy of dose calculation of TPS is fully verified during commissioning and, subsequently, during periodic quality control tests, checking the dosimetry prior to each treatment ensures that the correct source data is being used and that the source data has not been modified by software bugs either deliberately or incidentally. It is important to note here that the basic source data required for dose calculation by the verification code is entered into the computer independently and has no direct link with the data presented by the commercial planning system. This way the present code provides confidence to the user against any malfunction of the commercial TPS as the decay correction of the source is accounted for by the planning software automatically.
> Conclusion   
The code independently verifies the dose calculated by the TPS at selected points for HDR brachytherapy using AAPM TG43 parameters. The code is written to work with Varisource, but it can conceivably be modified for other sources also by using the fitted constant of the respective source.
> Acknowledgment   
We are thankful to Dr Y.S. Mayya, Head, RP and AD and Shri H.S. Kushwaha, Director, HS and EG, BARC for their encouragement and interest in this work.
> References   
1.  Atomic Energy Regulatory Board (AERB). Safety Code for Medical Applications of Ionizing Radiation, AERB SCMED 1 and 3, 2007. (Yet to be adopted) 
2.  Nath R, Anderson LL, Luxton G, Weaver KA, Williamson JF, Meigooni AS. Dosimetry of interstitial brachytherapy sources: Recommendations of the AAPM Radiation Therapy Committee Task Group No. 43. Med Phys 1995;22:20934. 
3.  Lliso F, PιrezCalatayud J, Carmona V, Ballester F, Puchades V, Granero D. Technical note: Fitted dosimetric parameters of high doserate 192Ir sources according to the AAPM TG43 formalism. Med Phys 2003;30:6514. 
4.  Rivard MJ. Refinements to the geometry factor used in the AAPM Task Group Report No 43 necessary for brachytherapy dosimetry calculations. Med Phys 1999;26:244550. 
5.  Angelopoulos A, Baras P, Sakelliou L, Karaiskos P, Sandilos P. Monte Carlo dosimetry of a new 192Ir high dose rate brachytherapy source. Med Phys 2000;27:25217. 
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2], [Table 3a], [Table 3b], [Table 3c], [Table 3ac]
This article has been cited by  1 
A dose verification tool for highdoserate interstitial brachytherapy treatment planning in accelerated partial breast irradiation 

 Mohamad Feras Marqa,JeanMichel Caudrelier,Nacim Betrouni   Brachytherapy. 2012; 11(5): 359   [Pubmed]  [DOI]   2 
The usefulness of an independent patientspecific treatment planning verification method using a benchmark plan in highdoserate intracavitary brachytherapy for carcinoma of the uterine cervix 

 Y. Takahashi,M. Koizumi,I. Sumida,F. Isohashi,T. Ogata,Y. Akino,Y. Yoshioka,S. Maruoka,S. Inoue,K. Konishi,K. Ogawa   Journal of Radiation Research. 2012; 53(6): 936   [Pubmed]  [DOI]   3 
The usefulness of an independent patientspecific treatment planning verification method using a benchmark plan in highdoserate intracavitary brachytherapy for carcinoma of the uterine cervix 

 Takahashi, Y. and Koizumi, M. and Sumida, I. and Isohashi, F. and Ogata, T. and Akino, Y. and Yoshioka, Y. and Maruoka, S. and Inoue, S. and Konishi, K. and Ogawa, K.   Journal of Radiation Research. 2012; 53(6): 936944   [Pubmed]   4 
A dose verification tool for highdoserate interstitial brachytherapy treatment planning in accelerated partial breast irradiation 

 Marqa, M.F. and Caudrelier, J.M. and Betrouni, N.   Brachytherapy. 2012; 11(5): 359368   [Pubmed]   5 
Physicsaspects of dose accuracy in high dose rate (HDR) brachytherapy: Source dosimetry, treatment planning, equipment performance and in vivo verification techniques 

 Palmer, A. and Bradley, D. and Nisbet, A.   Journal of Contemporary Brachytherapy. 2012; 4(2): 8191   [Pubmed]   6 
Primary gliosarcoma – clinical experience from a regional cancer centre in north India 

 Ahitagni Biswas, Narendra Kumar, Pankaj Kumar, Rakesh Kumar vasishta, Kirti Gupta, Suresh C. Sharma, Firuza Patel, Suresh Narain Mathuriya   British Journal of Neurosurgery. 2011; 25(6): 723   [VIEW]  [DOI]   7 
Sarcomatoid squamous cell carcinoma of uterine cervix 

 Young Park, Jee Hae Cho, Jong Sil Lee, Jeong Kyu Shin, Won Jun Choi, Jong Hak Lee, Won Young Paik   Korean Journal of Obstetrics and Gynecology. 2010; 53(10): 946   [VIEW]  [DOI]   8 
MUTomo: Independent dose validation software for helical tomotherapy 

 Papanikolaou, N., He, W., Luis, A.V.Q., Gutierrez, A., Stathakis, S., Alkhatib, H., Shi, C.   Journal of Cancer Science and Therapy. 2010; 2(5): 145152   [Pubmed]   9 
Sensitivity study to evaluate the dosimetric impact of offaxis ratio profiles misalignment on tomotherapy second dose validation 

 He, W., Luis, A., Vazquez, Q., Shi, C., Papanikolaou, N.   Technology in Cancer Research and Treatment. 2010; 9(5): 515522   [Pubmed]  



