Respiration-based sorting of dynamic MRI to derive representative 4D-MRI for radiotherapy planning

Retrospective reconstruction of four-dimensional magnetic resonance from interleaved cine imaging - A comparative study with four-dimensional computed tomography in the lung

Background and purpose

Imaging of respiration-induced anatomical changes is essential to ensure high accuracy in radiotherapy of lung cancer. We expanded here on methods for retrospective reconstruction of time-resolved volumetric magnetic resonance (4DMR) of the thoracic region and benchmarked the results against 4D computed tomography (4DCT).

Materials and method

MR data of six lung cancer patients were collected by interleaving cine-navigator images with 2D data frame images, acquired across the thorax. The data frame images have been stacked in volumes based on a similarity metric that considers the anatomical deformation of lungs, while addressing ambiguities in respiratory phase detection and interpolation of missing data. The resulting images were validated against cine-navigator images and compared to paired 4DCTs in terms of amplitude and period of motion, assessing differences in internal target volume (ITV) margin definition.

Results

4DMR-based motion amplitude was on average within 1.8 mm of that measured in the corresponding 2D cine-navigator images. In our dataset, the 4DCT motion and the 4DMR median amplitude were always within 3.8 mm. The median period was generally close to CT references, although deviations up to 24 % have been observed. These changes were reflected in the ITV, which was generally larger for MRI than for 4DCT (up to 39.7 %).

Conclusions

The proposed algorithm for retrospective reconstruction of time-resolved volumetric MR provided quality anatomical images with high temporal resolution for motion modelling and treatment planning. The potential for imaging organ motion variability makes 4DMR a valuable complement to standard 4DCT imaging.

Motion-resolved four-dimensional abdominal diffusion-weighted imaging using PROPELLER EPI (4D-DW-PROPELLER-EPI)

PURPOSE

To develop a distortion-free motion-resolved four-dimensional diffusion-weighted PROPELLER EPI (4D-DW-PROPELLER-EPI) technique for benefiting clinical abdominal radiotherapy (RT).

METHODS

An improved abdominal 4D-DWI technique based on 2D diffusion-weighted PROPELLER-EPI (2D-DW-PROPELLER-EPI), termed 4D-DW-PROPELLER-EPI, was proposed to improve the frame rate of repeated data acquisition and produce distortion-free 4D-DWI images. Since the radial or PROPELLER sampling with golden-angle rotation can achieve an efficient k-space coverage with a flexible time-resolved acquisition, the golden-angle multi-blade acquisition was used in the proposed 4D-DW-PROPELLER-EPI to improve the performance of data sorting. A new k-space and blade (K-B) amplitude binning method was developed for the proposed 4D-DW-PROPELLER-EPI to optimize the number of blades and the k-space uniformity before performing conventional PROPELLER-EPI reconstruction, by using two metrics to evaluate the adequacy of the acquired data. The proposed 4D-DW-PROPELLER-EPI was preliminarily evaluated in both simulation experiments and in vivo experiments with varying frame rates and different numbers of repeated acquisition.

RESULTS

The feasibility of achieving distortion-free 4D-DWI images by using the proposed 4D-DW-PROPELLER-EPI technique was demonstrated in both digital phantom and healthy subjects. Evaluation of the 4D completeness metrics shows that the K-B amplitude binning method could simultaneously improve the acquisition efficiency and data reconstruction performance for 4D-DW-PROPELLER-EPI.

CONCLUSION

4D-DW-PROPELLER-EPI with K-B amplitude binning is an advanced technique that can provide distortion-free 4D-DWI images for resolving respiratory motion, and may benefit the application of image-guided abdominal RT.

Transfer-learning is a key ingredient to fast deep learning-based 4D liver MRI reconstruction

Time-resolved volumetric magnetic resonance imaging (4D MRI) could be used to address organ motion in image-guided interventions like tumor ablation. Current 4D reconstruction techniques are unsuitable for most interventional settings because they are limited to specific breathing phases, lack temporal/spatial resolution, and have long prior acquisitions or reconstruction times. Deep learning-based (DL) 4D MRI approaches promise to overcome these shortcomings but are sensitive to domain shift. This work shows that transfer learning (TL) combined with an ensembling strategy can help alleviate this key challenge. We evaluate four approaches: pre-trained models from the source domain, models directly trained from scratch on target domain data, models fine-tuned from a pre-trained model and an ensemble of fine-tuned models. For that the data base was split into 16 source and 4 target domain subjects. Comparing ensemble of fine-tuned models (N = 10) with directly learned models, we report significant improvements (P < 0.001) of the root mean squared error (RMSE) of up to 12% and the mean displacement (MDISP) of up to 17.5%. The smaller the target domain data amount, the larger the effect. This shows that TL + Ens significantly reduces beforehand acquisition time and improves reconstruction quality, rendering it a key component in making 4D MRI clinically feasible for the first time in the context of 4D organ motion models of the liver and beyond.

Real-time motion monitoring using orthogonal cine MRI during MR-guided adaptive radiation therapy for abdominal tumors on 1.5T MR-Linac

BACKGROUND

Real-time motion monitoring (RTMM) is necessary for accurate motion management of intrafraction motions during radiation therapy (RT).

PURPOSE

Building upon a previous study, this work develops and tests an improved RTMM technique based on real-time orthogonal cine magnetic resonance imaging (MRI) acquired during magnetic resonance-guided adaptive RT (MRgART) for abdominal tumors on MR-Linac.

METHODS

A motion monitoring research package (MMRP) was developed and tested for RTMM based on template rigid registration between beam-on real-time orthogonal cine MRI and pre-beam daily reference 3D-MRI (baseline). The MRI data acquired under free-breathing during the routine MRgART on a 1.5T MR-Linac for 18 patients with abdominal malignancies of 8 liver, 4 adrenal glands (renal fossa), and 6 pancreas cases were used to evaluate the MMRP package. For each patient, a 3D mid-position image derived from an in-house daily 4D-MRI was used to define a target mask or a surrogate sub-region encompassing the target. Additionally, an exploratory case reviewed for an MRI dataset of a healthy volunteer acquired under both free-breathing and deep inspiration breath-hold (DIBH) was used to test how effectively the RTMM using the MMRP can address through-plane motion (TPM). For all cases, the 2D T2/T1-weighted cine MRIs were captured with a temporal resolution of 200 ms interleaved between coronal and sagittal orientations. Manually delineated contours on the cine frames were used as the ground-truth motion. Common visible vessels and segments of target boundaries in proximity to the target were used as anatomical landmarks for reproducible delineations on both the 3D and the cine MRI images. Standard deviation of the error (SDE) between the ground-truth and the measured target motion from the MMRP package were analyzed to evaluate the RTMM accuracy. The maximum target motion (MTM) was measured on the 4D-MRI for all cases during free-breathing.

RESULTS

The mean (range) centroid motions for the 13 abdominal tumor cases were 7.69 (4.71-11.15), 1.73 (0.81-3.05), and 2.71 (1.45-3.93) mm with an overall accuracy of <2 mm in the superior-inferior (SI), the left-right (LR), and the anterior-posterior (AP) directions, respectively. The mean (range) of the MTM from the 4D-MRI was 7.38 (2-11) mm in the SI direction, smaller than the monitored motion of centroid, demonstrating the importance of the real-time motion capture. For the remaining patient cases, the ground-truth delineation was challenging under free-breathing due to the target deformation and the large TPM in the AP direction, the implant-induced image artifacts, and/or the suboptimal image plane selection. These cases were evaluated based on visual assessment. For the healthy volunteer, the TPM of the target was significant under free-breathing which degraded the RTMM accuracy. RTMM accuracy of <2 mm was achieved under DIBH, indicating DIBH is an effective method to address large TPM.

CONCLUSIONS

We have successfully developed and tested the use of a template-based registration method for an accurate RTMM of abdominal targets during MRgART on a 1.5T MR-Linac without using injected contrast agents or radio-opaque implants. DIBH may be used to effectively reduce or eliminate TPM of abdominal targets during RTMM.

Predicting 4D liver MRI for MR-guided interventions

Organ motion poses an unresolved challenge in image-guided interventions like radiation therapy, biopsies or tumor ablation. In the pursuit of solving this problem, the research field of time-resolved volumetric magnetic resonance imaging (4D MRI) has evolved. However, current techniques are unsuitable for most interventional settings because they lack sufficient temporal and/or spatial resolution or have long acquisition times. In this work, we propose a novel approach for real-time, high-resolution 4D MRI with large fields of view for MR-guided interventions. To this end, we propose a network-agnostic, end-to-end trainable, deep learning formulation that enables the prediction of a 4D liver MRI with respiratory states from a live 2D navigator MRI. Our method can be used in two ways: First, it can reconstruct high quality fast (near real-time) 4D MRI with high resolution (209×128×128 matrix size with isotropic 1.8mm voxel size and 0.6s/volume) given a dynamic interventional 2D navigator slice for guidance during an intervention. Second, it can be used for retrospective 4D reconstruction with a temporal resolution of below 0.2s/volume for motion analysis and use in radiation therapy. We report a mean target registration error (TRE) of 1.19±0.74mm, which is below voxel size. We compare our results with a state-of-the-art retrospective 4D MRI reconstruction. Visual evaluation shows comparable quality. We compare different network architectures within our formulation. We show that small training sizes with short acquisition times down to 2 min can already achieve promising results and 24 min are sufficient for high quality results. Because our method can be readily combined with earlier time reducing methods, acquisition time can be further decreased while also limiting quality loss. We show that an end-to-end, deep learning formulation is highly promising for 4D MRI reconstruction.

Contrast-enhanced 4D-MRI for internal target volume generation in treatment planning for liver tumors

BACKGROUND

Liver tumors are often invisible on four-dimensional commuted tomography (4D-CT). Imperfect imaging surrogates are used to estimate the tumor motion. Here, we assessed multiple 4D magnetic resonance (MR) binning algorithms for directly visualizing liver tumor motion for radiotherapy planning.

METHODS

Patients were simulated using a 3 Tesla MR and CT scanner. Three prototype binning algorithms (phase, amplitude, and two-directional) were applied to the 4D-MRIs, and the image quality was assessed using a qualitative clarity score and quantitative sharpness score. Radiation plans were generated for internal target volumes (ITVs) derived using 4D-MRI and 4D-CT, and the dosimetry of targets were compared. Paired t-tests were used to compare sharpness scores and dosimetric data.

RESULTS

Twelve patients with 17 liver tumors were scanned between May and November 2021. Compared to phase binning, two-directional demonstrated equal or better clarity and sharpness scores (end-expiration: 0.33 vs. 0.38, p=0.018, end-inspiration: 0.28 vs. 0.31, p=0.010). Compared to amplitude binning, two-directional binning captured hysteresis of ≥3 mm in 35% of patients. Evaluation of dosimetry CT-optimized plans revealed that PTV coverage of MR-derived targets were significantly lower than CT-derived targets (PTV receiving 90% of prescription: 75.56% vs. 89.38%, p=0.002).

CONCLUSION

Using contrast-enhanced 4D-MRI is feasible for directly delineating liver tumors throughout the respiratory cycle. The current standard of using radiation plans optimized for 4D-CT-derived targets achieved lower coverage of directly visualized MRI targets, suggesting that adopting MRI for motion management may improve radiation treatment of liver lesions and reduce the risk of marginal misses.

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.

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.

First-in-human imaging using a MR-compatible e4D ultrasound probe for motion management of radiotherapy

PURPOSE

Respiration-induced tumor or organ positional changes can impact the accuracy of external beam radiotherapy. Motion management strategies are used to account for these changes during treatment. The authors report on the development, testing, and first-in-human evaluation of an electronic 4D (e4D) MR-compatible ultrasound probe that was designed for hands-free operation in a MR and linear accelerator (LINAC) environment.

METHODS

Ultrasound components were evaluated for MR compatibility. Electromagnetic interference (EMI) shielding was used to enclose the entire probe and a factory-fabricated cable shielded with copper braids was integrated into the probe. A series of simultaneous ultrasound and MR scans were acquired and analyzed in five healthy volunteers.

RESULTS

The ultrasound probe led to minor susceptibility artifacts in the MR images immediately proximal to the ultrasound probe at a depth of <10 mm. Ultrasound and MR-based motion traces that were derived by tracking the salient motion of endogenous target structures in the superior-inferior (SI) direction demonstrated good concordance (Pearson correlation coefficients of 0.95-0.98) between the ultrasound and MRI datasets.

CONCLUSION

We have demonstrated that our hands-free, e4D probe can acquire ultrasound images during a MR acquisition at frame rates of approximately 4 frames per second (fps) without impacting either the MR or ultrasound image quality. This use of this technology for interventional procedures (e.g. biopsies and drug delivery) and motion compensation during imaging are also being explored.

Slice-stacking T2-weighted MRI for fast determination of internal target volume for liver tumor

Background

To investigate the feasibility of generating maximum intensity projection (MIP) images to determine internal target volume (ITV) using slice-stacking MRI (SS-MRI) technique.

Methods

Slice-stacking is a technique which applies a multi-slice MRI acquisition to generate a 3D MIP for ITV contouring, without reconstructing 4D-MRI. 4D digital extended cardiac-torso (XCAT) phantom was used to generate MIP images with sequential 2D HASTE sequence, with different tumor diameters (10, 30 and 50 mm) and with simulated regular and irregular (patient) breathing motions. A reference MIP was generated using all acquisition images. Consecutive repetitions were then used to generate MIP to analyze the relationship between Dice's similarity coefficient (DSC) and the number of repetitions, and the relationship between the relative ITV volume difference and the number of repetitions. Images from XCAT phantom and from three hepatic carcinoma patients were collected in this study to demonstrate the feasibility of this technique.

Results

For both regular and irregular breathing motion, the average DSC of ITV is >0.94 and the average relative ITV volume difference is <10% (approximately 0.15 cm³) when using 5 repeated scanning images to reconstruct MIP for tumor diameter of 10 mm. As tumor diameter increases, the DSC of ITV is >0.97 and the relative ITV volume difference is <5% for regular breathing motion, and the DSC of ITV is >0.97 and the relative ITV volume difference is <5.5% for irregular breathing motion when using 5 repeated scanning images to reconstruct MIP. In patient image study, the mean relative ITV volume difference is <3% and the mean DSC is 0.99 when using 5 repeated scanning images to reconstruct MIP.

Conclusions

The number of scans required to generate tumor ITV for slice-stacking method (5-7 repetition) is 3-4 times less than that of 4D-MRI (15-20 repetitions). It is feasible to generate a fast clinically acceptable ITV using slice-stacking method with sequential 2D MR images.

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.

Free-breathing liver fat and R 2 ∗ quantification using motion-corrected averaging based on a nonlocal means algorithm

PURPOSE

To propose a motion-robust chemical shift-encoded (CSE) method with high signal-to-noise (SNR) for accurate quantification of liver proton density fat fraction (PDFF) and R 2 ∗ .

METHODS

A free-breathing multi-repetition 2D CSE acquisition with motion-corrected averaging using nonlocal means (NLM) was proposed. PDFF and R 2 ∗ quantified with 2D CSE-NLM were compared to two alternative 2D techniques: direct averaging and single acquisition (2D 1ave) in a digital phantom. Further, 2D NLM was compared in patients to 3D techniques (standard breath-hold, free-breathing and navigated), and the alternative 2D techniques. A reader study and quantitative analysis (Bland-Altman, correlation analysis, paired Student's t-test) were performed to evaluate the image quality and assess PDFF and R 2 ∗ measurements in regions of interest.

RESULTS

In simulations, 2D NLM resulted in lower standard deviations (STDs) of PDFF (2.7%) and R 2 ∗ (8.2  s - 1 ) compared to direct averaging (PDFF: 3.1%, R 2 ∗ : 13.6  s - 1 ) and 2D 1ave (PDFF: 8.7%, R 2 ∗ : 33.2  s - 1 ). In patients, 2D NLM resulted in fewer motion artifacts than 3D free-breathing and 3D navigated, less signal loss than 2D direct averaging, and higher SNR than 2D 1ave. Quantitatively, the STDs of PDFF and R 2 ∗ of 2D NLM were comparable to those of 2D direct averaging (p>0.05). 2D NLM reduced bias, particularly in R 2 ∗ (-5.73 to -0.36  s - 1 ) that arises in direct averaging (-3.96 to 11.22  s - 1 ) in the presence of motion.

CONCLUSIONS

2D CSE-NLM enables accurate mapping of PDFF and R 2 ∗ in the liver during free-breathing.

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.

The noise navigator: a surrogate for respiratory-correlated 4D-MRI for motion characterization in radiotherapy

Respiratory-correlated 4D-MRI can characterize respiratory-induced motion of tumors and organs-at-risk for radiotherapy treatment planning and is a necessity for image guidance of moving tumors treated on an MRI-linac. Essential for 4D-MRI generation is a robust respiratory surrogate signal. We investigated the feasibility of the noise navigator as respiratory surrogate signal for 4D-MRI generation. The noise navigator is based on the respiratory-induced modulation of the thermal noise variance measured by the receive coils during MR acquisition and thus is inherently present and synchronized with MRI data acquisition. Additionally, the noise navigator can be combined with any rectilinear readout strategy (e.g. radial and cartesian) and is independent of MR image contrast and imaging orientation. For radiotherapy applications, the noise navigator provides a robust respiratory signal for patients scanned with an elevated coil setup. This is particularly attractive for widely used cartesian sequences where currently a non-interfering self-navigation means is lacking for MRI-based simulation and MRI-guided radiotherapy. The feasibility of 4D-MRI generation with the noise navigator as respiratory surrogate signal was demonstrated for both cartesian and radial readout strategies in radiotherapy setup on four healthy volunteers and two radiotherapy patients on a dedicated 1.5 T MRI scanner and two radiotherapy patients on a 1.5 T MRI-linac system. Moreover, the respiratory-correlated 4D-MR images showed liver motion comparable to a reference 2D cine MRI series for the volunteers. For 2D cartesian cine MRI acquisitions, both the noise navigator and respiratory bellows were benchmarked against an image navigator. Respiratory phase detection based on the noise navigator agreed 1.4 times better with the image navigator than the respiratory bellows did. For a 3D Stack-of-Stars acquisitions, the noise navigator was compared to radial self-navigation and a 1.7 times higher respiratory phase detection agreement was observed than for the respiratory bellows compared to radial self-navigation.

Optimizing 4D abdominal MRI: image denoising using an iterative back-projection approach

4D-MRI is a promising tool for organ exploration, target delineation and treatment planning. Intra-scan motion artifacts may be greatly reduced by increasing the imaging frame rate. However, poor signal-to-noise ratios (SNR) are observed when increasing spatial and/or frame number per physiological cycle, in particular in the abdomen. In the current work, the proposed 4D-MRI method favored spatial resolution, frame number, isotropic voxels and large field-of-view (FOV) during MR-acquisition. The consequential SNR penalty in the reconstructed data is addressed retrospectively using an iterative back-projection (IBP) algorithm. Practically, after computing individual spatial 3D deformations present in the images using a deformable image registration (DIR) algorithm, each 3D image is individually enhanced by fusing several successive frames in its local temporal neighborood, these latter being likely to cover common independent informations. A tuning parameter allows one to freely readjust the balance between temporal resolution and precision of the 4D-MRI. The benefit of the method was quantitatively evaluated on the thorax of 6 mice under free breathing using a clinically acceptable duration. Improved 4D cardiac imaging was also shown in the heart of 1 mice. Obtained results are compared to theoretical expectations and discussed. The proposed implementation is easily parallelizable and optimized 4D-MRI could thereby be obtained with a clinically acceptable duration.

Effects of age-related breathing characteristics on the performance of four-dimensional magnetic resonance imaging reconstructed by prospective gating for radiation therapy planning

Background and purpose

Four-dimensional magnetic resonance imaging (4D MRI) has advanced recently by incorporating prospective gating, but its performance on pediatric populations has not been investigated. This study aimed to determine the age-related performance of prospective gating, as compared with retrospective sorting.

Materials and methods

Prospectively gated 4D MRI scans were acquired on a motion phantom driven by real respiratory waveforms obtained from 23 pediatric and young adult patients (aged 5–24 years). The correlations between patient-specific breathing characteristics and the performance of 4D MRI were comparatively evaluated against retrospective sorting for the same scan time. For six patients who underwent both 4D MRI and 4D CT, the internal target volumes (ITVs) determined by the two modalities were compared.

Results

Longer scan time and greater sorting error were most highly correlated (P < 0.001) with breathing irregularity and extent of diaphragm motion, but age was not a strong covariate because of interindividual variation. Prospective gating was more accurate than retrospective sorting except for those patients with severe breathing irregularity (peak-to-peak coefficient of variation >30%). The ITVs of 4D MRI and 4D CT were comparable (Dice similarity: >90%) unless the breathing characteristics changed between the two imaging sessions.

Conclusions

For most patients analyzed in this study, prospective gating provided more accurate 4D MRI (95th percentile of deviation: <1.5 mm) than did retrospective sorting within a clinically feasible scan time (median: 5.9 min). The 4D MRI tended to take longer and to give larger sorting errors with deeper and irregular breathers.

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.

A post-processing method based on interphase motion correction and averaging to improve image quality of 4D magnetic resonance imaging: a clinical feasibility study

METHODS:

Nine patients (seven pancreas, one liver, and one lung) were recruited. 4D-MRI was performed using two prototype k-space sorted techniques, stack-of-stars (SOS) and koosh-ball (KB) acquisitions. Post-processing using MoCoAve was implemented for both methods. Image quality score, apparent SNR (aSNR), sharpness, motion trajectory and standard deviation (σ_GTV) of the gross tumor volumes were compared between original and MoCoAve image sets.

RESULTS:

All subjects successfully underwent 4D-MRI scans and MoCoAve was performed on all data sets. Significantly higher image quality scores (2.64 ± 0.39 vs 1.18 ± 0.34, p = 0.001) and aSNR (37.6 ± 15.3 vs 18.1 ± 5.7, p = 0.001) was observed in the MoCoAve images when compared to the original images. High correlation in tumor motion trajectories in the superoinferior direction (SI: 0.91 ± 0.08) and weaker in the anteroposterior (AP: 0.51 ± 0.44) and mediolateral (ML: 0.37 ± 0.23) directions, similar image sharpness (0.367 ± 0.068 vs 0.369 ± 0.072, p = 0.805), and minimal average absolute difference (0.47 ± 0.34  mm) of the motion trajectory profiles was found between the two image sets. The σ_GTV in pancreas patients was significantly (p = 0.039) lower in MoCoAve images (1.48 ± 1.35  cm3) than in the original images (2.17 ± 1.31  cm3).

CONCLUSION:

MoCoAve using interphase motion correction and averaging has shown promise as a post-processing method for improving k-space sorted (SOS and KB) 4D-MRI image quality in thoracic and abdominal cancer patients.

ADVANCES IN KNOWLEDGE:

The proposed method is an image based post-processing method that could be applied to many k-space sorted 4D-MRI methods for improved image quality and signal-to-noise ratio while preserving image sharpness and respiratory motion fidelity. It is a useful technique for the radiotherapy planning community who are interested in using 4D-MRI but aren't satisfied with their current MR image quality.

A clustering approach to 4D MRI retrospective sorting for the investigation of different surrogates

PURPOSE

In retrospective 4-Dimensional Magnetic Resonance Imaging (4D MRI) sorting, respiratory surrogate selection affects the image quality of reconstructed volumes. We propose a method for retrospective 4D MRI sorting based on clustering, which allowed us to compare the performance of single or multiple internal surrogates vs. a conventional external signal.

METHODS

A k-medoids clustering algorithm was exploited for sorting 2D MRI into 4D MRI, relying on (A) multiple or (B) single automatically tracked internal landmarks or (C) respiratory belt signal. 4D MRI reconstructions for seven liver cancer patients were compared to those of the state-of-the-art mutual information (MI) approach. Sorting artifacts were measured by the root mean square error (RMSE) between the diaphragm profile and a fitted second order curve. Diaphragm and tumor motions were evaluated.

RESULTS

The median RMSEs ranged 0.97-1.66 mm, 1.24-1.89 mm, 1.43-2.27 mm, 1.74-3.72 mm for the MI, (A), (B) and (C) methods, respectively. Significant differences (Friedman, α = 5%) were found between (C) and all other methods, and between (B) and MI approaches. The discrepancies between (A) and MI approaches ranged 1.1-6.2 mm and 0.7-5.3 mm respectively in diaphragm and tumor motions. Methods (A) and (B) showed similar ranges of motion.

CONCLUSION

With multiple internal points, our method yielded the description of a higher range of motion and similar image quality with respect to the MI approach. The single point method led to more artifacts, suggesting the superior suitability of multiple internal surrogates for retrospective 4D MRI sorting. Considering internal rather than external information favored superior performance.

An image interpolation approach for acquisition time reduction in navigator-based 4D MRI

Navigated 2D multi-slice dynamic Magnetic Resonance (MR) imaging enables high contrast 4D MR imaging during free breathing and provides in-vivo observations for treatment planning and guidance. Navigator slices are vital for retrospective stacking of 2D data slices in this method. However, they also prolong the acquisition sessions. Temporal interpolation of navigator slices can be used to reduce the number of navigator acquisitions without degrading specificity in stacking. In this work, we propose a convolutional neural network (CNN) based method for temporal interpolation, with motion field prediction as an intermediate step. The proposed formulation incorporates the prior knowledge that a motion field underlies changes in the image intensities over time. Previous approaches that interpolate directly in the intensity space are prone to produce blurry images or even remove structures in the images. Our method avoids such problems and faithfully preserves the information in the image. Further, an important advantage of our formulation is that it provides an unsupervised estimation of bi-directional motion fields. These motion fields can potentially be used to halve the number of registrations required during 4D reconstruction, thus substantially reducing the reconstruction time. These advantages are achieved while preserving 4D reconstruction quality as compared to that with the true navigators.

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.

MRI-guidance for motion management in external beam radiotherapy: current status and future challenges

High precision conformal radiotherapy requires sophisticated imaging techniques to aid in target localisation for planning and treatment, particularly when organ motion due to respiration is involved. X-ray based imaging is a well-established standard for radiotherapy treatments. Over the last few years, the ability of magnetic resonance imaging (MRI) to provide radiation-free images with high-resolution and superb soft tissue contrast has highlighted the potential of this imaging modality for radiotherapy treatment planning and motion management. In addition, these advantageous properties motivated several recent developments towards combined MRI radiation therapy treatment units, enabling in-room MRI-guidance and treatment adaptation. The aim of this review is to provide an overview of the state-of-the-art in MRI-based image guidance for organ motion management in external beam radiotherapy. Methodological aspects of MRI for organ motion management are reviewed and their application in treatment planning, in-room guidance and adaptive radiotherapy described. Finally, a roadmap for an optimal use of MRI-guidance is highlighted and future challenges are discussed.

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.

A Spatiotemporal-Constrained Sorting Method for Motion-Robust 4D-MRI: A Feasibility Study

PURPOSE

To develop a spatiotemporal-constrained sorting technique for motion-robust 4 dimensional-magnetic resonance imaging.

METHODS AND MATERIALS

This sorting method implemented 2 new approaches for 4-dimensional imaging: (1) an optimized sparse k-space acquisition trajectory with self-gating signal derivation, and (2) a retrospective k-space sorting for reconstruction using a novel spatiotemporal-constrained strategy to minimize breathing variation-induced motion artifacts. Such sorting was regularized by a spatiotemporal index. Volumetric reconstruction was implemented iteratively with a secnd-order total generalized variation penalty. The proposed method was evaluated and compared with the conventional phase-sorting and amplitude-sorting methods in 2 studies. In a computer simulation study, 6 abdominal motion scenarios, including 2 cosine and 4 patient breathing motion patterns, were studied. Reconstruction accuracy was evaluated quantitatively in reference to the ground truth by average image relative error (IRE) in 10 phases and target Dice similarity coefficients (DSCs) in end-of-exhalation/inhalation phases. In addition, the proposed method was evaluated using a custom-made motion phantom. Reconstruction accuracy was evaluated by motion range measurement and image quality comparison in both fast and slow breathing motions.

RESULTS

In the simulation study, stitching motion artifacts in restricted images were lessened using the proposed method compared with those using the conventional methods. The average IRE and target DSC (end-of-exhalation/inhalation) were 0.031 and 0.95/0.94, respectively, suggesting better motion reconstruction accuracy than the phase-sorted method (IRE, 0.057; DSC, 0.89/0.89) and the amplitude-sorted method (IRE, 0.048; DSC, 0.91/0.88). In the phantom study, the moving target reconstructed by the proposed method demonstrated better rendering with less edge blurring. With fast breathing motion, the range measured using the proposed method was more accurate than that of the phase-sorted method and was comparable to the result of amplitude-sorted method and ground truths.

CONCLUSIONS

Preliminary results suggested that the proposed sorting technique could reconstruct high-quality images and accurate motion estimation with reduced artifacts in 4 dimensional-magnetic resonance imaging.

Comparison of four dimensional computed tomography and magnetic resonance imaging in abdominal radiotherapy planning

Background and Purpose

Four-dimensional (4D) computed tomography (CT) is widely used in radiotherapy (RT) planning and remains the current standard for motion evaluation. We assess a 4D magnetic resonance imaging (MRI) sequence in terms of motion and image quality in a phantom, healthy volunteers and patients undergoing RT.

Materials and Methods

The 4D-MRI sequence is a prototype T1-weighted 3D gradient echo with radial acquisition with self-gating. The accuracy of the 4D-MRI respiratory sorting based method was assessed using a MRI-CT compatible respiratory simulation phantom. In volunteers, abdominal viscera were evaluated for artefact, noise, structure delineation and overall image quality using a previously published four-point scoring system. In patients undergoing abdominal RT, the tumour (or a surrogate) was utilized to assess the range of motion on both 4D-CT and 4D-MRI. Furthermore, imaging quality was evaluated for both 4D-CT and 4D-MRI.

Results

In phantom studies 4D-MRI demonstrated amplitude of motion error of less than 0.2 mm for five, seven and ten bins. 4D-MRI provided excellent image quality for liver, kidney and pancreas. In patients, the median amplitude of motion seen on 4D-CT and 4D-MRI was 11.2 mm (range 2.8–20.3 mm) and 10.1 mm (range 0.7–20.7 mm) respectively. The median difference in amplitude between 4D-CT and 4D-MRI was −0.6 mm (range −3.4–5.2 mm). 4D-MRI demonstrated superior edge detection (median score 3 versus 1) and overall image quality (median score 2 versus 1) compared to 4D-CT.

Conclusions

The prototype 4D-MRI sequence demonstrated promising results and may be used in abdominal targeting, motion gating, and towards implementing MRI-based adaptive RT.

Reducing Navigators in Free-Breathing Abdominal MRI via Temporal Interpolation Using Convolutional Neural Networks

Navigated 2-D multi-slice dynamic magnetic resonance imaging (MRI) acquisitions are essential for MR guided therapies. This technique yields time-resolved volumetric images during free-breathing, which are ideal for visualizing and quantifying breathing induced motion. To achieve this, navigated dynamic imaging requires acquiring multiple navigator slices. Reducing the number of navigator slices would allow for acquiring more data slices in the same time, and hence, increasing through-plane resolution or alternatively the overall acquisition time can be reduced while keeping resolution unchanged. To this end, we propose temporal interpolation of navigator slices using convolutional neural networks (CNNs). Our goal is to acquire fewer navigators and replace the missing ones with interpolation. We evaluate the proposed method on abdominal navigated dynamic MRI sequences acquired from 14 subjects. Investigations with several CNN architectures and training loss functions show favorable results for cost and a simple feed-forward network with no skip connections. When compared with interpolation by non-linear registration, the proposed method achieves higher interpolation accuracy on average as quantified in terms of root mean square error and residual motion. Analysis of the differences shows that the better performance is due to more accurate interpolation at peak exhalation and inhalation positions. Furthermore, the CNN-based approach requires substantially lower execution times than that of the registration-based method. At last, experiments on dynamic volume reconstruction reveal minimal differences between reconstructions with acquired and interpolated navigator slices.

Ultrasound-driven 4D MRI

We present an ultrasound-driven 4D magnetic resonance imaging (US-4DMRI) method for respiratory motion imaging in the thorax and abdomen. The proposed US-4DMRI comes along with a high temporal resolution, and allows for organ motion imaging beyond a single respiratory cycle. With the availability of the US surrogate both inside and outside the MR bore, 4D MR images can be reconstructed for 4D treatment planning and online respiratory motion prediction during radiotherapy. US-4DMRI relies on simultaneously acquired 2D liver US images and abdominal 2D MR multi-slice scans under free respiration. MR volumes are retrospectively composed by grouping the MR slices corresponding to the most similar US images. We present two different US similarity metrics: an intensity-based approach, and a similarity measure relying on predefined fiducials which are being tracked over time. The proposed method is demonstrated on MR liver scans of eight volunteers acquired over a duration of 5.5 min each at a temporal resolution of 2.6 Hz with synchronous US imaging at 14 Hz-17 Hz. Visual inspection of the reconstructed MR volumes revealed satisfactory results in terms of continuity in organ boundaries and blood vessels. In quantitative leave-one-out experiments, both US similarity metrics reach the performance level of state-of-the-art navigator-based approaches.

A Self-Sorting Coronal 4D-MRI Method for Daily Image Guidance of Liver Lesions on an MR-LINAC

PURPOSE

Novel hybrid MR-LINAC devices provide MRI's superior soft-tissue contrast in the treatment room and thus have the potential to increase accuracy of liver stereotactic body radiation therapy (SBRT). Requirements for daily position verification using 4-dimensional MRI include tumor visibility and short acquisition-reconstruction time (preferably <5 min). The proposed method provides fast acquisition-reconstruction time and the flexibility to vary T1- and T2-weighting, using standard imaging sequences for straightforward implementation on an MR-LINAC.

METHODS AND MATERIALS

Images were acquired using a coronal 2-dimensional, multislice, single-shot turbo spin-echo (TSE) and turbo field-echo (TFE) sequence, which were repeated 30 times. An image-based self-sorting signal (ImS) was extracted from the data, and rigid registration of the diaphragm per slice position was performed and corrected for amplitude variation in the anteroposterior direction. Data were sorted into 10 bins according to amplitude and phase. ImS was validated in 4 healthy volunteers against a navigator signal. Positional variations within bins, missing data, and smoothness of the liver dome were compared between amplitude and phase binning in 10 volunteers. Tumor contrast and registration were investigated in 3 patients.

RESULTS

Each ImS was found to be in excellent agreement with the navigator signal with a correlation coefficient of >0.95 and binning differences of <1 bin. Better liver dome smoothness per bin in case of amplitude binning compared with that in phase binning (2.0-2.6 mm vs 2.4-3.7 mm, respectively) is a tradeoff for more missing data (3.5%-17.5% vs 3.5%-4.7%, respectively). Liver lesions were visible in almost all coronal TSE and TFE images, but the lesion boundary was better defined in the TSE images. Rigid registrations could be performed on the tumor area.

CONCLUSIONS

An efficient self-sorted 4-dimensional MRI method was developed and validated using standard sequences and fast reconstruction on a LINAC-integrated MRI scanner providing good tumor visibility for daily image-guided liver stereotactic body radiation therapy.

An efficient model to guide prospective T2-weighted 4D magnetic resonance imaging acquisition

PURPOSE

To establish a mathematical model to guide prospective T2-weighted four-dimensional magnetic resonance imaging (4DMRI) acquisition and to propose an efficient solution to expedite prospective T2-weighted 4DMRI acquisition.

METHODS

Prospective T2-weighted 4DMRI acquisition was characterized by a mathematical model with 4DMRI acquisition time as the objective function and completeness of the image set, acquisition timing, image contrast, and image artifacts as constraints. Given the irregular nature of human respiration, an efficient solution based on the greedy strategy (ESGS) was proposed. The efficiency of the ESGS method was validated using healthy human subjects. Comparisons were made with the previous 4DMRI method incorporating the prefixed-order respiratory state splitting (PO-RSS) technique.

RESULTS

4DMRI image sets acquired using the ESGS and PO-RSS methods had similar image quality. The average time to acquire a 4DMRI image set covering 60 slices at 10 respiratory states was reduced by 30%, from 13.1 min using the PO-RSS method to 9.0 min using the ESGS method. It was demonstrated that high-quality T2-weighted 4DMRI could be obtained within a reasonable amount of time and all slices within each of the three-dimensional volumes were indeed acquired at the same respiratory state.

CONCLUSIONS

The ESGS method substantially reduces the acquisition time for T2-weighted 4DMRI, making it ready for clinical evaluation to obtain abdominal tumor motion for radiotherapy treatment planning.

Geometric accuracy of the MR imaging techniques in the presence of motion

Magnetic Resonance Imaging (MRI) is increasingly being used for improving tumor delineation and tumor tracking in the presence of respiratory motion. The purpose of this work is to design and build an MR compatible motion platform and to use it for evaluating the geometric accuracy of MR imaging techniques during respiratory motion. The motion platform presented in this work is composed of a mobile base made up of a flat plate and four wheels. The mobile base is attached from one end and through a rigid rod to a synchrony motion table by Accuray^® placed at the end of the MRI table and from the other end to an elastic rod. The geometric accuracy was measured by placing a control point‐based phantom on top of the mobile base.

In‐house software module was used to automatically assess the geometric distortion. The blurring artifact was also assessed by measuring the Full Width Half Maximum (FWHM) of each control point. Our results were assessed for 50, 100, and 150 mm radial distances, with a mean geometric distortion during the superior–inferior motion of 0.27, 0.41, and 0.55 mm, respectively. Adding the anterior–posterior motion, the mean geometric distortions increased to 0.4, 0.6, and 0.8 mm. Blurring was observed during motion causing an increase in the FWHM of ≈30%. The platform presented in this work provides a valuable tool for the assessment of the geometric accuracy and blurring artifact for MR during motion. Although the main objective was to test the spatial accuracy of an MR system during motion, the modular aspect of the presented platform enables the use of any commercially available phantom for a full quality control of the MR system during motion.

An MR-Based Model for Cardio-Respiratory Motion Compensation of Overlays in X-Ray Fluoroscopy

In X-ray fluoroscopy, static overlays are used to visualize soft tissue. We propose a system for cardiac and respiratory motion compensation of these overlays. It consists of a 3-D motion model created from real-time magnetic resonance (MR) imaging. Multiple sagittal slices are acquired and retrospectively stacked to consistent 3-D volumes. Slice stacking considers cardiac information derived from the ECG and respiratory information extracted from the images. Additionally, temporal smoothness of the stacking is enhanced. Motion is estimated from the MR volumes using deformable 3-D/3-D registration. The motion model itself is a linear direct correspondence model using the same surrogate signals as slice stacking. In X-ray fluoroscopy, only the surrogate signals need to be extracted to apply the motion model and animate the overlay in real time. For evaluation, points are manually annotated in oblique MR slices and in contrast-enhanced X-ray images. The 2-D Euclidean distance of these points is reduced from 3.85 to 2.75 mm in MR and from 3.0 to 1.8 mm in X-ray compared with the static baseline. Furthermore, the motion-compensated overlays are shown qualitatively as images and videos.

The development of a 4D treatment planning methodology to simulate the tracking of central lung tumors in an MRI-linac

PURPOSE

Targeting and tracking of central lung tumors may be feasible on the Elekta MRI-linac (MRL) due to the soft-tissue visualization capabilities of MRI. The purpose of this work is to develop a novel treatment planning methodology to simulate tracking of central lung tumors with the MRL and to quantify the benefits in OAR sparing compared with the ITV approach.

METHODS

Full 4D-CT datasets for five central lung cancer patients were selected to simulate the condition of having 4D-pseudo-CTs derived from 4D-MRI data available on the MRL with real-time tracking capabilities. We used the MRL treatment planning system to generate two plans: (a) with a set of MLC-defined apertures around the target at each phase of the breathing ("4D-MRL" method); (b) with a fixed set of fields encompassing the maximum inhale and exhale of the breathing cycle ("ITV" method). For both plans, dose accumulation was performed onto a reference phase. To further study the potential benefits of a 4D-MRL method, the results were stratified by tumor motion amplitude, OAR-to-tumor proximity, and the relative OAR motion (ROM).

RESULTS

With the 4D-MRL method, the reduction in mean doses was up to 3.0 Gy and 1.9 Gy for the heart and the lung. Moreover, the lung's V12.5 Gy was spared by a maximum of 300 cc. Maximum doses to serial organs were reduced by up to 6.1 Gy, 1.5 Gy, and 9.0 Gy for the esophagus, spinal cord, and the trachea, respectively. OAR dose reduction with our method depended on the tumor motion amplitude and the ROM. Some OARs with large ROMs and in close proximity to the tumor benefited from tracking despite small tumor amplitudes.

CONCLUSIONS

We developed a novel 4D tracking methodology for the MRL for central lung tumors and quantified the potential dosimetric benefits compared with our current ITV approach.

Accelerating volumetric cine MRI (VC-MRI) using undersampling for real-time 3D target localization/tracking in radiation therapy: a feasibility study

PURPOSE

To accelerate volumetric cine MRI (VC-MRI) using undersampled 2D-cine MRI to provide real-time 3D guidance for gating/target tracking in radiotherapy.

METHODS

4D-MRI is acquired during patient simulation. One phase of the prior 4D-MRI is selected as the prior images, designated as MRIprior. The on-board VC-MRI at each time-step is considered a deformation of the MRIprior. The deformation field map is represented as a linear combination of the motion components extracted by principal component analysis from the prior 4D-MRI. The weighting coefficients of the motion components are solved by matching the corresponding 2D-slice of the VC-MRI with the on-board undersampled 2D-cine MRI acquired. Undersampled Cartesian and radial k-space acquisition strategies were investigated. The effects of k-space sampling percentage (SP) and distribution, tumor sizes and noise on the VC-MRI estimation were studied. The VC-MRI estimation was evaluated using XCAT simulation of lung cancer patients and data from liver cancer patients. Volume percent difference (VPD) and Center of Mass Shift (COMS) of the tumor volumes and tumor tracking errors were calculated.

RESULTS

For XCAT, VPD/COMS were 11.93  ±  2.37%/0.90  ±  0.27 mm and 11.53  ±  1.47%/0.85  ±  0.20 mm among all scenarios with Cartesian sampling (SP  =  10%) and radial sampling (21 spokes, SP  =  5.2%), respectively. When tumor size decreased, higher sampling rate achieved more accurate VC-MRI than lower sampling rate. VC-MRI was robust against noise levels up to SNR  =  20. For patient data, the tumor tracking errors in superior-inferior, anterior-posterior and lateral (LAT) directions were 0.46  ±  0.20 mm, 0.56  ±  0.17 mm and 0.23  ±  0.16 mm, respectively, for Cartesian-based sampling with SP  =  20% and 0.60  ±  0.19 mm, 0.56  ±  0.22 mm and 0.42  ±  0.15 mm, respectively, for radial-based sampling with SP  =  8% (32 spokes).

CONCLUSIONS

It is feasible to estimate VC-MRI from a single undersampled on-board 2D cine MRI. Phantom and patient studies showed that the temporal resolution of VC-MRI can potentially be improved by 5-10 times using a 2D cine image acquired with 10-20% k-space sampling.

MRI-guided lung SBRT: Present and future developments

Stereotactic body radiotherapy (SBRT) is rapidly becoming an alternative to surgery for the treatment of early-stage non-small cell lung cancer patients. Lung SBRT is administered in a hypo-fractionated, conformal manner, delivering high doses to the target. To avoid normal-tissue toxicity, it is crucial to limit the exposure of nearby healthy organs-at-risk (OAR). Current image-guided radiotherapy strategies for lung SBRT are mostly based on X-ray imaging modalities. Although still in its infancy, magnetic resonance imaging (MRI) guidance for lung SBRT is not exposure-limited and MRI promises to improve crucial soft-tissue contrast. Looking beyond anatomical imaging, functional MRI is expected to inform treatment decisions and adaptations in the future. This review summarises and discusses how MRI could be advantageous to the different links of the radiotherapy treatment chain for lung SBRT: diagnosis and staging, tumour and OAR delineation, treatment planning, and inter- or intrafractional motion management. Special emphasis is placed on a new generation of hybrid MRI treatment devices and their potential for real-time adaptive radiotherapy.

Improved vessel-tissue contrast and image quality in 3D radial sampling-based 4D-MRI

Purpose

In radiation treatment planning for thoracic and abdominal tumors, 4D‐MRI has shown promise in respiratory motion characterization with improved soft‐tissue contrast compared to clinical standard, 4D computed tomography (4D‐CT). This study aimed to further improve vessel–tissue contrast and overall image quality in 3D radial sampling‐based 4D‐MRI using a slab‐selective (SS) excitation approach.

Methods

The technique was implemented in a 3D radial sampling with self‐gating‐based k‐space sorting sequence. The SS excitation approach was compared to a non‐selective (NS) approach in six cancer patients and two healthy volunteers at 3T. Improvements in vessel–tissue contrast ratio (CR) and vessel signal‐to‐noise ratio (SNR) were analyzed in five of the eight subjects. Image quality was visually assessed in all subjects on a 4‐point scale (0: poor; 3: excellent). Tumor (patients) and pancreas (healthy) motion trajectories were compared between the two imaging approaches.

Results

Compared with NS‐4D‐MRI, SS‐4D‐MRI significantly improved the overall vessel–tissue CR (2.60 ± 3.97 vs. 1.03 ± 1.44, P < 0.05), SNR (63.33 ± 38.45 vs. 35.74 ± 28.59, P < 0.05), and image quality score (2.6 ± 0.5 vs. 1.4 ± 0.5, P = 0.02). Motion trajectories from the two approaches exhibited strong correlation in the superior–inferior (0.96 ± 0.06), but weaker in the anterior–posterior (0.78 ± 0.24) and medial–lateral directions (0.46 ± 0.44).

Conclusions

The proposed 4D‐MRI with slab‐selectively excited 3D radial sampling allows for improved blood SNR, vessel–tissue CR, and image quality.

T2-Weighted 4D Magnetic Resonance Imaging for Application in Magnetic Resonance-Guided Radiotherapy Treatment Planning

OBJECTIVES

The aim of this study was to develop and verify a method to obtain good temporal resolution T2-weighted 4-dimensional (4D-T2w) magnetic resonance imaging (MRI) by using motion information from T1-weighted 4D (4D-T1w) MRI, to support treatment planning in MR-guided radiotherapy.

MATERIALS AND METHODS

Ten patients with primary non-small cell lung cancer were scanned at 1.5 T axially with a volumetric T2-weighted turbo spin echo sequence gated to exhalation and a volumetric T1-weighted stack-of-stars spoiled gradient echo sequence with golden angle spacing acquired in free breathing. From the latter, 20 respiratory phases were reconstructed using the recently developed 4D joint MoCo-HDTV algorithm based on the self-gating signal obtained from the k-space center. Motion vector fields describing the respiratory cycle were obtained by deformable image registration between the respiratory phases and projected onto the T2-weighted image volume. The resulting 4D-T2w volumes were verified against the 4D-T1w volumes: an edge-detection method was used to measure the diaphragm positions; the locations of anatomical landmarks delineated by a radiation oncologist were compared and normalized mutual information was calculated to evaluate volumetric image similarity.

RESULTS

High-resolution 4D-T2w MRI was obtained. Respiratory motion was preserved on calculated 4D-T2w MRI, with median diaphragm positions being consistent with less than 6.6 mm (2 voxels) for all patients and less than 3.3 mm (1 voxel) for 9 of 10 patients. Geometrical positions were coherent between 4D-T1w and 4D-T2w MRI as Euclidean distances between all corresponding anatomical landmarks agreed to within 7.6 mm (Euclidean distance of 2 voxels) and were below 3.8 mm (Euclidean distance of 1 voxel) for 355 of 470 pairs of anatomical landmarks. Volumetric image similarity was commensurate between 4D-T1w and 4D-T2w MRI, as mean percentage differences in normalized mutual information (calculated over all respiratory phases and patients), between corresponding respiratory phases of 4D-T1w and 4D-T2w MRI and the tie-phase of 4D-T1w and 3-dimensional T2w MRI, were consistent to 0.41% ± 0.37%. Four-dimensional T2w MRI displayed tumor extent, structure, and position more clearly than corresponding 4D-T1w MRI, especially when mobile tumor sites were adjacent to organs at risk.

CONCLUSIONS

A methodology to obtain 4D-T2w MRI that retrospectively applies the motion information from 4D-T1w MRI to 3-dimensional T2w MRI was developed and verified. Four-dimensional T2w MRI can assist clinicians in delineating mobile lesions that are difficult to define on 4D-T1w MRI, because of poor tumor-tissue contrast.

Retrospective four-dimensional magnetic resonance imaging with image-based respiratory surrogate: a sagittal-coronal-diaphragm point of intersection motion tracking method

A four-dimensional magnetic resonance imaging (4-D-MRI) technique with Sagittal-Coronal-Diaphragm Point-of-Intersection (SCD-PoI) as a respiratory surrogate is proposed. To develop an image-based respiratory surrogate, the SCD-PoI motion tracking method is used for retrospective 4-D-MRI reconstruction. Single-slice sagittal MR cine was acquired at a location near the center of the diaphragmatic dome. Multiple-slice coronal MR cines were acquired for 4-D-MRI reconstruction. As a motion surrogate, the diaphragm motion was measured from the PoI among the sagittal MRI cine plane, coronal MRI cine planes, and the diaphragm surface. These points were defined as the SCD-PoI. This point is used as a one-dimensional diaphragmatic navigator in our study. The 4-D-MRI technique was evaluated on a 4-D digital extended cardiac-torso (XCAT) human phantom, a motion phantom, and seven human subjects (five healthy volunteers and two cancer patients). Motion trajectories of a selected region of interest were measured on 4-D-MRI and compared with the known XCAT motion that served as references. The mean absolute amplitude difference ([Formula: see text]) and the cross-correlation coefficient (CC) of the comparisons were determined. 4-D-MRI of the XCAT phantom demonstrated highly accurate motion information ([Formula: see text], [Formula: see text]). Motion trajectories of the motion phantom measured on 4-D-MRI matched well with the references ([Formula: see text], [Formula: see text]). 4-D-MRI of human subjects showed minimal artifacts and clearly revealed the respiratory motion of organs and tumor (mean [Formula: see text]; mean [Formula: see text]). A 4-D-MRI technique with image-based respiratory surrogate has been developed and tested on phantoms and human subjects.