An optimal radial profile order based on the Golden Ratio for time-resolved MRI

Real-time water/fat imaging at 0.55T with spiral out-in-out-in sampling

PURPOSE

To develop an efficient and flexible water/fat separated real-time MRI (RT-MRI) method using spiral out-in-out-in (OIOI) sampling and balanced SSFP (bSSFP) at 0.55T.

METHODS

A bSSFP sequence with golden-angle spiral OIOI readout was developed, capturing three echoes to allow water/fat separation. A low-latency reconstruction that combines all echoes was available for online visualization. An offline reconstruction provided water and fat RT-MRI in two steps: (1) image reconstruction with spatiotemporally constrained reconstruction (STCR) and (2) water/fat separation with hierarchical iterative decomposition of water and fat with echo asymmetry and least-squares estimation (HIDEAL). In healthy volunteers, spiral OIOI was acquired in the wrist during a radial-to-ulnar deviation maneuver, in the heart without breath-hold and cardiac gating, and in the lower abdomen during free-breathing for visualizing small bowel motility.

RESULTS

We demonstrate successful water/fat separated RT-MRI for all tested applications. In the wrist, resulting images provided clear depiction of ligament gaps and their interactions during the radial-to-ulnar deviation maneuver. In the heart, water/fat RT-MRI depicted epicardial fat, provided improved delineation of epicardial coronary arteries, and provided high blood-myocardial contrast for ventricular function assessment. In the abdomen, water-only RT-MRI captured small bowel mobility clearly with improved water-fat contrast.

CONCLUSIONS

We have demonstrated a novel and flexible bSSFP spiral OIOI sequence at 0.55T that can provide water/fat separated RT-MRI with a variety of application-specific temporal resolution and spatial resolution requirements.

Deep Learning for Retrospective Motion Correction in MRI: A Comprehensive Review

Motion represents one of the major challenges in magnetic resonance imaging (MRI). Since the MR signal is acquired in frequency space, any motion of the imaged object leads to complex artefacts in the reconstructed image in addition to other MR imaging artefacts. Deep learning has been frequently proposed for motion correction at several stages of the reconstruction process. The wide range of MR acquisition sequences, anatomies and pathologies of interest, and motion patterns (rigid vs. deformable and random vs. regular) makes a comprehensive solution unlikely. To facilitate the transfer of ideas between different applications, this review provides a detailed overview of proposed methods for learning-based motion correction in MRI together with their common challenges and potentials. This review identifies differences and synergies in underlying data usage, architectures, training and evaluation strategies. We critically discuss general trends and outline future directions, with the aim to enhance interaction between different application areas and research fields.

Ultra-rapid, Free-breathing, Real-time Cardiac Cine MRI Using GRASP Amplified with View Sharing and KWIC Filtering

Purpose To achieve ultra-high temporal resolution (approximately 20 msec) in free-breathing, real-time cardiac cine MRI using golden-angle radial sparse parallel (GRASP) reconstruction amplified with view sharing (VS) and k-space-weighted image contrast (KWIC) filtering. Materials and Methods Fourteen pediatric patients with congenital heart disease (mean age [SD], 9 years ± 2; 13 male) and 10 adult patients with arrhythmia (mean age, 62 years ± 8; nine male) who underwent both standard breath-hold cine and free-breathing real-time cine using GRASP were retrospectively identified. To achieve high temporal resolution, each time frame was reconstructed using six radial spokes, corresponding to acceleration factors ranging from 24 to 32. To compensate for loss in spatial resolution resulting from over-regularization in GRASP, VS and KWIC filtering were incorporated. The blur metric, visual image quality scores, and biventricular parameters were compared between clinical and real-time cine images. Results In pediatric patients, the incorporation of VS and KWIC into GRASP (ie, GRASP + VS + KWIC) produced significantly (P < .05) sharper x-y-t (blur metric: 0.36 ± 0.03, 0.41 ± 0.03, 0.48 ± 0.03, respectively) and x-y-f (blur metric: 0.28 ± 0.02, 0.31 ± 0.03, 0.37 ± 0.03, respectively) component images compared with GRASP + VS and conventional GRASP. Only the noise score differed significantly between GRASP + VS + KWIC and clinical cine; all visual scores were above the clinically acceptable (3.0) cutoff point. Biventricular volumetric parameters strongly correlated (R2 > 0.85) between clinical and real-time cine images reconstructed with GRASP + VS + KWIC and were in good agreement (relative error < 6% for all parameters). In adult patients, the visual scores of all categories were significantly lower (P < .05) for clinical cine compared with real-time cine with GRASP + VS + KWIC, except for noise (P = .08). Conclusion Incorporating VS and KWIC filtering into GRASP reconstruction enables ultra-high temporal resolution (approximately 20 msec) without significant loss in spatial resolution. Keywords: Cine, View Sharing, k-Space-weighted Image Contrast Filtering, Radial k-Space, Pediatrics, Arrhythmia, GRASP, Compressed Sensing, Real-Time, Free-Breathing Supplemental material is available for this article. © RSNA, 2024.

The future of CMR: All-in-one vs. real-time CMR (Part 2)

Cardiac Magnetic Resonance (CMR) protocols can be lengthy and complex, which has driven the research community to develop new technologies to make these protocols more efficient and patient-friendly. Two different approaches to improving CMR have been proposed, specifically "all-in-one" CMR, where several contrasts and/or motion states are acquired simultaneously, and "real-time" CMR, in which the examination is accelerated to avoid the need for breathholding and/or cardiac gating. The goal of this two-part manuscript is to describe these two different types of emerging rapid CMR protocols. To this end, the vision of all-in-one and real-time imaging are described, along with techniques which have been devised and tested along the pathway of clinical implementation. The pros and cons of the different methods are presented, and the remaining open needs of each are detailed. Part 1 tackles the "All-in-One" approaches, and Part 2 focuses on the "Real-Time" approaches along with an overall summary of these emerging methods.

Data-driven system matrix manipulation enabling fast functional imaging and intra-image nonrigid motion correction in tomography

Tomographic imaging modalities are described by large system matrices. Sparse sampling and tissue motion degrade system matrix and image quality. Various existing techniques improve the image quality without correcting the system matrices. Here, we compress the system matrices to improve computational efficiency (e.g., 42 times) using singular value decomposition and fast Fourier transform. Enabled by the efficiency, we propose (1) fast sparsely sampling functional imaging by incorporating a densely sampled prior image into the system matrix, which maintains the critical linearity while mitigating artifacts and (2) intra-image nonrigid motion correction by incorporating the motion as subdomain translations into the system matrix and reconstructing the translations together with the image iteratively. We demonstrate the methods in 3D photoacoustic computed tomography with significantly improved image qualities and clarify their applicability to X-ray CT and MRI or other types of imperfections due to the similarities in system matrices.

Synthetic multi-contrast late gadolinium enhancement imaging using post-contrast magnetic resonance fingerprinting

Late gadolinium enhancement (LGE) MRI is the non-invasive reference standard for identifying myocardial scar and fibrosis but has limitations, including difficulty delineating subendocardial scar and operator dependence on image quality. The purpose of this work is to assess the feasibility of generating multi-contrast synthetic LGE images from post-contrast T1 and T2 maps acquired using magnetic resonance fingerprinting (MRF). Fifteen consecutive patients with a history of prior ischemic cardiomyopathy (12 men; mean age 63  ±  13 years) were prospectively scanned at 1.5 T between Oct 2020 and May 2021 using conventional LGE and MRF after injection of gadolinium contrast. Three classes of synthetic LGE images were derived from MRF post-contrast T1 and T2 maps: bright-blood phase-sensitive inversion recovery (PSIR), black- and gray-blood T2 -prepared PSIR (T2 -PSIR), and a novel "tissue-optimized" image to enhance differentiation among scar, viable myocardium, and blood. Image quality was assessed on a 1-5 Likert scale by two cardiologists, and contrast was quantified as the mean absolute difference (MAD) in pixel intensities between two tissues, with different methods compared using Kruskal-Wallis with Bonferroni post hoc tests. Per-patient and per-segment scar detection rates were evaluated using conventional LGE images as reference. Image quality scores were highest for synthetic PSIR (4.0) and reference images (3.8), followed by synthetic tissue-optimized (3.3), gray-blood T2 -PSIR (3.0), and black-blood T2 -PSIR (2.6). Among synthetic images, PSIR yielded the highest myocardium/scar contrast (MAD = 0.42) but the lowest blood/scar contrast (MAD = 0.05), and vice versa for T2 -PSIR, while tissue-optimized images achieved a balance among all tissues (myocardium/scar MAD = 0.16, blood/scar MAD = 0.26, myocardium/blood MAD = 0.10). Based on reference mid-ventricular LGE scans, 13/15 patients had myocardial scar. The per-patient sensitivity/accuracy for synthetic images were the following: PSIR, 85/87%; black-blood T2 -PSIR, 62/53%; gray-blood T2 -PSIR, 100/93%; tissue optimized, 100/93%. Synthetic multi-contrast LGE images can be generated from post-contrast MRF data without additional scan time, with initial feasibility shown in ischemic cardiomyopathy patients.

MRI-based quantification of whole-organ renal metabolic rate of oxygen

During the early stages of diabetes, kidney oxygen utilization increases. The mismatch between oxygen demand and supply contributes to tissue hypoxia, a key driver of chronic kidney disease. Thus, whole-organ renal metabolic rate of oxygen (rMRO2 ) is a potentially valuable biomarker of kidney function. The key parameters required to determine rMRO2 include the renal blood flow rate (RBF) in the feeding artery and oxygen saturation in the draining renal vein (SvO2 ). However, there is currently no noninvasive method to quantify rMRO2 in absolute physiologic units. Here, a new MRI pulse sequence, Kidney Metabolism of Oxygen via T2 and Interleaved Velocity Encoding (K-MOTIVE), is described, along with evaluation of its performance in the human kidney in vivo. K-MOTIVE interleaves a phase-contrast module before a background-suppressed T2 -prepared balanced steady-state-free-precession (bSSFP) readout to measure RBF and SvO2 in a single breath-hold period of 22 s, yielding rMRO2 via Fick's principle. Variants of K-MOTIVE to evaluate alternative bSSFP readout strategies were studied. Kidney mass was manually determined from multislice gradient recalled echo images. Healthy subjects were recruited to quantify rMRO2 of the left kidney at 3-T field strength (N = 15). Assessments of repeat reproducibility and comparisons with individual measurements of RBF and SvO2 were performed, and the method's sensitivity was evaluated with a high-protein meal challenge (N = 8). K-MOTIVE yielded the following metabolic parameters: T2  = 157 ± 19 ms; SvO2  = 92% ± 6%; RBF = 400 ± 110 mL/min; and rMRO2  = 114 ± 117(μmol O2 /min)/100 g tissue. Reproducibility studies of T2 and RBF (parameters directly measured by K-MOTIVE) resulted in coefficients of variation less than 10% and intraclass correlation coefficients more than 0.75. The high-protein meal elicited an increase in rMRO2 , which was corroborated by serum biomarkers. The K-MOTIVE sequence measures SvO2 and RBF, the parameters necessary to quantify whole-organ rMRO2 , in a single breath-hold. The present work demonstrates that rMRO2 quantification is feasible with good reproducibility. rMRO2 is a potentially valuable physiological biomarker.

HyperSLICE: HyperBand optimized spiral for low-latency interactive cardiac examination

PURPOSE

Interactive cardiac MRI is used for fast scan planning and MR-guided interventions. However, the requirement for real-time acquisition and near-real-time visualization constrains the achievable spatio-temporal resolution. This study aims to improve interactive imaging resolution through optimization of undersampled spiral sampling and leveraging of deep learning for low-latency reconstruction (deep artifact suppression).

METHODS

A variable density spiral trajectory was parametrized and optimized via HyperBand to provide the best candidate trajectory for rapid deep artifact suppression. Training data consisted of 692 breath-held CINEs. The developed interactive sequence was tested in simulations and prospectively in 13 subjects (10 for image evaluation, 2 during catheterization, 1 during exercise). In the prospective study, the optimized framework-HyperSLICE- was compared with conventional Cartesian real-time and breath-hold CINE imaging in terms quantitative and qualitative image metrics. Statistical differences were tested using Friedman chi-squared tests with post hoc Nemenyi test (p < 0.05).

RESULTS

In simulations the normalized RMS error, peak SNR, structural similarity, and Laplacian energy were all statistically significantly higher using optimized spiral compared to radial and uniform spiral sampling, particularly after scan plan changes (structural similarity: 0.71 vs. 0.45 and 0.43). Prospectively, HyperSLICE enabled a higher spatial and temporal resolution than conventional Cartesian real-time imaging. The pipeline was demonstrated in patients during catheter pull back, showing sufficiently fast reconstruction for interactive imaging.

CONCLUSION

HyperSLICE enables high spatial and temporal resolution interactive imaging. Optimizing the spiral sampling enabled better overall image quality and superior handling of image transitions compared with radial and uniform spiral trajectories.

Dynamic cell tracking using time-lapse MRI with variable temporal resolution Cartesian sampling

PURPOSE

Temporal resolution of time-lapse MRI to track individual iron-labeled cells is limited by the required data-acquisition time to fill k-space and to reach sufficient SNR. Although motion of slowly patrolling monocytes can be resolved, detection of fast-moving immune cells requires improved acquisition and reconstruction strategies.

THEORY AND METHODS

For accelerated MRI cell tracking, a Cartesian sampling scheme was designed, in which the fully sampled and undersampled k-space data for different acceleration factors were acquired simultaneously, and multiple undersampling ratios could be chosen retrospectively. Compressed-sensing reconstruction was applied using dictionary learning and low-rank constraints. Detection of iron-labeled monocytes was evaluated with simulations, rotating phantom experiments and in vivo mouse brain measurements at 9.4 T.

RESULTS

Fully sampled and 2.4-times and 4.8-times accelerated images were reconstructed and had sufficient contrast-to-noise ratio (CNR) for single cells to be resolved and followed dynamically. The phantom experiments showed an improvement in CNR of 6.1% per μm/s in the 4.8-times undersampled images. Geometric distortion of cells caused by motion was visibly reduced in the accelerated images, which enabled detection of moving cells with velocities of up to 7.0 μm/s. In vivo, additional cells were resolved in the accelerated images due to the improved temporal resolution.

CONCLUSION

The easy-to-implement flexible Cartesian sampling scheme with compressed-sensing reconstruction permits simultaneous acquisition of both fully sampled and high temporal resolution images. The CNR of moving cells is effectively improved, enabling the recovery of high velocity cells with sufficient contrast at virtually no cost.

Design and development of a novel flexible ultra-short echo time (FUSE) sequence

PURPOSE

To present the validation of a new Flexible Ultra-Short Echo time (FUSE) pulse sequence using a short-T2 phantom.

METHODS

FUSE was developed to include a range of RF excitation pulses, trajectories, dimensionalities, and long-T2 suppression techniques, enabling real-time interchangeability of acquisition parameters. Additionally, we developed an improved 3D deblurring algorithm to correct for off-resonance artifacts. Several experiments were conducted to validate the efficacy of FUSE, by comparing different approaches for off-resonance artifact correction, variations in RF pulse and trajectory combinations, and long-T2 suppression techniques. All scans were performed on a 3 T system using an in-house short-T2 phantom. The evaluation of results included qualitative comparisons and quantitative assessments of the SNR and contrast-to-noise ratio.

RESULTS

Using the capabilities of FUSE, we demonstrated that we could combine a shorter readout duration with our improved deblurring algorithm to effectively reduce off-resonance artifacts. Among the different RF and trajectory combinations, the spiral trajectory with the regular half-inc pulse achieves the highest SNRs. The dual-echo subtraction technique delivers better short-T2 contrast and superior suppression of water and agar signals, whereas the off-resonance saturation method successfully suppresses water and lipid signals simultaneously.

CONCLUSION

In this work, we have validated the use of our new FUSE sequence using a short T2 phantom, demonstrating that multiple UTE acquisitions can be achieved within a single sequence. This new sequence may be useful for acquiring improved UTE images and the development of UTE imaging protocols.

Data-efficient Bayesian learning for radial dynamic MR reconstruction

BACKGROUND

Cardiac MRI has become the gold-standard imaging technique for assessing cardiovascular morphology and function. In spite of this, its slow data acquisition process presents imaging challenges due to the motion from heartbeats, respiration, and blood flow. In recent studies, deep learning (DL) algorithms have shown promising results for the task of image reconstruction. However, there have been instances where they have introduced artifacts that may be misinterpreted as pathologies or may obscure the detection of pathologies. Therefore, it is important to obtain a metric, such as the uncertainty of the network output, that identifies such artifacts. However, this can be quite challenging for large-scale image reconstruction problems such as dynamic multi-coil non-Cartesian MRI.

PURPOSE

To efficiently quantify uncertainties of a physics-informed DL-based image reconstruction method for a large-scale accelerated 2D multi-coil dynamic radial MRI reconstruction problem, and demonstrate the benefits of physics-informed DL over model-agnostic DL in reducing uncertainties while at the same time improving image quality.

METHODS

We extended a recently proposed physics-informed 2D U-Net that learns spatio-temporal slices (named XT-YT U-Net), and employed it for the task of uncertainty quantification (UQ) by using Monte Carlo dropout and a Gaussian negative log-likelihood loss function. Our data comprised 2D dynamic MR images acquired with a radial balanced steady-state free precession sequence. The XT-YT U-Net, which allows for training with a limited amount of data, was trained and validated on a dataset of 15 healthy volunteers, and further tested on data from four patients. An extensive comparison between physics-informed and model-agnostic neural networks (NNs) concerning the obtained image quality and uncertainty estimates was performed. Further, we employed calibration plots to assess the quality of the UQ.

RESULTS

The inclusion of the MR-physics model of data acquisition as a building block in the NN architecture led to higher image quality (NRMSE:- 33 ± 8.2 % $-33 \pm 8.2 \%$ , PSNR:6.3 ± 1.3 % $6.3 \pm 1.3 \%$ , and SSIM:1.9 ± 0.96 % $1.9 \pm 0.96 \%$ ), lower uncertainties (- 46 ± 8.7 % $-46 \pm 8.7 \%$ ), and, based on the calibration plots, an improved UQ compared to its model-agnostic counterpart. Furthermore, the UQ information can be used to differentiate between anatomical structures (e.g., coronary arteries, ventricle boundaries) and artifacts.

CONCLUSIONS

Using an XT-YT U-Net, we were able to quantify uncertainties of a physics-informed NN for a high-dimensional and computationally demanding 2D multi-coil dynamic MR imaging problem. In addition to improving the image quality, embedding the acquisition model in the network architecture decreased the reconstruction uncertainties as well as quantitatively improved the UQ. The UQ provides additional information to assess the performance of different network approaches.

MR Fingerprinting for Contrast Agent-free and Quantitative Characterization of Focal Liver Lesions

Purpose To evaluate the feasibility of liver MR fingerprinting (MRF) for quantitative characterization and diagnosis of focal liver lesions. Materials and Methods This single-site, prospective study included 89 participants (mean age, 62 years ± 15 [SD]; 45 women, 44 men) with various focal liver lesions who underwent MRI between October 2021 and August 2022. The participants underwent routine clinical MRI, non-contrast-enhanced liver MRF, and reference quantitative MRI with a 1.5-T MRI scanner. The bias and repeatability of the MRF measurements were assessed using linear regression, Bland-Altman plots, and coefficients of variation. The diagnostic capability of MRF-derived T1, T2, T2*, proton density fat fraction (PDFF), and a combination of these metrics to distinguish benign from malignant lesions was analyzed according to the area under the receiver operating characteristic curve (AUC). Results Liver MRF measurements showed moderate to high agreement with reference measurements (intraclass correlation = 0.94, 0.77, 0.45, and 0.61 for T1, T2, T2*, and PDFF, respectively), with underestimation of T2 values (mean bias in lesion = -0.5%, -29%, 5.8%, and -8.2% for T1, T2, T2*, and PDFF, respectively). The median coefficients of variation for repeatability of T1, T2, and T2* values were 2.5% (IQR, 3.6%), 3.1% (IQR, 5.6%), and 6.6% (IQR, 13.9%), respectively. After considering multicollinearity, a combination of MRF measurements showed a high diagnostic performance in differentiating benign from malignant lesions (AUC = 0.92 [95% CI: 0.86, 0.98]). Conclusion Liver MRF enabled the quantitative characterization of various focal liver lesions in a single breath-hold acquisition. Keywords: MR Imaging, Abdomen/GI, Liver, Imaging Sequences, Technical Aspects, Tissue Characterization, Technology Assessment, Diagnosis, Liver Lesions, MR Fingerprinting, Quantitative Characterization Supplemental material is available for this article. © RSNA, 2023.

Joint radial trajectory correction for accelerated T2 * mapping on an MR-Linac

BACKGROUND

T2 * mapping can characterize tumor hypoxia, which may be associated with resistance to therapy. Acquiring T2 * maps during MR-guided radiotherapy could inform treatment adaptation by, for example, escalating the dose to resistant sub-volumes.

PURPOSE

The purpose of this work is to demonstrate the feasibility of the accelerated T2 * mapping technique using model-based image reconstruction with integrated trajectory auto-correction (TrACR) for MR-guided radiotherapy on an MR-Linear accelerator (MR-Linac).

MATERIALS AND METHODS

The proposed method was validated in a numerical phantom, where two T2 * mapping approaches (sequential and joint) were compared for different noise levels (0,0.1,0.5,1) and gradient delays ([1, -1] and [1, -2] in units of dwell time for x- and y-axis, respectively). Fully sampled k-space was retrospectively undersampled using two different undersampling patterns. Root mean square errors (RMSEs) were calculated between reconstructed T2 * maps and ground truth. In vivo data was acquired twice weekly in one prostate and one head and neck cancer patient undergoing treatment on a 1.5 T MR-Linac. Data were retrospectively undersampled and T2 * maps reconstructed, with and without trajectory corrections were compared.

RESULTS

Numerical simulations demonstrated that, for all noise levels, T2 * maps reconstructed with a joint approach demonstrated less error compared to an uncorrected and sequential approach. For a noise level of 0.1, uniform undersampling and gradient delay [1, -1] (in units of dwell time for x- and y-axis, respectively), RMSEs for sequential and joint approaches were 13.01 and 9.32 ms, respectively, which reduced to 10.92 and 5.89 ms for a gradient delay of [1, 2]. Similarly, for alternate undersampling and gradient delay [1, -1], RMSEs for sequential and joint approaches were 9.80 and 8.90 ms, respectively, which reduced to 9.10 and 5.40 ms for gradient delay [1, 2]. For in vivo data, T2 * maps reconstructed with our proposed approach resulted in less artifacts and improved visual appearance compared to the uncorrected approach. For both prostate and head and neck cancer patients, T2 * maps reconstructed from different treatment fractions showed changes within the planning target volume (PTV).

CONCLUSION

Using the proposed approach, a retrospective data-driven gradient delay correction can be performed, which is particularly relevant for hybrid devices, where full information on the machine configuration is not available for image reconstruction. T2 * maps were acquired in under 5 min and can be integrated into MR-guided radiotherapy treatment workflows, which minimizes patient burden and leaves time for additional imaging for online adaptive radiotherapy on an MR-Linac.

Free-breathing high isotropic resolution quantitative susceptibility mapping (QSM) of liver using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction

PURPOSE

To enable free-breathing and high isotropic resolution liver quantitative susceptibility mapping (QSM) using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction.

METHODS

Using 3D multi-echo UTE cones MRI, a respiratory motion was estimated from the k-space center of the imaging data. After sorting the k-space data with estimated motion, respiratory motion state-resolved reconstruction was performed for multi-echo data followed by nonlinear least-squares fitting for proton density fat fraction (PDFF),R 2 * $$ {\mathrm{R}}_2^{\ast } $$ , and fat-corrected B0 field maps. PDFF and B0 field maps were subsequently used for QSM reconstruction. The proposed method was compared with motion-averaged (gridding) reconstruction and conventional 3D multi-echo Cartesian MRI in moving gadolinium phantom and in vivo studies. Region of interest (ROI)-based linear regression analysis was performed on these methods to investigate correlations between gadolinium concentration and QSM in the phantom study and betweenR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM in in vivo study.

RESULTS

Cones with motion-resolved reconstruction showed sharper image quality compared to motion-averaged reconstruction with a substantial reduction of motion artifacts in both moving phantom and in vivo studies. For ROI-based linear regression analysis of the phantom study, susceptibility values from cones with motion-resolved reconstruction (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$  = 0.31 × gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.05,R 2 $$ {R}^2 $$  = 0.999) and Cartesian without motion (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$  = 0.32× gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.04,R 2 $$ {R}^2 $$  = 1.000) showed linear relationships with gadolinium concentrations and showed good agreement with each other. For in vivo, motion-resolved reconstruction showed higher goodness of fit (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$  = 0.00261 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.524,R 2 $$ {R}^2 $$  = 0.977) compared to motion-averaged reconstruction (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.0021 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.572,R 2 $$ {R}^2 $$ = 0.723) in ROI-based linear regression analysis betweenR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM.

CONCLUSION

Feasibility of free-breathing liver QSM was demonstrated with motion-resolved 3D multi-echo UTE cones MRI, achieving high isotropic resolution currently unachievable in conventional Cartesian MRI.

Optimal Temporal Resolution to Achieve Good Image Quality and Perform Pharmacokinetic Analysis in Free-breathing Dynamic Contrast-enhanced MR Imaging of the Pancreas

PURPOSE

The optimal temporal resolution for free-breathing dynamic contrast-enhanced MRI (FBDCE-MRI) of the pancreas has not been determined. This study aimed to evaluate the appropriate temporal resolution to achieve good image quality and to perform pharmacokinetic analysis in FBDCE-MRI of the pancreas using golden-angle radial sparse parallel (GRASP).

METHODS

Sixteen participants (53 ± 15 years, eight females) undergoing FBDCE-MRI were included in this prospective study. Images were retrospectively reconstructed at four temporal resolutions (1.8, 3.0, 4.8, and 7.8s). Two radiologists (5 years of experience) evaluated the image quality of each reconstructed image by assessing the visualization of the celiac artery (CEA), the common hepatic artery, the splenic artery, each area of the pancreas, and artifacts using a 5-point scale. Using Tissue-4D, pharmacokinetic parameters were calculated for each area in the reconstructed images at each temporal resolution for 16 examinations, excluding two with errors in the pharmacokinetic modeling analysis. Friedman and Bonferroni tests were used for analysis. A P value < 0.05 was considered statistically significant.

RESULTS

During vascular assessment, only scores for the CEA at 7.8s were significantly lower than the other temporal resolutions. Scores of all pancreatic regions and artifacts were significantly lower at 1.8s than at 4.8s and 7.8s. In the pharmacokinetic analysis, all volume transfer coefficients (Ktrans), rate constants (Kep), and the initial area under the concentration curve (iAUC) in the pancreatic head and tail were significantly lower at 4.8s and 7.8s than at 1.8s. iAUC in the pancreatic body and extracellular extravascular volume fraction (Ve) in the pancreatic head were significantly lower at 7.8s than at 1.8s.

CONCLUSION

A temporal resolution of 3.0s is appropriate to achieve image quality and perform pharmacokinetic analysis in FBDCE-MRI of the pancreas using GRASP.

Comparative study of sub-second temporal resolution 4D-MRI and 4D-CT for target motion assessment in a phantom model

To develop and investigate the feasibility of sub-second temporal resolution volumetric T1-weighted four-dimensional (4D-) MRI in comparison with 4D-CT for respiratory-correlated motion assessment using an MRI/CT-compatible phantom. Sub-second high temporal resolution (0.5 s) gradient-echo T1-weighted 4D-MRI was developed using a volumetric acquisition scheme with compressed sensing. An MRI/CT-compatible motion phantom (simulated liver tumor) with three sinusoidal movements of amplitudes and two respiratory patterns was introduced and imaged with 4D-MRI and 4D-CT to investigate the geometric accuracy of the target movement. The geometric accuracy, including centroid position, volume, similarity index of dice similarity coefficient (DSC), and Hausdorff distance (HD), was systematically evaluated. Proposed 4D-MRI achieved a similar geometric accuracy compared with 4D-CT regarding the centroid position, volume, and similarity index. The observed position differences of the absolute average centroid were within 0.08 cm in 4D-MRI and 0.03 cm in 4D-CT, less than the 1-pixel resolution for each modality. The observed volume difference in 4D-MRI/4D-CT was within 0.73 cm3 (4.5%)/0.29 cm3 (2.1%) for a large target and 0.06 cm3 (11.3%)/0.04 cm3 (11.6%) for a small target. The observed DSC values for 4D-MRI/4D-CT were at least 0.93/0.95 for the large target and 0.83/0.84 for the small target. The maximum HD values were 0.25 cm/0.31 cm for the large target and 0.21 cm/0.15 cm for the small target. Although 4D-CT potentially exhibit superior numerical accuracy in phantom studies, the proposed high temporal resolution 4D-MRI demonstrates sub-millimetre geometric accuracy comparable to that of 4D-CT. These findings suggest that the 4D-MRI technique is a viable option for characterizing motion and generating phase-dependent internal target volumes within the realm of radiotherapy.

Motion-compensated T1 mapping in cardiovascular magnetic resonance imaging: a technical review

T 1 mapping is becoming a staple magnetic resonance imaging method for diagnosing myocardial diseases such as ischemic cardiomyopathy, hypertrophic cardiomyopathy, myocarditis, and more. Clinically, most T 1 mapping sequences acquire a single slice at a single cardiac phase across a 10 to 15-heartbeat breath-hold, with one to three slices acquired in total. This leaves opportunities for improving patient comfort and information density by acquiring data across multiple cardiac phases in free-running acquisitions and across multiple respiratory phases in free-breathing acquisitions. Scanning in the presence of cardiac and respiratory motion requires more complex motion characterization and compensation. Most clinical mapping sequences use 2D single-slice acquisitions; however newer techniques allow for motion-compensated reconstructions in three dimensions and beyond. To further address confounding factors and improve measurement accuracy, T 1 maps can be acquired jointly with other quantitative parameters such as T 2 , T 2 ∗ , fat fraction, and more. These multiparametric acquisitions allow for constrained reconstruction approaches that isolate contributions to T 1 from other motion and relaxation mechanisms. In this review, we examine the state of the literature in motion-corrected and motion-resolved T 1 mapping, with potential future directions for further technical development and clinical translation.

Quantitative longitudinal mapping of radiation-treated prostate cancer using MR fingerprinting with radial acquisition and subspace reconstruction

MR fingerprinting (MRF) enables fast multiparametric quantitative imaging with a single acquisition and has been shown to improve diagnosis of prostate cancer. However, most prostate MRF studies were performed with spiral acquisitions that are sensitive to B0 inhomogeneities and consequent blurring. In this work, a radial MRF acquisition with a novel subspace reconstruction technique was developed to enable fast T1/T2 mapping in the prostate in under 4 min. The subspace reconstruction exploits the extensive temporal correlations in the MRF dictionary to pre-compute a low dimensional space for the solution and thus reduce the number of radial spokes to accelerate the acquisition. Iterative reconstruction with the subspace model and additional regularization of the signal representation in the subspace is performed to minimize the number of spokes and maintain matching quality and SNR. Reconstruction accuracy was assessed using the ISMRM NIST phantom. In-vivo validation was performed on two healthy subjects and two prostate cancer patients undergoing radiation therapy. The longitudinal repeatability was quantified using the concordance correlation coefficient (CCC) in one of the healthy subjects by repeated scans over 1 year. One prostate cancer patient was scanned at three time points, before initiating therapy and following brachytherapy and external beam radiation. Changes in the T1/T2 maps obtained with the proposed method were quantified. The prostate, peripheral and transitional zones, and visible dominant lesion were delineated for each study, and the statistics and distribution of the quantitative mapping values were analyzed. Significant image quality improvements compared with standard reconstruction methods were obtained with the proposed subspace reconstruction method. A notable decrease in the spread of the T1/T2 values without biasing the estimated mean values was observed with the subspace reconstruction and agreed with reported literature values. The subspace reconstruction enabled visualization of small differences in T1/T2 values in the tumor region within the peripheral zone. Longitudinal imaging of a volunteer subject yielded CCC of 0.89 for MRF T1, and 0.81 for MRF T2 in the prostate gland. Longitudinal imaging of the prostate patient confirmed the feasibility of capturing radiation treatment related changes. This work is a proof-of-concept for a high resolution and fast quantitative mapping using golden-angle radial MRF combined with a subspace reconstruction technique for longitudinal treatment response assessment in subjects undergoing radiation treatment.

Rapid 3D T1 mapping using deep learning-assisted Look-Locker inversion recovery MRI

PURPOSE

Conventional 3D Look-Locker inversion recovery (LLIR) T1 mapping requires multi-repetition data acquisition to reconstruct images at different inversion times for T1 fitting. To ensure B1 robustness, sufficient time of delay (TD) is needed between repetitions, which prolongs scan time. This work proposes a novel deep learning-assisted LLIR MRI approach for rapid 3D T1 mapping without TD.

THEORY AND METHODS

The proposed approach is based on the fact thatT 1 * $$ {\mathrm{T}}_1^{\ast } $$ , the effective T1 in LLIR imaging, is independent of TD and can be estimated from both LLIR imaging with and without TD, while accurate conversion ofT 1 * $$ {\mathrm{T}}_1^{\ast } $$ to T1 requires TD. Therefore, deep learning can be used to learn the conversion ofT 1 * $$ {\mathrm{T}}_1^{\ast } $$ to T1 , which eliminates the need for TD. This idea was implemented for inversion-recovery-prepared Golden-angel RAdial Sparse Parallel T1 mapping (GraspT1 ). 39 GraspT1 datasets with a TD of 6 s (GraspT1 -TD6) were used for training, which also incorporates additional anatomical images. The trained network was applied for T1 estimation in 14 GraspT1 datasets without TD (GraspT1 -TD0). The robustness of the trained network was also tested.

RESULTS

Deep learning-based T1 estimation from GraspT1 -TD0 is accurate compared to the reference. Incorporation of additional anatomical images improves the accuracy of T1 estimation. The technique is also robust against slight variation in spatial resolution, imaging orientation and scanner platform.

CONCLUSION

Our approach eliminates the need for TD in 3D LLIR imaging without affecting the T1 estimation accuracy. It represents a novel use of deep learning towards more efficient and robust 3D LLIR T1 mapping.

A low-rank deep image prior reconstruction for free-breathing ungated spiral functional CMR at 0.55 T and 1.5 T

OBJECTIVE

This study combines a deep image prior with low-rank subspace modeling to enable real-time (free-breathing and ungated) functional cardiac imaging on a commercial 0.55 T scanner.

MATERIALS AND METHODS

The proposed low-rank deep image prior (LR-DIP) uses two u-nets to generate spatial and temporal basis functions that are combined to yield dynamic images, with no need for additional training data. Simulations and scans in 13 healthy subjects were performed at 0.55 T and 1.5 T using a golden angle spiral bSSFP sequence with images reconstructed using [Formula: see text]-ESPIRiT, low-rank plus sparse (L + S) matrix completion, and LR-DIP. Cartesian breath-held ECG-gated cine images were acquired for reference at 1.5 T. Two cardiothoracic radiologists rated images on a 1-5 scale for various categories, and LV function measurements were compared.

RESULTS

LR-DIP yielded the lowest errors in simulations, especially at high acceleration factors (R [Formula: see text] 8). LR-DIP ejection fraction measurements agreed with 1.5 T reference values (mean bias - 0.3% at 0.55 T and - 0.2% at 1.5 T). Compared to reference images, LR-DIP images received similar ratings at 1.5 T (all categories above 3.9) and slightly lower at 0.55 T (above 3.4).

CONCLUSION

Feasibility of real-time functional cardiac imaging using a low-rank deep image prior reconstruction was demonstrated in healthy subjects on a commercial 0.55 T scanner.

Free-Breathing Liver Fat, R₂* and B₀ Field Mapping Using Multi-Echo Radial FLASH and Regularized Model-Based Reconstruction

This work introduced a stack-of-radial multi-echo asymmetric-echo MRI sequence for free-breathing liver volumetric acquisition. Regularized model-based reconstruction was implemented in Berkeley Advanced Reconstruction Toolbox (BART) to jointly estimate all physical parameter maps (water, fat, R2∗ , and B0 field inhomogeneity maps) and coil sensitivity maps from self-gated k -space data. Specifically, locally low rank and temporal total variation regularization were employed directly on physical parameter maps. The proposed free-breathing radial technique was tested on a water/fat & iron phantom, a young volunteer, and obesity/diabetes/hepatic steatosis patients. Quantitative fat fraction and R2∗ accuracy were confirmed by comparing our technique with the reference breath-hold Cartesian scan. The multi-echo radial sampling sequence achieves fast k -space coverage and is robust to motion. Moreover, the proposed motion-resolved model-based reconstruction allows for free-breathing liver fat and R2∗ quantification in multiple motion states. Overall, our proposed technique offers a convenient tool for non-invasive liver assessment with no breath holding requirement.

Optimization of 4D combined angiography and perfusion using radial imaging and arterial spin labeling

PURPOSE

To extend and optimize a non-contrast MRI technique to obtain whole head 4D (time-resolved 3D) qualitative angiographic and perfusion images from a single scan.

METHODS

4D combined angiography and perfusion using radial imaging and arterial spin labeling (CAPRIA) uses pseudocontinuous labeling with a 3D golden ratio ("koosh ball") readout to continuously image the blood water as it travels through the arterial system and exchanges into the tissue. High spatial/temporal resolution angiograms and low spatial/temporal resolution perfusion images can be flexibly reconstructed from the same raw k-space data. Constant and variable flip angle (CFA and VFA, respectively) excitation schedules were optimized through simulations and tested in healthy volunteers. A conventional sensitivity encoding (SENSE) reconstruction was compared against a locally low rank (LLR) reconstruction, which leverages spatiotemporal correlations. Comparison was also made with time-matched time-of-flight angiography and multi-delay EPI perfusion images. Differences in image quality were assessed through split-scan repeatability.

RESULTS

The optimized VFA schedule (2-9°) resulted in a significant (p < 0.001) improvement in image quality (up to 84% vs. CFA), particularly for the lower SNR perfusion images. The LLR reconstruction provided effective denoising without biasing the signal timecourses, significantly improving angiographic and perfusion image quality and repeatability (up to 143%, p < 0.001). 4D CAPRIA performed well compared with time-of-flight angiography and had better perfusion signal repeatability than the EPI-based approach (p < 0.001).

CONCLUSION

4D CAPRIA optimized using a VFA schedule and LLR reconstruction can yield high quality whole head 4D angiograms and perfusion images from a single scan.

Free-breathing half-radial dual-echo balanced steady-state free precession thoracic imaging with wobbling Archimedean spiral pole trajectories

PURPOSE

To demonstrate free-breathing thoracic MRI with a minimal-TR balanced steady-state free precession (bSSFP) technique using wobbling Archimedean spiral pole (WASP) trajectories.

METHODS

Phantom and free-breathing in vivo chest imaging in healthy volunteers was performed at 1.5T with a half-radial, dual-echo, bSSFP sequence, termed bSTAR. For maximum sampling efficiency, a single analog-to-digital converter window along the full bipolar readout was used. To ensure a homogeneous coverage of the k-space over multiple breathing cycles, radial k-space sampling followed short-duration Archimedean spiral interleaves that were randomly titled by a small polar angle and rotated by a golden angle about the polar axis; depticting a wobbling Archimedean spiral pole (WASP) trajectory. In phantom and in vivo experiments, WASP trajectories were compared to spiral phyllotaxis sampling in terms of eddy currents and were used to generate in vivo thorax images at different respiratory phases.

RESULTS

WASP trajectories provided artifact-free bSTAR imaging in both phantom and in vivo and respiratory self-gated reconstruction was successfully performed in all subjects. The amount of the acquired data allowed the reconstruction of 10 volumes at different respiratory levels with isotropic resolution of 1.77mm from a scan of 5.5minutes (using a TR of 1.32ms), and one high-resolution 1.16mm end-expiratory volume from a scan of 4.7minutes (using a TR of 1.42ms). The very short TR of bSTAR mitigated off-resonance artifacts despite the large field-of-view.

CONCLUSION

We have demonstrated the feasibility of high-resolution free-breathing thoracic imaging with bSTAR using the wobbling Archimedean spiral pole in healthy subjects at 1.5T.

Single-spoke binning: Reducing motion artifacts in abdominal radial stack-of-stars imaging

PURPOSE

To increase the effectiveness of respiratory gating in radial stack-of-stars MRI, particularly when imaging at high spatial resolutions or with multiple echoes.

METHODS

Free induction decay (FID) navigators were integrated into a three-dimensional gradient echo radial stack-of-stars pulse sequence. These navigators provided a motion signal with a high temporal resolution, which allowed single-spoke binning (SSB): each spoke at each phase encode step was sorted individually to the corresponding motion state of the respiratory signal. SSB was compared with spoke-angle binning (SAB), in which all phase encode steps of one projection angle were sorted without the use of additional navigator data. To illustrate the benefit of SSB over SAB, images of a motion phantom and of six free-breathing volunteers were reconstructed after motion-gating using either method. Image sharpness was quantitatively compared using image gradient entropies.

RESULTS

The proposed method resulted in sharper images of the motion phantom and free-breathing volunteers. Differences in gradient entropy were statistically significant (p = 0.03) in favor of SSB. The increased accuracy of motion-gating led to a decrease of streaking artifacts in motion-gated four-dimensional reconstructions. To consistently estimate respiratory signals from the FID-navigator data, specific types of gradient spoiler waveforms were required.

CONCLUSION

SSB allowed high-resolution motion-corrected MR imaging, even when acquiring multiple gradient echo signals or large acquisition matrices, without sacrificing accuracy of motion-gating. SSB thus relieves restrictions on the choice of pulse sequence parameters, enabling the use of motion-gated radial stack-of-stars MRI in a broader domain of clinical applications.

Time-encoded pseudo-continuous arterial spin labeling: Increasing SNR in ASL dynamic angiography

PURPOSE

Dynamic angiography using arterial spin labeling (ASL) can provide detailed hemodynamic information. However, the long time-resolved readouts require small flip angles to preserve ASL signal for later timepoints, limiting SNR. By using time-encoded ASL to generate temporal information, the readout can be shortened. Here, the SNR improvements from using larger flip angles, made possible by the shorter readout, are quantitatively investigated.

METHODS

The SNR of a conventional protocol with nine Look-Locker readouts and a 4 × $$ \times $$ 3 time-encoded protocol with three Look-Locker readouts (giving nine matched timepoints) were compared using simulations and in vivo data. Both protocols were compared using readouts with constant flip angles (CFAs) and variable flip angles (VFAs), where the VFA scheme was designed to produce a consistent ASL signal across readouts. Optimization of the background suppression to minimize physiological noise across readouts was also explored.

RESULTS

The time-encoded protocol increased in vivo SNR by 103% and 96% when using CFAs or VFAs, respectively. Use of VFAs improved SNR compared with CFAs by 25% and 21% for the conventional and time-encoded protocols, respectively. The VFA scheme also removed signal discontinuities in the time-encoded data. Preliminary data suggest that optimizing the background suppression could improve in vivo SNR by a further 16%.

CONCLUSIONS

Time encoding can be used to generate additional temporal information in ASL angiography. This enables the use of larger flip angles, which can double the SNR compared with a non-time-encoded protocol. The shortened time-encoded readout can also lead to improved background suppression, reducing physiological noise and further improving SNR.

Kz-accelerated variable-density stack-of-stars MRI

This work aimed to develop a modified stack-of-stars golden-angle radial sampling scheme with variable-density acceleration along the slice (kz) dimension (referred to as VD-stack-of-stars) and to test this new sampling trajectory with multi-coil compressed sensing reconstruction for rapid motion-robust 3D liver MRI. VD-stack-of-stars sampling implements additional variable-density undersampling along the kz dimension, so that slice resolution (or volumetric coverage) can be increased without prolonging scan time. The new sampling trajectory (with increased slice resolution) was compared with standard stack-of-stars sampling with fully sampled kz (with standard slice resolution) in both non-contrast-enhanced free-breathing liver MRI and dynamic contrast-enhanced MRI (DCE-MRI) of the liver in volunteers. For both sampling trajectories, respiratory motion was extracted from the acquired radial data, and images were reconstructed using motion-compensated (respiratory-resolved or respiratory-weighted) dynamic radial compressed sensing reconstruction techniques. Qualitative image quality assessment (visual assessment by experienced radiologists) and quantitative analysis (as a metric of image sharpness) were performed to compare images acquired using the new and standard stack-of-stars sampling trajectories. Compared to standard stack-of-stars sampling, both non-contrast-enhanced and DCE liver MR images acquired with VD-stack-of-stars sampling presented improved overall image quality, sharper liver edges and increased hepatic vessel clarity in all image planes. The results have suggested that the proposed VD-stack-of-stars sampling scheme can achieve improved performance (increased slice resolution or volumetric coverage with better image quality) over standard stack-of-stars sampling in free-breathing DCE-MRI without increasing scan time. The reformatted coronal and sagittal images with better slice resolution may provide added clinical value.

Single-heartbeat cardiac cine imaging via jointly regularized nonrigid motion-corrected reconstruction

The aim of the current study was to develop a novel approach for 2D breath-hold cardiac cine imaging from a single heartbeat, by combining cardiac motion-corrected reconstructions and nonrigidly aligned patch-based regularization. Conventional cardiac cine imaging is obtained via motion-resolved reconstructions of data acquired over multiple heartbeats. Here, we achieve single-heartbeat cine imaging by incorporating nonrigid cardiac motion correction into the reconstruction of each cardiac phase, in conjunction with a motion-aligned patch-based regularization. The proposed Motion-Corrected CINE (MC-CINE) incorporates all acquired data into the reconstruction of each (motion-corrected) cardiac phase, resulting in a better posed problem than motion-resolved approaches. MC-CINE was compared with iterative sensitivity encoding (itSENSE) and Extra-Dimensional Golden Angle Radial Sparse Parallel (XD-GRASP) in 14 healthy subjects in terms of image sharpness, reader scoring (range: 1-5) and reader ranking (range: 1-9) of image quality, and single-slice left ventricular assessment. MC-CINE was significantly superior to both itSENSE and XD-GRASP using 20 heartbeats, two heartbeats, and one heartbeat. Iterative SENSE, XD-GRASP, and MC-CINE achieved a sharpness of 74%, 74%, and 82% using 20 heartbeats, and 53%, 66%, and 82% with one heartbeat, respectively. The corresponding results for reader scoring were 4.0, 4.7, and 4.9 with 20 heartbeats, and 1.1, 3.0, and 3.9 with one heartbeat. The corresponding results for reader ranking were 5.3, 7.3, and 8.6 with 20 heartbeats, and 1.0, 3.2, and 5.4 with one heartbeat. MC-CINE using a single heartbeat presented nonsignificant differences in image quality to itSENSE with 20 heartbeats. MC-CINE and XD-GRASP at one heartbeat both presented a nonsignificant negative bias of less than 2% in ejection fraction relative to the reference itSENSE. It was concluded that the proposed MC-CINE significantly improves image quality relative to itSENSE and XD-GRASP, enabling 2D cine from a single heartbeat.

Deep Learning-Based Reconstruction for Cardiac MRI: A Review

Cardiac magnetic resonance (CMR) is an essential clinical tool for the assessment of cardiovascular disease. Deep learning (DL) has recently revolutionized the field through image reconstruction techniques that allow unprecedented data undersampling rates. These fast acquisitions have the potential to considerably impact the diagnosis and treatment of cardiovascular disease. Herein, we provide a comprehensive review of DL-based reconstruction methods for CMR. We place special emphasis on state-of-the-art unrolled networks, which are heavily based on a conventional image reconstruction framework. We review the main DL-based methods and connect them to the relevant conventional reconstruction theory. Next, we review several methods developed to tackle specific challenges that arise from the characteristics of CMR data. Then, we focus on DL-based methods developed for specific CMR applications, including flow imaging, late gadolinium enhancement, and quantitative tissue characterization. Finally, we discuss the pitfalls and future outlook of DL-based reconstructions in CMR, focusing on the robustness, interpretability, clinical deployment, and potential for new methods.

Ultra-wide band radar for prospective respiratory motion correction in the liver

Objective.T1 mapping of the liver is time consuming and can be challenging due to respiratory motion. Here we present a prospective slice tracking approach, which utilizes an external ultra-wide band radar signal and allows for efficient T1 mapping during free-breathing.Approach.The fast radar signal is calibrated to an MR-based motion signal to create a motion model. This motion model provides motion estimates, which are used to carry out slice tracking for any subsequent clinical scan. This approach was evaluated in simulations, phantom experiments andin vivoscans.Main results.Radar-based slice tracking was implemented on an MR system with a total latency of 77 ms. Moving phantom experiments showed accurate motion prediction with an error of 0.12 mm in anterior-posterior and 0.81 mm in head-feet direction. The model error remained stable for up to two hours.In vivoexperiments showed visible image improvement with a motion model error three times smaller than with a respiratory bellow. For T1 mapping during free-breathing the proposed approach provided similar results compared to reference T1 mapping during a breathhold.Significance.The proposed radar-based approach achieves accurate slice tracking and enables efficient T1 mapping of the liver during free-breathing. This motion correction approach is independent from scanning parameters and could also be used for applications like MR guided radiotherapy or MR Elastography.

MP-RAVE: IR-Prepared T1 -Weighted Radial Stack-of-Stars 3D GRE imaging with retrospective motion correction

PURPOSE

To describe an inversion-recovery T₁ -weighted radial stack-of-stars 3D gradient echo (GRE) sequence with comparable image quality to conventional MP-RAGE and to demonstrate how the radial acquisition scheme can be utilized for additional retrospective motion correction to improve robustness to head motion.

METHODS

The proposed sequence, named MP-RAVE, has been derived from a previously described radial stack-of-stars 3D GRE sequence (RAVE) and includes a 180° inversion recovery pulse that is generated once for every stack of radial views. The sequence is combined with retrospective 3D motion correction to improve robustness. The effectiveness has been evaluated in phantoms and healthy volunteers and compared to conventional MP-RAGE acquisition.

RESULTS

MP-RAGE and MP-RAVE anatomical images were rated "good" to "excellent" in overall image quality, with artifact level between "mild" and "no artifacts", and with no statistically significant difference between methods. During head motion, MP-RAVE showed higher inherent robustness with artifacts confined to local brain regions. In combination with motion correction, MP-RAVE provided noticeably improved image quality during different head motion and showed statistically significant improvement in image sharpness.

CONCLUSION

MP-RAVE provides comparable image quality and contrast to conventional MP-RAGE with improved robustness to head motion. In combination with retrospective 3D motion correction, MP-RAVE can be a useful alternative to MP-RAGE, especially in non-cooperative or pediatric patients.

Latest Advances in Image Acceleration: All Dimensions are Fair Game

Magnetic resonance imaging (MRI) is a versatile modality that can generate high-resolution images with a variety of tissue contrasts. However, MRI is a slow technique and requires long acquisition times, which increase with higher temporal and spatial resolution and/or when multiple contrasts and large volumetric coverage is required. In order to speedup MR data acquisition, several approaches have been introduced in the literature. Most of these techniques acquire less data than required and exploit intrinsic redundancies in the MR images to recover the information that was not sampled. This article presents a review of MR acquisition and reconstruction methods that have exploited redundancies in the temporal, spatial, and contrast/parametric dimensions to accelerate image data acquisition, focusing on cardiac and abdominal MR imaging applications. The review describes how each of these dimensions has been separately exploited for speeding up MR acquisition to then discuss more advanced techniques where multiple dimensions are exploited together for further reducing scan times. Finally, future directions for multidimensional image acceleration and remaining technical challenges are discussed. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: 1.

4D Golden-Angle Radial MRI at Subsecond Temporal Resolution

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

Performance of spiral UTE-MRI of the lung in post-COVID patients

Patients recovered from COVID-19 may develop long-COVID symptoms in the lung. For this patient population (post-COVID patients), they may benefit from longitudinal, radiation-free lung MRI exams for monitoring lung lesion development and progression. The purpose of this study was to investigate the performance of a spiral ultrashort echo time MRI sequence (Spiral-VIBE-UTE) in a cohort of post-COVID patients in comparison with CT and to compare image quality obtained using different spiral MRI acquisition protocols. Lung MRI was performed in 36 post-COVID patients with different acquisition protocols, including different spiral sampling reordering schemes (line in partition or partition in line) and different breath-hold positions (inspiration or expiration). Three experienced chest radiologists independently scored all the MR images for different pulmonary structures. Lung MR images from spiral acquisition protocol that received the highest image quality scores were also compared against corresponding CT images in 27 patients for evaluating diagnostic image quality and lesion identification. Spiral-VIBE-UTE MRI acquired with the line in partition reordering scheme in an inspiratory breath-holding position achieved the highest image quality scores (score range = 2.17-3.69) compared to others (score range = 1.7-3.29). Compared to corresponding chest CT images, three readers found that 81.5% (22 out of 27), 81.5% (22 out of 27) and 37% (10 out of 27) of the MR images were useful, respectively. Meanwhile, they all agreed that MRI could identify significant lesions in the lungs. The Spiral-VIBE-UTE sequence allows for fast imaging of the lung in a single breath hold. It could be a valuable tool for lung imaging without radiation and could provide great value for managing different lung diseases including assessment of post-COVID lesions.

Dynamic Cardiac Magnetic Resonance Fingerprinting During Vasoactive Breathing Maneuvers: First Results

BACKGROUND

Cardiac magnetic resonance fingerprinting (cMRF) enables simultaneous mapping of myocardial T1 and T2 with very short acquisition times. Breathing maneuvers have been utilized as a vasoactive stress test to dynamically characterize myocardial tissue in vivo. We tested the feasibility of sequential, rapid cMRF acquisitions during breathing maneuvers to quantify myocardial T1 and T2 changes.

METHODS

We measured T1 and T2 values using conventional T1 and T2-mapping techniques (modified look locker inversion [MOLLI] and T2-prepared balanced-steady state free precession), and a 15 heartbeat (15-hb) and rapid 5-hb cMRF sequence in a phantom and in 9 healthy volunteers. The cMRF_5-hb sequence was also used to dynamically assess T1 and T2 changes over the course of a vasoactive combined breathing maneuver.

RESULTS

In healthy volunteers, the mean myocardial T1 of the different mapping methodologies were: MOLLI 1,224 ± 81 ms, cMRF_15-hb 1,359 ± 97 ms, and cMRF_5-hb 1,357 ± 76 ms. The mean myocardial T2 measured with the conventional mapping technique was 41.7 ± 6.7 ms, while for cMRF_15-hb 29.6 ± 5.8 ms and cMRF_5-hb 30.5 ± 5.8 ms. T2 was reduced with vasoconstriction (post-hyperventilation compared to a baseline resting state) (30.15 ± 1.53 ms vs. 27.99 ± 2.07 ms, p = 0.02), while T1 did not change with hyperventilation. During the vasodilatory breath-hold, no significant change of myocardial T1 and T2 was observed.

CONCLUSIONS

cMRF_5-hb enables simultaneous mapping of myocardial T1 and T2, and may be used to track dynamic changes of myocardial T1 and T2 during vasoactive combined breathing maneuvers.

MR Fingerprinting for Liver Tissue Characterization: A Histopathologic Correlation Study

Background Liver MR fingerprinting (MRF) enables simultaneous quantification of T1, T2, T2*, and proton density fat fraction (PDFF) maps in single breath-hold acquisitions. Histopathologic correlation studies are desired for its clinical use. Purpose To compare liver MRF-derived metrics with separate reference quantitative MRI in participants with diffuse liver disease, evaluate scan-rescan repeatability of liver MRF, and validate MRF-derived measurements for histologic grading of liver biopsies. Materials and Methods This prospective study included participants with diffuse liver disease undergoing MRI from July 2021 to January 2022. Participants underwent two-dimensional single-section liver MRF and separate reference quantitative MRI. Linear regression, Bland-Altman plots, and coefficients of variation were used to assess the bias and repeatability of liver MRF measurements. For participants undergoing liver biopsy, the association between mapping and histologic grading was evaluated by using the Spearman correlation coefficient. Results Fifty-six participants (mean age, 59 years ± 15 [SD]; 32 women) were included to compare mapping techniques and 23 participants were evaluated with liver biopsy (mean age, 52.7 years ± 12.7; 14 women). The linearity of MRF with reference measurements in participants with diffuse liver disease (R² value) for T1, T2, T2*, and PDFF maps was 0.86, 0.88, 0.54, and 0.99, respectively. The overall coefficients of variation for repeatability in the liver were 3.2%, 5.5%, 7.1%, and 4.6% for T1, T2, T2*, and PDFF maps, respectively. MRF-derived metrics showed high diagnostic performance in differentiating moderate or severe changes from mild or no changes (area under the receiver operating characteristic curve for fibrosis, inflammation, steatosis, and siderosis: 0.62 [95% CI: 0.52, 0.62], 0.92 [95% CI: 0.88, 0.92], 0.97 [95% CI: 0.96, 0.97], and 0.74 [95% CI: 0.57, 0.74], respectively). Conclusion Liver MR fingerprinting provided repeatable T1, T2, T2*, and proton density fat fraction maps in high agreement with reference quantitative mapping and may correlate with pathologic grades in participants with diffuse liver disease. © RSNA, 2022 Online supplemental material is available for this article.

A novel multipurpose device for guided knee motion and loading during dynamic magnetic resonance imaging

INTRODUCTION

This work aimed to develop a novel multipurpose device for guided knee flexion-extension, both passively using a motorized pneumatic system and actively (muscle-driven) with the joint unloaded or loaded during dynamic MRI. Secondary objectives were to characterize the participant experience during device use, and present preliminary dynamic MRI data to demonstrate the different device capabilities.

MATERIAL AND METHODS

Self-reported outcomes were used to characterize the pain, physical exertion and discomfort levels experienced by 10 healthy male participants during four different active knee motion and loading protocols using the novel device. Knee angular data were recorded during the protocols to determine the maximum knee range of motion achievable. Dynamic MRI was acquired for three healthy volunteers during passive, unloaded knee motion using 2D Cartesian TSE, 2D radial GRE and 3D UTE sequences; and during active, unloaded and loaded knee motion using 2D radial GRE imaging. Because of the different MRI sequences used, spatial resolution was inherently lower for active knee motion than for passive motion acquisitions.

RESULTS

Depending on the protocol, some participants reported slight pain, mild discomfort and varying levels of physical exertion. On average, participants achieved ∼40° of knee flexion; loaded conditions can create knee moments up to 27Nm. High quality imaging data were obtained during different motion and loading conditions. Dynamic 3D data allowed to retrospectively extract arbitrarily oriented slices.

CONCLUSION

A novel multipurpose device for guided, physiologically relevant knee motion and loading during dynamic MRI was developed. Device use was well tolerated and suitable for acquiring high quality images during different motion and loading conditions. Different bone positions between loaded and unloaded conditions were likely due to out-of-plane motion, particularly because image registration was not performed. Ultimately, this device could be used to advance our understanding of physiological and pathological joint mechanics.

Recent Technical Developments in ASL: A Review of the State of the Art

This review article provides an overview of a range of recent technical developments in advanced arterial spin labeling (ASL) methods that have been developed or adopted by the community since the publication of a previous ASL consensus paper by Alsop et al. It is part of a series of review/recommendation papers from the International Society for Magnetic Resonance in Medicine Perfusion Study Group. Here, we focus on advancements in readouts and trajectories, image reconstruction, noise reduction, partial volume correction, quantification of nonperfusion parameters, fMRI, fingerprinting, vessel selective ASL, angiography, deep learning, and ultrahigh field ASL. We aim to provide a high level of understanding of these new approaches and some guidance for their implementation, with the goal of facilitating the adoption of such advances by research groups and by MRI vendors. Topics outside the scope of this article that are reviewed at length in separate articles include velocity selective ASL, multiple-timepoint ASL, body ASL, and clinical ASL recommendations.

B-Spline Parameterized Joint Optimization of Reconstruction and K-Space Trajectories (BJORK) for Accelerated 2D MRI

Optimizing k-space sampling trajectories is a promising yet challenging topic for fast magnetic resonance imaging (MRI). This work proposes to optimize a reconstruction method and sampling trajectories jointly concerning image reconstruction quality in a supervised learning manner. We parameterize trajectories with quadratic B-spline kernels to reduce the number of parameters and apply multi-scale optimization, which may help to avoid sub-optimal local minima. The algorithm includes an efficient non-Cartesian unrolled neural network-based reconstruction and an accurate approximation for backpropagation through the non-uniform fast Fourier transform (NUFFT) operator to accurately reconstruct and back-propagate multi-coil non-Cartesian data. Penalties on slew rate and gradient amplitude enforce hardware constraints. Sampling and reconstruction are trained jointly using large public datasets. To correct for possible eddy-current effects introduced by the curved trajectory, we use a pencil-beam trajectory mapping technique. In both simulations and in- vivo experiments, the learned trajectory demonstrates significantly improved image quality compared to previous model-based and learning-based trajectory optimization methods for 10× acceleration factors. Though trained with neural network-based reconstruction, the proposed trajectory also leads to improved image quality with compressed sensing-based reconstruction.

Metabolism of oxygen via T2 and interleaved velocity encoding: A rapid method to quantify whole-brain cerebral metabolic rate of oxygen

PURPOSE

Cerebral metabolic rate of oxygen (CMRO2 ) is an important biomarker of brain function. Key physiological parameters required to quantify CMRO2 include blood flow rate in the feeding arteries and venous oxygen saturation (SvO2 ) in the draining vein. Here, a pulse sequence, metabolism of oxygen via T2 and interleaved velocity encoding (MOTIVE), was developed to measure both parameters simultaneously and enable CMRO2 quantification in a single pass.

METHODS

The MOTIVE sequence interleaves a phase-contrast module between a nonselective saturation and a background-suppressed T2 -prepared EPI readout (BGS-EPI) to measure T2 of blood water protons and cerebral blood flow in 20 s or less. The MOTIVE and standalone BGS-EPI sequences were compared against TRUST ("T2 relaxation under spin tagging") in the brain in healthy subjects (N = 24). Variants of MOTIVE to enhance resolution or shorten scan time were explored. Intrasession and intersession reproducibility studies were performed.

RESULTS

MOTIVE experiments yielded an average SvO2 of 61 ± 6% in the superior sagittal sinus of the brain and an average cerebral blood flow of 56 ± 10 ml/min/100 g. The bias in SvO2 of MOTIVE and BGS-EPI to TRUST was +2 ± 4% and +1 ± 3%, respectively. The bias in cerebral blood flow of MOTIVE to Cartesian phase-contrast reference was +1 ± 6 ml/min/100 g.

CONCLUSIONS

The MOTIVE sequence is an advance over existing T2 -based oximetric methods. It does not require a control image and simultaneously measures SvO2 and flow velocity. The measurements agree well with TRUST and reference phase-contrast sequences. This noninvasive technique enables CMRO2 quantification in under 20 s and is reproducible for in vivo applications.

Dynamic Contrast-Enhanced MRI in the Abdomen of Mice with High Temporal and Spatial Resolution Using Stack-of-Stars Sampling and KWIC Reconstruction

Application of quantitative dynamic contrast-enhanced (DCE) MRI in mouse models of abdominal cancer is challenging due to the effects of RF inhomogeneity, image corruption from rapid respiratory motion and the need for high spatial and temporal resolutions. Here we demonstrate a DCE protocol optimized for such applications. The method consists of three acquisitions: (1) actual flip-angle B₁ mapping, (2) variable flip-angle T₁ mapping and (3) acquisition of the DCE series using a motion-robust radial strategy with k-space weighted image contrast (KWIC) reconstruction. All three acquisitions employ spoiled radial imaging with stack-of-stars sampling (SoS) and golden-angle increments between the views. This scheme is shown to minimize artifacts due to respiratory motion while simultaneously facilitating view-sharing image reconstruction for the dynamic series. The method is demonstrated in a genetically engineered mouse model of pancreatic ductal adenocarcinoma and yielded mean perfusion parameters of K^trans = 0.23 ± 0.14 min⁻¹ and v_e = 0.31 ± 0.17 (n = 22) over a wide range of tumor sizes. The SoS-sampled DCE method is shown to produce artifact-free images with good SNR leading to robust estimation of DCE parameters.

Real-time MRI motion estimation through an unsupervised k-space-driven deformable registration network (KS-RegNet)

Purpose. Real-time three-dimensional (3D) magnetic resonance (MR) imaging is challenging because of slow MR signal acquisition, leading to highly under-sampled k-space data. Here, we proposed a deep learning-based, k-space-driven deformable registration network (KS-RegNet) for real-time 3D MR imaging. By incorporating prior information, KS-RegNet performs a deformable image registration between a fully-sampled prior image and on-board images acquired from highly-under-sampled k-space data, to generate high-quality on-board images for real-time motion tracking. Methods. KS-RegNet is an end-to-end, unsupervised network consisting of an input data generation block, a subsequent U-Net core block, and following operations to compute data fidelity and regularization losses. The input data involved a fully-sampled, complex-valued prior image, and the k-space data of an on-board, real-time MR image (MRI). From the k-space data, under-sampled real-time MRI was reconstructed by the data generation block to input into the U-Net core. In addition, to train the U-Net core to learn the under-sampling artifacts, the k-space data of the prior image was intentionally under-sampled using the same readout trajectory as the real-time MRI, and reconstructed to serve an additional input. The U-Net core predicted a deformation vector field that deforms the prior MRI to on-board real-time MRI. To avoid adverse effects of quantifying image similarity on the artifacts-ridden images, the data fidelity loss of deformation was evaluated directly in k-space. Results. Compared with Elastix and other deep learning network architectures, KS-RegNet demonstrated better and more stable performance. The average (±s.d.) DICE coefficients of KS-RegNet on a cardiac dataset for the 5- , 9- , and 13-spoke k-space acquisitions were 0.884 ± 0.025, 0.889 ± 0.024, and 0.894 ± 0.022, respectively; and the corresponding average (±s.d.) center-of-mass errors (COMEs) were 1.21 ± 1.09, 1.29 ± 1.22, and 1.01 ± 0.86 mm, respectively. KS-RegNet also provided the best performance on an abdominal dataset. Conclusion. KS-RegNet allows real-time MRI generation with sub-second latency. It enables potential real-time MR-guided soft tissue tracking, tumor localization, and radiotherapy plan adaptation.

Ultrahigh Resolution fMRI at 7T Using Radial-Cartesian TURBINE Sampling

Purpose

We investigate the use of TURBINE, a 3D radial‐Cartesian acquisition scheme in which EPI planes are rotated about the phase‐encoding axis to acquire a cylindrical k‐space for high‐fidelity ultrahigh isotropic resolution fMRI at 7 Tesla with minimal distortion and blurring.

Methods

An improved, completely self‐navigated version of the TURBINE sampling scheme was designed for fMRI at 7 Telsa. To demonstrate the image quality and spatial specificity of the acquisition, thin‐slab visual and motor BOLD fMRI at 0.67 mm isotropic resolution (16 mm slab, TRvol = 2.32 s), and 0.8 × 0.8 × 2.0 mm (whole‐brain, TRvol = 2.4 s) data were acquired. To prioritize the high spatial fidelity, we employed a temporally regularized reconstruction to improve sensitivity without any spatial bias.

Results

TURBINE images provide high structural fidelity with almost no distortion, dropout, or T₂* blurring for the thin‐slab acquisitions compared to conventional 3D EPI owing to the radial sampling in‐plane and the short echo train used. This results in activation that can be localized to pre‐ and postcentral gyri in a motor task, for example, with excellent correspondence to brain structure measured by a T₁‐MPRAGE. The benefits of TURBINE (low distortion, dropout, blurring) are reduced for the whole‐brain acquisition due to the longer EPI train. We demonstrate robust BOLD activation at 0.67 mm isotropic resolution (thin‐slab) and also anisotropic 0.8 × 0.8 × 2.0 mm (whole‐brain) acquisitions.

Conclusion

TURBINE is a promising acquisition approach for high‐resolution, minimally distorted fMRI at 7 Tesla and could be particularly useful for fMRI in areas of high B₀ inhomogeneity.

Click here for author‐reader discussions

Cramér-Rao bound-informed training of neural networks for quantitative MRI

PURPOSE

To improve the performance of neural networks for parameter estimation in quantitative MRI, in particular when the noise propagation varies throughout the space of biophysical parameters.

THEORY AND METHODS

A theoretically well-founded loss function is proposed that normalizes the squared error of each estimate with respective Cramér-Rao bound (CRB)-a theoretical lower bound for the variance of an unbiased estimator. This avoids a dominance of hard-to-estimate parameters and areas in parameter space, which are often of little interest. The normalization with corresponding CRB balances the large errors of fundamentally more noisy estimates and the small errors of fundamentally less noisy estimates, allowing the network to better learn to estimate the latter. Further, proposed loss function provides an absolute evaluation metric for performance: A network has an average loss of 1 if it is a maximally efficient unbiased estimator, which can be considered the ideal performance. The performance gain with proposed loss function is demonstrated at the example of an eight-parameter magnetization transfer model that is fitted to phantom and in vivo data.

RESULTS

Networks trained with proposed loss function perform close to optimal, that is, their loss converges to approximately 1, and their performance is superior to networks trained with the standard mean-squared error (MSE). The proposed loss function reduces the bias of the estimates compared to the MSE loss, and improves the match of the noise variance to the CRB. This performance gain translates to in vivo maps that align better with the literature.

CONCLUSION

Normalizing the squared error with the CRB during the training of neural networks improves their performance in estimating biophysical parameters.

Cardio-respiratory motion-corrected 3D cardiac water-fat MRI using model-based image reconstruction

PURPOSE

Myocardial fat infiltrations are associated with a range of cardiomyopathies. The purpose of this study was to perform cardio-respiratory motion-correction for model-based water-fat separation to image fatty infiltrations of the heart in a free-breathing, non-cardiac-triggered high-resolution 3D MRI acquisition.

METHODS

Data were acquired in nine patients using a free-breathing, non-cardiac-triggered high-resolution 3D Dixon gradient-echo sequence and radial phase encoding trajectory. Motion correction was combined with a model-based water-fat reconstruction approach. Respiratory and cardiac motion models were estimated using a dual-mode registration algorithm incorporating both motion-resolved water and fat information. Qualitative comparisons of fat structures were made between 2D clinical routine reference scans and reformatted 3D motion-corrected images. To evaluate the effect of motion correction the local sharpness of epicardial fat structures was analyzed for motion-averaged and motion-corrected fat images.

RESULTS

The reformatted 3D motion-corrected reconstructions yielded qualitatively comparable fat structures and fat structure sharpness in the heart as the standard 2D breath-hold. Respiratory motion correction improved the local sharpness on average by 32% ± 24% with maximum improvements of 81% and cardiac motion correction increased the sharpness further by another 15% ± 11% with maximum increases of 31%. One patient showed a fat infiltration in the myocardium and cardio-respiratory motion correction was able to improve its visualization in 3D.

CONCLUSION

The 3D water-fat separated cardiac images were acquired during free-breathing and in a clinically feasible and predictable scan time. Compared to a motion-averaged reconstruction an increase in sharpness of fat structures by 51% ± 27% using the presented motion correction approach was observed for nine patients.

High-resolution CINE imaging of active guided knee motion using continuously acquired golden-angle radial MRI and rotary sensor information

To explore and extend on dynamic imaging of joint motion, an MRI-safe device guiding knee motion with an attached rotary encoder was used in MRI measurements of multiple knee flexion-extension cycles using radial gradient echo imaging with the golden-angle as azimuthal angle increment. Reproducibility of knee motion was investigated. Real-time and CINE mode anatomical images were reconstructed for different knee flexion angles by synchronizing the encoder information with the MRI data, and performing flexion angle selective gating across multiple motion cycles. When investigating the influence of the rotation angle window width on reconstructed CINE images, it was found that angle windows between 0.5° and 3° exhibited acceptable image sharpness without suffering from significant motion-induced blurring. Furthermore, due to flexible retrospective image reconstruction afforded by the radial golden-angle imaging, the number of motion cycles included in the reconstruction could be retrospectively reduced to investigate the corresponding influence of acquisition time on image quality. Finally, motion reproducibility between motion cycles and accuracy of the flexion angle selective gating were sufficient to acquire whole-knee 3D dynamic imaging with a retrospectively gated 3D cone UTE sequence.

A Self-Supervised Deep Learning Reconstruction for Shortening the Breathhold and Acquisition Window in Cardiac Magnetic Resonance Fingerprinting

The aim of this study is to shorten the breathhold and diastolic acquisition window in cardiac magnetic resonance fingerprinting (MRF) for simultaneous T1, T2, and proton spin density (M0) mapping to improve scan efficiency and reduce motion artifacts. To this end, a novel reconstruction was developed that combines low-rank subspace modeling with a deep image prior, termed DIP-MRF. A system of neural networks is used to generate spatial basis images and quantitative tissue property maps, with training performed using only the undersampled k-space measurements from the current scan. This approach avoids difficulties with obtaining in vivo MRF training data, as training is performed de novo for each acquisition. Calculation of the forward model during training is accelerated by using GRAPPA operator gridding to shift spiral k-space data to Cartesian grid points, and by using a neural network to rapidly generate fingerprints in place of a Bloch equation simulation. DIP-MRF was evaluated in simulations and at 1.5 T in a standardized phantom, 18 healthy subjects, and 10 patients with suspected cardiomyopathy. In addition to conventional mapping, two cardiac MRF sequences were acquired, one with a 15-heartbeat(HB) breathhold and 254 ms acquisition window, and one with a 5HB breathhold and 150 ms acquisition window. In simulations, DIP-MRF yielded decreased nRMSE compared to dictionary matching and a sparse and locally low rank (SLLR-MRF) reconstruction. Strong correlation (R2 &gt; 0.999) with T1 and T2 reference values was observed in the phantom using the 5HB/150 ms scan with DIP-MRF. DIP-MRF provided better suppression of noise and aliasing artifacts in vivo, especially for the 5HB/150 ms scan, and lower intersubject and intrasubject variability compared to dictionary matching and SLLR-MRF. Furthermore, it yielded a better agreement between myocardial T1 and T2 from 15HB/254 ms and 5HB/150 ms MRF scans, with a bias of −9 ms for T1 and 2 ms for T2. In summary, this study introduces an extension of the deep image prior framework for cardiac MRF tissue property mapping, which does not require pre-training with in vivo scans, and has the potential to reduce motion artifacts by enabling a shortened breathhold and acquisition window.

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.

Event-recurring multiband SWIFT functional MRI with 200-ms temporal resolution during deep brain stimulation and isoflurane-induced burst suppression in rat

PURPOSE

To develop a high temporal resolution functional MRI method for tracking repeating events in the brain.

METHODS

We developed a novel functional MRI method using multiband sweep imaging with Fourier transformation (SWIFT), termed event-recurring SWIFT (EVER-SWIFT). The method is able to image similar repeating events with subsecond temporal resolution. Here, we demonstrate the use of EVER-SWIFT for detecting functional MRI responses during deep brain stimulation of the medial septal nucleus and during spontaneous isoflurane-induced burst suppression in the rat brain at 9.4 T with 200-ms temporal resolution.

RESULTS

The EVER-SWIFT approach showed that the shapes and time-to-peak values of the response curves to deep brain stimulation significantly differed between downstream brain regions connected to the medial septal nucleus, resembling findings obtained with traditional 2-second temporal resolution. In contrast, EVER-SWIFT allowed for detailed temporal measurement of a spontaneous isoflurane-induced bursting activity pattern, which was not achieved with traditional temporal resolution.

CONCLUSION

The EVER-SWIFT technique enables subsecond 3D imaging of both stimulated and spontaneously recurring brain activities, and thus holds great potential for studying the mechanisms of neuromodulation and spontaneous brain activity.

Embedded Quantitative MRI T1ρ Mapping Using Non-Linear Primal-Dual Proximal Splitting

Quantitative MRI (qMRI) methods allow reducing the subjectivity of clinical MRI by providing numerical values on which diagnostic assessment or predictions of tissue properties can be based. However, qMRI measurements typically take more time than anatomical imaging due to requiring multiple measurements with varying contrasts for, e.g., relaxation time mapping. To reduce the scanning time, undersampled data may be combined with compressed sensing (CS) reconstruction techniques. Typical CS reconstructions first reconstruct a complex-valued set of images corresponding to the varying contrasts, followed by a non-linear signal model fit to obtain the parameter maps. We propose a direct, embedded reconstruction method for T1ρ mapping. The proposed method capitalizes on a known signal model to directly reconstruct the desired parameter map using a non-linear optimization model. The proposed reconstruction method also allows directly regularizing the parameter map of interest and greatly reduces the number of unknowns in the reconstruction, which are key factors in the performance of the reconstruction method. We test the proposed model using simulated radially sampled data from a 2D phantom and 2D cartesian ex vivo measurements of a mouse kidney specimen. We compare the embedded reconstruction model to two CS reconstruction models and in the cartesian test case also the direct inverse fast Fourier transform. The T1ρ RMSE of the embedded reconstructions was reduced by 37–76% compared to the CS reconstructions when using undersampled simulated data with the reduction growing with larger acceleration factors. The proposed, embedded model outperformed the reference methods on the experimental test case as well, especially providing robustness with higher acceleration factors.

Simultaneous 3D acquisition of 1 H MRF and 23 Na MRI

PURPOSE

To develop a 3D MR technique to simultaneously acquire proton multiparametric maps (T₁ , T₂ , and proton density) and sodium density weighted images over the whole brain.

METHODS

We implemented a 3D stack-of-stars MR pulse sequence which consists of interleaved proton (¹ H) and sodium (²³ Na) excitations, tailored slice encoding gradients that can encode the same slice for both nuclei, and simultaneous readout with different radial trajectories (¹ H, full-radial; ²³ Na, center-out radial). The receive chain of our 7T scanner was modified to enable simultaneous acquisition of ¹ H and ²³ Na signal. A heuristically optimized flip angle train was implemented for proton MR fingerprinting (MRF). The SNR and the accuracy of proton T₁ and T₂ were evaluated in phantoms. Finally, in vivo application of the method was demonstrated in five healthy subjects.

RESULTS

The SNR for the simultaneous measurement was almost identical to that for the single-nucleus measurements (<2% change). The proton T₁ and T₂ maps remained similar to the results from a reference 2D MRF technique (normalized RMS error in T₁ ≈ 4.2% and T₂ ≈ 11.3%). Measurements in healthy subjects corroborated these results and demonstrated the feasibility of our method for in vivo application. The in vivo T₁ values measured using our method were lower than the results measured by other conventional techniques.

CONCLUSIONS

With the 3D simultaneous implementation, we were able to acquire sodium and proton density weighted images in addition to proton T₁ , T₂ , and B1+ from ¹ H MRF that covers the whole brain volume within 21 min.