Novel Super-Resolution Approach to Time-Resolved Volumetric 4-Dimensional Magnetic Resonance Imaging With High Spatiotemporal Resolution for Multi-Breathing Cycle Motion Assessment

Ultra-fast multi-parametric 4D-MRI image reconstruction for real-time applications using a downsampling-invariant deformable registration (D2R) model

BACKGROUND AND PURPOSE

Motion estimation from severely downsampled 4D-MRI is essential for real-time imaging and tumor tracking. This simulation study developed a novel deep learning model for simultaneous MR image reconstruction and motion estimation, named the Downsampling-Invariant Deformable Registration (D2R) model.

MATERIALS AND METHODS

Forty-three patients undergoing radiotherapy for liver tumors were recruited for model training and internal validation. Five prospective patients from another center were recruited for external validation. Patients received 4D-MRI scans and 3D MRI scans. The 4D-MRI was retrospectively down-sampled to simulate real-time acquisition. Motion estimation was performed using the proposed D2R model. The accuracy and robustness of the proposed D2R model and baseline methods, including Demons, Elastix, the parametric total variation (pTV) algorithm, and VoxelMorph, were compared. High-quality (HQ) 4D-MR images were also constructed using the D2R model for real-time imaging feasibility verification. The image quality and motion accuracy of the constructed HQ 4D-MRI were evaluated.

RESULTS

The D2R model showed significantly superior and robust registration performance than all the baseline methods at downsampling factors up to 500. HQ T1-weighted and T2-weighted 4D-MR images were also successfully constructed with significantly improved image quality, sub-voxel level motion error, and real-time efficiency. External validation demonstrated the robustness and generalizability of the technique.

CONCLUSION

In this study, we developed a novel D2R model for deformation estimation of downsampled 4D-MR images. HQ 4D-MR images were successfully constructed using the D2R model. This model may expand the clinical implementation of 4D-MRI for real-time motion management during liver cancer treatment.

Feasibility of delivered dose reconstruction for MR-guided SBRT of pancreatic tumors with fast, real-time 3D cine MRI

BACKGROUND AND PURPOSE

In MR-guided SBRT of pancreatic cancer, intrafraction motion is typically monitored with (interleaved) 2D cine MRI. However, tumor surroundings are often not fully captured in these images, and motion might be distorted by through-plane movement. In this study, the feasibility of highly accelerated 3D cine MRI to reconstruct the delivered dose during MR-guided SBRT was assessed.

MATERIALS AND METHODS

A 3D cine MRI sequence was developed for fast, time-resolved 4D imaging, featuring a low spatial resolution that allows for rapid volumetric imaging at 430 ms. The 3D cines were acquired during the entire beam-on time of 23 fractions of online adaptive MR-guided SBRT for pancreatic tumors on a 1.5 T MR-Linac. A 3D deformation vector field (DVF) was extracted for every cine dynamic using deformable image registration. Next, these DVFs were used to warp the partial dose delivered in the time interval between consecutive cine acquisitions. The warped dose plans were summed to obtain a total delivered dose. The delivered dose was also calculated under various motion correction strategies. Key DVH parameters of the GTV, duodenum, small bowel and stomach were extracted from the delivered dose and compared to the planned dose. The uncertainty of the calculated DVFs was determined with the inverse consistency error (ICE) in the high-dose regions.

RESULTS

The mean (SD) relative (ratio delivered/planned) D_99% of the GTV was 0.94 (0.06), and the mean (SD) relative D_0.5cc of the duodenum, small bowel, and stomach were respectively 0.98 (0.04), 1.00 (0.07), and 0.98 (0.06). In the fractions with the lowest delivered tumor coverage, it was found that significant lateral drifts had occurred. The DVFs used for dose warping had a low uncertainty with a mean (SD) ICE of 0.65 (0.07) mm.

CONCLUSION

We employed a fast, real-time 3D cine MRI sequence for dose reconstruction in the upper abdomen, and demonstrated that accurate DVFs, acquired directly from these images, can be used for dose warping. The reconstructed delivered dose showed only a modest degradation of tumor coverage, mostly attainable to baseline drifts. This emphasizes the need for motion monitoring and development of intrafraction treatment adaptation solutions, such as baseline drift corrections.

4D Golden-Angle Radial MRI at Subsecond Temporal Resolution

Intraframe motion blurring, as a major challenge in free-breathing dynamic MRI, can be reduced if high temporal resolution can be achieved. To address this challenge, this work proposes a highly accelerated 4D (3D + time) dynamic MRI framework with subsecond temporal resolution that does not require explicit motion compensation. The method combines standard stack-of-stars golden-angle radial sampling and tailored GRASP-Pro (Golden-angle RAdial Sparse Parallel imaging with imProved performance) reconstruction. Specifically, 4D dynamic MRI acquisition is performed continuously without motion gating or sorting. The k-space centers in stack-of-stars radial data are organized to guide estimation of a temporal basis, with which GRASP-Pro reconstruction is employed to enforce joint low-rank subspace and sparsity constraints. This new basis estimation strategy is the new feature proposed for subspace-based reconstruction in this work to achieve high temporal resolution (e.g., subsecond/3D volume). It does not require sequence modification to acquire additional navigation data, it is compatible with commercially available stack-of-stars sequences, and it does not need an intermediate reconstruction step. The proposed 4D dynamic MRI approach was tested in abdominal motion phantom, free-breathing abdominal MRI, and dynamic contrast-enhanced MRI (DCE-MRI). Our results have shown that GRASP-Pro reconstruction with the new basis estimation strategy enables highly-accelerated 4D dynamic imaging at subsecond temporal resolution (with five spokes or less for each dynamic frame per image slice) for both free-breathing non-DCE-MRI and DCE-MRI. In the abdominal phantom, better image quality with lower root mean square error and higher structural similarity index was achieved using GRASP-Pro compared with standard GRASP. With the ability to acquire each 3D image in less than 1 s, intraframe respiratory blurring can be intrinsically reduced for body applications with our approach, which eliminates the need for explicit motion detection and motion compensation.

A Simulation Study of Tolerance of Breathing Amplitude Variations in Radiotherapy of Lung Cancer Using 4DCT and Time-Resolved 4DMRI

As patient breathing irregularities can introduce a large uncertainty in targeting the internal tumor volume (ITV) of lung cancer patients, and thereby affect treatment quality, this study evaluates dose tolerance of tumor motion amplitude variations in ITV-based volumetric modulated arc therapy (VMAT). A motion-incorporated planning technique was employed to simulate treatment delivery of 10 lung cancer patients' clinical VMAT plans using original and three scaling-up (by 0.5, 1.0, and 2.0 cm) motion waveforms from single-breath four-dimensional computed tomography (4DCT) and multi-breath time-resolved 4D magnetic resonance imaging (TR-4DMRI). The planning tumor volume (PTV = ITV + 5 mm margin) dose coverage (PTV D95%) was evaluated. The repeated waveforms were used to move the isocenter in sync with the clinical leaf motion and gantry rotation. The continuous VMAT arcs were broken down into many static beam fields at the control points (2°-interval) and the composite plan represented the motion-incorporated VMAT plan. Eight motion-incorporated plans per patient were simulated and the plan with the native 4DCT waveform was used as a control. The first (D95% ≤ 95%) and second (D95% ≤ 90%) plan breaching points due to motion amplitude increase were identified and analyzed. The PTV D95% in the motion-incorporated plans was 99.4 ± 1.0% using 4DCT, closely agreeing with the corresponding ITV-based VMAT plan (PTV D95% = 100%). Tumor motion irregularities were observed in TR-4DMRI and triggered D95% ≤ 95% in one case. For small tumors, 4 mm extra motion triggered D95% ≤ 95%, and 6-8 mm triggered D95% ≤ 90%. For large tumors, 14 mm and 21 mm extra motions triggered the first and second breaching points, respectively. This study has demonstrated that PTV D95% breaching points may occur for small tumors during treatment delivery. Clinically, it is important to monitor and avoid systematic motion increase, including baseline drift, and large random motion spikes through threshold-based beam gating.

Improving liver tumor image contrast and synthesizing novel tissue contrasts by adaptive multiparametric MRI fusion

Purpose

Multiparametric MRI contains rich and complementary anatomical and functional information, which is often utilized separately. This study aims to propose an adaptive multiparametric MRI (mpMRI) fusion method and examine its capability in improving tumor contrast and synthesizing novel tissue contrasts among liver cancer patients.

Methods

An adaptive mpMRI fusion method was developed with five components: image pre-processing, fusion algorithm, database, adaptation rules, and fused MRI. Linear-weighted summation algorithm was used for fusion. Weight-driven and feature-driven adaptations were designed for different applications. A clinical-friendly graphic-user-interface (GUI) was developed in Matlab and used for mpMRI fusion. Twelve liver cancer patients and a digital human phantom were included in the study. Synthesis of novel image contrast and enhancement of image signal and contrast were examined in patient cases. Tumor contrast-to-noise ratio (CNR) and liver signal-to-noise ratio (SNR) were evaluated and compared before and after mpMRI fusion.

Results

The fusion platform was applicable in both XCAT phantom and patient cases. Novel image contrasts, including enhancement of soft-tissue boundary, vertebral body, tumor, and composition of multiple image features in a single image were achieved. Tumor CNR improved from -1.70 ± 2.57 to 4.88 ± 2.28 (p < 0.0001) for T1-w, from 3.39 ± 1.89 to 7.87 ± 3.47 (p < 0.01) for T2-w, and from 1.42 ± 1.66 to 7.69 ± 3.54 (p < 0.001) for T2/T1-w MRI. Liver SNR improved from 2.92 ± 2.39 to 9.96 ± 8.60 (p < 0.05) for DWI. The coefficient of variation (CV) of tumor CNR lowered from 1.57, 0.56, and 1.17 to 0.47, 0.44, and 0.46 for T1-w, T2-w and T2/T1-w MRI, respectively.

Conclusion

A multiparametric MRI fusion method was proposed and a prototype was developed. The method showed potential in improving clinically relevant features such as tumor contrast and liver signal. Synthesis of novel image contrasts including the composition of multiple image features into single image set was achieved.

Advances and potential of optical surface imaging in radiotherapy

This article reviews the recent advancements and future potential of optical surface imaging (OSI) in clinical applications as a four-dimensional (4D) imaging modality for surface-guided radiotherapy (SGRT), including OSI systems, clinical SGRT applications, and OSI-based clinical research. The OSI is a non-ionizing radiation imaging modality, offering real-time 3D surface imaging with a large field of view (FOV), suitable for in-room interactive patient setup, and real-time motion monitoring at any couch rotation during radiotherapy. So far, most clinical SGRT applications have focused on treating superficial breast cancer or deep-seated brain cancer in rigid anatomy, because the skin surface can serve as tumor surrogates in these two clinical scenarios, and the procedures for breast treatments in free-breathing (FB) or at deep-inspiration breath-hold (DIBH), and for cranial stereotactic radiosurgery (SRS) and radiotherapy (SRT) are well developed. When using the skin surface as a body-position surrogate, SGRT promises to replace the traditional tattoo/laser-based setup. However, this requires new SGRT procedures for all anatomical sites and new workflows from treatment simulation to delivery. SGRT studies in other anatomical sites have shown slightly higher accuracy and better performance than a tattoo/laser-based setup. In addition, radiographical image-guided radiotherapy (IGRT) is still necessary, especially for stereotactic body radiotherapy (SBRT). To go beyond the external body surface and infer an internal tumor motion, recent studies have shown the clinical potential of OSI-based spirometry to measure dynamic tidal volume as a tumor motion surrogate, and Cherenkov surface imaging to guide and assess treatment delivery. As OSI provides complete datasets of body position, deformation, and motion, it offers an opportunity to replace fiducial-based optical tracking systems. After all, SGRT has great potential for further clinical applications. In this review, OSI technology, applications, and potential are discussed since its first introduction to radiotherapy in 2005, including technical characterization, different commercial systems, and major clinical applications, including conventional SGRT on top of tattoo/laser-based alignment and new SGRT techniques attempting to replace tattoo/laser-based setup. The clinical research for OSI-based tumor tracking is reviewed, including OSI-based spirometry and OSI-guided tumor tracking models. Ongoing clinical research has created more SGRT opportunities for clinical applications beyond the current scope.

Real-Time 2D MR Cine From Beam Eye's View With Tumor-Volume Projection to Ensure Beam-to-Tumor Conformality for MR-Guided Radiotherapy of Lung Cancer

Purpose

To minimize computation latency using a predictive strategy to retrieve and project tumor volume onto 2D MR beam eye’s view (BEV) cine from time-resolved four-dimensional magnetic resonance imaging (TR-4DMRI) libraries (inhalation/exhalation) for personalized MR-guided intensity-modulated radiotherapy (IMRT) or volumetric-modulated arc therapy (VMAT).

Methods

Two time-series forecasting algorithms, autoregressive (AR) modeling and deep-learning-based long short-term memory (LSTM), were applied to predict the diaphragm position in the next 2D BEV cine to identify a motion-matched and hysteresis-accounted image to retrieve the tumor volume from the inhalation/exhalation TR-4DMRI libraries. Three 40-s TR-4DMRI (2 Hz, 3 × 80 images) per patient of eight lung cancer patients were used to create patient-specific inhalation/exhalation 4DMRI libraries, extract diaphragmatic waveforms, and interpolate them to f = 4 and 8 Hz to match 2D cine frame rates. Along a (40•f)-timepoint waveform, 30•f training timepoints were moved forward to produce 3×(10•f-1) predictions. The accuracy of position prediction was assessed against the waveform ground truth. The accuracy of tumor volume projections was evaluated using the center-of-mass difference (∆COM) and Dice similarity index against the TR-4DMRI ground truth for both IMRT (six beam angles, 30° interval) and VMAT (240/480 beam angles, 1.5°/0.75° interval, at 4/8 Hz, respectively).

Results

The accuracy of the first-timepoint prediction is 0.36 ± 0.10 mm (AR) and 0.62 ± 0.21 mm (LSTM) at 4 Hz and 0.06 ± 0.02 mm (AR) and 0.18 ± 0.06 mm (LSTM) at 8 Hz. A 10%–20% random error in prediction-library matching increases the overall uncertainty slightly. For both IMRT and VMAT, the accuracy of projected tumor volume contours on 2D BEV cine is ∆COM = 0.39 ± 0.13 mm and DICE = 0.97 ± 0.02 at 4 Hz and ∆COM = 0.10 ± 0.04 mm and DICE = 1.00 ± 0.00 at 8Hz.

Conclusion

This study demonstrates the feasibility of accurately predicting respiratory motion during 2D BEV cine imaging, identifying a motion-matched and hysteresis-accounted tumor volume, and projecting tumor volume contour on 2D BEV cine for real-time assessment of beam-to-tumor conformality, promising for optimal personalized MR-guided radiotherapy.

Stability and Reliability of Enhanced External-Internal Motion Correlation via Dynamic Phase-Shift Corrections Over 30-min Timeframe for Respiratory-Gated Radiotherapy

To assess the stability of patient-specific phase shifts between external- and internal-respiratory motion waveforms, the reliability of enhanced external–internal correlation with phase-shift correction, and the feasibility of guiding respiratory-gated radiotherapy (RGRT) over 30 min. In this clinical feasibility investigation, external bellows and internal-navigator waveforms were simultaneously and prospectively acquired along with two four-dimensional magnetic resonance imaging (4DMRI) scans (6–15 m each) with 15–20 m intervals in 10 volunteers. A bellows was placed 5 cm inferior to the xiphoid to monitor abdominal motion, and an MR navigator was used to track the diaphragmatic motion. The mean phase-domain (MPD) method was applied, which combines three individual phase-calculating methods: phase-space oval fitting, principal component analysis, and analytic signal analysis, weighted by the reciprocal of their residual errors (RE) excluding outliers (RE >2σ). The time-domain cross-correlation (TCC) analysis was applied for comparison. Dynamic phase-shift correction was performed based on the phase shift detected on the fly within two 10 s moving datasets. Simulating bellows-triggered gating, the median and 95% confidence interval for the navigator's position at beam-on/beam-off and %harm (percentage of beam-on time outside the safety margin) were calculated. Averaged across all subjects, the mean phase shifts are found indistinguishable (p > .05) between scan 1 (55˚ ± 9˚) and scan 2 (59˚ ± 11˚). Using the MPD method the averaged correlation increases from 0.56 ± 0.22 to 0.85 ± 0.11 for scan 1 and from 0.47 ± 0.30 to 0.84 ± 0.08 for scan 2. The TCC correction results in similar results. After phase-shift correction, the number of cases that were suitable for amplitude gating (with <10%harm) increased from 2 to 17 out of 20 cases. A patient-specific, stable phase-shift between the external and internal motions was observed and corrected using the MPD and TCC methods, producing long-lasting enhanced motion correlation over 30m. Phase-shift correction offers a feasible strategy for improving the accuracy of tumor-motion prediction during RGRT.

Multi-contrast four-dimensional magnetic resonance imaging (MC-4D-MRI): Development and initial evaluation in liver tumor patients

PURPOSE

To develop a novel multi-contrast four-dimensional magnetic resonance imaging (MC-4D-MRI) technique that expands single image contrast 4D-MRI to a spectrum of native and synthetic image contrasts and to evaluate its feasibility in liver tumor patients.

METHODS AND MATERIALS

The MC-4D-MRI technique integrates multi-parametric MRI fusion, 4D-MRI, and deformable image registration (DIR) techniques. The fusion technique consists of native MRI as input, image pre-processing, fusion algorithm, adaptation, and fused multi-contrast MRI as output. Four-dimensional deformation vector fields (4D-DVF) were generated from an original T2/T1-w 4D-MRI by deforming end-of-inhalation (EOI) to nine other phase volumes via DIR. The 4D-DVF were applied to multi-contrast MRI to generate a spectrum of 4D-MRI in different image contrasts. The MC-4D-MRI technique was evaluated in five liver tumor patients on tumor contrast-to-noise ratio (CNR), internal target volume (ITV) contouring consistency, diaphragm motion range, and tumor motion trajectory; and in digital anthropomorphic phantoms on 4D-DIR introduced errors in tumor motion range, centroid location, extent, and volume.

RESULTS

MC-4D-MRI consisting of 4D-MRIs in native image contrasts (T1-w, T2-w, and T2/T1-w) and synthetic image contrasts, such as tumor-enhanced contrast (TEC) were generated in five liver tumor patients. Patient tumor CNR increased from 2.6 ± 1.8 in the T2/T1-w MRI, to -4.4 ± 2.4, 6.6 ± 3.0, and 9.6 ± 3.9 in the T1-w, T2-w, and TEC MRI, respectively. Patient ITV inter-observer mean Dice similarity coefficient (mDSC) increased from 0.65 ± 0.10 in the original T2/T1-w 4D-MRI, to 0.76 ± 0.14, 0.77 ± 0.12, and 0.86 ± 0.05 in the T1-w, T2-w, and TEC 4D-MRI, respectively. Patient diaphragm motion range absolute differences between the three new 4D-MRIs and original T2/T1-w 4D-MRI were 1.2 ± 1.3, 0.3 ± 0.7, and 0.5 ± 0.5 mm, respectively. Patient tumor displacement phase-averaged absolute differences between the three 4D-MRIs and the original 4D-MRI were 0.72 ± 0.33, 0.62 ± 0.54, and 0.74 ± 0.43 mm in the superior-inferior (SI) direction, and 0.59 ± 0.36, 0.51 ± 0.30, and 0.50 ± 0.24 mm in the anterior-posterior (AP) direction, respectively. In the digital phantoms, phase-averaged absolute tumor centroid shift caused by the 4D-DIR were at or below 0.5 mm in SI, AP, and left-right (LR) directions.

CONCLUSION

We developed an MC-4D-MRI technique capable of expanding single image contrast 4D-MRI along a new dimension of image contrast. Initial evaluations in liver tumor patients showed enhancements in image contrast variety, tumor contrast, and ITV contouring consistencies using MC-4D-MRI. The technique might offer new perspectives on the image contrast of MRI and 4D-MRI in MR-guided radiotherapy.

Evaluation of super-resolution on 50 pancreatic cancer patients with real-time cine MRI from 0.35T MRgRT

MR-guided radiotherapy (MRgRT) systems provide excellent soft tissue imaging immediately prior to and in real time during radiation delivery for cancer treatment. However, 2D cine MRI often has limited spatial resolution due to high temporal resolution. This work applies a super resolution machine learning framework to 3.5 mm pixel edge length, low resolution (LR), sagittal 2D cine MRI images acquired on a MRgRT system to generate 0.9 mm pixel edge length, super resolution (SR), images originally acquired at 4 frames per second (FPS). LR images were collected from 50 pancreatic cancer patients treated on a ViewRay MR-LINAC. SR images were evaluated using three methods. 1) The first method utilized intrinsic image quality metrics for evaluation. 2) The second used relative metrics including edge detection and structural similarity index (SSIM). 3) Finally, automatically generated tumor contours were created on both low resolution and super resolution images to evaluate target delineation and compared with DICE and SSIM. Intrinsic image quality metrics all had statistically significant improvements for SR images versus LR images, with mean (±1 SD) BRISQUE scores of 29.65 ± 2.98 and 42.48 ± 0.98 for SR and LR, respectively. SR images showed good agreement with LR images in SSIM evaluation, indicating there was not significant distortion of the images. Comparison of LR and SR images with paired high resolution (HR) 3D images showed that SR images had a mean (±1 SD) SSIM value of 0.633 ± 0.063 and LR a value of 0.587 ± 0.067 (p ≪ 0.05). Contours generated on SR images were also more robust to noise addition than those generated on LR images. This study shows that super resolution with a machine learning framework can generate high spatial resolution images from 4fps low spatial resolution cine MRI acquired on the ViewRay MR-LINAC while maintaining tumor contour quality and without significant acquisition or post processing delay.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Enhanced super-resolution reconstruction of T1w time-resolved 4DMRI in low-contrast tissue using 2-step hybrid deformable image registration

Purpose

Deformable image registration (DIR) in low‐contrast tissues is often suboptimal because of low visibility of landmarks, low driving‐force to deform, and low penalty for misalignment. We aim to overcome the shortcomings for improved reconstruction of time‐resolved four‐dimensional magnetic resonance imaging (TR‐4DMRI).

Methods and Materials

Super‐resolution TR‐4DMRI reconstruction utilizes DIR to combine high‐resolution (highR:2x2x2mm³) breath‐hold (BH) and low‐resolution (lowR:5x5x5mm³) free‐breathing (FB) 3D cine (2Hz) images to achieve clinically acceptable spatiotemporal resolution. A 2‐step hybrid DIR approach was developed to segment low‐dynamic‐range (LDR) regions: low‐intensity lungs and high‐intensity “bodyshell” (=body‐lungs) for DIR refinement after conventional DIR. The intensity in LDR regions was renormalized to the full dynamic range (FDR) to enhance local tissue contrast. A T1‐mapped 4D XCAT digital phantom was created, and seven volunteers and five lung cancer patients were scanned with two BH and one 3D cine series per subject to compare the 1‐step conventional and 2‐step hybrid DIR using: (a) the ground truth in the phantom, (b) highR‐BH references, which were used to simulate 3D cine images by down‐sampling and Rayleigh‐noise‐adding, and (c) cross‐verification between two TR‐4DMRI images reconstructed from two BHs. To assess DIR improvement, 8‐17 blood vessel bifurcations were used in volunteers, and lung tumor position, size, and shape were used in phantom and patients, together with the voxel intensity correlation (VIC), structural similarity (SSIM), and cross‐consistency check (CCC).

Results

The 2‐step hybrid DIR improves contrast and DIR accuracy. In volunteers, it improves low‐contrast alignment from 6.5 ± 1.8 mm to 3.3 ± 1.0 mm. In phantom, it improves tumor center of mass alignment (COM = 1.3 ± 0.2 mm) and minimizes DIR directional difference. In patients, it produces almost‐identical tumor COM, size, and shape (dice> 0.85) as the reference. The VIC and SSIM are significantly increased and the number of CCC outliers are reduced by half.

Conclusion

The 2‐step hybrid DIR improves low‐contrast‐tissue alignment and increases lung tumor fidelity. It is recommended to adopt the 2‐step hybrid DIR for TR‐4DMRI reconstruction.

A modern review of the uncertainties in volumetric imaging of respiratory-induced target motion in lung radiotherapy

Radiotherapy has become a critical component for the treatment of all stages and types of lung cancer, often times being the primary gateway to a cure. However, given that radiation can cause harmful side effects depending on how much surrounding healthy tissue is exposed, treatment of the lung can be particularly challenging due to the presence of moving targets. Careful implementation of every step in the radiotherapy process is absolutely integral for attaining optimal clinical outcomes. With the advent and now widespread use of stereotactic body radiation therapy (SBRT), where extremely large doses are delivered, accurate, and precise dose targeting is especially vital to achieve an optimal risk to benefit ratio. This has largely become possible due to the rapid development of image-guided technology. Although imaging is critical to the success of radiotherapy, it can often be plagued with uncertainties due to respiratory-induced target motion. There has and continues to be an immense research effort aimed at acknowledging and addressing these uncertainties to further our abilities to more precisely target radiation treatment. Thus, the goal of this article is to provide a detailed review of the prevailing uncertainties that remain to be investigated across the different imaging modalities, as well as to highlight the more modern solutions to imaging motion and their role in addressing the current challenges.

A super-resolution framework for the reconstruction of T2-weighted (T2w) time-resolved (TR) 4DMRI using T1w TR-4DMRI as the guidance

PURPOSE

The purpose of this study was to develop T2-weighted (T2w) time-resolved (TR) four-dimensional magnetic resonance imaging (4DMRI) reconstruction technique with higher soft-tissue contrast for multiple breathing cycle motion assessment by building a super-resolution (SR) framework using the T1w TR-4DMRI reconstruction as guidance.

METHODS

The multi-breath T1w TR-4DMRI was reconstructed by deforming a high-resolution (HR: 2 × 2 × 2 mm³ ) volumetric breath-hold (BH, 20s) three-dimensional magnetic resonance imaging (3DMRI) image to a series of low-resolution (LR: 5 × 5 × 5 mm³ ) 3D cine images at a 2Hz frame rate in free-breathing (FB, 40 s) using an enhanced Demons algorithm, namely [T1^BH →FB] reconstruction. Within the same imaging session, respiratory-correlated (RC) T2w 4DMRI (2 × 2 × 2 mm³ ) was acquired based on an internal navigator to gain HR T2w (T2^HR ) in three states (full exhalation and mid and full inhalation) in ~5 min. Minor binning artifacts in the RC-4DMRI were automatically identified based on voxel intensity correlation (VIC) between consecutive slices as outliers (VIC < VIC^mean -σ) and corrected by deforming the artifact slices to interpolated slices from the adjacent slices iteratively until no outliers were identified. A T2^HR image with minimal deformation (<1 cm at the diaphragm) from the T1^BH image was selected for multi-modal B-Spline deformable image registration (DIR) to establish the T2^HR -T1^BH voxel correspondence. Two approaches to reconstruct T2w TR-4DMRI were investigated: (A) T2^HR →[T1^BH →FB]: to deform T2w HR to T1w BH only as T1w TR-4DMRI was reconstructed, and combine the two displacement vector fields (DVFs) to reconstruct T2w TR-4DMRI, and (B) [T2^HR ←T1^BH ]→FB: to deform T1w BH to T2w HR first and apply the deformed T1w BH to reconstruct T2w TR-4DMRI. The reconstruction times were similar, 8-12 min per volume. To validate the two methods, T2w- and T1w-mapped 4D XCAT digital phantoms were utilized with three synthetic spherical tumors (ϕ = 2.0, 3.0, and 4.0 cm) in the lower or mid lobes as the ground truth to evaluate the tumor location (the center of mass, COM), size (volume ratio, %V), and shape (Dice index). Six lung cancer patients were scanned under an IRB-approved protocol and the T2w TR-4DMRI images reconstructed from the two methods were compared based on the preservation of the three tumor characteristics. The local tumor-contained image quality was also characterized using the VIC and structure similarity (SSIM) indexes.

RESULTS

In the 4D digital phantom, excellent tumor alignment after T2^HR -T1^HR DIR is achieved: ∆COM = 0.8 ± 0.5 mm, %V = 1.06 ± 0.02, and Dice = 0.91 ± 0.03, in both deformation directions using the DIR-target image as the reference. In patients, binning artifacts are corrected with improved image quality: average VIC increases from 0.92 ± 0.03 to 0.95 ± 0.01. Both T2w TR-4DMRI reconstruction methods produce similar tumor alignment errors ∆COM = 2.9 ± 0.6 mm. However, method B ([T2^HR ←T1^BH ]→FB) produces superior results in preserving more T2w tumor features with a higher %V = 0.99 ± 0.03, Dice = 0.81 ± 0.06, VIC = 0.85 ± 0.06, and SSIM = 0.65 ± 0.10 in the T2w TR-4DMRI images.

CONCLUSIONS

This study has demonstrated the feasibility of T2w TR-4DMRI reconstruction with high soft-tissue contrast and adequately-preserved tumor position, size, and shape in multiple breathing cycles. The T2w-centric DIR (method B) produces a superior solution for the SR-based framework of T2w TR-4DMRI reconstruction with highly preserved tumor characteristics and local image features, which are useful for tumor delineation and motion management in radiation therapy.

Clinical evaluation of 4D MRI in the delineation of gross and internal tumor volumes in comparison with 4DCT

Purpose

To evaluate clinical utility of respiratory‐correlated (RC) four‐dimensional magnetic resonance imaging (4DMRI) for lung tumor delineation and motion assessment, in comparison with the current clinical standard of 4D computed tomography (4DCT).

Methods and Materials

A prospective T2‐weighted (T2w) RC‐4DMRI technique was applied to acquire coronal 4DMRI images for 14 lung cancer patients (16 lesions) during free breathing (FB) under an IRB‐approved protocol, together with a breath‐hold (BH) T1w 3DMRI and axial 4DMRI. Clinical simulation CT and 4DCT were acquired within 2 h. An internal navigator was applied to trigger amplitude‐binned 4DMRI acquisition whereas a bellows or real‐time position management (RPM) was used in the 4DCT reconstruction. Six radiation oncologists manually delineated the gross and internal tumor volumes (GTV and ITV) in 399 3D images using programmed clinical workflows under a tumor delineation guideline. The ITV was the union of GTVs within the breathing cycle without margin. Average GTV and motion range were assessed and ITV variation between 4DMRI and 4DCT was evaluated using the Dice similarity index, mean distance agreement (MDA), and volume difference.

Results

The mean tumor volume is similar between 4DCT (GTV^4DCT = 1.0, as the reference) and T2w‐4DMRI (GTV^T2wMR = 0.97), but smaller in T1w MRI (GTV^T1wMR = 0.76), suggesting possible peripheral edema around the tumor. Average GTV variation within the breathing cycle (22%) in 4DMRI is slightly greater than 4DCT (17%). GTV motion variation (−4 to 12 mm) and ITV variation (∆V^ITV=−25 to 95%) between 4DCT and 4DMRI are large, confirmed by relatively low ITV similarity (Dice = 0.72 ± 0.11) and large MDA = 2.9 ± 1.5 mm.

Conclusion

Average GTVs are similar between T2w‐4DMRI and 4DCT, but smaller by 25% in T1w BH MRI. Physician training and breathing coaching may be necessary to reduce ITV variability between 4DMRI and 4DCT. Four‐dimensional magnetic resonance imaging is a promising and viable technique for clinical lung tumor delineation and motion assessment.

An image regression motion prediction technique for MRI-guided radiotherapy evaluated in single-plane cine imaging

PURPOSE

To develop and evaluate a novel motion prediction method for magnetic resonance image (MRI)-guided radiotherapy applications. This method, which we deem "image regression," predicts future tissue motion based on a weighted combination of previously observed motion states. Motion predictions are derived from a sliding window of recent motion states which are defined by a temporal sequence of images. A key advantage of this method compared to other motion prediction methods is that its computational complexity scales weekly with the number of spatial points predicted. Applications of gating latency reduction and improvement in deformable registration-based target tracking are demonstrated.

METHODS

The image regression (IR) motion prediction method was developed and evaluated using 26.9 h of real-time imaging acquired from eight healthy volunteers and 13 patients using a 0.35 T MRI-guided radiotherapy system. Motion predictions were performed 0.25-0.33 s into the future using a weighted sum of previously observed motion states with image similarity-derived weights. The set of previously observed motion states were continuously updated to incorporate the changes in breathing patterns. The accuracy of the predicted radiotherapy gating decision, beam-on positive predictive value (PPV), and predicted vs ground-truth target centroid position errors are reported. The IR technique was compared against no prediction, linear extrapolation, and an established autoregressive linear prediction algorithm. The usage of IR to initialize the deformable registration and enhance the target tracking was demonstrated in the healthy volunteer studies. Deformable registration with IR initialization was compared to the initialization performed by current clinical software: no initialization, previous image registration initialization and linear motion extrapolation initialization.

RESULTS

The average IR-predicted radiation beam gating decision accuracy was 95.8%, with a PPV of 95.7%, and median and 95th percentile centroid position errors of 0.63 and 2.08 mm, respectively. Compared to the autoregressive linear prediction method, gating accuracy was 1.15% greater, PPV was 1.61% greater, and median and 95th percentile centroid distances were 0.21 and 0.23 mm smaller. The IR-initialized registration on average converged within 0.50 mm of the ground-truth position in fewer than 10 iterations whereas the next best initialization method required more than 25 iterations.

CONCLUSIONS

Image regression motion prediction has the potential to reduce the gating latencies and improve the speed and accuracy of deformable registration-based target tracking in MRI-guided radiotherapy.

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.

Enhancement of Long-Term External-Internal Correlation by Phase-Shift Detection and Correction Based on Concurrent External Bellows and Internal Navigator Signals

Purpose

The purpose of this study was to enhance the correlation between external and internal respiratory motions by dynamically determining and correcting the patient-specific phase shift between external and internal respiratory waveforms acquired concurrently during respiratory-correlated 4-dimensional magnetic resonance imaging scans.

Methods and Materials

Internal-navigator and external-bellows waveforms were acquired simultaneously during 6- to 15-minute respiratory-correlated 4-dimensional magnetic resonance imaging scans in 10 healthy participants under an institutional review board–approved protocol. The navigator was placed at the right lung–diaphragm interface, and the bellows were placed ∼5 cm inferior to the sternum. Three segments of each respiratory waveform, at the beginning, middle, and end of a scan, were analyzed. Three phase-domain methods were employed to estimate the phase shift, including analytical signal analysis, phase-space oval fitting, and principal component analysis. A robust strategy for estimating the phase shift was realized by combining these methods in a weighted average and by eliminating outliers (>2 σ) caused by breathing irregularities. Whether phase-shift correction affects the external-internal correlation was evaluated. The cross-correlation between the 2 waveforms in the time domain provided an independent check of the correlation enhancement.

Results

Phase-shift correction significantly enhanced the external-internal correlation in all participants across the entire 6- to 15-minute scans. On average, the correlation increased from 0.45 ± 0.28 to 0.85 ± 0.15 for the combined method. The combined method exhibited a 99.5% success rate and revealed that the phase of the external waveform leads that of the internal waveform in all 10 participants by 57 ^o ± 17^o (1.6 ± 0.5 bins) on average. Seven participants exhibited highly reproducible phase shifts over time, evidenced by standard deviations (σ) < 4^o, whereas 8^o < σ < 12^o in the remaining 3 participants. Regardless, phase-shift correction significantly improved the correlation in all participants.

Conclusions

Correcting the phase shift estimated by the phase-domain methods provides a new approach for enhancing the correlation between external and internal respiratory motions. This strategy holds promise for improving the accuracy of respiratory-gated radiation therapy.

Multislice motion modeling for MRI-guided radiotherapy gating

PURPOSE

On-board magnetic resonance imaging (MRI) greatly enhances real-time target tracking capability during radiotherapy treatments. However, multislice and volumetric MRI techniques are frame rate limited and introduce unacceptable latency between the target moving out of position and the beam being turned off. We present a technique to estimate continuous volumetric tissue motion using motion models built from a repeated acquisition of a stack of MR slices. Applications including multislice target visualization and out-of-slice motion estimation during MRI-guided radiotherapy are demonstrated.

METHODS

Eight healthy volunteer studies were performed using a 0.35 T MRI-guided radiotherapy system. Images were acquired at three frames per second in an interleaved fashion across ten adjacent sagittal slice positions covering 4.5 cm using a balanced steady-state-free precession sequence. A previously published five-dimensional (5D) linear motion model used for MRI-guided radiotherapy gating was extended to include multiple slices. This model utilizes an external respiratory bellows signal recorded during imaging to simultaneously estimate motion across all imaged slices. For comparison to an image-based approach, the manifold learning technique local linear embedding (LLE) was used to derive a respiratory surrogate for motion modeling. Manifolds for every slice were aligned during LLE in a group-wise fashion, enabling motion estimation outside the current imaged slice using a motion model, a process which we denote as mSGA. Additionally, a method is developed to evaluate out-of-slice motion estimates. The multislice motion model was evaluated in a single slice with each newly acquired image using a leave-one-out approach. Model-generated gating decision accuracy and beam-on positive predictive value (PPV) are reported along with the median and 95th percentile distance between model and ground truth target centroids.

RESULTS

The average model gating decision accuracy and PPV across all volunteer studies was 93.7% and 92.8% using the 5D model, and 96.8% and 96.1% using the mSGA model, respectively. The median and 95th percentile distance between model and ground truth target centroids was 0.91 and 2.90 mm, respectively, using the 5D model and 0.58 and 1.49 mm using the mSGA model, averaged over all eight subjects. The mSGA motion model provided a statistically significant improvement across all evaluation metrics compared to the external surrogate-based 5D model.

CONCLUSION

The proposed techniques for out-of-slice target motion estimation demonstrated accuracy likely sufficient for clinical use. Results indicate the mSGA model may provide higher accuracy, however, the external surrogate-based model allows for unbiased in vivo accuracy evaluation.

Super-resolution T2-weighted 4D MRI for image guided radiotherapy

BACKGROUND AND PURPOSE

The superior soft-tissue contrast of 4D-T2w MRI motivates its use for delineation in radiotherapy treatment planning. We address current limitations of slice-selective implementations, including thick slices and artefacts originating from data incompleteness and variable breathing.

MATERIALS AND METHODS

A method was developed to calculate midposition and 4D-T2w images of the whole thorax from continuously acquired axial and sagittal 2D-T2w MRI (1.5 × 1.5 × 5.0 mm3). The method employed image-derived respiratory surrogates, deformable image registration and super-resolution reconstruction. Volunteer imaging and a respiratory motion phantom were used for validation. The minimum number of dynamic acquisitions needed to calculate a representative midposition image was investigated by retrospectively subsampling the data (10-30 dynamic acquisitions).

RESULTS

Super-resolution 4D-T2w MRI (1.0 × 1.0 × 1.0 mm3, 8 respiratory phases) did not suffer from data incompleteness and exhibited reduced stitching artefacts compared to sorted multi-slice MRI. Experiments using a respiratory motion phantom and colour-intensity projection images demonstrated a minor underestimation of the motion range. Midposition diaphragm differences in retrospectively subsampled acquisitions were <1.1 mm compared to the full dataset. 10 dynamic acquisitions were found sufficient to generate midposition MRI.

CONCLUSIONS

A motion-modelling and super-resolution method was developed to calculate high quality 4D/midposition T2w MRI from orthogonal 2D-T2w MRI.

Model-Interpolated Gating for Magnetic Resonance Image-Guided Radiation Therapy

PURPOSE

To develop and validate a technique for radiation therapy gating using slow (≤1 frame per second) magnetic resonance imaging (MRI) and a motion model. Proposed uses of the technique include radiation therapy gating using T2-weighted images and conducting additional imaging studies during gated treatments.

METHODS AND MATERIALS

The technique uses a physiologically guided breathing motion model to interpolate deformed target position between 2-dimensional (2D) MRI images acquired every 1 to 3 seconds. The model is parameterized by a 1-dimensional respiratory bellows surrogate and is continuously updated with the most recently acquired 2D images. A phantom and 8 volunteers were imaged with a 0.35T MRI-guided radiation therapy system. A balanced steady-state free precession sequence with a 2D frame rate of 3 frames per second was used to evaluate the technique. The accuracy and beam-on positive predictive value (PPV) of the model-based gating decisions were evaluated using the gating decisions derived from imaging as a ground truth. A T2-weighted gating offline proof-of-concept study using a half-Fourier, single-shot, turbo-spin echo sequence is reported.

RESULTS

Model-interpolated gating accuracy, beam-on PPV, and median absolute distances between model and image-tracked target centroids were, on average, 98.3%, 98.4%, and 0.33 mm, respectively, in the balanced steady-state free precession phantom studies and 93.7%, 92.1%, and 0.86 mm, respectively, in the volunteer studies. T2 model-interpolated gating in 6 volunteers yielded an average accuracy and PPV of 94.3% and 92.5%, respectively, and the mean absolute median distance between modeled and imaged target centroids was 0.86 mm.

CONCLUSIONS

This work demonstrates the concept of model-interpolated gating for MRI-guided radiation therapy. The technique was found to be potentially sufficiently accurate for clinical use. Further development is needed to accommodate out-of-plane motion and the use of an internal MR-based respiratory surrogate.

Introduction of a pseudo demons force to enhance deformation range for robust reconstruction of super-resolution time-resolved 4DMRI

PURPOSE

The purpose of this study was to enhance the deformation range of demons-based deformable image registration (DIR) for large respiration-induced organ motion in the reconstruction of time-resolved four-dimensional magnetic resonance imaging (TR-4DMRI) for multi-breath motion simulation.

METHODS

A demons-based DIR algorithm was modified to enhance the deformation range for TR-4DMRI reconstruction using the super-resolution approach. A pseudo demons force was introduced to accelerate the coarse deformation in a multi-resolution (n = 3) DIR approach. The intensity gradient of a voxel was applied to its neighboring (5 × 5 × 5) voxels with a weight of Gaussian probability profile (σ = 1 voxel) to extend the demons force, especially on those voxels that have little intensity gradience but high-intensity difference. A digital 4DMRI phantom with 3-8 cm diaphragmatic motions was used for DIR comparison. Six volunteers were scanned with two high-resolution (highR: 2 × 2 × 2 mm³ ) breath-hold (BH) 3DMR images at full inhalation (BHI) and full exhalation (BHE) and low-resolution (lowR: 5 × 5 × 5 mm³ ) free-breathing (FB) 3DMR cine images (2 Hz) under an IRB-approved protocol. A cross-consistency check (CCC) (BHI→FB←BHE), with voxel intensity correlation (VIC) and inverse consistency error (ICE), was introduced for cross-verification of TR-4DMRI reconstruction.

RESULTS

Using the digital phantom, the maximum deformable magnitude is doubled using the modified DIR from 3 to 6 cm at the diaphragm. In six human subjects, the first 15-iteration DIR using the pseudo force deforms 200 ± 150% more than the original force, and succeeds in all 12 cases, whereas the original demons-based DIR failed in 67% of tested cases. Using the pseudo force, high VIC (>0.9) and small ICE (1.6 ± 0.6 mm) values are observed for DIR of BHI&BHE, BHI→FB, and BHE→FB. The CCC identifies four questionable cases, in which two cases need further DIR refinement, without missing true negative.

CONCLUSIONS

The introduction of a pseudo demons force enhances the largest deformation magnitude up to 6 cm. The cross-consistency check ensures the quality of TR-4DMRI reconstruction. Further investigation is ongoing to fully characterize TR-4DMRI for potential multi-breathing-cycle radiotherapy simulation.

Evaluation of automatic contour propagation in T2-weighted 4DMRI for normal-tissue motion assessment using internal organ-at-risk volume (IRV)

Purpose

The purpose of this study was to evaluate the quality of automatically propagated contours of organs at risk (OARs) based on respiratory‐correlated navigator‐triggered four‐dimensional magnetic resonance imaging (RC‐4DMRI) for calculation of internal organ‐at‐risk volume (IRV) to account for intra‐fractional OAR motion.

Methods and Materials

T2‐weighted RC‐4DMRI images were of 10 volunteers acquired and reconstructed using an internal navigator‐echo surrogate and concurrent external bellows under an IRB‐approved protocol. Four major OARs (lungs, heart, liver, and stomach) were delineated in the 10‐phase 4DMRI. Two manual‐contour sets were delineated by two clinical personnel and two automatic‐contour sets were propagated using free‐form deformable image registration. The OAR volume variation within the 10‐phase cycle was assessed and the IRV was calculated as the union of all OAR contours. The OAR contour similarity between the navigator‐triggered and bellows‐rebinned 4DMRI was compared. A total of 2400 contours were compared to the most probable ground truth with a 95% confidence level (S95) in similarity, sensitivity, and specificity using the simultaneous truth and performance level estimation (STAPLE) algorithm.

Results

Visual inspection of automatically propagated contours finds that approximately 5–10% require manual correction. The similarity, sensitivity, and specificity between manual and automatic contours are indistinguishable (P > 0.05). The Jaccard similarity indexes are 0.92 ± 0.02 (lungs), 0.89 ± 0.03 (heart), 0.92 ± 0.02 (liver), and 0.83 ± 0.04 (stomach). Volume variations within the breathing cycle are small for the heart (2.6 ± 1.5%), liver (1.2 ± 0.6%), and stomach (2.6 ± 0.8%), whereas the IRV is much larger than the OAR volume by: 20.3 ± 8.6% (heart), 24.0 ± 8.6% (liver), and 47.6 ± 20.2% (stomach). The Jaccard index is higher in navigator‐triggered than bellows‐rebinned 4DMRI by 4% (P < 0.05), due to the higher image quality of navigator‐based 4DMRI.

Conclusion

Automatic and manual OAR contours from Navigator‐triggered 4DMRI are not statistically distinguishable. The navigator‐triggered 4DMRI image provides higher contour quality than bellows‐rebinned 4DMRI. The IRVs are 20–50% larger than OAR volumes and should be considered in dose estimation.