|Year : 2019 | Volume
| Issue : 3 | Page : 207-212
Calibration of gafchromic EBT film using the Microtek ScanMaker 9800XL plus flatbed scanner with a modified one red-channel after three-channel method
Liyun Chang1, Sheng-Yow Ho2, Hueisch-Jy Ding1, Shyh-An Yeh1, Pang-Yu Chen3
1 Department of Medical Imaging and Radiological Sciences, I-Shou University, Kaohsiung, Taiwan
2 Department of Nursing, Chang Jung Christian University; Department of Radiation Oncology, Chi Mei Medical Center, Tainan, Taiwan
3 Department of Radiation Oncology, Sinlau Christian Hospital, Tainan, Taiwan
|Date of Submission||21-May-2019|
|Date of Decision||17-Aug-2019|
|Date of Acceptance||17-Aug-2019|
|Date of Web Publication||13-Sep-2019|
Dr. Liyun Chang
Department of Medical Imaging and Radiological Sciences, I-Shou University, Kaohsiung 82445
Dr. Pang-Yu Chen
Department of Radiation Oncology, Sinlau Christian Hospital, Tainan 70142
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Purpose: Using the Microtek ScanMaker 9800XL Plus (9800XL+) flatbed scanner, a method is presented to accurately calibrate EBT film, which cannot be calibrated simply using a general three-channel method because of the nonhomogeneous scanning. Materials and Methods: Through the percentage-depth-dose method, 6-MV photon beams with two different monitor units were delivered to eight EBT2 films, each of which was tightly sandwiched in a 30-cm cubic polystyrene phantom and positioned parallel to the central axis of the beam. Before and after irradiation, all films were scanned using the Microtek 9800XL+ scanner and the pixel values (PVs) were measured along the central axis of the beam on the film and fitted to the corresponding depth doses. Before calibration, the irradiated film image was first modified using a template matrix, which was generated using the prescanned background images. Then, a modified one red-channel after three-channel method was used to calibrate the film. Results: Without a template matrix, the three-channel method cannot be used because the PVs do not correspond to a rational fitting form. Using the proposed method, the difference between the fitted dose and the delivered dose is <2%. The green channel, and not the red, is found to have the largest dynamic range. Conclusion: The proposed technique allows the use of the three-channel method to calibrate film using a Microtek 9800XL+ scanner.
Keywords: EBT film, Microtek scanner, template matrix, three-channel calibration
|How to cite this article:|
Chang L, Ho SY, Ding HJ, Yeh SA, Chen PY. Calibration of gafchromic EBT film using the Microtek ScanMaker 9800XL plus flatbed scanner with a modified one red-channel after three-channel method. J Med Phys 2019;44:207-12
|How to cite this URL:|
Chang L, Ho SY, Ding HJ, Yeh SA, Chen PY. Calibration of gafchromic EBT film using the Microtek ScanMaker 9800XL plus flatbed scanner with a modified one red-channel after three-channel method. J Med Phys [serial online] 2019 [cited 2019 Oct 18];44:207-12. Available from: http://www.jmp.org.in/text.asp?2019/44/3/207/266848
| Introduction|| |
Gafchromic EBT film is widely used to determine the differences of dose distribution in the treatment volume for a patient between planning and delivery, largely because it is self-developing, has a near dose-to-water equivalence, gives high spatial resolution, is re-readable, features a relatively uniform dose-response across a wide range of photon energies ,, and is an inexpensive technique that uses commercially available flatbed document scanners.
Several generations of Gafchromic film have been developed, but only EBT2 and EBT3 film are recommended by Ashland company for the verification of beam-modulated techniques  because spatial nonhomogeneity is corrected by the yellow marker dye,,,,,,, it has a lower sensitivity to the visible spectrum and it can be used for repeated scans., Using a matte polyester substrate to avoid the formation of Newton's rings,,,, EBT3 film has a similar composition for the active layer and dosimetric properties as EBT2 and either side of the film can be used. The Ashland report states that the respective effective atomic numbers for EBT2 and EBT3 films are around 6.8 and 7.3, respectively, which is approximately equivalent to that for water. Therefore, they are eminently suited to patient dosimetry.,
The red, green, and blue color levels of EBT films can be analyzed for calibration. The red level calibration is sensitive to the dose of daily treatments within a range of 0–5 Gy ,, and the green and blue level calibrations are better for larger doses.,, Multichannel techniques generally give a better pass rate than a single-red-channel calibration for the verification of treatment planning using the local gamma index. Using the three-channel calibration method, Micke et al. proposed parameters that are dependent and independent of dose in the optical density (OD) domain to model and reduce deviations in dose-independent (DI) nonuniformities. van Hoof et al. noted  that nonuniformity corrections must be examined for larger fields. To further reduce nonuniformities, Chang et al. proposed a new parameter: channel-dependent (CD) scanner-nonuniformity in the OD domain. To reduce the calibration uncertainties, a power-function fitting process using only the red channel OD was used for the three-channel technique.
All of these methods mentioned use Epson scanners. An alternative cheaper choice is the Microtek 9800XL Plus (9800XL +) scanner, which is an A3 size scanner with a dynamic range of 48-bit. However, if the film dose is calibrated using a Microtek 9800XL +, the multi-channel method cannot be used, except the presented template matrix is used first.
| Materials and Methods|| |
Gafchromic™ EBT2 films (Ashland Inc., USA) of 8“×10” and a Microtek 9800XL + (MICROTEK™ company, Taiwan) scanner were used to generate 127 dpi tiff images before and after irradiation, for the purpose of calibration. The films were scanned using Scan Wizard Pro software (MICROTEK™ edition V7.26, Taiwan) and the functional mode, “averaging multiple-sampling of 2 lines,” was used.
To ensure that the film was oriented identically on the scanning bed for each scan, a homemade 43 cm × 36 cm × 0.2 cm acrylic transparent frame with a film-fitted hole was placed on the scanning platform in the portrait orientation. The scanner was warmed up by scanning a waste film eight times. The percentage-depth-dose (PDD) method  was used to calibrate the film using the 6 MV photon beam from an Elekta Synergy ® accelerator, which was previously calibrated at a depth of 5 cm according to the AAPM TG reports.,,,
Before the dose was delivered, the film was sandwiched in a 30 cm × 30 cm × 30 cm PTW™ (PTW-NEW YORK Corporation) RW3 polystyrene phantom and the midline that longitudinally separates the film into two equal parts was oriented to be coincident with the central beam axis. The film was also oriented parallel to the central beam axis and the upper edge of the film was parallel to the gantry rotation axis, which is conventionally the Y-axis in the literature.,,,, To measure the reference dose, a 0.6 cc Farmer chamber was located at a depth of 31 cm, 10 cm thick backup plates were placed under the entire 30 cm × 30 cm × 30 cm phantom, and the source-surface-distance was set to 100 cm. The reference dose measurement using the chamber at depth 31 cm may not be necessary if the monthly dosimetry calibration was performed just before the film-dose calibration. With field size 20 cm × 20 cm, four films were irradiated using a 290 monitor unit (MU) with delivered doses (Dd) around 320, 304, 280, 225, 177, and 108 cGy at respective depths, dmax, 3, 5, 10, 15, and 25 cm, and another four using an 88 MU with Dd around 97, 92, 85, 68, 54, and 33 cGy at respective depths, dmax, 3, 5, 10, 15, and 25 cm. The absorbed doses then ranged from around 30 cGy to 300 cGy on the midlines of the films. The doses are calculated using the verified PDDs and the measurements from the Farmer chamber.
After around 48 h, each film was rescanned with the same 127 dpi and “2-lines-average scanning (each color 16 bit).” All of the tiff format images were analyzed using MathWorks ® Matlab software (version R2017b). The inverse transmittance (IT) for the U channel (one of the R, G, and B channels), IU, for each image pixel is expressed as:
IU= (216 − 1)/PU, (1)
where PU is the pixel value (PV) for the U channel. However, equation (1) must be modified with a template matrix, TU, because the IT values do not correspond to a rational fitting form (equation 4) and the calibration cannot be correctly performed, as shown in [Figure 1]. Therefore, it is not possible to use the three-channel method  if TU is not used. Multiplying IU by the template matrix, the modified IT, NU, is written as:
|Figure 1: The inverse transmittance versus the delivered dose: the film was scanned using Microtek 9800+ without the template matrix TU|
Click here to view
NU= IU× TU (2)
The template matrix of channel U is given by:
TU= SU/MU (3)
where MU is a 2-dimensional PV matrix of channel U, using the average for the eight prescanned background films, and SU is a single value that is the average of all values for MU.
The delivered reference dose Dd as a function of NU along the central beam axis is fitted using the rational fitting form:
DU= (AU− BU× NU)/(NU− CU) (4)
where AU, BU, and CU are the fitting parameters for the calibration of the U channel and DU is the fitted dose for channel U. The inverse function is written as:
NU= (AU+ CU× DU)/(BU+ DU) (5)
Using the “one red-channel after three-channel (R-3C) method,” the OD of the U channel, OU, is the product of Odd, U, Odi and OL,U, which is the OD for CD and dose-dependent (DD), the OD for channel-independent and DI  and the OD for the CD scanner nonuniformity effect (the lateral effect), respectively. Therefore, the DD component of NU, which is Ndd,U, is expressed as:
where Ndd,U is equal to 10°dd, u. DR should be equal to DG and DB, but this is not the case because the Odi exists, but not the OL, U. This is because the depth doses are extracted from the middle line of the film, where no correction for lateral effect is required. Therefore, Odi is determined by minimizing K:
K = (DR– DG)2+ (DG– DB)2 + (DB–DR)2, (7)
by using equation (4) and (6) with the command “fminbnd” in MATLAB software. When Odi is eliminated, AU, BU, and CU are calculated again using equation (4) and are, respectively, renamed AnU, BnU, and CnU.
The profile doses for the 290-MU films with depths of 5, 10, and 20 cm (the request depths for inputting profiles into the planning system), which are related to the doses, D of 280, 225, and 140 cGy, or to PDDs of 87.8, 70.4, and 43.6, are used to calculate the OD for the lateral effect, OL,U, using the following equation:
where RD(y) is the profile dose that is obtained from the scan of the water phantom, normalized to that of the central beam axis at the depth of dose D (280, 225, or 140 cGy). The left part of the equation also applies to the normalized profile doses but is calculated using equation (4) with AnU, BnU, and CnU. The y-axis is parallel to the gantry-rotation axis with its origin at the isocenter. OL, U is calculated using the fitted results for a third-order polynomial form. When OL, U is obtained, the value of Ndd,U for the entire film is calculated using the equation.
Calibration in the red channel is more sensitive to a range of daily doses,,, and hence, the red channel of Ndd, U, Ndd, R, was obtained and this calibration is named “the red IT with the template.” To decrease the uncertainty further, a conventional fitting process with a power function was used with the equation:
where Odd, R= log10(Ndd, R); FR is the fitted dose for fitting parameters, A', B', C', and D'. The fitting process was repeated twice, with C'first bounded between 1 and 3, and then with C' and D' fixed to the value of the first fit.
For the dose calculation using Ndd, R, the total standard uncertainty, which is the combined uncertainty due to the experiment and the fitting processes, is calculated using equation (4) and equations in reference 26 with the fitted dose in percentage terms and the uncorrelated input quantities Ndd, R, AnR, BnR, and CnR. The calculation is performed based on the guidelines of the Joint Committee for Guides in Metrology (JCGM).
Similarly, the relative total standard uncertainty for the dose calculation using Odd,R, is calculated using equation (9) and equations in reference 26 according to the JCGM guidelines, using the uncorrelated input quantities Odd, R, A' and B'. The accuracy of the fitted dose and the uncertainty calculation using the traditional red channel net optical density (NOD) and the PDD method  is also calculated and compared with that for the R-3C method using the template.
| Results and Discussion|| |
In contrast to [Figure 1], when a template matrix is used and the DI OD (Odi) is eliminated using equations (3), (6), and (7), NU can be neatly fitted to the Dd using equation (4) for all three channels, as shown in [Figure 2]. The results demonstrate that the template matrix (equation 3) is required for the rational fitting form of equation (4). If it is not used, the difference between the fitted dose and the Dd could be >50%.
|Figure 2: A comparison of the fitting results for the inverse transmittances NU using the red-channel after three-channel method with a template matrix with those for the fitting data for IU in Figure 1|
Click here to view
If an Epson scanner is used, the red channel has a greater dynamic range than the green channel. In contrast, if the Microtek scanner is used, the green channel has a greater dynamic range than the red channel [Figure 2].
The fitted doses that are calculated using equations (4) and (6) in [Figure 3] better match the Dd than those that are calculated using the NOD of the red channel and the PDD method. All differences using the R-3C with a template are <2%, which are similar to the calibration results using an Epson scanner.,
|Figure 3: Differences between the fitted dose and the delivered dose, calculated using the red-channel after three-channel method with a template matrix and using the existing percentage-depth-dose method with red-channel net optical density|
Click here to view
[Figure 4] shows the OD calibration for lateral effect. It is seen that all of the calibration factors are approximately between 0.99 and 1.01, because the template matrix regulates inhomogeneity in the scanner. This calibration factor gives a more refined dose calculation.
|Figure 4: The optical density of lateral effect, OL, U, using equation (8) for doses of 280, 225, and 140 cGy: (a) The red channel, (b) The green channel, and (c) The blue channel|
Click here to view
When the DI OD and lateral effect are eliminated, the dose profiles for 280, 225, and 140 cGy are normalized to the dose at the central axis and then compared to the profile measured in the water phantom. The difference is <2% for all channels [Figure 5].
|Figure 5: Difference between the normalized dose profiles calculated after eliminating dose-independent and optical density for the lateral effect and the profile measured from the water phantom for the three channels: note that no data points are obscured by the legend|
Click here to view
The results of the uncertainty calculation are shown in [Figure 6]. It is seen that the total standard uncertainty relative to the fitted dose for the R-3C method using the template and the NOD calibration for the red channel is <10%. This value is high, when compared with the value of 5% for the Epson 10000XL scanner. However, the total uncertainty for red IT using the template is even higher and can be higher than 20% for Dd around 50 cGy.
|Figure 6: The total standard deviation, calculated using the red-channel after three-channel method and a template, the red-channel net optical density method and red inverse transmittance with a template|
Click here to view
| Conclusion|| |
This study demonstrates that a template matrix is required for the rational fitting process using equation (4). If it is not used, there is a significant error and the fitting cannot be performed. Using the proposed technique, the differences between the fitted dose and the Dd are <2%. The difference for the traditional one red-channel method is >3%. The calibration accuracy is similar to that for an Epson 10000XL scanner. The lateral effect can also be ignored because this calibration uses a template matrix.
This method means that the three-channel method can be used to calibrate film using a Microtek 9800XL + scanner and the template matrix can be useful for other brand scanners. The greater dynamic range of the green channel is a subject for future study. The green channel may replace the red channel for the verification of clinical daily treatment using a Microtek scanner.
This work was supported in part by the Ministry of Science and Technology of Taiwan (MOST 104-2221-E-214-014), the Research Support Scheme of the E-DA hospital (EDAHP106034) and Research Support Scheme of the I-Shou University (ISU-107-01-14A).
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Fuss M, Sturtewagen E, De Wagter C, Georg D. Dosimetric characterization of gafChromic EBT film and its implication on film dosimetry quality assurance. Phys Med Biol 2007;52:4211-25.
Niroomand-Rad A, Blackwell CR, Coursey BM, Gall KP, Galvin JM, McLaughlin WL, et al.
Radiochromic film dosimetry: Recommendations of AAPM radiation therapy committee task group 55. American Association of Physicists in Medicine. Med Phys 1998;25:2093-115.
Arjomandy B, Tailor R, Anand A, Sahoo N, Gillin M, Prado K, et al.
Energy dependence and dose response of gafchromic EBT2 film over a wide range of photon, electron, and proton beam energies. Med Phys 2010;37:1942-7.
Sutherland JG, Rogers DW. Monte Carlo calculated absorbed-dose energy dependence of EBT and EBT2 film. Med Phys 2010;37:1110-6.
Devic S, Tomic N, Soares CG, Podgorsak EB. Optimizing the dynamic range extension of a radiochromic film dosimetry system. Med Phys 2009;36:429-37.
Pérez Azorín JF, Ramos García LI, Martí-Climent JM. A method for multichannel dosimetry with EBT3 radiochromic films. Med Phys 2014;41:062101.
Ohuchi H. High sensitivity radiochromic film dosimetry using an optical common-mode rejection and a reflective-mode flatbed color scanner. Med Phys 2007;34:4207-12.
Aland T, Kairn T, Kenny J. Evaluation of a gafchromic EBT2 film dosimetry system for radiotherapy quality assurance. Australas Phys Eng Sci Med 2011;34:251-60.
Andrés C, del Castillo A, Tortosa R, Alonso D, Barquero R. A comprehensive study of the gafchromic EBT2 radiochromic film. A comparison with EBT. Med Phys 2010;37:6271-8.
McCaw TJ, Micka JA, Dewerd LA. Characterizing the marker-dye correction for gafchromic(®) EBT2 film: A comparison of three analysis methods. Med Phys 2011;38:5771-7.
Kairn T, Aland T, Kenny J. Local heterogeneities in early batches of EBT2 film: A suggested solution. Phys Med Biol 2010;55:L37-42.
Richley L, John AC, Coomber H, Fletcher S. Evaluation and optimization of the new EBT2 radiochromic film dosimetry system for patient dose verification in radiotherapy. Phys Med Biol 2010;55:2601-17.
Hartmann B, Martisiková M, Jäkel O. Homogeneity of gafchromic EBT2 film. Med Phys 2010;37:1753-6.
Chang L, Ho SY, Lee TF, Yeh SA, Ding HJ, Chen PY. The suitable dose range for the calibration of EBT2 film by the PDD method with a comparison of two curve fitting algorithms. Nucl Instrum Methods Phys Res A 2015;777:85-90.
Moylan R, Aland T, Kairn T. Dosimetric accuracy of gafchromic EBT2 and EBT3 film for in vivo
dosimetry. Australas Phys Eng Sci Med 2013;36:331-7.
Brown TA, Hogstrom KR, Alvarez D, Matthews KL 2nd
, Ham K, Dugas JP. Dose-response curve of EBT, EBT2, and EBT3 radiochromic films to synchrotron-produced monochromatic X-ray beams. Med Phys 2012;39:7412-7.
Desroches J, Bouchard H, Lacroix F. Potential errors in optical density measurements due to scanning side in EBT and EBT2 gafchromic film dosimetry. Med Phys 2010;37:1565-70.
Lewis D, Micke A, Yu X, Chan MF. An efficient protocol for radiochromic film dosimetry combining calibration and measurement in a single scan. Med Phys 2012;39:6339-50.
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:5257-62.
Marroquin EY, Herrera González JA, Camacho López MA, Barajas JE, García-Garduño OA. Evaluation of the uncertainty in an EBT3 film dosimetry system utilizing net optical density. J Appl Clin Med Phys 2016;17:466-81.
Fiandra C, Ricardi U, Ragona R, Anglesio S, Giglioli FR, Calamia E, et al.
Clinical use of EBT model gafchromic film in radiotherapy. Med Phys 2006;33:4314-9.
Hupe O, Brunzendorf J. A novel method of radiochromic film dosimetry using a color scanner. Med Phys 2006;33:4085-94.
Mathot M, Sobczak S, Hoornaert MT. Gafchromic film dosimetry: Four years experience using filmQA Pro software and Epson flatbed scanners. Phys Med 2014;30:871-7.
Micke A, Lewis DF, Yu X. Multichannel film dosimetry with nonuniformity correction. Med Phys 2011;38:2523-34.
van Hoof SJ, Granton PV, Landry G, Podesta M, Verhaegen F. Evaluation of a novel triple-channel radiochromic film analysis procedure using EBT2. Phys Med Biol 2012;57:4353-68.
Chang L, Ho SY, Lee TF, Yeh SA, Ding HJ, Chen PY. Calibration of EBT2 film using a red-channel PDD method in combination with a modified three-channel technique. Med Phys 2015;42:5838-47.
Chang L, Ho SY, Ding HJ, Hwang IM, Chen PY, Lee TF. Evaluation of multiple-sampling function used with a microtek flatbed scanner for Radiation dosimetry calibration of EBT2 film. Nucl Instrum Methods Phys Res A 2016;832:179-83.
Chang L, Chui CS, Ding HJ, Hwang IM, Ho SY. Calibration of EBT2 film by the PDD method with scanner non-uniformity correction. Phys Med Biol 2012;57:5875-87.
Klein EE, Hanley J, Bayouth J, Yin FF, Simon W, Dresser S, et al.
Task group 142 report: Quality assurance of medical accelerators. Med Phys 2009;36:4197-212.
Kutcher GJ, Coia L, Gillin M, Hanson WF, Leibel S, Morton RJ, et al.
Comprehensive QA for radiation oncology: Report of AAPM radiation therapy committee task group 40. Med Phys 1994;21:581-618.
Schulz RJ, Almond PR, Cunningham JR, Holt JG, Loevinger R, Suntharalingam NI, et al
. A protocol for the determination of absorbed dose from high-energy photon and electron beams. Med Phys 1983;10:741-71.
Almond PR, Biggs PJ, Coursey BM, Hanson WF, Huq MS, Nath R, et al.
AAPM's TG-51 protocol for clinical reference dosimetry of high-energy photon and electron beams. Med Phys 1999;26:1847-70.
Chang L, Ho SY, Ding HJ, Lee TF, Chen PY. Dependency of EBT2 film calibration curve on postirradiation time. Med Phys 2014;41:021726.
Joint Committee of Guides in Metrology. Evaluation of Measurement Data – Guide to the Expression of Uncertainty in Measurement. Vol. 100. Joint Committee of Guides in Metrology; 2008.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6]