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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
Annunziata S, Rabe M, Vai A, Molinelli S, Nakas A, Meschini G, Pella A, Vitolo V, Barcellini A, Imparato S, Ciocca M, Orlandi E, Landry G, Kamp F, Kurz C, Baroni G, Riboldi M, Paganelli C. Virtual 4DCT generated from 4DMRI in gated particle therapy: phantom validation and application to lung cancer patients. Phys Med Biol 2023. [PMID: 37321258 DOI: 10.1088/1361-6560/acdec5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
OBJECTIVE Respiration negatively affects the outcome of a radiation therapy treatment, with potentially severe effects especially in particle therapy (PT). If compensation strategies are not applied, accuracy cannot be achieved. To support the clinical practice based on 4D Computed Tomography (CT), 4D Magnetic Resonance Imaging (MRI) acquisitions can be exploited. The purpose of this study was to validate a method for virtual 4DCT generation from 4DMRI data for lung cancers on a porcine lung phantom, and to apply it to lung cancer patients in PT.
Approach: Deformable image registration was used to register each respiratory phase of the 4DMRI to a reference phase. Then, a static 3DCT was registered to this reference MR image set, and the virtual 4DCT was generated by warping the registered CT according to previously obtained deformation fields. The method was validated on a physical phantom for which a ground truth 4DCT was available and tested on lung tumor patients, treated with gated PT at end-exhale, by comparing the virtual 4DCT with a re-evaluation 4DCT. The geometric and dosimetric evaluation was performed for both proton and carbon ion treatment plans.
Main results: The phantom validation exhibited a geometrical accuracy within the maximum resolution of the MRI and mean dose deviations, with respect to the prescription dose, up to 3.2% for target D95%, with a mean gamma pass rate of 98%. For patients, the virtual and re-evaluation 4DCTs showed good correspondence, with errors on target D95% up to 2% within the gating window. For one patient, dose variations up to 10% at end-exhale were observed due to relevant inter-fraction anatomo-pathological changes that occurred between the planning and re-evaluation CTs. 
Significance: Results obtained on phantom data showed that the virtual 4DCT method was accurate, allowing its application on patient data for testing within a clinical scenario.

.
Collapse
Affiliation(s)
- Simone Annunziata
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Via G. Colombo 40, Milano, 20133, ITALY
| | - Moritz Rabe
- Radiation Oncology, Klinikum der Universität München, Ziemssenstraße 1, Munchen, 80336, GERMANY
| | - Alessandro Vai
- National Centre of Oncological Hadrontherapy, Strada Campeggi, Pavia, Lombardia, 27100, ITALY
| | - Silvia Molinelli
- National Centre of Oncological Hadrontherapy, Strada Campeggi, Pavia, Lombardia, 27100, ITALY
| | - Anestis Nakas
- Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Via Giuseppe Colombo 40, Milano, Lombardia, 20133, ITALY
| | - Giorgia Meschini
- Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Via Giuseppe Colombo 40, Milano, Lombardia, 20133, ITALY
| | - Andrea Pella
- Bioengineering, Clinical Department, National Centre of Oncological Hadrontherapy, Strada Campeggi, Pavia, Lombardia, 27100, ITALY
| | - Viviana Vitolo
- National Centre of Oncological Hadrontherapy, Strada Campeggi, Pavia, Lombardia, 27100, ITALY
| | - Amelia Barcellini
- National Centre of Oncological Hadrontherapy, Strada Campeggi, Pavia, Lombardia, 27100, ITALY
| | - Sara Imparato
- National Centre of Oncological Hadrontherapy, Strada Campeggi, Pavia, Lombardia, 27100, ITALY
| | - Mario Ciocca
- Medical Physics, National Centre of Oncological Hadrontherapy, Strada Privata Campeggi, 53, Pavia, Lombardia, 27100, ITALY
| | - Ester Orlandi
- National Centre of Oncological Hadrontherapy, Strada Campeggi, Pavia, Lombardia, 27100, ITALY
| | - Guillaume Landry
- Radiation Oncology, University Hospital Munich Campus Grosshadern, Marchioninistr. 15, Munchen, 81377, GERMANY
| | - Florian Kamp
- Klinik und Poliklinik für Radioonkologie, Cyberknife und Strahlentherapie, Uniklinik Köln, Kerpener Str. 62, Koln, Nordrhein-Westfalen, 50937, GERMANY
| | - Christopher Kurz
- Radiation Oncology, University Hospital Munich, Ziemssenstraße 1, Munich, Bavaria, 80336, GERMANY
| | - Guido Baroni
- Department of Electronics Information and Bioengineering, Politecnico di Milano, P.zza Leonardo da Vinci 32,, Milano, Lombardia, 20133, ITALY
| | - Marco Riboldi
- Ludwig-Maximilians-Universität München, Geschwister-Scholl-Platz 1, Munchen, Bayern, 80539, GERMANY
| | - Chiara Paganelli
- Dipartimento di Elettronica Informazione e Bioingegneria, Politecnico di Milano, via G. Colombo 40, Milano, Lombardia, 20133, ITALY
| |
Collapse
|
3
|
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: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [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
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
molinelli S, Vai A, Russo S, Loap P, Meschini G, Magro G, Paganelli C, Barcellini A, Vitolo V, Ciocca M, Orlandi E. MO-0635 The need for treatment adaptation in carbon ion radiotherapy of pancreatic cancer. Radiother Oncol 2022. [DOI: 10.1016/s0167-8140(22)02393-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
Meschini G, Calabrese D, De Mori Bajolin F, Vai A, Fontana G, Molinelli S, Pella A, Imparato S, Vitolo V, Barcellini A, Orlandi E, Paganelli C, Baroni G. PO-1660 Investigating the generation of synthetic CT for abdominal tumors treated with particle therapy. Radiother Oncol 2021. [DOI: 10.1016/s0167-8140(21)08111-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
8
|
Meschini G, Paganelli C, Vai A, Fontana G, Molinelli S, Pella A, Vitolo V, Barcellini A, Orlandi E, Ciocca M, Riboldi M, Baroni G. An MRI framework for respiratory motion modelling validation. J Med Imaging Radiat Oncol 2021; 65:337-344. [PMID: 33773081 PMCID: PMC8251859 DOI: 10.1111/1754-9485.13175] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Revised: 01/27/2021] [Accepted: 03/07/2021] [Indexed: 12/20/2022]
Abstract
Introduction Respiratory motion models establish a correspondence between respiratory‐correlated (RC) 4‐dimensional (4D) imaging and respiratory surrogates, to estimate time‐resolved (TR) 3D breathing motion. To evaluate the performance of motion models on real patient data, a validation framework based on magnetic resonance imaging (MRI) is proposed, entailing the use of RC 4DMRI to build the model, and on both (i) TR 2D cine‐MRI and (ii) additional 4DMRI data for testing intra‐/inter‐fraction breathing motion variability. Methods Repeated MRI data were acquired in 7 patients with abdominal lesions. The considered model relied on deformable image registration (DIR) for building the model and compensating for inter‐fraction baseline variations. Both 2D and 3D validation were performed, by comparing model estimations with the ground truth 2D cine‐MRI and 4DMRI respiratory phases, respectively. Results The median DIR error was comparable to the voxel size (1.33 × 1.33 × 5 mm3), with higher values in the presence of large inter‐fraction motion (median value: 2.97 mm). In the 2D validation, the median estimation error on anatomical landmarks’ position resulted below 4 mm in every scenario, whereas in the 3D validation it was 1.33 mm and 4.21 mm when testing intra‐ and inter‐fraction motion, respectively. The range of motion described in the cine‐MRI was comparable to the motion of the building 4DMRI, being always above the estimation error. Overall, the model performance was dependent on DIR error, presenting reduced accuracy when inter‐fraction baseline variations occurred. Conclusions Results suggest the potential of the proposed framework in evaluating global motion models for organ motion management in MRI‐guided radiotherapy.
Collapse
Affiliation(s)
- Giorgia Meschini
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milan, Italy
| | - Chiara Paganelli
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milan, Italy
| | - Alessandro Vai
- National Centre for Oncological Hadrontherapy (CNAO), Pavia, Italy
| | - Giulia Fontana
- National Centre for Oncological Hadrontherapy (CNAO), Pavia, Italy
| | - Silvia Molinelli
- National Centre for Oncological Hadrontherapy (CNAO), Pavia, Italy
| | - Andrea Pella
- National Centre for Oncological Hadrontherapy (CNAO), Pavia, Italy
| | - Viviana Vitolo
- National Centre for Oncological Hadrontherapy (CNAO), Pavia, Italy
| | | | - Ester Orlandi
- National Centre for Oncological Hadrontherapy (CNAO), Pavia, Italy
| | - Mario Ciocca
- National Centre for Oncological Hadrontherapy (CNAO), Pavia, Italy
| | - Marco Riboldi
- Department of Medical Physics, Ludwig-Maximilians-Universität (LMU), Garching bei München, Germany
| | - Guido Baroni
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milan, Italy.,National Centre for Oncological Hadrontherapy (CNAO), Pavia, Italy
| |
Collapse
|
9
|
Elisei G, Pella A, Ricotti R, Via R, Fiore MR, Calvi G, Mastella E, Paganelli C, Tagaste B, Bello F, Fontana G, Meschini G, Buizza G, Valvo F, Orlandi E, Ciocca M, Baroni G. Development and validation of a new set-up simulator dedicated to ocular proton therapy at CNAO. Phys Med 2021; 82:228-239. [PMID: 33657472 DOI: 10.1016/j.ejmp.2021.01.070] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/16/2020] [Revised: 11/27/2020] [Accepted: 01/14/2021] [Indexed: 10/22/2022] Open
Abstract
An Eye Tracking System (ETS) is used at CNAO for providing a stable and reproducible ocular proton therapy (OPT) set-up, featuring a fixation light (FL) and monitoring stereo-cameras embedded in a rigid case. The aim of this work is to propose an ETS set-up simulation algorithm, that automatically provides the FL positioning in space, according to patient-specific gaze direction and avoiding interferences with patient, beam and collimator. Two configurations are provided: one in the CT room for acquiring images required for treatment planning with the patient lying on a couch, and one related to the treatment room with the patient sitting in front of the beam. Algorithm validation was performed reproducing ETS simulation (CT) and treatment (room) set-up for 30 patients previously treated at CNAO. The positioning accuracy of the device was quantified through a set of 14 control points applied to the ETS case and localizable both in the CT volume and in room X-ray images. Differences between the position of ETS reference points estimated by the algorithm and those measured by imaging systems are reported. The corresponding gaze direction deviation is on average 0.2° polar and 0.3° azimuth for positioning in CT room and 0.1° polar and 0.4° azimuth in the treatment room. The simulation algorithm was embedded in a clinically usable software application, which we assessed as capable of ensuring ETS positioning with an average accuracy of 2 mm in CT room and 1.5 mm in treatment room, corresponding to gaze direction deviations consistently lower than 1°.
Collapse
Affiliation(s)
- G Elisei
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department-Bioengineering Unit, Pavia, Italy
| | - A Pella
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department-Bioengineering Unit, Pavia, Italy.
| | - R Ricotti
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department-Bioengineering Unit, Pavia, Italy
| | - R Via
- Center of Proton Therapy, Paul Scherrer Institut, 5232 Villigen, PSI, Switzerland
| | - M R Fiore
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department, Pavia, Italy
| | - G Calvi
- Centro Nazionale di Adroterapia Oncologica CNAO, Particle Accelerator Department, Pavia, Italy
| | - E Mastella
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department - Medical Physics Unit, Pavia, Italy
| | - C Paganelli
- Politecnico di Milano, Department of Electronics Information and Bioengineering, Milano, Italy
| | - B Tagaste
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department-Bioengineering Unit, Pavia, Italy
| | - F Bello
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department-Bioengineering Unit, Pavia, Italy
| | - G Fontana
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department-Bioengineering Unit, Pavia, Italy
| | - G Meschini
- Politecnico di Milano, Department of Electronics Information and Bioengineering, Milano, Italy
| | - G Buizza
- Politecnico di Milano, Department of Electronics Information and Bioengineering, Milano, Italy
| | - F Valvo
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department, Pavia, Italy
| | - E Orlandi
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department, Pavia, Italy
| | - M Ciocca
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department - Medical Physics Unit, Pavia, Italy
| | - G Baroni
- Centro Nazionale di Adroterapia Oncologica CNAO, Clinical Department-Bioengineering Unit, Pavia, Italy; Politecnico di Milano, Department of Electronics Information and Bioengineering, Milano, Italy
| |
Collapse
|
10
|
Dumlu HS, Meschini G, Kurz C, Kamp F, Baroni G, Belka C, Paganelli C, Riboldi M. Dosimetric impact of geometric distortions in an MRI-only proton therapy workflow for lung, liver and pancreas. Z Med Phys 2020; 32:85-97. [PMID: 33168274 PMCID: PMC9948883 DOI: 10.1016/j.zemedi.2020.10.002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2020] [Revised: 09/02/2020] [Accepted: 10/01/2020] [Indexed: 12/25/2022]
Abstract
In a radiation therapy workflow based on Magnetic Resonance Imaging (MRI), dosimetric errors may arise due to geometric distortions introduced by MRI. The aim of this study was to quantify the dosimetric effect of system-dependent geometric distortions in an MRI-only workflow for proton therapy applied at extra-cranial sites. An approach was developed, in which computed tomography (CT) images were distorted using an MRI displacement map, which represented the MR distortions in a spoiled gradient-echo sequence due to gradient nonlinearities and static magnetic field inhomogeneities. A retrospective study was conducted on 4DCT/MRI digital phantoms and 18 4DCT clinical datasets of the thoraco-abdominal site. The treatment plans were designed and separately optimized for each beam in a beam specific Planning Target Volume on the distorted CT, and the final dose distribution was obtained as the average. The dose was then recalculated in undistorted CT using the same beam geometry and beam weights. The analysis was performed in terms of Dose Volume Histogram (DVH) parameters. No clinically relevant dosimetric impact was observed on organs at risk, whereas in the target structure, geometric distortions caused statistically significant variations in the planned dose DVH parameters and dose homogeneity index (DHI). The dosimetric variations in the target structure were smaller in abdominal cases (ΔD2%, ΔD98%, and ΔDmean all below 0.1% and ΔDHI below 0.003) compared to the lung cases. Indeed, lung patients with tumors isolated inside lung parenchyma exhibited higher dosimetric variations (ΔD2%≥0.3%, ΔD98%≥15.9%, ΔDmean≥3.3% and ΔDHI≥0.102) than lung patients with tumor close to soft tissue (ΔD2%≤0.4%, ΔD98%≤5.6%, ΔDmean≤0.9% and ΔDHI≤0.027) potentially due to higher density variations along the beam path. Results suggest the potential applicability of MRI-only proton therapy, provided that specific analysis is applied for isolated lung tumors.
Collapse
Affiliation(s)
- Hatice Selcen Dumlu
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Via Ponzio 34/5, 20133 Milano, Italy; Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Am Coulombwall 1, 85748 Garching bei München, Germany
| | - Giorgia Meschini
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Via Ponzio 34/5, 20133 Milano, Italy
| | - Christopher Kurz
- Department of Radiation Oncology, University Hospital, LMU Munich, Marchioninistraße 15, 81377 München, Germany
| | - Florian Kamp
- Department of Radiation Oncology, University Hospital, LMU Munich, Marchioninistraße 15, 81377 München, Germany
| | - Guido Baroni
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Via Ponzio 34/5, 20133 Milano, Italy; Centro Nazionale di Adroterapia Oncologica, Strada Campeggi 53, 27100 Pavia, Italy
| | - Claus Belka
- Department of Radiation Oncology, University Hospital, LMU Munich, Marchioninistraße 15, 81377 München, Germany; German Cancer Consortium (DKTK) partner site Munich, Germany and German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, 69120 Heidelberg, Germany
| | - Chiara Paganelli
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Via Ponzio 34/5, 20133 Milano, Italy
| | - Marco Riboldi
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Am Coulombwall 1, 85748 Garching bei München, Germany.
| |
Collapse
|
11
|
Meschini G, Vai A, Paganelli C, Molinelli S, Maestri D, Fontana G, Pella A, Vitolo V, Valvo F, Ciocca M, Baroni G. Investigating the use of virtual 4DCT from 4DMRI in gated carbon ion radiation therapy of abdominal tumors. Z Med Phys 2020; 32:98-108. [PMID: 33069586 PMCID: PMC9948849 DOI: 10.1016/j.zemedi.2020.08.005] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2020] [Revised: 07/27/2020] [Accepted: 08/31/2020] [Indexed: 12/22/2022]
Abstract
PURPOSE To generate virtual 4DCT from 4DMRI with field of view (FOV) extended to the entire involved patient anatomy, in order to evaluate its use in carbon ion radiation therapy (CIRT) of the abdominal site in a clinical scenario. MATERIALS AND METHODS The virtual 4DCT was generated by deforming a reference CT in order to (1) match the anatomy depicted in the 4DMRI within its FOV, by calculating deformation fields with deformable image registration to describe inter-fractional and breathing motion, and (2) obtain physically plausible deformation outside of the 4DMRI FOV, by propagating and modulating the previously obtained deformation fields. The implemented method was validated on a digital anthropomorphic phantom, for which a ground truth (GT) 4DCT was available. A CIRT treatment plan was optimized at the end-exhale reference CT and the RBE-weighted dose distribution was recalculated on both the virtual and GT 4DCTs. The method estimation error was quantified by comparing the virtual and GT 4DCTs and the corresponding recomputed doses. The method was then evaluated on 8 patients with pancreas or liver tumors treated with CIRT using respiratory gating at end-exhale. The clinical treatment plans adopted at the National Center for Oncological Hadrontherapy (CNAO, Pavia, Italy) were considered and the dose distribution was recomputed on all respiratory phases of the planning and virtual 4DCTs. By comparing the two datasets and the corresponding dose distributions, the geometrical and dosimetric impact of organ motion was assessed. RESULTS For the phantom, the error outside of the 4DMRI FOV was up to 4.5mm, but it remained sub-millimetric in correspondence to the target within the 4DMRI FOV. Although the impact of motion on the target D95% resulted in variations ranging from 22% to 90% between the planned dose and the doses recomputed on the GT 4DCT phases, the corresponding estimation error was ≤2.2%. In the patient cases, the variation of the baseline tumor position between the planning and the virtual end-exhale CTs presented a median (interquartile range) value of 6.0 (4.9) mm. For baseline variations larger than 5mm, the tumor D95% variation between the plan and the dose recomputed on the end-exhale virtual CT resulted larger than 10%. Median variations higher than 10% in the target D95% and gastro-intestinal OARs D2% were quantified at the end-inhale, whereas close to the end-exhale phase, limited variations of relevant dose metrics were found for both tumor and OARs. CONCLUSIONS The negligible impact of the geometrical inaccuracy in the estimated anatomy outside of the 4DMRI FOV on the overall dosimetric accuracy suggests the feasibility of virtual 4DCT with extended FOV in CIRT of the abdominal site. In the analyzed patient group, inter-fractional variations such as baseline variation and breathing variability were quantified, demonstrating the method capability to support treatment planning in gated CIRT of the abdominal site.
Collapse
Affiliation(s)
- Giorgia Meschini
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milano 20133, Italy.
| | - Alessandro Vai
- Centro Nazionale di Adroterapia Oncologica, Pavia 27100, Italy
| | - Chiara Paganelli
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milano 20133, Italy
| | | | - Davide Maestri
- Centro Nazionale di Adroterapia Oncologica, Pavia 27100, Italy
| | - Giulia Fontana
- Centro Nazionale di Adroterapia Oncologica, Pavia 27100, Italy
| | - Andrea Pella
- Centro Nazionale di Adroterapia Oncologica, Pavia 27100, Italy
| | - Viviana Vitolo
- Centro Nazionale di Adroterapia Oncologica, Pavia 27100, Italy
| | - Francesca Valvo
- Centro Nazionale di Adroterapia Oncologica, Pavia 27100, Italy
| | - Mario Ciocca
- Centro Nazionale di Adroterapia Oncologica, Pavia 27100, Italy
| | - Guido Baroni
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milano 20133, Italy,Centro Nazionale di Adroterapia Oncologica, Pavia 27100, Italy
| |
Collapse
|
12
|
Kalantzopoulos C, Meschini G, Paganelli C, Fontana G, Vai A, Preda L, Vitolo V, Valvo F, Baroni G. Organ motion quantification and margins evaluation in carbon ion therapy of abdominal lesions. Phys Med 2020; 75:33-39. [PMID: 32485596 DOI: 10.1016/j.ejmp.2020.05.014] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Revised: 04/03/2020] [Accepted: 05/17/2020] [Indexed: 12/11/2022] Open
Abstract
PURPOSE In image-guided particle radiotherapy of abdominal lesions, respiratory motion hinders treatment accuracy. In this study, 2D cineMRI data were used to quantify the tumor (GTV) motion and to evaluate the clinical approach based on deriving an internal target volume (ITV) from a planning 4DCT for gating treatments. METHODS Seven patients with abdominal lesions were treated with carbon-ion therapy at the National Centre of Oncological Hadron-therapy (Italy). The MR scan was performed on the same day of the 4DCT acquisition. For four patients, an additional MR was acquired approximately after 1 week. The cineMRI combined with deformable image registration algorithm was used to quantify tumor motion. Afterwards, two ITVs were defined considering (1) all phases (ITVFB) and (2) only phases within the gating window (ITVG), and then compared with the clinical (4DCT-derived) ITVs (ITVCG and ITVCFB). RESULTS Tumor residual motion estimated by cineMRI data in the two MRI sessions resulted not significantly different from 4DCT, although cineMRI accounted for cycle-to-cycle variations. The ITV normalized for the GTV median values were higher for ITVFB with respect to ITVG, ITVCFB and ITVCG. The Hausdorff distances with respect to the GTV were up to 10.55 mm, 3.13 mm, 5.56 mm and 2.51 mm, for ITVFB, ITVG, ITVCFB and ITVCG, respectively. According to both metrics, ITVCG and ITVG were not found significantly different. CONCLUSIONS CineMRI acquisitions allowed to quantify organ motion without delivering additional dose to the patient and to verify treatment margins in gated carbon-ion therapy of abdominal lesions.
Collapse
Affiliation(s)
| | - Giorgia Meschini
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, 20133 Milan, Italy
| | - Chiara Paganelli
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, 20133 Milan, Italy
| | - Giulia Fontana
- Centro Nazionale di Adroterapia Oncologica, Str. Campeggi, 53, 27100 Pavia, Italy
| | - Alessandro Vai
- Centro Nazionale di Adroterapia Oncologica, Str. Campeggi, 53, 27100 Pavia, Italy
| | - Lorenzo Preda
- Centro Nazionale di Adroterapia Oncologica, Str. Campeggi, 53, 27100 Pavia, Italy
| | - Viviana Vitolo
- Centro Nazionale di Adroterapia Oncologica, Str. Campeggi, 53, 27100 Pavia, Italy
| | - Francesca Valvo
- Centro Nazionale di Adroterapia Oncologica, Str. Campeggi, 53, 27100 Pavia, Italy
| | - Guido Baroni
- Centro Nazionale di Adroterapia Oncologica, Str. Campeggi, 53, 27100 Pavia, Italy; Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, 20133 Milan, Italy
| |
Collapse
|
13
|
Meschini G, Kamp F, Hofmaier J, Reiner M, Sharp G, Paganetti H, Belka C, Wilkens JJ, Carlson DJ, Parodi K, Baroni G, Riboldi M. Modeling RBE-weighted dose variations in irregularly moving abdominal targets treated with carbon ion beams. Med Phys 2020; 47:2768-2778. [PMID: 32162332 DOI: 10.1002/mp.14135] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2019] [Revised: 03/09/2020] [Accepted: 03/09/2020] [Indexed: 01/01/2023] Open
Abstract
PURPOSE To model four-dimensional (4D) relative biological effectiveness (RBE)-weighted dose variations in abdominal lesions treated with scanned carbon ion beam in case of irregular breathing motion. METHODS The proposed method, referred to as bioWED method, combines the simulation of tumor motion in a patient- and beam-specific water equivalent depth (WED)-space with RBE modeling, aiming at the estimation of RBE-weighted dose changes due to respiratory motion. The method was validated on a phantom, simulating gated and free breathing dose delivery, and on a patient case, for which free breathing irradiation was assumed and both amplitude and baseline breathing irregularities were simulated through a respiratory motion model. We quantified (a) the effect of motion on the equivalent uniform dose (EUD) and the RBE-weighted dose-volume histograms (DVH), by comparing the planned dose distribution with "ground truth" 4D RBE-weighted doses computed using 4D computed tomography data, and (ii) the estimation error, by comparing the doses estimated with the bioWED method to "ground truth" 4D RBE-weighted doses. RESULTS In the phantom validation, the estimation error on the EUD was limited with respect to the motion effect and the median estimation error on relevant RBE-weighted DVH metrics remained within 5%. In the patient study, the estimation error as computed on the EUD was smaller than the corresponding motion effect, exhibiting the largest values in the baseline irregularity simulation. However, the median estimation error over all simulations was below 3.2% considering relevant DVH metrics. CONCLUSIONS In the evaluated cases, the bioWED method showed proper accuracy when compared to deformable image registration-based 4D dose calculation. Therefore, it can be seen as a tool to test treatment plan robustness against irregular breathing motion, although its accuracy decreases as a function of increasing soft tissue deformation and should be evaluated on a larger patient dataset.
Collapse
Affiliation(s)
- Giorgia Meschini
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milan, Italy
| | - Florian Kamp
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
| | - Jan Hofmaier
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
| | - Michael Reiner
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
| | - Gregory Sharp
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Harald Paganetti
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Claus Belka
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
| | - Jan J Wilkens
- Department of Radiation Oncology, School of Medicine, Technical University of Munich, Klinikum rechts der Isar, Munich, Germany
| | - David J Carlson
- Yale University, New Haven, CT, USA.,University of Pennsylvania, Philadelphia, PA, USA
| | - Katia Parodi
- Department of Experimental Physics -Medical Physics, Ludwig-Maximilians-Universität München (LMU Munich), Munich, Germany
| | - Guido Baroni
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milan, Italy.,Centro Nazionale di Adroterapia Oncologica, Pavia, Italy
| | - Marco Riboldi
- Department of Experimental Physics -Medical Physics, Ludwig-Maximilians-Universität München (LMU Munich), Munich, Germany
| |
Collapse
|
14
|
Meschini G, Vai A, Paganelli C, Molinelli S, Fontana G, Pella A, Preda L, Vitolo V, Valvo F, Ciocca M, Riboldi M, Baroni G. Virtual 4DCT from 4DMRI for the management of respiratory motion in carbon ion therapy of abdominal tumors. Med Phys 2020; 47:909-916. [PMID: 31880819 DOI: 10.1002/mp.13992] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2019] [Revised: 12/17/2019] [Accepted: 12/17/2019] [Indexed: 12/12/2022] Open
Abstract
PURPOSE To evaluate a method for generating virtual four-dimensional computed tomography (4DCT) from four-dimensional magnetic resonance imaging (4DMRI) data in carbon ion radiotherapy with pencil beam scanning for abdominal tumors. METHODS Deformable image registration is used to: (a) register each respiratory phase of the 4DMRI to the end-exhale MRI; (b) register the reference end-exhale CT to the end-exhale MRI volume; (c) generate the virtual 4DCT by warping the registered CT according to the obtained deformation fields. A respiratory-gated carbon ion treatment plan is optimized on the planning 4DCT and the corresponding dose distribution is recalculated on the virtual 4DCT. The method was validated on a digital anthropomorphic phantom and tested on eight patients (18 acquisitions). For the phantom, a ground truth dataset was available to assess the method performances from the geometrical and dosimetric standpoints. For the patients, the virtual 4DCT was compared with the planning 4DCT. RESULTS In the phantom, the method exhibits a geometrical accuracy within the voxel size and Dose Volume Histograms deviations up to 3.3% for target V95% (mean dose difference ≤ 0.2% of the prescription dose, gamma pass rate > 98%). For patients, the virtual and the planning 4DCTs show good agreement at end-exhale (3% median D95% difference), whereas other respiratory phases exhibit moderate motion variability with consequent dose discrepancies, confirming the need for motion mitigation strategies during treatment. CONCLUSIONS The virtual 4DCT approach is feasible to evaluate treatment plan robustness against intra- and interfraction motion in carbon ion therapy delivered at the abdominal site.
Collapse
Affiliation(s)
- Giorgia Meschini
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milano, 20133, Italy
| | - Alessandro Vai
- Centro Nazionale di Adroterapia Oncologica, Pavia, 27100, Italy
| | - Chiara Paganelli
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milano, 20133, Italy
| | | | - Giulia Fontana
- Centro Nazionale di Adroterapia Oncologica, Pavia, 27100, Italy
| | - Andrea Pella
- Centro Nazionale di Adroterapia Oncologica, Pavia, 27100, Italy
| | - Lorenzo Preda
- Centro Nazionale di Adroterapia Oncologica, Pavia, 27100, Italy.,Department of Clinical, Surgical, Diagnostic and Pediatric Sciences, University of Pavia, Pavia, 27100, Italy
| | - Viviana Vitolo
- Centro Nazionale di Adroterapia Oncologica, Pavia, 27100, Italy
| | - Francesca Valvo
- Centro Nazionale di Adroterapia Oncologica, Pavia, 27100, Italy
| | - Mario Ciocca
- Centro Nazionale di Adroterapia Oncologica, Pavia, 27100, Italy
| | - Marco Riboldi
- Chair of Experimental Physics - Medical Physics, Ludwig-Maximilians-Universität (LMU), Munich, 80539, Germany
| | - Guido Baroni
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milano, 20133, Italy.,Centro Nazionale di Adroterapia Oncologica, Pavia, 27100, Italy
| |
Collapse
|
15
|
Paganelli C, Portoso S, Garau N, Meschini G, Via R, Buizza G, Keall P, Riboldi M, Baroni G. Time-resolved volumetric MRI in MRI-guided radiotherapy: an in silico comparative analysis. Phys Med Biol 2019; 64:185013. [PMID: 31323645 DOI: 10.1088/1361-6560/ab33e5] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
MRI-treatment units enable 2D cine-MRI centred in the tumour for motion detection in radiotherapy, but they lack 3D information due to spatio-temporal limits. To derive time-resolved 3D information, different approaches have been proposed in the literature, but a rigorous comparison among these strategies has not yet been performed. The goal of this study is to quantitatively investigate five published strategies that derive time-resolved volumetric MRI in MRI-guided radiotherapy: Propagation, out-of-plane motion compensation, Fayad model, ROI-based model and Stemkens model. Comparisons were performed using an MRI digital phantom generated with six different patient-derived motion signals and tumour-shapes. An average 4D cycle was generated as well as 2D cine-MRI data with corresponding 3D in-room ground truth. Quantitative analysis was performed by comparing the estimated 3D volume to the ground truth available for each 2D cine-MRI sample. A grouped patient statistical analysis was performed to evaluate the performance of the selected methods, in case of tumour tracking or motion estimation of the whole anatomy. Analyses were also performed based on patient characteristics. Quantitative ranking of the investigated methods highlighted that Propagation and ROI-based model strategies achieved an overall median tumour centre of mass 3D distance from the ground truth of 1.1 mm and 1.3 mm, respectively, and a diaphragm distance below 1.6 mm. Higher errors and variabilities were instead obtained for other methods, which lack the ability to compensate for in-room variations and to account for regional changes. These results were especially evident when further analysing patient characteristics, where errors above 2 mm/5 mm in tumour/diaphragm were found for more irregular breathing patterns in case of out-of-plane motion compensation, Fayad and Stemkens models. These findings suggest the potential of the proposed in silico framework to develop and compare strategies to estimate time-resolved 3DMRI in MRI-guided radiotherapy.
Collapse
Affiliation(s)
- C Paganelli
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milano, Italy. Both authors contributed equally. Author to whom any correspondence should be addressed
| | | | | | | | | | | | | | | | | |
Collapse
|
16
|
Meschini G, Seregni M, Molinelli S, Vai A, Phillips J, Sharp GC, Pella A, Valvo F, Ciocca M, Riboldi M, Paganetti H, Baroni G. Validation of a model for physical dose variations in irregularly moving targets treated with carbon ion beams. Med Phys 2019; 46:3663-3673. [PMID: 31206718 DOI: 10.1002/mp.13662] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2018] [Revised: 06/07/2019] [Accepted: 06/07/2019] [Indexed: 12/14/2022] Open
Abstract
PURPOSE In particle therapy, conventional treatment planning systems rely on an imaging representation of the irradiated region to compute the dose. For irregular breathing, when an imaging dataset describing the actual motion is not available, a different approach for dose estimation is needed. To this aim, we validate a method for the estimation of physical dose variations in gated carbon ion treatments, providing also a demonstration of the feasibility of physical dose metrics to assess the method performance. Finally, we describe a sample use case, in which this method is used to assess plan robustness with respect to undetected irregular tumor motion. METHODS The method entails the definition of a patient- and beam-specific water equivalent depth (WED) space, the simulation of motion as a translation equal to tumor displacement, and the reconstruction of the altered dose. We validated the approach using four-dimensional computed tomographies (4DCTs) and clinical plans in 12 patients, treated with respiratory gated carbon ion beams at the National Centre for Oncological Hadrontherapy (Pavia, Italy). Using the end-exhale CT and dose distribution as a reference, the physical dose delivered at the end-inhale tumor position was estimated and compared to the ground-truth dose recalculation on the end-inhale CT. Biologically effective and physical dose variations between the plan and the recalculation were compared as well. As a use case, we evaluated dose changes caused by simulated irregular tumor motion, that is, linear and nonlinear baseline shifts and/or amplitude variations with hysteresis. RESULTS The ratio between biologically effective and physical equivalent uniform dose (EUD) variations due to end-exhale to end-inhale motion was less than one for 96% of investigated structures. In the validation study, we found a median error corresponding to a 14% EUD overestimation for the tumor and 4% EUD underestimation for a subgroup of organs at risk, together with a high EUD variation due to motion [median 352% EUD variation between end-exhale and end-inhale doses in the planning tumor volume (PTV)]. Considering relevant dose-volume histogram (DVH) metrics, the median difference between estimated and ground truth doses was ≤ 4%. Gamma analysis between estimated and recalculated dose distributions resulted in a pass rate > 80% for 83% of the target volumes. For the two patients selected for the sample use case, a patient-specific assessment of the method performance was performed on the 4DCT and it was possible to relate EUD variations of both tumor and organs at risk to the simulated target motion. CONCLUSIONS The physical dose distribution was found to be more sensitive to motion with respect to the biologically effective one, suggesting the suitability of the physical dose metrics for the WED-space method validation. We showed that the method can compensate for intra-fractional tumor motion with proper accuracy in the selected patient group, although its use is recommended when limited deformations are expected. In conclusion, the WED-space method can provide simulations of dose alteration due to irregular breathing when imaging data are lacking, and, once integrated with relative biological effectiveness (RBE) modeling, it would be useful in evaluating the robustness of carbon ion treatment plans.
Collapse
Affiliation(s)
| | | | | | - Alessandro Vai
- Centro Nazionale Adroterapia Oncologica, Pavia, 27100, Italy
| | - Justin Phillips
- Alexian Brothers Medical Center, Elk Grove Village, IL, 60007, USA
| | | | - Andrea Pella
- Centro Nazionale Adroterapia Oncologica, Pavia, 27100, Italy
| | - Francesca Valvo
- Centro Nazionale Adroterapia Oncologica, Pavia, 27100, Italy
| | - Mario Ciocca
- Centro Nazionale Adroterapia Oncologica, Pavia, 27100, Italy
| | - Marco Riboldi
- Ludwig-Maximilians-Universität, Munich, 80539, Germany
| | | | - Guido Baroni
- Politecnico di Milano, Milan, 20133, Italy.,Centro Nazionale Adroterapia Oncologica, Pavia, 27100, Italy
| |
Collapse
|
17
|
Vai A, Meschini G, Molinelli S, Paganelli C, Maestri D, Magro G, Mastella E, Mairani A, Mirandola A, Russo S, Preda L, Viselener G, Barcellini A, Vitolo V, Mancin A, Fontana G, Baroni G, Ciocca M. EP-1968 Respiratory-gated carbon-ion beam treatments of abdominal targets: clinical introduction of 4DMRI. Radiother Oncol 2019. [DOI: 10.1016/s0167-8140(19)32388-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
18
|
Meschini G, Paganelli C, Fontana G, Pella A, Mancin A, Vai A, Riboldi M, Valvo F, Baroni G. PO-0976 Validation of respiratory motion modeling through repeated 4DMRI in the abdomen: preliminary results. Radiother Oncol 2019. [DOI: 10.1016/s0167-8140(19)31396-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
19
|
Meschini G, Paganelli C, Gianoli C, Summers P, Bellomi M, Baroni G, Riboldi M. A clustering approach to 4D MRI retrospective sorting for the investigation of different surrogates. Phys Med 2019; 58:107-113. [PMID: 30824141 DOI: 10.1016/j.ejmp.2019.02.003] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/12/2018] [Revised: 01/19/2019] [Accepted: 02/06/2019] [Indexed: 12/25/2022] Open
Abstract
PURPOSE In retrospective 4-Dimensional Magnetic Resonance Imaging (4D MRI) sorting, respiratory surrogate selection affects the image quality of reconstructed volumes. We propose a method for retrospective 4D MRI sorting based on clustering, which allowed us to compare the performance of single or multiple internal surrogates vs. a conventional external signal. METHODS A k-medoids clustering algorithm was exploited for sorting 2D MRI into 4D MRI, relying on (A) multiple or (B) single automatically tracked internal landmarks or (C) respiratory belt signal. 4D MRI reconstructions for seven liver cancer patients were compared to those of the state-of-the-art mutual information (MI) approach. Sorting artifacts were measured by the root mean square error (RMSE) between the diaphragm profile and a fitted second order curve. Diaphragm and tumor motions were evaluated. RESULTS The median RMSEs ranged 0.97-1.66 mm, 1.24-1.89 mm, 1.43-2.27 mm, 1.74-3.72 mm for the MI, (A), (B) and (C) methods, respectively. Significant differences (Friedman, α = 5%) were found between (C) and all other methods, and between (B) and MI approaches. The discrepancies between (A) and MI approaches ranged 1.1-6.2 mm and 0.7-5.3 mm respectively in diaphragm and tumor motions. Methods (A) and (B) showed similar ranges of motion. CONCLUSION With multiple internal points, our method yielded the description of a higher range of motion and similar image quality with respect to the MI approach. The single point method led to more artifacts, suggesting the superior suitability of multiple internal surrogates for retrospective 4D MRI sorting. Considering internal rather than external information favored superior performance.
Collapse
Affiliation(s)
- Giorgia Meschini
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci, 32, 20133 Milan, Italy.
| | - Chiara Paganelli
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci, 32, 20133 Milan, Italy
| | - Chiara Gianoli
- Chair of Experimental Physics - Medical Physics, Ludwig-Maximilians-Universität München, Am Coulombwall 1, 85748 Garching bei München, Germany
| | - Paul Summers
- Department of Imaging and Radiological Science, European Institute of Oncology, Via Giuseppe Ripamonti, 435, 20141 Milan, Italy
| | - Massimo Bellomi
- Department of Imaging and Radiological Science, European Institute of Oncology, Via Giuseppe Ripamonti, 435, 20141 Milan, Italy; Department of Oncology and Emato-oncology, University of Milan, Via Festa del Perdono, 7, 20122, Italy
| | - Guido Baroni
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci, 32, 20133 Milan, Italy; Bioengineering Unit, CNAO Foundation, Str. Campeggi, 53, 27100 Pavia, Italy
| | - Marco Riboldi
- Chair of Experimental Physics - Medical Physics, Ludwig-Maximilians-Universität München, Am Coulombwall 1, 85748 Garching bei München, Germany
| |
Collapse
|
20
|
Garau N, Via R, Meschini G, Lee D, Keall P, Riboldi M, Baroni G, Paganelli C. A ROI-based global motion model established on 4DCT and 2D cine-MRI data for MRI-guidance in radiation therapy. Phys Med Biol 2019; 64:045002. [PMID: 30625459 DOI: 10.1088/1361-6560/aafcec] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
In-room magnetic resonance imaging (MRI) allows the acquisition of fast 2D cine-MRI centered in the tumor for advanced motion management in radiotherapy. To achieve 3D information during treatment, patient-specific motion models can be considered the most viable solution. However, conventional global motion models are built using a single motion surrogate, independently from the anatomical location. In this work, we present a novel motion model based on regions of interest (ROIs) established on 4D computed tomography (4DCT) and 2D cine-MRI, aiming at accurately compensating for changes during treatment. In the planning phase, a motion model is built on a 4DCT dataset, through 3D deformable image registration (DIR). ROIs are then defined and correlated with motion fields derived by 2D DIR between CT slices centered in the tumor. In the treatment phase, the model is applied to in-room cine-MRI data to compensate for organ motion in a multi-modal framework, aiming at estimating a time-resolved 3DCT. The method is validated on a digital phantom and tested on two lung patients. Analysis is performed by considering different anatomical planes (coronal, sagittal and a combination of the two) and evaluating the performance of the method on tumor and diaphragm. For the phantom study, the ROI-based model results in a uniform median error on both diaphragm and tumor below 1.5 mm. For what concerns patients, median errors on both diaphragm and tumor are around 2 mm (maximum patient resolution), confirming the capability of the method to regionally compensate for motion. A novel ROI-based motion model is proposed as an integral part of an envisioned clinical MRI-guided workflow aiming at enhanced image guidance compared to conventional strategies.
Collapse
Affiliation(s)
- Noemi Garau
- Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milano, Italy. Author to whom any correspondence should be addressed
| | | | | | | | | | | | | | | |
Collapse
|
21
|
Paganelli C, Meschini G, Molinelli S, Riboldi M, Baroni G. “Patient-specific validation of deformable image registration in radiation therapy: Overview and caveats”. Med Phys 2018; 45:e908-e922. [DOI: 10.1002/mp.13162] [Citation(s) in RCA: 58] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2017] [Revised: 07/30/2018] [Accepted: 08/24/2018] [Indexed: 12/26/2022] Open
Affiliation(s)
- Chiara Paganelli
- Dipartimento di Elettronica, Informazione e Bioingegneria; Politecnico di Milano; Milano 20133 Italy
| | - Giorgia Meschini
- Dipartimento di Elettronica, Informazione e Bioingegneria; Politecnico di Milano; Milano 20133 Italy
| | | | - Marco Riboldi
- Department of Medical Physics; Ludwig-Maximilians-Universitat Munchen; Munich 80539 Germany
| | - Guido Baroni
- Dipartimento di Elettronica, Informazione e Bioingegneria; Politecnico di Milano; Milano 20133 Italy
- Centro Nazionale di Adroterapia Oncologica; Pavia 27100 Italy
| |
Collapse
|
22
|
Meschini G, Vai A, Riboldi M, Pella A, Vitolo V, Caivano D, Ciocca M, Baroni G, Molinelli S. [P203] WEPL-corrected margins for moving targets in carbon-ion gated treatments: Preliminary results. Phys Med 2018. [DOI: 10.1016/j.ejmp.2018.06.497] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/28/2022] Open
|
23
|
Garau N, Paganelli C, Meschini G, Via R, Riboldi M, Baroni G. OC-0188: A ROI-based global motion model for MRI-guidance in radiation therapy: a phantom study. Radiother Oncol 2018. [DOI: 10.1016/s0167-8140(18)30498-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
24
|
Meschini G, Seregni M, Pella A, Baroni G, Riboldi M. SU-F-J-80: Deformable Image Registration for Residual Organ Motion Estimation in Respiratory Gated Treatments with Scanned Carbon Ion Beams. Med Phys 2016. [DOI: 10.1118/1.4955988] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|