1
|
Rabe M, Paganelli C, Schmitz H, Meschini G, Riboldi M, Hofmaier J, Nierer-Kohlhase L, Dinkel J, Reiner M, Parodi K, Belka C, Landry G, Kurz C, Kamp F. Continuous time-resolved estimated synthetic 4D-CTs for dose reconstruction of lung tumor treatments at a 0.35 T MR-linac. Phys Med Biol 2023; 68:235008. [PMID: 37669669 DOI: 10.1088/1361-6560/acf6f0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Accepted: 09/05/2023] [Indexed: 09/07/2023]
Abstract
Objective.To experimentally validate a method to create continuous time-resolved estimated synthetic 4D-computed tomography datasets (tresCTs) based on orthogonal cine MRI data for lung cancer treatments at a magnetic resonance imaging (MRI) guided linear accelerator (MR-linac).Approach.A breathing porcine lung phantom was scanned at a CT scanner and 0.35 T MR-linac. Orthogonal cine MRI series (sagittal/coronal orientation) at 7.3 Hz, intersecting tumor-mimicking gelatin nodules, were deformably registered to mid-exhale 3D-CT and 3D-MRI datasets. The time-resolved deformation vector fields were extrapolated to 3D and applied to a reference synthetic 3D-CT image (sCTref), while accounting for breathing phase-dependent lung density variations, to create 82 s long tresCTs at 3.65 Hz. Ten tresCTs were created for ten tracked nodules with different motion patterns in two lungs. For each dataset, a treatment plan was created on the mid-exhale phase of a measured ground truth (GT) respiratory-correlated 4D-CT dataset with the tracked nodule as gross tumor volume (GTV). Each plan was recalculated on the GT 4D-CT, randomly sampled tresCT, and static sCTrefimages. Dose distributions for corresponding breathing phases were compared in gamma (2%/2 mm) and dose-volume histogram (DVH) parameter analyses.Main results.The mean gamma pass rate between all tresCT and GT 4D-CT dose distributions was 98.6%. The mean absolute relative deviations of the tresCT with respect to GT DVH parameters were 1.9%, 1.0%, and 1.4% for the GTVD98%,D50%, andD2%, respectively, 1.0% for the remaining nodulesD50%, and 1.5% for the lungV20Gy. The gamma pass rate for the tresCTs was significantly larger (p< 0.01), and the GTVD50%deviations with respect to the GT were significantly smaller (p< 0.01) than for the sCTref.Significance.The results suggest that tresCTs could be valuable for time-resolved reconstruction and intrafractional accumulation of the dose to the GTV for lung cancer patients treated at MR-linacs in the future.
Collapse
Affiliation(s)
- Moritz Rabe
- Department of Radiation Oncology, LMU University Hospital, LMU Munich, Munich, Germany
| | - Chiara Paganelli
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milano, Italy
| | - Henning Schmitz
- Department of Radiation Oncology, LMU University Hospital, LMU Munich, Munich, Germany
| | - Giorgia Meschini
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milano, Italy
| | - Marco Riboldi
- Department of Medical Physics, Ludwig-Maximilians-Universität München (LMU Munich), Garching (Munich), Germany
| | - Jan Hofmaier
- Department of Radiation Oncology, LMU University Hospital, LMU Munich, Munich, Germany
| | - Lukas Nierer-Kohlhase
- Department of Radiation Oncology, LMU University Hospital, LMU Munich, Munich, Germany
| | - Julien Dinkel
- Department of Radiology, LMU University Hospital, LMU Munich, Munich, Germany
- Comprehensive Pneumology Center, German Center for Lung Research (DZL), Munich, Germany
| | - Michael Reiner
- Department of Radiation Oncology, LMU University Hospital, LMU Munich, Munich, Germany
| | - Katia Parodi
- Department of Medical Physics, Ludwig-Maximilians-Universität München (LMU Munich), Garching (Munich), Germany
| | - Claus Belka
- Department of Radiation Oncology, LMU University Hospital, LMU Munich, Munich, Germany
- German Cancer Consortium (DKTK), partner site Munich, a partnership between DKFZ and LMU University Hospital Munich, Germany
- Bavarian Cancer Research Center (BZKF), Munich, Germany
| | - Guillaume Landry
- Department of Radiation Oncology, LMU University Hospital, LMU Munich, Munich, Germany
| | - Christopher Kurz
- Department of Radiation Oncology, LMU University Hospital, LMU Munich, Munich, Germany
| | - Florian Kamp
- Department of Radiation Oncology, LMU University Hospital, LMU Munich, Munich, Germany
- Department of Radiation Oncology, University Hospital Cologne, Cologne, Germany
| |
Collapse
|
2
|
Parrella G, Vai A, Nakas A, Garau N, Meschini G, Camagni F, Molinelli S, Barcellini A, Pella A, Ciocca M, Vitolo V, Orlandi E, Paganelli C, Baroni G. Synthetic CT in Carbon Ion Radiotherapy of the Abdominal Site. Bioengineering (Basel) 2023; 10:bioengineering10020250. [PMID: 36829745 PMCID: PMC9951997 DOI: 10.3390/bioengineering10020250] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2023] [Revised: 02/06/2023] [Accepted: 02/09/2023] [Indexed: 02/17/2023] Open
Abstract
The generation of synthetic CT for carbon ion radiotherapy (CIRT) applications is challenging, since high accuracy is required in treatment planning and delivery, especially in an anatomical site as complex as the abdomen. Thirty-nine abdominal MRI-CT volume pairs were collected and a three-channel cGAN (accounting for air, bones, soft tissues) was used to generate sCTs. The network was tested on five held-out MRI volumes for two scenarios: (i) a CT-based segmentation of the MRI channels, to assess the quality of sCTs and (ii) an MRI manual segmentation, to simulate an MRI-only treatment scenario. The sCTs were evaluated by means of similarity metrics (e.g., mean absolute error, MAE) and geometrical criteria (e.g., dice coefficient). Recalculated CIRT plans were evaluated through dose volume histogram, gamma analysis and range shift analysis. The CT-based test set presented optimal MAE on bones (86.03 ± 10.76 HU), soft tissues (55.39 ± 3.41 HU) and air (54.42 ± 11.48 HU). Higher values were obtained from the MRI-only test set (MAEBONE = 154.87 ± 22.90 HU). The global gamma pass rate reached 94.88 ± 4.9% with 3%/3 mm, while the range shift reached a median (IQR) of 0.98 (3.64) mm. The three-channel cGAN can generate acceptable abdominal sCTs and allow for CIRT dose recalculations comparable to the clinical plans.
Collapse
Affiliation(s)
- Giovanni Parrella
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milan, Italy
- Correspondence: ; Tel.: +39-02-2399-18-9022
| | - Alessandro Vai
- Medical Physics Unit, National Center of Oncological Hadrontherapy (CNAO), Strada Campeggi, 53, 27100 Pavia, Italy
| | - Anestis Nakas
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milan, Italy
| | - Noemi Garau
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milan, Italy
| | - Giorgia Meschini
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milan, Italy
| | - Francesca Camagni
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milan, Italy
| | - Silvia Molinelli
- Medical Physics Unit, National Center of Oncological Hadrontherapy (CNAO), Strada Campeggi, 53, 27100 Pavia, Italy
| | - Amelia Barcellini
- Radiotherapy Unit, National Center of Oncological Hadrontherapy (CNAO), Strada Campeggi, 53, 27100 Pavia, Italy
- Department of Internal Medicine and Medical Therapy, University of Pavia, 27100 Pavia, Italy
| | - Andrea Pella
- Bioengineering Unit, National Center of Oncological Hadrontherapy (CNAO), Strada Campeggi, 53, 27100 Pavia, Italy
| | - Mario Ciocca
- Medical Physics Unit, National Center of Oncological Hadrontherapy (CNAO), Strada Campeggi, 53, 27100 Pavia, Italy
| | - Viviana Vitolo
- Radiotherapy Unit, National Center of Oncological Hadrontherapy (CNAO), Strada Campeggi, 53, 27100 Pavia, Italy
| | - Ester Orlandi
- Clinical Unit, National Center of Oncological Hadrontherapy (CNAO), Strada Campeggi, 53, 27100 Pavia, Italy
| | - Chiara Paganelli
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milan, Italy
| | - Guido Baroni
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milan, Italy
| |
Collapse
|
3
|
Duetschler A, Prendi J, Safai S, Weber DC, Lomax AJ, Zhang Y. Limitations of phase-sorting based pencil beam scanned 4D proton dose calculations under irregular motion. Phys Med Biol 2022; 68. [PMID: 36571234 DOI: 10.1088/1361-6560/aca9b6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2022] [Accepted: 12/07/2022] [Indexed: 12/12/2022]
Abstract
Objective.4D dose calculation (4DDC) for pencil beam scanned (PBS) proton therapy is typically based on phase-sorting of individual pencil beams onto phases of a single breathing cycle 4DCT. Understanding the dosimetric limitations and uncertainties of this approach is essential, especially for the realistic treatment scenario with irregular free breathing motion.Approach.For three liver and three lung cancer patient CTs, the deformable multi-cycle motion from 4DMRIs was used to generate six synthetic 4DCT(MRI)s, providing irregular motion (11/15 cycles for liver/lung; tumor amplitudes ∼4-18 mm). 4DDCs for two-field plans were performed, with the temporal resolution of the pencil beam delivery (4-200 ms) or with 8 phases per breathing cycle (500-1000 ms). For the phase-sorting approach, the tumor center motion was used to determine the phase assignment of each spot. The dose was calculated either using the full free breathing motion or individually repeating each single cycle. Additionally, the use of an irregular surrogate signal prior to 4DDC on a repeated cycle was simulated. The CTV volume with absolute dose differences >5% (Vdosediff>5%) and differences in CTVV95%andD5%-D95%compared to the free breathing scenario were evaluated.Main results.Compared to 4DDC considering the full free breathing motion with finer spot-wise temporal resolution, 4DDC based on a repeated single 4DCT resulted inVdosediff>5%of on average 34%, which resulted in an overestimation ofV95%up to 24%. However, surrogate based phase-sorting prior to 4DDC on a single cycle 4DCT, reduced the averageVdosediff>5%to 16% (overestimationV95%up to 19%). The 4DDC results were greatly influenced by the choice of reference cycle (Vdosediff>5%up to 55%) and differences due to temporal resolution were much smaller (Vdosediff>5%up to 10%).Significance.It is important to properly consider motion irregularity in 4D dosimetric evaluations of PBS proton treatments, as 4DDC based on a single 4DCT can lead to an underestimation of motion effects.
Collapse
Affiliation(s)
- A Duetschler
- Center for Proton Therapy, Paul Scherrer Institute, 5232 Villigen PSI, CH, Switzerland.,Department of Physics, ETH Zürich, 8092 Zürich, CH, Switzerland
| | - J Prendi
- Center for Proton Therapy, Paul Scherrer Institute, 5232 Villigen PSI, CH, Switzerland.,Department of Physics, University of Basel, 4056 Basel, CH, Switzerland
| | - S Safai
- Center for Proton Therapy, Paul Scherrer Institute, 5232 Villigen PSI, CH, Switzerland
| | - D C Weber
- Center for Proton Therapy, Paul Scherrer Institute, 5232 Villigen PSI, CH, Switzerland.,Department of Radiation Oncology, University Hospital of Zürich, 8091 Zürich, CH, Switzerland.,Department of Radiation Oncology, Inselspital, Bern University Hospital, University of Bern, 3010 Bern, CH, Switzerland
| | - A J Lomax
- Center for Proton Therapy, Paul Scherrer Institute, 5232 Villigen PSI, CH, Switzerland.,Department of Physics, ETH Zürich, 8092 Zürich, CH, Switzerland
| | - Ye Zhang
- Center for Proton Therapy, Paul Scherrer Institute, 5232 Villigen PSI, CH, Switzerland
| |
Collapse
|
4
|
Molinelli S, Vai A, Russo S, Loap P, Meschini G, Paganelli C, Barcellini A, Vitolo V, Orlandi E, Ciocca M. The role of multiple anatomical scenarios in plan optimization for carbon ion radiotherapy of pancreatic cancer. Radiother Oncol 2022; 176:1-8. [PMID: 36113776 DOI: 10.1016/j.radonc.2022.09.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Revised: 08/17/2022] [Accepted: 09/07/2022] [Indexed: 12/14/2022]
Abstract
PURPOSE /OBJECTIVE To quantify benefits of robust optimization on multiple 4DCT acquisitions combined with off-line treatment adaptation for neoadjuvant carbon ion therapy (CIRT) of pancreatic cancer. MATERIAL/METHODS For 10 previously treated patients, 4DCTs were acquired around -15 (CTPlan), -5 (RE1), -1 (RE2) and +6 (RE3) days from RT start. Treatment plans were newly optimized to a dose prescription of 38.4 Gy(RBE) (8 fractions) with a constraint of 38 Gy(RBE) to 1% of the gastrointestinal organs at risk volume (D1%). Three strategies were tested: (A) robust optimization on CTPlan maximum exhale (0Ex) with 3 mm set-up, 3% range uncertainty, including 30%-inhale; (B) addition of the RE1-0Ex scenario; (C) plan recalculation at each REi and adaptation (RPi) according to deviation thresholds from clinical goals. The cumulative variation of target coverage and GI-OARs doses was evaluated. Duodenum contours of all 4DCTs of each patient were registered on CTPlan-0Ex. The capacity of pre-RT acquisitions to predict duodenum position was investigated by computing the intersection of contours at CTplan, RE1, or their union, with respect to subsequent 4DCTs and the CTV, coupled with increasing margin. RESULTS (A) No recalculation exceeded the D1% constraint. (B) The inclusion of RE1-0Ex in the optimization problem improved inter-fraction robustness on a patient-specific basis, but was non-significant on average. (C) Half of the plans would be re-optimized to recover target coverage and/or minimize duodenum dose, at least once. A significant difference was observed between pre-RT duodenum contours when intersecting subsequent contours, either with a margin expansion. CONCLUSION Anatomical variations highlighted at multiple REi proved that a fast and efficient online adaptation is essential to optimize treatment quality of CIRT for pancreatic cancer.
Collapse
Affiliation(s)
| | | | | | - Pierre Loap
- Institut Curie, Department of Radiation Oncology, Paris, France
| | - Giorgia Meschini
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milano, Italy
| | - Chiara Paganelli
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milano, Italy
| | | | | | | | - Mario Ciocca
- Dipartimento Clinico, Fondazione CNAO, Pavia, Italy
| |
Collapse
|
5
|
Bäumer C, Frakulli R, Kohl J, Nagaraja S, Steinmeier T, Worawongsakul R, Timmermann B. Adaptive Proton Therapy of Pediatric Head and Neck Cases Using MRI-Based Synthetic CTs: Initial Experience of the Prospective KiAPT Study. Cancers (Basel) 2022; 14:cancers14112616. [PMID: 35681594 PMCID: PMC9179385 DOI: 10.3390/cancers14112616] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2022] [Revised: 05/20/2022] [Accepted: 05/20/2022] [Indexed: 12/07/2022] Open
Abstract
BACKGROUND AND PURPOSE Interfractional anatomical changes might affect the outcome of proton therapy (PT). We aimed to prospectively evaluate the role of Magnetic Resonance Imaging (MRI) based adaptive PT for children with tumors of the head and neck and base of skull. METHODS MRI verification images were acquired at half of the treatment course. A synthetic computed tomography (CT) image was created using this MRI and a deformable image registration (DIR) to the reference MRI. The methodology was verified with in-silico phantoms and validated using a clinical case with a shrinking cystic hygroma on the basis of dosimetric quantities of contoured structures. The dose distributions on the verification X-ray CT and on the synthetic CT were compared with a gamma-index test using global 2 mm/2% criteria. RESULTS Regarding the clinical validation case, the gamma-index pass rate was 98.3%. Eleven patients were included in the clinical study. The most common diagnosis was rhabdomyosarcoma (73%). Craniofacial tumor site was predominant in 64% of patients, followed by base of skull (18%). For one individual case the synthetic CT showed an increase in the median D2 and Dmax dose on the spinal cord from 20.5 GyRBE to 24.8 GyRBE and 14.7 GyRBE to 25.1 GyRBE, respectively. Otherwise, doses received by OARs remained relatively stable. Similarly, the target volume coverage seen by D95% and V95% remained unchanged. CONCLUSIONS The method of transferring anatomical changes from MRIs to a synthetic CTs was successfully implemented and validated with simple, commonly available tools. In the frame of our early results on a small cohort, no clinical relevant deterioration for neither PTV coverage nor an increased dose burden to OARs occurred. However, the study will be continued to identify a pediatric patient cohort, which benefits from adaptive treatment planning.
Collapse
Affiliation(s)
- Christian Bäumer
- West German Proton Therapy Centre Essen, 45147 Essen, Germany; (R.F.); (J.K.); (S.N.); (T.S.); (R.W.); (B.T.)
- University Hospital Essen, 45147 Essen, Germany
- West German Cancer Center (WTZ), 45147 Essen, Germany
- German Cancer Consortium (DKTK), 69120 Heidelberg, Germany
- Department of Physics, Technische Universität Dortmund, 44227 Dortmund, Germany
- Correspondence:
| | - Rezarta Frakulli
- West German Proton Therapy Centre Essen, 45147 Essen, Germany; (R.F.); (J.K.); (S.N.); (T.S.); (R.W.); (B.T.)
- University Hospital Essen, 45147 Essen, Germany
- West German Cancer Center (WTZ), 45147 Essen, Germany
- Department of Particle Therapy, 45147 Essen, Germany
| | - Jessica Kohl
- West German Proton Therapy Centre Essen, 45147 Essen, Germany; (R.F.); (J.K.); (S.N.); (T.S.); (R.W.); (B.T.)
- University Hospital Essen, 45147 Essen, Germany
- West German Cancer Center (WTZ), 45147 Essen, Germany
| | - Sindhu Nagaraja
- West German Proton Therapy Centre Essen, 45147 Essen, Germany; (R.F.); (J.K.); (S.N.); (T.S.); (R.W.); (B.T.)
- University Hospital Essen, 45147 Essen, Germany
- West German Cancer Center (WTZ), 45147 Essen, Germany
- Department of Particle Therapy, 45147 Essen, Germany
| | - Theresa Steinmeier
- West German Proton Therapy Centre Essen, 45147 Essen, Germany; (R.F.); (J.K.); (S.N.); (T.S.); (R.W.); (B.T.)
- University Hospital Essen, 45147 Essen, Germany
- West German Cancer Center (WTZ), 45147 Essen, Germany
- Department of Particle Therapy, 45147 Essen, Germany
| | - Rasin Worawongsakul
- West German Proton Therapy Centre Essen, 45147 Essen, Germany; (R.F.); (J.K.); (S.N.); (T.S.); (R.W.); (B.T.)
- University Hospital Essen, 45147 Essen, Germany
- West German Cancer Center (WTZ), 45147 Essen, Germany
- Department of Particle Therapy, 45147 Essen, Germany
- Radiation Oncology Unit, Department of Diagnostic and Therapeutic Radiology, Ramathibodi Hospital, Mahidol University, Nakhon 73170, Thailand
| | - Beate Timmermann
- West German Proton Therapy Centre Essen, 45147 Essen, Germany; (R.F.); (J.K.); (S.N.); (T.S.); (R.W.); (B.T.)
- University Hospital Essen, 45147 Essen, Germany
- West German Cancer Center (WTZ), 45147 Essen, Germany
- German Cancer Consortium (DKTK), 69120 Heidelberg, Germany
- Department of Particle Therapy, 45147 Essen, Germany
| |
Collapse
|
6
|
Meschini G, Vai A, Barcellini A, Fontana G, Molinelli S, Mastella E, Pella A, Vitolo V, Imparato S, Orlandi E, Ciocca M, Baroni G, Paganelli C. Time-resolved MRI for off-line treatment robustness evaluation in carbon-ion radiotherapy of pancreatic cancer. Med Phys 2022; 49:2386-2395. [PMID: 35124811 PMCID: PMC9306947 DOI: 10.1002/mp.15510] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2021] [Revised: 09/28/2021] [Accepted: 01/20/2022] [Indexed: 11/25/2022] Open
Abstract
PURPOSE In this study, we investigate the use of magnetic resonance imaging (MRI) for the clinical evaluation of gating treatment robustness in carbon-ion radiotherapy (CIRT) of pancreatic cancer. Indeed, MRI allows radiation-free repeated scans and fast dynamic sequences for time-resolved (TR) imaging (cine-MRI), providing information on inter- and intra-fraction cycle-to-cycle variations of respiratory motion. MRI can therefore support treatment planning and verification, overcoming the limitations of the current clinical standard, that is, four-dimensional computed tomography (4DCT), which describes an "average" breathing cycle neglecting breathing motion variability. METHODS We integrated a technique to generate a virtual CT (vCT) from 3D MRI with a method for 3D reconstruction from 2D cine-MRI, to produce TR vCTs for dose recalculations. For eight patients, the method allowed evaluating inter-fraction variations at end-exhale and intra-fraction cycle-to-cycle variability within the gating window in terms of tumor displacement and dose to the target and organs at risk. RESULTS The median inter-fraction tumor motion was in the range 3.33-12.16 mm, but the target coverage was robust (-0.4% median D95% variation). Concerning cycle-to-cycle variations, the gating technique was effective in limiting tumor displacement (1.35 mm median gating motion) and corresponding dose variations (-3.9% median D95% variation). The larger exposure of organs at risk (duodenum and stomach) was caused by inter-fraction motion, whereas intra-fraction cycle-to-cycle dose variations were limited. CONCLUSIONS This study proposed a method for the generation of TR vCTs from MRI, which enabled an off-line evaluation of gating treatment robustness and suggested its feasibility to support treatment planning of pancreatic tumors in CIRT.
Collapse
Affiliation(s)
- Giorgia Meschini
- Department of Electronics, Information and BioengineeringPolitecnico di MilanoMilanItaly
| | - Alessandro Vai
- Medical Physics UnitNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Amelia Barcellini
- Clinical DepartmentNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Giulia Fontana
- Clinical Bioengineering UnitNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Silvia Molinelli
- Medical Physics UnitNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Edoardo Mastella
- Medical Physics UnitNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Andrea Pella
- Clinical Bioengineering UnitNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Viviana Vitolo
- Clinical DepartmentNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Sara Imparato
- Radiology UnitNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Ester Orlandi
- Clinical DepartmentNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Mario Ciocca
- Medical Physics UnitNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Guido Baroni
- Department of Electronics, Information and BioengineeringPolitecnico di MilanoMilanItaly
- Clinical Bioengineering UnitNational Center for Oncological Hadrontherapy (Fondazione CNAO)PaviaItaly
| | - Chiara Paganelli
- Department of Electronics, Information and BioengineeringPolitecnico di MilanoMilanItaly
| |
Collapse
|
7
|
Duetschler A, Bauman G, Bieri O, Cattin PC, Ehrbar S, Engin-Deniz G, Giger A, Josipovic M, Jud C, Krieger M, Nguyen D, Persson GF, Salomir R, Weber DC, Lomax AJ, Zhang Y. Synthetic 4DCT(MRI) lung phantom generation for 4D radiotherapy and image guidance investigations. Med Phys 2022; 49:2890-2903. [PMID: 35239984 PMCID: PMC9313613 DOI: 10.1002/mp.15591] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2021] [Revised: 12/26/2021] [Accepted: 02/24/2022] [Indexed: 11/14/2022] Open
Abstract
Purpose Respiratory motion is one of the major challenges in radiotherapy. In this work, a comprehensive and clinically plausible set of 4D numerical phantoms, together with their corresponding “ground truths,” have been developed and validated for 4D radiotherapy applications. Methods The phantoms are based on CTs providing density information and motion from multi‐breathing‐cycle 4D Magnetic Resonance imagings (MRIs). Deformable image registration (DIR) has been utilized to extract motion fields from 4DMRIs and to establish inter‐subject correspondence by registering binary lung masks between Computer Tomography (CT) and MRI. The established correspondence is then used to warp the CT according to the 4DMRI motion. The resulting synthetic 4DCTs are called 4DCT(MRI)s. Validation of the 4DCT(MRI) workflow was conducted by directly comparing conventional 4DCTs to derived synthetic 4D images using the motion of the 4DCTs themselves (referred to as 4DCT(CT)s). Digitally reconstructed radiographs (DRRs) as well as 4D pencil beam scanned (PBS) proton dose calculations were used for validation. Results Based on the CT image appearance of 13 lung cancer patients and deformable motion of five volunteer 4DMRIs, synthetic 4DCT(MRI)s with a total of 871 different breathing cycles have been generated. The 4DCT(MRI)s exhibit an average superior–inferior tumor motion amplitude of 7 ± 5 mm (min: 0.5 mm, max: 22.7 mm). The relative change of the DRR image intensities of the conventional 4DCTs and the corresponding synthetic 4DCT(CT)s inside the body is smaller than 5% for at least 81% of the pixels for all studied cases. Comparison of 4D dose distributions calculated on 4DCTs and the synthetic 4DCT(CT)s using the same motion achieved similar dose distributions with an average 2%/2 mm gamma pass rate of 90.8% (min: 77.8%, max: 97.2%). Conclusion We developed a series of numerical 4D lung phantoms based on real imaging and motion data, which give realistic representations of both anatomy and motion scenarios and the accessible “ground truth” deformation vector fields of each 4DCT(MRI). The open‐source code and motion data allow foreseen users to generate further 4D data by themselves. These numeric 4D phantoms can be used for the development of new 4D treatment strategies, 4D dose calculations, DIR algorithm validations, as well as simulations of motion mitigation and different online image guidance techniques for both proton and photon radiation therapy.
Collapse
Affiliation(s)
- Alisha Duetschler
- Center for Proton Therapy, Paul Scherrer Institute, Villigen PSI, 5232, Switzerland.,Department of Physics, ETH Zurich, Zurich, 8092, Switzerland
| | - Grzegorz Bauman
- Department of Biomedical Engineering, University of Basel, Allschwil, 4123, Switzerland.,Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, 4031, Switzerland
| | - Oliver Bieri
- Department of Biomedical Engineering, University of Basel, Allschwil, 4123, Switzerland.,Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, 4031, Switzerland
| | - Philippe C Cattin
- Department of Biomedical Engineering, University of Basel, Allschwil, 4123, Switzerland.,Center for medical Image Analysis & Navigation, University of Basel, Allschwil, 4123, Switzerland
| | - Stefanie Ehrbar
- Department of Radiation Oncology, University Hospital of Zurich, Zurich, 8091, Switzerland.,University of Zurich, Zurich, 8006, Switzerland
| | - Georg Engin-Deniz
- Center for Proton Therapy, Paul Scherrer Institute, Villigen PSI, 5232, Switzerland.,Department of Physics, ETH Zurich, Zurich, 8092, Switzerland
| | - Alina Giger
- Department of Biomedical Engineering, University of Basel, Allschwil, 4123, Switzerland.,Center for medical Image Analysis & Navigation, University of Basel, Allschwil, 4123, Switzerland
| | - Mirjana Josipovic
- Department of Oncology, Rigshospitalet Copenhagen University Hospital, Copenhagen, 2100, Denmark
| | - Christoph Jud
- Department of Biomedical Engineering, University of Basel, Allschwil, 4123, Switzerland.,Center for medical Image Analysis & Navigation, University of Basel, Allschwil, 4123, Switzerland
| | - Miriam Krieger
- Center for Proton Therapy, Paul Scherrer Institute, Villigen PSI, 5232, Switzerland.,Department of Physics, ETH Zurich, Zurich, 8092, Switzerland
| | - Damien Nguyen
- Department of Biomedical Engineering, University of Basel, Allschwil, 4123, Switzerland.,Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, 4031, Switzerland
| | - Gitte F Persson
- Department of Oncology, Rigshospitalet Copenhagen University Hospital, Copenhagen, 2100, Denmark.,Department of Oncology, Herlev-Gentofte Hospital Copenhagen University Hospital, Herlev, 2730, Denmark.,Department of Clinical Medicine, Faculty of Medical Sciences, University of Copenhagen, Copenhagen, 2100, Denmark
| | - Rares Salomir
- Image Guided Interventions Laboratory (949), Faculty of Medicine, University of Geneva, Geneva, 1211, Switzerland.,Radiology Division, University Hospitals of Geneva, Geneva, 1205, Switzerland
| | - Damien C Weber
- Center for Proton Therapy, Paul Scherrer Institute, Villigen PSI, 5232, Switzerland.,Department of Radiation Oncology, University Hospital of Zurich, Zurich, 8091, Switzerland.,Department of Radiation Oncology, University of Bern, Bern, 3010, Switzerland
| | - Antony J Lomax
- Center for Proton Therapy, Paul Scherrer Institute, Villigen PSI, 5232, Switzerland.,Department of Physics, ETH Zurich, Zurich, 8092, Switzerland
| | - Ye Zhang
- Center for Proton Therapy, Paul Scherrer Institute, Villigen PSI, 5232, Switzerland
| |
Collapse
|
8
|
Abstract
Protons and carbon ions (hadrons) have useful properties for the treatments of patients affected by oncological pathologies. They are more precise than conventional X-rays and possess radiobiological characteristics suited for treating radio-resistant or inoperable tumours. This paper gives an overview of the status of hadron therapy around the world. It focusses on the Italian National Centre for Oncological Hadron therapy (CNAO), introducing operation procedures, system performance, expansion projects, methodologies and modelling to build individualized treatments. There is growing evidence that supports safety and effectiveness of hadron therapy for a variety of clinical situations. However, there is still a lack of high-level evidence directly comparing hadron therapy with modern conventional radiotherapy techniques. The results give an overview of pre-clinical and clinical research studies and of the treatments of 3700 patients performed at CNAO. The success and development of hadron therapy is strongly associated with the creation of networks among hadron therapy facilities, clinics, universities and research institutions. These networks guarantee the growth of cultural knowledge on hadron therapy, favour the efficient recruitment of patients and present available competences for R&D (Research and Development) programmes.
Collapse
|
9
|
Extension of RBE-weighted 4D particle dose calculation for non-periodic motion. Phys Med 2021; 91:62-72. [PMID: 34715550 DOI: 10.1016/j.ejmp.2021.10.009] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/08/2021] [Revised: 10/01/2021] [Accepted: 10/05/2021] [Indexed: 01/02/2023] Open
Abstract
PURPOSE Highly conformal scanned Carbon Ion Radiotherapy (CIRT) might permit dose escalation and improved local control in advanced stage thoracic tumors, but is challenged by target motion. Dose calculation algorithms typically assume a periodically repeating, regular motion. To assess the effect of realistic, irregular motion, new algorithms of validated accuracy are needed. METHODS We extended an in-house treatment planning system to calculate RBE-weighted dose distributions in CIRT on non-periodic CT image sequences. Dosimetric accuracy was validated experimentally on a moving, time-resolved ionization chamber array. Log-file based dose reconstructions were compared by gamma analysis and correlation to measurements at every intermediate detector frame during delivery. The impact of irregular motion on treatment quality was simulated on a virtual 4DCT thorax phantom. Periodic motion was compared to motion with varying amplitude and period ± baseline drift. Rescanning as a mitigation strategy was assessed on all scenarios. RESULTS In experimental validation, average gamma pass rates were 99.89+-0.30% for 3%/3 mm and 88.2+-2.2% for 2%/2 mm criteria. Average correlation for integral dose distributions was 0.990±0.002. Median correlation for single 200 ms frames was 0.947±0.006. In the simulations, irregular motion deteriorated V95 target coverage to 81.2%, 76.6% and 79.0% for regular, irregular motion and irregular motion with base-line drift, respectively. Rescanning restored V95 to >98% for both scenarios without baseline drift, but not with additional baseline drift at 83.7%. CONCLUSIONS The validated algorithm permits to study the effects of irregular motion and to develop and adapt appropriate motion mitigation techniques.
Collapse
|
10
|
Target motion management in breast cancer radiation therapy. Radiol Oncol 2021; 55:393-408. [PMID: 34626533 PMCID: PMC8647788 DOI: 10.2478/raon-2021-0040] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Accepted: 08/04/2021] [Indexed: 12/25/2022] Open
Abstract
Background Over the last two decades, breast cancer remains the main cause of cancer deaths in women. To treat this type of cancer, radiation therapy (RT) has proved to be efficient. RT for breast cancer is, however, challenged by intrafractional motion caused by respiration. The problem is more severe for the left-sided breast cancer due to the proximity to the heart as an organ-at-risk. While particle therapy results in superior dose characteristics than conventional RT, due to the physics of particle interactions in the body, particle therapy is more sensitive to target motion. Conclusions This review highlights current and emerging strategies for the management of intrafractional target motion in breast cancer treatment with an emphasis on particle therapy, as a modern RT technique. There are major challenges associated with transferring real-time motion monitoring technologies from photon to particles beams. Surface imaging would be the dominant imaging modality for real-time intrafractional motion monitoring for breast cancer. The magnetic resonance imaging (MRI) guidance and ultra high dose rate (FLASH)-RT seem to be state-of-the-art approaches to deal with 4D RT for breast cancer.
Collapse
|
11
|
Mastella E, Mirandola A, Russo S, Vai A, Magro G, Molinelli S, Barcellini A, Vitolo V, Orlandi E, Ciocca M. High-dose hypofractionated pencil beam scanning carbon ion radiotherapy for lung tumors: Dosimetric impact of different spot sizes and robustness to interfractional uncertainties. Phys Med 2021; 85:79-86. [PMID: 33984821 DOI: 10.1016/j.ejmp.2021.05.004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Revised: 03/18/2021] [Accepted: 05/03/2021] [Indexed: 02/08/2023] Open
Abstract
PURPOSE The robustness against setup and motion uncertainties of gated four-dimensional restricted robust optimization (4DRRO) was investigated for hypofractionated carbon ion radiotherapy (CIRT) of lung tumors. METHODS CIRT plans of 9 patients were optimized using 4DRRO strategy with 3 mm setup errors, 3% density errors and 3 breathing phases related to the gate window. The prescription was 60 Gy(RBE) in 4 fractions. Standard spots (SS) were compared to big spots (BS). Plans were recalculated on multiple 4DCTs acquired within 3 weeks from treatment simulation and rigidly registered with planning images using bone matching. Warped dose distributions were generated using deformable image registration and accumulated on the planning 4DCTs. Target coverage (D98%, D95% and V95%) and dose to lung were evaluated in the recalculated and accumulated dose distributions. RESULTS Comparable target coverage was obtained with both spot sizes (p = 0.53 for D95%). The mean lung dose increased of 0.6 Gy(RBE) with BS (p = 0.0078), still respecting the dose constraint of a 4-fraction stereotactic treatment for the risk of radiation pneumonitis. Statistically significant differences were found in the recalculated and accumulated D95% (p = 0.048 and p = 0.024), with BS showing to be more robust. Using BS, the average degradations of the D98%, D95% and V95% in the accumulated doses were -2.7%, -1.6% and -1.5%. CONCLUSIONS Gated 4DRRO was highly robust against setup and motion uncertainties. BS increased the dose to healthy tissues but were more robust than SS. The selected optimization settings guaranteed adequate target coverage during the simulated treatment course with acceptable risk of toxicity.
Collapse
Affiliation(s)
- Edoardo Mastella
- CNAO, National Center for Oncological Hadrontherapy, Strada Campeggi 53, I-27100 Pavia, Italy.
| | - Alfredo Mirandola
- CNAO, National Center for Oncological Hadrontherapy, Strada Campeggi 53, I-27100 Pavia, Italy
| | - Stefania Russo
- CNAO, National Center for Oncological Hadrontherapy, Strada Campeggi 53, I-27100 Pavia, Italy
| | - Alessandro Vai
- CNAO, National Center for Oncological Hadrontherapy, Strada Campeggi 53, I-27100 Pavia, Italy
| | - Giuseppe Magro
- CNAO, National Center for Oncological Hadrontherapy, Strada Campeggi 53, I-27100 Pavia, Italy
| | - Silvia Molinelli
- CNAO, National Center for Oncological Hadrontherapy, Strada Campeggi 53, I-27100 Pavia, Italy
| | - Amelia Barcellini
- CNAO, National Center for Oncological Hadrontherapy, Strada Campeggi 53, I-27100 Pavia, Italy
| | - Viviana Vitolo
- CNAO, National Center for Oncological Hadrontherapy, Strada Campeggi 53, I-27100 Pavia, Italy
| | - Ester Orlandi
- CNAO, National Center for Oncological Hadrontherapy, Strada Campeggi 53, I-27100 Pavia, Italy
| | - Mario Ciocca
- CNAO, National Center for Oncological Hadrontherapy, Strada Campeggi 53, I-27100 Pavia, Italy
| |
Collapse
|
12
|
Rabe M, Paganelli C, Riboldi M, Bondesson D, Jörg Schneider M, Chmielewski T, Baroni G, Dinkel J, Reiner M, Landry G, Parodi K, Belka C, Kamp F, Kurz C. Porcine lung phantom-based validation of estimated 4D-MRI using orthogonal cine imaging for low-field MR-Linacs. Phys Med Biol 2021; 66:055006. [PMID: 33171458 DOI: 10.1088/1361-6560/abc937] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Real-time motion monitoring of lung tumors with low-field magnetic resonance imaging-guided linear accelerators (MR-Linacs) is currently limited to sagittal 2D cine magnetic resonance imaging (MRI). To provide input data for improved intrafractional and interfractional adaptive radiotherapy, the 4D anatomy has to be inferred from data with lower dimensionality. The purpose of this study was to experimentally validate a previously proposed propagation method that provides continuous time-resolved estimated 4D-MRI based on orthogonal cine MRI for a low-field MR-Linac. Ex vivo porcine lungs were injected with artificial nodules and mounted in a dedicated phantom that allows for the simulation of periodic and reproducible breathing motion. The phantom was scanned with a research version of a commercial 0.35 T MR-Linac. Respiratory-correlated 4D-MRI were reconstructed and served as ground truth images. Series of interleaved orthogonal slices in sagittal and coronal orientation, intersecting the injected targets, were acquired at 7.3 Hz. Estimated 4D-MRI at 3.65 Hz were created in post-processing using the propagation method and compared to the ground truth 4D-MRI. Eight datasets at different breathing frequencies and motion amplitudes were acquired for three porcine lungs. The overall median (95[Formula: see text] percentile) deviation between ground truth and estimated deformation vector fields was 2.3 mm (5.7 mm), corresponding to 0.7 (1.6) times the in-plane imaging resolution (3.5 × 3.5 mm2). Median (95[Formula: see text] percentile) estimated nodule position errors were 1.5 mm (3.8 mm) for nodules intersected by orthogonal slices and 2.1 mm (7.1 mm) for nodules located more than 2 cm away from either of the orthogonal slices. The estimation error depended on the breathing phase, the motion amplitude and the location of the estimated position with respect to the orthogonal slices. By using the propagation method, the 4D motion within the porcine lung phantom could be accurately and robustly estimated. The method could provide valuable information for treatment planning, real-time motion monitoring, treatment adaptation, and post-treatment evaluation of MR-guided radiotherapy treatments.
Collapse
Affiliation(s)
- Moritz Rabe
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
| | - Chiara Paganelli
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milano, Italy
| | - Marco Riboldi
- Department of Medical Physics, Ludwig-Maximilians-Universität München (LMU Munich), Garching (Munich), Germany
| | - David Bondesson
- Department of Radiology, University Hospital, LMU Munich, Munich, Germany.,Comprehensive Pneumology Center, German Center for Lung Research (DZL), Munich, Germany
| | - Moritz Jörg Schneider
- Department of Radiology, University Hospital, LMU Munich, Munich, Germany.,Comprehensive Pneumology Center, German Center for Lung Research (DZL), Munich, Germany
| | | | - Guido Baroni
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milano, Italy.,Bioengineering Unit, Centro Nazionale di Adroterapia Oncologica, Pavia, Italy
| | - Julien Dinkel
- Department of Radiology, University Hospital, LMU Munich, Munich, Germany.,Comprehensive Pneumology Center, German Center for Lung Research (DZL), Munich, Germany
| | - Michael Reiner
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
| | - Guillaume Landry
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany.,Department of Medical Physics, Ludwig-Maximilians-Universität München (LMU Munich), Garching (Munich), Germany
| | - Katia Parodi
- Department of Medical Physics, Ludwig-Maximilians-Universität München (LMU Munich), Garching (Munich), Germany
| | - Claus Belka
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany.,German Cancer Consortium (DKTK), Partner Site Munich, Munich, Germany
| | - Florian Kamp
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
| | - Christopher Kurz
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany.,Department of Medical Physics, Ludwig-Maximilians-Universität München (LMU Munich), Garching (Munich), Germany
| |
Collapse
|
13
|
Mori S, Hirai R, Sakata Y. Simulated four-dimensional CT for markerless tumor tracking using a deep learning network with multi-task learning. Phys Med 2020; 80:151-158. [PMID: 33189045 DOI: 10.1016/j.ejmp.2020.10.023] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/02/2020] [Revised: 10/16/2020] [Accepted: 10/24/2020] [Indexed: 10/23/2022] Open
Abstract
INTRODUCTION Our markerless tumor tracking algorithm requires 4DCT data to train models. 4DCT cannot be used for markerless tracking for respiratory-gated treatment due to inaccuracies and a high radiation dose. We developed a deep neural network (DNN) to generate 4DCT from 3DCT data. METHODS We used 2420 thoracic 4DCT datasets from 436 patients to train a DNN, designed to export 9 deformation vector fields (each field representing one-ninth of the respiratory cycle) from each CT dataset based on a 3D convolutional autoencoder with shortcut connections using deformable image registration. Then 3DCT data at exhale were transformed using the predicted deformation vector fields to obtain simulated 4DCT data. We compared markerless tracking accuracy between original and simulated 4DCT datasets for 20 patients. Our tracking algorithm used a machine learning approach with patient-specific model parameters. For the training stage, a pair of digitally reconstructed radiography images was generated using 4DCT for each patient. For the prediction stage, the tracking algorithm calculated tumor position using incoming fluoroscopic image data. RESULTS Diaphragmatic displacement averaged over 40 cases for the original 4DCT were slightly higher (<1.3 mm) than those for the simulated 4DCT. Tracking positional errors (95th percentile of the absolute value of displacement, "simulated 4DCT" minus "original 4DCT") averaged over the 20 cases were 0.56 mm, 0.65 mm, and 0.96 mm in the X, Y and Z directions, respectively. CONCLUSIONS We developed a DNN to generate simulated 4DCT data that are useful for markerless tumor tracking when original 4DCT is not available. Using this DNN would accelerate markerless tumor tracking and increase treatment accuracy in thoracoabdominal treatment.
Collapse
Affiliation(s)
- Shinichiro Mori
- Research Center for Charged Particle Therapy, National Institute of Radiological Sciences, Inage-ku, Chiba 263-8555, Japan.
| | - Ryusuke Hirai
- Corporate Research and Development Center, Toshiba Corporation, Kanagawa 212-8582, Japan
| | - Yukinobu Sakata
- Corporate Research and Development Center, Toshiba Corporation, Kanagawa 212-8582, Japan
| |
Collapse
|
14
|
Shan J, Yang Y, Schild SE, Daniels TB, Wong WW, Fatyga M, Bues M, Sio TT, Liu W. Intensity-modulated proton therapy (IMPT) interplay effect evaluation of asymmetric breathing with simultaneous uncertainty considerations in patients with non-small cell lung cancer. Med Phys 2020; 47:5428-5440. [PMID: 32964474 DOI: 10.1002/mp.14491] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2020] [Revised: 07/14/2020] [Accepted: 09/11/2020] [Indexed: 12/15/2022] Open
Abstract
PURPOSE Intensity-modulated proton therapy (IMPT) is sensitive to uncertainties from patient setup and proton beam range, as well as interplay effect. In addition, respiratory motion may vary from cycle to cycle, and also from day to day. These uncertainties can severely degrade the original plan quality and potentially affect patient's outcome. In this work, we developed a new tool to comprehensively consider the impact of all these uncertainties and provide plan robustness evaluation under them. METHODS We developed a comprehensive plan robustness evaluation tool that considered both uncertainties from patient setup and proton beam range, as well as respiratory motion simultaneously. To mimic patients' respiratory motion, the time spent in each phase was randomly sampled based on patient-specific breathing pattern parameters as acquired during the four-dimensional (4D)-computed tomography (CT) simulation. Spots were then assigned to one specific phase according to the temporal relationship between spot delivery sequence and patients' respiratory motion. Dose in each phase was calculated by summing contributions from all the spots delivered in that phase. The final 4D dynamic dose was obtained by deforming all doses in each phase to the maximum exhalation phase. Three hundred (300) scenarios (10 different breathing patterns with 30 different setup and range uncertainty scenario combinations) were calculated for each plan. The dose-volume histograms (DVHs) band method was used to assess plan robustness. Benchmarking the tool as an application's example, we compared plan robustness under both three-dimensional (3D) and 4D robustly optimized IMPT plans for 10 nonrandomly selected patients with non-small cell lung cancer. RESULTS The developed comprehensive plan robustness tool had been successfully applied to compare the plan robustness between 3D and 4D robustly optimized IMPT plans for 10 lung cancer patients. In the presence of interplay effect with uncertainties considered simultaneously, 4D robustly optimized plans provided significantly better CTV coverage (D95% , P = 0.002), CTV homogeneity (D5% -D95% , P = 0.002) with less target hot spots (D5% , P = 0.002), and target coverage robustness (CTV D95% bandwidth, P = 0.004) compared to 3D robustly optimized plans. Superior dose sparing of normal lung (lung Dmean , P = 0.020) favoring 4D plans and comparable normal tissue sparing including esophagus, heart, and spinal cord for both 3D and 4D plans were observed. The calculation time for all patients included in this study was 11.4 ± 2.6 min. CONCLUSION A comprehensive plan robustness evaluation tool was successfully developed and benchmarked for plan robustness evaluation in the presence of interplay effect, setup and range uncertainties. The very high efficiency of this tool marks its clinical adaptation, highly practical and versatile nature, including possible real-time intra-fractional interplay effect evaluation as a potential application for future use.
Collapse
Affiliation(s)
- Jie Shan
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, 85054, USA
| | - Yunze Yang
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, 85054, USA
| | - Steven E Schild
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, 85054, USA
| | - Thomas B Daniels
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, 85054, USA
| | - William W Wong
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, 85054, USA
| | - Mirek Fatyga
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, 85054, USA
| | - Martin Bues
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, 85054, USA
| | - Terence T Sio
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, 85054, USA
| | - Wei Liu
- Department of Radiation Oncology, Mayo Clinic, Phoenix, AZ, 85054, USA
| |
Collapse
|
15
|
Hoffmann A, Oborn B, Moteabbed M, Yan S, Bortfeld T, Knopf A, Fuchs H, Georg D, Seco J, Spadea MF, Jäkel O, Kurz C, Parodi K. MR-guided proton therapy: a review and a preview. Radiat Oncol 2020; 15:129. [PMID: 32471500 PMCID: PMC7260752 DOI: 10.1186/s13014-020-01571-x] [Citation(s) in RCA: 73] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2019] [Accepted: 05/17/2020] [Indexed: 02/14/2023] Open
Abstract
Background The targeting accuracy of proton therapy (PT) for moving soft-tissue tumours is expected to greatly improve by real-time magnetic resonance imaging (MRI) guidance. The integration of MRI and PT at the treatment isocenter would offer the opportunity of combining the unparalleled soft-tissue contrast and real-time imaging capabilities of MRI with the most conformal dose distribution and best dose steering capability provided by modern PT. However, hybrid systems for MR-integrated PT (MRiPT) have not been realized so far due to a number of hitherto open technological challenges. In recent years, various research groups have started addressing these challenges and exploring the technical feasibility and clinical potential of MRiPT. The aim of this contribution is to review the different aspects of MRiPT, to report on the status quo and to identify important future research topics. Methods Four aspects currently under study and their future directions are discussed: modelling and experimental investigations of electromagnetic interactions between the MRI and PT systems, integration of MRiPT workflows in clinical facilities, proton dose calculation algorithms in magnetic fields, and MRI-only based proton treatment planning approaches. Conclusions Although MRiPT is still in its infancy, significant progress on all four aspects has been made, showing promising results that justify further efforts for research and development to be undertaken. First non-clinical research solutions have recently been realized and are being thoroughly characterized. The prospect that first prototype MRiPT systems for clinical use will likely exist within the next 5 to 10 years seems realistic, but requires significant work to be performed by collaborative efforts of research groups and industrial partners.
Collapse
Affiliation(s)
- Aswin Hoffmann
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden-Rossendorf, Dresden, Germany.,Helmholtz-Zentrum Dresden-Rossendorf, Institute of Radiooncology-OncoRay, Dresden, Germany.,Department of Radiotherapy and Radiation Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Dresden, Germany
| | - Bradley Oborn
- Centre for Medical Radiation Physics, University of Wollongong, Wollongong, Australia.,Illawarra Cancer Care Centre, Wollongong Hospital, Wollongong, Australia
| | - Maryam Moteabbed
- Division of Radiation Biophysics, Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, USA
| | - Susu Yan
- Division of Radiation Biophysics, Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, USA
| | - Thomas Bortfeld
- Division of Radiation Biophysics, Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, USA
| | - Antje Knopf
- Department of Radiation Oncology, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands
| | - Herman Fuchs
- Department of Radiation Oncology, Medical University of Vienna/AKH, Vienna, Austria.,Christian Doppler Laboratory for Medical Radiation Research for Radiation Oncology, Medical University of Vienna, Vienna, Austria
| | - Dietmar Georg
- Department of Radiation Oncology, Medical University of Vienna/AKH, Vienna, Austria.,Christian Doppler Laboratory for Medical Radiation Research for Radiation Oncology, Medical University of Vienna, Vienna, Austria
| | - Joao Seco
- Biomedical Physics in Radiation Oncology, Deutsches Krebsforschungszentrum DKFZ, Heidelberg, Germany.,Department of Physics and Astronomy, Heidelberg University, Heidelberg, Germany
| | - Maria Francesca Spadea
- Biomedical Physics in Radiation Oncology, Deutsches Krebsforschungszentrum DKFZ, Heidelberg, Germany.,Department of Experimental and Clinical Medicine, Magna Graecia University, Catanzaro, Italy
| | - Oliver Jäkel
- Medical Physics in Radiation Oncology, Deutsches Krebsforschungszentrum DKFZ and Heidelberg Ion-Beam Therapy Center at the University Medical Center, Heidelberg, Germany
| | - Christopher Kurz
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany.,Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Garching, Germany
| | - Katia Parodi
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Garching, Germany.
| |
Collapse
|
16
|
4D strategies for lung tumors treated with hypofractionated scanning proton beam therapy: Dosimetric impact and robustness to interplay effects. Radiother Oncol 2020; 146:213-220. [DOI: 10.1016/j.radonc.2020.02.025] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2020] [Revised: 02/26/2020] [Accepted: 02/28/2020] [Indexed: 12/13/2022]
|