 Research
 Open Access
 Published:
Fast 3D dosimetric verifications based on an electronic portal imaging device using a GPU calculation engine
Radiation Oncology volume 10, Article number: 85 (2015)
Abstract
Purpose
To use a graphic processing unit (GPU) calculation engine to implement a fast 3D pretreatment dosimetric verification procedure based on an electronic portal imaging device (EPID).
Methods
The GPU algorithm includes the deconvolution and convolution method for the fluencemap calculations, the collapsedcone convolution/superposition (CCCS) algorithm for the 3D dose calculations and the 3D gamma evaluation calculations. The results of the GPUbased CCCS algorithm were compared to those of Monte Carlo simulations. The planned and EPIDbased reconstructed dose distributions in overriddentowater phantoms and the original patients were compared for 6 MV and 10 MV photon beams in intensitymodulated radiation therapy (IMRT) treatment plans based on dose differences and gamma analysis.
Results
The total singlefield dose computation time was less than 8 s, and the gamma evaluation for a 0.1cm grid resolution was completed in approximately 1 s. The results of the GPUbased CCCS algorithm exhibited good agreement with those of the Monte Carlo simulations. The gamma analysis indicated good agreement between the planned and reconstructed dose distributions for the treatment plans. For the target volume, the differences in the mean dose were less than 1.8%, and the differences in the maximum dose were less than 2.5%. For the critical organs, minor differences were observed between the reconstructed and planned doses.
Conclusions
The GPU calculation engine was used to boost the speed of 3D dose and gamma evaluation calculations, thus offering the possibility of true realtime 3D dosimetric verification.
Background
As the frequency of dose prescriptions and the complexity of radiotherapy techniques increase, so too do the demands for accurate and efficient methods of verifying the doses delivered to patients. Hence, dosimetric verification is a prerequisite to ensure correct treatment planning and delivery. In recent years, an increasing number of reports have begun to focus on the 3D dosimetric verification of intensitymodulated radiation therapy (IMRT) treatment plans [13]. Most planar dosimetric verifications are based on uniform phantoms, but this method cannot be used to determine the algorithm error of a treatment planning system (TPS) with respect to inhomogeneities. Nelms et al. have reported that there is a lack of correlation between conventional IMRT quality assurance (QA) performance metrics (gamma pass rate) and dose errors in anatomical regions of interest. There is no clear relation between the gamma pass rate and the error in dose–volume histograms (DVHs) [4]. 3D dosimetric verifications can yield more intuitive anatomybased 3D dose differences for clinical verifications and other detailed evaluation indices, such as the DVH.
Because of their submillimeter resolution, realtime response and lower workload, electronic portal imaging devices (EPIDs) are often utilized for pretreatment and in vivo measurements [1,2,57]. There are two approaches to using an EPID for 3D dosimetric verifications. One technique is to calculate the 3D dose distributions based on the EPID measurement using a TPS. Steciw et al. have reported the use of a TPS to calculate dose distributions based on the fluences of treatment fields that were measured using an EPID [8]. This method can be easily performed in clinical practice because no effort is required for the independent dosecalculation algorithm and model commissioning. However, the effect of the TPS calculation errors cannot be estimated. Moreover, not all commercial TPSs have functions to export fluence maps or to calculate doses based on input fluence maps. The other technique is to use an independent dosecalculation algorithm. Wendling et al. have reported the use of the “backprojection” algorithm for this purpose [9]. This method is simple and has been used clinically in their institution [6]. The model is “waterequivalent,” as the parameters are determined under homogeneous conditions. For verifications of lung treatments in which inhomogeneities are important, a correction method known as the “in aqua vivo” method is applied [10]. Elmpt et al. have implemented 3D dose verifications based on Monte Carlo simulations. If implemented correctly, such simulations can provide, in principle, the highest accuracy that is currently available [2]. Recent publications have reported the implementation of a pencilbeam algorithm for fast 3D dosimetric verifications [1113]. Because of the simplicity of this pencilbeam algorithm, its computation speed is very fast. However, its limitations cause it to be unsuitable for application to inhomogeneous tissues [14].
In clinical applications, whether for pretreatment or for in vivo dosimetric verifications, two aspects must be considered: first, the computation time must be considered, with the potential goal of realtime verification, and second, the accuracy in inhomogeneous tissues must be examined. The collapsedcone convolution/superposition (CCCS) [15] algorithm is often used in radiotherapy TPSs. The calculation accuracy of the CCCS algorithm for inhomogeneous tissues is acceptable. A recent report demonstrated that the results produced by the CCCS dose engine that is used in the TomoTherapy® TPS are consistent with Monte Carlo results within 1%/1 mm in water and 2%/2 mm in the lung, except in cases involving a small field (1.25 cm × 1.25 cm) and a tissue–air interface [16]. Thus, the use of the CCCS algorithm is a compromise between the computation time and the accuracy in inhomogeneous tissues.
Although the CCCS algorithm incorporates the collapsedcone approximation, which allows the calculation complexity to be reduced in comparison with fullvolume pointtopoint convolution, it is still computationally demanding, especially for IMRT, which involves many fields or segments. Graphic processing units (GPUs) are increasingly being used for scientific applications, including 3D dosimetric calculations, that can be accelerated by exploiting the parallel architecture and unprecedented computing power density of GPUs [1720]. Quan et al. have reported that a GPU algorithm using exponential cumulative–cumulative kernels (CCKs) is 1000–3000 times faster than a highly optimized singlethreaded CPU implementation using a tabulated CCK [19]. Computation times on the order of seconds can be achieved. Thus, a GPUbased CCCS algorithm demonstrates good prospects for achieving true realtime 3D dosimetric verifications.
An EPID can be used to determine the delivery errors of treatment plans, and an accurate independent dosecalculation algorithm can be used to estimate the errors of TPS dose calculations, including scenarios involving the treatment of inhomogeneous sites such as lung tissues and oralnasal cavities, and to crosscheck the TPS modeling, which is strongly dependent on the input measurement data. The computation speed can be boosted via GPU implementation. The GPU implementation of fast EPIDbased 3D pretreatment dosimetric verifications is presented in this paper. The GPU algorithm includes the deconvolution and convolution method for the fluencemap calculations, the CCCS algorithm for the 3D dose calculations and the 3D gamma evaluations. The combination of an EPIDbased approach and a GPUbased approach can effectively reduce the measurement workload and calculation time that are required for pretreatment verifications. Furthermore, realtime IMRT dosimetric verifications and doseguided radiation therapy pose stringent requirements in terms of measurement and calculationtime capabilities. The combination of features provided by EPIDs and GPUs can satisfy these requirements.
Methods
EPID image acquisition
A Trilogy linear accelerator system (Varian Medical Systems, Palo Alto, CA, USA) with an aS1000 EPID (Varian Medical Systems, Palo Alto, CA, USA) was employed in this study. Both 6 MV and 10 MV photon beams were used for all measurements. The physical size of the EPID was as follows: the EPID had a sensitive area of 40 cm × 30 cm in size, and the effective pixel size was 0.04 cm × 0.04 cm. The pretreatment EPID measurement was performed in air without a phantom/patient. Because the EPID images were used as the input for subsequent dose calculations, it was necessary to completely capture the treatment fields in the crossplane direction. The effective sourcetodetector distance was set to 140 cm. Image acquisition was performed in integrated mode, and offset correction, gain correction and pixel correction were performed for each image. Because of the robotic arm that is located directly beneath the sensitive area of the Varian EPID, the backscatter from the arm can have a deleterious effect when the EPID is used for dosimetric purposes [2123]. The backscatter correction kernel was modeled as a radially symmetric Gaussian and applied to the raw EPID data [24]. To ensure positioning accuracy during image acquisition, a gantryangledependent correction associated with the panel position displacement was applied [25]. The EPID displacement was measured by comparing the image acquired at each gantry angle to the 0° image in the 10 cm × 10 cm field.
Fluencemap calculation
The pixel values of the EPID images, which captured the treatment fields during air delivery, were reconstructed into fluence maps at a source–axis distance (SAD) of 100 cm using the deconvolution and convolution method.
In our previous studies [26], the EPID, which consists of several layers, was treated as a uniform waterequivalent phantom. The relationship between the pixel values of the EPID images and the dose values on the central plane of the 20cmthick virtual water phantom, which was placed at a sourcetosurface distance (SSD) of 90 cm, can be represented as follows:
where ⊗ and ⊗ ^{− 1} denote the convolution and deconvolution operators, respectively. Every pixel in the EPID image is referred to by its indices i and j; P _{ ij } is the pixel value of the EPID image that is projected to the phantom level, D _{ ij } is the dose at the phantom level, c _{ad} is the absolute dose calibration factor, \( {K}_{ij}^{\mathrm{EPID}} \) is the EPID scatter kernel that is used to obtain the incident fluence map, and \( {K}_{ij}^{\mathrm{Water}} \) is the convolution kernel that is used to calculate the dose distribution in the water phantom.
In this study, this model was used to reconstruct the EPID image into a fluence map; thus, equation (1) was simplified as follows:
with
Here, d _{ ij } is the distance from the center of the kernel to a pixel ij. μ _{S} and μ _{L} are the attenuation coefficients, which depend on the energy and the material. There are two terms in this equation. The first term, which contains μ _{S}, is the shortrange primary dose component. It describes the energy per unit mass imparted by the charged particles released as a result of the first interaction of a primary photon. The second term, which contains μ _{L}, is the longrange scatter dose component. It describes the energy imparted by charged particles set in motion by secondary photons, i.e., by scattered, bremsstrahlung and annihilation photons.
Because of the broad distribution of photon energies within the phantom and the EPID, there will be an effective attenuation coefficient for the energy spectrum. Instead of the use of multiple monoenergetic kernels based on the beam’s energy spectrum, μ _{S} and μ _{L} were treated as free parameters in the fit procedure. For each nominal energy of the linear accelerator, the parameters of \( {K}_{ij}^{\mathrm{EPID}} \) are different. These parameters were determined using the central point dos × es, which were measured using an ion chamber in air for fields of 3 cm × 3 cm to 20 cm × 20 cm. The parameters of the \( {K}_{ij}^{\mathrm{EPID}} \) kernel that yielded the minimum variance could thus be determined via the goldensection search algorithm [27]. Because the structure of an EPID is invariant, these parameters were constants once the fit procedure had been completed.
A simple exponential kernel was introduced to reduce the noise arising from the application of the deconvolution operator, and a Gaussian function was used to correct the “horns” from the beam profile.
To improve the calculation speed, a fast Fourier transform (FFT) was used. The FFT was applied using the CUFFT library from the NVIDIA® CUDA™ library, which is highly optimized for computing parallel FFTs using an NVIDIA GPU.
The EPIDbased fluence maps for fields of 3 cm × 3 cm to 20 cm × 20 cm were compared to the ionchamber results. The S_{c} factor and the dose profiles in air were obtained using a 0.13cc ion chamber (IBA Dosimetry GmbH, Schwarzenbruck, Germany) with a buildup cap by scanning an empty 3D water tank (Blue Phantom, IBA Dosimetry GmbH, Schwarzenbruck, Germany) at an SAD of 100 cm.
3D dose reconstruction using a GPU
3D dosecalculation algorithm
The fluencemapbased CCCS algorithm includes two independent components: the calculation of the total energy released per mass (TERMA) and the convolution/superposition (CS) calculation of the energy deposition [15].
The TERMA calculation models the primary photons that interact with the material. Because the CCCS algorithm was used for EPIDbased 3D dosimetric verifications in this study, a measured source model was used for the TERMA calculation instead of a theoretically calculated source model, which must account for all elements of the system, including the flattening filter, the collimator system, the multileaf system, etc.[12]. The EPIDbased fluence maps, which included the intensity of the beam emitted from the head of the linear accelerator, were taken as the input. The spectrum of the beam was obtained from the phasespace files that were simulated for the Varian Trilogy 6MV and 10MV photon beams using BEAMnrc [28], and the geometry and material parameters of the linear accelerator were obtained from the manufacturer. The TERMA calculation was performed based on the electron density and the linear attenuation coefficient of water during the ray tracing, the latter of which was obtained from the National Institute of Standards and Technology (NIST). The calculation had a complexity of O (N^{3}) for N^{3} voxels. Each voxel was calculated by one GPU thread.
Three types of kernel expressions can be used in the discretization of the CCCS algorithm: the differential kernel (DK), the cumulative kernel (CK) and the cumulative–cumulative kernel (CCK) [29]. Lu et al. reported that the CCK is the most accurate because of its inherent voxel integration; therefore, the CCK was used in this study. The CCK was precalculated from discrete Monte Carlo kernels, which were simulated using the edknrc program [30] for uniform water. The lookup method based on tabulated kernels was used in the calculation. The kernels were scaled in accordance with the electron densities between the interaction point and the receiving points to allow the water kernels to be applied to arbitrary inhomogeneous media [15,29]. Considering that the kernel intensity had an angular distribution, the collapsedcone directions were divided into nonuniformly sampled zenith angles to improve the calculation speed. The collapsedcone directions were divided as follows: 22 uniformly sampled angles from 0° to 44°, 1 angle from 44° to 50°, 9 uniformly sampled angles from 50° to 95°, 1 angle from 95° to 120° and 3 uniformly sampled angles from 120° to 180°. The calculation also included 8 uniformly sampled azimuthal angles. In total, 280 collapsedcone directions were used. Each voxel was calculated in one kernel direction by one GPU thread and then looped through all directions. This method is well suited to GPU implementation because adjacent threads can be used to access adjacent memory locations [19].
Algorithm heterogeneity evaluation
A slab phantom similar to that used by Ahnesjo [15] was employed in this study. This phantom consisted of slabs of simulated tissues (adipose (A, density = 0.920 g cm^{−3}), muscle (M, density = 1.040 g cm^{−3}), bone (B, density = 1.850 g cm^{−3}) and lung (L, density = 0.250 g cm^{−3})). The thickness of the phantom was 30 cm, and its structure is illustrated in Figure 1. The sourcetosurface distance (SSD) was 85 cm. The sizes of the evaluation fields were 3 cm × 3 cm, 10 cm × 10 cm and 20 cm × 20 cm. The calculation methods included (1) the GPUbased CCCS algorithm, (2) a full Monte Carlo simulation and (3) the analytical anisotropic algorithm (AAA). The size of the calculation grid for all methods was 3 mm × 3 mm × 3 mm. The EGSnrc user code DOSRZnrc [31], with the phasespace files mentioned in the previous section, was used to simulate the dose delivered to the material. Because the treatment plans were created and optimized using Eclipse v10.0.28 (Varian Medical Systems, Palo Alto, CA, USA), which uses the AAA for dose calculations, the AAA was also considered in the evaluation.
3D dosimetric verifications for example clinical treatment plans
Treatment planning
One headandneck IMRT treatment plan and one lung IMRT treatment plan were chosen for 6 MV dosimetric verification. One lung IMRT treatment plan and one prostate IMRT treatment plan were chosen for 10 MV dosimetric verification. The treatment plans were created and optimized using Eclipse v10.0.28, and the AAA was used as the dosecalculation algorithm. The size of the calculation grid was 3 mm × 3 mm × 3 mm. Nine radiation fields were used for the headandneck and prostate treatment plans, and five radiation fields were used for the lung treatment plan. The optimization objective functions of the targets were generally set such that the planning target volume (PTV) received at least 95% of the prescribed dose. The optimization objective functions of the other normal tissues were set following clinical requirements. The prescribed doses for the 6 MV headandneck, 6 MV lung, 10 MV lung and 10 MV prostate plans were 68.2 Gy, 66.0 Gy, 66.0 Gy and 67.5 Gy, respectively.
Dosimetric verifications
The 3D dosereconstruction results that were obtained based on the EPID were compared to the TPS results. Because a measured source model was used in our dose calculation, it was difficult to directly quantify the differences in the dosecalculation algorithms and the sources of delivery error for these two different methods. It is well known that dosecalculation algorithms, including the CCCS algorithm and the AAA, are accurate in a uniform water phantom [15,29,32]. Therefore, the material within the patient’s body structure was first overridden to water (mass density 1.0 g cm^{3}). This aided in excluding any inherent difference between the dosecalculation algorithms when quantifying the effect of the delivery error. The TPS was used to recalculate the dose distributions in the overridden water phantom for the verification plans, and the EPIDbased fluence maps were used to reconstruct the dose distributions in the same phantom. The comparisons were performed with respect to the following: 1) Dose differences, including the differences in the mean doses (D_{mean}) and the maximum doses (D_{max}) to the targets and the organs at risk (OARs). The percent dose differences were calculated using the following equations:
2) 3D gamma evaluations. The same areas that were used for the dose comparisons were treated as the regions of interest (ROIs). The dose grid size was resampled to 1 mm × 1 mm × 1 mm before the gamma evaluation was performed. We set 3% of the maximum dose as the dose criterion and 3 mm as the distance criterion (henceforth referred to as the 3%/3 mm criterion), and a tighter 2%/2 mm criterion was also evaluated. A statistical evaluation was performed to determine the point at which the dose was greater than 10% of the maximum dose.
The 3D dosereconstruction results and the TPS results from the original patient CTs were also considered for comparison. These comparisons were performed with respect to the following: 1) Dose differences, including those in the D_{mean} and D_{max} values of the targets and OARs, and 2) 3D gamma evaluations. The evaluation criteria were the same as those listed above.
One disadvantage of the 3D gamma evaluation is its long computation time [33,34]. With the dose grid size resampled to 1 mm × 1 mm × 1 mm, the computation time using a fast algorithm for the body gamma evaluation would be approximately 1 minute [34]. However, applying a GPU implementation of the gamma evaluation can reduce the computation time [35]. Because GPU texture memory provides highefficiency trilinear interpolation for 3D data, the dose data were bound to the GPU texture memory for the resampling process. The gamma index of each voxel was calculated by a separate GPU thread. Voxels with dose values lower than the threshold were excluded from the calculation process.
All calculation algorithms discussed above were implemented in Visual C++. All computations were performed on a notebook with an Intel Core i74700MQ 2.40 GHz CPU and an NVIDIA GTX765 M video card.
Results
Evaluations of EPIDbased fluence maps
Figure 2 shows the centralpoint relative dose values in air that were measured using the ion chamber and the EPID for square fields of various sizes. When our deconvolution model correction was applied to the raw EPID data, the maximum difference between the EPID results and the ionchamber measurements was less than 0.5% for both the 6 MV and 10MV photon beams. Figure 3 shows the fluence profiles on the central axis that were measured using the ion chamber and the EPID after correction. The agreement was good for both the crossplane direction and the inplane direction.
Evaluations of dose algorithm heterogeneity
Figure 4 presents the axialdepth–dose curves in the heterogeneous slab phantom that were obtained using the various algorithms for the 6 MV photon beam. The percent dose differences presented in this section were normalized to the Monte Carlo simulated doses for a 2 cm depth in the case of the 6 MV beam and for a 5 cm depth in the case of the 10 MV beam; these depths correspond to positions in the muscle slab near the depth at which the maximum dose occurred. The CCCS results agreed well with the results of the Monte Carlo simulation. The dose difference was less than 1.5%, except at the interfaces of different materials. For the AAA, the dose results for the lowdensity lung region were lower than the Monte Carlo results. The largest dose difference was 7.5%. The dose results for the highdensity bone region were higher than the Monte Carlo results, and the largest dose difference was 15%. Figure 5 presents the axialdepth–dose curves in the heterogeneous slab phantom that were obtained using the various algorithms for the 10 MV photon beam. In this case, the CCCS results again agreed well with the results of the Monte Carlo simulation except at the interfaces of different materials. For the AAA, the results were superior to those obtained in the 6MV case. The largest dose difference in the lowdensity lung region was 6.4%, and that in the highdensity bone region was 9.0%.
3D dosimetric verifications
Table 1 presents the gamma pass rates of the gross target volume (GTV), the clinical target volume (CTV), the planning target volume (PTV) and the entire body for the treatment plans. The results for the patient body overridden to water and for the original material are shown. The gamma statistics are within the 10% isodose contours. Good results were obtained using the 3%/3 mm criterion for the planned and EPIDbased reconstructed dose distributions in both cases. However, for the tight criterion (2%/2 mm), the pass rates of the original material were decreased compared with the overriddentowater results; the maximum difference was 22%.
The detailed statistics of the dose differences for the targets and the OARs of the patient body overridden to water are summarized in Table 2, and the results for the original material are summarized in Table 3. Acceptable agreement was observed between the EPIDbased results and the planned results for both the target regions and the OARs. However, the dose differences in the case of the original material were slightly larger than those in the overriddentowater case. In the case of the patient body overridden to water, the average differences in the mean doses and the maximum doses to the target regions were 0.11 ± 0.52% and 0.32 ± 1.17%, respectively, and the corresponding average values for the OARs were 0.34 ± 1.35% and 0.12 ± 1.74%, respectively. In the case of the original material, the average differences in the mean doses and the maximum point doses to the target regions were 0.57 ± 0.71% and 1.34 ± 0.66%, respectively, and the corresponding average values for the OARs were 0.80 ± 2.05% and 0.41 ± 2.32%, respectively.
Computation times for dosimetric verifications
Table 4 lists the computation times associated with the various components of the calculation for the GPU implementation, including the initialization of the CUDA subsystem, the fluencemap calculation, the CCCS calculation of one field and the gamma evaluation of the body. According to these results, the total dose computation time for a single field was less than 8 s, and the gamma evaluation was completed in approximately 1 s.
Discussion
Treatment verification using a 3D dose reconstruction based on information acquired in the treatment room is feasible. EPIDbased reconstruction can be used to determine the delivery errors of a treatment plan and to estimate the errors of a TPS dosecalculation algorithm. During treatment, the output of the linear accelerator may differ from the planned delivery determined using the TPS. For example, the positioning accuracy of the multileaf collimator may lead to differences in the delivered and planned dose distributions. Because of its high resolution and linear dose response, an EPID can effectively detect the delivery errors of a linear accelerator [36,37]. Accurate and reliable dosecalculation algorithms that consider density inhomogeneities are necessary, especially for the head, neck and thorax [38,39]. For this purpose, the calculation accuracy of the CCCS algorithm for inhomogeneous tissues is acceptable [16,29]. An independent dosecalculation algorithm based on the CCCS algorithm was implemented to successfully estimate the errors of a TPS dosecalculation algorithm. In this study, fluence maps were obtained based on EPID measurements. These fluence maps were used to reconstruct the doses delivered to overriddentowater phantoms and the original patients. According to the gamma analysis results, good results were obtained for both cases using the 3%/3 mm criterion. The average gamma pass rates were greater than 94%. However, when a tighter criterion (2%/2 mm) was used, the pass rates of the original patients were decreased compared with the pass rates of the overriddentowater phantoms. The maximum decreases for the headandneck, lung and prostate treatment plans were approximately 15%, 20% and 8%, respectively. For both 6 MV treatment plans (headandneck and lung), the EPIDbased doses to the spinal cord were found to be lower than the planned doses; based on the heterogeneity evaluation, this result can be attributed to the presence of highdensity bone material around the spinal cord, which caused the dose results calculated using the AAA to be higher than those calculated using the CCCS algorithm. Fogliata et al. have reported similar results [40]. These results indicate that compared with the uniformphantombased dosimetric verification method, the complete 3D dosimetric verification method using an independent dosecalculation algorithm can provide more detailed evaluation information.
The advantage of the GPU implementation of the CCCS algorithm is that the 3D dose computation time is significantly reduced, while the necessary accuracy of the dose calculation in inhomogeneous tissues is maintained. These advantages are important for 3D dosimetric techniques that are intended to be applied in laborintensive clinical verifications. For either pretreatment or in vivo dosimetric verifications, GPUbased dose reconstruction can achieve true realtime calculations. Chen et al. have reported that the singlefield computation time can reach 0.5 s when two higherperformance GPUs and a fully optimized CCCS algorithm are used [19]. In this study, for the three different sites considered (which contained 160 × 100 × 72, 180 × 50 × 76 and 136 × 120 × 106 voxels), the singlefield dosereconstruction times were 7 s, 6 s and 8 s, respectively. Thus, our study indicates that there is room for further improvement: 1) Highperformance hardware can effectively improve the computation speed. An NVIDIA GTX765 M notebook video card with 768 CUDA cores was used to perform the computations reported in this study. When we performed the computations on an NVIDIA GTXTITAN, which contains 2688 CUDA cores, the computation speed was 56 times faster than that achieved using the GTX765 M. 2) Optimization of the CCCS algorithm can also improve the computation speed. Here, we used the lookup method based on tabulated kernels for the calculation. The computation time could be reduced by using the exponential kernel model and the builtin NVIDIA special function unit.
Not only the dose calculation but also the fluencemap reconstruction and the gammaindex calculation were implemented using the GPU. The computation time required to obtain a fluence map from EPID images of 1024 × 768 pixels was reduced to 0.3 s by computing FFTs in parallel on the NVIDIA GPU. By contrast, the computation time that was required when a single CPU was used was approximately 8 s. Wendling et al. have presented a fast algorithm for gamma calculation, for which the computation time for a 38.5 cm × 35.3 cm × 25.3 cm dose dataset at a 0.1cm grid resolution was found to be 28 s [34]. This algorithm still represents a nonnegligible time cost for realtime clinical verifications. With the gamma calculation implemented using a GPU, the computation time, including data resampling and gamma calculation, for a 48.0 cm × 30.0 cm × 21.6 cm dose dataset at a 0.1cm grid resolution was found to be less than 1.8 s.
Conclusion
A 3D pretreatment dosimetric verification method based on an EPID was presented. A GPU calculation engine can boost the speed of 3D dose and gamma evaluation calculations, thereby bringing the calculation and evaluation times to the level required for realtime verifications. Although this investigation focused on pretreatment dosimetric verification, it is our intent to extend the methodology to true realtime 3D dosimetric verification using EPID transit measurements and a GPU calculation engine.
Consent
Written informed consent was obtained from the patient for the publication of this report and any accompanying images.
Abbreviations
 GPU:

Graphic processing unit
 EPID:

Electronic portal imaging device
 IMRT:

Intensitymodulated radiation therapy
 CCCS:

Collapsedcone convolution/superposition
 TERMA:

Total energy released per mass
 AAA:

Analytical anisotropic algorithm
 OAR:

Organ at risk
 DVH:

Dosevolume histogram
References
 1.
van Zijtveld M, Dirkx ML, de Boer HC, Heijmen BJ. 3D dose reconstruction for clinical evaluation of IMRT pretreatment verification with an EPID. Radiother Oncol. 2007;82:201–7.
 2.
van Elmpt W, Petit S, De Ruysscher D, Lambin P, Dekker A. 3D dose delivery verification using repeated conebeam imaging and EPID dosimetry for stereotactic body radiotherapy of nonsmall cell lung cancer. J European Soc Ther Radiology Oncol. 2010;94:188–94.
 3.
Nelms BE, Opp D, Robinson J, Wolf TK, Zhang G, Moros E, et al. VMAT QA: measurementguided 4D dose reconstruction on a patient. Med Phys. 2012;39:4228–38.
 4.
Nelms BE, Zhen H, Tome WA. Perbeam, planar IMRT QA passing rates do not predict clinically relevant patient dose errors. Med Phys. 2011;38:1037–44.
 5.
PecharrománGallego R, Mans A, Sonke JJ, Stroom JC, OlacireguiRuiz I, van Herk M, et al. Simplifying EPID dosimetry for IMRT treatment verification. Med Phys. 2011;38:983.
 6.
Mans A, Remeijer P, OlacireguiRuiz I, Wendling M, Sonke JJ, Mijnheer B, et al. 3D Dosimetric verification of volumetricmodulated arc therapy by portal dosimetry. Radiother Oncol. 2010;94:181–7.
 7.
van Elmpt W, McDermott L, Nijsten S, Wendling M, Lambin P, Mijnheer B. A literature review of electronic portal imaging for radiotherapy dosimetry. Radiother Oncol. 2008;88:289–309.
 8.
Steciw S, Warkentin B, Rathee S, Fallone BG. Threedimensional IMRT verification with a flatpanel EPID. Med Phys. 2005;32:600–12.
 9.
Wendling M, McDermott LN, Mans A, Sonke JJ, van Herk M, Mijnheer BJ. A simple backprojection algorithm for 3D in vivo EPID dosimetry of IMRT treatments. Med Phys. 2009;36:3310–21.
 10.
Wendling M, NM L, Mans A, OlacireguiRuiz I, PecharromanGallego R, Sonke JJ, et al. Med Phys. 2012;39:367.
 11.
Renner WD, Norton K, Holmes T. A method for deconvolution of integrated electronic portal images to obtain incident fluence for dose reconstruction. J Appl Clin Med Phys. 2005;6:22–39.
 12.
Renner WD. 3D dose reconstruction to insure correct external beam treatment of patients. Med Dosim. 2007;32:157–65.
 13.
Liang Y, Kim GY, Pawlicki T, Mundt AJ, Mell LK. Feasibility study on dosimetry verification of volumetricmodulated arc therapybased total marrow irradiation. J Appl Clin Med Phys. 2013;14:3852.
 14.
Knoos T, Ahnesjo A, Nilsson P, Weber L. Limitations of a pencil beam approach to photon dose calculations in lung tissue. Phys Med Biol. 1995;40:1411–20.
 15.
Ahnesjo A. Collapsed cone convolution of radiant energy for photon dose calculation in heterogeneous media. Med Phys. 1989;16:577–92.
 16.
Sterpin E, Salvat F, Olivera G, Vynckier S. Monte Carlo evaluation of the convolution/superposition algorithm of HiArt tomotherapy in heterogeneous phantoms and clinical cases. Med Phys. 2009;36:1566–75.
 17.
Hissoiny S, Ozell B, Despres P. Fast convolutionsuperposition dose calculation on graphics hardware. Med Phys. 2009;36:1998–2005.
 18.
Hissoiny S, Ozell B, Despres P. A convolutionsuperposition dose calculation engine for GPUs. Med Phys. 2010;37:1029–37.
 19.
Chen Q, Chen M, Lu W. Ultrafast convolution/superposition using tabulated and exponential kernels on GPU. Med Phys. 2011;38:1150–61.
 20.
Jacques R, Wong J, Taylor R, McNutt T. Realtime dose computation: GPUaccelerated source modeling and superposition/convolution. Med Phys. 2011;38:294–305.
 21.
Ko L, Kim JO, Siebers JV. Investigation of the optimal backscatter for an aSi electronic portal imaging device. Phys Med Biol. 2004;49:1723–38.
 22.
Moore JA, Siebers JV. Verification of the optimal backscatter for an aSi electronic portal imaging device. Phys Med Biol. 2005;50:2341–50.
 23.
Berry SL, Polvorosa CS, Wuu CS. A field size specific backscatter correction algorithm for accurate EPID dosimetry. Med Phys. 2010;37:2425–34.
 24.
Rowshanfarzad P, McCurdy BM, Sabet M, Lee C, O’Connor DJ, Greer PB. Measurement and modeling of the effect of support arm backscatter on dosimetry with a varian EPID. Med Phys. 2010;37:2269–78.
 25.
Jin GH, Zhu JH, Chen LX, Deng XW, Huang BT. Yuan K. Gantry angledependent correction of dose detection error due to panel position displacement in IMRT dose verification using EPIDs. Phys Med: Liu XW; 2013.
 26.
Jinhan CLx ZHU, Guanghua JIN, Wufei CAO, Xiaowu DENG, Xiaowei LIU. The study of twodimensional dosimetric verifications of volumetric intensitymodulated arc therapy based on electronic portal imaging devices. Chinese J Radiation Oncol. 2012;21:550–3.
 27.
Kiefer J. Sequential Minimax Search for a Maximum. Proceed The American Math Soc. 1953;4:502–2.
 28.
Rogers DW, Faddegon BA, Ding GX, Ma CM, We J, Mackie TR. BEAM: a Monte Carlo code to simulate radiotherapy treatment units. Med Phys. 1995;22:503–24.
 29.
Lu W, Olivera GH, Chen ML, Reckwerdt PJ, Mackie TR. Accurate convolution/superposition for multiresolution dose calculation using cumulative tabulated kernels. Phys Med Biol. 2005;50:655–80.
 30.
MainegraHing E, Rogers DW, Kawrakow I. Calculation of photon energy deposition kernels and electron dose point kernels in water. Med Phys. 2005;32:685–99.
 31.
Walters IK B, Rogers DWO. DOSXYZnrc users manual. Ionizing radiation standards National Research Council of Canada. 2011.
 32.
Tillikainen L, Helminen H, Torsti T, Siljamaki S, Alakuijala J, Pyyry J, et al. A 3D pencilbeambased superposition algorithm for photon dose calculation in heterogeneous media. Phys Med Biol. 2008;53:3821–39.
 33.
Low DA, Dempsey JF. Evaluation of the gamma dose distribution comparison method. Med Phys. 2003;30:2455–64.
 34.
Wendling M, Zijp LJ, McDermott LN, Smit EJ, Sonke JJ, Mijnheer BJ, et al. A fast algorithm for gamma evaluation in 3D. Med Phys. 2007;34:1647–54.
 35.
Gu X, Jia X, Jiang SB. GPUbased fast gamma index calculation. Phys Med Biol. 2011;56:1431–41.
 36.
Wendling M, Louwe RJ, McDermott LN, Sonke JJ, van Herk M, Mijnheer BJ. Accurate twodimensional IMRT verification using a backprojection EPID dosimetry method. Med Phys. 2006;33:259–73.
 37.
Zhu J, Chen L, Jin G, Luo G, Cao W, Liu X. A comparison of VMAT dosimetric verifications between fixed and rotating gantry positions. Phys Med Biol. 2013;58:1315–22.
 38.
Nikos Papanikolaou JJB, Boyer AL, Constantin K, Eric K, Rock Mackie T, Michael S, et al. Tissue inhomogeneity corrections for megavoltage photon beams. Madison,WI: Task Group No. 65 of the Radiation Therapy Committee of the American Association of Physicists in Medicine. 2004.
 39.
Tommy K, Elinore W, Luca C, Carsten B, Antonella F, Dirk A, et al. Comparison of dose calculation algorithms for treatment planning in external photon beam therapy for clinical situations. Phys Med Biol. 2006;51:5785.
 40.
Fogliata A, Nicolini G, Clivio A, Vanetti E, Cozzi L. Dosimetric evaluation of Acuros XB Advanced dose calculation algorithm in heterogeneous media. Radiat Oncol. 2011;6:82.
Acknowledgements
This study was supported by the National Natural Science Foundation of China (8127248) and a Combination of Guangdong Province and Ministry of Education research project (Grant ID 2012B091000144).
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors contributed to the research, data collection and analysis. JZ and LC contributed equally to this work. JZ developed the GPU calculation program and performed the data analysis. LC and XL conceived of the study, including the study design and coordination. AC and GL assisted with the data measurements. XD assisted in the study design and the data analysis. All authors read and approved the final manuscript.
Jinhan Zhu and Lixin Chen contributed equally to this work.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhu, J., Chen, L., Chen, A. et al. Fast 3D dosimetric verifications based on an electronic portal imaging device using a GPU calculation engine. Radiat Oncol 10, 85 (2015). https://doi.org/10.1186/s1301401503877
Received:
Accepted:
Published:
Keywords
 3D dosimetric verification
 EPID
 GPU