Volumetric coronary endothelial function assessment: a feasibility study exploiting stack-of-stars 3D cine MRI and image-based respiratory self-gating.

Abnormal coronary endothelial function (CEF), manifesting as depressed vasoreactive responses to endothelial-specific stressors, occurs early in atherosclerosis, independently predicts cardiovascular events, and responds to cardioprotective interventions. CEF is spatially heterogeneous along a coronary artery in patients with atherosclerosis, and thus recently developed and tested non-invasive 2D MRI techniques to measure CEF may not capture the extent of changes in CEF in a given coronary artery. The purpose of this study was to develop and test the first volumetric coronary 3D MRI cine method for assessing CEF along the proximal and mid-coronary arteries with isotropic spatial resolution and in free-breathing. This approach, called 3D-Stars, combines a 6 min continuous, untriggered golden-angle stack-of-stars acquisition with a novel image-based respiratory self-gating method and cardiac and respiratory motion-resolved reconstruction. The proposed respiratory self-gating method agreed well with respiratory bellows and center-of-k-space methods. In healthy subjects, 3D-Stars vessel sharpness was non-significantly different from that by conventional 2D radial in proximal segments, albeit lower in mid-portions. Importantly, 3D-Stars detected normal vasodilatation of the right coronary artery in response to endothelial-dependent isometric handgrip stress in healthy subjects. Coronary artery cross-sectional areas measured using 3D-Stars were similar to those from 2D radial MRI when similar thresholding was used. In conclusion, 3D-Stars offers good image quality and shows feasibility for non-invasively studying vasoreactivity-related lumen area changes along the proximal coronary artery in 3D during free-breathing.

Precision assessment of bowel motion quantification using 3D cine-MRI for radiotherapy.

Objective. The bowel is an important organ at risk for toxicity during pelvic and abdominal radiotherapy. Identifying regions of high and low bowel motion with MRI during radiotherapy may help to understand the development of bowel toxicity, but the acquisition time of MRI is rather long. The aim of this study is to retrospectively evaluate the precision of bowel motion quantification and to estimate the minimum MRI acquisition time.Approach. We included 22 gynaecologic cancer patients receiving definitive radiotherapy with curative intent. The 10 min pre-treatment 3D cine-MRI scan consisted of 160 dynamics with an acquisition time of 3.7 s per volume. Deformable registration of consecutive images generated 159 deformation vector fields (DVFs). We defined two motion metrics, the 50th percentile vector lengths (VL50) of the complete set of DVFs was used to measure median bowel motion. The 95th percentile vector lengths (VL95) was used to quantify high motion of the bowel. The precision of these metrics was assessed by calculating their variation (interquartile range) in three different time frames, defined as subsets of 40, 80, and 120 consecutive images, corresponding to acquisition times of 2.5, 5.0, and 7.5 min, respectively.Main results. For the full 10 min scan, the minimum motion per frame of 50% of the bowel volume (M50%) ranged from 0.6-3.5 mm for the VL50 motion metric and 2.3-9.0 mm for the VL95 motion metric, across all patients. At 7.5 min scan time, the variation in M50% was less than 0.5 mm in 100% (VL50) and 95% (VL95) of the subsets. A scan time of 5.0 and 2.5 min achieved a variation within 0.5 mm in 95.2%/81% and 85.7%/57.1% of the subsets, respectively.Significance. Our 3D cine-MRI technique quantifies bowel loop motion with 95%-100% confidence with a precision of 0.5 mm variation or less, using a 7.5 min scan time.

3D cine-magnetic resonance imaging using spatial and temporal implicit neural representation learning (STINR-MR).

Objective. 3D cine-magnetic resonance imaging (cine-MRI) can capture images of the human body volume with high spatial and temporal resolutions to study anatomical dynamics. However, the reconstruction of 3D cine-MRI is challenged by highly under-sampled k-space data in each dynamic (cine) frame, due to the slow speed of MR signal acquisition. We proposed a machine learning-based framework, spatial and temporal implicit neural representation learning (STINR-MR), for accurate 3D cine-MRI reconstruction from highly under-sampled data.Approach. STINR-MR used a joint reconstruction and deformable registration approach to achieve a high acceleration factor for cine volumetric imaging. It addressed the ill-posed spatiotemporal reconstruction problem by solving a reference-frame 3D MR image and a corresponding motion model that deforms the reference frame to each cine frame. The reference-frame 3D MR image was reconstructed as a spatial implicit neural representation (INR) network, which learns the mapping from input 3D spatial coordinates to corresponding MR values. The dynamic motion model was constructed via a temporal INR, as well as basis deformation vector fields (DVFs) extracted from prior/onboard 4D-MRIs using principal component analysis. The learned temporal INR encodes input time points and outputs corresponding weighting factors to combine the basis DVFs into time-resolved motion fields that represent cine-frame-specific dynamics. STINR-MR was evaluated using MR data simulated from the 4D extended cardiac-torso (XCAT) digital phantom, as well as two MR datasets acquired clinically from human subjects. Its reconstruction accuracy was also compared with that of the model-based non-rigid motion estimation method (MR-MOTUS) and a deep learning-based method (TEMPEST).Main results. STINR-MR can reconstruct 3D cine-MR images with high temporal (<100 ms) and spatial (3 mm) resolutions. Compared with MR-MOTUS and TEMPEST, STINR-MR consistently reconstructed images with better image quality and fewer artifacts and achieved superior tumor localization accuracy via the solved dynamic DVFs. For the XCAT study, STINR reconstructed the tumors to a mean ± SD center-of-mass error of 0.9 ± 0.4 mm, compared to 3.4 ± 1.0 mm of the MR-MOTUS method. The high-frame-rate reconstruction capability of STINR-MR allows different irregular motion patterns to be accurately captured.Significance. STINR-MR provides a lightweight and efficient framework for accurate 3D cine-MRI reconstruction. It is a 'one-shot' method that does not require external data for pre-training, allowing it to avoid generalizability issues typically encountered in deep learning-based methods.

Inter-scanner super-resolution of 3D cine MRI using a transfer-learning network for MRgRT.

Objective. Deep-learning networks for super-resolution (SR) reconstruction enhance the spatial-resolution of 3D magnetic resonance imaging (MRI) for MR-guided radiotherapy (MRgRT). However, variations between MRI scanners and patients impact the quality of SR for real-time 3D low-resolution (LR) cine MRI. In this study, we present a personalized super-resolution (psSR) network that incorporates transfer-learning to overcome the challenges in inter-scanner SR of 3D cine MRI.Approach: Development of the proposed psSR network comprises two-stages: (1) a cohort-specific SR (csSR) network using clinical patient datasets, and (2) a psSR network using transfer-learning to target datasets. The csSR network was developed by training on breath-hold and respiratory-gated high-resolution (HR) 3D MRIs and their k-space down-sampled LR MRIs from 53 thoracoabdominal patients scanned at 1.5 T. The psSR network was developed through transfer-learning to retrain the csSR network using a single breath-hold HR MRI and a corresponding 3D cine MRI from 5 healthy volunteers scanned at 0.55 T. Image quality was evaluated using the peak-signal-noise-ratio (PSNR) and the structure-similarity-index-measure (SSIM). The clinical feasibility was assessed by liver contouring on the psSR MRI using an auto-segmentation network and quantified using the dice-similarity-coefficient (DSC).Results. Mean PSNR and SSIM values of psSR MRIs were increased by 57.2% (13.8-21.7) and 94.7% (0.38-0.74) compared to cine MRIs, with the reference 0.55 T breath-hold HR MRI. In the contour evaluation, DSC was increased by 15% (0.79-0.91). Average time consumed for transfer-learning was 90 s, psSR was 4.51 ms per volume, and auto-segmentation was 210 ms, respectively.Significance. The proposed psSR reconstruction substantially increased image and segmentation quality of cine MRI in an average of 215 ms across the scanners and patients with less than 2 min of prerequisite transfer-learning. This approach would be effective in overcoming cohort- and scanner-dependency of deep-learning for MRgRT.

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) D99% of the GTV was 0.94 (0.06), and the mean (SD) relative D0.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.

3D cine-magnetic resonance imaging using spatial and temporal implicit neural representation learning (STINR-MR).

Objective: 3D cine-magnetic resonance imaging (cine-MRI) can capture images of the human body volume with high spatial and temporal resolutions to study the anatomical dynamics. However, the reconstruction of 3D cine-MRI is challenged by highly undersampled k-space data in each dynamic (cine) frame, due to the slow speed of MR signal acquisition. We proposed a machine learning-based framework, spatial and temporal implicit neural representation learning (STINR-MR), for accurate 3D cine-MRI reconstruction from highly undersampled data. Approach: STINR-MR used a joint reconstruction and deformable registration approach to achieve a high acceleration factor for cine volumetric imaging. It addressed the ill-posed spatiotemporal reconstruction problem by solving a reference-frame 3D MR image and a corresponding motion model which deforms the reference frame to each cine frame. The reference-frame 3D MR image was reconstructed as a spatial implicit neural representation (INR) network, which learns the mapping from input 3D spatial coordinates to corresponding MR values. The dynamic motion model was constructed via a temporal INR, as well as basis deformation vector fields (DVFs) extracted from prior/onboard 4D-MRIs using principal component analysis (PCA). The learned temporal INR encodes input time points and outputs corresponding weighting factors to combine the basis DVFs into time-resolved motion fields that represent cine-frame-specific dynamics. STINR-MR was evaluated using MR data simulated from the 4D extended cardiac-torso (XCAT) digital phantom, as well as MR data acquired clinically from a healthy human subject. Its reconstruction accuracy was also compared with that of the model-based non-rigid motion estimation method (MR-MOTUS). Main results: STINR-MR can reconstruct 3D cine-MR images with high temporal (<100 ms) and spatial (3 mm) resolutions. Compared with MR-MOTUS, STINR-MR consistently reconstructed images with better image quality and fewer artifacts and achieved superior tumor localization accuracy via the solved dynamic DVFs. For the XCAT study, STINR reconstructed the tumors to a mean±S.D. center-of-mass error of 1.0±0.4 mm, compared to 3.4±1.0 mm of the MR-MOTUS method. The high-frame-rate reconstruction capability of STINR-MR allows different irregular motion patterns to be accurately captured. Significance: STINR-MR provides a lightweight and efficient framework for accurate 3D cine-MRI reconstruction. It is a 'one-shot' method that does not require external data for pre-training, allowing it to avoid generalizability issues typically encountered in deep learning-based methods.

A 3D cine-MRI acquisition technique and image analysis framework to quantify bowel motion demonstrated in gynecological cancer patients.

PURPOSE: Magnetic resonance imaging (MRI) is increasingly used in radiation oncology for target delineation and radiotherapy treatment planning, for example, in patients with gynecological cancers. As a consequence of pelvic radiotherapy, a part of the bowel is irradiated, yielding risk of bowel toxicity. Existing dose-effect models predicting bowel toxicity are inconclusive and bowel motion might be an important confounding factor. The exact motion of the bowel and dosimetric effects of its motion are yet uncharted territories in radiotherapy. In diagnostic radiology methods on the acquisition of dynamic MRI sequences were developed for bowel motility visualization and quantification. Our study aim was to develop an imaging technique based on three-dimensional (3D) cine-MRI to visualize and quantify bowel motion and demonstrate it in a cohort of gynecological cancer patients. METHODS: We developed an MRI acquisition suitable for 3D bowel motion quantification, namely a balanced turbo field echo sequence (TE = 1.39 ms, TR = 2.8 ms), acquiring images in 3.7 s (dynamic) with a 1.25 × 1.25 × 2.5 mm3 resolution, yielding a field of view of 200 × 200 × 125 mm3 . These MRI bowel motion sequences were acquired in 22 gynecological patients. During a 10-min scan, 160 dynamics were acquired. Subsequent dynamics were deformably registered using a B-spline transformation model, resulting in 159 3D deformation vector fields (DVFs) per MRI set. From the 159 DVFs, the average vector length was calculated per voxel to generate bowel motion maps. Quality assurance was performed on all 159 DVFs per MRI, using the Jacobian Determinant and the Harmonic Energy as deformable image registration error metrics. In order to quantify bowel motion, we introduced the concept of cumulative motion-volume histogram (MVH) of the bowel bag volume. Finally, interpatient variation of bowel motion was analyzed using the MVH parameters M10%, M50%, and M90%. The M10%/M50%/M90% represents the minimum bowel motion per frame of 10%/50%/90% of the bowel bag volume. RESULTS: The motion maps resulted in a visualization of areas with small and large movements within the bowel bag. After applying quality assurance, the M10%, M50%, and M90% were 4.4 (range 2.2-7.6) mm, 2.2 (range 0.9-4.1) mm, and 0.5 (range 0.2-1.4) mm per frame, on average over all patients, respectively. CONCLUSION: We have developed a method to visualize and quantify 3D bowel motion with the use of bowel motion specific MRI sequences in 22 gynecological cancer patients. This 3D cine-MRI-based quantification tool and the concept of MVHs can be used in further studies to determine the effect of radiotherapy on bowel motion and to find the relation with dose effects to the small bowel. In addition, the developed technique can be a very interesting application for bowel motility assessment in diagnostic radiology.