|
Journal of Cancer Research and Therapeutics, Vol. 4, No. 4, October-December, 2008, pp. 173-177 Original Article A dose verification method for high-dose-rate brachytherapy treatment plans Kumar Rajesh, Sharma SD, Vijaykumar C, Deshpande Sudesh, Sharma PK, Vandana S, Philomena A, Chilkulwar RaviH Radiological Physics and Advisory Division, Bhabha Atomic Research Centre, CT and CRS, Anushaktinagar, Mumbai-400 094 Code Number: cr08051 Abstract Aim: To evolve a fast dose verification method for high-dose-rate (HDR) brachytherapy treatment plans and to demonstrate its applicability in different clinical cases.
Keywords: Brachytherapy, dose verification, dosimetry, TG-43 formalism, treatment planning Introduction High-dose-rate (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 TG-43 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 TG-43 formalism [2] was used for dose calculation in the verification code:
where Sk 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 k(r) = k1rk2 + k3r + k4; 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 3a-c] 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 patient-specific quality assurance. A simple, fast, and accurate method of dose verification is required to fulfil this requirement. The AAPM TG-43 formalism-based 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 TG-43 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
Copyright 2008 - Journal of Cancer Research and Therapeutics The following images related to this document are available:Photo images[cr08051f2.jpg] [cr08051t3.jpg] [cr08051t2.jpg] [cr08051f4.jpg] [cr08051f1.jpg] [cr08051f3.jpg] [cr08051t1.jpg] |
|