- Open Access
Potential dosimetric benefits of adaptive tumor tracking over the internal target volume concept for stereotactic body radiation therapy of pancreatic cancer
Radiation Oncologyvolume 12, Article number: 175 (2017)
Radiotherapy for pancreatic cancer has two major challenges: (I) the tumor is adjacent to several critical organs and, (II) the mobility of both, the tumor and its surrounding organs at risk (OARs). A treatment planning study simulating stereotactic body radiation therapy (SBRT) for pancreatic tumors with both the internal target volume (ITV) concept and the tumor tracking approach was performed. The two respiratory motion-management techniques were compared in terms of doses to the target volume and organs at risk.
Methods and Materials
Two volumetric-modulated arc therapy (VMAT) treatment plans (5 × 5 Gy) were created for each of the 12 previously treated pancreatic cancer patients, one using the ITV concept and one the tumor tracking approach. To better evaluate the overall dose delivered to the moving tumor volume, 4D dose calculations were performed on four-dimensional computed tomography (4DCT) scans. The resulting planning target volume (PTV) size for each technique was analyzed. Target and OAR dose parameters were reported and analyzed for both 3D and 4D dose calculation.
Tumor motion ranged from 1.3 to 11.2 mm. Tracking led to a reduction of PTV size (max. 39.2%) accompanied with significant better tumor coverage (p<0.05, paired Wilcoxon signed rank test) both in 3D and 4D dose calculations and improved organ at risk sparing. Especially for duodenum, stomach and liver, the mean dose was significantly reduced (p<0.05) with tracking for 3D and 4D dose calculations.
By using an adaptive tumor tracking approach for respiratory-induced pancreatic motion management, a significant reduction in PTV size can be achieved, which subsequently facilitates treatment planning, and improves organ dose sparing. The dosimetric benefit of tumor tracking is organ and patient-specific.
Pancreatic cancer is one of the most aggressive tumors with high mortality rates. According to cancer mortality statistics, for 2015 in European Union (EU) a total of 85300 pancreatic cancer deaths are predicted (42700 men and 42600 women), making it the fourth leading cause of cancer-related deaths for both genders in the EU .
Surgery is generally considered the standard treatment for patients with resectable pancreatic cancer (PC) while standard-dose chemoradiation therapy and chemotherapy alone are suitable in case of unresectable disease . Stereotactic body radiation therapy (SBRT) allows for higher doses per fraction to the tumor than the conventional radiotherapy with sharp dose gradients in a shorter overall treatment time, and better dose sparing of adjacent organs at risk. It has been successful in the treatment of thoracic tumors and has been used for intra-abdominal sites, including pancreas [2–5]. Besides the patient comfort due to significantly shorter treatment times, additional advantages of SBRT include the possibility of being better integrated into multidisciplinary treatment concepts, significantly less bone marrow toxicity and shorter interruption of full-dose chemotherapy . SBRT for pancreatic tumors has been either delivered by traditional linear accelerator [6, 7] or a robotic arm-mounted linear accelerator [6, 8].
Radiotherapy for pancreatic cancer has two major challenges: (i) the tumor is proximal to several critical organs such as duodenum, stomach, small intestine, kidneys, liver and spinal cord, and (ii) the mobility of both the tumor and its surrounding organs at risk (OARs). The duodenum, which often lies immediately adjacent to pancreatic tumors, is most commonly the organ that limits the maximum dose that can be delivered to the target .
As radiotherapy treatments become more and more conformal, without a full understanding of the individualized tumor/organ motion, the target volume could be under-dosed while the nearby OARs could be over-dosed, resulting in the degradation of the overall treatment effectiveness and accuracy. Several studies have addressed the important aspect of pancreatic mobility caused by respiration using different imaging modalities such as magnetic resonance imaging (MRI) [9, 10], four-dimensional computed tomography (4DCT) [11–14], cone-beam computed tomography , fluoroscopy [16, 17] and ultrasound . They have shown that, on average, the magnitude of pancreatic motion is largest in the cranial-caudal (superior-inferior) direction and the largest reported value among them was 24±16 mm .
The implementation of multi-fraction SBRT has demonstrated promising results with favorable tumor control and toxicity rates  which might be further improved using advanced motion management strategies which lead to margin reduction and therefore PTV size reduction.
Different strategies have been proposed to account for respiratory motion [20, 21]: active breath-hold techniques , abdominal compression techniques , gating during end-expiration , motion-encompassing methods (Internal Target Volume (ITV) concept) , and respiration-synchronized techniques (tumor tracking approach) . Tumor tracking is an advanced motion management strategy, which adjusts dynamically the treatment beam [26–28] or the patient position [29, 30] to the tumor motion.
Currently, most clinics use the ITV concept, which accounts for the whole tumor motion resulting in broad safety margins. However, among the aforementioned motion management methods, the tumor tracking is recognized as an advanced method of managing respiratory motion reducing the size of the planning target volume (PTV). This might lead to improved targeting and better tumor control while minimizes the radiation-induced toxicity to normal tissue. Real-time tracking of the pancreatic target has been mainly performed by a robotic arm-mounted linear accelerator  which utilizes prediction and correlation models [32–34] in order to relate the external breathing signals to 3D tumor motion (Cyberknife Synchrony, Saunnyvale, CA).
The aim of this study was to investigate the possible dosimetric advantages of using a tumor tracking approach compared with the ITV concept to account for respiratory-induced pancreatic tumor motion during gantry-based SBRT. For a better estimation of the dose in the moving tumor, we implemented 4D dose calculation of SBRT treatments with the ITV [25, 35] and the tumor tracking approach [25, 36, 37]. This allows to include respiratory motion effects into the evaluation of both motion managements strategies, which has not been done before so far for pancreatic tumors. Additionally, the differences between the 3D and 4D dose calculations were evaluated for both motion management strategies.
Methods and materials
Patients and four-dimensional computed tomography data
Four-dimensional computed tomography (4DCT) data, from 12 consecutive patients with pancreatic cancer, was included in this retrospective study. Among them, three had undergone total duodenectomy. Each patient had received a 4DCT scan using the SOMATOM Definition AS Open (Siemens AG, Germany) CT scanner. Scans were acquired under free breathing without coaching and with respiratory monitoring (RPM, Varian Medical Systems, Palo Alto CA). An average CT and ten respiratory phase CTs were reconstructed. For all but two patients amplitude-sorted 4DCT scans had been obtained, while for two patients phase-sorted 4DCT scans. The data sets were transferred to Eclipse Treatment Planning System (Varian Medical Systems, Palo Alto). Slice thickness of the 4DCT scans was 2 mm and a pitch of 0.09 was used.
Organ and target segmentation
The gross tumor volume (GTV) and multiple abdominal organs were contoured manually by the radiation oncologists of the Radiation Oncology Department of University Hospital Zurich on one respiratory phase of the 4DCT image set. This contouring phase, which is approximately between inhale and exhale, was considered to better approximate the average tumor position during the full breathing cycle . The following abdominal organs were considered as organs at risk (OAR): duodenum, stomach, bowel, left and right kidneys, liver and spinal cord.
Two individual structure sets were generated to simulate treatment planning for dynamic motion compensation with tumor tracking and for the conventional motion encompassing ITV concept. For tracking, the structures on the contouring phase were used with the GTV as primary target volume (GTV track ). For the ITV concept, structure propagation was performed to all respiration phases of the 4DCT image set by deformable registration using MIM Maestro (v6.1, MIM Software Inc., Cleveland, OH). The individual structures from all phases were combined to a total structure, which was then used for treatment planning. The motion envelope enclosing all GTVs was defined as ITV. The OARs were also created from their motion envelope over all respiratory 4DCT phases. An ITV-to-PTV margin of 5 mm was added for the PTV itv . For tracking, a fixed 5mm-margin was added to GTV track to form the PTV track .
The center-of-volume motion for GTV at each respiratory phase along the superior-inferior (SI), anterior-posterior (AP), and left-right axes was measured using MIM software. The tumor motion per direction was determined as the maximal displacement. Three dimensional tumor motion vector was defined as: 3D tumor motion= (LR 2+AP 2+SI 2)1/2.
Treatment planning - three dimensional (3D) dose calculation
For both concepts, ITV and tracking, stereotactic body radiotherapy (SBRT) was planned based on the concept-specific structure sets. SBRT was planned using a nominal energy of 10 MV photon beam in flattening-filter-free mode. Two full arcs with collimator rotations of 5° and 355° were used. Treatment plan optimizations for volumetric-modulated arc therapy (VMAT) and dose calculations using the Analytical Anisotropic Algorithm (AAA) 13.6.23 were performed in Eclipse Treatment Planning System. The treatment plans were calculated on the average reconstruction of the 4DCT data set for the ITV concept and on the contouring phase for the tumor tracking concept. The derived dose distribution is referred from now on as 3D dose distribution.
A dose of 25 Gy in 5 fractions was prescribed to the 60%-isodose surrounding the PTV, which is now widely used [6, 8], corresponding to a maximum allowed dose of up to 41.5 Gy (166% of the prescribed dose). Treatment plans were created with more than 95% of the PTV receiving the prescribed dose. A risk adapted prescription similar to the Simultaneous Integrated Protection (SIP) concept of Gkika et al. was used . Parts of the PTV overlapping with OARs received between a minimum of 25 Gy and a maximum of 27.5 Gy. However, in the parts of the PTV not overlapping with OAR, we were aiming for a minimum dose of 27.5 Gy escalating the dose in the GTV as high as 41.5 Gy. The constraints for OARs were D max <25 Gy on the spinal cord, D max <27.5 Gy for the bowel, the duodenum, the stomach and D mean <10 Gy for the kidneys.
Four dimensional (4D) dose calculation
Inter-fractional and intra-fractional anatomical variations can introduce significant errors in dose delivered by radiation therapy. The impact of respiratory motion should not be ignored as respiration can induce both rigid body translation/rotation and organ deformation. The 3D dose distribution does not fully represent the dose accumulated by a moving tumor . 4D dose calculation considers the changing beam aperture and the changing anatomy. It strongly depends on the patient’s respiratory cycle and therefore should be performed explicitly on a patient-specific base. In this study, for the 4D dose calculation, we used the concept presented by Ehrbar et al. .
The respiratory cycles of all the patients were divided into 10 breathing phases, according to the 4DCT. The original treatment plans (for ITV and tumor tracking concept) were divided into angular segments, which were temporally assigned to the breathing phases. This procedure was done by using in-house MATLAB scripts [25, 39].The created sub-plans were imported back to the treatment planning system and calculated on the different breathing phase. For the tumor tracking concept, the sub-plans were calculated on the different phase CT image set while the beam isocenter was shifted according to the current center of volume position for the GTV to simulate adaptive motion compensation. Subsequently, 4D dose accumulation was performed in MIM Software.
The dose accumulation was performed against the initial contouring phase and the other CT phases were registrated to this one phase using deformable image registration. Following, the dose distributions at other phases were summed up to the dose distribution at reference CT to make 4D dose accumulation . The resulting accumulated 4D dose distributions were recorded for the GTV volume and the OARs of the reference phase.
PTV size and dosimetric comparisons
The two motion-management techniques were compared regarding the resulting PTV size, dose coverage of the tumor, and sparing of organs at risk. The dosimetric comparisons were performed in 3D and 4D. For the 3D dose comparison, the dose parameters were evaluated on the structures used for treatment planning. Dose parameters for the target (GTV track and ITV, respectively) include the minimum (D min ), maximum (D max ), mean dose (D mean ), D 2, D 95 and D 98 - the dose received by 2, 95 and 98% of the target volume. For OARs, two dose parameters were recorded: the D mean and the D 0.1cc , which corresponds to the dose received by 0.1 cc of the volume. For the comparison of 4D dose distributions, the same dose parameters were evaluated within the structures delineated on the reference phase for the dose accumulation. This means that for both motion-management techniques, the 4D target dose parameters were evaluated on the GTV.
Dosimetric comparisons between ITV concept and tumor tracking concept were performed with a two-sided Wilcoxon paired signed rank test, with p-values < 0.05 considered significant. The degree of association between the differences in dose parameters and the 3D tumor motion was calculated by Spearman’s rank correlation and a significance level of 5% was used (p-value < 0.05). Statistical analysis was performed using R v.3.2.4 .
Tumor motion and PTV size
The respiratory-induced pancreatic motion was studied in 12 patients (Table 1). For all patients, the largest amplitude of pancreatic motion was found in the SI direction. The 3D motion ranged from 1.3 to 11.2 mm. With an adaptive tumor tracking approach, the PTV size could be significantly (p<0.01) reduced with an average reduction of 16.8% (Table 1). Figure 1 shows the linear regression between PTV reduction and 3D motion. The values are moderately correlated (rho=−0.46, p=0.131).
Treatment plan quality
The reduction of PTV size allowed for a decrease of the overlap between the PTV and duodenum. This improved not only OAR dose sparing but also facilitated treatment planning optimization. For five patients (3, 9, 10, 11 and 12) not all the optimization objectives could be satisfied in the ITV treatment plan. For these five cases we had to comprise on the dose to the overlapping part of ITV and OARs (ITV-OAR), where we aimed for a dose of 37.5 Gy to 95% of the overlapping volume. Instead, this volume was compromised to values between 83.5% and 93.9% for these five patients. However, the optimization objectives were fulfilled for all patients using the tumor tracking concept.
Comparisons of 3D and 4D dose calculations
In Table 2, the two motion management techniques are compared using dosimetric parameters of the 3D and 4D dose calculations. The target dose parameters are additionally shown in Fig. 2. There is a clear dosimetric advantage for tumor tracking, regarding the tumor coverage (D 95% ) in both, 3D and 4D dose calculation. Additionally, many of the OAR dose parameters such as duodenum, stomach, and liver D mean showed a significant improvement with tracking in 3D and 4D dose. There were also cases in which significant OAR dose reduction was only observed in the 4D dose distribution (bowel, left and right kidney D mean ), or the 3D dose distribution (bowel D 0.1cc ). Furthermore, highly correlated (rho ≥ 0.71), in absolute values, and significant (p<0.05) relationships were observed between the 3D dose difference (3D Tracking −3D ITV ) and the 3D tumor motion for the bowel D mean (rho=−0.73) and the target volume Dmean (rho = 0.73). No significant results were found for the correlation coefficients investigating the 4D dose differences.
Figure 3 shows dose-volume histograms (DVHs) for 3D and 4D dose calculations. The data is presented for two typical patients, patient 1 with pancreatic motion 11.2 mm and patient 7 with pancreatic motion 2.1 mm. For patient 1 there is an obvious gain in the dose coverage of the tumor and OAR dose sparing with tumor tracking approach. However, for patient 7 minimal improvements were observed regarding the tumor coverage and the OAR dose reduction. Additionally, using the ITV concept for patient 1 the 3D dose results in an overestimation for all OARs in comparison with 4D dose calculation. On the other hand, for the same patient using the tumor tracking approach almost no difference between the two dose calculations is noticed with the exception for the dose to the duodenum.
In the present study we investigated the potential dosimetric benefits of using the tumor tracking concept instead of the ITV concept for the respiratory-induced motion of pancreatic tumors. Tumor motion ranged between 2 and 11 mm. Planning target volumes were reduced using the tracking approach with increasing tumor motion. For tumor tracking all 3D plans fulfilled the constraints whereas for the ITV concept we had to compromise the PTV coverage in five cases. Mean and maximum doses to most organs at risk were significantly reduced using tumor tracking. However, using 4D calculation this significant reduction could only be confirmed for the mean doses to the organs at risk.
Our motion results have shown that the magnitude of motion is largest in the SI direction and are in agreement with previous studies with 4DCT. Tai et al. , Lens et al. , and Hallman et al.  found a mean tumor motion between 5 and 8 mm in SI direction. However, in the literature larger pancreas motion has been also reported. Particularly, Minn et al.  calculated on a planning 4DCT the range of centroid movement in SI direction between 0.9 to 28.8 m (from the maximum inspiration to the maximum expiration). Similar large motion amplitudes were reported by Bussels et al.  and Feng et al.  using MRI imaging. Observed differences in the pancreatic tumor motion between the current work and the previous studies may be due to patient heterogeneity and variation of inter-fractional pancreatic motion .
When the ITV concept is used, a large volume is unnecessarily irradiated since the whole envelope of motion is treated. This can be avoided if the tumor tracking is considered. Our findings with the tumor tracking approach led to a mean PTV volume reduction of 16.8%±8.4% (range: −4.7 to −39.2%). This significant reduction of the PTV is in agreement with the result found by Lens et al. using the midventilation approach as a motion management technique for the pancreatic tumors. They found a mean absolute PTV volume reduction of 13.9%. However, it is worth mentioning that this single-phase tracking is able to deposit higher dose to target volume than the midventilation approach.
Tumor tracking allowed fulfilling all planning objectives whereas for five ITV plans we had to make a compromise on the PTV coverage in order to spare the organs at risk adequately. Therefore, we consider tracking as a promising tool for complicated cases for which the objectives cannot be fulfilled using the ITV concept. Similarly, for tracking of lung tumors Ehrbar et al.  showed a significant reduction in the mean dose to the lung compared to the ITV concept.
The 4D dose calculation can be considered as an appropriate tool for better understanding of the differences between a static 3D dose distribution and taking the respiratory motion into account. 4D dose calculation has been also performed by other authors [25, 35–37, 45] but mostly for lung [25, 37, 45] and not for pancreas. The results of our study have shown that the 3D dose distribution adequately represented the actual 4D dose for the mean tumor dose as well as the coverage of the tumor for both concepts. This is different compared to investigations for the lung, where the 3D dose calculation underestimates the dose to the tumor compared to the actual 4D calculation for the ITV concept [25, 45]. This difference is due to the changing density in the lung, which is not that pronounced in the pancreatic region.
In our study a single 4DCT was acquired per patient and based on this, we performed our treatment planning. However, it has been shown that the respiratory-induced tumor motion in pancreatic cancer patients given by a single 4DCT is not necessarily representative of the amplitude during treatment . Furthermore, during the treatment of pancreatic tumors, non-rigid shape and position variations occur frequently both in the tumor and the OARs . These effects were not considered in this work as structure delineation and 4D dose calculations were based on a single 4DCT before the treatment. Respiratory patterns reproducibility during the imaging and delivery is an important issue that needs further investigation . Such changing breathing and motion patterns are not an issue for adaptive tracking approaches which use real-time position feedback of the tumor, since they are able to adapt to the actual tumor position during delivery and this improves the daily both tumor coverage and OAR sparing.
Another limitation of the study was that we assumed perfect tumor tracking. We wanted to show what is possible with a perfect system rather than evaluating possible current systems. We assumed that generating the PTV with a 5 mm safety margin on the ITV/GTV, may be adequate to account for tracking errors, residual organ deformation and rotation . However, we performed 4D dose calculation which deals with the main contributors of tracking inaccuracies such as rotation and deformation .
Our results indicate that there is a possible benefit of using tumor tracking approach instead of ITV concept in SBRT treatment for pancreatic cancer patients. However, further investigation is needed in order to relate the present results to clinical outcome in terms of toxicity. This study showed that the tumor motion differs between patients, indicating that tumor motion should be assessed individually. Furthermore, the tumor tracking approach resulted in better tumor coverage and a significant reduction of PTV facilitating subsequently treatment planning and decreasing mean doses of OARs.
Taking the pancreatic tumor motion into account is important and critical in SBRT. If this movement is neglected, it could offset the potential benefits of SBRT. In this study two motion management techniques, ITV concept and tumor tracking concept, were evaluated and compared to each other using 3D and 4D dose calculations.
The dosimetric comparisons showed that the dose benefit is patient and OAR specific. Tumor tracking might be considered for patients in whom the planning objectives with ITV concept cannot be met. However, tumor tracking requires a 4D dose calculation to have an accurate estimation of the doses to the organs at risk and additionally, a suitable position information and motion compensation system.
Four-dimensional computed tomography
Analytical anisotropic algorithm
Anterior - posterior
- D2 :
Dose received by the 2% of target volume
- D95 :
Dose received by the 2% of target volume
- D98 :
Dose received by the 2% of target volume
- D max :
- D mean :
- D min :
Gross tumor volume
Internal target volume
Magnetic resonance imaging
Organ at risk
Planning target volume
Stereotactic body radiation therapy
Superior - inferior
Simultaneous integrated protection
Volumetric-modulated arc therapy
Malvezzi M, Bertuccio P, Rosso T, Rota M, Levi F, La Vecchia C, Negri E. European cancer mortality predictions for the year 2015: does lung cancer have the highest death rate in EU women? Ann Oncol. 2015; 26(4):779–86.
Petrelli F, Comito T, Ghidini A, Torri V, Scorsetti M, Barni S. Stereotactic body radiation therapy for locally advanced pancreatic cancer: a systematic review and pooled analysis of 19 trials. Int J Radiat Oncol Biol Phys. 2017; 97:313–22.
Trakul N, Koong AC, Chang DT. Stereotactic body radiotherapy in the treatment of pancreatic cancer. Semin Radiat Oncol. 2014; 24:140–7.
Casamassima F, Cavedon C, Francescon P, Stancanello J, Avanzo M, Cora S, et al.Use of motion tracking in stereotactic body radiotherapy: Evaluation of uncertainty in off-target dose distribution and optimization strategies. Acta Oncol. 2006; 45:943–7.
Panje C, Andratschke N, Brunner TB, Niyazi M, Guckenberger M. Stereotactic body radiotherapy for renal cell cancer and pancreatic cancer: literature review and practice recommendations of the DEGRO working group on stereotactic radiotherapy. Strahlenther Onkol. 2016; 192:875–5.
Pollom EL, Alagappan M, von Eyben R, Kunz PL, Fisher GA, Ford JA, et al. Single- versus multifraction stereotactic body radiation therapy for pancreatic adenocarcinoma: outcomes and toxicity. Int J Radiat Oncol Biol Phys. 2014; 90:918–25.
Mellon EA, Hoffe SE, Springett GM, Frakes JM, Strom TJ, Hodul PJ, et al. Long-term outcomes of induction chemotherapy and neoadjuvant stereotactic body radiotherapy for borderline resectable and locally advanced pancreatic adenocarcinoma. Acta Oncol. 2015; 69:979–85.
Chuong MD, Springett GM, Freilich JM, Park CK, Weber JM, Mellon EA, et al. Stereotactic body radiation therapy for locally advanced and borderline resectable pancreatic cancer is effective and well tolerated. Int J Radiat Oncol Biol Phys. 2013; 86:516–22.
Feng M, Balter JM, Normolle D, Adusumili S, Cao Y, Chenevert TL, et al. Characterization of pancreatic tumor motion using cine MRI: surrogates for tumor position should be used with caution. Int J Radiat Oncol Biol Phys. 2009; 74:884–91.
Bussels B, Goethals L, Feron M, Bielen D, Dymarkowski S, Suetens P, et al. Respiration-induced movement of the upper abdominal organs: a pitfall for the three-dimensional conformal radiation treatment of pancreatic cancer. Radiother Oncol. 2003; 68:69–74.
Goldstein SD, Ford EC, Duhon M, McNutt T, Wong J, Herman JM. Use of respiratory-correlated four-dimensional computed tomography to determine acceptable treatment margins for locally advanced pancreatic adenocarcinoma. Int J Radiat Oncol Biol Phys. 2010; 76:597–602.
Cattaneo GM, Passoni P, Sangalli G, Slim N, Longobardi B, Mancosu P, et al. Internal target volume defined by contrast-enhanced 4D-CT scan in unresectable pancreatic tumour: Evaluation and reproducibility. Radiother Oncol. 2010; 97:525–29.
Shiinoki T, Shibuya K, Nakamura M, Nakamura A, Matsuo Y, Sawada A, et al. Interfractional reproducibility in pancreatic position based on four-dimensional computed tomography. Radiat Oncol Biol Phys. 2011; 80:1567–72.
Hallman JL, Mori S, Sharp GC, Lu HM, Hong TS, Chen GT. A four-dimensional computed tomography analysis of multiorgan abdominal motion. Int J Radiat Oncol Biol Phys. 2012; 83:435–1.
Ohira S, Isono M, Ueda Y, Hirata T, Ashida R, Takahashi H, et al. Assessment with cone-beam computed tomography of intrafractional motion and interfractional position changes of resectable and borderline resectable pancreatic tumours with implanted fiducial marker. Br J Radiol. 2017; 90:20160815.
Gierga DP, Chen GT, Kung JH, Betke M, Lombardi J, Willett CG. Quantification of respiration-induced abdominal tumor motion and its impact on imrt dose distributions. Int J Radiat Oncol Biol Phys. 2004; 58:1584–95.
Bhasin DK, Rana SS, Jahagirdar S, Nagi B. Does the pancreas move with respiration?Gastroenterol Hepatol. 2006; 21:1424–27.
Suramo I, Päivänsalo M, Myllylä V. Cranio-caudal movements of the liver, pancreas and kidneys in respiration. Acta Radiol Diagn (Stockh). 1984; 25:129–31.
Moningi S, Marciscano AE, Rosati L, Sook Ng K, Teboh Forbang R, Jackson J, et al. Stereotactic body radiation therapy in pancreatic cancer: the new frontier. Expert Rev Anticancer Ther. 2014; 14:1461–75.
Keall PJ, Mageras GS, Balter JM, Emery RS, Forster KM, Jiang SB, et al. The management of respiratory motion in radiation oncology report of AAPM task group 76. Med Phys. 2006; 33:3874–900.
Seuntjens J, Lartigau EF, Cora S, Ding GX, Goetsch S, Nuyttens J, et al.Prescribing, recording, and reporting of stereotactic treatments with small photon beams. J ICRU. 2014; 14:1–160. https://academic.oup.com/jicru/article/14/2/1/4035832/Report-91.
Dawson LA, Brock KK, Kazanjian S, Fitch D, McGinn CJ, Lawrence TS, et al. The reproducibility of organ position using active breathing control (ABC) during liver radiotherapy. Int J Radiat Oncol Biol Phys. 2001; 51:1410–21.
Murray B, Forster K, Timmerman R. Frame-based immobilization and targeting for stereotactic body radiation therapy. Med Dosim. 2007; 32:86–91.
Taniguchi CM, Murphy JD, Eclov N, Atwood TF, Kielar KN, Christman-Skieller C, et al. Dosimetric analysis of organs at risk during expiratory gating in stereotactic body radiation therapy for pancreatic cancer. Int J Radiat Oncol Biol Phys. 2013; 85:1090–5.
Ehrbar S, Jöhl A, Tartas A, Stark LS, Riesterer O, Klöck S, et al. Itv, mid-ventilation, gating or couch tracking - a comparison of respiratory motion-management techniques based on 4D dose calculations. Radiother Oncol. 2017; 124:80–8.
Kilby W, Dooley JR, Kuduvalli G, Sayeh S, Maurer C R J. The CyberKnife robotic radiosurgery system in 2010. Technol Cancer Res Treat. 2010; 9:433–52.
Depuydt T, Verellen D, Haas O, Gevaert T, Linthout N, Duchateau M, et al. Geometric accuracy of a novel gimbals based radiation therapy tumor tracking system. Radiother Oncol. 2011; 98:365–72.
Keall PJ, Colvill E, O’Brien R, Ng JA, Poulsen PR, Eade T, et al. The first clinical implementation of electromagnetic transponder-guided MLC tracking. Med Phys. 2014; 41:0207702.
Lang S, Zeimetz J, Ochsner G, Schmid Daners M, Riesterer O, Klöck S. Development and evaluation of a prototype tracking system using the treatment couch. Med Phys. 2014; 41:021720.
D’Souza WD, Naqvi SA, Yu CX. Real-time intra-fraction-motion tracking using the treatment couch: a feasibility study. Phys med Biol. 2005; 50:4021–33.
Papalazarou C, Klop GJ, Milder MTW, Marijnissen JPA, Gupta V, Heijmen BJM, et al. CyberKnife with integrated CT-on-rails: System description and first clinical application for pancreas SBRT. Med Phys. 2017; 44:4816–27.
Malinowski K, McAvoy TJ, George R, Dieterich S, D’Souza WD. Online monitoring and error detection of real-time tumor displacement prediction accuracy using control limits on respiratory surrogate statistics. Med Phys. 2012; 39:2042–8.
Malinowski K, McAvoy TJ, George R, Dieterich S, D’Souza WD. Incidence of changes in respiration-induced tumor motion and its relationship with respiratory surrogates during individual treatment fractions. Int J Radiat Oncol Biol Phys. 2012; 82:1665–73.
Zhang H, Zhao G, Djajaputra D, Xie Y. Determination of acquisition frequency for intrafractional motion of pancreas in CyberKnife radiotherapy. Sci World J. 2014; 2014:408019. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4053084/.
Rietzel E, Chen GT, Choi NC, Willet CG. Four-dimensional image-based treatment planning: target volume segmentation and dose calculation in the presence of respiratory motion. Int J Radiat Oncol Biol Phys. 2005; 61:1535–50.
Schlaefer A, Fisseler J, Dieterich S, Shiomi H, Cleary K, Schweikard A. Feasibility of four-dimensional conformal planning for robotic radiosurgery. Med Phys. 2005; 32:3786–92.
Chan MK, Kwong DL, Ng SC, Tam EK, Tong AS. Investigation of four-dimensional (4D) Monte Carlo dose calculation in real-time tumor tracking stereotactic body radiotherapy for lung cancers. Med Phys. 2012; 90:5479–87.
Gkika E, Adebahr S, Kirste S, Schimek-Jasch T, Wiehle R, Claus R, et al. Stereotactic body radiotherapy (SBRT) in recurrent or oligometastatic pancreatic cancer: A toxicity review of simultaneous integrated protection (SIP) versus conventional SBRT. Strahlenther Onkol. 2017; 193:433–3.
Ehrbar S, Lang S, Stieb S, Riesterer O, Stark LS, Guckenberger M, et al. Three-dimensional versus four-dimensional dose calculation for volumetric modulated arc therapy of hypo-fractionated treatments. J Med Phys. 2016; 26:45–53.
R Development Core Team. R: A language and environment for statistical computing. Vienna; 2017. https://www.r-project.org.
Tai A, Liang Z, Erickson B, Li XA. Management of respiration-induced motion with 4-dimensional computed tomography (4DCT) for pancreas irradiation. Int J Radiat Oncol Biol Phys. 2013; 86:908–13.
Lens E, Van der Horst A, Kroon PS, van Hooft JE, Fajardo RD, Fockens P, et al. Differences in respiratory induced pancreatic tumor motion between 4D treatment planning CT and daily cone beam CT, measured using intratumoral fiducials. Acta Oncol. 2014; 53:1257–64.
Minn AY, Schellenberg D, Maxim P, Suh Y, McKenna S, Cox B, et al. Pancreatic tumor motion on a single planning 4D-CT does not correlate with intrafraction tumor motion during treatment. Am J Clin Oncol. 2009; 32:364–8.
Lens E, van der Horst A, Versteijne E, van Tienhoven G, Bel A. Dosimetric advantages of midventilation compared with internal target volume for radiation therapy of pancreatic cancer. Int J Radiat Oncol Biol Phys. 2015; 92:675–82.
Guckenberger M, Wilbert J, Krieger T, Richter A, Baier K, Meyer J, et al. Four-dimensional treatment planning for stereotactic body radiotherapy. Int J Radiat Oncol Biol Phys. 2007; 69:276–85.
Chan M, Grehn M, Cremers F, Siebert FA, Wurster S, Huttenlocher S, et al. Dosimetric implications of residual tracking errors during robotic SBRT of liver metastases. Int J Radiat Oncol Biol Phys. 2017; 97:839–48.
Chan MK, Werner R, Ayadi M, Blanck O. Comparison of 3D and 4D Monte Carlo optimization in robotic tracking stereotactic body radiotherapy of lung cancer. Strahlenther Onkol. 2015; 191:161–71.
Availability of data and materials
The datasets used in the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
For this study anonymized patient data was used. According to our local ethics committee this does not require ethics approval.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.