Comparison of 3D anatomical dose verification and 2D phantom dose verification of IMRT/VMAT treatments for nasopharyngeal carcinoma
© Lin et al.; licensee BioMed Central Ltd. 2014
Received: 16 December 2013
Accepted: 2 March 2014
Published: 7 March 2014
The two-dimensional phantom dose verification (2D-PDV) using hybrid plan and planar dose measurement has been widely used for IMRT treatment QA. Due to the lack of information about the correlations between the verification results and the anatomical structure of patients, it is inadequate in clinical evaluation. A three-dimensional anatomical dose verification (3D-ADV) method was used in this study to evaluate the IMRT/VMAT treatment delivery for nasopharyngeal carcinoma and comparison with 2D-PDV was analyzed.
Twenty nasopharyngeal carcinoma (NPC) patients treated with IMRT/VMAT were recruited in the study. A 2D ion-chamber array was used for the 2D-PDV in both single-gantry-angle composite (SGAC) and multi-gantry-angle composite (MGAC) verifications. Differences in the gamma pass rate between the 2 verification methods were assessed. Based on measurement of irradiation dose fluence, the 3D dose distribution was reconstructed for 3D-ADV in the above cases. The reconstructed dose homogeneity index (HI), conformity index (CI) of the planning target volume (PTV) were calculated. Gamma pass rate and deviations in the dose-volume histogram (DVH) of each PTV and organ at risk (OAR) were analyzed.
In 2D-PDV, the gamma pass rate (3%, 3 mm) of SGAC (99.55% ± 0.83%) was significantly higher than that of MGAC (92.41% ± 7.19%). In 3D-ADV, the gamma pass rates (3%, 3 mm) were 99.75% ± 0.21% in global, 83.82% ± 16.98% to 93.71% ± 6.22% in the PTVs and 45.12% ± 32.78% to 98.08% ± 2.29% in the OARs. The maximum HI increment in PTVnx was 19.34%, while the maximum CI decrement in PTV1 and PTV2 were -32.45% and -6.93%, respectively. Deviations in dose volume of PTVs were all within ±5%. D2% of the brainstem, spinal cord, left/right optic nerves, and the mean doses to the left/right parotid glands maximally increased by 3.5%, 6.03%, 31.13%/26.90% and 4.78%/4.54%, respectively.
The 2D-PDV and global gamma pass rate might be insufficient to provide an accurate assessment for the complex NPC IMRT operation. In contrast, the 3D-ADV is superior in clinic-related quality assurance offering evaluation of organ specific pass rate and dose-volume deviations.
KeywordsIMRT VMAT 3D anatomical dose 2D phantom dose Dosimetry verification
Intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) techniques are able to provide very high dose conformity for cancer radiotherapy; thus, the surrounding normal tissue and organs can be well protected when high-dose radiation is delivered to the target volume. However, many uncertainties exist in the treatment planning and operation process that can lead to deviations of the IMRT or VMAT dose distribution. Therefore, it is necessary to verify the irradiation dose distribution that is delivered by the accelerator before such kinds of treatments . Until recently, a hybrid plan has been adopted most often in 2-dimensional (2D) planar dose measurement verifications. The precision of irradiating doses is evaluated and verified by comparing the planned dose distribution, calculated by the treatment planning system (TPS), to the measured results . However, due to the lack of information about the correlations between the verification measurement results and the anatomical structure of patients, as well as the resulting lack of information about the actual irradiation doses to different target volumes and organs at risk (OARs), it is difficult to identify the geometric locations where dose errors occur during plan implementation, thus leading to inadequate clinical evaluation information . Recently, some 3-dimensional (3D) dose verification tools that provide patient anatomical structure information were applied clinically. These tools can provide important information such as the dose deviations, the pass rates and the locations of the dose deviations in the patients’ target volumes and organs, as well as identification of the error origins [4, 5]. In this study, we adopted a 3D anatomical dose verification (3D-ADV) based on measurements of delivered dose fluence and patients’ anatomical images. Meanwhile, the traditional 2D phantom verification (2D-PDV) using an ionization chamber array with angular response correction applied an in-house software, were used to compare the efficacies of dose verification for IMRT and VMAT of nasopharyngeal carcinoma (NPC). As a result, the differences between the two verification methods and their clinical significances in evaluations of irradiating dose deviations were analyzed and clarified.
Materials and methods
Dose verification tools
A commercialized 3D dose verification system  (COMPASS, IBA Dosimetry, Schwarzenbruck, Germany) was used. The system included an online 2-dimensional ion-chamber array (2D-IC array) and dose reconstruction software based on a beam model describing the characteristics of the accelerator (e.g., energy spectrum, lateral beam quality variations) and the collapsed-cone convolution/superposition (CC) algorithm, which computed the radiation dose distribution on to the patient CT image set. A strict commissioning of the whole system, including the validation of accuracy for 2D-IC array measurement, beam modeling and dose reconstruction, was performed in advance according to the same standards as the clinic used TPS. In this previously commissioning job, which had been published , phantom plans of regular, irregular fields and IMRT were selected tested. All test plans were implemented and the dose distributions were measured using thimble ion-chamber and the 2D-IC array, the accuracy of the 3D-ADV system were then evaluated by comparing the corresponding measurement results. The average deviation of the 3D-ADV system was less then 1% with the largest difference of 2.12% comparing to the thimble ion-chamber measurements. In the comparison of IMRT planning calculation and the 3D-ADV system computation, the global and PTV gamma pass rates were better than 99% and 98%, respectively.
The 2D-IC array (MatriXX, IBA Dosimetry) used in this study comprised 1020 plane-parallel ion-chambers arranged with a distance of 0.762 cm between the chamber centers. A build-up of 3.3-mm equivalent water thickness was placed on top of the chambers, and a 22-mm-thick RW3 backscattering phantom (composition: 98% Polystyrol, 2% TiO2 and density: 1.045 g/cm3) was placed on the back. The array was inserted into an additional phantom with exterior dimensions of 31.4 cm × 34 cm × 22 cm (MULTICube Lite, IBA Dosimetry) for multi-gantry-angle composite 2D dose verification, equipped with an angular position sensor for angular response corrections.
Planning and delivery system
The IMRT and VMAT plans measured in this study were designed by a 3D inverse treatment planning system (Monaco 3.0, Elekta, Stockholm, Sweden). The planned dose distribution was calculated using a fast Monte Carlo algorithm (XVMC) with the variance rate of 3% and grid size of 3 mm, as same as the reported study of Boggula et al. . A 6MV X-ray irradiation was delivered via a linear accelerator (Synergy VMAT, Elekta) that supported step-and-shoot IMRT and VMAT therapies.
Methods of dose verification
The above-mentioned MatriXX 2D-IC array with the MULTICube phantom was used to perform the 2D-PDV for 20 NPC treatment cases, 10 cases treated with IMRT and 10 cases with VMAT. The tests were performed in single-gantry-angle composite (0°, SGAC) and multi-gantry-angle composite (MGAC) methods, respectively . Different from the published method of angle dependence correction for MatriXX , an in-house correction software program was used to compute the angular correction for each independent ion-chamber of the 2D-IC array in the MGAC measurements. A comparison was performed to identify the differences in the gamma pass rates (3%, 3 mm)  at the isocenter plane from the 2 composite dose verifications.
Angle correction of the ion-chamber detector
Within a gantry angle range from 0° to 360° with a step size of 5° (a step size of 1° within the range of 90° ± 5° and 270° ± 5°), the dose at the point of central detector was calculated by TPS and measured when the gantry was maintained at the same position. The correction factor of the central detector at each angle was then obtained by comparing the difference between the measurement and TPS calculation, and was normalized to the one at zero degree gantry angle. These correction factors can be applied to all the other ion-chambers in the array assuming the angular response coincidence is high among all of them. Due to the differences in the incident angles of the central and the other ion-chambers, the incident angle of each chamber was calculated with the following formula. For an angle other then the above tested ones, a linear interpolation method was adopted to get the corresponding angular correction factor. Visual C++ language was used to program this angle correction for the measurement results.
3D-ADV test was done for the same 20 NPC IMRT/VMAT plans as well, with the 3D dose verification system that passed the commissioning test in advance as briefly described above, similar with that reported by Boggula et al. . During the plan delivery, the irradiated fluence of every radiation beam was measured by the 2D-IC array mounted on a holder attached to the radiation head, perpendicular to the radiation beam. According to the measured results, the dose distribution calculation was corrected and the final dose reconstruction (RDD) was obtained for each plan case; subsequently, the RDD results were compared to the original treatment plan. Parameters such as the homogeneity index (HI) and the conformity index (CI) of the planning target volumes for primary nasopharynx tumor (PTVnx), high-risk subclinical region (PTV1), and the preventive irradiation region (PTV2), the dose volume parameters of every PTV and OAR and the gamma pass rate (3%, 3 mm) for global and for each organ were analyzed.
D2%, D50% and D98% indicate the doses that covered 2% (near-maximum dose), 50%, and 98% (near-minimum dose) of the PTV, respectively .
SPSS 18.0 software (Chicago, IL, USA) was used for 1 sample T-tests of the dose verification data from the 20 cases of IMRT/VMAT test plans. A value of P < 0.05 was considered statistically significant.
Results from the SGAC/MGAC 2D-PDV
The gamma pass rates (3%, 3 mm) of the 0° SGAC and MGAC in 2D-PDV
Matrixx build-in ang corr
In-house ang corr
Mean pass rate ± σ
99.55% ± 0.83%
89.22% ± 10.61%
92.41% ± 7.19%
Comparison of 3D-ADV and TPS planning
The mean gamma pass rates (3%, 3 mm) of global volume and selected organs in 3D-ADV for the 20 NPC IMRT/VMAT plans (Std indicated the standard deviation)
Mean pass rate ± Std.
99.75 ± 0.21
83.82 ± 16.98
90.68 ± 9.34
93.71 ± 6.22
96.93 ± 4.58
65.69 ± 20.54
97.33 ± 3.72
98.08 ± 2.29
The deviations in target volume HI and CI between the RDD and the plan value
-0.94 ± 9.24
-3.21 ± 10.53
-1.92 ± 3.00
The relative deviations in dose-volume histogram (DVH) parameters between those obtained from the RDD and those expected by the therapy plan
0.93 ± 1.88
0.93 ± 1.52
1.07 ± 1.85
0.44 ± 1.72
0.97 ± 2.17
1.07 ± 1.44
1.01 ± 2.06
0.04 ± 0.35
0.21 ± 1.43
1.04 ± 1.57
0.54 ± 1.48
-0.05 ± 0.19
-0.90 ± 2.22
1.56 ± 2.03
-0.21 ± 11.80
Optic nerve (Left)
0.34 ± 9.49
Optic nerve (Right)
-0.38 ± 9.08
Parotid gland (Left)
-0.55 ± 2.71
Parotid gland (Right)
0.40 ± 2.00
An effective evaluation of the treatment operation in clinical practice should be able to reflect the true delivery condition of the treatment and any errors occurred in the planned parameters. The traditional 2D-PDV QA, especially the SGAC measurement which used only a single fixed incident angle to avoid the existence of angular response errors, can result in inconsistencies between the QA and clinical therapy conditions , and not able to provide information regarding the relationship between the dose error and the anatomical structures. This leads to decreasing of QA abilities in clinical evaluations . In our study of 2D-PDV of 20 NPC IMRT/VMAT plans, the gamma pass rates obtained from the SGAC verification were significantly higher than the results from the MGAC measurements done with the same gantry angle of therapy. This might be explained by the fact that the SGAC verification cannot reflect the effects on the radiation dose caused by gravity-induced changes in the multi-leaf collimator (MLC) position, or the output dose angle dependence from the accelerator under different gantry angles. The results indicated that the pass rates corrected by MGAC verification with the MatriXX build-in angle correction factor were lower than those corrected by our in-house correction software. The reason could be that the correction factor of the former was determined according to the correction angle of the central detector and then used to correct all of the other chambers without considering the differences in the incident angles of each detector. The angle correction from the in-house software considered the influence of different incident angles of each detector in the 2D-IC array; therefore, the pass rates in the verification were significantly superior to the results without the independent incident angle correction.
3D-ADV can provide us with information such as the pass rates (the global, each target and OAR volumes in the measured area), the statistical results of deviations in the dose-volume histogram (DVH) parameters (including the dose volume and the volume dose) of each organ and the anatomical positions that correspond to the dose deviations. The 3D-ADV results from the 20 nasopharyngeal carcinoma patients, who received IMRT/VMAT irradiation, revealed that the mean global gamma pass rate was 99.75% ± 0.21%. Lee et al.  reported that the mean gamma pass rates (3%, 3 mm) of 2D-PDV were 98.2% ± 1.3% and 98.5% ± 1.3% in nasopharyngeal carcinoma patients treated with IMRT and VMAT, respectively. Our results from the 2D-PDV with SGAC and MGAC on a similar group of patients were 99.55% ± 0.83% and 92.41% ± 7.19%, respectively. The global gamma pass rate from the 3D dose verification was similar with that from the 2D-PDV of the SGAC method but was higher than that from the 2D-PDV of MGAC measurement. The reason for the higher global pass rate in the 3D dose verification could be that, compared to 2D dose verification, the dose pixels evaluated in the 3D dose verification included all points within the CT scanning area, thus resulting in a relatively lower ratio of pixel numbers at which dose deviations occurred to the overall pixel number. Additionally, during the implementation of MGAC 2D-PDV, the pass rates might decrease due to the non-uniform phantom density that can result from the detector arrangement when the incident angles are parallel to the detector plane . The gamma pass rates of each PTV and OAR decreased at different levels in all cases when compared to the global pass rate, indicating that higher gamma pass rates might be caused by improper evaluation strategies in which some errors in the delivered dose distribution were disguised and ignored due to the use of the global pass rate in the evaluation.
In this study of 3D-ADV, no statistically significant difference was found in the HI of the PTVnx between the measurement based dose reconstruction and the planned value of the TPS. However, there were obvious alterations in the verification results for each patient that the HI deviation ranged from -15.62% to 19.34%, indicating that there were significant individual differences in the irradiation results. Some patients had greater HI values in the PTVnx with irradiated doses than the planned values, indicating decreased homogeneity of the dose distribution to the target volumes. Even greater deviations from the planned value were observed in the CI of the PTV1, which were the high-risk lymph nodular target volume, ranged from -32.45% to 13.97%, indicating that the dose conformity of PTV1 in some individual cases decreased a lot after the plan delivery, which might have led to inadequate dose coverage in this target volume or to elevated doses in the surrounding normal tissues. The CI of the PTV2 decreased by an average of approximately 2% (P < 0.05), thus showing a total reduction of the PTV2 dose conformity during the implementation of the plan.
Elmpt et al.  reported a work of similar 3D-ADV, using the planning CT image to reconstruct the dose distribution in combination with a Monte Carlo calculation and the energy fluence of the actual treatment beams measured pre-treatment with the electronic portal image device (EPID). In their study on head and neck cancers IMRT cases, most deviations between the reconstructed delivered dose and the planned value in the PTV, including D5%, Dmean and D95%, were less than 3%, while the mean dose in the parotid gland decreased by 3.2% ± 1.2% and the maximum spinal cord dose increased by 3.1% ± 1.9%. Our 3D verification data showed that the deviations in the D2%, D98% and D95% of the PTVs and the deviations in the V100% and V95% were all within ±5%, which were similar to the previous report. In contrast, the maximum deviations of D2% in the optic nerves and the optic chiasm were 31% and 27%, respectively, suggesting that extra attention should be given to future plan verifications with regard to whether the delivered dose to these organs will exceed the clinical limits. In the reported study of Stasi  and Carrasco , a weak correlations or even no correlation has been found between the gamma index and the clinical impact of a delivery dose discrepancy, like the deviation on DVH for PTV and OAR volumes. The acceptance criteria for which we had the highest frequency of correlations were (3%, 3 mm), however, this criterion hid relevant clinical dose metric differences which is not clinically acceptable. Our study also shown that there might be clinically unacceptable dose discrepancy in some cases even the (3%, 3 mm) gamma pass rate was very high.
Traditional 2D-PDV is relatively simple and easy to perform. However, this method cannot provide information regarding anatomical positions and the dose volumes that correspond to the dose deviations; thus, it can only serve as a basic quality assurance tool for IMRT. using online measured results, 3D-ADV is able to reconstruct the dose distributions from patients' anatomical images and provide us with more clinical relative information, verify the delivery deviation in both dosimetric and geometric parameters in the same way of plan evaluation. In the other hand, unlike 2D-PDV devices, the 3D-ADV utilized an independent calculation algorithm in the dose reconstruction, which differed from that used in the planning system. Therefore, it might lead to additional discrepancy in the verification results, if the dose reconstruction computation of the 3D-ADV was not accurate enough. It is important that a very strict pre-commissioning and proper evaluation has been performed when the 3D-ADV system is used for the clinic treatment plan QA.
In this study of 20 nasopharyngeal carcinoma patients who were treated with IMRT or VMAT, comparing the results of the delivered dose distributions from traditional 2D-PDV and the 3D-ADV, it was confirmed that a relatively large local dose deviation might exist in the delivered dose distributions, even when the global gamma pass rate is very high in both verification. A 3D-ADV, providing structure by structure volumetric dose evaluation is suggested to be used as further clinical quality assurances for IMRT/VMAT therapy of complex cases like nasopharyngeal carcinoma.
- Hendee WR, Herman MG: Improving patient safety in radiation oncology. Med Phys 2011,38(1):78-82. 10.1118/1.3522875View ArticlePubMedGoogle Scholar
- Low DA, Moran JM, Dempsey JF, Dong L, Oldham M: Dosimetry tools and techniques for IMRT. Med Phys 2011,38(3):1313-1338. 10.1118/1.3514120View ArticlePubMedGoogle Scholar
- Nelms BE, Zhen HM, Tome WA: Per-beam planar IMRT QA passing rates do not predict clinically relevant patient dose errors. Med Phys 2011,38(2):1037-1043. 10.1118/1.3544657PubMed CentralView ArticlePubMedGoogle Scholar
- Boggula R, Jahnke L, Wertz H, Lohr F, Wenz F: Patient-specific 3D pretreatment and potential 3D online dose verification of Monte Carlo-calculated IMRT prostate treatment plans. Int J Radiat Oncol Biol Phys 2011,81(4):1168-1175. 10.1016/j.ijrobp.2010.09.010View ArticlePubMedGoogle Scholar
- McCurdy , Mueller L, Backman E, Venkataraman S, Fleming E, Asuni G, Jensen M, Ur-Rehman F, Pistorius S: Commissioning and validation of a novel measurement-based IMRT QA method, incorporating dose recalculation on patient CT data. Med Phys 2008,35(6):2760-2760.View ArticleGoogle Scholar
- Godart J, Korevaar EW, Visser R, Wauben DJL, Van't Veld AA: Reconstruction of high-resolution 3D dose from matrix measurements: error detection capability of the COMPASS correction kernel method. Phys Med Biol 2011,56(15):5029-5043. 10.1088/0031-9155/56/15/023View ArticlePubMedGoogle Scholar
- Lin HL, Huang SM, Deng XW, Jin GH, Sun WZ, Yao XH, Zhang DD, Chen LX: A clinical test and application research of IMRT dose verification system based on patient's anatomical structure and on-line dosimetry. Chin J Radiat Oncol 2012,21(3):271-275.Google Scholar
- Nelms BE, Simon JA: A survey of planar IMRT QA analysis. J Appl Clin Med Phys 2007,8(3):1-15.View ArticleGoogle Scholar
- Shimohigashi Y, Araki F, Tominaga H, Sakata J, Kawasaki K, Kanetake N, Iwashita Y, Yoshimura S, Kawakami T, Ishihara T, Okuda T, Kogo K: Angular dependence correction of MatriXX and its application to composite dose verification. J Appl Clin Med Phys 2012,13(5):198-214.Google Scholar
- Low DA, Harms WB, Mutic S, Purdy JA: A technique for the quantitative evaluation of dose distributions. Med Phys 1998,25(5):656-661. 10.1118/1.598248View ArticlePubMedGoogle Scholar
- Boggula R, Lorenz F, Mueller L, Birkner M, Wertz H, Stieler F, Steil V, Lohr F, Wenz F: Experimental validation of a commercial 3D dose verification system for intensity-modulated arc therapies. Phys Med Biol 2010,55(19):5619-5633. 10.1088/0031-9155/55/19/001View ArticlePubMedGoogle Scholar
- Feuvret L, Noël G, Mazeron JJ, Bey P: Conformity index: a review. Int J Radiat Oncol Biol Phys 2006,64(2):333-342. 10.1016/j.ijrobp.2005.09.028View ArticlePubMedGoogle Scholar
- ICRU report 83: Prescribing, recording, and reporting intensity-modulated photon-beam therapy. J ICRU 2010,10(1):1-106.View ArticleGoogle Scholar
- Wolfsberger L, Wagar M, Nitsch P, Bhagwat M, Zygmanski P: Angular dose dependency of matrixx TM and its calibration. J Appl Clin Med Phys 2010,11(1):241-251.Google Scholar
- Lee TF, Chao PJ, Ting HM: Comparative analysis of SmartArc-based dual arc volumetric-modulated arc radiotherapy (VMAT) versus intensity-modulated radiotherapy (IMRT) for nasopharyngeal carcinoma. J Appl Clin Med Phy 2011,12(4):158-174.Google Scholar
- Li QL, Deng XW, Chen LX, Huang XY, Huang SM: The angular dependence of a 2-dimensional diode array and the feasibility of its application in verifying the composite dose distribution of intensity-modulated radiation therapy. Chin J Canc 2010,29(6):617-620. 10.5732/cjc.009.10592View ArticleGoogle Scholar
- Elmpt WV, Nijsten S, Mijnheer B, Dekker A, Lambin P: The next step in patient-specific QA: 3D dose verification of conformal and intensity-modulated RT based on EPID dosimetry and Monte Carlo dose calculations. Radiother Oncol 2008,86(1):86-92. 10.1016/j.radonc.2007.11.007View ArticlePubMedGoogle Scholar
- Stasi M, Bresciani S, Miranti A, Maggio A, Sapino V, Gabriele P: Pretreatment patient-specific IMRT quality assurance: a correlation study between gamma index and patient clinical dose volume histogram. Med Phys 2012,39(12):7626-7634. 10.1118/1.4767763View ArticlePubMedGoogle Scholar
- Carrasco P, Jornet N, Latorre A, Eudaldo T, Ruiz A, Ribas M: 3D DVH-based metric analysis versus per-beam planar analysis in IMRT pretreatment verification. Med Phys 2012,39(8):5040-5049. 10.1118/1.4736949View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.