- Research
- Open access
- Published:
Dose calibration of Gafchromic EBT3 film for Ir-192 brachytherapy source using 3D-printed PLA and ABS plastics
3D Printing in Medicine volume 5, Article number: 3 (2019)
Abstract
3D printing technology has allowed the creation of custom applicators for high dose rate (HDR) brachytherapy, especially for complex anatomy. With conformal therapy comes the need for advanced dosimetric verification. It is important to demonstrate how dose to 3D printed materials can be related to dose to water. This study aimed to determine dose differences and uncertainties using 3D printed PLA and ABS plastics for Radiochromic film calibration in HDR brachytherapy.
Gafchromic EBT3 film pieces were irradiated in water with an Ir-192 source at calculated dose levels ranging from 0 to 800 cGy, to create the control calibration curve. Similarly, film was placed below 3D printed PLA and ABS blocks and irradiated at the same dose levels calculated for water, ranging from 0 to 800 cGy. After a 72-h development time, film pieces were scanned on a flatbed scanner and the median pixel value was recorded in the region of highest dose. This value was converted to net optical density (NOD). A rational function was used to fit a calibration curve in water that relates NOD to dose for red, green, and blue color channels. Based on this fitted curve, ABS and PLA NOD values were used to estimate dose in 3D printed plastics.
From the fitted calibration curve, mean residual error between measured and planned dose to water was less than 1% for each color channel at high dose levels. At high dose levels, ABS and PLA mean residual errors were about 6.9 and 7.8% in the red channel, while 5.2 and 5.7% in the green channel. Combined uncertainties measured to be about 6.9% at high dose levels. This study demonstrated dose differences and uncertainties using 3D printed applicators for HDR Ir-192 brachytherapy.
Introduction
Custom applicators are useful for treating superficial tumors with high dose rate (HDR) brachytherapy. They allow highly conformal dose delivery, and are especially beneficial for oblique surfaces, such as the face [1, 2]. The Freiburg Flap is a commonly used customizable applicator for treating cutaneous lesions at depths less than 5 mm3. While this applicator is useful for treating lesions located on relatively uniform surfaces, it has difficulty conforming to irregularly-shaped structures such as the nose or ear. Recently, the proliferation of 3D printing technology has enabled the creation of customized applicators for the delivery of highly conformal radiotherapy treatment, even in the setting of complex anatomic geometries [1].
With advancement of applicators and conformal therapy comes the need for accurate dosimetry methods. Film dosimetry has proved to be a high-resolution tool for radiotherapy treatment verification. Radiochromic film is useful for HDR brachytherapy quality assurance and is more practical than use of radiographic film or ion chambers [4, 5].
Radiochromic film darkens with radiation exposure and dose, measured by changes in optical density without the need for chemical processing. Radiochromic film dosimetry allows the conversion of net optical density (NOD) to dose based on a source-specific calibration curve in water, or solid water [4]. To create a calibration curve, film is exposed in water at increasing prescribed dose levels. Film response (darkening) is measured in a standard flatbed scanner as NOD. The NOD can be related to dose by a rational function. This fitted function, known as a calibration curve, can then be used as a dosimetry tool for future dose measurements made with Radiochromic film.
Common 3D printed plastics such as acrylic butadiene styrene (ABS) and polymethylmethacrylate (PLA) have been shown to be near-water equivalent [6]. Currently, film calibration for the HDR source is done in water. The process of film calibration in water can be laborious and time consuming. In a busy clinical setting, the use of 3D printed PLA and ABS plastics could offer a faster alternative for obtaining film calibration curves. The goal of this work was to obtain the relationship of doses and uncertainties in 3D printed PLA and ABS plastics when used in lieu of water for film calibration for HDR Iridium-192 brachytherapy sources. Furthermore, as part of a large-scale study to design custom 3D-printed superficial HDR applicators, it is important to demonstrate how dose in 3D printed materials can be related to dose to water.
Materials and methods
Source
An HDR remote afterloader (Nucletron, Elekta, Stockholm, Sweden) was used with a MICROSELECTRON V2 Iridium-192 source (Alpha-Omega Services, Inc., Edgerly, LA).
Film
Gafchromic EBT3 film was used in this study and is the most current model of Radiochromic film available. EBT3 film is near tissue equivalent with an effective atomic number of 7.26 [7]. It has two polyester substrate layers (125 μm) surrounding an active layer (28 μm), creating a symmetric structure that eliminates side orientation dependence. EBT3 film has been shown to be energy independent at doses above 100 cGy, and at energies greater than 100 keV photons [8], making it useful for Ir-192 sources, which have an average photon energy of 380 keV. In addition, EBT3 characteristics such as uniformity, orientation, and energy dependence have been widely studied [9, 10].
All film was obtained from the same lot number, 06131702. Sheets were cut into fourths, creating 10.2 × 12.7 cm2 pieces and marked to maintain side and orientation dependence. Film was handled and stored in accordance with American Association of Physicists in Medicine (AAPM) Task Group 55 recommendations and guidelines [11].
Dose to water measurements
A custom 3D printed holder (made of ABS) was used to take film measurements in water (see Fig. 1). Two opposed catheters were placed 5 cm from the film's center to deliver a total dose of 50, 100, 200, 400, 600, and 800 cGy. Treatment planning software (Oncentra, Elekta, Stockholm, Sweden) calculated doses based on AAPM Task Group 43 (which does not correct for inhomogeneities) [12, 13]. A non-irradiated film sample was used as a 0 cGy background measurement.
3D-printed material measurements
A 4x4x7 cm3 ABS block was 3D printed using a commercially available 3D printer (Taz 6, Aleph Objects Inc., Loveland, CO) with 100% infill and a 0.2 mm layer height. The center of the catheter channel was offset by 1 cm from midline with the film placed along the surface of the block 3 cm from the catheter (see Figs. 2, 3 and 4). The treatment planning software (Oncentra, Elekta, Stockholm, Sweden) was again used to measure dose to water ranging 50–800 cGy to the film. This procedure was repeated for a similar block made of PLA. ABS and PLA were assumed to be near water equivalent in the treatment planning system (TPS) [6].
Scanning the film
A flatbed scanner (Epson Expression 11000XL, Seiko Epson Corp., Tokyo, Japan) and associated EPSON SCAN, were used to scan the film approximately 72 h post-irradiation. Previous findings have shown minimal variation in NOD beyond a 24 h development time [10]. A foam board positioning template was used to achieve reproducible and uniform position in the center of the scanner. Films were scanned three times each to determine scanning consistency. RBG-positive images were acquired with a spatial resolution of 72 dpi, and a depth of 48 bits (16 bits per color channel). Images were saved in tiff format. The ImageJ software platform (National Institutes of Health, Bethseda, MD) was used to analyze the film with a 9 × 9 pixel region of interest (ROI) selected at the highest point of exposure. The median pixel value (related to image’s intensity) from the ROI was then recorded for each color channel (red/blue/green). Median pixel value can then be converted to NOD by the following formula:
Creating and using a calibration curve
The relationship between film response and dose was determined with a calibration curve based on water measurements. A rational function was used to fit the NOD data with the expected dose from the TPS [14]:
Using coefficients, a, b, and c to create calibration curve, dose to the can be measured based on changing film response.
Uncertainty analysis
In this study, uncertainty was determined based on precision of measurements and calculations made. Based on AAPM Task Group 43 uncertainties, both random (type A), and systematic (type B), are measured [12, 13]. Uncertainties were divided into two parts; determining dose from NOD, and film exposure. An estimate of combined uncertainty was calculated using a square root of the sum of squared individual uncertainty components, as recommended by AAPM [13].
Results
Calibration curve fit
The parameters shown in Table 1 were fit for red, green, and blue channels, based on the relationship between planned dose and NOD. The dosimetric error was measured as a percent difference between the measured dose with film and the expected dose from the TPS. An absolute mean error was determined to quantify the sensitivity of each channel. Since EBT3 film is less accurate at low doses, the absolute mean error was divided into low dose (≤100 cGy) and high dose (> 100 cGy).
For water measurements all channels showed high dose mean error less than 1%. At low doses, there were inaccuracies in the calibration curve fit, as shown in Table 2. The blue channel showed the highest error, especially at the 50 cGy dose level. Since the red and green channels showed little dosimetric error at both high (≤0.1%) and low (< 5.0%) dose levels, they will be used for the remainder of the analysis. Previous studies emphasize that red and green channels were most sensitive for high dose levels [4, 10].
Dose measurements in PLA and ABS
The red channel calibration curve from water (in Table 1) was used to fit NOD values for ABS and PLA film of increasing dose levels. Table 3 shows the measured dose and the percent error from the expected dose based on the red channel fit. Again, mean error was divided into low and high dose to account for inaccurate film response for doses at and below 100 cGy. A mean error of 6.9% and 7.8%, for ABS and PLA respectively for high doses using the red channel as seen in Table 3. The green channel measured a high dose mean error of 5.2% and 5.7% for ABS and PLA respectively. Table 4 displays green channel dose measurements. Figures 5 and 6 demonstrate how error changes with increasing dose for water, ABS, and PLA.
Discussion
Uncertainty analysis
Uncertainties reported in AAPM Task Group 43 were both statistical, random (type A) and systematic, nonrandom (type B) [13]. Also suggested by AAPM Task Group 43, the combined uncertainty was estimated using a simple root sum of squares of individual components. Our uncertainty analysis was broken into two parts; obtaining dose from NOD, and film exposure. Timer error and calibration curve fit are dose dependent errors, therefore both low and high dose uncertainties were estimated. Table 5 summarizes sources of uncertainty and combined uncertainty for high and lose dose measurements in film.
Determining dose from NOD
Scanning consistency added minimal statistical uncertainty, less than 0.1%. Scanning consistency was determined from the difference between six separate scans of the same film piece. Film uniformity was determined by measuring difference between four quadrants of a background film sample. An additional source to consider was the ROI size. Decreasing the ROI size would reduce the amount of data sampled and may overestimate the optical density. This error was again more prominent at lower dose levels, with a low signal to noise ratio. Previous literature has found little uncertainty for a ROI 4 × 4 mm2 or larger [15, 16]. Our 9 × 9 pixel ROI was equivalent to a 4.1 × 4.1 mm2 ROI. Compared to a 2.3 × 2.3 mm2 ROI, net optical density for our larger ROI differed by about 1.5% for red channel water measurements.
With the obtained fitting parameters, the water calibration curves gave residual error for dose, especially at low dose levels as outlined in Table 2. For this reason, a systematic uncertainty error was added to the analysis for low and high dose levels. Combined uncertainty from Dose to OD procedure was estimated to be 5.1 and 2.0% for low and high dose, respectively.
Film exposure
Positional accuracy between the film and source was a large cause of dosimetric error. In water measurements, a slight bend in the film may cause the source to be incorrectly placed up to 1 mm. Although a double catheter setup was used to minimize film position error, it was still present. Source film position error was determined by relating exposure at 5 cm to exposure at 5.1 cm away by the inverse square relationship. For example, a 1 mm uncertainty would lead to a 4% error, for water measurements at 5 cm.
Scatter equilibrium is an important consideration in superficial brachytherapy. At the skin surface there are less photons scattered from surrounding tissue, causing a smaller dose to the target. For ABS and PLA measurements, there was only 10 cm of solid water surrounding the 3D printed block/catheter (Fig. 4). Previous literature has suggested that full scatter equilibrium is achieved with at least 40 cm of water surrounding an Ir-192 source, otherwise the radial dose function can differ by 5–10%, thus under-dosing the target [17]. With the knowledge that our ABS and PLA film exposures did not achieve full scatter equilibrium, 5% systematic uncertainty is added. While this adds large uncertainty, it is a clinically relevant consideration. In practice full scatter equilibrium is not achievable for superficial HDR brachytherapy treatments within 5 mm of the skin surface [18]. Future studies should be completed to address loss of scatter conditions in superficial brachytherapy dosimetry.
Timer error added minimal uncertainty, but still was considered in our analysis. Low doses required smaller dwell times; thus, a higher proportion of timer error. For example, a 50 cGy dose to water required a dwell time of 50 s. A 0.4 s timer uncertainty would lead to a 0.8% timer error, while an 800 cGy dose and 800 s dwell time would produce an error of only 0.05%. Total uncertainty for low and high dose measurement in film was 8.4 and 6.9%, respectively.
Clinical relevance of results
In HDR brachytherapy, most prescribed fractionation schemes require doses of 500 cGy per fraction, or more [2, 3]. Our results showed that at doses beyond 200 cGy, dosimetric error was near the expected level of uncertainty. The green channel proved to have the most sensitivity by measuring dose within 6% of the expected dose at high dose levels. This finding was in alignment with previous literature [4, 10]. The red channel showed less sensitivity but still within 8% at high dose levels for both PLA and ABS.
Conclusions
A calibration curve was accurately obtained for water, and doses to water were compared to doses to 3D printed PLA and ABS plastics. For the calibration curve fit in water at low dose levels (100 cGy and below), there was an error up to 4.8% in the red channel. This error in the calibration curve fit is minimized for doses beyond 200 cGy. At high dose levels (200 cGy and higher), the calibration curve fit in water presented a mean error of 0.8, 0.3, and 0.7% in respective red, green, and blue channels. An accurate setup, with high positional accuracy is necessary to ensure that a quality calibration curve is acquired.
The results estimated the uncertainties and errors in measured doses when using 3D printed PLA and ABS plastics as summarized in Table 6. Distance and lack of scatter equilibrium were the largest source of uncertainty in our measurements. Measured doses in ABS and PLA were within the expected range of uncertainty.
These results have shown that 3D printed PLA and ABS plastics can be used in lieu of water for film calibration for HDR Iridium-192 brachytherapy sources, however systematic uncertainties related to source position, and scatter equilibrium conditions due to the phantom size should be carefully accounted for.
Abbreviations
- AAPM:
-
American association of physicists in medicine
- ABS:
-
Acrylic butadiene styrene
- ADCL:
-
Accredited dosimetric calibration laboratory
- HDR:
-
High dose rate
- NOD:
-
Net optical density
- PLA:
-
Polymethylmethacrylate
- ROI:
-
Region of interest
- TPS:
-
Treatment planning system
References
Jones EL, Baldion AT, Thomas C, Burrows T, Byrne N, Newton V, Aldridge S. Introduction of novel 3D-printed superficial applicators for high-dose-rate skin brachytherapy. Brachy. 2017;16(2):409–14.
Ouhib Z, Kasper M, Perez-Calatayud J, Rodriguez S, Bhatnagar A, Pai S, Strasswimmer J. Aspects of dosimetry and clinical practice of skin brachytherapy: the American brachytherapy society working group report. Brachy. 2015;14(6):840–58.
Guinot JL, Rembielak A, Perez-Calatayud J, Kovacs G. GEC-ESTRO ACROP recommendations in skin brachytherapy. Rad Onc. 2018;126:377–85.
Ayoobian N, Asl AS, Poorbaygi H, Javanshir MR. Gafchromic film dosimetry of a new HDR Ir-192 brachytherapy source. J Appl Clin Med Phys. 2015;17:194–205.
Mathot M, Sobczak S, Hoornaert MT. Gafchromic film dosimetry: four years experience using FilmQA pro software and Epson flatbed scanners. Physca Medica. 2014;30:871–7.
Craft DF, Kry SF, Balter P, Salehpour M, Woodward W, Howell RM. Material matters: analysis of density uncertainty in 3D printing and its consequences for radiation oncology. Med Phys. 2018. https://doi.org/10.1002/mp.12839.
León-Marroquín EY, Herrera-González JA, Camacho-López MA, Villarreal-Barajas JE, García-Garduño OA. Evaluation of the uncertainty in a EBT3 film dosimetry system utilizing net optical density. J Appl Clin Med Physv17. 2016:466–81.
Massillon JL, Chiu-Tsao ST, Domingo-Munoz I, Chan MF. Energy dependence of new Gafchromic EBT3 film: dose response curves for 50 kV, 6 and 15 MV X-ray beams. IJMPCERO. 2012;1:60–5.
Reinhardt S, Hillbrand M, Wilkens JJ, Assmann W. Comparison of Gafchromic EBT2 and EBT3 films for clinical photon and proton beams. Med Phys. 2012;39(8):5257–62.
Casanova Borca V, Pasquino M, Russo G, Grosso, Cante D, Sciacero P, Girelli G, La Porta MR, Tofani S (2013) GAFCHROMIC EBT3 film for IMRT dose verification. J Appl Clin Med Phys.
Niroomand-Rad A, Blackwell CR, Coursey BM, Gall KP, Galvin JM, McLaughlin WL, Meigooni AS, Nath R, Rodgers JE, Soares CG. Radiochromic film dosimetry: recommendations of AAPM radiation therapy committee task group 55. Med Phys. 1998;25:2093–115.
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:209–34.
Rivard MJ, Coursey BM, DeWerd LA, Hanson WF, Huq MS, Ibbott GS, Mitch MG, Nath R, Williamson JF. Update of AAPM task group no. 43 report: a revised AAPM protocol for brachytherapy dose calculations. Med Phys. 2004;31:633–74.
Micke A, Lewis DF, Yu X. Multichannel film dosimetry with nonuniformity correction. Med Phys. 2011;38(5):2523–34.
Bouchard H, Lacroix F, Beaudoin G, Kawrakow I. On the characterization and uncertainty analysis of Radiochromic film dosimetry. Med Phys. 2009;36:1931–46.
Sendani NG, Karimian A, Ferreira C, Alaei P. Technical note: impact of region of interest size and location in Gafchromic film dosimetry. Med Phys. 2018;45(5):2329–36.
Perez- Calatayud J, Granero D, Ballester F. Phantom size in brachytherapy source dosimetric studies. Med Phys. 2004;31(7).
Melhus CS, Rivard MJ. Approaches to calculating AAPM TG-43 brachytherapy dosimetry parameters for 137Cs, 125I, 192Ir, 103Pd, and 169Yb sources. Med Phys. 2006;33(6):1729–37.
Acknowledgements
The article processing charge for this publication has been funded generously by an unrestricted educational grant from Stratasys.
Ethics and consent to participate
Not applicable.
Funding
This research was partly funded by the 2017 Radiation Oncology Innovation Initiative Award, at the University of Minnesota.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Author information
Authors and Affiliations
Contributions
CO contributed to HDR setup, delivery, and film measurements, as well the analysis. CF also contributed to the HDR setup, delivery, film analysis and is the one who suggested the need for this research topic. CW and EE 3D designed and printed plastics. In addition, EE assisted with the linear model calibration and dosimetric analysis. DS aided HDR planning and delivery, as well as error analysis. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Consent for publication
Not applicable.
Competing interests
None of the authors has any relationship with Stratasys. The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Oare, C., Wilke, C., Ehler, E. et al. Dose calibration of Gafchromic EBT3 film for Ir-192 brachytherapy source using 3D-printed PLA and ABS plastics. 3D Print Med 5, 3 (2019). https://doi.org/10.1186/s41205-019-0040-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s41205-019-0040-4