- Open Access
Comparison of treatment plans between IMRT with MR-linac and VMAT for lung SABR
Radiation Oncology volume 14, Article number: 105 (2019)
The aim of this study was to compare the plan quality of magnetic-resonance image-based intensity modulated radiation therapy (MRI-based-IMRT) with the MRIdian Linac system to that of volumetric modulated arc therapy (VMAT) with the TrueBeam STx system for lung stereotactic ablative radiotherapy (SABR).
A total of 22 patients with tumors located in the lower lobe were retrospectively selected for the study. For each patient, both the MRI-based-IMRT and VMAT plans were generated using an identical CT image set and identical structures with the exception of the planning target volume (PTV). The PTVs of the MRI-based-IMRT were generated by adding an isotropic margin of 3 mm from the gross tumor volume, whereas those of VMAT were generated by adding an isotropic margin of 5 mm from the internal target volume. For both the MRI-based-IMRT and VMAT, the prescription doses to the PTVs were 60 Gy in four fractions.
The average PTV volume of the MRI-based-IMRT was approximately 4-times smaller than that of VMAT (p < 0.001). The maximum dose to the bronchi for the MRI-based-IMRT was smaller than that for the VMAT (20.4 Gy versus 24.2 Gy, p < 0.001). In addition, V40Gy of the rib for the MRI-based-IMRT was smaller than that for the VMAT (1.8 cm3 versus 7.7 cm3, p = 0.008). However, the maximum doses to the skin and spinal cord for the MRI-based-IMRT (33.0 Gy and 14.5 Gy, respectively) were larger than those for the VMAT (27.8 Gy and 11.0 Gy, respectively) showing p values of less than 0.02. For the ipsilateral lung, the mean dose, V20Gy, V10Gy, and V5Gy for the MRI-based-IMRT were smaller than those for the VMAT (all with p < 0.05). For the contralateral lung, V5Gy, V10Gy, D1500cc, and D1000cc for the MRI-based-IMRT were larger than those for the VMAT (all with p < 0.05). The mean dose and V50% of the whole body for the MRI-based-IMRT were smaller than those for the VMAT (0.9 Gy versus 1.2 Gy, and 78.7 cm3 versus 103.5 cm3, respectively, all at p < 0.001).
The MRI-based-IMRT using the MRIdian Linac system could reduce doses to bronchi, rib, ipsilateral lung, and whole body compared to VMAT for lung SABR when the tumor was located in the lower lobe.
Since the release of the MRIdian® Linac system (ViewRay Inc., Oakwood Village, OH), linear accelerators with a magnetic resonance imaging system (MR-linac) have become clinically available. The MRIdian Linac system can generate 6 MV flattening filter free (FFF) beams with a maximum dose rate of 6 Gy/min using an S-band standing wave linear accelerator (linac) . Therefore, the MRIdian Linac system generates treatment beams with higher penetrability and smaller penumbrae than those of the ViewRay® system (ViewRay Inc., Oakwood Village, OH), which is a previous model for magnetic resonance image-guided radiation therapy (MR-IGRT) that uses Co-60 radioisotopes to generate treatment beams [2, 3]. Although the MRIdian Linac system is equipped with a linac, 3D volumetric and 2D planar magnetic resonance images (MRI) can be acquired using a 0.35 T magnetic field, similar with the ViewRay system . The mutual interference that occurs between the linac and MR imaging systems of the MRIdian Linac system were eliminated using patented technology from ViewRay, Inc. . In addition to the successful integration of the linac and MR imaging systems, a notable feature of the MRIdian Linac system is a double-stacked and double-focused multi-leaf collimator (MLC) system that can project field sizes of 0.2 cm × 0.4 cm to 27.4 cm × 24.1 cm at the isoplane located 90 cm from the source . Using this MLC system, the MRIdian Linac system can conduct step-and-shoot intensity modulated radiation therapy (IMRT) with fluences at a finer resolution than those of the ViewRay system, thereby potentially improving the quality of the treatment plan [4, 5].
One of the indications of MR-IGRT is lung cancer owing to the substantial respiratory motions of the lung tumors . Because the MR-IGRT technique can be used in respiratory gated radiation therapy based on the application of real-time cine planar MR images during treatment without additional imaging doses, the margins of the lung target volumes can be significantly reduced, which results in substantially smaller doses to normal tissue . This advantage can be maximized in stereotactic ablative radiotherapy (SABR) for lung cancer because the fractional doses of lung SABR are much larger than those of conventional fractionated radiation therapy and should be delivered to patients accurately for lung SABR . Several previous studies have compared the dose distributions of MR-IGRT using the ViewRay system to those of linac-based IMRT or linac-based volumetric modulated arc therapy (VMAT) and reported that the plan quality of the ViewRay system is clinically acceptable but no better than that of linac-based IMRT and VMAT owing to the inferior beam quality of the Co-60 source and the large MLC leaf width of the ViewRay system [2, 3, 8,9,10,11]. Specifically for lung target volumes of smaller than 10 cm3, the plan quality of internal target volume (ITV) based VMAT is statistically significantly better than that of gross tumor volume (GTV) based IMRT using the ViewRay system despite the margin reduction capability of the ViewRay system .
Because the MRIdian Linac system overcomes the disadvantages of the ViewRay system, including the poor characteristics of the Co-60 source beams and the large MLC leaf width, the plan quality of the MRIdian Linac system might be better than that of a conventional IMRT or VMAT by virtue of its margin reduction capability. In contrast, the MRIdian Linac system can only generate coplanar step-and-shoot IMRT plans whose plan quality might be inferior to that of VMAT. It is therefore unclear whether the lung SABR plan quality of the MRIdian Lianc system is better than that of VMAT. Thus, in the present study, we compared lung SABR IMRT plans using the MRIdian Linac system based on the GTV to those of VMAT plans using TrueBeam STx™ (Varian Medical Systems, Palo Alto, CA) based on the ITV for a total of 22 patients with lung cancer.
Patient selection and simulation
After receiving approval from the institutional review board, a total of 22 patients treated for lung cancer with the SABR VMAT technique were retrospectively selected for this study. Because the respiratory motion of lung tumors located in the lower lobe is known to be large, we selected only patients with lung tumors located in the lower lobe . All patients were scanned using the Brilliance CT Big Bore™ (Phillips, Cleveland, OH). Ten phase 4D CT images were acquired using the Real-time Position Management™ (RPM, Varian Medical Systems, Palo Alto, CA) system for each patient to define the ITV for VMAT planning. The slice thickness of the CT images was 2 mm. When acquiring the CT images, all patients were immobilized using the Body Pro-Lok™ (CIVICO, Orange City, IA) system to reduce the respiratory motion of the lung tumors.
Treatment planning of VMAT using TrueBeam STx
The planning target volume (PTV) was defined by adding an isotropic margin of 5 mm from the ITV in the 80% phase CT image set . The prescription dose to the PTV was 60 Gy in four fractions for all patients. VMAT plans were generated using a 6 MV FFF beam of TrueBeam STx with a High-Definition 120 MLC™ (Varian Medical Systems, Palo Alto, CA). For the planning, the Eclipse™ system (Varian Medical Systems, Palo Alto, CA) was used. To minimize the doses to the contralateral lung, two coplanar half arcs near the ipsilateral lung were used. One half arc was rotated from 180° to 0°, whereas the other half arc was rotated from 0° to 180°. All VMAT plans were optimized using the progressive resolution optimizer (PRO3, ver. 13, Varian Medical Systems, Palo Alto, CA). The dose-volume constraints during optimization were set according to the National Comprehensive Cancer Network Clinical Practice Guidelines in Oncology (NCCN Guidelines®, National Comprehensive Cancer Network, Fort Washington, PA). After optimization, the dose distributions were calculated using the Acuros® XB algorithm (AXB, ver. 13, Varian Medical Systems, Palo Alto, CA). The size of the dose calculation grid was 1 mm, which is the finest resolution available with the Eclipse system. After a dose calculation, all plans were normalized to cover 95% of the PTV volume with 100% of the prescription dose.
Treatment planning of IMRT using MRIdian Linac system
Because this is a comparative planning study, CT images identical to those of VMAT plans were used for the IMRT planning with the MRIdian Linac system. The PTV was defined by adding an isotropic margin of 3 mm from the GTV in the 80% phase CT image set because the MRIdian Linac system is capable of anatomy-based gated radiation therapy. The PTV was used as the gating window during treatment, i.e., the treatment beams were delivered to a patient only when the GTV was inside the PTV. We added a 3-mm isotropic PTV margin empirically to consider the daily variations in lung tumor movements of the patients. The prescription dose to the PTV was identical to that of the VMAT plans, which was 60 Gy in the four fractions. Coplanar step-and-shoot IMRT plans were generated for each patient with a 6 MV FFF photon beam of the MRIdian Linac system using double-stacked and double-focused MLCs. For the planning, the treatment planning system (TPS) of the MRIdian Linac system which is the MRIdian system (ViewRay Inc., Oakwood Village, OH) was used. Eight or nine fields with various gantry angles were used for each patient according to the tumor locations. The isocenter was not moved laterally to avoid collisions between the patient body and bore of the MRIdian Linac system. All plans were optimized to meet the dose-volumetric constraints based on the same NCCN guidelines as used in the VMAT planning. The sequencing number of segments was set to 60 . The IMRT efficiency was set to 3.00, which determines the smoothness of the fluence map . The dose distributions were calculated using a dose calculation grid of 2 mm. When calculating the dose distributions, the number of histories per cm2 was set to 25,000. After the dose calculation, the calculated dose distributions were imported from the MRIdian system to the Eclipse system in order to calculate the dose-volumetric parameters and dose volume histograms (DVHs) . All plans were normalized to cover 100% of the PTV volume with 99.99% of the prescription dose because the PTV of the MRIdian Linac system was generated without considering the internal target motion through respiration.
Evaluation of the treatment plans
For each plan, the total monitor unit (MU) and calculated beam-on time in the TPS were acquired. For the IMRT plans using the MRIdian Linac system, the average number of segments for each plan was acquired.
For the PTVs, the maximum, minimum, and mean doses of the PTV and PTV volumes were calculated. The minimum doses to 2% volume of the PTV (D2%), D5%, D95%, and D98% were also calculated. The homogeneity index (HI) was calculated as follows .
For organs at risk (OARs), the maximum doses to the bronchi, heart, rib, skin, and spinal cord were calculated. The percent volume received at least 34.8 Gy (V34.8Gy) for the bronchi, V30Gy for the esophagus, V34Gy for the heart, V40Gy for the rib, V36Gy for the skin, and V26Gy for the spinal cord were also calculated. For both the ipsilateral and contralateral lungs, V20Gy, V10Gy, V5Gy, D1500cc, D1000cc, and mean doses were calculated.
For whole body of a patient, mean dose to whole body as well as the percent volume received at least 50% of the prescription dose (V50%) were acquired.
A paired t-test was conducted to examine the statistical significances of the differences in the parameters between VMAT and IMRT using the MRIdian Linac system. We regarded differences with p values of less than 0.05 as statistically significant.
Differences in the plan parameters
The plan parameters of both VMAT and IMRT using an MR-linac are shown in Table 1. The MU of the MRI-based IMRT was twice that of VMAT on average (8373 MU for MRI-based IMRT versus 4173 MU for VMAT, p < 0.001). The calculated beam-on time of the MRI-based IMRT was approximately 4-times larger than that of VMAT on average (18.4 min for MRI-based IMRT versus 4.8 min for VMAT, p < 0.001).
Dose-volumetric parameters of PTV
The dose-volumetric parameters of the PTVs of both the VMAT and MRI-based IMRT plans are shown in Table 2. The dose distributions of VMAT and MRI-based IMRT for the representative patients are shown in Fig. 1. The DVHs of the representative patients are shown in Fig. 2.
The PTV volume of the MRI-based IMRT is approximately 4-times smaller than that of VMAT on average (p < 0.001). The average values of D98%, D95%, and the minimum dose of MRI-based IMRT were larger than those of VMAT (all at p < 0.02) while the average values of D5%, D2%, and the maximum dose of MRI-based IMRT were smaller than those of VMAT (all at p < 0.04). Therefore, HI indicates that the dose homogeneity inside the PTV volume of MRI-based IMRT was better than that of VMAT (0.062 for MRI-based IMRT versus 0.099 for VMAT, p < 0.001).
Dose-volumetric parameters of OARs
The dose-volumetric parameters of OARs are shown in Table 3.
The maximum dose to the bronchi of MRI-based IMRT was smaller than that of VMAT with a statistical significance (20.4 Gy for MRI-based IMRT versus 24.2 Gy for VMAT, p < 0.001). Both the maximum dose and V34Gy for the heart were smaller in the MRI-based IMRT plans than in the VMAT plans; however, the differences were not statistically significant (all at p > 0.05). The maximum dose and V40Gy of the rib of MRI-based IMRT (43.0 Gy and 1.8 cm3, respectively) were smaller than those of VMAT (47.2 Gy and 7.7 cm3, respectively) showing p values less than 0.04, whereas the maximum doses to the skin and spinal cord of MRI-based IMRT (33.0 and 14.5 Gy, respectively) were larger than those of VMAT (27.8 and 11.0 Gy, respectively) showing p values of less than 0.02.
For the ipsilateral lung, the mean dose, V20Gy, V10Gy, and V5Gy of the MRI-based IMRT were smaller than those of VMAT (all at p < 0.05), whereas D1500cc and D1000cc of the MRI-based IMRT were larger than those of VMAT (all at p < 0.02). For the contralateral lung, V10Gy, V5Gy, D1500cc, and D1000cc of the MRI-based IMRT were larger than those of VMAT with a statistical significance (all at p < 0.05).
Dose-volumetric parameters of the whole body
The dose-volumetric parameters of the whole body are shown in Table 3. The mean dose and V50% consistently indicated that the irradiation of normal tissue of MRI-based IMRT was lower than that of VMAT with a statistical significance (all at p < 0.001).
In the present study, we compared GTV-based IMRT plans generated with the MR-linac to the ITV-based VMAT plans for lung SABR. We retrospectively selected patients with tumors located at lower lobes of the lungs, of which the respiratory motions were large . Therefore, the PTV sizes of VMAT, including margins compensating for the respiratory motion of the lung tumors, were much larger than those of MRI-based IMRT (approximately 4.3-times larger), which did not include the margins for respiratory motion. By doing so, we were able to maximize the advantage of MRI-based IMRT, which is internal anatomy-based gated radiotherapy with real-time cine MR images. Despite the smaller PTV sizes of MR-IGRT than those of conventional ITV-based radiotherapy, a previous study demonstrated that the lung SABR plan quality of the ViewRay system was not better than that of VMAT owing to the poor characteristics of Co-60 beams (large penumbrae and low penetrability) and the large MLC width of the ViewRay system . Specifically for target volumes of smaller than 10 cm3, no benefits were observed when utilizing the ViewRay system for lung SABR as compared to VMAT . However, in the present study, the MR-IGRT plans with the MR-linac demonstrated superiority to the VMAT plans even for target volumes of smaller than 10 cm3, showing less irradiation of normal tissue with intermediate and high doses than in the VMAT plans. As the results show, the volumes of body irradiated by 50% of the prescription dose of MRI-based IMRT were approximately 76% that of VMAT. In addition, the MRI-based IMRT was generally able to reduce the dose-volumetric parameters of the bronchi, rib, heart, and ipsilateral lung as compared to VMAT. This better plan quality was attributed to the margin reduction capability of MRI-based IMRT combined with the improved photon beam characteristics generated using the linac and small leaf widths of the MLCs.
Although most of the clinically relevant dose-volumetric parameters of MRI-based IMRT were better than those of VMAT, this was not always the case. Regarding the D1500cc and D1000cc of the ipsilateral lung, and all dose-volumetric parameters of the contralateral lung examined in this study, the maximum dose to the skin, and the maximum dose to the spinal cord from the MR-IGRT plans were worse than those of the VMAT plans. This might be attributed to the different characteristics of the IMRT and VMAT technique as well as the different beam arrangements of their plans, as described in the present study [13, 14]. For D1500cc and D1000cc of the ipsilateral lung, the average differences between MRI-based IMRT and VMAT were less than only 0.15 Gy, which is insignificant. For the contralateral lung, VMAT used half arcs proximal to the ipsilateral lung, and therefore no beams were delivered through the contralateral lung. In contrast, for MRI-based IMRT plans, beams were occasionally delivered to the target volume through the contralateral lung to acquire the acceptable target conformity, as shown in Fig. 1 (b). This resulted in relatively higher dose-volumetric parameters of the contralateral lung in the MRI-based IMRT plans than in the VMAT plans. In the case of skin, because the VMAT technique delivers beams through arcs, i.e., numerous beams in various directions are utilized, the doses can be distributed all over the skin but not concentrated within specific regions . This can reduce the maximum doses to the skin for VMAT plans. In the case of the spinal cord, because MRI-based IMRT occasionally uses a beam through the spinal cord according to the target volume location, as shown in Fig. 1 (d), this could increase the maximum dose to the spinal cord.
The required MUs and beam-on times of MRI-based IMRT were much larger than those of VMAT. According to a simple calculation, a beam-on time of approximately 4-times longer on average was required for MRI-based IMRT as compared to VMAT, which increases the patient discomfort. In addition, the respiratory gated beam delivery of MRI-based IMRT would increase the treatment time. However, the patient setup of MRI-based IMRT is relatively simpler than that of VMAT for lung SABR because rigorous immobilization devices used to suppress the respiratory motion of the patient to reduce the ITV margins are not necessarily applied to a patient in the case of MRI-based IMRT. This simple setup could increase the degree of patient comfort. Therefore, MRI-based IMRT has both advantages and disadvantages in terms of patient comfort. Besides the patient comfort issue, for the long beam-on time of MRI-based IMRT, there might be a concern regarding an inaccurate treatment when the treatment time is long. However, a previous study reported that a long treatment time is not correlated with an inaccurate treatment for lung SABR . Moreover, because the treatment beams of MRI-based IMRT are delivered to patients when monitoring their internal anatomy motion in real time, the long treatment time of MRI-based IMRT will not hamper the treatment accuracy.
Several studies comparing MRI-based IMRT using the ViewRay system to linac-based IMRT or VMAT for lung SABR have reported that the former are clinically acceptable but generally worse than those of linac-based IMRT or VMAT [2, 3, 8,9,10,11]. However, no studies have yet been conducted to compare MRI-based IMRT plans using the MRIdian Linac system to VMAT plans for lung SABR. In the present study, we conducted a planning study comparing MRI-based IMRT plans using the MRIdian Linac system to VMAT plans for lung SABR when the tumors were located in the lower lob. We presented that the MRI-based IMRT could reduce doses to bronchi, rib, ipsilateral lung, and whole body compared to VMAT and this was attributed to the combination of 6 MV FFF photon beams generated from the MR-linac, a double-stacked and double-focus MLC system with small leaf widths, and the margin reduction capability using a real-time anatomy-based gating technique.
In the present study, we identified the MRI-based IMRT with the MRIdian Linac system could improve the dose homogeneity inside the target volume as well as reduce doses to bronchi, rib, ipsilateral lung, and whole body compared to VMAT for lung SABR when the tumors were located in the lower lobe. In addition, MRI-based IMRT with the MRIdian Linac system is capable of internal anatomy-based gated radiotherapy by monitoring real-time tumor motion, which guarantees an accurate treatment, but currently cannot be conducted using the VMAT technique. Therefore, it seems appropriate to utilize MRI-based IMRT for lung SABR when the tumors were located in the lower lobe.
Dose volume histogram
Flattening filter free
Gross tumor volume
Intensity modulated radiation therapy
Internal target volume
Magnetic resonance image
Magnetic resonance image-guided radiation therapy
Linear accelerator with magnetic resonance imaging system
National comprehensive cancer network
Organ at risk
Progressive resolution optimizer 3
Planning target volume
Real-time Position Management
Stereotactic ablative radiotherapy
Treatment planning system
Volumetric modulated arc therapy
Ramey SJ, Padgett KR, Lamichhane N, Neboori HJ, Kwon D, Mellon EA, Brown K, Duffy M, Victoria J, Dogan N, Portelance L. Dosimetric analysis of stereotactic body radiation therapy for pancreatic cancer using MR-guided tri-60Co unit, MR-guided LINAC, and conventional LINAC-based plans. Pract Radiat Oncol. 2018;8:e312–21.
Park JM, Park SY, Kim HJ, Wu HG, Carlson J, Kim JI. A comparative planning study for lung SABR between tri-co-60 magnetic resonance image guided radiation therapy system and volumetric modulated arc therapy. Radiother Oncol. 2016;120:279–85.
Choi CH, Park SY, Kim JI, Kim JH, Kim K, Carlson J, Park JM. Quality of tri-co-60 MR-IGRT treatment plans in comparison with VMAT treatment plans for spine SABR. Br J Radiol. 2017;90(1070):20160652.
Rodal J, Sovik A, Malinen E. Influence of MLC leaf width on biologically adapted IMRT plans. Acta Oncol. 2010;49:1116–23.
Park JM, Park SY, Kim JH, Carlson J, Kim JI. The effect of extremely narrow MLC leaf width on the plan quality of VMAT for prostate cancer. Radiat Oncol. 2016;11:85.
Seppenwoolde Y, Shirato H, Kitamura K, Shimizu S, van Herk M, Lebesque JV, Miyasaka K. Precise and real-time measurement of 3D tumor motion in lung due to breathing and heartbeat, measured during radiotherapy. Int J Radiat Oncol Biol Phys. 2002;53:822–34.
Yeo SG, Kim ES. Efficient approach for determining four-dimensional computed tomography-based internal target volume in stereotactic radiotherapy of lung cancer. Radiat Oncol J. 2013;31:247–51.
Saenz DL, Paliwal BR, Bayouth JE. A dose homogeneity and conformity evaluation between ViewRay and pinnacle-based linear accelerator IMRT treatment plans. J Med Phys. 2014;39:64–70.
Kishan AU, Cao M, Wang PC, Mikaeilian AG, Tenn S, Rwigema JC, Sheng K, Low DA, Kupelian PA, Steinberg ML, Lee P. Feasibility of magnetic resonance imaging-guided liver stereotactic body radiation therapy: a comparison between modulated tri-cobalt-60 teletherapy and linear accelerator-based intensity modulated radiation therapy. Pract Radiat Oncol. 2015;5:330–7.
Wooten HO, Green O, Yang M, DeWees T, Kashani R, Olsen J, Michalski J, Yang D, Tanderup K, Hu Y, et al. Quality of intensity modulated radiation therapy treatment plans using a 60Co magnetic resonance image guidance radiation therapy system. Int J Radiat Oncol Biol Phys. 2015;92:771–8.
Merna C, Rwigema JC, Cao M, Wang PC, Kishan AU, Michailian A, Lamb J, Sheng K, Agazaryan N, Low DA, et al. A treatment planning comparison between modulated tri-cobalt-60 teletherapy and linear accelerator-based stereotactic body radiotherapy for central early-stage non-small cell lung cancer. Med Dosim. 2016;41:87–91.
Kim JI, Han JH, Choi CH, An HJ, Wu HG, Park JM. Discrepancies in dose-volume histograms henerated from different treatment planning systems. J Radiat Prot and Res. 2018;43(2):59–65.
Park JM, Kim K, Chie EK, Choi CH, Ye SJ, Ha SW. RapidArc vs intensity-modulated radiation therapy for hepatocellular carcinoma: a comparative planning study. Br J Radiol. 2012;85:e323–9.
Holt A, van Vliet-Vroegindeweij C, Mans A, Belderbos JS, Damen EM. Volumetric-modulated arc therapy for stereotactic body radiotherapy of lung tumors: a comparison with intensity-modulated radiotherapy techniques. Int J Radiat Oncol Biol Phys. 2011;81:1560–7.
Li W, Purdie TG, Taremi M, Fung S, Brade A, Cho BC, Hope A, Sun A, Jaffray DA, Bezjak A, Bissonnette JP. Effect of immobilization and performance status on intrafraction motion for stereotactic lung radiotherapy: analysis of 133 patients. Int J Radiat Oncol Biol Phys. 2011;81:1568–75.
Availability of the data and materials
All data generated or analyzed during the current study are available from the corresponding author upon reasonable request.
This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIP) (No. 2017M2A2A7A02020640, 2017M2A2A7A02020641, and 2017M2A2A7A02020643).
Ethics approval and consent to participate
Approval for this study was obtained from the institutional review board of Seoul National University Hospital (IRB No. 1901–059-1002).
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Park, J.M., Wu, H., Kim, H.J. et al. Comparison of treatment plans between IMRT with MR-linac and VMAT for lung SABR. Radiat Oncol 14, 105 (2019). https://doi.org/10.1186/s13014-019-1314-0
- Planning study