469 research article Quality assurance procedures based on dosimetric, gamma analysis as a fast reliable tool for commissioning brachytherapy treatment planning systems Grzegorz Zwierzchowski12, Grzegorz Bieleda12, Janusz Skowronek13 1 Poznan University of Medical Sciences, Faculty of Health Sciences, Poznana, Poland 2 Greater Poland Cancer Centre, Medical Physics Department, Poznan, Poland 3 Greater Poland Cancer Centre, Brachytherapy Department, Poznan, Poland Radiol Oncol 2017; 51(4): 469-474. Received 25 November 2016 Accepted 10 October 2017 Correspondence to: Dr. Grzegorz Zwierzchowski, Greater Poland Cancer Centre, Medical Physics Department / Brachytherapy Department, Garbary 15 Str, PL 61-855 Poznan, Poland. E-mail: grzegorz.zwierzchowski@gmail.com Disclosure: No potential conflicts of interest were disclosed. Background. Fast and easily repeatable methods for commissioning procedures for brachytherapy (BT) treatment planning systems (TPS) are needed. Radiochromic film dosimetry with gamma analysis is widely used in external beam quality assurance (QA) procedures and planar film dosimetry is also increasingly used for verification of the dose distribution in BT applications. Using the gamma analysis method for comparing calculated and measured dose data could be used for commissioning procedures of the newly developed TG-186 and MBDCA calculation algorithms. The aim of this study was dosimetric verification of the calculation algorithm used in TPS when the CT/MRI ring applicator is used. Materials and methods. Ring applicators with 26 and 30 mm diameters and a 60 mm intra-uterine tube with 60° angle were used for verification. Gafchromic® EBT films were used as dosimetric media. Dose grids, corresponding to each plane (dosimetric film location), were exported from the TPS as a raw data. Gafchromic® films were digitized after irradiation. gamma analysis of the data were performed using the OMNI Pro I'mRT® system, as recommended by the AAPM TG-119 rapport criterion for gamma analysis of 3%, 3 mm and a level of 95%. Results. For the 26 mm and 30 mm rings, the average gamma ranged, respectively, from 0.1 to 0.44 and from 0.1 to 0.27. In both cases, 99% of the measured points corresponded with the calculated data. Conclusions. This analysis showed excellent agreement between the dose distribution calculated with the TPS and the doses measured by Gafchromic films. This finding confirms the viability of using film dosimetry in BT. Key words: brachytherapy; quality assurance; film dosimetry; ring applicator Introduction Radiochromic film dosimetry with gamma analysis is widely used in quality assurance procedures for external beam radiotherapy (EBRT). The main advantage of this method is the ability to collect dosimetric data with very high planar resolution in contrast to other methods that are based on point dose measurements. Introducing fast and easily repeatable methods for commissioning procedures for brachytherapy (BT) treatment planning systems (TPS) is a complex undertaking in most cases.12 The relatively low energy range of BT sources induce very high spatial dose gradients in irradiated volumes; as a result, it is difficult to use point dose measurement in high dose volumes because small errors in detector positioning can induce large uncertainties in the measured values. Planar film dosimetry is increasingly used to verify dose distributions in BT applications. When the stepping source is used, the overall dose distri- Radiol Oncol 2017; 51(4): 447-454. doi: 10.1515/raon-2017-0024 470 Zwierzchowski G et al. / Film dosimetry procedures for quality assurances in brachytherapy bution in the medium is a product of the contribution from each source position and the modulated step time, which are governed by optimization routines.3,4 An advantage of film dosimetry is that it offers the possibility of collecting the dose data in the planar area and using techniques of additive measurements.5,6 Self-developing dosimetry films have long been used to successfully verify dose distributions.7,8 Although this approach still needs to be customized — in terms of its technical possibilities and BT TPS realization regime in particular facilities. Most TPS used for treatment plan preparation are still based on the TG-43 recommendations. These algorithms and optimization routines have been verified by many authors in homogenous conditions.9,10 The Sivert integral and modular dose calculation models allow calculations of dose rates under the assumption that the elemental source position is surrounded by a homogenous water en-vironment.11,12 The TG-186 recommendations use the MBDCA (Model Based Dose Calculation Algorithm). Because these recommendations deliver accurate tissue segmentation and take into account the elemental composition of the structure, their use continues to make inroads because they provide better accuracy in BT dosimetry.13,14,15 Ring applicators are widely used in high dose rate (HDR) BT applications for patients with cervix cancer. This model is convenient to use due to its fixed geometry and availability in versions com- plane A plane B © X 1 T J ^ plan* FIGURE 1. Geometry of the measurements, tandem ring applicator and Gafchromic® film located in planes (A), (B) and (C), respectively. patible with computed tomography (CT)/magnetic resonance imaging (MRI) and many variations including additional interstitial needles. The most common observed problems with setup are associated with source positioning uncertainties during treatment compared to the dwell position placements used to calculate the dose distribution. This problem occurs when source path models from libraries are used and can also occur when the path is modeled manually.16 To our knowledge, no previous studies have attempted to perform dosimetric verification of the calculation algorithm used in TPS when the CT/ MRI ring applicator is used. For this reason, we conducted the present study, in which we also sought to develop treatment planning commissioning procedures for current and further use based on the planar film dosimetry. Materials and methods For the preparation of the treatment plans Oncentra Brachy® 4.3 were used as this is the main treatment planning software used in authors department. This version of the system was equipped with TG43 based calculation algorithms only. Two treatment plans were prepared for evaluation purposes. Standard CT/MRI ring applicator sets were used. The setup was based on 26 and 30 mm diameter rings and a 60 mm intra-uterine tube with 60° bent angle. Reconstruction of the application geometry was based on the applicator library module in Oncentra Brachy® 4.3. The reference dose of 3 Gy was prescribed to the standard Manchester A points using a HDR iridium source. Gafchromic® EBT films were used as the do-simetric media. Irradiation setup was based on a PMMA phantom. The phantom was build using large PMMA blocks with prepared cavity for intrautherine probe. Blocks of PMMA with different thickness were used to prepare repeatable setup and for assuring dosimetric media flatness and proper (parallel and perpendicular) positioning relative to the applicator. The dose to water in PMMA (calibration) and dose to water in water (measurements and TG43-based dose calculations) are not equivalent. However, D and D differ only by ' w,w w,PMMA j j 0.8% for distances of 25 mm from an 192Ir source. For this reason, relative character of the performed measurements and that PMMA is commonly available in radiotherapy facilities this material was chosen and used. Radiol Oncol 2017; 51(4): 469-474. 471 Zwierzchowski G et al. / Film dosimetry procedures for quality assurances in brachytherapy Dosimetric films (Gafchromic® EBT) were then placed in three planes. The first plane (plane A) was located at the surface of the ring part of the applicator with a round hole for the intra-uterine probe. Plane B was located on the surface of the probe in parallel orientation to the plane where Manchester A points were located. The third film (plane C) was placed on the surface of the probe perpendicular to the previous (plane B) localization. Geometry of the measurements is presented on the Figure 1. After preparation of the treatment plan dose grids corresponding to each plane (dosimetric film localization), these were exported from the TPS as a raw data. The smallest available dose grid resolution was 1 mm x 1 mm, although under this condition only a small part of the data collected from film analysis could be used for comparisons. In the plane A area that included the hole data were removed from the analysis to avoid obvious differences in the no-film area. Additionally, the dose limit in the export module of the TPS was set to 400% of the reference dose and any dose values exceeding this limit were exported as 400% value in the dose grid. After 72 hours, the irradiated Gafchromic® films were digitized with a flat table scanner (Epson® Perfection V750 Pro), all with the same orientation. Bitmap representation of the digitized films was then converted to the dose data and then local gamma analysis of the data were performed using the OMNI Pro I'mRT ® package with 3% and 3 mm criteria. Calibration data for the films were collected by separately irradiating 14 sheets (20 mm x 30 mm) of Gafchromic® EBT (Lot #: 47207-031, ISP) films with doses ranging from 0.25 Gy to 8.0 Gy using an HDR Ir-192 (192-Ir-mHDR-v2) source. To assure homogenous dose distribution, the films were placed between two blocks of 25 mm thick PMMA and two catheters were placed above and below the films at a distance of 25 mm. The doses were prescribed to the dose points in the centre of the film. After 72 hours, the films were digitized with a flat table scanner (Epson® Perfection V750 Pro) with light source on the one side and the detector on the other side of the film, all with the same orientation. Mean values from the most homogenous central part of the film (10 mm x 5 mm) were calculated using the VeriSoft® package. In region of interest of 10 mm x 5 mm the dose variation was estimated below 5%. Then the calibration curve was prepared and used recalculate the optical density (analog to digital conversion value; ADC) to the doses.17 FIGURE 2. Digitized dosimetric films as a visual input data for two applicators setup and for all three analyzed dose planes. (I) Applicator ring diameter 26 mm, Gafchromic® film located in planes (A), (B) and (C), respectively. (II) Applicator ring diameter 30 mm: Gafchromic® film located in planes (A), (B) and (C), respectively. Results Co-registration of data imported from TPS and dose distribution from scanned film were guided by centre of the intrauterine probe. Dose was normalized to 400% as maximum measured dose by OmniPro® software for each analyzed film. The films were scanned at 250 DPI (Dots Per Inch) resolution, always with the same orientation. Treatment planning system export files contains dose data with 1 mm resolution therefore the data from film have to be downsampled by analyzing software. The dose inside the hole in the films - prepared for insertion intrauterine probe was manually changed to 0 Gy, in both - film data and TPS exported ones. We assumed that dose distributions (planned and measured) were consistent if the count of pixel values from 0.00 to 1.00 (blue and green representation) in the gamma analysis was over 95% and the average pixel value was lower than 1.00. Pixel values higher than 1.00 (yellow and red representation) show the regions of inconsistencies (Figures 3-8). The gamma analysis showed that all measured dose distributions were consistent with the planned distributions. Table 1 presents average and maximum gammas and percentage of the analyzed dose points that met the 3% and 3 mm criteria. It's difficult to clearly state the uncertainties during realizing this study (based on gamma analysis), the main source of possible error is the mechani- Radiol Oncol 2017; 51(4): 469-474. 472 Zwierzchowski G et al. / Film dosimetry procedures for quality assurances in brachytherapy 100%=700,00 cGy II I1111!1111!'........... 100%-700.00 cGy I" 2 2i I li 1 li 2 ■ 3 H I I im 360Ï 3m sat 2m 20ÖV 169V I 120S I 3m- 2«* I ' 1ED5 I 1SR ;--200% [ II 1005È=7Q0.0a cGy III .100% ■ 150% ■ 200% ■ 100% 1507« ■200% IV FIGURE 3. Data analysis for the 26 mm ring in plane A. (I)- dose distribution obtained from film, (II) - dose distribution obtained from TPS, (III) - isodoses obtained from film (dots) and TPS (lines), (IV) - gamma factor map. FIGURE 4. Data analysis for the 26 mm ring in plane B. (I) - dose distribution obtained from film, (II) - dose distribution obtained from TPS, (III) - isodoses obtained from film (dots) and TPS (lines), (IV) - gamma factor map. FIGURE 5. Data analysis for the 26 mm ring in plane C. (I) - dose distribution obtained from film, (II) - dose distribution obtained from TPS, (III) - isodoses obtained from film (dots) and TPS (lines), (IV) - gamma factor map. FIGURE 6. Data analysis for the 30 mm ring in plane A. (I) - dose distribution obtained from film, (II) - dose distribution obtained from TPS, (III) - isodoses obtained from film (dots) and TPS (lines), (IV) - gamma factor map. FIGURE 7. Data analysis for the 30 mm ring in plane B. (I) - dose distribution obtained from film, (II) - dose distribution obtained from TPS, (III) - isodoses obtained from film (dots) and TPS (lines), (IV) - gamma factor map. FIGURE 8. Data analysis for the 30 mm ring in plane C. (I) - dose distribution obtained from film, (II) - dose distribution obtained from TPS, (III) - isodoses obtained from film (dots) and TPS (lines), (IV) - gamma factor map. cal positioning of the films in phantom. Another problem could be find in the spatial aligning of the two types of analyzed data. Authors decided to manually register two data series using intrauter-ine probe position assuming that the film was flat and properly placed. Authors estimated possible positioning errors at less than 1 mm level. During the calibration - homogenous irradiation of film detector with a point source in one position isn't achievable, therefore a special arrangement Radiol Oncol 2017; 51(4): 469-474. 473 Zwierzchowski G et al. / Film dosimetry procedures for quality assurances in brachytherapy of irradiation based on Khushdeep Singh's18 work with two catheters was used, which allowed to deposit not less than 95% with standard deviation of 0,82% of reference dose on a defined film area. Discussion In the study, we have attempted to use the typical BT setup with tandem ring applicator as a test platform to develop a commissioning procedure to verify the MBDCA algorithm recommended in TG-186. The main finding was that self-developing flat film dosimetry is a fast and reliable commissioning method in TG-43 conditions that could easily be adopted to almost any clinical setup in which point dose dosimetry is difficult to use and cannot provide valuable information. The data packages obtained from digitizing the Gafchromic® films were prepared for comparison with the planar dose distribution data exported from the TPS. The OMNI Pro I'mRT® software typically used for IMRT was used to perform gamma analysis of the data. In this approach, a commercially available method, typically used for external beam radiotherapy (EBRT), was used to compare the calculated and measured dose distribution. Overall dose distributions from EBDRT and BT is significantly different in terms of the dose gradients.19 For brachytherapy dose distributions we observe much larger dose ranges in small volumes. As a result, the use of point dose detectors to make measurements is highly difficult due to positioning errors. In contrast, film dosimetry seems to be offer a fast and reliable method for commissioning calculation algorithms, and planning and treatment delivery procedure.22021 The common rule for the gamma criterion when BT verification using gafrchromic films is performed has not been established. For purposes of this study, we adopted the AAPM TG-119 rapport criterion for gamma analysis: 3%, 3 mm and a level of 95%. The results of the analysis were acceptable for two applicator size on all three planes at the 95% level and above. This results confirms the correctness of this measuring method and is a positive reference for further analysis in more complex cases with more complicated density and geometry setup. Another important issue is the use of shields and commissioning of the calculation result for shielded applicators, where high density material is placed in the close proximity of the source.22 TABLE 1. Results of Gamma analysis for the 26 mm ring and 30 mm ring applicator Ring, 26 mm Gamma Gamma Gamma average max (0.00 - 1.00) plane A 0.22 1.42 99.04% plane B 0.10 1.29 99.31% plane C 0.44 1.75 98.88% Ring, 30 mm Gamma Gamma Gamma average max (0.00 - 1.00) plane A 0.25 1.96 98.11% plane B 0.27 2.00 97.94% plane C 0.10 1.22 99.54% In most TPS in current use, the elemental composition of the structures is not taken into account during the calculation (atomic number - Z) and dose distribution is based on geometric models of the shields. The TG-186 recommendations with MBDCA introduce the need to develop reliable verification methods that are more convenient and accurate, which can be performed in a more repeat-able manner than point dose methods.2324 Better accuracy in BT dosimetry appears to be a common need when the benefits from accurate tissue segmentation and the structure's elemental composition are considered as an important step up.13,14,25 It bears mentioning that for the analysis performed in the present study, we only used a small amount of the data from the digitized dosimetric films. The planar resolution of the films is very high and limited only by the chemical structure of the dosimetric media itself and physical resolution of the scanner. On the other hand, the maximum resolution of the exported dose grid was only 1 mm. But the most important factor for QA in BT is the possibility to export data exactly from the plane where the film was located during the measurements. This allows for the design of very effective phantoms with convenient and repeatable geom-etry.2,26 The film dosimetry used to verify the dose distribution and also for the direct reconstruction of the source path should be considered the method of choice for commissioning these newly-designed applicators. Precise reconstruction of the real source positions allows the dose distributions calculated by the planning system to be checked and leads to more conformal treatment planning.27 It should be noted that introducing effective, repeat-able and easy-to-use QA procedures for BT is es- Radiol Oncol 2017; 51(4): 469-474. 474 Zwierzchowski G et al. / Film dosimetry procedures for quality assurances in brachytherapy sential. Numerous factors, including relatively high doses, dose gradients, limited number of fractions, and advanced optimization algorithms with extensive dwell time modulation, make it nearly impossible to apply in-vivo dosimetry, thus leaving virtually no room to implement plan corrections during treatment. In these conditions, developing convenient "offline" methods for commissioning calculation routines, applicators and also for more frequently used QA procedures is very important and merits more study. Conclusions The analysis performed in this study showed excellent agreement between the dose distributions calculated using TPS and the doses measured by Gafchromic films. This confirms the viability of using film dosimetry in brachytherapy. The proposed commissioning procedure for further use with the MBDCA algorithm was established for use at the authors' facility, and it seems likely that the same procedure could be replicated at other centers. References 1. Malicki J. The importance of accurate treatment planning, delivery, and dose verification. Rep Pract Oncol Radiother 2012; 17: 63-5. doi: 10.1016/j. rpor.2012.02.001 2. Palmer AL, Lee C, Ratcliffe AJ, Bradley D, Nisbet A. Design and implementation of a film dosimetry audit tool for comparison of planned and delivered dose distributions in high dose rate (HDR) brachytherapy. Phys Med Biol 2013; 58: 6623-40. doi: 10.1088/0031-9155/58/19/6623 3. Adamczyk M, Zwierzchowski G, Malicki J, Skowronek J. Evaluation of clinical benefits achievable by using different optimization algorithms during real time prostate brachytherapy. Phys Medica 2013; 29: 111-7. doi: 10.1016/j. ejmp.2011.12.005 4. Rivard MJ, Beaulieu L, Mourtada F. Enhancements to commissioning techniques and quality assurance of brachytherapy treatment planning systems that use model-based dose calculation algorithms. Med Phys 2010; 37: 2645-58. doi: 10.1118/1.3429131 5. Granero D, Vijande J, Ballester F, Rivard MJ. Dosimetry revisited for the HDR 192Ir brachytherapy source model mHDR-v2. Med Phys 2011; 38: 487-94. doi: 10.1118/1.3531973 6. Perez-Calatayud J, Ballester F, Das RK, Dewerd LA, Ibbott GS, Meigooni AS, et al. Dose calculation for photon-emitting brachytherapy sources with average energy higher than 50 keV: Report of the AAPM and ESTRO. Med Phys 2012; 39: 2904-29. doi: 10.1118/1.3703892 7. Brown T, Hogstrom RK. Dose-response curve of EBT, EBT2 and EBT3 radi-ochromic films to synchrotron-produced monochromatic x-ray beams, Med Phys 2012; 39: 7412-7. doi: 10.1118/1.4767770 8. Chiu-Tsao ST, Ho Y, Shankar R, Wang L, Harrison LB. Energy dependence of response of new high sensitivity radiochromic films for megavolt-age and kilovoltage radiation energies, Med Phys 2005; 32: 3350-4. doi: 10.1118/1.20165467 9. DeWerd LA, Ibbott GS, Meigooni AS. A dosimetric uncertainty analysis for photon-emitting brachytherapy sources: report of AAPM Task Group No. 138 and GEC-ESTRO. Med Phys 2011; 38: 782-801. doi: 10.1118/1.3533720 10. Hyer DE, Sheybani A, Jacobson GM, Kim Y. The dosimetric impact of heterogeneity corrections in high-dose rate (192)Ir brachytherapy for cervical cancer: investigation of both conventional Point-A and volume-optimized plans. Brachytherapy 2012; 11: 515-20. doi: 10.1016/j.brachy.2012.01.011 11. de Oliveira CA, Lopes MD, Matos A. Commissioning of brachytherapy module of Oncentra MasterPlan treatment planning system. J Contemp Brachyther 2009; 1: 224-30. 12. Palmer A, Bradley D, Nisbet A. Physics-aspects of dose accuracy in high dose rate (HDR) brachytherapy: source dosimetry, treatment planning, equipment performance and in vivo verification techniques. J Contemp Brachytherapy 2012; 4: 81-91. doi: 10.5114/jcb.2012.32558 13. Beaulieu L, Carlsson Tedgren A, Carrier JF, Davis SD, Mourtada F, Rivard MJ, et al. Report of the Task Group 186 on model-based dose calculation methods in brachytherapy beyond the TG-43 formalism: current status and recommendations for clinical implementation. Med Phys 2012; 39: 620836. doi: 10.1118/1.4747264 14. Mikell JK, Klopp AH, Gonzalez GM, Kisling KD, Price MJ, Berner PA, et al. Impact of heterogeneity-based dose calculation using a deterministic grid-based Boltzmann equation solver for intracavitary brachytherapy. Int J Radiat Oncol Biol Phys 2012; 83: 417-22. doi: 10.1016/j.ijrobp.2011.12.074 15. Oliveira S, Teixeira N, Fernandes L, Teles P, Vaz P. Dosimetric effect of tissue heterogeneity for 125I prostate implants. Rep Pract Oncol Radiother 2014; 19: 392-8. 16. Awunor OA, Dixon B, Walker C. Direct reconstruction and associated uncertainties of 192Ir source dwell positions in ring applicators using gafchromic film in the treatment planning of HDR brachytherapy cervix patients. Phys Med Biol 2013; 58: 3207-25. doi: 10.1088/0031-9155/58/10/3207 17. Zwierzchowski G, Bielçda G, Skowronek J, Mazur M. Film based verification of calculation algorithms used for brachytherapy planning-getting ready for upcoming challenges of MBDCA. J Contemp Brachytherapy 2016; 8: 326-35. doi: 10.5114/jcb.2016.61828 18. Singh K. Investigation of the energy response of EBT-2 gafchromic film model. Montreal: Department of Medical Physics, McGill University; 2011. 19. Manikandan A, Biplab S, David PA, Holla R, Vivek TR, Sujatha N. Relative dosimetrical verification in high dose rate brachytherapy using two-dimensional detector array IMatriXX. J Med Phys 2011; 36: 171-5. doi: 10.4103/0971-6203.83491 20. Vijande J, Ballester F, Ouhib Z, Granero D, Pujades-Claumarchirant MC, Perez-Calatayud J. Dosimetry comparison between TG-43 and Monte Carlo calculations using the Freiburg flap for skin high-dose-rate brachytherapy. Brachytherapy 2012; 11: 528-35. doi: 10.1016/j.brachy.2011.11.005 21. Rivard MJ, Venselaar JL, Beaulieu L. The evolution of brachytherapy treatment planning. Med Phys 2009; 36: 2136-53. doi: 10.1118/1.3125136 22. Hira M, Podgorsak MB, Jaggernauth W, Malhotra HK. Measurement of dose perturbation around shielded ovoids in high-dose-rate brachytherapy. Brachytherapy 2011; 10: 232-41. doi: 10.1016/j.brachy.2010.08.008 23. Uniyal SC, Sharma SD, Naithani UC. A dosimetry method in the transverse plane of HDR Ir-192 brachytherapy source using gafchromic EBT2 film. Phys Med 2012; 28: 129-33. doi: 10.1016/j.ejmp.2011.03.005 24. Landry G, Reniers B, Murrer L, Lutgens L, Gurp EB, Pignol JP, et al. Sensitivity of low energy brachytherapy Monte Carlo dose calculations to uncertainties in human tissue composition. Med Phys 2010; 37: 5188-98. doi: 10.1118/1.3477161 25. Kirisits C, Rivard MJ, Baltas D, Ballester F, De Brabandere M, van der Laarse R, et al. Review of clinical brachytherapy uncertainties: analysis guidelines of GEC-ESTRO and the AAPM. Radiother Oncol 2014; 110: 199-212. doi: 10.1016/j.radonc.2013.11.002 26. Niatsetski Y, Fekkes S, Vreeken H. Source path measurements for ring applicators. [Abstract]. Radiother Oncol 2011; 99: 278. Poster No. 695 27. Berger D, Dimopoulos J, Georg P, Georg D, Pötter R, Kirisits C. Uncertainties in assessment of the vaginal dose for intracavitary brachytherapy of cervical cancer using a tandem-ring applicator. Int J Radiat Oncol Biol Phys 2007; 67: 1451-9. doi: 10.1016/j.ijrobp.2006.11.021 Radiol Oncol 2017; 51(4): 469-474.