Respiratory motion-resolved, self-gated 4D-MRI using Rotating Cartesian K-space (ROCK): Initial clinical experience on an MRI-guided radiotherapy system

Intra-frame motion deterioration effects and deep-learning-based compensation in MR-guided radiotherapy

BACKGROUND

Current commercially available hybrid magnetic resonance linear accelerators (MR-Linac) use 2D+t cine MR imaging to provide intra-fractional motion monitoring. However, given the limited temporal resolution of cine MR imaging, target intra-frame motion deterioration effects, resulting in effective time latency and motion artifacts in the image domain, can be appreciable, especially in the case of fast breathing.

PURPOSE

The aim of this work is to investigate intra-frame motion deterioration effects in MR-guided radiotherapy (MRgRT) by simulating the motion-corrupted image acquisition, and to explore the feasibility of deep-learning-based compensation approaches, relying on the intra-frame motion information which is spatially and temporally encoded in the raw data (k-space).

METHODS

An intra-frame motion model was defined to simulate motion-corrupted MR images, with 4D anthropomorphic digital phantoms being exploited to provide ground truth 2D+t cine MR sequences. A total number of 10 digital phantoms were generated for lung cancer patients, with randomly selected eight patients for training or validation and the remaining two for testing. The simulation code served as the data generator, and a dedicated motion pattern perturbation scheme was proposed to build the intra-frame motion database, where three degrees of freedom were designed to guarantee the diversity of intra-frame motion trajectories, enabling a thorough exploration in the domain of the potential anatomical structure positions. U-Nets with three types of loss functions: L1 or L2 loss defined in image or Fourier domain, referred to as NNImgLoss-L1 , NNFloss-L1 and NNL2-Loss were trained to extract information from the motion-corrupted image and used to estimate the ground truth final-position image, corresponding to the end of the acquisition. Images before and after compensation were evaluated in terms of (i) image mean-squared error (MSE) and mean absolute error (MAE), and (ii) accuracy of gross tumor volume (GTV) contouring, based on optical-flow image registration.

RESULTS

Image degradation caused by intra-frame motion was observed: for a linearly and fully acquired Cartesian readout k-space trajectory, intra-frame motion resulted in an imaging latency of approximately 50% of the acquisition time; in comparison, the motion artifacts exhibited only a negligible contribution to the overall geometric errors. All three compensation models led to a decrease in image MSE/MAE and GTV position offset compared to the motion-corrupted image. In the investigated testing dataset for GTV contouring, the average dice similarity coefficients (DSC) improved from 88% to 96%, and the 95th percentile Hausdorff distance (HD95 ) dropped from 4.8 mm to 2.1 mm. Different models showed slight performance variations across different intra-frame motion amplitude categories: NNImgLoss-L1 excelled for small/medium amplitudes, whereas NNFloss-L1 demonstrated higher DSC median values at larger amplitudes. The saliency maps of the motion-corrupted image highlighted the major contribution of the later acquired k-space data, as well as the edges of the moving anatomical structures at their final positions, during the model inference stage.

CONCLUSIONS

Our results demonstrate the deep-learning-based approaches have the potential to compensate for intra-frame motion by utilizing the later acquired data to drive the convergence of the earlier acquired k-space components.

The Use of Magnetic Resonance Imaging in Radiation Therapy Treatment Simulation and Planning

Ever since its introduction as a diagnostic imaging tool the potential of magnetic resonance imaging (MRI) in radiation therapy (RT) treatment simulation and planning has been recognized. Recent technical advances have addressed many of the impediments to use of this technology and as a result have resulted in rapid and growing adoption of MRI in RT. The purpose of this article is to provide a broad review of the multiple uses of MR in the RT treatment simulation and planning process, identify several of the most used clinical scenarios in which MR is integral to the simulation and planning process, highlight existing limitations and provide multiple unmet needs thereby highlighting opportunities for the diagnostic MR imaging community to contribute and collaborate with our oncology colleagues. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 5.

Intrafraction Motion Management With MR-Guided Radiation Therapy

High quality radiation therapy requires highly accurate and precise dose delivery. MR-guided radiotherapy (MRgRT), integrating an MRI scanner with a linear accelerator, offers excellent quality images in the treatment room without subjecting patient to ionizing radiation. MRgRT therefore provides a powerful tool for intrafraction motion management. This paper summarizes different sources of intrafraction motion for different disease sites and describes the MR imaging techniques available to visualize and quantify intrafraction motion. It provides an overview of MR guided motion management strategies and of the current technical capabilities of the commercially available MRgRT systems. It describes how these motion management capabilities are currently being used in clinical studies, protocols and provides a future outlook.

Respiratory motion management using a single rapid MRI scan for a 0.35 T MRI-Linac system

BACKGROUND

MRI has a rapidly growing role in radiation therapy (RT) for treatment planning, real-time image guidance, and beam gating (e.g., MRI-Linac). Free-breathing 4D-MRI is desirable in respiratory motion management for therapy. Moreover, high-quality 3D-MRIs without motion artifacts are needed to delineate lesions. Existing MRI methods require multiple scans with lengthy acquisition times or are limited by low spatial resolution, contrast, and signal-to-noise ratio.

PURPOSE

We developed a novel method to obtain motion-resolved 4D-MRIs and motion-integrated 3D-MRI reconstruction using a single rapid (35-45 s scan on a 0.35 T MRI-Linac.

METHODS

Golden-angle radial stack-of-stars MRI scans were acquired from a respiratory motion phantom and 12 healthy volunteers (n = 12) on a 0.35 T MRI-Linac. A self-navigated method was employed to detect respiratory motion using 2000 (acquisition time = 5-7 min) and the first 200 spokes (acquisition time = 35-45 s). Multi-coil non-uniform fast Fourier transform (MCNUFFT), compressed sensing (CS), and deep-learning Phase2Phase (P2P) methods were employed to reconstruct motion-resolved 4D-MRI using 2000 spokes (MCNUFFT2000) and 200 spokes (CS200 and P2P200). Deformable motion vector fields (MVFs) were computed from the 4D-MRIs and used to reconstruct motion-corrected 3D-MRIs with the MOtion Transformation Integrated forward-Fourier (MOTIF) method. Image quality was evaluated quantitatively using the structural similarity index measure (SSIM) and the root mean square error (RMSE), and qualitatively in a blinded radiological review.

RESULTS

Evaluation using the respiratory motion phantom experiment showed that the proposed method reversed the effects of motion blurring and restored edge sharpness. In the human study, P2P200 had smaller inaccuracy in MVFs estimation than CS200. P2P200 had significantly greater SSIMs (p < 0.0001) and smaller RMSEs (p < 0.001) than CS200 in motion-resolved 4D-MRI and motion-corrected 3D-MRI. The radiological review found that MOTIF 3D-MRIs using MCNUFFT2000 exhibited the highest image quality (scoring > 8 out of 10), followed by P2P200 (scoring > 5 out of 10), and then motion-uncorrected (scoring < 3 out of 10) in sharpness, contrast, and artifact-freeness.

CONCLUSIONS

We have successfully demonstrated a method for respiratory motion management for MRI-guided RT. The method integrated self-navigated respiratory motion detection, deep-learning P2P 4D-MRI reconstruction, and a motion integrated reconstruction (MOTIF) for 3D-MRI using a single rapid MRI scan (35-45 s) on a 0.35 T MRI-Linac system.

MRI-guided Radiotherapy (MRgRT) for treatment of Oligometastases: Review of clinical applications and challenges

PURPOSE

Early clinical results on the application of magnetic resonance imaging (MRI) coupled with a linear accelerator to deliver MR-guided radiation therapy (MRgRT) have demonstrated feasibility for safe delivery of stereotactic body radiotherapy (SBRT) in treatment of oligometastatic disease. Here we set out to review the clinical evidence and challenges associated with MRgRT in this setting.

METHODS AND MATERIALS

We performed a systematic review of the literature pertaining to clinical experiences and trials on the use of MRgRT primarily for the treatment of oligometastatic cancers. We reviewed the opportunities and challenges associated with the use of MRgRT.

RESULTS

Benefits of MRgRT pertaining to superior soft-tissue contrast, real-time imaging and gating, and online adaptive radiotherapy facilitate safe and effective dose escalation to oligometastatic tumors while simultaneously sparing surrounding healthy tissues. Challenges concerning further need for clinical evidence and technical considerations related to planning, delivery, quality assurance (QA) of hypofractionated doses, and safety in the MRI environment must be considered.

CONCLUSIONS

The promising early indications of safety and effectiveness of MRgRT for SBRT-based treatment of oligometastatic disease in multiple treatment locations should lead to further clinical evidence to demonstrate the benefit of this technology.

Dosimetry impact of gating latency in cine magnetic resonance image guided breath-hold pancreatic cancer radiotherapy

Objective.We investigated dosimetry effect of gating latency in cine magnetic resonance image (cine MRI) guided breath-hold pancreatic cancer radiotherapy.Approach.The gating latency was calculated based on cine MRI obtained from 17 patients who received MRI guided radiotherapy. Because of the cine MRI-related latency, beam overshoot occurs when beam remains on while the tracking target already moves out of the target boundary. The number of beam on/off events was calculated from the cine MRI data. We generated both IMRT and VMAT plans for all 17 patients using 33 Gy prescription, and created motion plans by applying isocenter shift that corresponds to motion-induced tumor displacement. The GTV and PTV coverage and dose to nearby critical structures were compared between the motion and original plan to evaluate the dosimetry change caused by cine MRI latency.Main results.The time ratio of cine MRI imaging latency over the treatment duration is 6.6 ± 3.1%, the mean and median percentage of beam-on events <4 s are 67.0 ± 14.3% and 66.6%. When a gating boundary of 4 mm and a target-out threshold of 5% is used, there is no significant difference for GTV V33Gy between the motion and original plan (p = 0.861 and 0.397 for IMRT and VMAT planning techniques, respectively). However, the PTV V33Gy and stomach Dmax for the motion plans are significantly lower; duodenum V12.5 Gy and V18Gy are significantly higher when compared with the original plans, for both IMRT and VMAT planning techniques.Significance.The cine MRI gating latency can significantly decrease the dose delivered to the PTV, and increase the dose to the nearby critical structures. However, no significant difference is observed for the GTV coverage. The dosimetry impact can be mitigated by implementing additional beam-on control techniques which reduces unnecessary beam on events and/or by using faster cine MRI sequences which reduces the latency period.

Abdominal motion tracking with free-breathing XD-GRASP acquisitions using spatio-temporal geodesic trajectories

Free-breathing external beam radiotherapy remains challenging due to the complex elastic or irregular motion of abdominal organs, as imaging moving organs leads to the creation of motion blurring artifacts. In this paper, we propose a radial-based MRI reconstruction method from 3D free-breathing abdominal data using spatio-temporal geodesic trajectories, to quantify motion during radiotherapy. The prospective study was approved by the institutional review board and consent was obtained from all participants. A total of 25 healthy volunteers, 12 women and 13 men (38 years ± 12 [standard deviation]), and 11 liver cancer patients underwent imaging using a 3.0 T clinical MRI system. The radial acquisition based on golden-angle sparse sampling was performed using a 3D stack-of-stars gradient-echo sequence and reconstructed using a discretized piecewise spatio-temporal trajectory defined in a low-dimensional embedding, which tracks the inhale and exhale phases, allowing the separation between distinct motion phases. Liver displacement between phases as measured with the proposed radial approach based on the deformation vector fields was compared to a navigator-based approach. Images reconstructed with the proposed technique with 20 motion states and registered with the multiscale B-spline approach received on average the highest Likert scores for the overall image quality and visual SNR score 3.2 ± 0.3 (mean ± standard deviation), with liver displacement errors varying between 0.1 and 2.0 mm (mean 0.8 ± 0.6 mm). When compared to navigator-based approaches, the proposed method yields similar deformation vector field magnitudes and angle distributions, and with improved reconstruction accuracy based on mean squared errors. Schematic illustration of the proposed 4D-MRI reconstruction method based on radial golden-angle acquisitions and a respiration motion model from a manifold embedding used for motion tracking. First, data is extracted from the center of k-space using golden-angle sampling, which is then mapped onto a low-dimensional embedding, describing the relationship between neighboring samples in the breathing cycle. The trained model is then used to extract the respiratory motion signal for slice re-ordering. The process then improves the image quality through deformable image registration. Using a reference volume, the deformation vector field (DVF) of sequential motion states are extracted, followed by deformable registrations. The output is a 4DMRI which allows to visualize and quantify motion during free-breathing.

Technical Challenges of Real-Time Adaptive MR-Guided Radiotherapy

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

Practical Safety Considerations for Integration of Magnetic Resonance Imaging in Radiation Therapy

Interest in integrating magnetic resonance imaging (MRI) in radiation therapy (RT) practice has increased dramatically in recent years owing to its unique advantages such as excellent soft tissue contrast and capability of measuring biological properties. Continuous real-time imaging for intrafractional motion tracking without ionizing radiation serves as a particularly attractive feature for applications in RT. Despite its many advantages, the integration of MRI in RT workflows is not straightforward, with many unmet needs. MR safety remains one of the key challenges and concerns in the clinical implementation of MR simulators and MR-guided radiation therapy systems in radiation oncology. Most RT staff are not accustomed to working in an environment with a strong magnetic field. There are specific requirements in RT that are different from diagnostic applications. A large variety of implants and devices used in routine RT practice do not have clear MR safety labels. RT-specific imaging pulse sequences focusing on fast acquisition, high spatial integrity, and continuous, real-time acquisition require additional MR safety testing and evaluation. This article provides an overview of MR safety tailored toward RT staff, followed by discussions on specific requirements and challenges associated with MR safety in the RT environment. Strategies and techniques for developing an MR safety program specific to RT are presented and discussed.

Artificial intelligence in radiotherapy: a technological review

Radiation therapy (RT) is widely used to treat cancer. Technological advances in RT have occurred in the past 30 years. These advances, such as three-dimensional image guidance, intensity modulation, and robotics, created challenges and opportunities for the next breakthrough, in which artificial intelligence (AI) will possibly play important roles. AI will replace certain repetitive and labor-intensive tasks and improve the accuracy and consistency of others, particularly those with increased complexity because of technological advances. The improvement in efficiency and consistency is important to manage the increasing cancer patient burden to the society. Furthermore, AI may provide new functionalities that facilitate satisfactory RT. The functionalities include superior images for real-time intervention and adaptive and personalized RT. AI may effectively synthesize and analyze big data for such purposes. This review describes the RT workflow and identifies areas, including imaging, treatment planning, quality assurance, and outcome prediction, that benefit from AI. This review primarily focuses on deep-learning techniques, although conventional machine-learning techniques are also mentioned.

4D MRI: Robust sorting of free breathing MRI slices for use in interventional settings

Purpose

We aim to develop a robust 4D MRI method for large FOVs enabling the extraction of irregular respiratory motion that is readily usable with all MRI machines and thus applicable to support a wide range of interventional settings.

Method

We propose a 4D MRI reconstruction method to capture an arbitrary number of breathing states. It uses template updates in navigator slices and search regions for fast and robust vessel cross-section tracking. It captures FOVs of 255 mm x 320 mm x 228 mm at a spatial resolution of 1.82 mm x 1.82 mm x 4mm and temporal resolution of 200ms. A total of 37 4D MRIs of 13 healthy subjects were reconstructed to validate the method. A quantitative evaluation of the reconstruction rate and speed of both the new and baseline method was performed. Additionally, a study with ten radiologists was conducted to assess the subjective reconstruction quality of both methods.

Results

Our results indicate improved mean reconstruction rates compared to the baseline method (79.4% vs. 45.5%) and improved mean reconstruction times (24s vs. 73s) per subject. Interventional radiologists perceive the reconstruction quality of our method as higher compared to the baseline (262.5 points vs. 217.5 points, p = 0.02).

Conclusions

Template updates are an effective and efficient way to increase 4D MRI reconstruction rates and to achieve better reconstruction quality. Search regions reduce reconstruction time. These improvements increase the applicability of 4D MRI as a base for seamless support of interventional image guidance in percutaneous interventions.

MRSIGMA: Magnetic Resonance SIGnature MAtching for real-time volumetric imaging

PURPOSE

To propose a real-time 3D MRI technique called MR SIGnature MAtching (MRSIGMA) for high-resolution volumetric imaging and motion tracking with very low imaging latency.

METHODS

MRSIGMA consists of two steps: (1) offline learning of a database of possible 3D motion states and corresponding motion signature ranges and (2) online matching of new motion signatures acquired in real time with prelearned motion states. Specifically, the offline learning step (non-real-time) reconstructs motion-resolved 4D images representing different motion states and assigns a unique motion range to each state. The online matching step (real-time) acquires motion signatures only and selects one of the prelearned 3D motion states for each newly acquired signature, which generates 3D images efficiently in real time. The MRSIGMA technique was evaluated on 15 golden-angle stack-of-stars liver data sets, and the performance of respiratory motion tracking with the online-generated real-time 3D MRI was compared with the corresponding 2D projections acquired in real time.

RESULTS

The total latency of generating each 3D image during online matching was about 300 ms, including acquisition of the motion signature data (~138 ms) and corresponding matching process (~150 ms). Linear correlation assessment suggested excellent correlation (R²  = 0.948) between motion displacement measured from the online-generated real-time 3D images and the 2D real-time projections.

CONCLUSION

This proof-of-concept study demonstrates the feasibility of MRSIGMA for high-resolution real-time volumetric imaging, which shifts the acquisition and reconstruction burden to an offline learning step and leaves fast online matching for online imaging with very low imaging latency. The MRSIGMA technique can potentially be used for real-time motion tracking in MRI-guided radiation therapy.

Virtual 4DCT from 4DMRI for the management of respiratory motion in carbon ion therapy of abdominal tumors

PURPOSE

To evaluate a method for generating virtual four-dimensional computed tomography (4DCT) from four-dimensional magnetic resonance imaging (4DMRI) data in carbon ion radiotherapy with pencil beam scanning for abdominal tumors.

METHODS

Deformable image registration is used to: (a) register each respiratory phase of the 4DMRI to the end-exhale MRI; (b) register the reference end-exhale CT to the end-exhale MRI volume; (c) generate the virtual 4DCT by warping the registered CT according to the obtained deformation fields. A respiratory-gated carbon ion treatment plan is optimized on the planning 4DCT and the corresponding dose distribution is recalculated on the virtual 4DCT. The method was validated on a digital anthropomorphic phantom and tested on eight patients (18 acquisitions). For the phantom, a ground truth dataset was available to assess the method performances from the geometrical and dosimetric standpoints. For the patients, the virtual 4DCT was compared with the planning 4DCT.

RESULTS

In the phantom, the method exhibits a geometrical accuracy within the voxel size and Dose Volume Histograms deviations up to 3.3% for target V_95% (mean dose difference ≤ 0.2% of the prescription dose, gamma pass rate > 98%). For patients, the virtual and the planning 4DCTs show good agreement at end-exhale (3% median D_95% difference), whereas other respiratory phases exhibit moderate motion variability with consequent dose discrepancies, confirming the need for motion mitigation strategies during treatment.

CONCLUSIONS

The virtual 4DCT approach is feasible to evaluate treatment plan robustness against intra- and interfraction motion in carbon ion therapy delivered at the abdominal site.

Respiratory-Correlated (RC) vs. Time-Resolved (TR) Four-Dimensional Magnetic Resonance Imaging (4DMRI) for Radiotherapy of Thoracic and Abdominal Cancer

Recent technological and clinical advancements of both respiratory-correlated (RC) and time-resolved (TR) four-dimensional magnetic resonance imaging (4DMRI) techniques are reviewed in light of tumor/organ motion simulation, monitoring, and assessment in radiotherapy. For radiotherapy of thoracic and abdominal cancer, respiratory-induced tumor motion, and motion variation due to breathing irregularities are the major uncertainties in treatment. RC-4DMRI is developed to assess tumor motion for treatment planning, whereas TR-4DMRI is developed to assess both motion and motion variation for treatment planning, delivery and assessment. RC-4DMRI is reconstructed to provide one-breathing-cycle motion, similar to 4D computed tomography (4DCT), the current clinical standard, but with higher soft-tissue contrast, no ionizing radiation, and less binning artifacts due to the use of an internal respiratory surrogate. Recent studies have shown that its spatial resolution has reached or exceeded that of 4DCT and scanning time becomes clinically acceptable. TR-4DMRI is recently developed with an adequate spatiotemporal resolution to assess tumor motion and motion variations for treatment simulation, delivery and assessment. The super-resolution approach is most promising since it can image any organ/body motion, whereas RC-4D MRI are limited to resolve only respiration-induced motion and some TR-4DMRI approaches may more or less depend on RC-4DMRI. TR-4DMRI provides multi-breath motion data that are useful not only in MR-guided radiotherapy but also for building a patient-specific motion model to guide radiotherapy treatment using an non-MR-equipped linear accelerator. Based on 4DMRI motion data, motion-corrected dynamic contrast imaging and diffusion-weighted imaging have also been reported, aiming to facilitate tumor delineation for more accurate radiotherapy targeting. Both RC- and TR-4DMRI have been evaluated for potential clinical applications, such as delineation of tumor volumes, where sufficiently high spatial resolution and large field-of-view are required. The 4DMRI techniques are promising to play a role in motion assessment in radiotherapy treatment planning, delivery, assessment, and adaptation.

"Après Mois, Le Déluge": Preparing for the Coming Data Flood in the MRI-Guided Radiotherapy Era

Magnetic resonance imaging provides a sea of quantitative and semi-quantitative data. While radiation oncologists already navigate a pool of clinical (semantic) and imaging data, the tide will swell with the advent of hybrid MRI/linear accelerator devices and increasing interest in MRI-guided radiotherapy (MRIgRT), including adaptive MRIgRT. The variety of MR sequences (of greater complexity than the single parameter Hounsfield unit of CT scanning routinely used in radiotherapy), the workflow of adaptive fractionation, and the sheer quantity of daily images acquired are challenges for scaling this technology. Biomedical informatics, which is the science of information in biomedicine, can provide helpful insights for this looming transition. Funneling MRIgRT data into clinically meaningful information streams requires committing to the flow of inter-institutional data accessibility and interoperability initiatives, standardizing MRIgRT dosimetry methods, streamlining MR linear accelerator workflow, and standardizing MRI acquisition and post-processing. This review will attempt to conceptually ford these topics using clinical informatics approaches as a theoretical bridge.

A fast volumetric 4D-MRI with sub-second frame rate for abdominal motion monitoring and characterization in MRI-guided radiotherapy

Background

To propose a fast volumetric 4D-MRI based on 3D pulse sequence acquisition for abdominal motion monitoring and characterization in MRI-guided radiotherapy (MRgRT).

Methods

A 3D spoiled gradient echo sequence volumetric interpolated breath-hold examination (VIBE) [repetition time/echo time (TR/TE) =0.53/1.57 ms, flip-angle =5°, receiver bandwidth (RBW) =1,400 Hz/voxel] based 4D-MRI acquisition, accelerated by 4-fold controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA), named CAIPIRINHA-VIBE 4D-MRI, was implemented on a 1.5T MRI simulator (MR-sim) and applied for abdominal imaging of nine healthy volunteers under free breathing. One hundred and forty-four dynamics of the entire abdomen volume (56 slices), in total 8,064 (144×56) images with a voxel size of 2.7×2.7×4.0 mm³, were acquired in 89 s for 4D-MRI. This CAIPIRINHA-VIBE 4D-MRI was qualitatively compared with a 2D half-Fourier acquisition single-shot turbo spin-echo (2D-HASTE) based 4D-MRI. The motions of liver dome, kidney and spleen were analyzed using the CAIPIRINHA-VIBE 4D-MRI data. The kidney motion was quantitatively characterized in terms of motion range and the correlations between left and right kidneys.

Results

CAIPIRINHA-VIBE 4D-MRI was successfully conducted in all subjects. CAIPIRINHA-VIBE 4D-MRI exhibited much higher effective volumetric temporal resolution (0.615 vs. ~5 s/volume) and better reconstructed volume consistency than 2D-HASTE 4D-MRI. CAIPIRINHA-VIBE 4D-MRI was able to characterize the respiratory motion of abdominal organs simultaneously in three orthogonal directions, and could potentially be used for whole abdomen deformable motion tracking. Renal motion range was most pronounced in superior-inferior (SI) direction (L: 10.03±2.65 mm; R: 10.38±2.80 mm), significantly larger (P<0.001) than that in anterior-posterior (AP) and the least in left-right (LR) directions. Right kidney had significantly larger mobility (4.18±2.19 vs. 2.32±1.34 mm, P=0.045) than left kidney in AP, but not in LR and SI directions. The Pearson correlation coefficients r between left and right kidney motion were 0.5063 (P=0.164), 0.6624 (P=0.052) and 0.5752 (P=0.105) in LR, AP and SI correspondingly. The correlation of renal motion in SI and AP was found significant in right kidney (r=0.843, P=0.004) but not in left kidney (r=0.467, P=0.205).

Conclusions

A fast volumetric 4D-MRI was implemented for abdominal motion monitoring in MRgRT. A sub-second volumetric temporal resolution of 0.615 s, covering the entire abdomen, was demonstrated for respiratory motion monitoring and characterization. This technique holds potentials for MRgRT applications.