MRI-only treatment planning: benefits and challenges

The Promise of Magnetic Resonance Imaging in Radiation Oncology Practice in the Management of Brain, Prostate, and GI Malignancies

Magnetic resonance imaging (MRI) has a key role to play at multiple steps of the radiotherapy (RT) treatment planning and delivery process. Development of high-precision RT techniques such as intensity-modulated RT, stereotactic ablative RT, and particle beam therapy has enabled oncologists to escalate RT dose to the target while restricting doses to organs at risk (OAR). MRI plays a critical role in target volume delineation in various disease sites, thus ensuring that these high-precision techniques can be safely implemented. Accurate identification of gross disease has also enabled selective dose escalation as a means to widen the therapeutic index. Morphological and functional MRI sequences have also facilitated an understanding of temporal changes in target volumes and OAR during a course of RT, allowing for midtreatment volumetric and biological adaptation. The latest advancement in linear accelerator technology has led to the incorporation of an MRI scanner in the treatment unit. MRI-guided RT provides the opportunity for MRI-only workflow along with online adaptation for either target or OAR or both. MRI plays a key role in post-treatment response evaluation and is an important tool for guiding decision making. In this review, we briefly discuss the RT-related applications of MRI in the management of brain, prostate, and GI malignancies.

Development of an anthropomorphic multimodality pelvic phantom for quantitative evaluation of a deep-learning-based synthetic computed tomography generation technique

PURPOSE

The objective of this study was to fabricate an anthropomorphic multimodality pelvic phantom to evaluate a deep-learning-based synthetic computed tomography (CT) algorithm for magnetic resonance (MR)-only radiotherapy.

METHODS

Polyurethane-based and silicone-based materials with various silicone oil concentrations were scanned using 0.35 T MR and CT scanner to determine the tissue surrogate. Five tissue surrogates were determined by comparing the organ intensity with patient CT and MR images. Patient-specific organ modeling for three-dimensional printing was performed by manually delineating the structures of interest. The phantom was finally fabricated by casting materials for each structure. For the quantitative evaluation, the mean and standard deviations were measured within the regions of interest on the MR, simulation CT (CT_sim ), and synthetic CT (CT_syn ) images. Intensity-modulated radiation therapy plans were generated to assess the impact of different electron density assignments on plan quality using CT_sim and CT_syn . The dose calculation accuracy was investigated in terms of gamma analysis and dose-volume histogram parameters.

RESULTS

For the prostate site, the mean MR intensities for the patient and phantom were 78.1 ± 13.8 and 86.5 ± 19.3, respectively. The mean intensity of the synthetic image was 30.9 Hounsfield unit (HU), which was comparable to that of the real CT phantom image. The original and synthetic CT intensities of the fat tissue in the phantom were -105.8 ± 4.9 HU and -107.8 ± 7.8 HU, respectively. For the target volume, the difference in D_95% was 0.32 Gy using CT_syn with respect to CT_sim values. The V_65Gy values for the bladder in the plans using CT_sim and CT_syn were 0.31% and 0.15%, respectively.

CONCLUSION

This work demonstrated that the anthropomorphic phantom was physiologically and geometrically similar to the patient organs and was employed to quantitatively evaluate the deep-learning-based synthetic CT algorithm.

Evaluation of Multisource Adaptive MRI Fusion for Gross Tumor Volume Delineation of Hepatocellular Carcinoma

Purpose

Tumor delineation plays a critical role in radiotherapy for hepatocellular carcinoma (HCC) patients. The incorporation of MRI might improve the ability to correctly identify tumor boundaries and delineation consistency. In this study, we evaluated a novel Multisource Adaptive MRI Fusion (MAMF) method in HCC patients for tumor delineation.

Methods

Ten patients with HCC were included in this study retrospectively. Contrast-enhanced T1-weighted MRI at portal-venous phase (T1W_PP), contrast-enhanced T1-weighted MRI at 19-min delayed phase (T1W_DP), T2-weighted (T2W), and diffusion-weighted MRI (DWI) were acquired on a 3T MRI scanner and imported to in-house-developed MAMF software to generate synthetic MR fusion images. The original multi-contrast MR image sets were registered to planning CT by deformable image registration (DIR) using MIM. Four observers independently delineated gross tumor volumes (GTVs) on the planning CT, four original MR image sets, and the fused MRI for all patients. Tumor contrast-to-noise ratio (CNR) and Dice similarity coefficient (DSC) of the GTVs between each observer and a reference observer were measured on the six image sets. Inter-observer and inter-patient mean, SD, and coefficient of variation (CV) of the DSC were evaluated.

Results

Fused MRI showed the highest tumor CNR compared to planning CT and original MR sets in the ten patients. The mean ± SD tumor CNR was 0.72 ± 0.73, 3.66 ± 2.96, 4.13 ± 3.98, 4.10 ± 3.17, 5.25 ± 2.44, and 9.82 ± 4.19 for CT, T1W_PP, T2W, DWI, T1W_DP, and fused MRI, respectively. Fused MRI has the minimum inter-observer and inter-patient variations as compared to original MR sets and planning CT sets. GTV delineation inter-observer mean DSC across the ten patients was 0.81 ± 0.09, 0.85 ± 0.08, 0.88 ± 0.04, 0.89 ± 0.08, 0.90 ± 0.04, and 0.95 ± 0.02 for planning CT, T1W_PP, T2W, DWI, T1W_DP, and fused MRI, respectively. The patient mean inter-observer CV of DSC was 3.3%, 3.2%, 1.7%, 2.6%, 1.5%, and 0.9% for planning CT, T1W_PP, T2W, DWI, T1W_DP, and fused MRI, respectively.

Conclusion

The results demonstrated that the fused MRI generated using the MAMF method can enhance tumor CNR and improve inter-observer consistency of GTV delineation in HCC as compared to planning CT and four commonly used MR image sets (T1W_PP, T1W_DP, T2W, and DWI). The MAMF method holds great promise in MRI applications in HCC radiotherapy treatment planning.

Target localization accuracy in frame‐based stereotactic radiosurgery: Comparison between MR‐only and MR/CT co‐registration approaches

Purpose

In frame‐based Gamma Knife (GK) stereotactic radiosurgery two treatment planning workflows are commonly employed; one based solely on magnetic resonance (MR) images and the other based on magnetic resonance/computed tomography (MR/CT) co‐registered images. In both workflows, target localization accuracy (TLA) can be deteriorated due to MR‐related geometric distortions and/or MR/CT co‐registration uncertainties. In this study, the overall TLA following both clinical workflows is evaluated for cases of multiple brain metastases.

Methods

A polymer gel‐filled head phantom, having the Leksell stereotactic headframe attached, was CT‐imaged and irradiated by a GK Perfexion unit. A total of 26 4‐mm shots were delivered at 26 locations directly defined in the Leksell stereotactic space (LSS), inducing adequate contrast in corresponding T2‐weighted (T2w) MR images. Prescribed shot coordinates served as reference locations. An additional MR scan was acquired to implement the “mean image” distortion correction technique. The TLA for each workflow was assessed by comparing the radiation‐induced target locations, identified in MR images, with corresponding reference locations. Using T1w MR and CT images of 15 patients (totaling 81 lesions), TLA in clinical cases was similarly assessed, considering MR‐corrected data as reference. For the MR/CT workflow, both global and region of interest (ROI)‐based MR/CT registration approaches were studied.

Results

In phantom measurements, the MR‐corrected workflow demonstrated unsurpassed TLA (median offset of 0.2 mm) which deteriorated for MR‐only and MR/CT workflows (median offsets of 0.8 and 0.6 mm, respectively). In real‐patient cases, the MR‐only workflow resulted in offsets that exhibit a significant positive correlation with the distance from the MR isocenter, reaching 1.1 mm (median 0.6 mm). Comparable results were obtained for the MR/CT‐global workflow, although a maximum offset of 1.4 mm was detected. TLA was improved with the MR/CT‐ROI workflow resulting in median/maximum offsets of 0.4 mm/1.1 mm.

Conclusions

Subpixel TLA is achievable in all workflows. For the MR/CT workflow, a ROI‐based MR/CT co‐registration approach could considerably increase TLA and should be preferred instead of a global registration.

Multi-sequence MR generated sCT is promising for HNC MR-only RT: a comprehensive evaluation of previously developed sCT generation networks

PURPOSE

To verify the feasibility of our in-house developed multi-sequence magnetic resonance (MR)-generated synthetic computed tomography (sCT) for the accurate dose calculation and fractional positioning for head and neck MR-only radiation therapy (RT).

MATERIALS AND METHODS

Forty-five patients with nasopharyngeal carcinoma were retrospectively studied. By applying our previously in-house developed network, a patient's sCT can rapidly be generated with respect to feeding the sole T1 image, T1C image, T1DixonC image, T2 image, and their combination respectively (five pipelines in total). A k(5)-fold strategy was implemented during model establishment. Dose recalculation was performed for each pipeline generation to attain a dosimetric feasibility evaluation. Fractional positioning evaluation was performed by calculating the digitally reconstructed radiograph (DRR) of the sCT and planning CT and their offset to the portal image.

RESULTS

The dose mean absolute error values are (0.47±0.16)%, (0.48±0.15)% (p<0.05), (0.50±0.16)% (p<0.05), (0.50±0.15)% (p<0.05), and (0.45±0.16)% (p<0.05) for the T1, T1C, T1Dixon C, T2, and 4-channel generated sCT to the prescription dose, respectively. The 4-channel-generated sCT outperforms any other single-sequence pipelines. Among the single-sequence MR imaging-generated sCTs, the T1-generated shows the most accurate HU image quality and provide a reliable dose result. Quantified positioning errors with calculation of the difference to the planning CT offsets are (-0.26±0.50)mm, (-0.58±0.52)mm (p<0.05), (-0.27±0.57)mm (p>0.05), (-0.31±0.44)mm (p>0.05), and (-0.19±0.37)mm (p>0.05) at LNG and (0.34±0.53)mm, (0.48±0.56)mm (p>0.05), (0.55±0.56)mm (p>0.05), (0.37±0.61)mm (p>0.05), and (0.24±0.43)mm (p>0.05) at LAT of the anterior-posterior direction for the five pipelines.

CONCLUSION

Multi-sequence MR-generated sCT allows for accurate dose calculation and fractional positioning for head and neck MR-only RT. This article is protected by copyright. All rights reserved.

Head and neck cancer patient positioning using synthetic CT data in MRI-only radiation therapy

Purpose

The accuracy and precision of patient positioning is crucial in radiotherapy; however, there are no publications available using synthetic computed tomography (sCT) that evaluate rotations in head and neck (H&N) patients positioning or the effect of translation and rotation combined. The aim of this work was to evaluate the differences between using sCT with the CT for 2D‐ and 3D‐patient positioning in a magnetic resonance imaging (MRI)‐only workflow.

Methods

This study included 14 H&N cancer patients, with generated sCT data (MRI Planner v2.2) and the CT deformably registered to the MRI. Patient positioning was evaluated by comparing sCT against CT data: 3D cone beam CT (CBCT) was registered to the deformed CT (dCT) and sCT in six degrees of freedom (DoF) with a rigid auto‐registration algorithm and bone threshold, and 2D deformed digital reconstructed radiographs (dDRR) and synthetic DRRs (sDRR) were manually registered to orthogonal projections in five DoF by six blinded observers. The difference in displacement in all DoF were calculated for dCT and sCT, as well as for dDRR and sDRR. The interobserver variation was evaluated by separate application of the paired dDRR and sDRR registration matrices to the original coordinates of the planning target volume (PTV) structures and calculation of the Euclidean distance between the corresponding points. The Dice similarity coefficient (DSC) was calculated between dDRR/sDRR‐registered PTVs.

Results

The mean difference in patient positioning using CBCT was <0.7 mm and <0.3° and using orthogonal projections <0.4 mm and <0.2° in all directions. The maximum Euclidean distance was 5.1 mm, the corresponding mean (1SD) Euclidean distance and mean DSC were 3.5 ± 0.7 mm and 0.93, respectively.

Conclusions

This study shows that the sCT‐based patient positioning gives a comparable result with that based on CT images, allowing sCT to replace CT as reference for patient treatment positioning.

Prospective Clinical Feasibility Study for MRI-Only Brain Radiotherapy

OBJECTIVES

MRI-only radiotherapy (RT) provides a workflow to decrease the geometric uncertainty introduced by the image registration process between MRI and CT data and to streamline the RT planning. Despite the recent availability of validated synthetic CT (sCT) methods for the head region, there are no clinical implementations reported for brain tumors. Based on a preceding validation study of sCT, this study aims to investigate MRI-only brain RT through a prospective clinical feasibility study with endpoints for dosimetry and patient setup.

MATERIAL AND METHODS

Twenty-one glioma patients were included. MRI Dixon images were used to generate sCT images using a CE-marked deep learning-based software. RT treatment plans were generated based on MRI delineated anatomical structures and sCT for absorbed dose calculations. CT scans were acquired but strictly used for sCT quality assurance (QA). Prospective QA was performed prior to MRI-only treatment approval, comparing sCT and CT image characteristics and calculated dose distributions. Additional retrospective analysis of patient positioning and dose distribution gamma evaluation was performed.

RESULTS

Twenty out of 21 patients were treated using the MRI-only workflow. A single patient was excluded due to an MRI artifact caused by a hemostatic substance injected near the target during surgery preceding radiotherapy. All other patients fulfilled the acceptance criteria. Dose deviations in target were within ±1% for all patients in the prospective analysis. Retrospective analysis yielded gamma pass rates (2%, 2 mm) above 99%. Patient positioning using CBCT images was within ± 1 mm for registrations with sCT compared to CT.

CONCLUSION

We report a successful clinical study of MRI-only brain radiotherapy, conducted using both prospective and retrospective analysis. Synthetic CT images generated using the CE-marked deep learning-based software were clinically robust based on endpoints for dosimetry and patient positioning.

MRI classification using semantic random forest with auto-context model

BACKGROUND

It is challenging to differentiate air and bone on MR images of conventional sequences due to their low contrast. We propose to combine semantic feature extraction under auto-context manner into random forest to improve reasonability of the MRI segmentation for MRI-based radiotherapy treatment planning or PET attention correction.

METHODS

We applied a semantic classification random forest (SCRF) method which consists of a training stage and a segmentation stage. In the training stage, patch-based MRI features were extracted from registered MRI-CT training images, and the most informative elements were selected via feature selection to train an initial random forest. The rest sequence of random forests was trained by a combination of MRI feature and semantic feature under an auto-context manner. During segmentation, the MRI patches were first fed into these random forests to derive patch-based segmentation. By using patch fusion, the final end-to-end segmentation was obtained.

RESULTS

The Dice similarity coefficient (DSC) for air, bone and soft tissue classes obtained via proposed method were 0.976±0.007, 0.819±0.050 and 0.932±0.031, compared to 0.916±0.099, 0.673±0.151 and 0.830±0.083 with random forest (RF), and 0.942±0.086, 0.791±0.046 and 0.917±0.033 with U-Net. SCRF also outperformed the competing methods in sensitivity and specificity for all three structure types.

CONCLUSIONS

The proposed method accurately segmented bone, air and soft tissue. It is promising in facilitating advanced MR application in diagnosis and therapy.

The benefit of MR-only radiotherapy treatment planning for anal and rectal cancers: A planning study

INTRODUCTION

Limited evidence exists showing the benefit of magnetic resonance (MR)-only radiotherapy treatment planning for anal and rectal cancers. This study aims to assess the impact of MR-only planning on target volumes (TVs) and treatment plan doses to organs at risks (OARs) for anal and rectal cancers versus a computed tomography (CT)-only pathway.

MATERIALS AND METHODS

Forty-six patients (29 rectum and 17 anus) undergoing preoperative or radical external beam radiotherapy received CT and T2 MR simulation. TV and OARs were delineated on CT and MR, and volumetric arc therapy treatment plans were optimized independently (53.2 Gy/28 fractions for anus, 45 Gy/25 fractions for rectum). Further treatment plans assessed gross tumor volume (GTV) dose escalation. Differences in TV volumes and OAR doses, in terms of Vx Gy (organ volume (%) receiving x dose (Gy)), were assessed.

RESULTS

MR GTV and primary planning TV (PTV) volumes systematically reduced by 13 cc and 98 cc (anus) and 44 cc and 109 cc (rectum) respectively compared to CT volumes. Statistically significant OAR dose reductions versus CT were found for bladder and uterus (rectum) and bladder, penile bulb, and genitalia (anus). With GTV boosting, statistically significant dose reductions were found for sigmoid, small bowel, vagina, and penile bulb (rectum) and vagina (anus).

CONCLUSION

Our findings provide evidence that the introduction of MR (whether through MR-only or CT-MR pathways) to radiotherapy treatment planning for anal and rectal cancers has the potential to improve treatments. MR-related OAR dose reductions may translate into less treatment-related toxicity for patients or greater ability to dose escalate.

Dose accumulation of daily adaptive plans to decide optimal plan adaptation strategy for head-and-neck patients treated with MR-Linac

Advances in magnetic resonance linear accelerators (MR-Linacs) allow for superior visualization of soft tissue to guide online adaptive replanning for precise radiotherapy delivery. Elekta Unity MR-Linacs (Elekta AB, Stockholm, Sweden) provides 2 plan adaptation approaches, adapt-to-position (ATP), plan reoptimization based on the reference CT with the iso-shift measured from daily MR scans, and adapt-to-shape (ATS), full plan reoptimization based on the re-contoured daily MR scans. Our study aims to close the gap in knowledge regarding the use of the ATP technique in the treatment of head and neck (HN) cancers through the analysis of accumulated dose of daily ATP plans to organs at risk (OARs). Daily accumulated doses of 8 HN patients using deformable registration were analyzed to estimate the actual delivered dose versus the planned dose to evaluate the impact from daily anatomical changes and setup uncertainties. This process was completed through the collection of doses to OARs which were chosen based on the rigidity and size of the organ and the substantial dose it received. Results showed that the actual dose delivered to some OARs was significantly higher than the originally planned dose and was more pronounced in structures that were within the high-dose gradient for some subdisease sites. These findings suggest that the ATS approach should be used for plan adaptation in some specific HN diseases where OARs receive substantial dose with anatomy changes that could not be accounted for by the ATP approach. We also investigated the possibility of predicting the actual delivered dose at an early stage of the treatment course, with the intention of exploring a possibly more optimal alternative for planning through the combination of ATP and ATS approaches throughout treatment.

Technical overview of magnetic resonance fingerprinting and its applications in radiation therapy

Magnetic resonance fingerprinting (MRF) is an emerging imaging technique for rapid and simultaneous quantification of multiple tissue properties. The technique has been developed for quantitative imaging of different organs. The obtained quantitative measures have the potential to improve multiple steps of a typical radiotherapy workflow and potentially further improve integration of magnetic resonance imaging guided clinical decision making. In this review paper, we first provide a technical overview of the MRF method from data acquisition to postprocessing, along with recent development in advanced reconstruction methods. We further discuss critical aspects that could influence its usage in radiation therapy, such as accuracy and precision, repeatability and reproducibility, geometric distortion, and motion robustness. Finally, future directions for MRF application in radiation therapy are discussed.

Deep learning methods to generate synthetic CT from MRI in radiotherapy: A literature review

PURPOSE

In radiotherapy, MRI is used for target volume and organs-at-risk delineation for its superior soft-tissue contrast as compared to CT imaging. However, MRI does not provide the electron density of tissue necessary for dose calculation. Several methods of synthetic-CT (sCT) generation from MRI data have been developed for radiotherapy dose calculation. This work reviewed deep learning (DL) sCT generation methods and their associated image and dose evaluation, in the context of MRI-based dose calculation.

METHODS

We searched the PubMed and ScienceDirect electronic databases from January 2010 to March 2021. For each paper, several items were screened and compiled in figures and tables.

RESULTS

This review included 57 studies. The DL methods were either generator-only based (45% of the reviewed studies), or generative adversarial network (GAN) architecture and its variants (55% of the reviewed studies). The brain and pelvis were the most commonly investigated anatomical localizations (39% and 28% of the reviewed studies, respectively), and more rarely, the head-and-neck (H&N) (15%), abdomen (10%), liver (5%) or breast (3%). All the studies performed an image evaluation of sCTs with a diversity of metrics, with only 36 studies performing dosimetric evaluations of sCT.

CONCLUSIONS

The median mean absolute errors were around 76 HU for the brain and H&N sCTs and 40 HU for the pelvis sCTs. For the brain, the mean dose difference between the sCT and the reference CT was <2%. For the H&N and pelvis, the mean dose difference was below 1% in most of the studies. Recent GAN architectures have advantages compared to generator-only, but no superiority was found in term of image or dose sCT uncertainties. Key challenges of DL-based sCT generation methods from MRI in radiotherapy is the management of movement for abdominal and thoracic localizations, the standardization of sCT evaluation, and the investigation of multicenter impacts.

Deep learning based synthetic-CT generation in radiotherapy and PET: A review

Recently,deep learning (DL)-based methods for the generation of synthetic computed tomography (sCT) have received significant research attention as an alternative to classical ones. We present here a systematic review of these methods by grouping them into three categories, according to their clinical applications: (i) to replace computed tomography in magnetic resonance (MR) based treatment planning, (ii) facilitate cone-beam computed tomography based image-guided adaptive radiotherapy, and (iii) derive attenuation maps for the correction of positron emission tomography. Appropriate database searching was performed on journal articles published between January 2014 and December 2020. The DL methods' key characteristics were extracted from each eligible study, and a comprehensive comparison among network architectures and metrics was reported. A detailed review of each category was given, highlighting essential contributions, identifying specific challenges, and summarizing the achievements. Lastly, the statistics of all the cited works from various aspects were analyzed, revealing the popularity and future trends and the potential of DL-based sCT generation. The current status of DL-based sCT generation was evaluated, assessing the clinical readiness of the presented methods.

Harnessing benefit of highly conformal RT techniques for lymphoma patients

This status article describes current state-of-the-art radiotherapy for lymphomas and new emerging techniques. Current state-of-the-art radiotherapy is sophisticated, individualised, CT-based, intensity-modulated treatment, using PET/CT to define the target. The concept of involved site radiotherapy should be used, delineating the target using the exact same principles as for solid tumours. The optimal treatment delivery includes motion management and online treatment verification systems, which reduce intra- and interfractional anatomical variation. Emerging radiotherapy techniques in lymphomas include adaptive radiotherapy in MR- and CT-based treatment systems and proton therapy. The next generation linear accelerators have the capability to deliver adaptive treatment and allow relatively quick online adaptation to the daily variations of the anatomy. The computer systems use machine leaning to facilitate rapid automatic contouring of the target and organs-at-risk. Moreover, emerging MR-based planning and treatment facilities allow target definition directly from MR scans and allow intra-fractional tracking of structures recognisable on MR. Proton facilities are now being widely implemented. The benefits of proton therapy are due to the physical properties of protons, which in many cases allow sparing of normal tissue. The variety of techniques in modern radiotherapy means that the radiation oncologist must be able to choose the right technique for each patient. The choice is mainly based on experience and standard protocols, but new systems calculating risks for the patients with a specific treatment plan and also systems integrating clinical factors and risk factors into the planning process itself are emerging.

Toward PET/MRI as one-stop shop for radiotherapy planning in cervical cancer patients

BACKGROUND

Radiotherapy (RT) planning for cervical cancer patients entails the acquisition of both Computed Tomography (CT) and Magnetic Resonance Imaging (MRI). Further, molecular imaging by Positron Emission Tomography (PET) could contribute to target volume delineation as well as treatment response monitoring. The objective of this study was to investigate the feasibility of a PET/MRI-only RT planning workflow of patients with cervical cancer. This includes attenuation correction (AC) of MRI hardware and dedicated positioning equipment as well as evaluating MRI-derived synthetic CT (sCT) of the pelvic region for positioning verification and dose calculation to enable a PET/MRI-only setup.

MATERIAL AND METHODS

16 patients underwent PET/MRI using a dedicated RT setup after the routine CT (or PET/CT), including eight pilot patients and eight cervical cancer patients who were subsequently referred for RT. Data from 18 patients with gynecological cancer were added for training a deep convolutional neural network to generate sCT from Dixon MRI. The mean absolute difference between the dose distributions calculated on sCT and a reference CT was measured in the RT target volume and organs at risk. PET AC by sCT and a reference CT were compared in the tumor volume.

RESULTS

All patients completed the examination. sCT was inferred for each patient in less than 5 s. The dosimetric analysis of the sCT-based dose planning showed a mean absolute error (MAE) of 0.17 ± 0.12 Gy inside the planning target volumes (PTV). PET images reconstructed with sCT and CT had no significant difference in quantification for all patients.

CONCLUSIONS

These results suggest that multiparametric PET/MRI can be successfully integrated as a one-stop-shop in the RT workflow of patients with cervical cancer.

Task group 284 report: magnetic resonance imaging simulation in radiotherapy: considerations for clinical implementation, optimization, and quality assurance

The use of dedicated magnetic resonance simulation (MR-SIM) platforms in Radiation Oncology has expanded rapidly, introducing new equipment and functionality with the overall goal of improving the accuracy of radiation treatment planning. However, this emerging technology presents a new set of challenges that need to be addressed for safe and effective MR-SIM implementation. The major objectives of this report are to provide recommendations for commercially available MR simulators, including initial equipment selection, siting, acceptance testing, quality assurance, optimization of dedicated radiation therapy specific MR-SIM workflows, patient-specific considerations, safety, and staffing. Major contributions include guidance on motion and distortion management as well as MRI coil configurations to accommodate patients immobilized in the treatment position. Examples of optimized protocols and checklists for QA programs are provided. While the recommendations provided here are minimum requirements, emerging areas and unmet needs are also highlighted for future development.

Dose Calculation Algorithms for External Radiation Therapy: An Overview for Practitioners

Radiation therapy (RT) is a constantly evolving therapeutic technique; improvements are continuously being introduced for both methodological and practical aspects. Among the features that have undergone a huge evolution in recent decades, dose calculation algorithms are still rapidly changing. This process is propelled by the awareness that the agreement between the delivered and calculated doses is of paramount relevance in RT, since it could largely affect clinical outcomes. The aim of this work is to provide an overall picture of the main dose calculation algorithms currently used in RT, summarizing their underlying physical models and mathematical bases, and highlighting their strengths and weaknesses, referring to the most recent studies on algorithm comparisons. This handy guide is meant to provide a clear and concise overview of the topic, which will prove useful in helping clinical medical physicists to perform their responsibilities more effectively and efficiently, increasing patient benefits and improving the overall quality of the management of radiation treatment.

CT synthesis from MRI using multi-cycle GAN for head-and-neck radiation therapy

Magnetic Resonance Imaging (MRI) guided Radiation Therapy is a hot topic in the current studies of radiotherapy planning, which requires using MRI to generate synthetic Computed Tomography (sCT). Despite recent progress in image-to-image translation, it remains challenging to apply such techniques to generate high-quality medical images. This paper proposes a novel framework named Multi-Cycle GAN, which uses the Pseudo-Cycle Consistent module to control the consistency of generation and the domain control module to provide additional identical constraints. Besides, we design a new generator named Z-Net to improve the accuracy of anatomy details. Extensive experiments show that Multi-Cycle GAN outperforms state-of-the-art CT synthesis methods such as Cycle GAN, which improves MAE to 0.0416, ME to 0.0340, PSNR to 39.1053.

Head-to-Head Comparison of Neck 18F-FDG PET/MR and PET/CT in the Diagnosis of Differentiated Thyroid Carcinoma Patients after Comprehensive Treatment

Simple Summary

The most advanced positron emission tomography–magnetic resonance (PET/MR) combines the high soft tissue contrast of MRI with the high functional/metabolic sensitivity of PET and has the potential to achieve the highest level of diagnostic performance for refractory malignancies in differentiated thyroid cancer (DTC) patients. The utility of PET/MR in the postoperative follow-up of DTC patients has been relatively ambiguous. This retrospective study compared ¹⁸F-fluorodeoxyglucose neck PET/MR with PET/CT head-to-head, in order to evaluate the diagnostic efficacy of PET/MR in assessment malignancy in DTC patients after comprehensive treatment. We determined that PET/MR presented better detection rates, image conspicuity, and sensitivity than PET/CT in recurrent DTC lesions and cervical lymph node metastases. The addition of neck PET/MR scan after whole-body PET/CT may provide more favorable diagnostic information.

Abstract

We explored the clinical value of ¹⁸F-FDG PET/MR in a head-to-head comparison with PET/CT in loco-regional recurrent and metastatic cervical lymph nodes of differentiated thyroid carcinoma (DTC) patients after comprehensive treatment. ¹⁸F-FDG PET/CT and neck PET/MR scans that were performed in DTC patients with suspected recurrence or cervical lymph node metastasis after comprehensive treatment were retrospectively analyzed. Detection rates, diagnostic efficacy, image conspicuity, and measured parameters were compared between ¹⁸F-FDG PET/CT and PET/MR. The gold standard was histopathological diagnosis or clinical and imaging follow-up results for more than 6 months. Among the 37 patients enrolled, no suspicious signs of tumor were found in 10 patients, 24 patients had lymph node metastasis, and 3 patients had both recurrence and lymph node metastases. A total of 130 lesions were analyzed, including 3 malignant and 6 benign thyroid nodules, as well as 74 malignant and 47 benign cervical lymph nodes. Compared with PET/CT, PET/MR presented better detection rates (91.5% vs. 80.8%), image conspicuity (2.74 ± 0.60 vs. 1.9 ± 0.50, p < 0.001, especially in complex level II), and sensitivity (80.5% vs. 61.0%). SUVmax differed in benign and malignant lymph nodes in both imaging modalities (p < 0.05). For the same lesion, the SUVmax, SUVmean, and diameters measured by PET/MR and PET/CT were consistent and had significant correlation. In conclusion, compared with ¹⁸F-FDG PET/CT, PET/MR was more accurate in determining recurrent and metastatic lesions, both from a patient-based and from a lesion-based perspective. Adding local PET/MR after whole-body PET/CT may be recommended to provide more precise diagnostic information and scope of surgical resection without additional ionizing radiation. Further scaling-up prospective studies and economic benefit analysis are expected.

Effects of MR imaging time reduction on substitute CT generation for prostate MRI-only treatment planning

The introduction of MRI linear accelerators (MR-linacs) and the increased use of MR imaging in radiotherapy, requires improved approaches to MRI-only radiotherapy. MRI provides excellent soft tissue visualisation but does not provide any electron density information required for radiotherapy dose calculation, instead MRI is registered to CT images to enable dose calculations. MRI-only radiotherapy eliminates registration errors and reduces patient discomfort, workload and cost. Electron density requirements may be addressed in different ways, from manually applying bulk density corrections, to more computationally intensive methods to produce substitute CT datasets (sCT), requiring additional sequences, increasing overall imaging time. Reducing MR imaging time would reduce potential artefacts from intrafraction motion and patient discomfort. The aim of this study was to assess the impact of reducing MR imaging time on a hybrid atlas-voxel sCT conversion for prostate MRI-only treatment planning, considering both anatomical and dosimetric parameters. 10 volunteers were scanned on a Siemens Skyra 3T MRI. Sequences included the 3D T2-weighted (T2-w) SPACE sequence used for sCT conversion as previously validated against CT, along with variations to this sequence in repetition time (TR), turbo factor, and combinations of these to reduce the imaging time. All scans were converted to sCT and were compared to the sCT from the original SPACE sequence, evaluating for anatomical changes and dosimetric differences for a standard prostate VMAT plan. Compared to the previously validated T2-w SPACE sequence, scan times were reduced by up to 80%. The external volume and bony anatomy were compared, with all but one sequence meeting a DICE coefficient of 0.9 or better, with the largest variations occurring at the edges of the external body volume. The generated sCT agreed with the gold standard sCT within an isocentre dose of 1% and a gamma pass rate of 99% for a 1%/1 mm gamma tolerance for all but one sequence. This study demonstrates that the MR imaging sequence time was able to be reduced by approximately 80% with similar dosimetric results.

Patient position verification in magnetic-resonance imaging only radiotherapy of anal and rectal cancers

BACKGROUND AND PURPOSE

Magnetic resonance (MR)-only treatment pathways require either the MR-simulation or synthetic-computed tomography (sCT) as an alternative reference image for cone beam computed tomography (CBCT) patient position verification. This study assessed whether using T2 MR or sCT as CBCT reference images introduces systematic registration errors as compared to CT for anal and rectal cancers.

MATERIALS AND METHODS

A total of 32 patients (18 rectum,14 anus) received pre-treatment CT- and T2 MR- simulation. Routine treatment CBCTs were acquired. sCTs were generated using a validated research model. The local clinical registration protocol, using a grey-scale registration algorithm, was performed for 216 CBCTs using CT, MR and sCT as the reference image. Linear mixed effects modelling identified systematic differences between modalities.

RESULTS

Systematic translation and rotation differences to CT for MR were -0.3 to + 0.3 mm and -0.1 to 0.4° for anal cancers and -0.4 to 0.0 mm and 0.0 to 0.1° for rectal cancers, and for sCT were -0.4 to + 0.8 mm, -0.1 to 0.2° for anal cancers and -0.6 to + 0.2 mm, -0.1 to + 0.1° for rectal cancers.

CONCLUSIONS

T2 MR or sCT can successfully be used as reference images for anal and rectal cancer CBCT position verification with systematic differences to CT <±1 mm and <±0.5°. Clinical enabling of alternative modalities as reference images by vendors is required to reduce challenges associated with their use.

Feasibility of using a single MRI acquisition for fiducial marker localization and synthetic CT generation towards MRI-only prostate radiation therapy treatment planning

Purpose.To investigate the feasibility of using a single MRI acquisition for fiducial marker identification and synthetic CT (sCT) generation towards MRI-only treatment planning for prostate external beam radiation therapy (EBRT).Methods.Seven prostate cancer patients undergoing EBRT, each with three implanted gold fiducial markers, participated in this study. In addition to the planning CT scan, all patients were scanned on a 3 T MR scanner with a 3D double-echo gradient echo (GRE) sequence. Quantitative susceptibility mapping (QSM) was performed for marker localization. QSM-derived marker positions were compared to those from CT. The bulk density assignment technique for sCT generation was adopted. The magnitude GRE images were segmented into muscle, bone, fat, and air using a combination of unsupervised intensity-based classification of soft tissue and convolutional neural networks (CNN) for bone segmentation.Results.All implanted markers were visualized and accurately identified (average error: 0.7 ± 0.5 mm). QSM generated distinctive contrast for hemorrhage, calcifications, and gold fiducial markers. The estimated susceptibility/HU values on QSM/CT for gold and calcifications were 31.5 ± 2.9 ppm/1220 ± 100 HU and 14.6 ± 0.9 ppm/440 ± 100 HU, respectively. The intensity-based soft tissue classification resulted in an average Dice score of 0.97 ± 0.02; bone segmentation using CNN resulted in an average Dice score of 0.93 ± 0.03.Conclusion.This work indicates the feasibility of simultaneous fiducial marker identification and sCT generation using a single MRI acquisition. Future works includes evaluation of the proposed method in a large cohort of patients with optimized acquisition parameters as well as dosimetric evaluations.

Three-Dimensional Printed Anatomic Models Derived From Magnetic Resonance Imaging Data: Current State and Image Acquisition Recommendations for Appropriate Clinical Scenarios

Three-dimensional (3D) printing technologies have been increasingly utilized in medicine over the past several years and can greatly facilitate surgical planning thereby improving patient outcomes. Although still much less utilized compared to computed tomography (CT), magnetic resonance imaging (MRI) is gaining traction in medical 3D printing. The purpose of this study was two-fold: 1) to determine the prevalence in the existing literature of using MRI to create 3D printed anatomic models for surgical planning and 2) to provide image acquisition recommendations for appropriate clinical scenarios where MRI is the most suitable imaging modality. The workflow for creating 3D printed anatomic models from medical imaging data is complex and involves image segmentation of the regions of interest and conversion of that data into 3D surface meshes, which are compatible with printing technologies. CT is most commonly used to create 3D printed anatomic models due to the high image quality and relative ease of performing image segmentation from CT data. As compared to CT datasets, 3D printing using MRI data offers advantages since it provides exquisite soft tissue contrast needed for accurate organ segmentation and it does not expose patients to unnecessary ionizing radiation. MRI, however, often requires complicated imaging techniques and time-consuming postprocessing procedures to generate high-resolution 3D anatomic models needed for 3D printing. Despite these challenges, 3D modeling and printing from MRI data holds great clinical promises thanks to emerging innovations in both advanced MRI imaging and postprocessing techniques. EVIDENCE LEVEL: 2 TECHNICAL EFFICATCY: 5.

Evaluation of a commercial synthetic computed tomography generation solution for magnetic resonance imaging-only radiotherapy

Purpose

To evaluate the Siemens solution generating Synthetic computed tomography (sCT) for magnetic resonance imaging (MRI)‐only radiotherapy (RT).

Method

A retrospective study was conducted on 47 patients treated with external beam RT for brain or prostate cancer who underwent both MRI and CT for treatment planning. sCT images were generated from MRI using automatic bulk densities segmentation. The geometric accuracy of the sCT was assessed by comparing the Hounsfield Units (HU) difference between sCT and CT for bone structures, soft‐tissue, and full body contour. VMAT plans were computed on the CT for treatment preparation and then copied and recalculated with the same monitor units on the sCT using the AcurosXB algorithm. A 1%‐1mm gamma analysis was performed and DVH metrics for the Planning Target Volume (PTV) like the D_mean and the D_98% were compared. In addition, we evaluate the usability of sCT for daily position verification with cone beam computed tomography (CBCT) for 14 prostate patients by comparing sCT/CBCT registration results to CT/CBCT.

Results

Mean HU differences were small except for the skull (207 HU) and right femoral head of four patients where significant aberrations were found. The mean gamma pass rate was 73.2% for the brain and 84.7% for the prostate and D_mean were smaller than 0.5%. Large differences for the D_98% of the prostate group could be correlated to low Dice index of the PTV. The mean difference of translations and rotations were inferior to 3.5 mm and 0.2° in all directions with a major difference in the anterior‐posterior direction.

Conclusion

The performances of the software were shown to be similar to other sCT generation algorithms in terms of HU difference, dose comparison and daily image localization.

Magnetic resonance biomarkers in radiation oncology: The report of AAPM Task Group 294

A magnetic resonance (MR) biologic marker (biomarker) is a measurable quantitative characteristic that is an indicator of normal biological and pathogenetic processes or a response to therapeutic intervention derived from the MR imaging process. There is significant potential for MR biomarkers to facilitate personalized approaches to cancer care through more precise disease targeting by quantifying normal versus pathologic tissue function as well as toxicity to both radiation and chemotherapy. Both of which have the potential to increase the therapeutic ratio and provide earlier, more accurate monitoring of treatment response. The ongoing integration of MR into routine clinical radiation therapy (RT) planning and the development of MR guided radiation therapy systems is providing new opportunities for MR biomarkers to personalize and improve clinical outcomes. Their appropriate use, however, must be based on knowledge of the physical origin of the biomarker signal, the relationship to the underlying biological processes, and their strengths and limitations. The purpose of this report is to provide an educational resource describing MR biomarkers, the techniques used to quantify them, their strengths and weakness within the context of their application to radiation oncology so as to ensure their appropriate use and application within this field.

Artificial Intelligence in magnetic Resonance guided Radiotherapy: Medical and physical considerations on state of art and future perspectives

Over the last years, technological innovation in Radiotherapy (RT) led to the introduction of Magnetic Resonance-guided RT (MRgRT) systems. Due to the higher soft tissue contrast compared to on-board CT-based systems, MRgRT is expected to significantly improve the treatment in many situations. MRgRT systems may extend the management of inter- and intra-fraction anatomical changes, offering the possibility of online adaptation of the dose distribution according to daily patient anatomy and to directly monitor tumor motion during treatment delivery by means of a continuous cine MR acquisition. Online adaptive treatments require a multidisciplinary and well-trained team, able to perform a series of operations in a safe, precise and fast manner while the patient is waiting on the treatment couch. Artificial Intelligence (AI) is expected to rapidly contribute to MRgRT, primarily by safely and efficiently automatising the various manual operations characterizing online adaptive treatments. Furthermore, AI is finding relevant applications in MRgRT in the fields of image segmentation, synthetic CT reconstruction, automatic (on-line) planning and the development of predictive models based on daily MRI. This review provides a comprehensive overview of the current AI integration in MRgRT from a medical physicist's perspective. Medical physicists are expected to be major actors in solving new tasks and in taking new responsibilities: their traditional role of guardians of the new technology implementation will change with increasing emphasis on the managing of AI tools, processes and advanced systems for imaging and data analysis, gradually replacing many repetitive manual tasks.

An Inter-observer Study to Determine Radiotherapy Planning Target Volumes for Recurrent Gynaecological Cancer Comparing Magnetic Resonance Imaging Only With Computed Tomography-Magnetic Resonance Imaging

AIMS

Target delineation uncertainty is arguably the largest source of geometric uncertainty in radiotherapy. Several factors can affect it, including the imaging modality used for delineation. It is accounted for by applying safety margins to the target to produce a planning target volume (PTV), to which treatments are designed. To determine the margin, the delineation uncertainty is measured as the delineation error, and then a margin recipe used. However, there is no published evidence of such analysis for recurrent gynaecological cancers (RGC). The aims of this study were first to quantify the delineation uncertainty for RGC gross tumour volumes (GTVs) and to calculate the associated PTV margins and then to quantify the difference in GTV, delineation uncertainty and PTV margin, between a computed tomography-magnetic resonance imaging (CT-MRI) and MRI workflow.

MATERIALS AND METHODS

Seven clinicians delineated the GTV for 20 RGC tumours on co-registered CT and MRI datasets (CT-MRI) and on MRI alone. The delineation error, the standard deviation of distances from each clinician's outline to a reference, was measured and the required PTV margin determined. Differences between using CT-MRI and MRI alone were assessed.

RESULTS

The overall delineation error and the resulting margin were 3.1 mm and 8.5 mm, respectively, for CT-MRI, reducing to 2.5 mm and 7.1 mm, respectively, for MRI alone. Delineation errors and therefore the theoretical margins, varied widely between patients. MRI tumour volumes were on average 15% smaller than CT-MRI tumour volumes.

DISCUSSION

This study is the first to quantify delineation error for RGC tumours and to calculate the corresponding PTV margin. The determined margins were larger than those reported in the literature for similar patients, bringing into question both current margins and margin calculation methods. The wide variation in delineation error between these patients suggests that applying a single population-based margin may result in PTVs that are suboptimal for many. Finally, the reduced tumour volumes and safety margins suggest that patients with RGC may benefit from an MRI-only treatment workflow.

Synthetic CT for single-fraction neoadjuvant partial breast irradiation on an MRI-linac

A synthetic computed tomography (sCT) is required for daily plan optimization on an MRI-linac. Yet, only limited information is available on the accuracy of dose calculations on sCT for breast radiotherapy. This work aimed to (1) evaluate dosimetric accuracy of treatment plans for single-fraction neoadjuvant partial breast irradiation (PBI) on a 1.5 T MRI-linac calculated on a) bulk-density sCT mimicking the current MRI-linac workflow and b) deep learning-generated sCT, and (2) investigate the number of bulk-density levels required. For ten breast cancer patients we created three bulk-density sCTs of increasing complexity from the planning-CT, using bulk-density for: (1) body, lungs, and GTV (sCT_BD1); (2) volumes for sCT_BD1plus chest wall and ipsilateral breast (sCT_BD2); (3) volumes for sCT_BD2plus ribs (sCT_BD3); and a deep learning-generated sCT (sCT_DL) from a 1.5 T MRI in supine position. Single-fraction neoadjuvant PBI treatment plans for a 1.5 T MRI-linac were optimized on each sCT and recalculated on the planning-CT. Image evaluation was performed by assessing mean absolute error (MAE) and mean error (ME) in Hounsfield Units (HU) between the sCTs and the planning-CT. Dosimetric evaluation was performed by assessing dose differences, gamma pass rates, and dose-volume histogram (DVH) differences. The following results were obtained (median across patients for sCT_BD1/sCT_BD2/sCT_BD3/sCT_DLrespectively): MAE inside the body contour was 106/104/104/75 HU and ME was 8/9/6/28 HU, mean dose difference in the PTV_GTVwas 0.15/0.00/0.00/-0.07 Gy, median gamma pass rate (2%/2 mm, 10% dose threshold) was 98.9/98.9/98.7/99.4%, and differences in DVH parameters were well below 2% for all structures except for the skin in the sCT_DL. Accurate dose calculations for single-fraction neoadjuvant PBI on an MRI-linac could be performed on both bulk-density and deep learning sCT, facilitating further implementation of MRI-guided radiotherapy for breast cancer. Balancing simplicity and accuracy, sCT_BD2showed the optimal number of bulk-density levels for a bulk-density approach.

Clinical validation of a commercially available deep learning software for synthetic CT generation for brain

BACKGROUND

Most studies on synthetic computed tomography (sCT) generation for brain rely on in-house developed methods. They often focus on performance rather than clinical feasibility. Therefore, the aim of this work was to validate sCT images generated using a commercially available software, based on a convolutional neural network (CNN) algorithm, to enable MRI-only treatment planning for the brain in a clinical setting.

METHODS

This prospective study included 20 patients with brain malignancies of which 14 had areas of resected skull bone due to surgery. A Dixon magnetic resonance (MR) acquisition sequence for sCT generation was added to the clinical brain MR-protocol. The corresponding sCT images were provided by the software MRI Planner (Spectronic Medical AB, Sweden). sCT images were rigidly registered and resampled to CT for each patient. Treatment plans were optimized on CT and recalculated on sCT images for evaluation of dosimetric and geometric endpoints. Further analysis was also performed for the post-surgical cases. Clinical robustness in patient setup verification was assessed by rigidly registering cone beam CT (CBCT) to sCT and CT images, respectively.

RESULTS

All sCT images were successfully generated. Areas of bone resection due to surgery were accurately depicted. Mean absolute error of the sCT images within the body contour for all patients was 62.2 ± 4.1 HU. Average absorbed dose differences were below 0.2% for parameters evaluated for both targets and organs at risk. Mean pass rate of global gamma (1%/1 mm) for all patients was 100.0 ± 0.0% within PTV and 99.1 ± 0.6% for the full dose distribution. No clinically relevant deviations were found in the CBCT-sCT vs CBCT-CT image registrations. In addition, mean values of voxel-wise patient specific geometric distortion in the Dixon images for sCT generation were below 0.1 mm for soft tissue, and below 0.2 mm for air and bone.

CONCLUSIONS

This work successfully validated a commercially available CNN-based software for sCT generation. Results were comparable for sCT and CT images in both dosimetric and geometric evaluation, for both patients with and without anatomical anomalies. Thus, MRI Planner is feasible to use for radiotherapy treatment planning of brain tumours.

Initial Feasibility and Clinical Implementation of Daily MR-Guided Adaptive Head and Neck Cancer Radiation Therapy on a 1.5T MR-Linac System: Prospective R-IDEAL 2a/2b Systematic Clinical Evaluation of Technical Innovation

PURPOSE

This prospective study is, to our knowledge, the first report of daily adaptive radiation therapy (ART) for head and neck cancer (HNC) using a 1.5T magnetic resonance imaging-linear accelerator (MR-linac) with particular focus on safety and feasibility and dosimetric results of an online rigid registration-based adapt to position (ATP) workflow.

METHODS AND MATERIALS

Ten patients with HNC received daily ART on a 1.5T/7MV MR-linac, 6 using ATP only and 4 using ATP with 1 offline adapt-to-shape replan. Setup variability with custom immobilization masks was assessed by calculating the mean systematic error (M), standard deviation of the systematic error (Σ), and standard deviation of the random error (σ) of the isocenter shifts. Quality assurance was performed with a cylindrical diode array using 3%/3 mm γ criteria. Adaptive treatment plans were summed for each patient to compare the delivered dose with the planned dose from the reference plan. The impact of dosimetric variability between adaptive fractions on the summation plan doses was assessed by tracking the number of optimization constraint violations at each individual fraction.

RESULTS

The random errors (mm) for the x, y, and z isocenter shifts, respectively, were M = -0.3, 0.7, 0.1; Σ = 3.3, 2.6, 1.4; and σ = 1.7, 2.9, 1.0. The median (range) γ pass rate was 99.9% (90.9%-100%). The differences between the reference and summation plan doses were -0.61% to 1.78% for the clinical target volume and -11.74% to 8.11% for organs at risk (OARs), although an increase greater than 2% in OAR dose only occurred in 3 cases, each for a single OAR. All cases had at least 2 fractions with 1 or more constraint violations. However, in nearly all instances, constraints were still met in the summation plan despite multiple single-fraction violations.

CONCLUSIONS

Daily ART on a 1.5T MR-linac using an online ATP workflow is safe and clinically feasible for HNC and results in delivered doses consistent with planned doses.

Magnetic Resonance Imaging-Guided Adaptive Radiotherapy for Colorectal Liver Metastases

Technological advances have enabled well tolerated and effective radiation treatment for small liver metastases. Stereotactic ablative radiation therapy (SABR) refers to ablative dose delivery (>100 Gy BED) in five fractions or fewer. For larger tumors, the safe delivery of SABR can be challenging due to a more limited volume of healthy normal liver parenchyma and the proximity of the tumor to radiosensitive organs such as the stomach, duodenum, and large intestine. In addition to stereotactic treatment delivery, controlling respiratory motion, the use of image guidance, adaptive planning and increasing the number of radiation fractions are sometimes necessary for the safe delivery of SABR in these situations. Magnetic Resonance (MR) image-guided adaptive radiation therapy (MRgART) is a new and rapidly evolving treatment paradigm. MR imaging before, during and after treatment delivery facilitates direct visualization of both the tumor target and the adjacent normal healthy organs as well as potential intrafraction motion. Real time MR imaging facilitates non-invasive tumor tracking and treatment gating. While daily adaptive re-planning permits treatment plans to be adjusted based on the anatomy of the day. MRgART therapy is a promising radiation technology advance that can overcome many of the challenges of liver SABR and may facilitate the safe tumor dose escalation of colorectal liver metastases.

Synthetic dual-energy CT for MRI-only based proton therapy treatment planning using label-GAN

MRI-only treatment planning is highly desirable in the current proton radiation therapy workflow due to its appealing advantages such as bypassing MR-CT co-registration, avoiding x-ray CT exposure dose and reduced medical cost. However, MRI alone cannot provide stopping power ratio (SPR) information for dose calculations. Given that dual energy CT (DECT) can estimate SPR with higher accuracy than conventional single energy CT, we propose a deep learning-based method in this study to generate synthetic DECT (sDECT) from MRI to calculate SPR. Since the contrast difference between high-energy and low-energy CT (LECT) is important, and in order to accurately model this difference, we propose a novel label generative adversarial network-based model which can not only discriminate the realism of sDECT but also differentiate high-energy CT (HECT) and LECT from DECT. A cohort of 57 head-and-neck cancer patients with DECT and MRI pairs were used to validate the performance of the proposed framework. The results of sDECT and its derived SPR maps were compared with clinical DECT and the corresponding SPR, respectively. The mean absolute error for synthetic LECT and HECT were 79.98 ± 18.11 HU and 80.15 ± 16.27 HU, respectively. The corresponding SPR maps generated from sDECT showed a normalized mean absolute error as 5.22% ± 1.23%. By comparing with the traditional Cycle GANs, our proposed method significantly improves the accuracy of sDECT. The results indicate that on our dataset, the sDECT image form MRI is close to planning DECT, and thus shows promising potential for generating SPR maps for proton therapy.

Technical Challenges of Real-Time Adaptive MR-Guided Radiotherapy

In the past few years, radiotherapy (RT) has experienced a major technological innovation with the development of hybrid machines combining magnetic resonance (MR) imaging and linear accelerators. This new technology for MR-guided cancer treatment has the potential to revolutionize the field of adaptive RT due to the opportunity to provide high-resolution, real-time MR imaging before and during treatment application. However, from a technical point of view, several challenges remain which need to be tackled to ensure safe and robust real-time adaptive MR-guided RT delivery. In this manuscript, several technical challenges to MR-guided RT are discussed. Starting with magnetic field strength tradeoffs, the potential and limitations for purely MR-based RT workflows are discussed. Furthermore, the current status of real-time 3D MR imaging and its potential for real-time RT are summarized. Finally, the potential of quantitative MR imaging for future biological RT adaptation is highlighted.

Technical Note: Validation of an automatic ACR phantom quality assurance tool for an MR-guided radiotherapy system

PURPOSE

To modify and evaluate an automatic American College of Radiology (ACR) phantom analysis toolbox for ACR quality assurance (QA) on a low-field MR-guided radiotherapy system (ViewRay).

METHODS

An open-source toolbox was modified for ACR QA of a 0.35T MRI system (ViewRay MRIdian). A total of 17 ACR datasets were evaluated, including 10 datasets acquired from different systems across the world, and seven datasets acquired at our center between 2014 and 2020. All required ACR tests, geometric accuracy (GA), high-contrast spatial resolution (HCSR), slice thickness accuracy (ST), slice position accuracy (SP), percent integral uniformity (PIU), percentage signal ghosting (PSG), and low-contrast object detectability (LCOD), were assessed manually and using the toolbox automatically. Measurements between manual and automatic analysis were compared. Precision, recall, and accuracy were calculated, where the manual results were used as the ground truth.

RESULTS

The software took less than 2 min to complete all seven tests, which usually requires 40 min or more if analyzed manually. Overall, the automatic measurement was consistent with the manual result. The absolute differences between the two measurements were 0.72 ± 0.66 (mm), 0.01 ± 0.03, 0.50 ± 0.60 (mm), 0.39 ± 0.41 (mm), 1.01 ± 1.00 (%), 0.0016 ± 0.0019, and 2.79 ± 2.29 for the seven tests. The precision of the automatic toolbox was 100% for all tests, indicating that a test would 100% pass a manual analysis if it had passed the automatic analysis. Recall and accuracy were ≥ 96% for GA, HCSR, SP, PIU, PSG, and LCOD tests, and 91% for the ST test.

CONCLUSION

An automatic ACR QA tool was adopted and evaluated for the low-field MR-guided radiotherapy (MRgRT) system. Overall, the toolbox provided comparable results as manual analysis, and reduced the processing time from over 40 min to <2 min. This toolbox holds the potential to be widely adopted either as a second check tool or partially replace human measurement for MRgRT programs using the same system.

Multicentre, deep learning, synthetic-CT generation for ano-rectal MR-only radiotherapy treatment planning

BACKGROUND AND PURPOSE

Comprehensive dosimetric analysis is required prior to the clinical implementation of pelvic MR-only sites, other than prostate, due to the limited number of site specific synthetic-CT (sCT) dosimetric assessments in the literature. This study aims to provide a comprehensive assessment of a deep learning-based, conditional generative adversarial network (cGAN) model for a large ano-rectal cancer cohort. The following challenges were investigated; T2-SPACE MR sequences, patient data from multiple centres and the impact of sex and cancer site on sCT quality.

METHOD

RT treatment position CT and T2-SPACE MR scans, from two centres, were collected for 90 ano-rectal patients. A cGAN model trained using a focal loss function, was trained and tested on 46 and 44 CT-MR ano-rectal datasets, paired using deformable registration, respectively. VMAT plans were created on CT and recalculated on sCT. Dose differences and gamma indices assessed sCT dosimetric accuracy. A linear mixed effect (LME) model assessed the impact of centre, sex and cancer site.

RESULTS

A mean PTV D95% dose difference of 0.1% (range: -0.5% to 0.7%) was found between CT and sCT. All gamma index (1%/1 mm threshold) measurements were >99.0%. The LME model found the impact of modality, cancer site, sex and centre was clinically insignificant (effect ranges: -0.4% and 0.3%). The mean dose difference for all OAR constraints was 0.1%.

CONCLUSION

Focal loss cGAN models using T2-SPACE MR sequences from multiple centres can produce generalisable, dosimetrically accurate sCTs for ano-rectal cancers.

Accuracy of an MR-only workflow for prostate radiotherapy using semi-automatically burned-in fiducial markers

BACKGROUND

The benefit of MR-only workflow compared to current CT-based workflow for prostate radiotherapy is reduction of systematic errors in the radiotherapy chain by 2-3 mm. Nowadays, MRI is used for target delineation while CT is needed for position verification. In MR-only workflows, MRI based synthetic CT (sCT) replaces CT. Intraprostatic fiducial markers (FMs) are used as a surrogate for the position of the prostate improving targeting. However, FMs are not visible on sCT. Therefore, a semi-automatic method for burning-in FMs on sCT was developed. Accuracy of MR-only workflow using semi-automatically burned-in FMs was assessed and compared to CT/MR workflow.

METHODS

Thirty-one prostate cancer patients receiving radiotherapy, underwent an additional MR sequence (mDIXON) to create an sCT for MR-only workflow simulation. Three sources of accuracy in the CT/MR- and MR-only workflow were investigated. To compare image registrations for target delineation, the inter-observer error (IOE) of FM-based CT-to-MR image registrations and soft-tissue-based MR-to-MR image registrations were determined on twenty patients. Secondly, the inter-observer variation of the resulting FM positions was determined on twenty patients. Thirdly, on 26 patients CBCTs were retrospectively registered on sCT with burned-in FMs and compared to CT-CBCT registrations.

RESULTS

Image registration for target delineation shows a three times smaller IOE for MR-only workflow compared to CT/MR workflow. All observers agreed in correctly identifying all FMs for 18 out of 20 patients (90%). The IOE in CC direction of the center of mass (COM) position of the markers was within the CT slice thickness (2.5 mm), the IOE in AP and RL direction were below 1.0 mm and 1.5 mm, respectively. Registrations for IGRT position verification in MR-only workflow compared to CT/MR workflow were equivalent in RL-, CC- and AP-direction, except for a significant difference for random error in rotation.

CONCLUSIONS

MR-only workflow using sCT with burned-in FMs is an improvement compared to the current CT/MR workflow, with a three times smaller inter observer error in CT-MR registration and comparable CBCT registration results between CT and sCT reference scans. Trial registry Medical Research Involving Human Subjects Act (WMO) does apply to this study and was approved by the Medical Ethics review Committee of the Academic Medical Center. Registration number: NL65414.018.18. Date of registration: 21-08-2018.

Feasibility of MR-guided radiotherapy using beam-eye-view 2D-cine with tumor-volume projection

PURPOSE

Current magnetic resonance imaging (MRI) guided radiotherapy (MRgRT) applies sagittal/coronal 2D-cine to monitor major tumor motions, however, the beam eye's view (BEV) with volumetric tumor projection would be the best measure for radiation beam conformality, independent of tumor through-plane motion. The goal is to assess the feasibility, accuracy, and performance of the BEV approach.

METHODS

Beam-specific BEV 2D-cine with volume-projected tumor contours were simulated to establish a 2D/3D tumor match against a tumor-motion library based on multi-breath time-resolved (TR) 4DMRI images. Two BEV-library-matching methods were developed: (1) fast screening with tumor center-of-mass (∆COM), in-plane area ratio, and DICE similarity, and finalizing with the highest DICE score and (2) DICE screening for top-3 candidates and finalizing with rigid registration. A 4D-XCAT digital phantom and 8 lung-cancer patients were used for assessment. For each patient, 3 sets of 40 s TR-4DMRI were acquired at 2 Hz and 6 representative BEV were created with the isocenter set at tumor COM in mid-respiration. One TR-4DMRI set (40 × 2 = 80-images) was used to simulate BEV 2D-cine and the other two (160-images) were used to create a library. The matching result was validated against the ground truth within the test set. Using a leave-one-out strategy, the success rate, accuracy, and speed of tumor matching were assessed for volume-projected tumors over 11520 time-points (=8patients•3sets•80images•6BEVs).

RESULTS

Volume-projected tumor contour area on the 6 BEVs varies by 60% ± 8% and [Formula: see text] (in-plane/volume-projected) varies by 82% ± 9%. The [Formula: see text] changes with tumor shape, orientation, and through-plane motion. Method-1 produces 96% matching success (ΔCOM = 0.7 ± 0.2 mm, [Formula: see text]=1.01 ± 0.02, Dice=0.92 ± 0.02) with the computational time of 15 ± 1 ms/match, while method-2 produces 94% ± 1% success (ΔCOM = 0.2 ± 0.1 mm, [Formula: see text]=1.00 ± 0.01, Dice = 0.94 ± 0.02) with 223 ± 13 ms/match.

CONCLUSION

This study has demonstrated the feasibility, accuracy, and benefits of BEV 2D-cine imaging with tumor-volume projection, allowing real-time tumor motion monitoring and beam conformality checking. Further clinical evaluation is necessary before MRgRT applications.

An MR-only acquisition and artificial intelligence based image-processing protocol for photon and proton therapy using a low field MR

OBJECTIVE

Recent developments on synthetically generated CTs (sCT), hybrid MRI linacs and MR-only simulations underlined the clinical feasibility and acceptance of MR guided radiation therapy. However, considering clinical application of open and low field MR with a limited field of view can result in truncation of the patient's anatomy which further affects the MR to sCT conversion. In this study an acquisition protocol and subsequent MR image stitching is proposed to overcome the limited field of view restriction of open MR scanners, for MR-only photon and proton therapy.

MATERIAL AND METHODS

12 prostate cancer patients scanned with an open 0.35T scanner were included. To obtain the full body contour an enhanced imaging protocol including two repeated scans after bilateral table movement was introduced. All required structures (patient contour, target and organ at risk) were delineated on a post-processed combined transversal image set (stitched MRI). The postprocessed MR was converted into a sCT by a pretrained neural network generator. Inversely planned photon and proton plans (VMAT and SFUD) were designed using the sCT and recalculated for rigidly and deformably registered CT images and compared based on D2%, D50%, V70Gy for organs at risk and based on D2%, D50%, D98% for the CTV and PTV. The stitched MRI and the untruncated MRI were compared to the CT, and the maximum surface distance was calculated. The sCT was evaluated with respect to delineation accuracy by comparing on stitched MRI and sCT using the DICE coefficient for femoral bones and the whole body.

RESULTS

Maximum surface distance analysis revealed uncertainties in lateral direction of 1-3mm on average. DICE coefficient analysis confirms good performance of the sCT conversion, i.e. 92%, 93%, and 100% were obtained for femoral bone left and right and whole body. Dose comparison resulted in uncertainties below 1% between deformed CT and sCT and below 2% between rigidly registered CT and sCT in the CTV for photon and proton treatment plans.

DISCUSSION

A newly developed acquisition protocol for open MR scanners and subsequent Sct generation revealed good acceptance for photon and proton therapy. Moreover, this protocol tackles the restriction of the limited FOVs and expands the capacities towards MR guided proton therapy with horizontal beam lines.

Synthetic computed tomography data allows for accurate absorbed dose calculations in a magnetic resonance imaging only workflow for head and neck radiotherapy

•

The geometry of the synthetic CT is comparable to the CT in the H&N region.

•

Synthetic CT in the H&N region provides similar absorbed dose calculation as the CT.

•

Absorbed dose calculations in the dental region could benefit from using synthetic CT.

Background and purpose

Few studies on magnetic resonance imaging (MRI) only head and neck radiation treatment planning exist, and none using a generally available software. The aim of this study was to evaluate the accuracy of absorbed dose for head and neck synthetic computed tomography data (sCT) generated by a commercial convolutional neural network-based algorithm.

Materials and methods

For 44 head and neck cancer patients, sCT were generated and the geometry was validated against computed tomography data (CT). The clinical CT based treatment plan was transferred to the sCT and recalculated without re-optimization, and differences in relative absorbed dose were determined for dose-volume-histogram (DVH) parameters and the 3D volume.

Results

For overall body, the results of the geometric validation were (Mean ± 1sd): Mean error −5 ± 10 HU, mean absolute error 67 ± 14 HU, Dice similarity coefficient 0.98 ± 0.05, and Hausdorff distance difference 4.2 ± 1.7 mm. Water equivalent depth difference for region Th1-C7, mid mandible and mid nose were −0.3 ± 3.4, 1.1 ± 2.0 and 0.7 ± 3.8 mm respectively. The maximum mean deviation in absorbed dose for all DVH parameters was 0.30% (0.12 Gy). The absorbed doses were considered equivalent (p-value < 0.001) and the mean 3D gamma passing rate was 99.4 (range: 95.7–99.9%).

Conclusions

The convolutional neural network-based algorithm generates sCT which allows for accurate absorbed dose calculations for MRI-only head and neck radiation treatment planning. The sCT allows for statistically equivalent absorbed dose calculations compared to CT based radiotherapy.

Magnetic resonance imaging-based synthetic computed tomography of the lumbar spine for surgical planning: a clinical proof-of-concept

OBJECTIVE

Computed tomography scanning of the lumbar spine incurs a radiation dose ranging from 3.5 mSv to 19.5 mSv as well as relevant costs and is commonly necessary for spinal neuronavigation. Mitigation of the need for treatment-planning CT scans in the presence of MRI facilitated by MRI-based synthetic CT (sCT) would revolutionize navigated lumbar spine surgery. The authors aim to demonstrate, as a proof of concept, the capability of deep learning-based generation of sCT scans from MRI of the lumbar spine in 3 cases and to evaluate the potential of sCT for surgical planning.

METHODS

Synthetic CT reconstructions were made using a prototype version of the "BoneMRI" software. This deep learning-based image synthesis method relies on a convolutional neural network trained on paired MRI-CT data. A specific but generally available 4-minute 3D radiofrequency-spoiled T1-weighted multiple gradient echo MRI sequence was supplemented to a 1.5T lumbar spine MRI acquisition protocol.

RESULTS

In the 3 presented cases, the prototype sCT method allowed voxel-wise radiodensity estimation from MRI, resulting in qualitatively adequate CT images of the lumbar spine based on visual inspection. Normal as well as pathological structures were reliably visualized. In the first case, in which a spiral CT scan was available as a control, a volume CT dose index (CTDIvol) of 12.9 mGy could thus have been avoided. Pedicle screw trajectories and screw thickness were estimable based on sCT findings.

CONCLUSIONS

The evaluated prototype BoneMRI method enables generation of sCT scans from MRI images with only minor changes in the acquisition protocol, with a potential to reduce workflow complexity, radiation exposure, and costs. The quality of the generated CT scans was adequate based on visual inspection and could potentially be used for surgical planning, intraoperative neuronavigation, or for diagnostic purposes in an adjunctive manner.

MRI-guided Radiation Therapy: An Emerging Paradigm in Adaptive Radiation Oncology

Radiation therapy (RT) continues to be one of the mainstays of cancer treatment. Considerable efforts have been recently devoted to integrating MRI into clinical RT planning and monitoring. This integration, known as MRI-guided RT, has been motivated by the superior soft-tissue contrast, organ motion visualization, and ability to monitor tumor and tissue physiologic changes provided by MRI compared with CT. Offline MRI is already used for treatment planning at many institutions. Furthermore, MRI-guided linear accelerator systems, allowing use of MRI during treatment, enable improved adaptation to anatomic changes between RT fractions compared with CT guidance. Efforts are underway to develop real-time MRI-guided intrafraction adaptive RT of tumors affected by motion and MRI-derived biomarkers to monitor treatment response and potentially adapt treatment to physiologic changes. These developments in MRI guidance provide the basis for a paradigm change in treatment planning, monitoring, and adaptation. Key challenges to advancing MRI-guided RT include real-time volumetric anatomic imaging, addressing image distortion because of magnetic field inhomogeneities, reproducible quantitative imaging across different MRI systems, and biologic validation of quantitative imaging. This review describes emerging innovations in offline and online MRI-guided RT, exciting opportunities they offer for advancing research and clinical care, hurdles to be overcome, and the need for multidisciplinary collaboration.

Evaluation of patient-specific MR distortion correction schemes for improved target localization accuracy in SRS

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.

Integration of spatial distortion effects in a 4D computational phantom for simulation studies in extra-cranial MRI-guided radiation therapy: Initial results

PURPOSE

Spatial distortions in magnetic resonance imaging (MRI) are mainly caused by inhomogeneities of the static magnetic field, nonlinearities in the applied gradients, and tissue-specific magnetic susceptibility variations. These factors may significantly alter the geometrical accuracy of the reconstructed MR image, thus questioning the reliability of MRI for guidance in image-guided radiation therapy. In this work, we quantified MRI spatial distortions and created a quantitative model where different sources of distortions can be separated. The generated model was then integrated into a four-dimensional (4D) computational phantom for simulation studies in MRI-guided radiation therapy at extra-cranial sites.

METHODS

A geometrical spatial distortion phantom was designed in four modules embedding laser-cut PMMA grids, providing 3520 landmarks in a field of view of (345 × 260 × 480) mm³ . The construction accuracy of the phantom was verified experimentally. Two fast MRI sequences for extra-cranial imaging at 1.5 T were investigated, considering axial slices acquired with online distortion correction, in order to mimic practical use in MRI-guided radiotherapy. Distortions were separated into their sources by acquisition of images with gradient polarity reversal and dedicated susceptibility calculations. Such a separation yielded a quantitative spatial distortion model to be used for MR imaging simulations. Finally, the obtained spatial distortion model was embedded into an anthropomorphic 4D computational phantom, providing registered virtual CT/MR images where spatial distortions in MRI acquisition can be simulated.

RESULTS

The manufacturing accuracy of the geometrical distortion phantom was quantified to be within 0.2 mm in the grid planes and 0.5 mm in depth, including thickness variations and bending effects of individual grids. Residual spatial distortions after MRI distortion correction were strongly influenced by the applied correction mode, with larger effects in the trans-axial direction. In the axial plane, gradient nonlinearities caused the main distortions, with values up to 3 mm in a 1.5 T magnet, whereas static field and susceptibility effects were below 1 mm. The integration in the 4D anthropomorphic computational phantom highlighted that deformations can be severe in the region of the thoracic diaphragm, especially when using axial imaging with 2D distortion correction. Adaptation of the phantom based on patient-specific measurements was also verified, aiming at increased realism in the simulation.

CONCLUSIONS

The implemented framework provides an integrated approach for MRI spatial distortion modeling, where different sources of distortion can be quantified in time-dependent geometries. The computational phantom represents a valuable platform to study motion management strategies in extra-cranial MRI-guided radiotherapy, where the effects of spatial distortions can be modeled on synthetic images in a virtual environment.

Dosimetric impact of geometric distortions in an MRI-only proton therapy workflow for lung, liver and pancreas

In a radiation therapy workflow based on Magnetic Resonance Imaging (MRI), dosimetric errors may arise due to geometric distortions introduced by MRI. The aim of this study was to quantify the dosimetric effect of system-dependent geometric distortions in an MRI-only workflow for proton therapy applied at extra-cranial sites. An approach was developed, in which computed tomography (CT) images were distorted using an MRI displacement map, which represented the MR distortions in a spoiled gradient-echo sequence due to gradient nonlinearities and static magnetic field inhomogeneities. A retrospective study was conducted on 4DCT/MRI digital phantoms and 18 4DCT clinical datasets of the thoraco-abdominal site. The treatment plans were designed and separately optimized for each beam in a beam specific Planning Target Volume on the distorted CT, and the final dose distribution was obtained as the average. The dose was then recalculated in undistorted CT using the same beam geometry and beam weights. The analysis was performed in terms of Dose Volume Histogram (DVH) parameters. No clinically relevant dosimetric impact was observed on organs at risk, whereas in the target structure, geometric distortions caused statistically significant variations in the planned dose DVH parameters and dose homogeneity index (DHI). The dosimetric variations in the target structure were smaller in abdominal cases (ΔD_2%, ΔD_98%, and ΔD_mean all below 0.1% and ΔDHI below 0.003) compared to the lung cases. Indeed, lung patients with tumors isolated inside lung parenchyma exhibited higher dosimetric variations (ΔD_2%≥0.3%, ΔD_98%≥15.9%, ΔD_mean≥3.3% and ΔDHI≥0.102) than lung patients with tumor close to soft tissue (ΔD_2%≤0.4%, ΔD_98%≤5.6%, ΔD_mean≤0.9% and ΔDHI≤0.027) potentially due to higher density variations along the beam path. Results suggest the potential applicability of MRI-only proton therapy, provided that specific analysis is applied for isolated lung tumors.

Investigating the use of virtual 4DCT from 4DMRI in gated carbon ion radiation therapy of abdominal tumors

PURPOSE

To generate virtual 4DCT from 4DMRI with field of view (FOV) extended to the entire involved patient anatomy, in order to evaluate its use in carbon ion radiation therapy (CIRT) of the abdominal site in a clinical scenario.

MATERIALS AND METHODS

The virtual 4DCT was generated by deforming a reference CT in order to (1) match the anatomy depicted in the 4DMRI within its FOV, by calculating deformation fields with deformable image registration to describe inter-fractional and breathing motion, and (2) obtain physically plausible deformation outside of the 4DMRI FOV, by propagating and modulating the previously obtained deformation fields. The implemented method was validated on a digital anthropomorphic phantom, for which a ground truth (GT) 4DCT was available. A CIRT treatment plan was optimized at the end-exhale reference CT and the RBE-weighted dose distribution was recalculated on both the virtual and GT 4DCTs. The method estimation error was quantified by comparing the virtual and GT 4DCTs and the corresponding recomputed doses. The method was then evaluated on 8 patients with pancreas or liver tumors treated with CIRT using respiratory gating at end-exhale. The clinical treatment plans adopted at the National Center for Oncological Hadrontherapy (CNAO, Pavia, Italy) were considered and the dose distribution was recomputed on all respiratory phases of the planning and virtual 4DCTs. By comparing the two datasets and the corresponding dose distributions, the geometrical and dosimetric impact of organ motion was assessed.

RESULTS

For the phantom, the error outside of the 4DMRI FOV was up to 4.5mm, but it remained sub-millimetric in correspondence to the target within the 4DMRI FOV. Although the impact of motion on the target D_95% resulted in variations ranging from 22% to 90% between the planned dose and the doses recomputed on the GT 4DCT phases, the corresponding estimation error was ≤2.2%. In the patient cases, the variation of the baseline tumor position between the planning and the virtual end-exhale CTs presented a median (interquartile range) value of 6.0 (4.9) mm. For baseline variations larger than 5mm, the tumor D_95% variation between the plan and the dose recomputed on the end-exhale virtual CT resulted larger than 10%. Median variations higher than 10% in the target D_95% and gastro-intestinal OARs D_2% were quantified at the end-inhale, whereas close to the end-exhale phase, limited variations of relevant dose metrics were found for both tumor and OARs.

CONCLUSIONS

The negligible impact of the geometrical inaccuracy in the estimated anatomy outside of the 4DMRI FOV on the overall dosimetric accuracy suggests the feasibility of virtual 4DCT with extended FOV in CIRT of the abdominal site. In the analyzed patient group, inter-fractional variations such as baseline variation and breathing variability were quantified, demonstrating the method capability to support treatment planning in gated CIRT of the abdominal site.

Fully automated multiorgan segmentation in abdominal magnetic resonance imaging with deep neural networks

PURPOSE

Segmentation of multiple organs-at-risk (OARs) is essential for magnetic resonance (MR)-only radiation therapy treatment planning and MR-guided adaptive radiotherapy of abdominal cancers. Current practice requires manual delineation that is labor-intensive, time-consuming, and prone to intra- and interobserver variations. We developed a deep learning (DL) technique for fully automated segmentation of multiple OARs on clinical abdominal MR images with high accuracy, reliability, and efficiency.

METHODS

We developed Automated deep Learning-based abdominal multiorgan segmentation (ALAMO) technique based on two-dimensional U-net and a densely connected network structure with tailored design in data augmentation and training procedures such as deep connection, auxiliary supervision, and multiview. The model takes in multislice MR images and generates the output of segmentation results. 3.0-Tesla T1 VIBE (Volumetric Interpolated Breath-hold Examination) images of 102 subjects were used in our study and split into 66 for training, 16 for validation, and 20 for testing. Ten OARs were studied, including the liver, spleen, pancreas, left/right kidneys, stomach, duodenum, small intestine, spinal cord, and vertebral bodies. An experienced radiologist manually labeled each OAR, followed by reediting, if necessary, by a senior radiologist, to create the ground-truth. The performance was measured using volume overlapping and surface distance.

RESULTS

The ALAMO technique generated segmentation labels in good agreement with the manual results. Specifically, among the ten OARs, nine achieved high dice similarity coefficients (DSCs) in the range of 0.87-0.96, except for the duodenum with a DSC of 0.80. The inference completed within 1 min for a three-dimensional volume of 320 × 288 × 180. Overall, the ALAMO model matched the state-of-the-art techniques in performance.

CONCLUSION

The proposed ALAMO technique allows for fully automated abdominal MR segmentation with high accuracy and practical memory and computation time demands.