1
|
Goodburn RJ, Philippens MEP, Lefebvre TL, Khalifa A, Bruijnen T, Freedman JN, Waddington DEJ, Younus E, Aliotta E, Meliadò G, Stanescu T, Bano W, Fatemi‐Ardekani A, Wetscherek A, Oelfke U, van den Berg N, Mason RP, van Houdt PJ, Balter JM, Gurney‐Champion OJ. The future of MRI in radiation therapy: Challenges and opportunities for the MR community. Magn Reson Med 2022; 88:2592-2608. [PMID: 36128894 PMCID: PMC9529952 DOI: 10.1002/mrm.29450] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2022] [Revised: 08/17/2022] [Accepted: 08/22/2022] [Indexed: 01/11/2023]
Abstract
Radiation therapy is a major component of cancer treatment pathways worldwide. The main aim of this treatment is to achieve tumor control through the delivery of ionizing radiation while preserving healthy tissues for minimal radiation toxicity. Because radiation therapy relies on accurate localization of the target and surrounding tissues, imaging plays a crucial role throughout the treatment chain. In the treatment planning phase, radiological images are essential for defining target volumes and organs-at-risk, as well as providing elemental composition (e.g., electron density) information for radiation dose calculations. At treatment, onboard imaging informs patient setup and could be used to guide radiation dose placement for sites affected by motion. Imaging is also an important tool for treatment response assessment and treatment plan adaptation. MRI, with its excellent soft tissue contrast and capacity to probe functional tissue properties, holds great untapped potential for transforming treatment paradigms in radiation therapy. The MR in Radiation Therapy ISMRM Study Group was established to provide a forum within the MR community to discuss the unmet needs and fuel opportunities for further advancement of MRI for radiation therapy applications. During the summer of 2021, the study group organized its first virtual workshop, attended by a diverse international group of clinicians, scientists, and clinical physicists, to explore our predictions for the future of MRI in radiation therapy for the next 25 years. This article reviews the main findings from the event and considers the opportunities and challenges of reaching our vision for the future in this expanding field.
Collapse
Affiliation(s)
- Rosie J. Goodburn
- Joint Department of PhysicsInstitute of Cancer Research and Royal Marsden NHS Foundation TrustLondonUnited Kingdom
| | | | - Thierry L. Lefebvre
- Department of PhysicsUniversity of CambridgeCambridgeUnited Kingdom
- Cancer Research UK Cambridge Research InstituteUniversity of CambridgeCambridgeUnited Kingdom
| | - Aly Khalifa
- Department of Medical BiophysicsUniversity of TorontoTorontoOntarioCanada
| | - Tom Bruijnen
- Department of RadiotherapyUniversity Medical Center UtrechtUtrechtNetherlands
| | | | - David E. J. Waddington
- Faculty of Medicine and Health, Sydney School of Health Sciences, ACRF Image X InstituteThe University of SydneySydneyNew South WalesAustralia
| | - Eyesha Younus
- Department of Medical Physics, Odette Cancer CentreSunnybrook Health Sciences CentreTorontoOntarioCanada
| | - Eric Aliotta
- Department of Medical PhysicsMemorial Sloan Kettering Cancer CenterNew YorkNew YorkUSA
| | - Gabriele Meliadò
- Unità Operativa Complessa di Fisica SanitariaAzienda Ospedaliera Universitaria Integrata VeronaVeronaItaly
| | - Teo Stanescu
- Department of Radiation Oncology, University of Toronto and Medical Physics, Princess Margaret Cancer CentreUniversity Health NetworkTorontoOntarioCanada
| | - Wajiha Bano
- Joint Department of PhysicsInstitute of Cancer Research and Royal Marsden NHS Foundation TrustLondonUnited Kingdom
| | - Ali Fatemi‐Ardekani
- Department of PhysicsJackson State University (JSU)JacksonMississippiUSA
- SpinTecxJacksonMississippiUSA
- Department of Radiation OncologyCommunity Health Systems (CHS) Cancer NetworkJacksonMississippiUSA
| | - Andreas Wetscherek
- Joint Department of PhysicsInstitute of Cancer Research and Royal Marsden NHS Foundation TrustLondonUnited Kingdom
| | - Uwe Oelfke
- Joint Department of PhysicsInstitute of Cancer Research and Royal Marsden NHS Foundation TrustLondonUnited Kingdom
| | - Nico van den Berg
- Department of RadiotherapyUniversity Medical Center UtrechtUtrechtNetherlands
| | - Ralph P. Mason
- Department of RadiologyUniversity of Texas Southwestern Medical CenterDallasTexasUSA
| | - Petra J. van Houdt
- Department of Radiation OncologyNetherlands Cancer InstituteAmsterdamNetherlands
| | - James M. Balter
- Department of Radiation OncologyUniversity of MichiganAnn ArborMichiganUSA
| | - Oliver J. Gurney‐Champion
- Imaging and Biomarkers, Cancer Center Amsterdam, Amsterdam UMCUniversity of AmsterdamAmsterdamNetherlands
| |
Collapse
|
2
|
Retif P, Djibo Sidikou A, Mathis C, Letellier R, Verrecchia-Ramos E, Dupres R, Michel X. Evaluation of the ability of the Brainlab Elements Cranial Distortion Correction algorithm to correct clinically relevant MRI distortions for cranial SRT. Strahlenther Onkol 2022; 198:907-918. [PMID: 35980455 DOI: 10.1007/s00066-022-01988-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2021] [Accepted: 07/10/2022] [Indexed: 11/29/2022]
Abstract
PURPOSE Cranial stereotactic radiotherapy (SRT) requires highly accurate lesion delineation. However, MRI can have significant inherent geometric distortions. We investigated how well the Elements Cranial Distortion Correction algorithm of Brainlab (Munich, Germany) corrects the distortions in MR image-sets of a phantom and patients. METHODS A non-distorted reference computed tomography image-set of a CIRS Model 603-GS (CIRS, Norfolk, VA, USA) phantom was acquired. Three-dimensional T1-weighted images were acquired with five MRI scanners and reconstructed with vendor-derived distortion correction. Some were reconstructed without correction to generate heavily distorted image-sets. All MR image-sets were corrected with the Brainlab algorithm relative to the computed tomography acquisition. CIRS Distortion Check software measured the distortion in each image-set. For all uncorrected and corrected image-sets, the control points that exceeded the 0.5-mm clinically relevant distortion threshold and the distortion maximum, mean, and standard deviation were recorded. Empirical cumulative distribution functions (eCDF) were plotted. Intraclass correlation coefficient (ICC) was calculated. The algorithm was evaluated with 10 brain metastases using Dice similarity coefficients (DSC). RESULTS The algorithm significantly reduced mean and standard deviation distortion in all image-sets. It reduced the maximum distortion in the heavily distorted image-sets from 2.072 to 1.059 mm and the control points with > 0.5-mm distortion fell from 50.2% to 4.0%. Before and especially after correction, the eCDFs of the four repeats were visually similar. ICC was 0.812 (excellent-good agreement). The algorithm increased the DSCs for all patients and image-sets. CONCLUSION The Brainlab algorithm significantly and reproducibly ameliorated MRI distortion, even with heavily distorted images. Thus, it increases the accuracy of cranial SRT lesion delineation. After further testing, this tool may be suitable for SRT of small lesions.
Collapse
Affiliation(s)
- Paul Retif
- Medical Physics Unit, CHR Metz-Thionville, Metz, France. .,Université de Lorraine, CNRS, CRAN, 54000, Nancy, France.
| | | | | | | | | | - Rémi Dupres
- Medical Imaging Department, CHR Metz-Thionville, Metz, France
| | - Xavier Michel
- Radiation Therapy Department, CHR Metz-Thionville, Metz, France
| |
Collapse
|
3
|
Theocharis S, Pappas EP, Seimenis I, Kouris P, Dellios D, Kollias G, Karaiskos P. Geometric distortion assessment in 3T MR images used for treatment planning in cranial Stereotactic Radiosurgery and Radiotherapy. PLoS One 2022; 17:e0268925. [PMID: 35605005 PMCID: PMC9126373 DOI: 10.1371/journal.pone.0268925] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Accepted: 05/10/2022] [Indexed: 12/31/2022] Open
Abstract
Magnetic Resonance images (MRIs) are employed in brain Stereotactic Radiosurgery and Radiotherapy (SRS/SRT) for target and/or critical organ localization and delineation. However, MRIs are inherently distorted, which also impacts the accuracy of the Magnetic Resonance Imaging/Computed Tomography (MRI/CT) co-registration process. In this phantom-based study, geometric distortion is assessed in 3T T2-weighted images (T2WIs), while the efficacy of an MRI distortion correction technique is also evaluated. A homogeneous polymer gel-filled phantom was CT-imaged before being irradiated with 26 4-mm Gamma Knife shots at predefined locations (reference control points). The irradiated phantom was MRI-scanned at 3T, implementing a T2-weighted protocol suitable for SRS/SRT treatment planning. The centers of mass of all shots were identified in the 3D image space by implementing an iterative localization algorithm and served as the evaluated control points for MRI distortion detection. MRIs and CT images were spatially co-registered using a mutual information algorithm. The inverse transformation matrix was applied to the reference control points and compared with the corresponding MRI-identified ones to evaluate the overall spatial accuracy of the MRI/CT dataset. The mean image distortion correction technique was implemented, and resulting MRI-corrected control points were compared against the corresponding reference ones. For the scanning parameters used, increased MRI distortion (>1mm) was detected at areas distant from the MRI isocenter (>5cm), while median radial distortion was 0.76mm. Detected offsets were slightly higher for the MRI/CT dataset (0.92mm median distortion). The mean image distortion correction improves geometric accuracy, but residual distortion cannot be considered negligible (0.51mm median distortion). For all three datasets studied, a statistically significant positive correlation between detected spatial offsets and their distance from the MRI isocenter was revealed. This work contributes towards the wider adoption of 3T imaging in SRS/SRT treatment planning. The presented methodology can be employed in commissioning and quality assurance programmes of corresponding treatment workflows.
Collapse
Affiliation(s)
- Stefanos Theocharis
- Medical Physics Laboratory, Medical School, National and Kapodistrian University of Athens, Athens, Greece
| | - Eleftherios P. Pappas
- Medical Physics Laboratory, Medical School, National and Kapodistrian University of Athens, Athens, Greece
| | - Ioannis Seimenis
- Medical Physics Laboratory, Medical School, National and Kapodistrian University of Athens, Athens, Greece
| | - Panagiotis Kouris
- Medical Physics Laboratory, Medical School, National and Kapodistrian University of Athens, Athens, Greece
| | - Dimitrios Dellios
- Medical Physics Laboratory, Medical School, National and Kapodistrian University of Athens, Athens, Greece
| | - Georgios Kollias
- Medical Physics and Gamma Knife Department, Hygeia Hospital, Marousi, Greece
| | - Pantelis Karaiskos
- Medical Physics Laboratory, Medical School, National and Kapodistrian University of Athens, Athens, Greece
- * E-mail:
| |
Collapse
|
4
|
Dellios D, Pappas EP, Seimenis I, Paraskevopoulou C, Lampropoulos KI, Lymperopoulou G, Karaiskos P. Evaluation of patient-specific MR distortion correction schemes for improved target localization accuracy in SRS. Med Phys 2020; 48:1661-1672. [PMID: 33230923 DOI: 10.1002/mp.14615] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2020] [Revised: 10/16/2020] [Accepted: 11/16/2020] [Indexed: 11/08/2022] Open
Abstract
PURPOSE This work aims at promoting target localization accuracy in cranial stereotactic radiosurgery (SRS) applications by focusing on the correction of sequence-dependent (also patient induced) magnetic resonance (MR) distortions at the lesion locations. A phantom-based quality assurance (QA) methodology was developed and implemented for the evaluation of three distortion correction techniques. The same approach was also adapted to cranial MR images used for SRS treatment planning purposes in single or multiple brain metastases cases. METHODS A three-dimensional (3D)-printed head phantom was filled with a 3D polymer gel dosimeter. Following treatment planning and dose delivery, volumes of radiation-induced polymerization served as hypothetical lesions, offering adequate MR contrast with respect to the surrounding unirradiated areas. T1-weighted (T1w) MR imaging was performed at 1.5 T using the clinical scanning protocol for SRS. Additional images were acquired to implement three distortion correction methods; the field mapping (FM), mean image (MI) and signal integration (SI) techniques. Reference lesion locations were calculated as the averaged centroid positions of each target identified in the forward and reverse read gradient polarity MRI scans. The same techniques and workflows were implemented for the correction of contrast-enhanced T1w MR images of 10 patients with a total of 27 brain metastases. RESULTS All methods employed in the phantom study diminished spatial distortion. Median and maximum distortion magnitude decreased from 0.7 mm (2.10 ppm) and 0.8 mm (2.36 ppm), respectively, to <0.2 mm (0.61 ppm) at all target locations, using any of the three techniques. Image quality of the corrected images was acceptable, while contrast-to-noise ratio slightly increased. Results of the patient study were in accordance with the findings of the phantom study. Residual distortion in corrected patient images was found to be <0.3 mm in the vast majority of targets. Overall, the MI approach appears to be the most efficient correction method from the three investigated. CONCLUSIONS In cranial SRS applications, patient-specific distortion correction at the target location(s) is feasible and effective, despite the expense of longer imaging time since additional MRI scan(s) need to be performed. A phantom-based QA methodology was developed and presented to reassure efficient implementation of correction techniques for sequence-dependent spatial distortion.
Collapse
Affiliation(s)
- Dimitrios Dellios
- Medical Physics Laboratory, Medical School, National and Kapodistrian University of Athens, Athens, 115 27, Greece
| | - Eleftherios P Pappas
- Medical Physics Laboratory, Medical School, National and Kapodistrian University of Athens, Athens, 115 27, Greece
| | - Ioannis Seimenis
- Medical Physics Laboratory, Medical School, National and Kapodistrian University of Athens, Athens, 115 27, Greece
| | | | - Kostas I Lampropoulos
- Medical Physics and Gamma Knife Department, Hygeia Hospital, Marousi, 151 23, Greece
| | - Georgia Lymperopoulou
- 1st Department of Radiology, Medical School, National and Kapodistrian University of Athens, Athens, 115 28, Greece
| | - Pantelis Karaiskos
- Medical Physics Laboratory, Medical School, National and Kapodistrian University of Athens, Athens, 115 27, Greece.,Medical Physics and Gamma Knife Department, Hygeia Hospital, Marousi, 151 23, Greece
| |
Collapse
|
5
|
Neppl S, Landry G, Kurz C, Hansen DC, Hoyle B, Stöcklein S, Seidensticker M, Weller J, Belka C, Parodi K, Kamp F. Evaluation of proton and photon dose distributions recalculated on 2D and 3D Unet-generated pseudoCTs from T1-weighted MR head scans. Acta Oncol 2019; 58:1429-1434. [PMID: 31271093 DOI: 10.1080/0284186x.2019.1630754] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
Introduction: The recent developments of magnetic resonance (MR) based adaptive strategies for photon and, potentially for proton therapy, require a fast and reliable conversion of MR images to X-ray computed tomography (CT) values. CT values are needed for photon and proton dose calculation. The improvement of conversion results employing a 3D deep learning approach is evaluated. Material and methods: A database of 89 T1-weighted MR head scans with about 100 slices each, including rigidly registered CTs, was created. Twenty-eight validation patients were randomly sampled, and four patients were selected for application. The remaining patients were used to train a 2D and a 3D U-shaped convolutional neural network (Unet). A stack size of 32 slices was used for 3D training. For all application cases, volumetric modulated arc therapy photon and single-field uniform dose pencil-beam scanning proton plans at four different gantry angles were optimized for a generic target on the CT and recalculated on 2D and 3D Unet-based pseudoCTs. Mean (absolute) error (MAE/ME) and a gradient sharpness estimate were used to quantify the image quality. Three-dimensional gamma and dose difference analyses were performed for photon (gamma criteria: 1%, 1 mm) and proton dose distributions (gamma criteria: 2%, 2 mm). Range (80% fall off) differences for beam's eye view profiles were evaluated for protons. Results: Training 36 h for 1000 epochs in 3D (6 h for 200 epochs in 2D) yielded a maximum MAE of 147 HU (135 HU) for the application patients. Except for one patient gamma pass rates for photon and proton dose distributions were above 96% for both Unets. Slice discontinuities were reduced for 3D training at the cost of sharpness. Conclusions: Image analysis revealed a slight advantage of 2D Unets compared to 3D Unets. Similar dose calculation performance was reached for the 2D and 3D network.
Collapse
Affiliation(s)
- Sebastian Neppl
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München (LMU Munich), Garching bei München, Germany
| | - Guillaume Landry
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München (LMU Munich), Garching bei München, 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 (LMU Munich), Garching bei München, Germany
| | - David C. Hansen
- Department of Medical Physics, Aarhus University Hospital, Aarhus, Denmark
| | - Ben Hoyle
- University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität München (LMU Munich), Munich, Germany
| | - Sophia Stöcklein
- Department of Radiology, University Hospital, LMU Munich, Munich, Germany
| | - Max Seidensticker
- Department of Radiology, University Hospital, LMU Munich, Munich, Germany
| | - Jochen Weller
- University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität München (LMU Munich), Munich, Germany
- Optical and Interpretative Astronomy, Max Planck Institute for Extraterrestrial Physics, Garching bei München, Germany
| | - Claus Belka
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
- German Cancer Consortium (DKTK), Partner site Munich, Munich, Germany
| | - Katia Parodi
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München (LMU Munich), Garching bei München, Germany
| | - Florian Kamp
- Department of Radiation Oncology, University Hospital, LMU Munich, Munich, Germany
| |
Collapse
|
6
|
|
7
|
Lei Y, Harms J, Wang T, Liu Y, Shu HK, Jani AB, Curran WJ, Mao H, Liu T, Yang X. MRI-only based synthetic CT generation using dense cycle consistent generative adversarial networks. Med Phys 2019; 46:3565-3581. [PMID: 31112304 DOI: 10.1002/mp.13617] [Citation(s) in RCA: 146] [Impact Index Per Article: 29.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2019] [Revised: 05/14/2019] [Accepted: 05/14/2019] [Indexed: 02/07/2023] Open
Abstract
PURPOSE Automated synthetic computed tomography (sCT) generation based on magnetic resonance imaging (MRI) images would allow for MRI-only based treatment planning in radiation therapy, eliminating the need for CT simulation and simplifying the patient treatment workflow. In this work, the authors propose a novel method for generation of sCT based on dense cycle-consistent generative adversarial networks (cycle GAN), a deep-learning based model that trains two transformation mappings (MRI to CT and CT to MRI) simultaneously. METHODS AND MATERIALS The cycle GAN-based model was developed to generate sCT images in a patch-based framework. Cycle GAN was applied to this problem because it includes an inverse transformation from CT to MRI, which helps constrain the model to learn a one-to-one mapping. Dense block-based networks were used to construct generator of cycle GAN. The network weights and variables were optimized via a gradient difference (GD) loss and a novel distance loss metric between sCT and original CT. RESULTS Leave-one-out cross-validation was performed to validate the proposed model. The mean absolute error (MAE), peak signal-to-noise ratio (PSNR), and normalized cross correlation (NCC) indexes were used to quantify the differences between the sCT and original planning CT images. For the proposed method, the mean MAE between sCT and CT were 55.7 Hounsfield units (HU) for 24 brain cancer patients and 50.8 HU for 20 prostate cancer patients. The mean PSNR and NCC were 26.6 dB and 0.963 in the brain cases, and 24.5 dB and 0.929 in the pelvis. CONCLUSION We developed and validated a novel learning-based approach to generate CT images from routine MRIs based on dense cycle GAN model to effectively capture the relationship between the CT and MRIs. The proposed method can generate robust, high-quality sCT in minutes. The proposed method offers strong potential for supporting near real-time MRI-only treatment planning in the brain and pelvis.
Collapse
Affiliation(s)
- Yang Lei
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Joseph Harms
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Tonghe Wang
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Yingzi Liu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Hui-Kuo Shu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Ashesh B Jani
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Walter J Curran
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Hui Mao
- Department of Radiology and Imaging Sciences and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Tian Liu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| | - Xiaofeng Yang
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
| |
Collapse
|
8
|
Dai T, Li Q, Liu X, Dai Z, He P, Ma Y, Shen G, Chen W, Zhang H, Meng Q, Zhang X. Technical Note: Effect of magnetic fields on the microdosimetry of carbon-ion beams. Med Phys 2019; 46:3746-3750. [PMID: 31148177 DOI: 10.1002/mp.13633] [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/07/2019] [Revised: 05/10/2019] [Accepted: 05/24/2019] [Indexed: 11/06/2022] Open
Abstract
PURPOSE To investigate the influence of magnetic fields on the microdosimetry of carbon-ion beams and the scaling effect of tissue equivalent proportional counter (TEPC) defined as the change of energy deposition in a TEPC from that in a microscopic scale region of interest due to the presence of a magnetic field in combination with the TEPC larger physical dimensions. METHODS Geant4-based Monte Carlo simulations were conducted to calculate the microdosimetric quantities for carbon-ion beams with different initial energies (10-290 MeV/u) under magnetic fields of various strengths (0.5-3 T). The calculations were performed for a 1 μm spherical volume made of tissue, and for spherical TEPCs of 1 and 10 mm in diameter. Then, values of dose-averaged lineal energy (yD ) were acquired for the different scenarios to analyze the effect of magnetic fields on the microdosimetry of carbon-ion beams and the scaling effect of TEPC. RESULTS The yD values and lineal energy spectra in the 1 μm spherical tissue volume for the scenarios without magnetic field and with magnetic fields of different strengths and directions remained nearly the same for the various energy carbon-ion beams. However, compared with those of the 1 μm spherical tissue volume, an increase of of yD values and an obvious shift of the lineal energy spectra for the TEPCs of 1 and 10 mm in diameter under magnetic fields were found. CONCLUSIONS The application of magnetic fields under 3 T has no significant influence on the microdosimetric results of carbon-ion beams. However, there is definitely a scaling effect when using TEPC for microdosimetric study, which alters the reading of TEPC in the presence of magnetic fields. Novel methods to correct the reading of TEPC or scaling effect-resistant microdosimetric measurement detectors are urgently needed to perform experimental microdosimetric studies under magnetic fields.
Collapse
Affiliation(s)
- Tianyuan Dai
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Qiang Li
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Xinguo Liu
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Zhongying Dai
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Pengbo He
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Yuanyuan Ma
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Guosheng Shen
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Weiqiang Chen
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Hui Zhang
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Qianqian Meng
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Xiaofang Zhang
- Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou, 73000, China.,Key Laboratory of Heavy Ion Radiation Biology and Medicine of Chinese Academy of Science, Lanzhou, 730000, China.,Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, Gansu Province, Lanzhou, 730000, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| |
Collapse
|
9
|
Lei Y, Harms J, Wang T, Tian S, Zhou J, Shu HK, Zhong J, Mao H, Curran WJ, Liu T, Yang X. MRI-based synthetic CT generation using semantic random forest with iterative refinement. Phys Med Biol 2019; 64:085001. [PMID: 30818292 PMCID: PMC7778365 DOI: 10.1088/1361-6560/ab0b66] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Target delineation for radiation therapy treatment planning often benefits from magnetic resonance imaging (MRI) in addition to x-ray computed tomography (CT) due to MRI's superior soft tissue contrast. MRI-based treatment planning could reduce systematic MR-CT co-registration errors, medical cost, radiation exposure, and simplify clinical workflow. However, MRI-only based treatment planning is not widely used to date because treatment-planning systems rely on the electron density information provided by CTs to calculate dose. Additionally, air and bone regions are difficult to separategiven their similar intensities in MR imaging. The purpose of this work is to develop a learning-based method to generate patient-specific synthetic CT (sCT) from a routine anatomical MRI for use in MRI-only radiotherapy treatment planning. An auto-context model with patch-based anatomical features was integrated into a classification random forest to generate and improve semantic information. The semantic information along with anatomical features was then used to train a series of regression random forests based on the auto-context model. After training, the sCT of a new MRI can be generated by feeding anatomical features extracted from the MRI into the well-trained classification and regression random forests. The proposed algorithm was evaluated using 14 patient datasets withT1-weighted MR and corresponding CT images of the brain. The mean absolute error (MAE), peak signal-to-noise ratio (PSNR), and normalized cross correlation (NCC) were 57.45 ± 8.45 HU, 28.33 ± 1.68 dB, and 0.97 ± 0.01. We also compared the difference between dose maps calculated on the sCT and those on the original CT, using the same plan parameters. The average DVH differences among all patients are less than 0.2 Gy for PTVs, and less than 0.02 Gy for OARs. The sCT generation by the proposed method allows for dose calculation based MR imaging alone, and may be a useful tool for MRI-based radiation treatment planning.
Collapse
Affiliation(s)
- Yang Lei
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| | - Joseph Harms
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| | - Tonghe Wang
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| | - Sibo Tian
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| | - Jun Zhou
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| | - Hui-Kuo Shu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| | - Jim Zhong
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| | - Hui Mao
- Department of Radiology and Imaging Sciences and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| | - Walter J Curran
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| | - Tian Liu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| | - Xiaofeng Yang
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, United States of America
| |
Collapse
|
10
|
Yang X, Wang T, Lei Y, Higgins K, Liu T, Shim H, Curran WJ, Mao H, Nye JA. MRI-based attenuation correction for brain PET/MRI based on anatomic signature and machine learning. Phys Med Biol 2019; 64:025001. [PMID: 30524027 PMCID: PMC7773209 DOI: 10.1088/1361-6560/aaf5e0] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Deriving accurate attenuation maps for PET/MRI remains a challenging problem because MRI voxel intensities are not related to properties of photon attenuation and bone/air interfaces have similarly low signal. This work presents a learning-based method to derive patient-specific computed tomography (CT) maps from routine T1-weighted MRI in their native space for attenuation correction of brain PET. We developed a machine-learning-based method using a sequence of alternating random forests under the framework of an iterative refinement model. Anatomical feature selection is included in both training and predication stages to achieve optimal performance. To evaluate its accuracy, we retrospectively investigated 17 patients, each of which has been scanned by PET/CT and MR for brain. The PET images were corrected for attenuation on CT images as ground truth, as well as on pseudo CT (PCT) images generated from MR images. The PCT images showed mean average error of 66.1 ± 8.5 HU, average correlation coefficient of 0.974 ± 0.018 and average Dice similarity coefficient (DSC) larger than 0.85 for air, bone and soft tissue. The side-by-side image comparisons and joint histograms demonstrated very good agreement of PET images after correction by PCT and CT. The mean differences of voxel values in selected VOIs were less than 4%, the mean absolute difference of all active area is around 2.5%, and the mean linear correlation coefficient is 0.989 ± 0.017 between PET images corrected by CT and PCT. This work demonstrates a novel learning-based approach to automatically generate CT images from routine T1-weighted MR images based on a random forest regression with patch-based anatomical signatures to effectively capture the relationship between the CT and MR images. Reconstructed PET images using the PCT exhibit errors well below accepted test/retest reliability of PET/CT indicating high quantitative equivalence.
Collapse
Affiliation(s)
- Xiaofeng Yang
- Department of Radiation Oncology, Emory University, Atlanta, GA, United States of America
- Winship Cancer Institute, Emory University, Atlanta, GA, United States of America
| | - Tonghe Wang
- Department of Radiation Oncology, Emory University, Atlanta, GA, United States of America
| | - Yang Lei
- Department of Radiation Oncology, Emory University, Atlanta, GA, United States of America
| | - Kristin Higgins
- Department of Radiation Oncology, Emory University, Atlanta, GA, United States of America
- Winship Cancer Institute, Emory University, Atlanta, GA, United States of America
| | - Tian Liu
- Department of Radiation Oncology, Emory University, Atlanta, GA, United States of America
- Winship Cancer Institute, Emory University, Atlanta, GA, United States of America
| | - Hyunsuk Shim
- Department of Radiation Oncology, Emory University, Atlanta, GA, United States of America
- Department of Radiology and Imaging Sciences, Emory University, Atlanta, GA, United States of America
- Winship Cancer Institute, Emory University, Atlanta, GA, United States of America
| | - Walter J Curran
- Department of Radiation Oncology, Emory University, Atlanta, GA, United States of America
- Winship Cancer Institute, Emory University, Atlanta, GA, United States of America
| | - Hui Mao
- Department of Radiology and Imaging Sciences, Emory University, Atlanta, GA, United States of America
- Winship Cancer Institute, Emory University, Atlanta, GA, United States of America
| | - Jonathon A Nye
- Department of Radiology and Imaging Sciences, Emory University, Atlanta, GA, United States of America
| |
Collapse
|
11
|
Adjeiwaah M, Bylund M, Lundman JA, Söderström K, Zackrisson B, Jonsson JH, Garpebring A, Nyholm T. Dosimetric Impact of MRI Distortions: A Study on Head and Neck Cancers. Int J Radiat Oncol Biol Phys 2018; 103:994-1003. [PMID: 30496879 DOI: 10.1016/j.ijrobp.2018.11.037] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2018] [Revised: 11/13/2018] [Accepted: 11/19/2018] [Indexed: 10/27/2022]
Abstract
PURPOSE To evaluate the effect of magnetic resonance (MR) imaging (MRI) geometric distortions on head and neck radiation therapy treatment planning (RTP) for an MRI-only RTP. We also assessed the potential benefits of patient-specific shimming to reduce the magnitude of MR distortions for a 3-T scanner. METHODS AND MATERIALS Using an in-house Matlab algorithm, shimming within entire imaging volumes and user-defined regions of interest were simulated. We deformed 21 patient computed tomography (CT) images with MR distortion fields (gradient nonlinearity and patient-induced susceptibility effects) to create distorted CT (dCT) images using bandwidths of 122 and 488 Hz/mm at 3 T. Field parameters from volumetric modulated arc therapy plans initially optimized on dCT data sets were transferred to CT data to compute a new plan. Both plans were compared to determine the impact of distortions on dose distributions. RESULTS Shimming across entire patient volumes decreased the percentage of voxels with distortions of more than 2 mm from 15.4% to 2.0%. Using the user-defined region of interest (ROI) shimming strategy, (here the Planning target volume (PTV) was the chosen ROI volume) led to increased geometric for volumes outside the PTV, as such voxels within the spinal cord with geometric shifts above 2 mm increased from 11.5% to 32.3%. The worst phantom-measured residual system distortions after 3-dimensional gradient nonlinearity correction within a radial distance of 200 mm from the isocenter was 2.17 mm. For all patients, voxels with distortion shifts of more than 2 mm resulting from patient-induced susceptibility effects were 15.4% and 0.0% using bandwidths of 122 Hz/mm and 488 Hz/mm at 3 T. Dose differences between dCT and CT treatment plans in D50 at the planning target volume were 0.4% ± 0.6% and 0.3% ± 0.5% at 122 and 488 Hz/mm, respectively. CONCLUSIONS The overall effect of MRI geometric distortions on data used for RTP was minimal. Shimming over entire imaging volumes decreased distortions, but user-defined subvolume shimming introduced significant errors in nearby organs and should probably be avoided.
Collapse
Affiliation(s)
- Mary Adjeiwaah
- Department of Radiation Sciences, Umeå University, Umeå, Sweden.
| | - Mikael Bylund
- Department of Radiation Sciences, Umeå University, Umeå, Sweden
| | - Josef A Lundman
- Department of Radiation Sciences, Umeå University, Umeå, Sweden
| | | | | | | | | | - Tufve Nyholm
- Department of Radiation Sciences, Umeå University, Umeå, Sweden
| |
Collapse
|
12
|
Lei Y, Jeong JJ, Wang T, Shu HK, Patel P, Tian S, Liu T, Shim H, Mao H, Jani AB, Curran WJ, Yang X. MRI-based pseudo CT synthesis using anatomical signature and alternating random forest with iterative refinement model. J Med Imaging (Bellingham) 2018; 5:043504. [PMID: 30840748 PMCID: PMC6280993 DOI: 10.1117/1.jmi.5.4.043504] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2018] [Accepted: 11/12/2018] [Indexed: 12/20/2022] Open
Abstract
We develop a learning-based method to generate patient-specific pseudo computed tomography (CT) from routinely acquired magnetic resonance imaging (MRI) for potential MRI-based radiotherapy treatment planning. The proposed pseudo CT (PCT) synthesis method consists of a training stage and a synthesizing stage. During the training stage, patch-based features are extracted from MRIs. Using a feature selection, the most informative features are identified as an anatomical signature to train a sequence of alternating random forests based on an iterative refinement model. During the synthesizing stage, we feed the anatomical signatures extracted from an MRI into the sequence of well-trained forests for a PCT synthesis. Our PCT was compared with original CT (ground truth) to quantitatively assess the synthesis accuracy. The mean absolute error, peak signal-to-noise ratio, and normalized cross-correlation indices were 60.87 ± 15.10 HU , 24.63 ± 1.73 dB , and 0.954 ± 0.013 for 14 patients' brain data and 29.86 ± 10.4 HU , 34.18 ± 3.31 dB , and 0.980 ± 0.025 for 12 patients' pelvic data, respectively. We have investigated a learning-based approach to synthesize CTs from routine MRIs and demonstrated its feasibility and reliability. The proposed PCT synthesis technique can be a useful tool for MRI-based radiation treatment planning.
Collapse
Affiliation(s)
- Yang Lei
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
| | - Jiwoong Jason Jeong
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
| | - Tonghe Wang
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
| | - Hui-Kuo Shu
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
| | - Pretesh Patel
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
| | - Sibo Tian
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
| | - Tian Liu
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
| | - Hyunsuk Shim
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
- Emory University, Winship Cancer Institute, Department of Radiology and Imaging Sciences, Atlanta, Georgia, United States
| | - Hui Mao
- Emory University, Winship Cancer Institute, Department of Radiology and Imaging Sciences, Atlanta, Georgia, United States
| | - Ashesh B. Jani
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
| | - Walter J. Curran
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
| | - Xiaofeng Yang
- Emory University, Winship Cancer Institute, Department of Radiation Oncology, Atlanta, Georgia, United States
| |
Collapse
|
13
|
Wang T, Manohar N, Lei Y, Dhabaan A, Shu HK, Liu T, Curran WJ, Yang X. MRI-based treatment planning for brain stereotactic radiosurgery: Dosimetric validation of a learning-based pseudo-CT generation method. Med Dosim 2018; 44:199-204. [PMID: 30115539 DOI: 10.1016/j.meddos.2018.06.008] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2017] [Revised: 03/19/2018] [Accepted: 06/22/2018] [Indexed: 01/23/2023]
Abstract
Magnetic resonance imaging (MRI)-only radiotherapy treatment planning is attractive since MRI provides superior soft tissue contrast without ionizing radiation compared with computed tomography (CT). However, it requires the generation of pseudo CT from MRI images for patient setup and dose calculation. Our machine-learning-based method to generate pseudo CT images has been shown to provide pseudo CT images with excellent image quality, while its dose calculation accuracy remains an open question. In this study, we aim to investigate the accuracy of dose calculation in brain frameless stereotactic radiosurgery (SRS) using pseudo CT images which are generated from MRI images using the machine learning-based method developed by our group. We retrospectively investigated a total of 19 treatment plans from 14 patients, each of whom has CT simulation and MRI images acquired during pretreatment. The dose distributions of the same treatment plans were calculated on original CT simulation images as ground truth, as well as on pseudo CT images generated from MRI images. Clinically-relevant DVH metrics and gamma analysis were extracted from both ground truth and pseudo CT results for comparison and evaluation. The side-by-side comparisons on image quality and dose distributions demonstrated very good agreement of image contrast and calculated dose between pseudo CT and original CT. The average differences in Dose-volume histogram (DVH) metrics for Planning target volume (PTVs) were less than 0.6%, and no differences in those for organs at risk at a significance level of 0.05. The average pass rate of gamma analysis was 99%. These quantitative results strongly indicate that the pseudo CT images created from MRI images using our proposed machine learning method are accurate enough to replace current CT simulation images for dose calculation in brain SRS treatment. This study also demonstrates the great potential for MRI to completely replace CT scans in the process of simulation and treatment planning.
Collapse
Affiliation(s)
- Tonghe Wang
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322
| | - Nivedh Manohar
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322
| | - Yang Lei
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322
| | - Anees Dhabaan
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322
| | - Hui-Kuo Shu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322
| | - Tian Liu
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322
| | - Walter J Curran
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322
| | - Xiaofeng Yang
- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA 30322.
| |
Collapse
|
14
|
Kashani R, Olsen JR. Magnetic Resonance Imaging for Target Delineation and Daily Treatment Modification. Semin Radiat Oncol 2018; 28:178-184. [DOI: 10.1016/j.semradonc.2018.02.002] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
|
15
|
Pappas EP, Seimenis I, Dellios D, Kollias G, Lampropoulos KI, Karaiskos P. Assessment of sequence dependent geometric distortion in contrast-enhanced MR images employed in stereotactic radiosurgery treatment planning. ACTA ACUST UNITED AC 2018; 63:135006. [DOI: 10.1088/1361-6560/aac7bf] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
16
|
Yan Y, Yang J, Beddar S, Ibbott G, Wen Z, Court LE, Hwang KP, Kadbi M, Krishnan S, Fuller CD, Frank SJ, Yang J, Balter P, Kudchadker RJ, Wang J. A methodology to investigate the impact of image distortions on the radiation dose when using magnetic resonance images for planning. Phys Med Biol 2018. [PMID: 29528037 DOI: 10.1088/1361-6560/aab5c3] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
We developed a novel technique to study the impact of geometric distortion of magnetic resonance imaging (MRI) on intensity-modulated radiation therapy treatment planning. The measured 3D datasets of residual geometric distortion (a 1.5 T MRI component of an MRI linear accelerator system) was fitted with a second-order polynomial model to map the spatial dependence of geometric distortions. Then the geometric distortion model was applied to computed tomography (CT) image and structure data to simulate the distortion of MRI data and structures. Fourteen CT-based treatment plans were selected from patients treated for gastrointestinal, genitourinary, thoracic, head and neck, or spinal tumors. Plans based on the distorted CT and structure data were generated (as the distorted plans). Dose deviations of the distorted plans were calculated and compared with the original plans to study the dosimetric impact of MRI distortion. The MRI geometric distortion led to notable dose deviations in five of the 14 patients, causing loss of target coverage of up to 3.68% and dose deviations to organs at risk in three patients, increasing the mean dose to the chest wall by up to 6.19 Gy in a gastrointestinal patient, and increases the maximum dose to the lung by 5.17 Gy in a thoracic patient.
Collapse
Affiliation(s)
- Yue Yan
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, United States of America. Joint first authors
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
17
|
Emmerich J, Laun FB, Pfaffenberger A, Schilling R, Denoix M, Maier F, Sterzing F, Bostel T, Straub S. Technical Note: On the size of susceptibility-induced MR image distortions in prostate and cervix in the context of MR-guided radiation therapy. Med Phys 2018; 45:1586-1593. [DOI: 10.1002/mp.12785] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2017] [Revised: 12/11/2017] [Accepted: 01/14/2018] [Indexed: 11/06/2022] Open
Affiliation(s)
- Julian Emmerich
- Department of Medical Physics in Radiology; German Cancer Research Center (DKFZ); Heidelberg Germany
| | - Frederik B. Laun
- Department of Medical Physics in Radiology; German Cancer Research Center (DKFZ); Heidelberg Germany
- Institute of Radiology; University Hospital Erlangen; Erlangen Germany
| | - Asja Pfaffenberger
- Department of Medical Physics in Radiation Oncology; German Cancer Research Center (DKFZ); Heidelberg Germany
| | | | - Michael Denoix
- Department of Medical Physics in Radiology; German Cancer Research Center (DKFZ); Heidelberg Germany
| | - Florian Maier
- Department of Medical Physics in Radiology; German Cancer Research Center (DKFZ); Heidelberg Germany
| | - Florian Sterzing
- Clinical Cooperation Unit Radiation Oncology; German Cancer Research Center (DKFZ); Heidelberg Germany
- Department of Radiation Oncology; University Hospital Heidelberg; Heidelberg Germany
- National Center for Research in Radiation Oncology; Heidelberg Institute for Radiation Oncology (HIRO); Heidelberg Germany
| | - Tilman Bostel
- Clinical Cooperation Unit Radiation Oncology; German Cancer Research Center (DKFZ); Heidelberg Germany
- Department of Radiation Oncology; University Hospital Heidelberg; Heidelberg Germany
- National Center for Research in Radiation Oncology; Heidelberg Institute for Radiation Oncology (HIRO); Heidelberg Germany
| | - Sina Straub
- Department of Medical Physics in Radiology; German Cancer Research Center (DKFZ); Heidelberg Germany
| |
Collapse
|
18
|
Koivula L, Wee L, Korhonen J. Feasibility of MRI-only treatment planning for proton therapy in brain and prostate cancers: Dose calculation accuracy in substitute CT images. Med Phys 2017; 43:4634. [PMID: 27487880 DOI: 10.1118/1.4958677] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
PURPOSE Magnetic resonance imaging (MRI) is increasingly used for radiotherapy target delineation, image guidance, and treatment response monitoring. Recent studies have shown that an entire external x-ray radiotherapy treatment planning (RTP) workflow for brain tumor or prostate cancer patients based only on MRI reference images is feasible. This study aims to show that a MRI-only based RTP workflow is also feasible for proton beam therapy plans generated in MRI-based substitute computed tomography (sCT) images of the head and the pelvis. METHODS The sCTs were constructed for ten prostate cancer and ten brain tumor patients primarily by transforming the intensity values of in-phase MR images to Hounsfield units (HUs) with a dual model HU conversion technique to enable heterogeneous tissue representation. HU conversion models for the pelvis were adopted from previous studies, further extended in this study also for head MRI by generating anatomical site-specific conversion models (a new training data set of ten other brain patients). This study also evaluated two other types of simplified sCT: dual bulk density (for bone and water) and homogeneous (water only). For every clinical case, intensity modulated proton therapy (IMPT) plans robustly optimized in standard planning CTs were calculated in sCT for evaluation, and vice versa. Overall dose agreement was evaluated using dose-volume histogram parameters and 3D gamma criteria. RESULTS In heterogeneous sCTs, the mean absolute errors in HUs were 34 (soft tissues: 13, bones: 92) and 42 (soft tissues: 9, bones: 97) in the head and in the pelvis, respectively. The maximum absolute dose differences relative to CT in the brain tumor clinical target volume (CTV) were 1.4% for heterogeneous sCT, 1.8% for dual bulk sCT, and 8.9% for homogenous sCT. The corresponding maximum differences in the prostate CTV were 0.6%, 1.2%, and 3.6%, respectively. The percentages of dose points in the head and pelvis passing 1% and 1 mm gamma index criteria were over 91%, 85%, and 38% with heterogeneous, dual bulk, and homogeneous sCTs, respectively. There were no significant changes to gamma index pass rates for IMPT plans first optimized in CT and then calculated in heterogeneous sCT versus IMPT plans first optimized in heterogeneous sCT and then calculated on standard CT. CONCLUSIONS This study demonstrates that proton therapy dose calculations on heterogeneous sCTs are in good agreement with plans generated with standard planning CT. An MRI-only based RTP workflow is feasible in IMPT for brain tumors and prostate cancers.
Collapse
Affiliation(s)
- Lauri Koivula
- Department of Radiation Oncology, Comprehensive Cancer Center, Helsinki University Central Hospital, P.O. Box 180, Helsinki 00029 HUS, Finland and Department of Medical Physics, Oncology Services, Vejle Hospital, Kabbeltoft 25, Vejle DK-7100, Denmark
| | - Leonard Wee
- Department of Medical Physics, Oncology Services, Vejle Hospital, Kabbeltoft 25, Vejle DK-7100, Denmark and Danish Colorectal Cancer Centre South, Vejle Hospital, Kabbeltoft 25, Vejle DK-7100, Denmark
| | - Juha Korhonen
- Department of Radiation Oncology, Comprehensive Cancer Center, Helsinki University Central Hospital, P.O. Box 180, Helsinki 00029 HUS, Finland; Danish Colorectal Cancer Centre South, Vejle Hospital, Kabbeltoft 25, Vejle DK-7100, Denmark; and Department of Radiology, Helsinki University Central Hospital, P.O. Box 180, Helsinki 00029 HUS, Finland
| |
Collapse
|
19
|
Ito K, Kadoya N, Nakajima Y, Saito M, Sato K, Nagasaka T, Yamanaka K, Dobashi S, Takeda K, Matsushita H, Jingu K. Feasibility of a Direct-Conversion Method from Magnetic Susceptibility to Relative Electron Density for Radiation Therapy Treatment Planning. ACTA ACUST UNITED AC 2017. [DOI: 10.4236/ijmpcero.2017.63023] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
20
|
Pappas EP, Seimenis I, Moutsatsos A, Georgiou E, Nomikos P, Karaiskos P. Characterization of system-related geometric distortions in MR images employed in Gamma Knife radiosurgery applications. Phys Med Biol 2016; 61:6993-7011. [DOI: 10.1088/0031-9155/61/19/6993] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
21
|
Weygand J, Fuller CD, Ibbott GS, Mohamed ASR, Ding Y, Yang J, Hwang KP, Wang J. Spatial Precision in Magnetic Resonance Imaging-Guided Radiation Therapy: The Role of Geometric Distortion. Int J Radiat Oncol Biol Phys 2016; 95:1304-16. [PMID: 27354136 DOI: 10.1016/j.ijrobp.2016.02.059] [Citation(s) in RCA: 113] [Impact Index Per Article: 14.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2015] [Revised: 02/05/2016] [Accepted: 02/25/2016] [Indexed: 12/11/2022]
Abstract
Because magnetic resonance imaging-guided radiation therapy (MRIgRT) offers exquisite soft tissue contrast and the ability to image tissues in arbitrary planes, the interest in this technology has increased dramatically in recent years. However, intrinsic geometric distortion stemming from both the system hardware and the magnetic properties of the patient affects MR images and compromises the spatial integrity of MRI-based radiation treatment planning, given that for real-time MRIgRT, precision within 2 mm is desired. In this article, we discuss the causes of geometric distortion, describe some well-known distortion correction algorithms, and review geometric distortion measurements from 12 studies, while taking into account relevant imaging parameters. Eleven of the studies reported phantom measurements quantifying system-dependent geometric distortion, while 2 studies reported simulation data quantifying magnetic susceptibility-induced geometric distortion. Of the 11 studies investigating system-dependent geometric distortion, 5 reported maximum measurements less than 2 mm. The simulation studies demonstrated that magnetic susceptibility-induced distortion is typically smaller than system-dependent distortion but still nonnegligible, with maximum distortion ranging from 2.1 to 2.6 mm at a field strength of 1.5 T. As expected, anatomic landmarks containing interfaces between air and soft tissue had the largest distortions. The evidence indicates that geometric distortion reduces the spatial integrity of MRI-based radiation treatment planning and likely diminishes the efficacy of MRIgRT. Better phantom measurement techniques and more effective distortion correction algorithms are needed to achieve the desired spatial precision.
Collapse
Affiliation(s)
- Joseph Weygand
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas; The University of Texas Graduate School of Biomedical Sciences at Houston, Houston, Texas.
| | - Clifton David Fuller
- The University of Texas Graduate School of Biomedical Sciences at Houston, Houston, Texas; Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Geoffrey S Ibbott
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas; The University of Texas Graduate School of Biomedical Sciences at Houston, Houston, Texas
| | - Abdallah S R Mohamed
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas; Department of Clinical Oncology and Nuclear Medicine, Alexandria University, Alexandria, Egypt
| | - Yao Ding
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Jinzhong Yang
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas; The University of Texas Graduate School of Biomedical Sciences at Houston, Houston, Texas
| | - Ken-Pin Hwang
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Jihong Wang
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas; The University of Texas Graduate School of Biomedical Sciences at Houston, Houston, Texas
| |
Collapse
|
22
|
Walker A, Liney G, Holloway L, Dowling J, Rivest-Henault D, Metcalfe P. Continuous table acquisition MRI for radiotherapy treatment planning: distortion assessment with a new extended 3D volumetric phantom. Med Phys 2015; 42:1982-91. [PMID: 25832089 DOI: 10.1118/1.4915920] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE Accurate geometry is required for radiotherapy treatment planning (RTP). When considering the use of magnetic resonance imaging (MRI) for RTP, geometric distortions observed in the acquired images should be considered. While scanner technology and vendor supplied correction algorithms provide some correction, large distortions are still present in images, even when considering considerably smaller scan lengths than those typically acquired with CT in conventional RTP. This study investigates MRI acquisition with a moving table compared with static scans for potential geometric benefits for RTP. METHODS A full field of view (FOV) phantom (diameter 500 mm; length 513 mm) was developed for measuring geometric distortions in MR images over volumes pertinent to RTP. The phantom consisted of layers of refined plastic within which vitamin E capsules were inserted. The phantom was scanned on CT to provide the geometric gold standard and on MRI, with differences in capsule location determining the distortion. MRI images were acquired with two techniques. For the first method, standard static table acquisitions were considered. Both 2D and 3D acquisition techniques were investigated. With the second technique, images were acquired with a moving table. The same sequence was acquired with a static table and then with table speeds of 1.1 mm/s and 2 mm/s. All of the MR images acquired were registered to the CT dataset using a deformable B-spline registration with the resulting deformation fields providing the distortion information for each acquisition. RESULTS MR images acquired with the moving table enabled imaging of the whole phantom length while images acquired with a static table were only able to image 50%-70% of the phantom length of 513 mm. Maximum distortion values were reduced across a larger volume when imaging with a moving table. Increased table speed resulted in a larger contribution of distortion from gradient nonlinearities in the through-plane direction and an increased blurring of capsule images, resulting in an apparent capsule volume increase by up to 170% in extreme axial FOV regions. Blurring increased with table speed and in the central regions of the phantom, geometric distortion was less for static table acquisitions compared to a table speed of 2 mm/s over the same volume. Overall, the best geometric accuracy was achieved with a table speed of 1.1 mm/s. CONCLUSIONS The phantom designed enables full FOV imaging for distortion assessment for the purposes of RTP. MRI acquisition with a moving table extends the imaging volume in the z direction with reduced distortions which could be useful particularly if considering MR-only planning. If utilizing MR images to provide additional soft tissue information to the planning CT, standard acquisition sequences over a smaller volume would avoid introducing additional blurring or distortions from the through-plane table movement.
Collapse
Affiliation(s)
- Amy Walker
- Centre for Medical Radiation Physics, University of Wollongong, Wollongong, NSW 2522, Australia and Liverpool and Macarthur Cancer Therapy Centres and Ingham Institute for Applied Medical Research, Liverpool, NSW 2170, Australia
| | - Gary Liney
- Centre for Medical Radiation Physics, University of Wollongong, Wollongong, NSW 2522, Australia; Liverpool and Macarthur Cancer Therapy Centres and Ingham Institute for Applied Medical Research, Liverpool, NSW 2170, Australia; and South West Clinical School, University of New South Wales, Sydney, NSW 2170, Australia
| | - Lois Holloway
- Centre for Medical Radiation Physics, University of Wollongong, Wollongong, NSW 2522, Australia; Liverpool and Macarthur Cancer Therapy Centres and Ingham Institute for Applied Medical Research, Liverpool, NSW 2170, Australia; South West Clinical School, University of New South Wales, Sydney, NSW 2170, Australia; and Institute of Medical Physics, School of Physics, University of Sydney, Sydney, NSW 2006, Australia
| | - Jason Dowling
- Commonwealth Scientific and Industrial Research Organisation, Australian E-Health Research Centre, Herston, QLD 4029, Australia
| | - David Rivest-Henault
- Commonwealth Scientific and Industrial Research Organisation, Australian E-Health Research Centre, Herston, QLD 4029, Australia
| | - Peter Metcalfe
- Centre for Medical Radiation Physics, University of Wollongong, Wollongong, NSW 2522, Australia and Liverpool and Macarthur Cancer Therapy Centres and Ingham Institute for Applied Medical Research, Liverpool, NSW 2170, Australia
| |
Collapse
|
23
|
Abstract
The use of magnetic resonance imaging (MRI) in radiotherapy (RT) planning is rapidly expanding. We review the wide range of image contrast mechanisms available to MRI and the way they are exploited for RT planning. However a number of challenges are also considered: the requirements that MR images are acquired in the RT treatment position, that they are geometrically accurate, that effects of patient motion during the scan are minimized, that tissue markers are clearly demonstrated, that an estimate of electron density can be obtained. These issues are discussed in detail, prior to the consideration of a number of specific clinical applications. This is followed by a brief discussion on the development of real-time MRI-guided RT.
Collapse
Affiliation(s)
- Maria A Schmidt
- Cancer Research UK Cancer Imaging Centre, Royal Marsden Hospital and the Institute of Cancer Research, Downs Road, Sutton, Surrey, SM2 5PT, UK
| | | |
Collapse
|
24
|
Tumour Movement in Proton Therapy: Solutions and Remaining Questions: A Review. Cancers (Basel) 2015; 7:1143-53. [PMID: 26132317 PMCID: PMC4586762 DOI: 10.3390/cancers7030829] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2015] [Revised: 06/10/2015] [Accepted: 06/18/2015] [Indexed: 12/25/2022] Open
Abstract
Movement of tumours, mostly by respiration, has been a major problem for treating lung cancer, liver tumours and other locations in the abdomen and thorax. Organ motion is indeed one component of geometrical uncertainties that includes delineation and target definition uncertainties, microscopic disease and setup errors. At present, minimising motion seems to be the easiest to implement in clinical practice. If combined with adaptive approaches to correct for gradual anatomical variations, it may be a practical strategy. Other approaches such as repainting and tracking could increase the accuracy of proton therapy delivery, but advanced 4D solutions are needed. Moreover, there is a need to perform clinical studies to investigate which approach is the best in a given clinical situation. The good news is that existing and emerging technology and treatment planning systems as will without doubt lead in the forthcoming future to practical solutions to tackle intra-fraction motion in proton therapy. These developments may also improve motion management in photon therapy as well.
Collapse
|
25
|
Uh J, Merchant TE, Li Y, Li X, Hua C. MRI-based treatment planning with pseudo CT generated through atlas registration. Med Phys 2014; 41:051711. [PMID: 24784377 DOI: 10.1118/1.4873315] [Citation(s) in RCA: 127] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE To evaluate the feasibility and accuracy of magnetic resonance imaging (MRI)-based treatment planning using pseudo CTs generated through atlas registration. METHODS A pseudo CT, providing electron density information for dose calculation, was generated by deforming atlas CT images previously acquired on other patients. The authors tested 4 schemes of synthesizing a pseudo CT from single or multiple deformed atlas images: use of a single arbitrarily selected atlas, arithmetic mean process using 6 atlases, and pattern recognition with Gaussian process (PRGP) using 6 or 12 atlases. The required deformation for atlas CT images was derived from a nonlinear registration of conjugated atlas MR images to that of the patient of interest. The contrasts of atlas MR images were adjusted by histogram matching to reduce the effect of different sets of acquisition parameters. For comparison, the authors also tested a simple scheme assigning the Hounsfield unit of water to the entire patient volume. All pseudo CT generating schemes were applied to 14 patients with common pediatric brain tumors. The image similarity of real patient-specific CT and pseudo CTs constructed by different schemes was compared. Differences in computation times were also calculated. The real CT in the treatment planning system was replaced with the pseudo CT, and the dose distribution was recalculated to determine the difference. RESULTS The atlas approach generally performed better than assigning a bulk CT number to the entire patient volume. Comparing atlas-based schemes, those using multiple atlases outperformed the single atlas scheme. For multiple atlas schemes, the pseudo CTs were similar to the real CTs (correlation coefficient, 0.787-0.819). The calculated dose distribution was in close agreement with the original dose. Nearly the entire patient volume (98.3%-98.7%) satisfied the criteria of chi-evaluation (<2% maximum dose and 2 mm range). The dose to 95% of the volume and the percentage of volume receiving at least 95% of the prescription dose in the planning target volume differed from the original values by less than 2% of the prescription dose (root-mean-square, RMS < 1%). The PRGP scheme did not perform better than the arithmetic mean process with the same number of atlases. Increasing the number of atlases from 6 to 12 often resulted in improvements, but statistical significance was not always found. CONCLUSIONS MRI-based treatment planning with pseudo CTs generated through atlas registration is feasible for pediatric brain tumor patients. The doses calculated from pseudo CTs agreed well with those from real CTs, showing dosimetric accuracy within 2% for the PTV when multiple atlases were used. The arithmetic mean process may be a reasonable choice over PRGP for the synthesis scheme considering performance and computational costs.
Collapse
Affiliation(s)
- Jinsoo Uh
- Department of Radiological Sciences, St. Jude Children's Research Hospital, Memphis, Tennessee 38105
| | - Thomas E Merchant
- Department of Radiological Sciences, St. Jude Children's Research Hospital, Memphis, Tennessee 38105
| | - Yimei Li
- Department of Biostatistics, St. Jude Children's Research Hospital, Memphis, Tennessee 38105
| | - Xingyu Li
- Department of Biostatistics, St. Jude Children's Research Hospital, Memphis, Tennessee 38105
| | - Chiaho Hua
- Department of Radiological Sciences, St. Jude Children's Research Hospital, Memphis, Tennessee 38105
| |
Collapse
|
26
|
Moteabbed M, Schuemann J, Paganetti H. Dosimetric feasibility of real-time MRI-guided proton therapy. Med Phys 2014; 41:111713. [PMID: 25370627 PMCID: PMC4209014 DOI: 10.1118/1.4897570] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2014] [Revised: 08/06/2014] [Accepted: 09/15/2014] [Indexed: 01/13/2023] Open
Abstract
PURPOSE Magnetic resonance imaging (MRI) is a prime candidate for image-guided radiotherapy. This study was designed to assess the feasibility of real-time MRI-guided proton therapy by quantifying the dosimetric effects induced by the magnetic field in patients' plans and identifying the associated clinical consequences. METHODS Monte Carlo dose calculation was performed for nine patients of various treatment sites (lung, liver, prostate, brain, skull-base, and spine) and tissue homogeneities, in the presence of 0.5 and 1.5 T magnetic fields. Dose volume histogram (DVH) parameters such as D95, D5, and V20 as well as equivalent uniform dose were compared for the target and organs at risk, before and after applying the magnetic field. The authors further assessed whether the plans affected by clinically relevant dose distortions could be corrected independent of the planning system. RESULTS By comparing the resulting dose distributions and analyzing the respective DVHs, it was determined that despite the observed lateral beam deflection, for magnetic fields of up to 0.5 T, neither was the target coverage jeopardized nor was the dose to the nearby organs increased in all cases except for prostate. However, for a 1.5 T magnetic field, the dose distortions were more pronounced and of clinical concern in all cases except for spine. In such circumstances, the target was severely underdosed, as indicated by a decrease in D95 of up to 41% of the prescribed dose compared to the nominal situation (no magnetic field). Sites such as liver and spine were less affected due to higher tissue homogeneity, typically smaller beam range, and the choice of beam directions. Simulations revealed that small modifications to certain plan parameters such as beam isocenter (up to 19 mm) and gantry angle (up to 10°) are sufficient to compensate for the magnetic field-induced dose disturbances. The authors' observations indicate that the degree of required corrections strongly depends on the beam range and direction relative to the magnetic field. This method was also applicable to more heterogeneous scenarios such as skull-base tumors. CONCLUSIONS This study confirmed the dosimetric feasibility of real-time MRI-guided proton therapy and delivering a clinically acceptable dose to patients with various tumor locations within magnetic fields of up to 1.5 T. This work could serve as a guide and encouragement for further efforts toward clinical implementation of hybrid MRI-proton gantry systems.
Collapse
Affiliation(s)
- M Moteabbed
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts 02114
| | - J Schuemann
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts 02114
| | - H Paganetti
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts 02114
| |
Collapse
|
27
|
Lagendijk JJW, Raaymakers BW, Van den Berg CAT, Moerland MA, Philippens ME, van Vulpen M. MR guidance in radiotherapy. Phys Med Biol 2014; 59:R349-69. [PMID: 25322150 DOI: 10.1088/0031-9155/59/21/r349] [Citation(s) in RCA: 153] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Jan J W Lagendijk
- Department of Radiotherapy, University Medical Centre Utrecht, Heidelberglaan 100, Utrecht, The Netherlands
| | | | | | | | | | | |
Collapse
|
28
|
Li W, Li Y, Zhu W, Chen X. Changes in brain functional network connectivity after stroke. Neural Regen Res 2014; 9:51-60. [PMID: 25206743 PMCID: PMC4146323 DOI: 10.4103/1673-5374.125330] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/16/2013] [Indexed: 01/15/2023] Open
Abstract
Studies have shown that functional network connection models can be used to study brain network changes in patients with schizophrenia. In this study, we inferred that these models could also be used to explore functional network connectivity changes in stroke patients. We used independent component analysis to find the motor areas of stroke patients, which is a novel way to determine these areas. In this study, we collected functional magnetic resonance imaging datasets from healthy controls and right-handed stroke patients following their first ever stroke. Using independent component analysis, six spatially independent components highly correlated to the experimental paradigm were extracted. Then, the functional network connectivity of both patients and controls was established to observe the differences between them. The results showed that there were 11 connections in the model in the stroke patients, while there were only four connections in the healthy controls. Further analysis found that some damaged connections may be compensated for by new indirect connections or circuits produced after stroke. These connections may have a direct correlation with the degree of stroke rehabilitation. Our findings suggest that functional network connectivity in stroke patients is more complex than that in hea-lthy controls, and that there is a compensation loop in the functional network following stroke. This implies that functional network reorganization plays a very important role in the process of rehabilitation after stroke.
Collapse
Affiliation(s)
- Wei Li
- Key Laboratory of Image Processing and Intelligent Control, Ministry of Education, Wuhan, Hubei Province, China ; Department of Control Science and Engineering, Huazhong University of Science and Technology, Wuhan, Hubei Province, China
| | - Yapeng Li
- Key Laboratory of Image Processing and Intelligent Control, Ministry of Education, Wuhan, Hubei Province, China ; Department of Control Science and Engineering, Huazhong University of Science and Technology, Wuhan, Hubei Province, China
| | - Wenzhen Zhu
- Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei Province, China
| | - Xi Chen
- Key Laboratory of Image Processing and Intelligent Control, Ministry of Education, Wuhan, Hubei Province, China ; Department of Control Science and Engineering, Huazhong University of Science and Technology, Wuhan, Hubei Province, China
| |
Collapse
|
29
|
Matakos A, Balter J, Cao Y. Estimation of geometrically undistorted B(0) inhomogeneity maps. Phys Med Biol 2014; 59:4945-59. [PMID: 25109506 DOI: 10.1088/0031-9155/59/17/4945] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Geometric accuracy of MRI is one of the main concerns for its use as a sole image modality in precision radiation therapy (RT) planning. In a state-of-the-art scanner, system level geometric distortions are within acceptable levels for precision RT. However, subject-induced B0 inhomogeneity may vary substantially, especially in air-tissue interfaces. Recent studies have shown distortion levels of more than 2 mm near the sinus and ear canal are possible due to subject-induced field inhomogeneity. These distortions can be corrected with the use of accurate B0 inhomogeneity field maps. Most existing methods estimate these field maps from dual gradient-echo (GRE) images acquired at two different echo-times under the assumption that the GRE images are practically undistorted. However distortion that may exist in the GRE images can result in estimated field maps that are distorted in both geometry and intensity, leading to inaccurate correction of clinical images. This work proposes a method for estimating undistorted field maps from GRE acquisitions using an iterative joint estimation technique. The proposed method yields geometrically corrected GRE images and undistorted field maps that can also be used for the correction of images acquired by other sequences. The proposed method is validated through simulation, phantom experiments and applied to patient data. Our simulation results show that our method reduces the root-mean-squared error of the estimated field map from the ground truth by ten-fold compared to the distorted field map. Both the geometric distortion and the intensity corruption (artifact) in the images caused by the B0 field inhomogeneity are corrected almost completely. Our phantom experiment showed improvement in the geometric correction of approximately 1 mm at an air-water interface using the undistorted field map compared to using a distorted field map. The proposed method for undistorted field map estimation can lead to improved geometric distortion correction at air-tissue interfaces, especially in low readout-bandwidth acquisitions, thus making them suitable for clinical use in precision RT without increasing the treatment planning margin.
Collapse
Affiliation(s)
- A Matakos
- Department of Radiation Oncology, University of Michigan, Ann Arbor, MI 48109, USA
| | | | | |
Collapse
|
30
|
Kupelian P, Sonke JJ. Magnetic Resonance–Guided Adaptive Radiotherapy: A Solution to the Future. Semin Radiat Oncol 2014; 24:227-32. [DOI: 10.1016/j.semradonc.2014.02.013] [Citation(s) in RCA: 81] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
|
31
|
Crijns S, Raaymakers B. From static to dynamic 1.5T MRI-linac prototype: impact of gantry position related magnetic field variation on image fidelity. Phys Med Biol 2014; 59:3241-7. [DOI: 10.1088/0031-9155/59/13/3241] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
32
|
MRI distortion: considerations for MRI based radiotherapy treatment planning. AUSTRALASIAN PHYSICAL & ENGINEERING SCIENCES IN MEDICINE 2014; 37:103-13. [DOI: 10.1007/s13246-014-0252-2] [Citation(s) in RCA: 72] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/16/2013] [Accepted: 02/01/2014] [Indexed: 10/25/2022]
|
33
|
Bostel T, Nicolay NH, Grossmann JG, Mohr A, Delorme S, Echner G, Häring P, Debus J, Sterzing F. MR-guidance--a clinical study to evaluate a shuttle- based MR-linac connection to provide MR-guided radiotherapy. Radiat Oncol 2014; 9:12. [PMID: 24401489 PMCID: PMC3904210 DOI: 10.1186/1748-717x-9-12] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2013] [Accepted: 12/13/2013] [Indexed: 11/10/2022] Open
Abstract
Background The purpose of this clinical study is to investigate the clinical feasibility and safety of a shuttle-based MR-linac connection to provide MR-guided radiotherapy. Methods/Design A total of 40 patients with an indication for a neoadjuvant, adjuvant or definitive radiation treatment will be recruited including tumors of the head and neck region, thorax, upper gastrointestinal tract and pelvic region. All study patients will receive standard therapy, i.e. highly conformal radiation techniques like CT-guided intensity-modulated radiotherapy (IMRT) with or without concomitant chemotherapy or other antitumor medication, and additionally daily short MR scans in treatment position with the same immobilisation equipment used for irradiation for position verification and imaging of the anatomical and functional changes during the course of radiotherapy. For daily position control, skin marks and a stereotactic frame will be used for both imaging modalities. Patient transfer between the MR device and the linear accelerator will be performed with a shuttle system which uses an air-bearing patient platform for both procedures. The daily acquired MR and CT data sets will be digitally registrated, correlated with the planning CT and compared with each other regarding translational and rotational errors. Aim of this clinical study is to establish a shuttle-based approach for realising MR-guided radiotherapy for certain clinical situations. Second objectives are to compare MR-guided radiotherapy with the gold standard of CT image guidance for quality assurance of radiotherapy, to establish an appropiate MR protocol therefore, and to assess the possibility of using MR-based image guidance not only for position verification but also for adaptive strategies in radiotherapy. Discussion Compared to CT, MRI might offer the advantage of providing IGRT without delivering an additional radiation dose to the patients and the possibility of optimisation of adaptive therapy strategies due to its superior soft tissue contrast. However, up to now, hybrid MR-linac devices are still under construction and not clinically applicable. For the near future, a shuttle-based approach would be a promising alternative for providing MR-guided radiotherapy, so that the present study was initiated to determine feasibility and safety of such an approach. Besides positioning information, daily MR data under treatment offer the possibility to assess tumor regression and functional parameters, with a potential impact not only on adaptive therapy strategies but also on early assessment of treatment response.
Collapse
Affiliation(s)
- Tilman Bostel
- Clinical Cooperation Unit Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, 69120 Heidelberg, Germany.
| | | | | | | | | | | | | | | | | |
Collapse
|
34
|
|
35
|
Wang H, Balter J, Cao Y. Patient-induced susceptibility effect on geometric distortion of clinical brain MRI for radiation treatment planning on a 3T scanner. Phys Med Biol 2013; 58:465-77. [PMID: 23302471 DOI: 10.1088/0031-9155/58/3/465] [Citation(s) in RCA: 72] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Concerns about the geometric accuracy of MRI in radiation therapy (RT) have been present since its invention. Although modern scanners typically have system levels of geometric accuracy that meet requirements of RT, subject-specific distortion is variable, and methods to in vivo assess and control patient-induced geometric distortion are not yet resolved. This study investigated the nature and magnitude of the subject-induced susceptibility effect on geometric distortions in clinical brain MRI, and tested the feasibility of in vivo quality control using field inhomogeneity mapping. For 19 consecutive patients scanned on a dedicated 3T MR scanner, B0 field inhomogeneity maps were acquired and analyzed to determine subject-induced distortions. For 3D T1 weighted images frequency-encoded with a bandwidth of 180 Hz/pixel, 86.9% of the estimated displacements were <0.5 mm, 97.4% <1 mm, and only 0.1% of displacements > 2 mm. The maximum displacement was <4 mm. The greatest distortions were observed at the interfaces with air at the sinuses. Displacements decayed to less than 1 mm over a distance of 8 mm. Metal surgical wires generated smaller distortions, with an averaged maximum displacement of 0.76 mm. Repeat acquisition of the field maps in 17 patients revealed a within-subject standard deviation of 0.25 ppm, equivalent to 0.22 mm displacement in the frequency-encoding direction in the 3D T1 weighted images. Susceptibility-induced voxel displacements in the brain are generally small, but should be monitored for precision RT. These effects are manageable at 3T and lower fields, and the methods applied can be used to monitor for potential local errors in individual patients, as well as to correct for local distortions as needed.
Collapse
Affiliation(s)
- H Wang
- Department of Radiation Oncology, University of Michigan, Ann Arbor, MI 48109, USA
| | | | | |
Collapse
|
36
|
Stanescu T, Wachowicz K, Jaffray DA. Characterization of tissue magnetic susceptibility-induced distortions for MRIgRT. Med Phys 2012; 39:7185-93. [DOI: 10.1118/1.4764481] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
37
|
Stam MK, Crijns SPM, Zonnenberg BA, Barendrecht MM, van Vulpen M, Lagendijk JJW, Raaymakers BW. Navigators for motion detection during real-time MRI-guided radiotherapy. Phys Med Biol 2012; 57:6797-805. [PMID: 23032581 DOI: 10.1088/0031-9155/57/21/6797] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
An MRI-linac system provides direct MRI feedback and with that the possibility of adapting radiation treatments to the actual tumour position. This paper addresses the use of fast 1D MRI, pencil-beam navigators, for this feedback. The accuracy of using navigators was determined on a moving phantom. The possibility of organ tracking and breath-hold monitoring based on navigator guidance was shown for the kidney. Navigators are accurate within 0.5 mm and the analysis has a minimal time lag smaller than 30 ms as shown for the phantom measurements. The correlation of 2D kidney images and navigators shows the possibility of complete organ tracking. Furthermore the breath-hold monitoring of the kidney is accurate within 1.5 mm, allowing gated radiotherapy based on navigator feedback. Navigators are a fast and precise method for monitoring and real-time tracking of anatomical landmarks. As such, they provide direct MRI feedback on anatomical changes for more precise radiation delivery.
Collapse
Affiliation(s)
- Mette K Stam
- Department of Radiotherapy, University Medical Center Utrecht, Heidelberglaan 100, 3584 CX, Utrecht, The Netherlands.
| | | | | | | | | | | | | |
Collapse
|
38
|
Crijns SPM, Bakker CJG, Seevinck PR, de Leeuw H, Lagendijk JJW, Raaymakers BW. Towards inherently distortion-free MR images for image-guided radiotherapy on an MRI accelerator. Phys Med Biol 2012; 57:1349-58. [PMID: 22349351 DOI: 10.1088/0031-9155/57/5/1349] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
In MR-guided interventions, it is mandatory to establish a solid relationship between the imaging coordinate system and world coordinates. This is particularly important in image-guided radiotherapy (IGRT) on an MRI accelerator, as the interaction of matter with γ-radiation cannot be visualized. In conventional acquisitions, off-resonance effects cause discrepancies between coordinate systems. We propose to mitigate this by using only phase encoding and to reduce the longer acquisitions by under-sampling and regularized reconstruction. To illustrate the performance of this acquisition in the presence of off-resonance phenomena, phantom and in vivo images are acquired using spin-echo (SE) and purely phase-encoded sequences. Data are retrospectively under-sampled and reconstructed iteratively. We observe accurate geometries in purely phase-encoded images for all cases, whereas SE images of the same phantoms display image distortions. Regularized reconstruction yields accurate phantom images under high acceleration factors. In vivo images were reconstructed faithfully while using acceleration factors up to 4. With the proposed technique, inherently undistorted images with one-to-one correspondence to world coordinates can be obtained. It is a valuable tool in geometry quality assurance, treatment planning and online image guidance. Under-sampled acquisition combined with regularized reconstruction can be used to accelerate the acquisition while retaining geometrical accuracy.
Collapse
Affiliation(s)
- S P M Crijns
- Department of Radiotherapy, Imaging Division, University Medical Center Utrecht, Heidelberglaan 100, 3584 CX Utrecht, The Netherlands.
| | | | | | | | | | | |
Collapse
|
39
|
Raaymakers BW, de Boer JCJ, Knox C, Crijns SPM, Smit K, Stam MK, Bosch MRVD, Kok JGM, Lagendijk JJW. Integrated megavoltage portal imaging with a 1.5 T MRI linac. Phys Med Biol 2011; 56:N207-14. [DOI: 10.1088/0031-9155/56/19/n01] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
40
|
Crijns SPM, Kok JGM, Lagendijk JJW, Raaymakers BW. Towards MRI-guided linear accelerator control: gating on an MRI accelerator. Phys Med Biol 2011; 56:4815-25. [DOI: 10.1088/0031-9155/56/15/012] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
|