Spatial Precision in Magnetic Resonance Imaging-Guided Radiation Therapy: The Role of Geometric Distortion

Evaluating stereotactic accuracy with patient-specific MRI distortion corrections for frame-based radiosurgery

PURPOSE

This study examines how MRI distortions affect frame-based SRS treatments and assesses the need for clinical distortion corrections.

METHODS

The study included 18 patients with 80 total brain targets treated using frame-based radiosurgery. Distortion within patients' MRIs were corrected using Cranial Distortion Correction (CDC) software, which utilizes the patient's CT to alter planning MRIs to reduce inherent intra-cranial distortion. Distortion was evaluated by comparing the original planning target volumes (PTVORIG) to targets contoured on corrected MRIs (PTVCORR). To provide an internal control, targets were also re-contoured on uncorrected (PTVRECON) MRIs. Additional analysis was done to assess if 1 mm expansions to PTVORIG targets would compensate for patient-specific distortions. Changes in target volumes, DICE and JACCARD similarity coefficients, minimum PTV dose (Dmin), dose to 95% of the PTV (D95%), and normal tissue receiving 12 Gy (V12Gy), 10 Gy (V10Gy), and 5 Gy (V5Gy) were calculated and evaluated. Student's t-tests were used to determine if changes in PTVCORR were significantly different than intra-contouring variability quantified by PTVRECON.

RESULTS

PTVRECON and PTVCORR relative changes in volume were 6.19% ± 10.95% and 1.48% ± 32.92%. PTVRECON and PTVCORR similarity coefficients were 0.90 ± 0.08 and 0.73 ± 0.16 for DICE and 0.82 ± 0.12 and 0.60 ± 0.18 for JACCARD. PTVRECON and PTVCORR changes in Dmin were -0.88% ± 8.77% and -12.9 ± 17.3%. PTVRECON and PTVCORR changes in D95% were -0.34% ± 5.89 and -8.68% ± 13.21%. The 1 mm expanded PTVORIG targets did not entirely cover 14 of the 80 PTVCORR targets. Normal tissue changes (V12Gy, V10Gy, V5Gy) calculated with PTVRECON were (-0.09% ± 7.39%, -0.38% ± 5.67%, -0.08% ± 2.04%) and PTVCORR were (-2.14% ± 7.34%, -1.42% ± 5.45%, -0.61% ± 1.93%). Except for V10Gy, all PTVCORR changes were significantly different (p < 0.05) than PTVRECON.

CONCLUSION

MRIs used for SRS target delineation exhibit notable geometric distortions that may compromise optimal dosimetric accuracy. A uniform 1 mm expansion may result in geometric misses; however, the CDC algorithm provides a feasible solution for rectifying distortions, thereby enhancing treatment precision.

Geometrical accuracy of magnetic resonance imaging for ocular proton therapy planning

Background & purpose

Magnetic resonance imaging (MRI) is increasingly used in treatment preparation of ocular proton therapy, but its spatial accuracy might be limited by geometric distortions due to susceptibility artefacts. A correct geometry of the MR images is paramount since it defines where the dose will be delivered. In this study, we assessed the geometrical accuracy of ocular MRI.

Materials & methods

A dedicated ocular 3 T MRI protocol, with localized shimming and increased gradients, was compared to computed tomography (CT) and X-ray images in a phantom and in 15 uveal melanoma patients. The MRI protocol contained three-dimensional T2-weighted and T1-weighted sequences with an isotropic reconstruction resolution of 0.3-0.4 mm. Tantalum clips were identified by three observers and clip-clip distances were compared between T2-weighted and T1-weighted MRI, CT and X-ray images for the phantom and between MRI and X-ray images for the patients.

Results

Interobserver variability was below 0.35 mm for the phantom and 0.30(T1)/0.61(T2) mm in patients. Mean absolute differences between MRI and reference were below 0.27 ± 0.16 mm and 0.32 ± 0.23 mm for the phantom and in patients, respectively. In patients, clip-clip distances were slightly larger on MRI than on X-ray images (mean difference T1: 0.11 ± 0.38 mm, T2: 0.10 ± 0.44 mm). Differences did not increase at larger distances and did not correlate to interobserver variability.

Conclusions

A dedicated ocular MRI protocol can produce images of the eye with a geometrical accuracy below half the MRI acquisition voxel (<0.4 mm). Therefore, these images can be used for ocular proton therapy planning, both in the current model-based workflow and in proposed three-dimensional MR-based workflows.

Assessment of traumatic mandibular nerve using MR neurography sequence: a preliminary study

BACKGROUND

Iatrogenic mandibular nerve damage resulting from oral surgeries and dental procedures is painful and a formidable challenge for patients and oral surgeons alike, mainly because the absence of objective and quantitative methods for diagnosing nerve damage renders treatment and compensation ambiguous while often leading to medico-legal disputes. The aim of this study was to examine discriminating factors of traumatic mandibular nerve within a specific magnetic resonance imaging (MRI) protocol and to suggest tangible diagnostic criteria for peripheral trigeminal nerve injury.

METHODS

Twenty-six patients with ipsilateral mandibular nerve trauma underwent T2 Flex water, 3D short tau inversion recovery (STIR), and diffusion-weighted imaging (DWI) acquired by periodically rotating overlapping parallel lines with enhanced reconstruction (PROPELLER) pulse sequences; 26 injured nerves were thus compared with contra-lateral healthy nerves at anatomically corresponding sites. T2 Flex apparent signal to noise ratio (FSNR), T2 Flex apparent nerve-muscle contrast to noise ratio (FNMCNR) 3D STIR apparent signal to noise ratio (SSNR), 3D STIR apparent nerve-muscle contrast to noise ratio (SNMCNR), apparent diffusion coefficient (ADC) and area of cross-sectional nerve (Area) were evaluated.

RESULTS

Mixed model analysis revealed FSNR and FNMCNR to be the dual discriminators for traumatized mandibular nerve (p < 0.05). Diagnostic performance of both parameters was also determined with area under the receiver operating characteristic curve (AUC for FSNR = 0.712; 95% confidence interval [CI]: 0.5660, 0.8571 / AUC for FNMCNR = 0.7056; 95% confidence interval [CI]: 1.011, 1.112).

CONCLUSIONS

An increase in FSNR and FNMCNR within our MRI sequence seems to be accurate indicators of the presence of traumatic nerve. This prospective study may serve as a foundation for sophisticated model diagnosing trigeminal nerve trauma within large patient cohorts.

The role of artificial intelligence in cardiovascular magnetic resonance imaging

Cardiovascular magnetic resonance (CMR) imaging is the gold standard test for myocardial tissue characterization and chamber volumetric and functional evaluation. However, manual CMR analysis can be time-consuming and is subject to intra- and inter-observer variability. Artificial intelligence (AI) is a field that permits automated task performance through the identification of high-level and complex data relationships. In this review, we review the rapidly growing role of AI in CMR, including image acquisition, sequence prescription, artifact detection, reconstruction, segmentation, and data reporting and analysis including quantification of volumes, function, myocardial infarction (MI) and scar detection, and prediction of outcomes. We conclude with a discussion of the emerging challenges to widespread adoption and solutions that will allow for successful, broader uptake of this powerful technology.

Neuronavigation in endoscopic skull base surgery and the accuracy of different MRI sequences

OBJECTIVE

Neuronavigation is common technology used by skull base teams when performing endoscopic endonasal surgery. A common practice of MRI imagining is to obtain 3D isotopic gadolinium enhanced T1W magnetisation prepared rapid gradient echo (MPRAGE) sequences. These are prone to distortion when undertaken on 3 T magnets. The aim of this project is to compare the in vivo accuracy of MRI sequences between current and new high resolution 3D sequences. The goal is to determine if geometric distortion significantly affects neuronavigation accuracy.

METHODS

Patients were scanned with a 3D T1 MPRAGE sequence, 3D T1 SPACE sequence and a CT stereotactic localisation. Following general anaesthesia, patients were registered on the Stealth Station (Medtronic, USA) using a side mount emitter for Electromagnetic navigation. A variety of surgically relevant anatomical landmarks in the sagittal and coronal plane were selected with real and virtual data points measured.

RESULTS

A total of 10 patients agreed be enrolled in the study with datapoints collected during surgery. The distance between real and virtual datapoints trended to be lower in SPACE sequences compared to MPRAGE. Paired t test did not demonstrate a significant difference.

CONCLUSION

We have demonstrated that navigational accuracy is not significantly affected by the type of MRI sequence selected and that current corrective algorithms are sufficient. Navigational accuracy is affected by many factors, with registration error likely playing the most significant role. Further research involving real time imaging such as endoscopic ultrasound may hopefully address this potential error.

Geometric distortion caused by metallic femoral head prosthesis in prostate cancer imaging on an MR Linac: in-vivo measurements of spatial deformation

OBJECTIVES

Metallic implants cause artefacts and distortion on MRI. To ensure accurate dose delivery and plan adaptation on an MR Linac, there is a need to evaluate distortion caused.

METHODS

Participants were imaged on an MR Linac (Elekta Unity, Elekta AB Stockholm). Three sequences were evaluated. Two vendor supplied (T2W TSE 3D), and one T2W TSE 3D optimized to reduce metal artefact distortions. Images were rigidly registered to CT images by a single observer, using bony anatomy. Three coronal and three axial images were selected, and six paired, adjacent, bony landmarks were identified on each slice. Images bisecting treatment isocentre were included. Difference between landmark coordinates was taken to be measure of distortion.

RESULTS

Five observers participated. Thirty six pairs of bony landmarks were identified. Median difference in position of landmarks was ≤3 mm (range 0.3-4.4 mm). One-way analysis of variance (ANOVA) between observer means showed no significant variation between sequences or patients (P = 1.26 in plane, P = 0.11 through plane). Interobserver intra class correlation (ICC) was 0.70 in-plane and 0.78 through-plane. Intra-observer ICC for three observers was 0.76, 0.81, 0.83, showing moderate to good reliability on this small cohort.

CONCLUSIONS

This in-vivo feasibility study suggests distortion due to metallic hip prosthesis is not an obstacle for pelvic radiotherapy on an MR Linac. Research on the impact on plan quality is warranted.

ADVANCES IN KNOWLEDGE

This work supports feasibility of treating patients with metallic hip prosthesis on an MR Linac.

Assessing the impact of distortion correction on Gamma Knife radiosurgery for multiple metastasis: Volumetric and dosimetric analysis

Introduction

Magnetic resonance imaging (MRI) is a robust neuroimaging technique and is the preferred method for stereotactic radiosurgery (SRS) planning. However, MRI data always contain distortions caused by hardware and patient factors.

Research question

Can these distortions potentially compromise the effectiveness and safety of SRS treatments?

Material and methods

Twenty-six MR datasets with multiple metastatic brain tumors (METs) used for Gamma Knife radiosurgery (GKRS) were retrospectively evaluated. A commercially available software was used for distortion correction. Geometrical agreement between corrected and uncorrected tumor volumes was evaluated using MacDonald criteria, Euclidian distance, and Dice similarity coefficient (DSC). SRS plans were generated using uncorrected tumor volumes, which were assessed to determine their coverage of the corrected tumor volumes.

Results

The median target volume was 0.38 cm3 (range,0.01-12.38 cm3). A maximum displacement of METs of up to 2.87 mm and a median displacement of 0.55 mm (range,0.1-2.87 mm) were noted. The median DSC between uncorrected and corrected MRI was 0.92, and the most concerning case had a DSC of 0.46. Although all plans met the optimization criterion of at least 98% of the uncorrected tumor volume (median 99.55%, range 98.1-100%) receiving at least 100% of the prescription dose, the percent of the corrected tumor volume receiving the total prescription dose was a median of 95.45% (range,23.1-99.5%).

Discussion and conclusion

MRI distortion, though visually subtle, has significant implications for SRS planning. Regular utilization of corrected MRI is recommended for SRS planning as distortion is sometimes enough to cause a volumetric miss of SRS targets.

Commissioning and implementing a Quality Assurance program for dedicated radiation oncology MRI scanners

PURPOSE

ACR and AAPM task group's guidelines addressing commissioning for dedicated MR simulators were recently published. The goal of the current paper is to present the authors' 2-year experience regarding the commissioning and introduction of a QA program based on these guidelines and an associated automated workflow.

METHODS

All mandatory commissioning tests suggested by AAPM report 284 were performed and results are reported for two MRI scanners (MAGNETOM Sola and Aera). Visual inspection, vendor clinical or service platform, third-party software, or in-house python-based code were used. Automated QA and data analysis was performed via vendor, in-house or third-party software. QATrack+ was used for QA data logging and storage. 3D geometric distortion, B0 inhomogeneity, EPI, and parallel imaging performance were evaluated.

RESULTS

Contrasting with AAPM report 284 recommendations, homogeneity and RF tests were performed monthly. The QA program allowed us to detect major failures over time (shimming, gradient calibration and RF interference). Automated QA, data analysis, and logging allowed fast ACR analysis daily and monthly QA to be performed in 3 h. On the Sola, the average distortion is 1 mm for imaging radii of 250 mm or less. For radii of up to 200 mm, the maximum, average (standard deviation) distortion is 1.2  and 0.4 mm (0.3 mm). Aera values are roughly double the Sola for radii up to 200 mm. EPI geometric distortion, ghosting ratio, and long-term stability were found to be under the maximum recommended values. Parallel imaging SNR ratio was stable and close to the theoretical value (ideal g-factor). No major failures were detected during commissioning.

CONCLUSION

An automated workflow and enhanced QA program allowed to automatically track machine and environmental changes over time and to detect periodic failures and errors that might otherwise have gone unnoticed. The Sola is more geometrically accurate, with a more homogenous B0 field than the Aera.

mcTFI QSM MRI ABC/2 intracranial hemorrhage to noncontrast head CT volume measurement equivalence

BACKGROUND/OBJECTIVES

Intracranial hemorrhage (ICH) volume assessment is an important part of patient management and is routinely obtained by non-contrast head CT (NCHCT) using the validated ABC/2 measurement method. Because conventional MRI imaging sequences demonstrate variability in ICH appearance, volumetric analyses for MRI bleed volume in a standardized manner using ABC/2 is not possible. The recently introduced multiecho-complex total field inversion quantitative susceptibility mapping (mcTFI QSM) MRI technique, which maps brain tissue susceptibility to both depict brain tissue structures and quantify tissue susceptibility, may provide a viable alternative. In this study we evaluated mcTFI QSM ABC/2 ICH volume assessment relative to NCHCT.

METHODS

Patients with ICH who had undergone NCHCT and MRI brain scans within 48 h were recruited for this retrospective study. The ABC/2 method was applied to estimate the bleed volume for both NCHCT and MRI by a CAQ-certified neuroradiologist with 10 years of experience and a trained laboratory assistant. Results were analyzed via Bland-Altman (B-A) and linear regression.

RESULTS

54 patients (27 females) who had undergone NCHCT and MRI within 48 h (<24 h., n = 31, 24-48 h, n = 10) were enrolled. mcTFI QSM ICH volume measurement method showed a positive correlation (99.5%) compared to NCHCT. B-A plot comparing ABC/2 ICH volume on NCHCT and mcTFI MRI done for patients within 24 h demonstrates a bias of -0.09%.

CONCLUSIONS

ICH volume calculation using ABC/2 on mcTFI QSM showed a high correlation with NCHCT measurement. These results suggest mcTFI QSM is a promising MRI method for ABC/2 for bleed volume measurement.

MRI scarcity in low- and middle-income countries

Since the introduction of MRI as a sustainable diagnostic modality, global accessibility to its services has revealed a wide discrepancy between populations-leaving most of the population in LMICs without access to this important imaging modality. Several factors lead to the scarcity of MRI in LMICs; for example, inadequate infrastructure and the absence of a dedicated workforce are key factors in the scarcity observed. RAD-AID has contributed to the advancement of radiology globally by collaborating with our partners to make radiology more accessible for medically underserved communities. However, progress is slow and further investment is needed to ensure improved global access to MRI.

Use of a Commercial 7-T MRI Scanner for Clinical Brain Imaging: Indications, Protocols, Challenges, and Solutions-A Single-Center Experience

The first commercially available 7-T MRI scanner (Magnetom Terra) was approved by the FDA in 2017 for clinical imaging of the brain and knee. After initial protocol development and sequence optimization efforts in volunteers, the 7-T system, in combination with an FDA-approved 1-channel transmit/32-channel receive array head coil, can now be routinely used for clinical brain MRI examinations. The ultrahigh field strength of 7-T MRI has the advantages of improved spatial resolution, increased SNR, and increased CNR but also introduces an array of new technical challenges. The purpose of this article is to describe an institutional experience with the use of the commercially available 7-T MRI scanner for routine clinical brain imaging. Specific clinical indications for which 7-T MRI may be useful for brain imaging include brain tumor evaluation with possible perfusion imaging and/or spectroscopy, radiotherapy planning; evaluation of multiple sclerosis and other demyelinating diseases, evaluation of Parkinson disease and guidance of deep brain stimulator placement, high-detail intracranial MRA and vessel wall imaging, evaluation of pituitary pathology, and evaluation of epilepsy. Detailed protocols, including sequence parameters, for these various indications are presented, and implementation challenges (including artifacts, safety, and side effects) and potential solutions are explored.

MRI sequence optimisation methods to identify cranial nerve course for radiotherapy planning

INTRODUCTION

Magnetic resonance imaging (MRI) is being increasingly used to improve radiation therapy planning by allowing visualisation of organs at risk that cannot be well-defined on computed tomography (CT). Diagnostic sequences are increasingly being adapted for radiation therapy planning, such as the use of heavily T2-weighted 3D SPACE (Sampling Perfection with Application optimised Contrasts using different flip angle Evolution) sequence for cranial nerve identification in head and neck tumour treatment planning.

METHODS

A 3D isotropic T2 SPACE sequence used for cranial nerve identification was adapted for radiation therapy purposes. Distortion was minimised using a spin-echo-based sequence, 3D distortion correction, isocentre scanning and an increased readout bandwidth. Radiation therapy positioning was accounted for by utilising two small flex, 4-channel coils. The protocol was validated for cranial nerve identification in clinical applications and distortion minimisation using an MRI QA phantom.

RESULTS

Normal anatomy of the cranial nerves CI-CIX, were presented, along with a selection of clinical applications and abnormal anatomy. The usefulness of cranial nerve identification is discussed for several case studies, particularly in proximity to tumours extending into the base of skull region. In-house testing validated that higher bandwidths of 600 Hz resulted in minimal displacement well below 1 mm.

CONCLUSION

The use of MRI for radiation therapy planning allows for greater individualisation and prediction of patient outcomes. Dose reduction to cranial nerves can decrease late side effects such as cranial neuropathy. In addition to current applications, future directions include further applications of this technology for radiation therapy treatments.

Geometric distortions in clinical MRI sequences for radiotherapy: insights gained from a multicenter investigation

BACKGROUND

As magnetic resonance imaging (MRI) becomes increasingly integrated into radiotherapy (RT) for enhanced treatment planning and adaptation, the inherent geometric distortion in acquired MR images pose a potential challenge to treatment accuracy. This study aimed to evaluate the geometric distortion levels in the clinical MRI protocols used across Danish RT centers and discuss influence of specific sequence parameters. Based on the variety in geometric performance across centers, we assess if harmonization of MRI sequences is a relevant measure.

MATERIALS AND METHODS

Nine centers participated with 12 MRI scanners and MRI-Linacs (MRL). Using a travelling phantom approach, a reference MRI sequence was used to assess variation in baseline distortion level between scanners. The phantom was also scanned with local clinical MRI sequences for brain, head/neck (H/N), abdomen, and pelvis. The influence of echo time, receiver bandwidth, image weighting, and 2D/3D acquisition was investigated.

RESULTS

We found a large variation in geometric accuracy across 93 clinical sequences examined, exceeding the baseline variation found between MRI scanners (σ = 0.22 mm), except for abdominal sequences where the variation was lower. Brain and abdominal sequences showed lowest distortion levels ([0.22, 2.26] mm), and a large variation in performance was found for H/N and pelvic sequences ([0.19, 4.07] mm). Post hoc analyses revealed that distortion levels decreased with increasing bandwidth and a less clear increase in distortion levels with increasing echo time. 3D MRI sequences had lower distortion levels than 2D (median of 1.10 and 2.10 mm, respectively), and in DWI sequences, the echo-planar imaging read-out resulted in highest distortion levels.

CONCLUSION

There is a large variation in the geometric distortion levels of clinical MRI sequences across Danish RT centers, and between anatomical sites. The large variation observed makes harmonization of MRI sequences across institutions and adoption of practices from well-performing anatomical sites, a relevant measure within RT.

Treatment of Central Nervous System Tumors on Combination MR-Linear Accelerators: Review of Current Practice and Future Directions

Magnetic resonance imaging (MRI) provides excellent visualization of central nervous system (CNS) tumors due to its superior soft tissue contrast. Magnetic resonance-guided radiotherapy (MRgRT) has historically been limited to use in the initial treatment planning stage due to cost and feasibility. MRI-guided linear accelerators (MRLs) allow clinicians to visualize tumors and organs at risk (OARs) directly before and during treatment, a process known as online MRgRT. This novel system permits adaptive treatment planning based on anatomical changes to ensure accurate dose delivery to the tumor while minimizing unnecessary toxicity to healthy tissue. These advancements are critical to treatment adaptation in the brain and spinal cord, where both preliminary MRI and daily CT guidance have typically had limited benefit. In this narrative review, we investigate the application of online MRgRT in the treatment of various CNS malignancies and any relevant ongoing clinical trials. Imaging of glioblastoma patients has shown significant changes in the gross tumor volume over a standard course of chemoradiotherapy. The use of adaptive online MRgRT in these patients demonstrated reduced target volumes with cavity shrinkage and a resulting reduction in radiation dose to uninvolved tissue. Dosimetric feasibility studies have shown MRL-guided stereotactic radiotherapy (SRT) for intracranial and spine tumors to have potential dosimetric advantages and reduced morbidity compared with conventional linear accelerators. Similarly, dosimetric feasibility studies have shown promise in hippocampal avoidance whole brain radiotherapy (HA-WBRT). Next, we explore the potential of MRL-based multiparametric MRI (mpMRI) and genomically informed radiotherapy to treat CNS disease with cutting-edge precision. Lastly, we explore the challenges of treating CNS malignancies and special limitations MRL systems face.

Assessment of dual-energy computed tomography derived virtual monoenergetic imaging for target volume delineation of brain metastases

BACKGROUND

Objective and subjective assessment of image quality of brain metastases on dual-energy computed tomography (DECT) virtual monoenergetic imaging (VMI) and its impact on target volume delineation.

MATERIALS AND METHODS

26 patients with 37 brain metastases receiving Magnetic Resonance Imaging (MRI) and DECT for stereotactic radiotherapy planning were included in this retrospective analysis. Lesion contrast (LC), contrast-to-noise ratio (CNR) and signal-to-noise ratio (SNR) were assessed for reconstructed VMI at 63 keV and artificial 120 kV Computed Tomography (CT). Image contrast and demarcation of metastases between 120 kV CT, VMI and MRI were subjectively assessed. Brain metastases were delineated by four radiation oncologists on VMI with a fixed or free brain window and contours were compared to solely MRI-based delineation using the Dice similarity coefficient.

RESULTS

LC, CNR and SNR were significantly higher in VMI than in 120 kV CT (p < 0.0001). Image contrast and lesion demarcation were significantly better on VMI compared to 120 kV CT (p < 0.0001). Mean gross tumor volume (GTV)/planning target volume (PTV) Dice similarity coefficients were 0.87/0.9 for metastases without imaging uncertainties (no artifacts, calcification or impaired visibility with MRI) but worse for metastases with imaging uncertainties (0.71/0.74). Target volumes delineated on VMI were around 5-10% smaller compared to MRI.

CONCLUSION

Image quality of VMI is objectively and subjectively superior to conventional CT. VMI provides significant advantages in stereotactic radiotherapy planning with improved visibility of brain metastases and geometrically distortion-free representation of brain metastases. Beside a plausibility check of MRI-based target volume delineation, VMI might improve reliability and accuracy in target volume definition particularly in cases with imaging uncertainties with MRI.

Accurate, repeatable, and geometrically precise diffusion-weighted imaging on a 0.35 T magnetic resonance imaging-guided linear accelerator

Background and purpose

Diffusion weighted imaging (DWI) allows for the interrogation of tissue cellularity, which is a surrogate for cellular proliferation. Previous attempts to incorporate DWI into the workflow of a 0.35 T MR-linac (MRL) have lacked quantitative accuracy. In this study, accuracy, repeatability, and geometric precision of apparent diffusion coefficient (ADC) maps produced using an echo planar imaging (EPI)-based DWI protocol on the MRL system is illustrated, and in vivo potential for longitudinal patient imaging is demonstrated.

Materials and methods

Accuracy and repeatability were assessed by measuring ADC values in a diffusion phantom at three timepoints and comparing to reference ADC values. System-dependent geometric distortion was quantified by measuring the distance between 93 pairs of phantom features on ADC maps acquired on a 0.35 T MRL and a 3.0 T diagnostic scanner and comparing to spatially precise CT images. Additionally, for five sarcoma patients receiving radiotherapy on the MRL, same-day in vivo ADC maps were acquired on both systems, one of which at multiple timepoints.

Results

Phantom ADC quantification was accurate on the 0.35 T MRL with significant discrepancies only seen at high ADC. Average geometric distortions were 0.35 (±0.02) mm and 0.85 (±0.02) mm in the central slice and 0.66 (±0.04) mm and 2.14 (±0.07) mm at 5.4 cm off-center for the MRL and diagnostic system, respectively. In the sarcoma patients, a mean pretreatment ADC of 910x10-6 (±100x10-6) mm2/s was measured on the MRL.

Conclusions

The acquisition of accurate, repeatable, and geometrically precise ADC maps is possible at 0.35 T with an EPI approach.

Comprehensive MR imaging QA of 0.35 T MR-Linac using a multi-purpose large FOV phantom: A single-institution experience

PURPOSE

Magnetic resonance-guided radiotherapy (MRgRT) is desired for the treatment of diseases in the abdominothoracic region, which has a broad imaging area and continuous motion. To ensure accurate treatment delivery, an effective image quality assurance (QA) program, with a phantom that covers the field of view (FOV) similar to a human torso, is required. However, routine image QA for a large FOV is not readily available at many MRgRT centers. In this work, we present the clinical experience of the large FOV MRgRT Insight phantom for periodic daily and monthly comprehensive magnetic resonance imaging (MRI)-QA and its feasibility compared to the existing institutional routine MRI-QA procedures in 0.35 T MRgRT.

METHODS

Three phantoms; ViewRay cylindrical water phantom, Fluke 76-907 uniformity and linearity phantom, and Modus QA large FOV MRgRT Insight phantom, were imaged on the 0.35 T MR-Linac. The measurements were made in MRI mode with the true fast imaging with steady-state free precession (TRUFI) sequence. The ViewRay cylindrical water phantom was imaged in a single-position setup whereas the Fluke phantom and Insight phantom were imaged in three different orientations: axial, sagittal, and coronal. Additionally, the phased array coil QA was performed using the horizontal base plate of the Insight phantom by placing the desired coil around the base section which was compared to an in-house built Polyurethane foam phantom for reference.

RESULT

The Insight phantom captured image artifacts across the entire planar field of view, up to 400 mm, in a single image acquisition, which is beyond the FOV of the conventional phantoms. The geometric distortion test showed a similar distortion of 0.45 ± 0.01  and 0.41 ± 0.01 mm near the isocenter, that is, within 300 mm lengths for Fluke and Insight phantoms, respectively, but showed higher geometric distortion of 0.8 ± 0.4 mm in the peripheral region between 300 and 400 mm of the imaging slice for the Insight phantom. The Insight phantom with multiple image quality features and its accompanying software utilized the modulation transform function (MTF) to evaluate the image spatial resolution. The average MTF values were 0.35 ± 0.01, 0.35 ± 0.01, and 0.34 ± 0.03 for axial, coronal, and sagittal images, respectively. The plane alignment and spatial accuracy of the ViewRay water phantom were measured manually. The phased array coil test for both the Insight phantom and the Polyurethane foam phantoms ensured the proper functionality of each coil element.

CONCLUSION

The multifunctional large FOV Insight phantom helps in tracking MR imaging quality of the system to a larger extent compared to the routine daily and monthly QA phantoms currently used in our institute. Also, the Insight phantom is found to be more feasible for routine QA with easy setup.

Off-resonance artifact correction for MRI: A review

In magnetic resonance imaging (MRI), inhomogeneity in the main magnetic field used for imaging, referred to as off-resonance, can lead to image artifacts ranging from mild to severe depending on the application. Off-resonance artifacts, such as signal loss, geometric distortions, and blurring, can compromise the clinical and scientific utility of MR images. In this review, we describe sources of off-resonance in MRI, how off-resonance affects images, and strategies to prevent and correct for off-resonance. Given recent advances and the great potential of low-field and/or portable MRI, we also highlight the advantages and challenges of imaging at low field with respect to off-resonance.

Quality assurance assessment of intra-acquisition diffusion-weighted and T2-weighted magnetic resonance imaging registration and contour propagation for head and neck cancer radiotherapy

BACKGROUND/PURPOSE

Adequate image registration of anatomical and functional magnetic resonance imaging (MRI) scans is necessary for MR-guided head and neck cancer (HNC) adaptive radiotherapy planning. Despite the quantitative capabilities of diffusion-weighted imaging (DWI) MRI for treatment plan adaptation, geometric distortion remains a considerable limitation. Therefore, we systematically investigated various deformable image registration (DIR) methods to co-register DWI and T2-weighted (T2W) images.

MATERIALS/METHODS

We compared three commercial (ADMIRE, Velocity, Raystation) and three open-source (Elastix with default settings [Elastix Default], Elastix with parameter set 23 [Elastix 23], Demons) post-acquisition DIR methods applied to T2W and DWI MRI images acquired during the same imaging session in twenty immobilized HNC patients. In addition, we used the non-registered images (None) as a control comparator. Ground-truth segmentations of radiotherapy structures (tumour and organs at risk) were generated by a physician expert on both image sequences. For each registration approach, structures were propagated from T2W to DWI images. These propagated structures were then compared with ground-truth DWI structures using the Dice similarity coefficient and mean surface distance.

RESULTS

19 left submandibular glands, 18 right submandibular glands, 20 left parotid glands, 20 right parotid glands, 20 spinal cords, and 12 tumours were delineated. Most DIR methods took <30 s to execute per case, with the exception of Elastix 23 which took ∼458 s to execute per case. ADMIRE and Elastix 23 demonstrated improved performance over None for all metrics and structures (Bonferroni-corrected p < 0.05), while the other methods did not. Moreover, ADMIRE and Elastix 23 significantly improved performance in individual and pooled analysis compared to all other methods.

CONCLUSIONS

The ADMIRE DIR method offers improved geometric performance with reasonable execution time so should be favoured for registering T2W and DWI images acquired during the same scan session in HNC patients. These results are important to ensure the appropriate selection of registration strategies for MR-guided radiotherapy.

Stereotactic Magnetic Resonance-Guided Adaptive and Non-Adaptive Radiotherapy on Combination MR-Linear Accelerators: Current Practice and Future Directions

Stereotactic body radiotherapy (SBRT) is an effective radiation therapy technique that has allowed for shorter treatment courses, as compared to conventionally dosed radiation therapy. As its name implies, SBRT relies on daily image guidance to ensure that each fraction targets a tumor, instead of healthy tissue. Magnetic resonance imaging (MRI) offers improved soft-tissue visualization, allowing for better tumor and normal tissue delineation. MR-guided RT (MRgRT) has traditionally been defined by the use of offline MRI to aid in defining the RT volumes during the initial planning stages in order to ensure accurate tumor targeting while sparing critical normal tissues. However, the ViewRay MRIdian and Elekta Unity have improved upon and revolutionized the MRgRT by creating a combined MRI and linear accelerator (MRL), allowing MRgRT to incorporate online MRI in RT. MRL-based MR-guided SBRT (MRgSBRT) represents a novel solution to deliver higher doses to larger volumes of gross disease, regardless of the proximity of at-risk organs due to the (1) superior soft-tissue visualization for patient positioning, (2) real-time continuous intrafraction assessment of internal structures, and (3) daily online adaptive replanning. Stereotactic MR-guided adaptive radiation therapy (SMART) has enabled the safe delivery of ablative doses to tumors adjacent to radiosensitive tissues throughout the body. Although it is still a relatively new RT technique, SMART has demonstrated significant opportunities to improve disease control and reduce toxicity. In this review, we included the current clinical applications and the active prospective trials related to SMART. We highlighted the most impactful clinical studies at various tumor sites. In addition, we explored how MRL-based multiparametric MRI could potentially synergize with SMART to significantly change the current treatment paradigm and to improve personalized cancer care.

Rapid distortion correction enables accurate magnetic resonance imaging-guided real-time adaptive radiotherapy

Background and purpose

Magnetic resonance imaging (MRI)-Linac systems combine simultaneous MRI with radiation delivery, allowing treatments to be guided by anatomically detailed, real-time images. However, MRI can be degraded by geometric distortions that cause uncertainty between imaged and actual anatomy. In this work, we develop and integrate a real-time distortion correction method that enables accurate real-time adaptive radiotherapy.

Materials and methods

The method was based on the pre-treatment calculation of distortion and the rapid correction of intrafraction images. A motion phantom was set up in an MRI-Linac at isocentre (P₀ ), the edge (P ₁) and just outside (P ₂) the imaging volume. The target was irradiated and tracked during real-time adaptive radiotherapy with and without the distortion correction. The geometric tracking error and latency were derived from the measurements of the beam and target positions in the EPID images.

Results

Without distortion correction, the mean geometric tracking error was 1.3 mm at P ₁ and 3.1 mm at P ₂. When distortion correction was applied, the error was reduced to 1.0 mm at P ₁ and 1.1 mm at P ₂. The corrected error was similar to an error of 0.9 mm at P₀ where the target was unaffected by distortion indicating that this method has accurately accounted for distortion during tracking. The latency was 319 ± 12 ms without distortion correction and 335 ± 34 ms with distortion correction.

Conclusions

We have demonstrated a real-time distortion correction method that maintains accurate radiation delivery to the target, even at treatment locations with large distortion.

Empirical planning target volume modeling for high precision MRI guided intracranial radiotherapy

Purpose

Magnetic resonance image-guided radiotherapy for intracranial indications is a promising advance; however, uncertainties remain for both target localization after translation-only MR setup and intrafraction motion. This investigation quantified these uncertainties and developed a population-based planning target volume (PTV) model to explore target and organ-at-risk (OAR) volumetric coverage tradeoffs.

Methods

Sixty-six patients, 49 with a primary brain tumor and 17 with a post-surgical resection cavity, treated on a 1.5T-based MR-linac across 1329 fractions were included. At each fraction, patients were setup by translation-only fusion of the online T1 MRI to the planning image. Each fusion was independently repeated offline accounting for rotations. The six degree-of-freedom difference between fusions was applied to transform the planning CTV at each fraction (CTV_fx). A PTV model parameterized by volumetric CTV_fx coverage, proportion of fractions, and proportion of patients was developed. Intrafraction motion was quantified in a 412 fraction subset as the fusion difference between post- and pre-irradiation T1 MRIs.

Results

For the left-right/anterior-posterior/superior-inferior axes, mean ± SD of the rotational fusion differences were 0.1 ± 0.8/0.1 ± 0.8/-0.2 ± 0.9°. Covering 98 % of the CTV_fx in 95 % of fractions in 95 % of patients required a 3 mm PTV margin. Margin reduction decreased PTV-OAR overlap; for example, the proportion of optic chiasm overlapped by the PTV was reduced up to 23.5 % by margin reduction from 4 mm to 3 mm.

Conclusions

An evidence-based PTV model was developed for brain cancer patients treated on the MR-linac. Informed by this model, we have clinically adopted a 3 mm PTV margin for conventionally fractionated intracranial patients.

Assessment of MRI image distortion based on 6 consecutive years of annual QAs and measurements on 14 MRI scanners used for radiation therapy

PURPOSE

To determine the magnitude of MRI image distortion based on 6 consecutive years of annual quality assurances/measurements on 14 MRI scanners used for radiation therapy and to provide evidence for the inclusion of additional margin for treatment planning.

METHODS AND MATERIALS

We used commercial MRI image phantoms to quantitatively study the MRI image distortion over period of 6 years for up to 14 1.5 and 3 T MRI scanners that could potentially be used to provide MRI images for treatment planning. With the phantom images collected from 2016 to 2022, we investigated the MRI image distortion, the dependence of distortion on the distance from the imaging isocenter, and the possible causes of large distortion discovered.

RESULTS

MRI image distortion increases with the distance from the imaging isocenter. For a region of interest (ROI) with a radius of 100 mm centered at the isocenter, the mean magnitude of distortion for all MRI scanners is 0.44 ± 0.18 mm $0.44 \pm 0.18\;{\rm{mm}}$ , and the maximum distortion varies from 0.52 to 1.31 mm $0.52\;{\rm{to}}\;1.31\;{\rm{mm}}$ depending on MRI scanners. For an ROI with a radius of 200 mm centered at the isocenter, the mean magnitude of distortion increases to 0.84 ± 0.45 mm $0.84 \pm 0.45\;{\rm{mm}}$ , and the range of the maximum distortion increases to 1.92 - 5.03 mm $1.92 - 5.03\;{\rm{mm}}$ depending on MRI scanners. The distortion could reach 2 mm at 150 mm from the isocenter.

CONCLUSION

An additional margin to accommodate image distortion should be considered for treatment planning. Imaging with proper patient alignment to the isocenter is vital to reducing image distortion. We recommend performing image distortion checks annually and after major upgrade on MRI scanners.

Generation of synthetic megavoltage CT for MRI-only radiotherapy treatment planning using a 3D deep convolutional neural network

BACKGROUND

Megavoltage computed tomography (MVCT) has been implemented on many radiotherapy treatment machines for on-board anatomical visualization, localization, and adaptive dose calculation. Implementing an MR-only workflow by synthesizing MVCT from magnetic resonance imaging (MRI) would offer numerous advantages for treatment planning and online adaptation.

PURPOSE

In this work, we sought to synthesize MVCT (sMVCT) datasets from MRI using deep learning to demonstrate the feasibility of MRI-MVCT only treatment planning.

METHODS

MVCTs and T1-weighted MRIs for 120 patients treated for head-and-neck cancer were retrospectively acquired and co-registered. A deep neural network based on a fully-convolutional 3D U-Net architecture was implemented to map MRI intensity to MVCT HU. Input to the model were volumetric patches generated from paired MRI and MVCT datasets. The U-Net was initialized with random parameters and trained on a mean absolute error (MAE) objective function. Model accuracy was evaluated on 18 withheld test exams. sMVCTs were compared to respective MVCTs. Intensity-modulated volumetric radiotherapy (IMRT) plans were generated on MVCTs of four different disease sites and compared to plans calculated onto corresponding sMVCTs using the gamma metric and dose-volume-histograms (DVHs).

RESULTS

MAE values between sMVCT and MVCT datasets were 93.3 ± 27.5, 78.2 ± 27.5, and 138.0 ± 43.4 HU for whole body, soft tissue, and bone volumes, respectively. Overall, there was good agreement between sMVCT and MVCT, with bone and air posing the greatest challenges. The retrospective dataset introduced additional deviations due to sinus filling or tumor growth/shrinkage between scans, differences in external contours due to variability in patient positioning, or when immobilization devices were absent from diagnostic MRIs. Dose distributions of IMRT plans evaluated for four test cases showed close agreement between sMVCT and MVCT images when evaluated using DVHs and gamma dose metrics, which averaged to 98.9 ± 1.0% and 96.8 ± 2.6% analyzed at 3%/3 mm and 2%/2 mm, respectively.

CONCLUSIONS

MVCT datasets can be generated from T1-weighted MRI using a 3D deep convolutional neural network with dose calculation on a sample sMVCT in close agreement with the MVCT. These results demonstrate the feasibility of using MRI-derived sMVCT in an MR-only treatment planning workflow.

Practical guidelines of online MR-guided adaptive radiotherapy

The first magnetic resonance (MR)-guided radiotherapy system in Japan was installed in May 2017. Implementation of online MR-guided adaptive radiotherapy (MRgART) began in February 2018. Online MRgART offers greater treatment accuracy owing to the high soft-tissue contrast in MR-images (MRI), compared to that in X-ray imaging. The Japanese Society for Magnetic Resonance in Medicine (JSMRM), Japan Society of Medical Physics (JSMP), Japan Radiological Society (JRS), Japanese Society of Radiological Technology (JSRT), and Japanese Society for Radiation Oncology (JASTRO) jointly established the comprehensive practical guidelines for online MRgART. These guidelines propose the essential requirements for clinical implementation of online MRgART with respect to equipment, personnel, institutional environment, practice guidance, and quality assurance/quality control (QA/QC). The minimum requirements for related equipment and QA/QC tools, recommendations for safe operation of MRI system, and the implementation system are described. The accuracy of monitor chamber and detector in dose measurements should be confirmed because of the presence of magnetic field. The ionization chamber should be MR-compatible. Non-MR-compatible devices should be used in an area that is not affected by the static magnetic field (outside the five Gauss line), and their operation should be checked to ensure that they do not affect the MR image quality. Dose verification should be performed using an independent dose verification system that has been confirmed to be reliable through commissioning. This guideline proposes the checklists to ensure the safety of online MRgART. Successful clinical implementation of online MRgART requires close collaboration between physician, radiological technologist, nurse, and medical physicist.

Evaluation of the ability of the Brainlab Elements Cranial Distortion Correction algorithm to correct clinically relevant MRI distortions for cranial SRT

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.

Integrated MRI-guided radiotherapy - opportunities and challenges

MRI can help to categorize tissues as malignant or non-malignant both anatomically and functionally, with a high level of spatial and temporal resolution. This non-invasive imaging modality has been integrated with radiotherapy in devices that can differentially target the most aggressive and resistant regions of tumours. The past decade has seen the clinical deployment of treatment devices that combine imaging with targeted irradiation, making the aspiration of integrated MRI-guided radiotherapy (MRIgRT) a reality. The two main clinical drivers for the adoption of MRIgRT are the ability to image anatomical changes that occur before and during treatment in order to adapt the treatment approach, and to image and target the biological features of each tumour. Using motion management and biological targeting, the radiation dose delivered to the tumour can be adjusted during treatment to improve the probability of tumour control, while simultaneously reducing the radiation delivered to non-malignant tissues, thereby reducing the risk of treatment-related toxicities. The benefits of this approach are expected to increase survival and quality of life. In this Review, we describe the current state of MRIgRT, and the opportunities and challenges of this new radiotherapy approach.

Time-course assessment of 3D-image distortion on the 1.5 T Marlin/Elekta Unity MR-LINAC

PURPOSE

The efficacy of MR-guided radiotherapy on a MR-LINAC (MR-L) is dependent on the geometric accuracy of its MR images over clinically relevant Fields-of-View (FOVs). Our objectives were to: evaluate gradient non-linearity (GNL) on the Elekta Unity MR-L across time via 76 weekly measurements of 3D-distortion over concentrically larger diameter spherical volumes (DSVs); quantify distortion measurement error; and assess the temporal stability of spatial distortion using statistical process control (SPC).

METHODS

MR-image distortion was assessed using a large-FOV 3D-phantom containing 1932 markers embedded in seven parallel plates, spaced 25 mm × 25 mm in- and 55 mm through-plane. Automatically analyzed T1 images yielded distortions in 200, 300, 400 and 500 mm concentric DSVs. Distortion measurement error was evaluated using median absolute difference analysis of imaging repeatability tests.

RESULTS

Over the measurement period absolute time-averaged distortion varied between: dr = 0.30 - 0.49 mm, 0.53 - 0.80 mm, 1.0 - 1.4 mm and 2.28 - 2.37 mm, for DSVs 200, 300, 400 and 500 mm at the 98^th percentile level. Repeatability tests showed that imaging/repositioning introduces negligible error: mean ≤ 0.02 mm (max ≤ 0.3 mm). SPC analysis showed image distortion was stable across all DSVs; however, noticeable changes in GNL were observed following servicing at the one-year mark.

CONCLUSIONS

Image distortion on the MR-L is in the sub-millimeter range for DSVs ≤ 300 mm and stable across time, with SPC analysis indicating all measurements remain within control for each DSV.

Increased brain volumetric measurement precision from multi-site 3D T1-weighted 3 T magnetic resonance imaging by correcting geometric distortions

PURPOSE

Magnetic resonance imaging (MRI) scanner-specific geometric distortions may contribute to scanner induced variability and decrease volumetric measurement precision for multi-site studies. The purpose of this study was to determine whether geometric distortion correction increases the precision of brain volumetric measurements in a multi-site multi-scanner study.

METHODS

Geometric distortion variation was quantified over a one-year period at 10 sites using the distortion fields estimated from monthly 3D T1-weighted MRI geometrical phantom scans. The variability of volume and distance measurements were quantified using synthetic volumes and a standard quantitative MRI (qMRI) phantom. The effects of geometric distortion corrections on MRI derived volumetric measurements of the human brain were assessed in two subjects scanned on each of the 10 MRI scanners and in 150 subjects with cerebrovascaular disease (CVD) acquired across imaging sites.

RESULTS

Geometric distortions were found to vary substantially between different MRI scanners but were relatively stable on each scanner over a one-year interval. Geometric distortions varied spatially, increasing in severity with distance from the magnet isocenter. In measurements made with the qMRI phantom, the geometric distortion correction decreased the standard deviation of volumetric assessments by 35% and distance measurements by 42%. The average coefficient of variance decreased by 16% in gray matter and white matter volume estimates in the two subjects scanned on the 10 MRI scanners.

CONCLUSION

Geometric distortion correction using an up-to-date correction field is recommended to increase precision in volumetric measurements made from MRI images.

Geometric distortion assessment in 3T MR images used for treatment planning in cranial Stereotactic Radiosurgery and Radiotherapy

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.

Developing quality assurance tests for simultaneous Positron Emission Tomography – Magnetic Resonance imaging for radiotherapy planning

Background and purpose Simultaneous Positron Emission Tomography – Magnetic Resonance (PET-MR) imaging can potentially improve radiotherapy by enabling more accurate tumour delineation and dose painting. The use of PET-MR imaging for radiotherapy planning requires a comprehensive Quality Assurance (QA) programme to be developed. This study aimed to develop the QA tests required and assess their repeatability and stability. Materials and methods QA tests were developed for: MR image quality, MR geometric accuracy, electromechanical accuracy, PET-MR alignment accuracy, Diffusion Weighted (DW)-MR Apparent Diffusion Coefficient (ADC) accuracy and PET Standard Uptake Value (SUV) accuracy. Each test used a dedicated phantom and was analysed automatically or semi-automatically, with in–house software. Repeatability was evaluated by three same-day measurements with independent phantom positions. Stability was assessed through 12 monthly measurements. Results The repeatability Standard Deviations (SDs) of distortion for the MR geometric accuracy test were ⩽0.7mm. The repeatability SDs in ADC difference from reference were ⩽3% for the DW-MR accuracy test. The PET SUV difference from reference repeatability SD was 0.3%. The stability SDs agreed within 0.6mm, 1 percentage point and 1.4 percentage points of the repeatability SDs for the geometric, ADC and SUV accuracy tests respectively. There were no monthly trends apparent. These results were representative of the other tests. Conclusions QA Tests for radiotherapy planning PET-MR have been developed. The tests appeared repeatable and stable over a 12-month period. The developed QA tests could form the basis of a QA programme that enables high-quality, robust PET-MR imaging for radiotherapy 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.

Impact of geometric distortion on dose deviation for photon and proton treatment plans

We investigated the dose deviation related to geometric distortion and dose gradient on magnetic resonance‐only treatment planning for intensity‐modulated radiation therapy and proton therapy. The residual geometric distortion of two different magnetic resonance imaging (MRI) sequences (A) and (B) was applied in the computed tomography image and the structure set of each patient through a polynomial MRI geometric distortion model to simulate MRI‐based treatment planning. A 3D histogram was generated to specify the relationship of dose deviation to geometric distortion and dose gradient. When the dose gradient (G_d) approached zero, the maximum dose deviation reached 1.64% and 2.71% for photon plans of sequences A and B, respectively. For proton plans, the maximum dose deviation reached 3.15% and 4.89% for sequences A and B, respectively. When the geometric distortion (d) was close to zero, the maximum dose deviation was less than 0.8% for photon and proton plans of both sequences. Under extreme conditions (d = 2 mm and G_d = 4.5%/mm), the median value of dose deviation reached 3% and 3.49% for photon and proton plans, respectively for sequence A, and 2.93% and 4.55% for photon and proton plans, respectively, for sequence B. We demonstrate that the dose deviation is specific to MRI hardware parameters. Compared to the photon plan, the proton plan is more sensitive to the changes in geometric distortion. For typical clinical MRI geometric distortion (d ≤2 mm), the median dose deviation is expected to be within 3% and 5% for photon and proton plans, respectively.

A narrative review of MRI acquisition for MR-guided-radiotherapy in prostate cancer

Magnetic resonance guided radiotherapy (MRgRT), enabled by the clinical introduction of the integrated MRI and linear accelerator (MR-LINAC), is a novel technique for prostate cancer (PCa) treatment, promising to further improve clinical outcome and reduce toxicity. The role of prostate MRI has been greatly expanded from the traditional PCa diagnosis to also PCa screening, treatment and surveillance. Diagnostic prostate MRI has been relatively familiar in the community, particularly with the development of Prostate Imaging - Reporting and Data System (PI-RADS). But, on the other hand, the use of MRI in the emerging clinical practice of PCa MRgRT, which is substantially different from that in PCa diagnosis, has been so far sparsely presented in the medical literature. This review attempts to give a comprehensive overview of MRI acquisition techniques currently used in the clinical workflows of PCa MRgRT, from treatment planning to online treatment guidance, in order to promote MRI practice and research for PCa MRgRT. In particular, the major differences in the MRI acquisition of PCa MRgRT from that of diagnostic prostate MRI are demonstrated and explained. Limitations in the current MRI acquisition for PCa MRgRT are analyzed. The future developments of MRI in the PCa MRgRT are also discussed.

Geometrical Comparison and Quantitative Evaluation of 18F-FDG PET/CT- and DW-MRI-Based Target Delineation Before and During Radiotherapy for Esophageal Squamous Carcinoma

Background and Purpose

This study aimed to evaluate the geometrical differences in and metabolic parameters of ¹⁸F-fluorodeoxyglucose positron emission tomography–computed tomography (¹⁸F-FDG PET-CT) and diffusion-weighted magnetic resonance imaging (DW-MRI) performed before and during radiotherapy (RT) for patients with esophageal cancer based on the three-dimensional CT (3DCT) medium and explore whether the high signal area derived from DW-MRI can be used as a tool for an individualized definition of the volume in need of dose escalation for esophageal squamous cancer.

Materials and Methods

Thirty-two patients with esophageal squamous cancer sequentially underwent repeated 3DCT, ¹⁸F-FDG PET-CT, and enhanced MRI before the initiation of RT and after the 15th fraction. All images were fused with 3DCT images through deformable registration. The gross tumor volume (GTV) was delineated based on PET Edge on the first and second PET-CT images and defined as GTV_PETpre and GTV_PETdur, respectively. GTV_DWIpre and GTV_DWIdur were delineated on the first and second DWI and corresponding T₂-weighted MRI (T₂W-MRI)-fused images. The maximum, mean, and peak standardized uptake values (SUVs; SUV_max, SUV_mean, and SUV_peak, respectively); metabolic tumor volume (MTV); and total lesion glycolysis(TLG) and its relative changes were calculated automatically on PET. Similarly, the minimum and mean apparent diffusion coefficient (ADC; ADC_min and ADC_mean) and its relative changes were measured manually using ADC maps.

Results

The volume of GTV_CT exhibited a significant positive correlation with that of GTV_PET and GTV_DWI (both p < 0.001). Significant differences were observed in both ADCs and ¹⁸F-FDG PET metabolic parameters before and during RT (both p < 0.001). No significant correlation was observed between SUVs and ADCs before and during RT (p = 0.072–0.944) and between ∆ADCs and ∆SUVs (p = 0.238–0.854). The conformity index and degree of inclusion of GTV_PETpre to GTV_DWIpre were significantly higher than those of GTV_PETdur to GTV_DWIdur (both p < 0.001). The maximum diameter shrinkage rate (∆LD_DWI) (24%) and the tumor volume shrinkage rate (VRR_DWI) (60%) based on DW-MRI during RT were significantly greater than the corresponding PET-based ∆LD_PET (14%) and VRR_PET (41%) rates (p = 0.017 and 0.000, respectively).

Conclusion

Based on the medium of CT images, there are significant differences in spatial position, biometabolic characteristics, and the tumor shrinkage rate for GTVs derived from ¹⁸F-FDG PET-CT and DW-MRI before and during RT for esophageal squamous cancer. Further studies are needed to determine if DW-MRI will be used as tool for an individualized definition of the volume in need of dose escalation.

Phantom assessment of three-dimensional geometric distortion of a dedicated wide-bore MR-simulator for radiotherapy

This study evaluated the machine-dependent three-dimensional geometric distortion images acquired from a 1.5T 700mm-wide bore MR-simulator based on a large geometric accuracy phantom. With the consideration of radiation therapy (RT) application requirements, every sequence was examined in various combinations of acquisition-orientations and receiver-bandwidths with console-integrated distortion correction enabled. Distortion was repeatedly measured over a six-month period. The distortion measured from the images acquired at the beginning of this period was employed to retrospectively correct the distortion in the subsequent acquisitions. Geometric distortion was analyzed within the largest field-of-view allowed. Six sequences were examined for comprehensive distortion analysis - VIBE, SPACE, TSE, FLASH, BLADE and PETRA. Based on optimal acquisition parameters, their diameter-sphere-volumes (DSVs) of CT-comparable geometric fidelity (where 1mm distortion was allowed) were 333.6mm, 315.1mm, 316.0mm, 318.9mm, 306.2mm and 314.5mm respectively. This was a significant increase from 254.0mm, 245.5mm, 228.9mm, 256.6mm, 230.8mm and 254.2mm DSVs respectively, when images were acquired using un-optimized parameters. The longitudinal stability of geometric distortion and the efficacy of retrospective correction of console-corrected images, based on prior distortion measurements, were inspected using VIBE and SPACE. The retrospectively corrected images achieved over 500mm DSVs with 1mm distortion allowed. The median distortion was below 1mm after retrospective correction, proving that obtaining prior distortion map for subsequent retrospective distortion correction is beneficial. The systematic evaluation of distortion using various combinations of sequence-type, acquisition-orientation and receiver-bandwidth in a six-month time span would be a valuable guideline for optimizing sequence for various RT applications.

Radiomics and Digital Image Texture Analysis in Oncology (Review)

One of the most promising areas of diagnosis and prognosis of diseases is radiomics, a science combining radiology, mathematical modeling, and deep machine learning. The main concept of radiomics is image biomarkers (IBMs), the parameters characterizing various pathological changes and calculated based on the analysis of digital image texture. IBMs are used for quantitative assessment of digital imaging results (CT, MRI, ultrasound, PET). The use of IBMs in the form of “virtual biopsy” is of particular relevance in oncology.

The article provides the basic concepts of radiomics identifying the main stages of obtaining IBMs: data collection and preprocessing, tumor segmentation, data detection and extraction, modeling, statistical processing, and data validation. The authors have analyzed the possibilities of using IBMs in oncology, describing the currently known features and advantages of using radiomics and image texture analysis in the diagnosis and prognosis of cancer. The limitations and problems associated with the use of radiomics data are considered.

Although the novel effective tool for performing virtual biopsy of human tissue is at the development stage, quite a few projects have already been implemented, and medical software packages for radiomics analysis of digital images have been created.

Clinical impact of magnetic resonance imaging distortions on gamma knife radiosurgery

INTRODUCTION

Magnetic resonance imaging (MRI) is the preferred imaging modality for Leksell Gamma Knife® (LGK) stereotactic radiosurgery (SRS) treatment planning (TP) due to superior soft tissue definition compared to computed tomography (CT). However, inherent distortions in MRI can affect treatment accuracy. The aim of this study was to develop a model to visualise the effect of MRI distortion on LGK SRS target coverage.

METHODS

A model was developed using MR images of a QUASARTM GRID3D QA phantom. One hundred and twenty-five points were compared against known phantom geometry. Using linear interpolation, the model was applied retrospectively to 10 brain metastases patient data sets treated with LGK. The model estimated the corrected shot position accounting for distortion. A total of 44 metastases were investigated regarding the effects of MRI distortion on target coverage.

RESULTS

The model indicated significantly reduced mean error by 0.30 mm and variance by 0.09 mm (P = 0.008). After model application, 23 (53%) metastases showed reduced coverage. Six of the 23 metastases were deemed to be potentially clinically significant changes. Results indicated MRI distortion had a greater effect on smaller targets (mean 0.06cc) located further away from the image isocentre (mean 64.88 mm).

CONCLUSION

This study developed a model to visualise the effect of MRI distortion on LGK SRS target coverage. Results suggest that MRI distortion can affect target coverage and the developed model may be one method to assess its impact. These results indicate that MRI distortion may have a greater effect on smaller targets located at the image periphery.

3D T1-weighted turbo spin echo contrast-enhanced MRI at 1.5 T for frameless brain metastases radiotherapy

PURPOSE

Performance of 3D-T1W-TSE has been proven superior to 3D-MP-GRE at 3 T on brain metastases (BM) contrast-enhanced (CE) MRI. However, its performance at 1.5 T is largely unknown and sparsely reported. This study aims to assess image quality, lesion detectability and conspicuity of 1.5 T 3D-T1W-TSE on planning MRI of frameless BM radiotherapy.

METHODS

94 BM patients to be treated by frameless brain radiotherapy were scanned using 3D-T1W-TSE with immobilization on multi-vendor 1.5 T MRI-simulators. BMs were jointly diagnosed by 4 reviewers. Enhanced lesion conspicuity was quantitatively assessed by calculating contrast ratio (CR) and contrast-to-noise ratio (CNR). Signal-to-noise ratio (SNR) reduction of white matter due to the use of flexible coil was assessed. Lesion detectability and conspicuity were compared between 1.5 T planning MRI and 3 T diagnostic MRI by an oncologist and a radiologist in 10 patients.

RESULTS

497 BMs were jointly diagnosed. The CR and CNR were 75.2 ± 39.9% and 14.2 ± 8.1, respectively. SNR reduced considerably from 31.7 ± 8.3 to 21.9 ± 5.4 with the longer distance to coils. 3 T diagnostic MRI and 1.5 T planning MRI yielded exactly the same detection of 84 BMs. Qualitatively, lesion conspicuity at 1.5 T was not inferior to that at 3 T. Quantitatively, lower brain SNR and lesion CNR were found at 1.5 T, while lesion CR at 1.5 T was highly comparable to that at 3 T.

CONCLUSION

1.5 T 3D-T1W-TSE planning MRI of frameless BM radiotherapy was comprehensively assessed. Highly comparable BM detectability and conspicuity were achieved by 1.5 T planning MRI compared to 3 T diagnostic MRI. 1.5 T 3D-T1W-TSE should be valuable for frameless brain radiotherapy planning.

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.

Clinical assessment of geometric distortion for a 0.35T MR-guided radiotherapy system

Purpose

To estimate the overall spatial distortion on clinical patient images for a 0.35 T MR‐guided radiotherapy system.

Methods

Ten patients with head‐and‐neck cancer underwent CT and MR simulations with identical immobilization. The MR images underwent the standard systematic distortion correction post‐processing. The images were rigidly registered and landmark‐based analysis was performed by an anatomical expert. Distortion was quantified using Euclidean distance between each landmark pair and tagged by tissue interface: bone‐tissue, soft tissue, or air‐tissue. For baseline comparisons, an anthropomorphic phantom was imaged and analyzed.

Results

The average spatial discrepancy between CT and MR landmarks was 1.15 ± 1.14 mm for the phantom and 1.46 ± 1.78 mm for patients. The error histogram peaked at 0–1 mm. 66% of the discrepancies were <2 mm and 51% <1 mm. In the patient data, statistically significant differences (p‐values < 0.0001) were found between the different tissue interfaces with averages of 0.88 ± 1.24 mm, 2.01 ± 2.20 mm, and 1.41 ± 1.56 mm for the air/tissue, bone/tissue, and soft tissue, respectively. The distortion generally correlated with the in‐plane radial distance from the image center along the longitudinal axis of the MR.

Conclusion

Spatial distortion remains in the MR images after systematic distortion corrections. Although the average errors were relatively small, large distortions observed at bone/tissue interfaces emphasize the need for quantitative methods for assessing and correcting patient‐specific spatial distortions.

Integration of functional imaging in brachytherapy

Functional imaging allows the evaluation of numerous biological properties that could be considered at all steps of the therapeutic management of patients treated with brachytherapy. Indeed, it enables better initial staging of the disease, and some parameters may also be used as predictive biomarkers for treatment response, allowing better selection of patients eligible for brachytherapy. It may also improve the definition of target volumes with the aim of dose escalations by dose-painting. Finally, it could be useful during the follow-up to assess response to treatment. In this review, we report how functional imaging is integrated at the present time during the brachytherapy procedure, and what are its potential future contributions in the main tumour locations where brachytherapy is recommended. Functional imaging has great potential in the contact of brachytherapy, but still, several issues remain to be resolved before integrating it into clinical practice, especially as a biomarker or in dose painting strategies.

A Comparison of the Distortion in the Same Field MRI and MR-Linac System With a 3D Printed Phantom

PURPOSE

A 3D printed geometric phantom was developed that can be scanned with computed tomography (CT) and magnetic resonance imaging (MRI) to measure the geometric distortion and determine the relevant dose changes.

MATERIALS AND METHODS

A self-designed 3D printed photosensitive resin phantom was used, which adopts grid-like structures and has 822 1 cm2 squares. The scanning plan was delivered by three MRI scanners: the Elekta Unity MR-Linac 1.5T, GE Signa HDe 1.5T, and GE Discovery-sim 750 3.0T. The geometric distortion comparison was concentrated on two 1.5T MRI systems, whereas the 3.0T MRI was used as a supplemental experiment. The most central transverse images in each dataset were selected to demonstrate the plane distortion. Some mark points were selected to analyze the distortion in the 3D direction based on the plane geometric distortion. A treatment plan was created with the off-line Monaco system.

RESULTS

The distortion increases gradually from the center to the outside. The distortion range is 0.79 ± 0.40 mm for the Unity, 1.31 ± 0.56 mm for the GE Signa HDe, and 2.82 ± 1.48 mm for the GE Discovery-sim 750. Additionally, the geometric distortion slightly affects the actual planning dose of the radiotherapy.

CONCLUSION

Geometric distortion increases gradually from the center to the outside. The distortion values of the Unity were smaller than those of the GE Signa HDe, and the Unity has the smallest geometric distortion. Finally, the Unity's dose variation best matched with the standard treatment plan.

In phantom evaluation of targeting accuracy in MRI-based brain radiosurgery

PURPOSE

To determine the targeting accuracy of brain radiosurgery when planning procedures employing different MRI and MRI + CT combinations are adopted.

MATERIALS AND METHOD

A new phantom, the BrainTool, has been designed and realized to test image co-registration and targeting accuracy in a realistic anatomical situation. The phantom was created with a 3D printer and materials that mimic realistic brain MRI and CT contrast using a model extracted from a synthetic MRI study of a human brain. Eight markers distributed within the BrainTool provide for assessment of the accuracy of image registrations while two cavities that host an ionization chamber are used to perform targeting accuracy measurements with an iterative cross-scan method. Two procedures employing 1.5 T MRI-only or a combination of MRI (taken with 1.5 T or 3 T scanners) and CT to carry out Gamma Knife treatments were investigated. As distortions can impact targeting accuracy, MR images were preliminary evaluated to assess image deformation extent using GammaTool phantom.

RESULTS

MR images taken with both scanners showed average and maximum distortion of 0.3 mm and 1 mm respectively. The marker distances in co-registered images resulted below 0.5 mm for both MRI scans. The targeting mismatches obtained were 0.8 mm, 1.0 mm and 1.2 mm for MRI-only and MRI + CT (1,5T and 3 T), respectively.

CONCLUSIONS

Procedures using a combination of MR and CT images provide targeting accuracies comparable to those of MRI-only procedures. The BrainTool proved to be a suitable tool for carrying out co-registration and targeting accuracy of Gamma Knife brain radiosurgery treatments.

Determination of acceptance criteria for geometric accuracy of magnetic resonance imaging scanners used in radiotherapy planning

•

Geometric accuracy of MRI-scanners in radiotherapy planning must be evaluated.

•

Phantom acquisitions with standard and clinical sequences were performed.

•

Geometric distortions were determined in several volumes of interest.

•

We recommend acceptance criteria for MRI-scanners in radiotherapy planning.

•

Explicit and simple acceptance criteria enable effective regulatory inspections.

Background and Purpose

Magnetic resonance imaging is increasingly used in radiotherapy planning; yet, the performance of the utilized scanners is rarely regulated by any authority. The aim of this study was to determine the geometric accuracy of several magnetic resonance imaging scanners used for radiotherapy planning, and to establish acceptance criteria for such scanners.

Materials and Methods

The geometric accuracy of five different scanners was measured with three sequences using a commercial large-field-of-view phantom. The distortion magnitudes were determined in spherical volumes around the scanner isocenter and in cylindrical volumes along scanner z-axis. The repeatability of the measurements was determined on a single scanner with two quality assurance sequences with three single-setup and seven repeated-setup measurements.

Results

For all scanners and sequences except one, the mean and median distortion magnitude was <1 mm and <2 mm in spherical volumes with diameters of 400 mm and 500 mm, respectively. For all sequences maximum distortion was <2 mm in spherical volume with diameter of 300 mm. The mean standard deviation of marker-by-marker distortion magnitudes over repeated acquisitions was ≤0.6 mm with both tested sequences.

Conclusions

All tested scanners were geometrically accurate for their current use in radiotherapy planning. The acceptance criteria of geometric accuracy for regulatory inspections of a supervising authority could be set according to these results.

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.

A 3D printed phantom to assess MRI geometric distortion

Magnetic Resonance has become a standard imaging modality for target volume delineation and treatment planning in radiation oncology. Geometric distortions, however, have the potential to detrimentally affect both tumour definition and the dose delivered to the target volume. We report the design, fabrication and imaging of a 3D printed unibody MR distortion phantom along with quantitative image analysis.

METHODS

The internal cavity of the phantom is an orthogonal three-dimensional planar lattice, composed of 3mm diameter rods spaced equidistantly at a 20mm centre-centre offset repeating along the X, Y and Z axes. The phantom featured an overall length of 308.5 mm, a width of 246 mm and a height of 264 mm with lines on the external surface for phantom positioning matched to external lasers. The MR phantom was 3D printed in Nylon-12 using an advancement on traditional selective laser sintering (SLS) (HP Jet Fusion 3D - 4200 machine). The phantom was scanned on a Toshiba Aquilion CT scanner to check the integrity of the 3D print and to correct for any resultant issues. The phantom was then filled with NiSO4 solution and scanned on a 3T PET-MR Siemens scanner for selected T1 and T2 sequences, from which distortion vectors were generated and analysed using in-house software written in Python.

RESULTS

All deviations were less than 1 mm, with an average displacement of 0.228 mm. The majority of the deviations are smaller than the 0.692 mm pixel size for this dataset.

CONCLUSION

A cost-effective, 3D printed MRI-phantom was successfully printed and tested for assessing geometric distortion on MRI scanners. The custom phantom with markings for phantom alignment may be considered for radiotherapy departments looking to add MR scanners for simulation and image guidance.

Trimodality PET/CT/MRI and Radiotherapy: A Mini-Review

Computed tomography (CT) has revolutionized external radiotherapy by making it possible to visualize and segment the tumors and the organs at risk in a three-dimensional way. However, if CT is a now a standard, it presents some limitations, notably concerning tumor characterization and delineation. Its association with functional and anatomical images, that are positron emission tomography (PET) and magnetic resonance imaging (MRI), surpasses its limits. This association can be in the form of a trimodality PET/CT/MRI. The objective of this mini-review is to describe the process of performing this PET/CT/MRI trimodality for radiotherapy and its potential clinical applications. Trimodality can be performed in two ways, either a PET/MRI fused to a planning CT (possibly with a pseudo-CT generated from the MRI for the planning), or a PET/CT fused to an MRI and then registered to a planning CT (possibly the CT of PET/CT if calibrated for radiotherapy). These examinations should be performed in the treatment position, and in the second case, a patient transfer system can be used between the PET/CT and MRI to limit movement. If trimodality requires adapted equipment, notably compatible MRI equipment with high-performance dedicated coils, it allows the advantages of the three techniques to be combined with a synergistic effect while limiting their disadvantages when carried out separately. Trimodality is already possible in clinical routine and can have a high clinical impact and good inter-observer agreement, notably for head and neck cancers, brain tumor, prostate cancer, cervical cancer.

Evaluation of the influence of susceptibility-induced magnetic field distortions on the precision of contouring intracranial organs at risk for stereotactic radiosurgery

•

45 data sets (18 on a 1.5 T MR and 27 on a 3 T MR) were evaluated for susceptibility induced distortions.

•

Maximum distortions of up to 1.7 mm were found for organs at risk in standard diagnostic settings.

•

Median distortions ranged between 0.1 and 0.2 mm for all organs at risk.

•

Active shimming was estimated to reduce distortions by a factor of 2.3 to 2.9.

•

A safety margin of 1 mm would have encompassed 99.8% of the distortions.

Background and purpose

Magnetic resonance imaging (MRI) is a crucial factor in optimal treatment planning for stereotactic radiosurgery. To further the awareness of possible errors in MRI, this work aimed to investigate the magnitude of susceptibility induced MRI distortions for intracranial organs at risk (OARs) and test the effectiveness of actively shimming these distortions.

Materials and methods

Distortion maps for 45 exams of 42 patients (18 on a 1.5 T MRI scanner, 27 on a 3 T MRI scanner) were calculated based on a high-bandwidth double-echo gradient echo sequence. The investigated OARs were brainstem, chiasm, eyes, and optic nerves. The influence of active shimming was investigated by comparing unshimmed 1.5 T data with shimmed 3 T data and comparing the results to a model based prediction.

Results

The median distortion for the different OARs was found to be between 0.13 and 0.18 mm for 1.5 T and between 0.11 and 0.13 mm for 3 T. The maximum distortion was found to be between 1.3 and 1.7 mm for 1.5 T and between 1.1 and 1.4 mm for 3 T. The variation of values was much higher for 1.5 T than for 3 T across all investigated OARs. Active shimming was found to reduce distortions by a factor of 2.3 to 2.9 compared to the expected values.

Conclusions

Using a safety margin for OARs of 1 mm would have encompassed 99.8% of the distortions. Since distortions are inversely proportional to the readout bandwidth, they can be further reduced by increasing the bandwidth. Additional error sources like gradient nonlinearities need to be addressed separately.