Figures
Abstract
Background and Purpose
Carbon-ion radiotherapy of prostate cancer is challenging in patients with metal implants in one or both hips. Problems can be circumvented by using fields at oblique angles. To evaluate the influence of setup and range uncertainties accompanying oblique field angles, we calculated rectal dose changes with oblique orthogonal field angles, using a device with fixed fields at 0° and 90° and a rotating patient couch.
Material and Methods
Dose distributions were calculated at the standard angles of 0° and 90°, and then at 30° and 60°. Setup uncertainty was simulated with changes from −2 mm to +2 mm for fields in the anterior-posterior, left-right, and cranial-caudal directions, and dose changes from range uncertainty were calculated with a 1 mm water-equivalent path length added to the target isocenter in each angle. The dose distributions regarding the passive irradiation method were calculated using the K2 dose algorithm.
Results
The rectal volumes with 0°, 30°, 60°, and 90° field angles at 95% of the prescription dose were 3.4±0.9 cm3, 2.8±1.1 cm3, 2.2±0.8 cm3, and 3.8±1.1 cm3, respectively. As compared with 90° fields, 30° and 60° fields had significant advantages regarding setup uncertainty and significant disadvantages regarding range uncertainty, but were not significantly different from the 90° field setup and range uncertainties.
Citation: Kubota Y, Kawamura H, Sakai M, Tsumuraya R, Tashiro M, Yusa K, et al. (2016) Changes in Rectal Dose Due to Alterations in Beam Angles for Setup Uncertainty and Range Uncertainty in Carbon-Ion Radiotherapy for Prostate Cancer. PLoS ONE 11(4): e0153894. https://doi.org/10.1371/journal.pone.0153894
Editor: Shian-Ying Sung, Taipei Medical University, TAIWAN
Received: October 18, 2015; Accepted: April 5, 2016; Published: April 20, 2016
Copyright: © 2016 Kubota et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper.
Funding: The authors have no support or funding to report.
Competing interests: The authors have declared that no competing interests exist.
Introduction
As compared with photon beams, particle beams provide sharper dose distributions by taking advantage of the Bragg peak and a sharp lateral penumbra [1]. Avoiding excessive exposure to organs at risk (OARs) requires understanding of the influences of setup error and of beam range error.
In particle therapy for prostate cancer, 90° horizontal fields are often used to reduce the rectal dose. It is easy to reduce the rectal dose by collimation without considering changes in the specific range of the particle beam, because small changes in range do not significantly affect the dose. However, horizontal fields are contraindicated after hip surgery employing metal plates or prostheses because of the unpredictability of the beam’s path through the metal and the influence of artifacts. For example, Jäkel et al. reported that in the case of tungsten and steel, metal path range errors of −5% and −18%, respectively, were observed, together with 1% path range errors involving the artifacts from titanium and steel [2]. Although artifact unpredictability was small if the metal was light, path unpredictability through the metal was large. Therefore, it is preferable to use 0° vertical (perpendicular to the patient body surface) or oblique fields in these cases. The oblique field might have a different sensitivity to the horizontal field regarding inaccuracies in patient setup and beam range; however, its influence is not well defined. Tang et al. and Christodouleas et al. reported a comparison of dose distribution in the anterior-oriented fields used for proton therapy; but, they did not consider the uncertainties involved [3,4]. Inter/intra fractional motion changes of the prostate could have an effect. However, only the influences of inaccuracies in setup and beam range were evaluated in this study. Although cases involving cancer patients with metal implants are not frequent, it is important to determine their influence on the oblique fields, because this has the potential to reduce uncertainty regarding the rectal dose using the current treatment.
The polybinary calibration method between CT density value and effective density for particle beam radiotherapy has an accuracy of 99% [5,6]. The resulting beam range uncertainties cause dose deviations that can result in errors in dose to the clinical target volume (CTV) and OARs located along or near the beam path. OARs located lateral to the target may be exposed to higher doses as a result of setup errors. Although robust optimizations of treatment planning including setup and range uncertainty for proton therapy have been proposed [7,8], the influences of dose distribution per field angle have not been considered.
We evaluated the influence of setup and range uncertainties on the rectal and CTV dose distribution of oblique fields as compared with a horizontal (90°) field in prostate cancer. Although the bladder dose might also change for each field angle, the rectal dose was focused in our study to simplify the problem because the bladder is unlikely to be a clinical problem.
Materials and Methods
Patients
We retrospectively studied the data of ten prostate cancer patients aged 59–74 years with a median age of 69.5 years. Three patients had a titanium metal hip implant. These were located on the left side in two patients and on the right side in one patient; seven patients did not have hip implants. The CTV includes the prostate and proximal seminal vesicle (SV), and rectal volume measured from CT images was 18.0–97.2 cm3 and 48.5–84.7 cm3 with medians of 36.1 cm3 and 70.7 cm3, respectively. Patient information, CTVs, and rectal volumes are detailed in Table 1. This study was approved by the Institutional Review Board at Gunma University Hospital (approval number: 1310), and patient records/information were anonymized and de-identified prior to analysis.
CTV shows the clinical target volume, metal implant shows which side the patient has in or not.
Irradiation devices and treatment planning
The Gunma University Heavy Ion Medical Center (GHMC) provides carbon-ion therapy [9] using a heavy ion irradiation device (Mitsubishi Electric, Tokyo, Japan) with a passive irradiation method [10] and a treatment planning system (TPS) (XiO-N, Mitsubishi Electric). The passive irradiation field was generated using the scatterer and wobbling, and the field was collimated to the outside of the PTV using a multi-leaf collimator (MLC). X-ray CT (Acquilion LB, Self-Propelled, Toshiba Medical Systems, Tochigi, Japan) images were acquired with non-helical, 2.0 mm × 4 acquisitions, full reconstruction mode, and pixel spacing was 1.07 × 1.07 mm. The average number of CT slices for prostate cancer patients was approximately 140. XiO-N incorporates a dose engine for ion beam radiotherapy dose calculations (K2-Dose) [11–14]. The relative biological effectiveness (RBE) was included in the absorbed dose using a spread-out Bragg peak concept [15], and the clinical dose including this was defined as Gy(RBE). This RBE concept was incorporated into the XiO-N. The planning target volume (PTV) for prostate cancer was created by adding the anterior and lateral margins of 10 mm, cranial and caudal margins of 6 mm, and a posterior margin of 5 mm to the CTV, but lateral margins to the proximal SV were 10 mm. Carbon ion treatment plans were generated as each PTV was covered with 95% of the prescribed dose. In treating prostate cancer, we used five fields, and the number of fractions for each field was normally 3, 3, 3, 4, and 3 or 3, 3, 3, 3, and 4 (16 fractions total). Thus, for one fraction, we used 3.6 Gy(RBE);and the total dose was 3.6 × 16 = 57.6 Gy(RBE).
In this planning study, the two patterns of CT image sets shown in Fig 1 were used for the calculation of dose distribution to evaluate the influence of dose deviations in each field angle, and to evaluate the actual fields used for the treatment. The first pattern was seven CT datasets for patients that had no implants as shown in Fig 1(A), and three CT datasets for patients with hip implants but with the opposite side of the implant as shown in Fig 1(B). Four different field angles (0°, 30°, 60°, and 90°) in each image set were used, with the patient couch rotated accordingly; the beam parameters used in the planning for each field angle are detailed in Table 2. The second pattern was three CT image sets for patients with implants, using oblique fields as shown in Fig 1(C). The field angles used for P1, P2, and P3 in treatment planning were 60°, 67.8°, and −35°, respectively. A prescription dose in all of the directional fields shown in Fig 1(A), 1(B) and 1(C) was set to 10.8 Gy(RBE), corresponding to three fractions per field.
Arrows show the beam directions, blue regions show the CTV, and red regions show the metal implant. (a) Diagram of a patient with no implant and a beam that can enter from the left (negative angle, gray arrows) or the right (positive angle, white arrows). (b) Diagram of a patient with a hip implant, showing the field directions of 0°, 30°, 60°, and 90°. (c) Diagram of a patient with a hip implant, showing the oblique field avoiding the implant. (d) Field directions from −90° to 90°; 90° represents left horizontal, and −90° represents right horizontal.
Creating setup uncertainty
In patient positioning, both orthogonal (frontal and lateral) X-ray images and digitally reconstructed radiographs from CT images are used, with the bony structures for landmarks [16]. We employed a 2-mm setup tolerance [17]. The QA mode of the TPS was used to evaluate the results of setup uncertainty, calculating the dose distributions after moving the field center from -2 mm to 2 mm in anterior-posterior (AP), left-right (LR), and cranial-caudal (CC) directions. The dose calculations for the evaluation of setup uncertainty were performed in four field angles for seven patients as shown in Fig 1(A), in four field angles for three patients as shown in Fig 1(B), and in each field angle for three patients as shown in Fig 1(C).
Creating range uncertainty
The stopping power ratios of the planning volumes were calculated with the polybinary calibration method using the CT density measurement/stopping power ratio [5,6]. The K2 dose used the ratio for calculating dose. Given the 99% accuracy of this method, we evaluated the uncertainty range using the following equation: (1) where RBody is the alteration in path through the patient’s body estimated from the range uncertainty, and RBeam is the alteration in the path the carbon beam travels before hitting the patient’s body surface. In this planning study, RBody was set to 2% of a water equivalent path length from the patient surface to the isocenter (IC) and RBeam was set to 1 mm from the specifications of our accelerator. The dose distributions with the range uncertainty were recalculated by changing the parameters of the range shifter (RSF) in four field angles for seven patients as shown in Fig 1(A), in four field angles for three patients as shown in Fig 1(B), and in each field angle for three patients as shown in Fig 1(C).
Estimation of setup and range uncertainties
Setup and range uncertainties were simulated by simultaneously changing the field center (along the worst-case direction in the AP, LR, and CC directions) and RSF parameters to construct a worst-case scenario; their dose distributions were calculated in four field angles for seven patients as shown in Fig 1(A), in four field angles for three patients as shown in Fig 1(B), and in each field angle for three patients as shown in Fig 1(C). A calculated case without considering setup and range uncertainties was defined as a normal-case, the highest mean rectal dose was defined as the worst-case, and the lowest mean rectal dose was defined as the best-case in each combination of the setup and range uncertainties.
Evaluation method
To evaluate the influence of the dose deviations due to the uncertainties in each field angle, we used a mean dose increase ratio RInc defined as (2) where Dmean,N is the mean rectal dose in the normal scenario, and Dmean,W is the mean rectal dose in the worst-case scenario.
Additionally, to evaluate the rectal dose volumes as a result of the uncertainties in each field angle, we used rectal 10, 50, and 95% volumes regarding the prescription dose (defined as V10, V50, and V95) in the normal, best, and worst cases for each field angle.
The RInc results for setup uncertainty and for range uncertainty were analyzed using the Wilcoxon test, and both of the RInc results for setup and range uncertainties and the results of the rectal dose volumes in the normal case, were analyzed using the Shapiro–Wilk normality test to determine if data were normally distributed, and using Dunnett’s multiple test. The level of statistical significance in the Wilcoxon and Dunnett’s multiple tests was set to 5%.
Results
The dose distribution of one case with a right hip implant is shown in Fig 2. For this patient treatment was as follows: three fractions using vertical fields, three fractions using horizontal fields from the left, two fractions using −67.8° fields from the right, five fractions using horizontal boost fields from the left, and two fractions using −67.8° boost fields from the right. All doses were 3.6 Gy(RBE) per fraction.
The upper row shows CT images and the bottom row shows CT images together with the dose distribution. Left column shows axial images, middle column shows sagittal images, and right column shows coronal images. Red lines show the metal implant after hip replacement. Green line shows prostate, light yellow line shows PTV, magenta line shows rectum, and purple line shows the bladder.
Dose distributions for four field angles in one patient without the implant are shown in Fig 3. RInc graphs from the uncertainties for ten patients are shown in Fig 4.
Data on one patient without a metal implant, for four field angles: (a) 0° field, (b) 30° field, (c) 60° field, and (d) 90° field. Green line shows prostate, light yellow line shows PTV, and magenta line shows rectum. (i) Dose distribution in the normal case. (ii) Yellow line shows the 95% isodose line for the prescription dose in the normal case, blue line shows the 95% isodose line of the prescription dose in the worst case, and red line shows the 95% isodose line of the prescription dose in the best case.
(a) is the increasing ratio from setup uncertainty in anterior-posterior (AP), left-right (LR), and cranial-caudal (CC) directions, (b) is the ratio from range uncertainty, and (c) is the ratio from the setup and range uncertainties. The error bars represent the standard deviations for 10 patients. * in (a), and (b) show p < 0.05 using the Wilcoxon test, and * in (c) shows p < 0.05 using Dunnett’s multiple test.
Dose volume histogram (DVH) graphs for the rectal dose and the CTV from setup and range uncertainties for ten patients, and DVH graphs in three patient cases with hip implants (Fig 1(C)) are shown in Fig 5; V10, V50, and V95 in normal, best, and worst cases for each field angle are shown in Table 3. In the normal cases (P1, P2, and P3) with the implant (Fig 1(C)), the V10 was 16.2 cm3, 18.4 cm3, and 25.7 cm3, respectively; the corresponding V50 was 5.5 cm3, 8.2 cm3, and 7.1 cm3, respectively, and the corresponding V95 was 1.4 cm3, 2.6 cm3, and 3.6 cm3, respectively. Additionally, RInc from the setup and range uncertainties for P1, P2, and P3 with the implant shown in Fig 1(C) was 25%, 33.1%, and 24.1%, respectively.
Red lines are the DVHs of the CTV dose shown as relative volume (%) and blue lines are DVHs of the rectal dose shown as absolute volume (cm3). (i) Ten patients in each beam angle. The light blue error bars represent the standard deviations for 10 patients. (ii) (e) is patient 1 with a 60° field. (f) is patient 2 with a 68° field. (g) is patient 3 with a −35° field. Patients in (e), (f) and (g) have hip implants and all fields avoid the implants. The solid lines show normal cases, and the dashed lines show the best or worst cases for setup and range uncertainties.
The values are the mean and standard deviation for 10 patients.
Discussion
The influences of field angles on rectal dose
The dose profile of the 0° field is affected by the depth and direction of the field, the dose profile of the 90° field is affected by the lateral direction of the field, and the dose profile of the 30° and 60° fields are affected by both depth and lateral directions. Considering the 0° fields, the rectum posterior to the PTV, is affected by the distal dose fall-off of the spread-out Bragg peak. Therefore, V10 from the 0° field was significantly larger than V10 from the 90° field angles shown in Table 3 because of the distal tail. With the 90° fields, the rectum, lateral to the PTV for the field’s eye view, is affected by the lower lateral penumbra dose. The 30° and 60° fields increase the dose to the rectum by both effects. Therefore, the V50 from the 90° field is significantly larger than the V50 from the other field angles shown in Table 3 because of the lateral penumbra dose. Additionally, the V95 from the 90° field is larger than the V95 from the other field angles detailed in Table 3, because the 90° field cannot deform into a re-entrant form of PTV on the beam’s path through. However, there were no significant differences from the 90° field to the 0° and 30° fields, but there was a significant difference between the 90° and 60° fields. Rucinski et al. reported that the V70 and V90 in the 90° field were respectively 12.2±4.7 cm3 and 5.9±2.6 cm3 for carbon beams [18], and Weber et al. reported the V50 Gy in the 90° field was 19.3±3.1% for proton beams [19]. Our results in the 90° field were similar.
Tang et al. reported the rectal volume changes from the prescription dose in 0°, 30°, and 90° fields for proton beams [3]. Kraft and Bassler et al. reported that the lateral penumbra of the carbon beams is sharper than the lateral penumbra of the proton beams, and that the distal tail doses of the carbon beams are higher than the distal tail doses of the proton beams [1, 20]. Using these results, V10, V50, and V95 in each field angle are considered. As compared with the V10 for the carbon beams, V10 in the 90° field for the proton beams is higher than V10 in 0° and 30° fields. The cause is assumed to be that the distal tail of the proton beams is lower than the tail of the carbon beam. Both the V50 and V95 in the 90° field for the proton beams are lower than both the V50 and V95 in the 0° and 30° fields, similar to the fact that both the V50 and V95 in the 90° field for the carbon beams are lower than both the V50 and V95 in the 0° and 30° fields; however, the differences for the proton beams are larger than the differences for the carbon beams. The causes are assumed to be that the lateral penumbra for the proton beams is larger than the penumbra for the carbon beams, and the lateral penumbra for the proton beams in the 90° field beams results in increasing V50 and V95.
The influence of separate setup or range uncertainties on rectal dose
Considering RInc from the setup uncertainty presented in Fig 4, RInc in the AP direction for a 90° field is significantly higher than the ratio for 0°, 30°, or 60°, and RInc in the LR and CC directions is lower than the ratio in the AP direction for 60° and 90° field angles. These findings indicate that the 90° field is disadvantageous for setup uncertainty, and the worst case for the setup uncertainty in the 90° field angle might be concerned only with the setup error in the CC direction. Additionally, RInc on 90° from range uncertainty is significantly lower than the ratio at 0°, 30°, and 60°. This shows that the 90° field is advantageous regarding range uncertainty.
The influence of simultaneous setup and range uncertainties on rectal dose
Considering the influence of both setup and range uncertainties in Fig 4, RInc at 90° shows no statistically significant difference from the ratio for 30° or 60°. However, the 0° field was significantly lower than the ratio for 90°. This suggests that the 0° fields are less affected by uncertainty than fields at different angles. Meanwhile, there are some differences for the forms between DVHs; however, rectal dose increases for the 30° and 60° fields were almost the same as the dose increase for the 90° field. Therefore, the oblique fields can be safely used after checking the dose distribution and DVH. In particular, the V95 regarding the worst case for all angles was similar; however, the V50 was 0°<30°<60°<90°, and the V10 was 0°≈30°>60°>90° (Table 3). Therefore, the 0° field can be used to reduce the middle dose to the rectum, the 90° field can be used to reduce the low dose to the rectum, and the oblique fields can be used to reduce the middle dose and the low dose averagely. Additionally, the sensitivity of the applied irradiation procedure to the setup and range uncertainties is substantially limited because the standard deviations of V10, V50, and V95 were low.
However, in treatment planning, adjusting the rectal dose is achieved as follows: for the 0° field by altering the bolus; for the 90° field by altering the MLC; and for the 30° and 60° fields by altering both. Therefore, planning for the 30° and 60° fields is more complicated than for the 0° and 90° fields.
In the current study, the dose changes from setup and range uncertainties were evaluated using this simple realistic model, and the bladder dose was not evaluated. Although the bladder dose is not a clinical problem, changes in bladder capacity should be noted because they are sensitive to the beam range changes. It is valid not only for the passive irradiation method but also for an active irradiation method. It was useful for clinical treatment; however, we did not factor in dose changes during or between fractions [21–23]. If we assume that position changes regarding prostate intra/inter fraction motion contribute to setup uncertainty, we could use the setup uncertainty for the oblique field and the horizontal field accordingly. However, in the future a more extensive study will be necessary because the results of the current study were limited to a few patient cases and the study was performed using a specific treatment planning procedure.
The influence on CTV dose from simultaneous setup and range uncertainties
There was no influence on dose to the CTV from setup and range uncertainties. We set PTV margins to the CTV in each direction as described in the Materials and Methods section. For example, the posterior margin does not protect against the range of uncertainties when the horizontal field is applied, but it does for the vertical field. In contrast, the posterior margin does not protect from the setup uncertainties when the vertical field is applied, but it does for the horizontal field. Taking into account the various factors, margins in all directions to the CTV are necessary to guarantee sufficient CTV coverage. Additionally, the CTV coverage will be guaranteed when the applied PTV margins widely exceed the shifts in the AP, LR, and CC directions.
The evaluation of oblique fields avoiding the implants
The DVH forms shown in panels (e) and (f) in Fig 5(ii) are similar to the DVH form of averaging 60° shown in panel (c) in Fig 5(i); the DVH form shown in panel (g) in Fig 5(ii) is similar to the DVH form of averaging 30° in Fig 5(B)(i). Additionally, as compared with the V10, V50, and V95 for similar field angles of 30° (35.4±7.8 cm3, 6.7±1.3 cm3, and 2.8±1.1 cm3, respectively) or 60° (26.1±5.2 cm3, 8.3±1.3 cm3, and 2.2±0.8 cm3, respectively), the V10, V50, and V95 for P1, P2, and P3 (16.2 cm3, 5.5 cm3, and 1.4 cm3, respectively for P1; 18.4 cm3, 8.2 cm3, and 2.6 cm3, respectively for P2; and 25.7 cm3, 7.1 cm3, and 3.6 cm3, respectively for P3) were similar or lower, as detailed in Table 3. Furthermore, RInc values from the setup and range uncertainties were similar, as shown in Fig 4. Therefore, this demonstrates that oblique fields avoiding the implant could be safely used in the same way on the non-implant side. In treatment planning, the oblique fields were used to avoid the metal implant while remaining as close to the horizontal as possible. The oblique fields were as good as the horizontal fields in terms of the uncertainties. However, vertical fields might be better than oblique and horizontal fields, as shown in Fig 4.
Conclusion
The influences of setup and range uncertainties on dose deviations in vertical, horizontal, and oblique fields were evaluated in this study. For the basic effect on the rectal dose, it was found that the vertical field could reduce the middle dose to the rectum, the horizontal field could reduce the low dose to the rectum, and the oblique fields could reduce the middle dose and the low dose averagely relative to the other fields. Additionally, the rectal dose deviations from the uncertainties in oblique fields showed no significant difference from those of the horizontal fields; it was found that oblique fields avoiding metal implants could be safely employed because the deviations did not increase with increasingly oblique field angles. The dose to the CTV was preserved over all obliquities.
Because robust optimization methods for the correction of uncertainties have been developed in intensity modulated radiation therapy [24,25], similar methods are needed in particle beam therapy. We hope our findings are the beginning of that process.
Acknowledgments
The authors would like to thank the staff of the GHMC and the Accelerator Engineering Corporation, Chiba, Japan. The authors would also like to thank Dr. Anette Houweling for many helpful discussions.
Author Contributions
Conceived and designed the experiments: YK H. Kawamura. Performed the experiments: YK RT. Analyzed the data: YK MS. Contributed reagents/materials/analysis tools: YK H. Kawamura MS KY NK HS MK H. Katoh. Wrote the paper: YK. Discussed the results and implications: YK H. Kawamura MS MT TK TO TN.
References
- 1. Kraft G. Tumor therapy with heavy charged particles. Prog Part Nucl Phys 2000;45:S473–S544.
- 2. Jäkel O, Reiss P. The influence of metal artefacts on the range of ion beams. Phys Med Biol 2007;52:635–44. pmid:17228110
- 3. Tang S, Both S, Bentefour H, Paly JJ, Tochner Z, Efstathiou J, et al. Improvement of prostate treatment by anterior proton fields. Int J Radiat Oncol Biol Phys 2012;83:408–18. pmid:22133626
- 4. Christodouleas JP, Tang S, Susil RC, McNutt TR, Song DY, Bekelman J, et al. The effect of anterior proton beams in the setting of a prostate-rectum spacer. Med Dosim 2013;38:315–9. pmid:23578497
- 5. Kanematsu N. Quality control of CT system for treatment planning using the polybinary calibration method. Igaku Butsuri 2003;22:140–146 [in Japanese].
- 6. Kanematsu N, Matsufuji N, Kohno R, Minohara S, Kanai T. A CT calibration method based on the polybinary tissue model for radiotherapy treatment planning. Phys Med Biol 2003;48:1053–64. pmid:12741501
- 7. Inaniwa T, Kanematsu N, Furukawa T, Hasegawa A. A robust algorithm of intensity modulated proton therapy for critical tissue sparing and target coverage. Phys Med Biol 2011;56:4749–4770. pmid:21753233
- 8. Liu W, Zhang X, Li Y, Mohan R. Robust optimization of intensity modulated proton therapy. Med Phys 2012;39:1079–1091. pmid:22320818
- 9. Ohno T, Kanai T, Yamada S, Yusa K, Tashiro M, Shimada H, et al. Carbon ion radiotherapy at the Gunma University heavy ion medical center: new facility set-up. Cancers 2011;3:4046–60. pmid:24213124
- 10. Renner TR, Chu WT. Wobbler facility for biomedical experiments. Med Phys 1987;14:825–34. pmid:3683312
- 11. Kanematsu N, Akagi T, Takatani Y, Yonai S, Sakamoto H, Yamashita H. Extended collimator model for pencil-beam dose calculation in proton radiotherapy. Phys Med Biol 2006;51:4807–17. pmid:16985272
- 12. Kanematsu N, Torikoshi M, Mizota M, Kanai T. Secondary range shifting with range compensator for reduction of beam data library in heavy-ion radiotherapy. Med Phys 2007;34:1907–10. pmid:17654892
- 13. Kanematsu N, Yonai S, Ishizaki A. The grid-dose-spreading algorithm for dose distribution calculation in heavy charged particle radiotherapy. Med Phys 2008;35:602–7. pmid:18383681
- 14. Kanematsu N. Dose calculation algorithm of fast fine-heterogeneity correction for heavy charged particle radiotherapy. Physica Med 2011;27:97–102.
- 15. Kanai T, Endo M, Minohara S, Miyahara N, Koyama-ito H, Tomura H, et al. Biophysical characteristics of HIMAC clinical irradiation system for heavy-ion radiation therapy. Int J Radiat Oncol Biol Phys 1999;44:201–10. pmid:10219815
- 16. Kubota Y, Tashiro M, Shinohara A, Abe S, Souda S, Okada R, et al. Development of an automatic evaluation method for patient positioning error. J Appl Clin Med Phys 2015:16:5400–11. pmid:26219004
- 17. Tashiro T, Ishii T, Koya J, Okada R, Kurosawa Y, Arai K, et al. Technical approach to individualized respiratory-gated carbon-ion therapy for mobile organs. Radiol Phys Technol 2013;6:356–366. pmid:23568337
- 18. Rucinski A, Brons S, Richter D, Habl G, Debus J, Bert C, et al. Ion therapy of prostate cancer: daily rectal dose reduction by application of spacer gel. Radiat Oncol 2015;10:1–10.
- 19. Weber DC, Zilli T, Vallee JP, Rouzaud M, Miralbell R, Cozzi L. Intensity modulated proton and photon therapy for early prostate cancer with or without transperineal injection of a polyethylen glycol spacer: a treatment planning comparison study. Int J Radiat Oncol Biol Phys 2012;84:e311–8. pmid:22999271
- 20. Bassler N, Kantemiris I, Karaiskos P, Engelke J, Holzscheiter MH, Petersen JB. Comparison of optimized single and multifield irradiation plans of antiproton, proton and carbon ion beams. Radiother Oncol. 2010;95:87–93. pmid:20304512
- 21. Mah D, Freedman G, Milestone B, Hanlon A, Palacio E, Richardson T, et al. Measurement of intrafractional prostate motion using magnetic resonance imaging. Int J Radiat Oncol Biol Phys 2002;54:568–575. pmid:12243837
- 22. Kurosawa Y, Ishikawa H, Hoshino Y, Higuchi H, Ogano T, Kawamura H, et al. Intra-fractional set-up and organ motion errors in intensity-modulated radiation therapy for prostate cancer. Nihon Hoshasen Gijutsu Gakkai Zasshi 2012;68:290–98 [in Japanese]. pmid:22449906
- 23. Huang K, Palma DA, Scott D, McGregor D, Gaede S, Yartsev S, et al. Inter- and intrafraction uncertainty in prostate bed image-guided radiotherapy. Int J Radiat Oncol Biol Phys 2012;84:402–7. pmid:22381905
- 24. Chu M, Zinchenko Y, Henderson SG, Sharpe MB. Robust optimization for intensity modulated radiation therapy treatment planning under uncertainty. Phys Med Biol 2005;50:5463–5477. pmid:16306645
- 25. Chan TCY, Mišić VV. Adaptive and robust radiation therapy optimization for lung cancer. Eur J Oper Res 2013;231:745–756.