Improved magnetic resonance fingerprinting reconstruction with low-rank and subspace modeling

Efficient pulse sequence design framework for high-dimensional MR fingerprinting scans using systematic error index

PURPOSE

For effective optimization of MR fingerprinting (MRF) pulse sequences, estimating and minimizing errors from actual scan conditions are crucial. Although virtual-scan simulations offer an approximation to these errors, their computational demands become expensive for high-dimensional MRF frameworks, where interactions between more than two tissue properties are considered. This complexity makes sequence optimization impractical. We introduce a new mathematical model, the systematic error index (SEI), to address the scalability challenges for high-dimensional MRF sequence design.

METHODS

By eliminating the need to perform dictionary matching, the SEI model approximates quantification errors with low computational costs. The SEI model was validated in comparison with virtual-scan simulations. The SEI model was further applied to optimize three high-dimensional MRF sequences that quantify two to four tissue properties. The optimized scans were examined in simulations and healthy subjects.

RESULTS

The proposed SEI model closely approximated the virtual-scan simulation outcomes while achieving hundred- to thousand-times acceleration in the computational speed. In both simulation and in vivo experiments, the optimized MRF sequences yield higher measurement accuracy with fewer undersampling artifacts at shorter scan times than the heuristically designed sequences.

CONCLUSION

We developed an efficient method for estimating real-world errors in MRF scans with high computational efficiency. Our results illustrate that the SEI model could approximate errors both qualitatively and quantitatively. We also proved the practicality of the SEI model of optimizing sequences for high-dimensional MRF frameworks with manageable computational power. The optimized high-dimensional MRF scans exhibited enhanced robustness against undersampling and system imperfections with faster scan times.

Quantification of water exchange across the blood-brain barrier using noncontrast MR fingerprinting

PURPOSE

A method is proposed to quantify cerebral blood volume (v b $$ {v}_b $$ ) and intravascular water residence time (τ b $$ {\tau}_b $$ ) using MR fingerprinting (MRF), applied using a spoiled gradient echo sequence without the need for contrast agent.

METHODS

An in silico study optimized an acquisition protocol to maximize the sensitivity of the measurement tov b $$ {v}_b $$ andτ b $$ {\tau}_b $$ changes. Its accuracy in the presence of variations inT 1 , t $$ {\mathrm{T}}_{1,t} $$ ,T 1 , b $$ {\mathrm{T}}_{1,b} $$ , andB 1 $$ {\mathrm{B}}_1 $$ was evaluated. The optimized protocol (scan time of 19 min) was then tested in a exploratory healthy volunteer study (10 volunteers, mean age 24± $$ \pm $$ 3, six males) at 3 T with a repeat scan taken after repositioning to allow estimation of repeatability.

RESULTS

Simulations show that assuming literature values forT 1 , b $$ {\mathrm{T}}_{1,b} $$ andT 1 , t $$ {\mathrm{T}}_{1,t} $$ , no variation inB 1 $$ {\mathrm{B}}_1 $$ , while fitting onlyv b $$ {v}_b $$ andτ b $$ {\tau}_b $$ , leads to large errors in quantification ofv b $$ {v}_b $$ andτ b $$ {\tau}_b $$ , regardless of noise levels. However, simulations also show that matchingT 1 , t $$ {\mathrm{T}}_{1,t} $$ ,T 1 , b $$ {\mathrm{T}}_{1,b} $$ ,B 1 + $$ {\mathrm{B}}_1^{+} $$ ,v b $$ {v}_b $$ andτ b $$ {\tau}_b $$ , simultaneously is feasible at clinically achievable noise levels. Across the healthy volunteers, all parameter quantifications fell within the expected literature range. In addition, the maps show good agreement between hemispheres suggesting physiologically relevant information is being extracted. Expected differences between white and gray matterT 1 , t $$ {\mathrm{T}}_{1,t} $$ (p < 0.0001) andv b $$ {v}_b $$ (p < 0.0001) are observed,T 1 , b $$ {\mathrm{T}}_{1,b} $$ andτ b $$ {\tau}_b $$ show no significant differences, p = 0.4 and p = 0.6, respectively. Moderate to excellent repeatability was seen between repeat scans: mean intra-class correlation coefficient ofT 1 , t : 0 . 91 $$ {\mathrm{T}}_{1,t}:0.91 $$ ,T 1 , b : 0 . 58 $$ {\mathrm{T}}_{1,b}:0.58 $$ ,v b : 0 . 90 $$ {v}_b:0.90 $$ , andτ b : 0 . 96 $$ {\tau}_b:0.96 $$ .

CONCLUSION

We demonstrate that regional simultaneous quantification ofv b $$ {v}_b $$ ,τ b $$ {\tau}_b $$ ,T 1 , b , T 1 , t $$ {\mathrm{T}}_{1,b},{T}_{1,t} $$ , andB 1 + $$ {\mathrm{B}}_1^{+} $$ using MRF is feasible in vivo.

Bias-reduced neural networks for parameter estimation in quantitative MRI

PURPOSE

To develop neural network (NN)-based quantitative MRI parameter estimators with minimal bias and a variance close to the Cramér-Rao bound.

THEORY AND METHODS

We generalize the mean squared error loss to control the bias and variance of the NN's estimates, which involves averaging over multiple noise realizations of the same measurements during training. Bias and variance properties of the resulting NNs are studied for two neuroimaging applications.

RESULTS

In simulations, the proposed strategy reduces the estimates' bias throughout parameter space and achieves a variance close to the Cramér-Rao bound. In vivo, we observe good concordance between parameter maps estimated with the proposed NNs and traditional estimators, such as nonlinear least-squares fitting, while state-of-the-art NNs show larger deviations.

CONCLUSION

The proposed NNs have greatly reduced bias compared to those trained using the mean squared error and offer significantly improved computational efficiency over traditional estimators with comparable or better accuracy.

High-resolution myelin-water fraction and quantitative relaxation mapping using 3D ViSTa-MR fingerprinting

PURPOSE

This study aims to develop a high-resolution whole-brain multi-parametric quantitative MRI approach for simultaneous mapping of myelin-water fraction (MWF), T1, T2, and proton-density (PD), all within a clinically feasible scan time.

METHODS

We developed 3D visualization of short transverse relaxation time component (ViSTa)-MRF, which combined ViSTa technique with MR fingerprinting (MRF), to achieve high-fidelity whole-brain MWF and T1/T2/PD mapping on a clinical 3T scanner. To achieve fast acquisition and memory-efficient reconstruction, the ViSTa-MRF sequence leverages an optimized 3D tiny-golden-angle-shuffling spiral-projection acquisition and joint spatial-temporal subspace reconstruction with optimized preconditioning algorithm. With the proposed ViSTa-MRF approach, high-fidelity direct MWF mapping was achieved without a need for multicompartment fitting that could introduce bias and/or noise from additional assumptions or priors.

RESULTS

The in vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide fast multi-parametric mapping with high SNR and good quality. The in vivo results of 1 mm- and 0.66 mm-isotropic resolution datasets indicate that the MWF values measured by the proposed method are consistent with standard ViSTa results that are 30× slower with lower SNR. Furthermore, we applied the proposed method to enable 5-min whole-brain 1 mm-iso assessment of MWF and T1/T2/PD mappings for infant brain development and for post-mortem brain samples.

CONCLUSIONS

In this work, we have developed a 3D ViSTa-MRF technique that enables the acquisition of whole-brain MWF, quantitative T1, T2, and PD maps at 1 and 0.66 mm isotropic resolution in 5 and 15 min, respectively. This advancement allows for quantitative investigations of myelination changes in the brain.

Zero-DeepSub: Zero-shot deep subspace reconstruction for rapid multiparametric quantitative MRI using 3D-QALAS

PURPOSE

To develop and evaluate methods for (1) reconstructing 3D-quantification using an interleaved Look-Locker acquisition sequence with T2 preparation pulse (3D-QALAS) time-series images using a low-rank subspace method, which enables accurate and rapid T1 and T2 mapping, and (2) improving the fidelity of subspace QALAS by combining scan-specific deep-learning-based reconstruction and subspace modeling.

THEORY AND METHODS

A low-rank subspace method for 3D-QALAS (i.e., subspace QALAS) and zero-shot deep-learning subspace method (i.e., Zero-DeepSub) were proposed for rapid and high fidelity T1 and T2 mapping and time-resolved imaging using 3D-QALAS. Using an ISMRM/NIST system phantom, the accuracy and reproducibility of the T1 and T2 maps estimated using the proposed methods were evaluated by comparing them with reference techniques. The reconstruction performance of the proposed subspace QALAS using Zero-DeepSub was evaluated in vivo and compared with conventional QALAS at high reduction factors of up to nine-fold.

RESULTS

Phantom experiments showed that subspace QALAS had good linearity with respect to the reference methods while reducing biases and improving precision compared to conventional QALAS, especially for T2 maps. Moreover, in vivo results demonstrated that subspace QALAS had better g-factor maps and could reduce voxel blurring, noise, and artifacts compared to conventional QALAS and showed robust performance at up to nine-fold acceleration with Zero-DeepSub, which enabled whole-brain T1, T2, and PD mapping at 1 mm isotropic resolution within 2 min of scan time.

CONCLUSION

The proposed subspace QALAS along with Zero-DeepSub enabled high fidelity and rapid whole-brain multiparametric quantification and time-resolved imaging.

Motion Robust MR Fingerprinting Scan to Image Neonates With Prenatal Opioid Exposure

PURPOSE

To explore whether MR fingerprinting (MRF) scans provide motion-robust and quantitative brain tissue measurements for non-sedated infants with prenatal opioid exposure (POE).

STUDY TYPE

Prospective.

POPULATION

13 infants with POE (3 male; 12 newborns (age 7-65 days) and 1 infant aged 9-months).

FIELD STRENGTH/SEQUENCE

3T, 3D T1-weighted MPRAGE, 3D T2-weighted TSE and MRF sequences.

ASSESSMENT

The image quality of MRF and MRI was assessed in a fully crossed, multiple-reader, multiple-case study. Sixteen image quality features in three types-image artifacts, structure and myelination visualization-were ranked by four neuroradiologists (8, 7, 5, and 8 years of experience respectively), using a 3-point scale. MRF T1 and T2 values in 8 white matter brain regions were compared between babies younger than 1 month and babies between 1 and 2 months.

STATISTICAL TESTS

Generalized estimating equations model to test the significance of differences of regional T1 and T2 values of babies under 1 month and those older. MRI and MRF image quality was assessed using Gwet's second order auto-correlation coefficient (AC2) with confidence levels. The Cochran-Mantel-Haenszel test was used to assess the difference in proportions between MRF and MRI for all features and stratified by the type of features. A P value <0.05 was considered statistically significant.

RESULTS

The MRF of two infants were excluded in T1 and T2 value analysis due to severe motion artifact but were included in the image quality assessment. In infants under 1 month of age (N = 6), the T1 and T2 values were significantly higher compared to those between 1 and 2 months of age (N = 4). MRF images showed significantly higher image quality ratings in all three feature types compared to MRI images.

CONCLUSIONS

MR Fingerprinting scans have potential to be a motion-robust and efficient method for nonsedated infants.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1.

Self-calibrated subspace reconstruction for multidimensional MR fingerprinting for simultaneous relaxation and diffusion quantification

PURPOSE

To propose a new reconstruction method for multidimensional MR fingerprinting (mdMRF) to address shading artifacts caused by physiological motion-induced measurement errors without navigating or gating.

METHODS

The proposed method comprises two procedures: self-calibration and subspace reconstruction. The first procedure (self-calibration) applies temporally local matrix completion to reconstruct low-resolution images from a subset of under-sampled data extracted from the k-space center. The second procedure (subspace reconstruction) utilizes temporally global subspace reconstruction with pre-estimated temporal subspace from low-resolution images to reconstruct aliasing-free, high-resolution, and time-resolved images. After reconstruction, a customized outlier detection algorithm was employed to automatically detect and remove images corrupted by measurement errors. Feasibility, robustness, and scan efficiency were evaluated through in vivo human brain imaging experiments.

RESULTS

The proposed method successfully reconstructed aliasing-free, high-resolution, and time-resolved images, where the measurement errors were accurately represented. The corrupted images were automatically and robustly detected and removed. Artifact-free T1, T2, and ADC maps were generated simultaneously. The proposed reconstruction method demonstrated robustness across different scanners, parameter settings, and subjects. A high scan efficiency of less than 20 s per slice has been achieved.

CONCLUSION

The proposed reconstruction method can effectively alleviate shading artifacts caused by physiological motion-induced measurement errors. It enables simultaneous and artifact-free quantification of T1, T2, and ADC using mdMRF scans without prospective gating, with robustness and high scan efficiency.

High-fidelity intravoxel incoherent motion parameter mapping using locally low-rank and subspace modeling

PURPOSE

Intravoxel incoherent motion (IVIM) is a quantitative magnetic resonance imaging (MRI) method used to quantify perfusion properties of tissue non-invasively without contrast. However, clinical applications are limited by unreliable parameter estimates, particularly for the perfusion fraction (f) and pseudodiffusion coefficient (D*). This study aims to develop a high-fidelity reconstruction for reliable estimation of IVIM parameters. The proposed method is versatile and amenable to various acquisition schemes and fitting methods.

METHODS

To address current challenges with IVIM, we adapted several advanced reconstruction techniques. We used a low-rank approximation of IVIM images and temporal subspace modeling to constrain the magnetization dynamics of the bi-exponential diffusion signal decay. In addition, motion-induced phase variations were corrected between diffusion directions and b-values, facilitating the use of high SNR real-valued diffusion data. The proposed method was evaluated in simulations and in vivo brain acquisitions in six healthy subjects and six individuals with a history of SARS-CoV-2 infection and compared with the conventionally reconstructed magnitude data. Following reconstruction, IVIM parameters were estimated voxel-wise.

RESULTS

Our proposed method reduced noise contamination in simulations, resulting in a 60%, 58.9%, and 83.9% reduction in the NRMSE for D, f, and D*, respectively, compared to the conventional reconstruction. In vivo, anisotropic properties of D, f, and D* were preserved with the proposed method, highlighting microvascular differences in gray matter between individuals with a history of COVID-19 and those without (p = 0.0210), which wasn't observed with the conventional reconstruction.

CONCLUSION

The proposed method yielded a more reliable estimation of IVIM parameters with less noise than the conventional reconstruction. Further, the proposed method preserved anisotropic properties of IVIM parameter estimates and demonstrated differences in microvascular perfusion in COVID-affected subjects, which weren't observed with conventional reconstruction methods.

Bias-Reduced Neural Networks for Parameter Estimation in Quantitative MRI

Purpose

To develop neural network (NN)-based quantitative MRI parameter estimators with minimal bias and a variance close to the Cramér-Rao bound.

Theory and Methods

We generalize the mean squared error loss to control the bias and variance of the NN's estimates, which involves averaging over multiple noise realizations of the same measurements during training. Bias and variance properties of the resulting NNs are studied for two neuroimaging applications.

Results

In simulations, the proposed strategy reduces the estimates' bias throughout parameter space and achieves a variance close to the Cramér-Rao bound. In vivo, we observe good concordance between parameter maps estimated with the proposed NNs and traditional estimators, such as non-linear least-squares fitting, while state-of-the-art NNs show larger deviations.

Conclusion

The proposed NNs have greatly reduced bias compared to those trained using the mean squared error and offer significantly improved computational efficiency over traditional estimators with comparable or better accuracy.

Rapid quantitative magnetization transfer imaging: Utilizing the hybrid state and the generalized Bloch model

PURPOSE

To explore efficient encoding schemes for quantitative magnetization transfer (qMT) imaging with few constraints on model parameters.

THEORY AND METHODS

We combine two recently proposed models in a Bloch-McConnell equation: the dynamics of the free spin pool are confined to the hybrid state, and the dynamics of the semi-solid spin pool are described by the generalized Bloch model. We numerically optimize the flip angles and durations of a train of radio frequency pulses to enhance the encoding of three qMT parameters while accounting for all eight parameters of the two-pool model. We sparsely sample each time frame along this spin dynamics with a three-dimensional radial koosh-ball trajectory, reconstruct the data with subspace modeling, and fit the qMT model with a neural network for computational efficiency.

RESULTS

We extracted qMT parameter maps of the whole brain with an effective resolution of 1.24 mm from a 12.6-min scan. In lesions of multiple sclerosis subjects, we observe a decreased size of the semi-solid spin pool and longer relaxation times, consistent with previous reports.

CONCLUSION

The encoding power of the hybrid state, combined with regularized image reconstruction, and the accuracy of the generalized Bloch model provide an excellent basis for efficient quantitative magnetization transfer imaging with few constraints on model parameters.

Emerging Trends in Magnetic Resonance Fingerprinting for Quantitative Biomedical Imaging Applications: A Review

Magnetic resonance imaging (MRI) stands as a vital medical imaging technique, renowned for its ability to offer high-resolution images of the human body with remarkable soft-tissue contrast. This enables healthcare professionals to gain valuable insights into various aspects of the human body, including morphology, structural integrity, and physiological processes. Quantitative imaging provides compositional measurements of the human body, but, currently, either it takes a long scan time or is limited to low spatial resolutions. Undersampled k-space data acquisitions have significantly helped to reduce MRI scan time, while compressed sensing (CS) and deep learning (DL) reconstructions have mitigated the associated undersampling artifacts. Alternatively, magnetic resonance fingerprinting (MRF) provides an efficient and versatile framework to acquire and quantify multiple tissue properties simultaneously from a single fast MRI scan. The MRF framework involves four key aspects: (1) pulse sequence design; (2) rapid (undersampled) data acquisition; (3) encoding of tissue properties in MR signal evolutions or fingerprints; and (4) simultaneous recovery of multiple quantitative spatial maps. This paper provides an extensive literature review of the MRF framework, addressing the trends associated with these four key aspects. There are specific challenges in MRF for all ranges of magnetic field strengths and all body parts, which can present opportunities for further investigation. We aim to review the best practices in each key aspect of MRF, as well as for different applications, such as cardiac, brain, and musculoskeletal imaging, among others. A comprehensive review of these applications will enable us to assess future trends and their implications for the translation of MRF into these biomedical imaging applications.

Cardiac MR Fingerprinting: Overview, Technical Developments, and Applications

Cardiovascular magnetic resonance (CMR) is an established imaging modality with proven utility in assessing cardiovascular diseases. The ability of CMR to characterize myocardial tissue using T1 - and T2 -weighted imaging, parametric mapping, and late gadolinium enhancement has allowed for the non-invasive identification of specific pathologies not previously possible with modalities like echocardiography. However, CMR examinations are lengthy and technically complex, requiring multiple pulse sequences and different anatomical planes to comprehensively assess myocardial structure, function, and tissue composition. To increase the overall impact of this modality, there is a need to simplify and shorten CMR exams to improve access and efficiency, while also providing reproducible quantitative measurements. Multiparametric MRI techniques that measure multiple tissue properties offer one potential solution to this problem. This review provides an in-depth look at one such multiparametric approach, cardiac magnetic resonance fingerprinting (MRF). The article is structured as follows. First, a brief review of single-parametric and (non-Fingerprinting) multiparametric CMR mapping techniques is presented. Second, a general overview of cardiac MRF is provided covering pulse sequence implementation, dictionary generation, fast k-space sampling methods, and pattern recognition. Third, recent technical advances in cardiac MRF are covered spanning a variety of topics, including simultaneous multislice and 3D sampling, motion correction algorithms, cine MRF, synthetic multicontrast imaging, extensions to measure additional clinically important tissue properties (proton density fat fraction, T2 *, and T1ρ ), and deep learning methods for image reconstruction and parameter estimation. The last section will discuss potential clinical applications, concluding with a perspective on how multiparametric techniques like MRF may enable streamlined CMR protocols. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 1.

High-resolution myelin-water fraction and quantitative relaxation mapping using 3D ViSTa-MR fingerprinting

Purpose

This study aims to develop a high-resolution whole-brain multi-parametric quantitative MRI approach for simultaneous mapping of myelin-water fraction (MWF), T1, T2, and proton-density (PD), all within a clinically feasible scan time.

Methods

We developed 3D ViSTa-MRF, which combined Visualization of Short Transverse relaxation time component (ViSTa) technique with MR Fingerprinting (MRF), to achieve high-fidelity whole-brain MWF and T1/T2/PD mapping on a clinical 3T scanner. To achieve fast acquisition and memory-efficient reconstruction, the ViSTa-MRF sequence leverages an optimized 3D tiny-golden-angle-shuffling spiral-projection acquisition and joint spatial-temporal subspace reconstruction with optimized preconditioning algorithm. With the proposed ViSTa-MRF approach, high-fidelity direct MWF mapping was achieved without a need for multi-compartment fitting that could introduce bias and/or noise from additional assumptions or priors.

Results

The in-vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide fast multi-parametric mapping with high SNR and good quality. The in-vivo results of 1mm- and 0.66mm-iso datasets indicate that the MWF values measured by the proposed method are consistent with standard ViSTa results that are 30x slower with lower SNR. Furthermore, we applied the proposed method to enable 5-minute whole-brain 1mm-iso assessment of MWF and T1/T2/PD mappings for infant brain development and for post-mortem brain samples.

Conclusions

In this work, we have developed a 3D ViSTa-MRF technique that enables the acquisition of whole-brain MWF, quantitative T1, T2, and PD maps at 1mm and 0.66mm isotropic resolution in 5 and 15 minutes, respectively. This advancement allows for quantitative investigations of myelination changes in the brain.

Learned Tensor Low-CP-Rank and Bloch Response Manifold Priors for Non-Cartesian MRF Reconstruction

Magnetic resonance fingerprinting (MRF) can rapidly perform simultaneous imaging of multiple tissue parameters. However, the rapid acquisition schemes used in MRF inevitably introduce aliasing artifacts in the recovered tissue fingerprints, reducing the accuracy of the predicted parameter maps. Current regularized reconstruction methods are based on iterative procedures which are usually time-consuming. In addition, most of the current deep learning-based methods for MRF often lack interpretability owing to the black-box nature, and most deep learning-based methods are not applicable for non-Cartesian scenarios, which limits the practical applications. In this paper, we propose a joint reconstruction model incorporating MRF-physics prior and the data correlation constraint for non-Cartesian MRF reconstruction. To avoid time-consuming iterative procedures, we unroll the reconstruction model into a deep neural network. Specifically, we propose a learned CANDECOMP/PARAFAC (CP) decomposition module to exploit the tensor low-rank priors of high-dimensional MRF data, which avoids computationally burdensome singular value decomposition. Inspired by the MRF-physics, we also propose a Bloch response manifold module to learn the mapping between reconstructed MRF data and the multiple parameter maps. Numerical experiments show that the proposed network can reconstruct high-quality MRF data and multiple parameter maps within significantly reduced computational time.

Cramér-Rao Bound Optimized Subspace Reconstruction in Quantitative MRI

We extend the traditional framework for estimating subspace bases that maximize the preserved signal energy to additionally preserve the Cramér-Rao bound (CRB) of the biophysical parameters and, ultimately, improve accuracy and precision in the quantitative maps. To this end, we introduce an approximate compressed CRB based on orthogonalized versions of the signal's derivatives with respect to the model parameters. This approximation permits singular value decomposition (SVD)-based minimization of both the CRB and signal losses during compression. Compared to the traditional SVD approach, the proposed method better preserves the CRB across all biophysical parameters with negligible cost to the preserved signal energy, leading to reduced bias and variance of the parameter estimates in simulation. In vivo, improved accuracy and precision are observed in two quantitative neuroimaging applications, permitting the use of smaller basis sizes in subspace reconstruction and offering significant computational savings.

High-fidelity direct contrast synthesis from magnetic resonance fingerprinting

PURPOSE

This work was aimed at proposing a supervised learning-based method that directly synthesizes contrast-weighted images from the Magnetic Resonance Fingerprinting (MRF) data without performing quantitative mapping and spin-dynamics simulations.

METHODS

To implement our direct contrast synthesis (DCS) method, we deploy a conditional generative adversarial network (GAN) framework with a multi-branch U-Net as the generator and a multilayer CNN (PatchGAN) as the discriminator. We refer to our proposed approach as N-DCSNet. The input MRF data are used to directly synthesize T1-weighted, T2-weighted, and fluid-attenuated inversion recovery (FLAIR) images through supervised training on paired MRF and target spin echo-based contrast-weighted scans. The performance of our proposed method is demonstrated on in vivo MRF scans from healthy volunteers. Quantitative metrics, including normalized root mean square error (nRMSE), peak signal-to-noise ratio (PSNR), structural similarity (SSIM), learned perceptual image patch similarity (LPIPS), and Fréchet inception distance (FID), were used to evaluate the performance of the proposed method and compare it with others.

RESULTS

In-vivo experiments demonstrated excellent image quality with respect to that of simulation-based contrast synthesis and previous DCS methods, both visually and according to quantitative metrics. We also demonstrate cases in which our trained model is able to mitigate the in-flow and spiral off-resonance artifacts typically seen in MRF reconstructions, and thus more faithfully represent conventional spin echo-based contrast-weighted images.

CONCLUSION

We present N-DCSNet to directly synthesize high-fidelity multicontrast MR images from a single MRF acquisition. This method can significantly decrease examination time. By directly training a network to generate contrast-weighted images, our method does not require any model-based simulation and therefore can avoid reconstruction errors due to dictionary matching and contrast simulation (code available at:https://github.com/mikgroup/DCSNet).

Multiparametric MRI: From Simultaneous Rapid Acquisition Methods and Analysis Techniques Using Scoring, Machine Learning, Radiomics, and Deep Learning to the Generation of Novel Metrics

ABSTRACT

With the recent advancements in rapid imaging methods, higher numbers of contrasts and quantitative parameters can be acquired in less and less time. Some acquisition models simultaneously obtain multiparametric images and quantitative maps to reduce scan times and avoid potential issues associated with the registration of different images. Multiparametric magnetic resonance imaging (MRI) has the potential to provide complementary information on a target lesion and thus overcome the limitations of individual techniques. In this review, we introduce methods to acquire multiparametric MRI data in a clinically feasible scan time with a particular focus on simultaneous acquisition techniques, and we discuss how multiparametric MRI data can be analyzed as a whole rather than each parameter separately. Such data analysis approaches include clinical scoring systems, machine learning, radiomics, and deep learning. Other techniques combine multiple images to create new quantitative maps associated with meaningful aspects of human biology. They include the magnetic resonance g-ratio, the inner to the outer diameter of a nerve fiber, and the aerobic glycolytic index, which captures the metabolic status of tumor tissues.

Accelerated Black-Blood Cine MR Imaging with Low-Rank and Sparsity Constraints

Black-blood MRI is a promising imaging technique for assessing vascular diseases (e.g., stroke). Vessel wall dynamic characterization using black-blood cine MRI has been recognized as an effective tool for studying vascular diseases. However, acquiring time-resolved 3D vessel wall images often requires a long acquisition time, which limits its clinical utility. In this work, we develop a new method to achieve rapid, time-resolved 3D black-blood cine MRI. Specifically, the proposed method performs (k, t)-space undersampling to accelerate the volumetric data acquisition process. Moreover, it utilizes an image reconstruction method with low-rank and sparsity constraints to enable high-quality image reconstruction from highly-undersampled data. We validate the performance of the proposed method with 3D in vivo black-blood cine MRI experiments and show representative results to demonstrate the utility of the proposed method.

Low-rank motion correction for accelerated free-breathing first-pass myocardial perfusion imaging

PURPOSE

Develop a novel approach for accelerated 2D free-breathing myocardial perfusion via low-rank motion-corrected (LRMC) reconstructions.

METHODS

Myocardial perfusion imaging requires high spatial and temporal resolution, despite scan time constraints. Here, we incorporate LRMC models into the reconstruction-encoding operator, together with high-dimensionality patch-based regularization, to produce high quality, motion-corrected myocardial perfusion series from free-breathing acquisitions. The proposed framework estimates beat-to-beat nonrigid respiratory (and any other incidental) motion and the dynamic contrast subspace from the actual acquired data, which are then incorporated into the proposed LRMC reconstruction. LRMC was compared with iterative SENSitivity Encoding (SENSE) (itSENSE) and low-rank plus sparse (LpS) reconstruction in 10 patients based on image-quality scoring and ranking by two clinical expert readers.

RESULTS

LRMC achieved significantly improved results relative to itSENSE and LpS in terms of image sharpness, temporal coefficient of variation, and expert reader evaluation. Left ventricle image sharpness was approximately 75%, 79%, and 86% for itSENSE, LpS and LRMC, respectively, indicating improved image sharpness for the proposed approach. Corresponding temporal coefficient of variation results were 23%, 11% and 7%, demonstrating improved temporal fidelity of the perfusion signal with the proposed LRMC. Corresponding clinical expert reader scores (1-5, from poor to excellent image quality) were 3.3, 3.9 and 4.9, demonstrating improved image quality with the proposed LRMC, in agreement with the automated metrics.

CONCLUSION

LRMC produces motion-corrected myocardial perfusion in free-breathing acquisitions with substantially improved image quality when compared with iterative SENSE and LpS reconstructions.

Feasibility of joint mapping of triglyceride saturation and water longitudinal relaxation in a single breath hold applied to high fat-fraction adipose depots in the periclavicular anatomy

INTRODUCTION

Simultaneous mapping of triglyceride (TAG) saturation and tissue water relaxation may improve the characterization of the structure and function of anatomies with significant adipose tissue. While several groups have demonstrated in vivo TAG saturation imaging using MRI, joint mapping of relaxation and TAG saturation is understudied. Such mappings may avoid bias from physiological motion, if they can be done within a single breath-hold, and also account for static and applied magnetic field heterogeneity.

METHODS

We propose a transient-state/MR fingerprinting single breath-hold sequence at 3 T, a low-rank reconstruction, and a parameter estimation pipeline that jointly estimates the number of double bonds (NDB), number of methylene interrupted double bonds (NMIDB), and tissue water T1, while accounting for non-ideal radiofrequency transmit scaling and off-resonance effects. We test the proposed method in simulations, in phantom against MR spectroscopy (MRS), and in vivo regions in and around high fat fraction (FF) periclavicular adipose tissue. Partial volume and multi-peak transverse relaxation effects are explored.

RESULTS

The simulation results demonstrate accurate NDB, NMIDB, and water T1 estimates across a range of NDB, NMIDB, and T1 values. In phantoms, the proposed method's estimates of NDB and NMIDB correlate with those from MR spectroscopy (Pearson correlation ≥0.98), while the water T1 estimates are concordant with a standard phantom. The NDB and NMIDB are sensitive to partial volumes of water, showing increasing bias at FF < 40%. This bias is found to be due to noise and transverse relaxation effects. The in vivo periclavicular adipose tissue has high FF (>90%). The adipose tissue NDB and NMIDB, and muscle T1 estimates are comparable to those reported in the literature.

CONCLUSION

Robust estimation of NDB, NMIDB at high FF and water T1 across a broad range of FFs are feasible using the proposed methods. Further reduction of noise and model bias are needed to employ the proposed technique in low FF anatomies and pathologies.

Motion-resolved real-time 4D flow MRI with low-rank and subspace modeling

PURPOSE

To develop a new motion-resolved real-time four-dimensional (4D) flow MRI method, which enables the quantification and visualization of blood flow velocities with three-directional flow encodings and volumetric coverage without electrocardiogram (ECG) synchronization and respiration control.

METHODS

An integrated imaging method is presented for real-time 4D flow MRI, which encompasses data acquisition, image reconstruction, and postprocessing. The proposed method features a specialized continuous ( k , t ) $$ \left(\mathbf{k},t\right) $$ -space acquisition scheme, which collects two sets of data (i.e., training data and imaging data) in an interleaved manner. By exploiting strong spatiotemporal correlation of 4D flow data, it reconstructs time-series images from highly-undersampled ( k , t ) $$ \left(\mathbf{k},t\right) $$ -space measurements with a low-rank and subspace model. Through data-binning-based postprocessing, it constructs a five-dimensional dataset (i.e., x-y-z-cardiac-respiratory), from which respiration-dependent flow information is further analyzed. The proposed method was evaluated in aortic flow imaging experiments with ten healthy subjects and two patients with atrial fibrillation.

RESULTS

The proposed method achieves 2.4 mm isotropic spatial resolution and 34.4 ms temporal resolution for measuring the blood flow of the aorta. For the healthy subjects, it provides flow measurements in good agreement with those from the conventional 4D flow MRI technique. For the patients with atrial fibrillation, it is able to resolve beat-by-beat pathological flow variations, which cannot be obtained from the conventional technique. The postprocessing further provides respiration-dependent flow information.

CONCLUSION

The proposed method enables high-resolution motion-resolved real-time 4D flow imaging without ECG gating and respiration control. It is able to resolve beat-by-beat blood flow variations as well as respiration-dependent flow information.

Latent signal models: Learning compact representations of signal evolution for improved time-resolved, multi-contrast MRI

PURPOSE

To improve time-resolved reconstructions by training auto-encoders to learn compact representations of Bloch-simulated signal evolution and inserting the decoder into the forward model.

METHODS

Building on model-based nonlinear and linear subspace techniques, we train auto-encoders on dictionaries of simulated signal evolution to learn compact, nonlinear, latent representations. The proposed latent signal model framework inserts the decoder portion of the auto-encoder into the forward model and directly reconstructs the latent representation. Latent signal models essentially serve as a proxy for fast and feasible differentiation through the Bloch equations used to simulate signal. This work performs experiments in the context of T₂ -shuffling, gradient echo EPTI, and MPRAGE-shuffling. We compare how efficiently auto-encoders represent signal evolution in comparison to linear subspaces. Simulation and in vivo experiments then evaluate if reducing degrees of freedom by incorporating our proxy for the Bloch equations, the decoder portion of the auto-encoder, into the forward model improves reconstructions in comparison to subspace constraints.

RESULTS

An auto-encoder with 1 real latent variable represents single-tissue fast spin echo, EPTI, and MPRAGE signal evolution to within 0.15% normalized RMS error, enabling reconstruction problems with 3 degrees of freedom per voxel (real latent variable + complex scaling) in comparison to linear models with 4-8 degrees of freedom per voxel. In simulated/in vivo T₂ -shuffling and in vivo EPTI experiments, the proposed framework achieves consistent quantitative normalized RMS error improvement over linear approaches. From qualitative evaluation, the proposed approach yields images with reduced blurring and noise amplification in MPRAGE-shuffling experiments.

CONCLUSION

Directly solving for nonlinear latent representations of signal evolution improves time-resolved MRI reconstructions.

Optimized three-dimensional ultrashort echo time: Magnetic resonance fingerprinting for myelin tissue fraction mapping

Quantitative assessment of brain myelination has gained attention for both research and diagnosis of neurological diseases. However, conventional pulse sequences cannot directly acquire the myelin-proton signals due to its extremely short T2 and T2* values. To obtain the myelin-proton signals, dedicated short T2 acquisition techniques, such as ultrashort echo time (UTE) imaging, have been introduced. However, it remains challenging to isolate the myelin-proton signals from tissues with longer T2. In this article, we extended our previous two-dimensional ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) with dual-echo acquisition to three dimensional (3D). Given a relatively low proton density (PD) of myelin-proton, we utilized Cramér-Rao Lower Bound to encode myelin-proton with the maximal SNR efficiency for optimizing the MR fingerprinting design, in order to improve the sensitivity of the sequence to myelin-proton. In addition, with a second echo of approximately 3 ms, myelin-water component can be also captured. A myelin-tissue (myelin-proton and myelin-water) fraction mapping can be thus calculated. The optimized 3D UTE-MRF with dual-echo acquisition is tested in simulations, physical phantom and in vivo studies of both healthy subjects and multiple sclerosis patients. The results suggest that the rapidly decayed myelin-proton and myelin-water signal can be depicted with UTE signals of our method at clinically relevant resolution (1.8 mm isotropic) in 15 min. With its good sensitivity to myelin loss in multiple sclerosis patients demonstrated, our method for the whole brain myelin-tissue fraction mapping in clinical friendly scan time has the potential for routine clinical imaging.

DEEP FACTOR MODEL: A NOVEL APPROACH FOR MOTION COMPENSATED MULTI-DIMENSIONAL MRI

Recent quantitative parameter mapping methods including MR fingerprinting (MRF) collect a time series of images that capture the evolution of magnetization. The focus of this work is to introduce a novel approach termed as Deep Factor Model(DFM), which offers an efficient representation of the multi-contrast image time series. The higher efficiency of the representation enables the acquisition of the images in a highly undersampled fashion, which translates to reduced scan time in 3D high-resolution multi-contrast applications. The approach integrates motion estimation and compensation, making the approach robust to subject motion during the scan.

Multicomponent MR fingerprinting reconstruction using joint-sparsity and low-rank constraints

PURPOSE

To develop an efficient algorithm for multicomponent MR fingerprinting (MC-MRF) reconstructions directly from highly undersampled data without making prior assumptions about tissue relaxation times and expected number of tissues.

METHODS

The proposed method reconstructs MC-MRF maps from highly undersampled data by iteratively applying a joint-sparsity constraint to the estimated tissue components. Intermediate component maps are obtained by a low-rank multicomponent alternating direction method of multipliers (MC-ADMM) including the non-negativity of tissue weights as an extra regularization term. Over iterations, the used dictionary compression is adjusted. The proposed method (k-SPIJN) is compared with a two-step approach in which image reconstruction and multicomponent estimations are performed sequentially and tested in numerical simulations and in vivo by applying different undersampling factors in eight healthy volunteers. In the latter case, fully sampled data serves as the reference.

RESULTS

The proposed method shows improved precision and accuracy in simulations compared with a state-of-art sequential approach. Obtained in vivo magnetization fraction maps for different tissue types show reduced systematic errors and reduced noise-like effects. Root mean square errors in estimated magnetization fraction maps significantly reduce from 13.0% <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>±</mml:mo></mml:mrow> <mml:annotation>$$ \pm $$</mml:annotation></mml:semantics> </mml:math> 5.8% with the conventional, two-step approach to 9.6% <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>±</mml:mo></mml:mrow> <mml:annotation>$$ \pm $$</mml:annotation></mml:semantics> </mml:math> 3.9% and 9.6% <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>±</mml:mo></mml:mrow> <mml:annotation>$$ \pm $$</mml:annotation></mml:semantics> </mml:math> 3.2% with the proposed MC-ADMM and k-SPIJN methods, respectively. Mean standard deviation in homogeneous white matter regions reduced significantly from 8.6% to 2.9% (two step vs. k-SPIJN).

CONCLUSION

The proposed MC-ADMM and k-SPIJN reconstruction methods estimate MC-MRF maps from highly undersampled data resulting in improved image quality compared with the existing method.

3D magnetic resonance fingerprinting for rapid simultaneous T1, T2, and T1ρ volumetric mapping of human articular cartilage at 3 T

Quantitative MRI can detect early biochemical changes in cartilage; however, the conventional techniques only measure one parameter (e.g., T1 , T2 , and T1ρ ) at a time while also being comparatively slow. We implemented a 3D magnetic resonance fingerprinting (3D-MRF) technique for simultaneous, volumetric mapping of T1 , T2 , and T1ρ in knee articular cartilage in under 9 min. It is evaluated on 11 healthy volunteers (mean age: 53 ± 9 years), five mild knee osteoarthritis (OA) patients (Kellgren-Lawrence (KL) score: 2, mean age: 60 ± 4 years), and the National Institute of Standards and Technology (NIST)/International Society for Magnetic Resonance in Medicine (ISMRM) system phantom. Proton density image, and T1 , T2, T1ρ relaxation times, and B1 + were estimated in the NIST/ISMRM system phantom as well as in the human knee medial and lateral femur, medial and lateral tibia, and patellar cartilage. The repeatability and reproducibility of the proposed technique were assessed in the phantom using analysis of the Bland-Altman plots. The intrasubject repeatability was assessed with the coefficient of variation (CV) and root mean square CV (rmsCV). The Mann-Whitney U test was used to assess the difference between healthy subjects and mild knee OA patients. The Bland-Altman plots in the NIST/ISMRM phantom demonstrated an average difference of 0.001% ± 015%, 1.2% ± 7.1%, and 0.47% ± 3% between two scans from the same 3-T scanner (repeatability), and 0.002% ± 015%, 0.62% ± 10.5%, and 0.97% ± 14% between the scans acquired on two different 3-T scanners (reproducibility) for T1 , T2 , and T1ρ , respectively. The in vivo knee study showed excellent repeatability with rmsCV less than 1%, 2%, and 1% for T1 , T2 , and T1ρ , respectively. T1ρ relaxation time in the mild knee OA patients was significantly higher (p < 0.05) than in healthy subjects. The proposed 3D-MRF sequence is fast, reproducible, robust to B1 + inhomogeneity, and can simultaneously measure the T1 , T2 , T1ρ , and B1 + volumetric maps of the knee joint in a single scan within a clinically feasible scan time.

Spatiotemporal encoding MRI using subspace-constrained sampling and locally-low-rank regularization: Applications to diffusion weighted and diffusion kurtosis imaging of human brain and prostate

The benefits of performing locally low-rank (LLR) reconstructions on subsampled diffusion weighted and diffusion kurtosis imaging data employing spatiotemporal encoding (SPEN) methods, is investigated. SPEN allows for self-referenced correction of motion-induced phase errors in case of interleaved diffusion-oriented acquisitions, and allows one to overcome distortions otherwise observed along EPI's phase-encoded dimension. In combination with LLR-based reconstructions of the pooled imaging data and with a joint subsampling of b-weighted and interleaved images, additional improvements in terms of sensitivity as well as shortened acquisition times are demonstrated, without noticeable penalties. Details on how the LLR-regularized, subspace-constrained image reconstructions were adapted to SPEN are given; the improvements introduced by adopting these reconstruction frameworks for the accelerated acquisition of diffusivity and of kurtosis imaging data in both relatively homogeneous regions like the human brain and in more challenging regions like the human prostate, are presented and discussed within the context of similar efforts in the field.

Artificial intelligence in cardiac magnetic resonance fingerprinting

Magnetic resonance fingerprinting (MRF) is a fast MRI-based technique that allows for multiparametric quantitative characterization of the tissues of interest in a single acquisition. In particular, it has gained attention in the field of cardiac imaging due to its ability to provide simultaneous and co-registered myocardial T₁ and T₂ mapping in a single breath-held cardiac MRF scan, in addition to other parameters. Initial results in small healthy subject groups and clinical studies have demonstrated the feasibility and potential of MRF imaging. Ongoing research is being conducted to improve the accuracy, efficiency, and robustness of cardiac MRF. However, these improvements usually increase the complexity of image reconstruction and dictionary generation and introduce the need for sequence optimization. Each of these steps increase the computational demand and processing time of MRF. The latest advances in artificial intelligence (AI), including progress in deep learning and the development of neural networks for MRI, now present an opportunity to efficiently address these issues. Artificial intelligence can be used to optimize candidate sequences and reduce the memory demand and computational time required for reconstruction and post-processing. Recently, proposed machine learning-based approaches have been shown to reduce dictionary generation and reconstruction times by several orders of magnitude. Such applications of AI should help to remove these bottlenecks and speed up cardiac MRF, improving its practical utility and allowing for its potential inclusion in clinical routine. This review aims to summarize the latest developments in artificial intelligence applied to cardiac MRF. Particularly, we focus on the application of machine learning at different steps of the MRF process, such as sequence optimization, dictionary generation and image reconstruction.

Motion-corrected 3D-EPTI with efficient 4D navigator acquisition for fast and robust whole-brain quantitative imaging

PURPOSE

To develop a motion estimation and correction method for motion-robust three-dimensional (3D) quantitative imaging with 3D-echo-planar time-resolved imaging.

THEORY AND METHODS

The 3D-echo-planar time-resolved imaging technique was designed with additional four-dimensional navigator acquisition (x-y-z-echoes) to achieve fast and motion-robust quantitative imaging of the human brain. The four-dimensional-navigator is inserted into the relaxation-recovery deadtime of the sequence in every pulse TR (∼2 s) to avoid extra scan time, and to provide continuous tracking of the 3D head motion and B0 -inhomogeneity changes. By using an optimized spatiotemporal encoding combined with a partial-Fourier scheme, the navigator acquires a large central k-t data block for accurate motion estimation using only four small-flip-angle excitations and readouts, resulting in negligible signal-recovery reduction to the 3D-echo-planar time-resolved imaging acquisition. By incorporating the estimated motion and B0 -inhomogeneity changes into the reconstruction, multi-contrast images can be recovered with reduced motion artifacts.

RESULTS

Simulation shows the cost to the SNR efficiency from the added navigator acquisitions is <1%. Both simulation and in vivo retrospective experiments were conducted, that demonstrate the four-dimensional navigator provided accurate estimation of the 3D motion and B0 -inhomogeneity changes, allowing effective reduction of image artifacts in quantitative maps. Finally, in vivo prospective undersampling acquisition was performed with and without head motion, in which the motion corrupted data after correction show close image quality and consistent quantifications to the motion-free scan, providing reliable quantitative measurements even with head motion.

CONCLUSION

The proposed four-dimensional navigator acquisition provides reliable tracking of the head motion and B0 change with negligible SNR cost, equips the 3D-echo-planar time-resolved imaging technique for motion-robust and efficient quantitative imaging.

MR fingerprinting of the prostate

Multiparametric magnetic resonance imaging (mpMRI) has been adopted as the key tool for detection, localization, characterization, and risk stratification of patients suspected to have prostate cancer. Despite advantages over systematic biopsy, the interpretation of prostate mpMRI has limitations including a steep learning curve, leading to considerable interobserver variation. There is growing interest in clinical translation of quantitative imaging techniques for more objective lesion assessment. However, traditional mapping techniques are slow, precluding their use in the clinic. Magnetic resonance fingerprinting (MRF) is an efficient approach for quantitative maps of multiple tissue properties simultaneously. The T₁ and T₂ values obtained with MRF have been validated with phantom studies as well as in normal volunteers and patients. Studies have shown that MRF-derived T₁ and T₂ along with ADC values are all significant independent predictors in the differentiation between normal prostate tissue and prostate cancer, and hold promise in differentiating low and intermediate/high-grade cancers. This review seeks to introduce the basics of the prostate MRF technique, discuss the potential applications of prostate MRF for the characterization of prostate cancer, and describes ongoing areas of research.

Improved Balanced Steady-State Free Precession Based MR Fingerprinting with Deep Autoencoders

Magnetic Resonance (MR) Fingerprinting is an emerging transient-state imaging paradigm, which enables the quantization of multiple MR tissue parameters in a single experiment. Balanced steady-state free precession (bSSFP)-based MR Fingerprinting has excellent signal-to-noise characteristics and also allows for acquiring both tissue parameter maps and field inhomogeneity maps. However, field inhomogeneity often results in complex magnetization evolutions in bSSFP-based MR Fingerprinting, which creates significant challenges in image reconstruction. In this paper, we introduce a new method to address the image reconstruction problem. The proposed method incorporates a low-dimensional nonlinear manifold learned from an ensemble of magnetization evolutions using a deep autoencoder. It provides much better representation power than a low-dimensional linear subspace in capturing complex magnetization evolutions. We formulate the image reconstruction problem with this nonlinear model and solve the resulting optimization problem using an algorithm based on variable splitting and the alternating direction method of multipliers. We evaluate the performance of the proposed method using numerical experiments and demonstrate that it significantly outperforms the state-of-art method using a linear subspace model.

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.

Simultaneous bilateral T1 , T2 , and T1ρ relaxation mapping of the hip joint with magnetic resonance fingerprinting

Quantitative MRI can detect early biochemical changes in cartilage, but its bilateral use in clinical routines is challenging. The aim of this prospective study was to demonstrate the feasibility of magnetic resonance fingerprinting for bilateral simultaneous T₁ , T₂ , and T_1ρ mapping of the hip joint. The study population consisted of six healthy volunteers with no known trauma or pain in the hip. Monoexponential T₁ , T₂ , and T_1ρ relaxation components were assessed in femoral lateral, superolateral, and superomedial, and inferior, as well as acetabular, superolateral, and superomedial subregions in left and right hip cartilage. Aligned ranked nonparametric factorial analysis was used to assess the side's impact on the subregions. Kruskal-Wallis and Wilcoxon tests were used to compare subregions, and coefficient of variation to assess repeatability. Global averages of T₁ (676.0 ± 45.4 and 687.6 ± 44.5 ms), T₂ (22.5 ± 2.6 and 22.1 ± 2.5 ms), and T_1ρ (38.2 ± 5.5 and 38.2 ± 5.5 ms) were measured in the left and right hip, and articular cartilage, respectively. The Kruskal-Wallis test showed a significant difference between different subregions' relaxation times regardless of the hip side (p < 0.001 for T₁ , p = 0.012 for T₂ , and p < 0.001 for T_1ρ ). The Wilcoxon test showed that T₁ of femoral layers was significantly (p < 0.003) higher than that for acetabular cartilage. The experiments showed excellent repeatability with CV_rms of 1%, 2%, and 4% for T₁ , T₂ , and T_1ρ, respectively. It was concluded that bilateral T₁ , T₂ , and T_1ρ relaxation times, as well as B₁⁺ maps, can be acquired simultaneously from hip joints using the proposed MRF sequence.

Rapid MR relaxometry using deep learning: An overview of current techniques and emerging trends

Quantitative mapping of MR tissue parameters such as the spin-lattice relaxation time (T₁ ), the spin-spin relaxation time (T₂ ), and the spin-lattice relaxation in the rotating frame (T_1ρ ), referred to as MR relaxometry in general, has demonstrated improved assessment in a wide range of clinical applications. Compared with conventional contrast-weighted (eg T₁ -, T₂ -, or T_1ρ -weighted) MRI, MR relaxometry provides increased sensitivity to pathologies and delivers important information that can be more specific to tissue composition and microenvironment. The rise of deep learning in the past several years has been revolutionizing many aspects of MRI research, including image reconstruction, image analysis, and disease diagnosis and prognosis. Although deep learning has also shown great potential for MR relaxometry and quantitative MRI in general, this research direction has been much less explored to date. The goal of this paper is to discuss the applications of deep learning for rapid MR relaxometry and to review emerging deep-learning-based techniques that can be applied to improve MR relaxometry in terms of imaging speed, image quality, and quantification robustness. The paper is comprised of an introduction and four more sections. Section 2 describes a summary of the imaging models of quantitative MR relaxometry. In Section 3, we review existing "classical" methods for accelerating MR relaxometry, including state-of-the-art spatiotemporal acceleration techniques, model-based reconstruction methods, and efficient parameter generation approaches. Section 4 then presents how deep learning can be used to improve MR relaxometry and how it is linked to conventional techniques. The final section concludes the review by discussing the promise and existing challenges of deep learning for rapid MR relaxometry and potential solutions to address these challenges.

3D MR fingerprinting using Seiffert spirals

PURPOSE

Seiffert spirals were recently explored as an efficient way to traverse 3D k-space compared to traditional 3D techniques. Several studies have shown the ability of 3D MR fingerprinting (MRF) techniques to acquire T₁ and T₂ relaxation maps in a short period of time. However, these sequences do not sample across a large region of 3D k-space every TR, especially in the way that Seiffert trajectories can.

METHODS

A 3D MRF sequence was designed using 8 Seiffert spirals rotated in 3D k-space, with flip angle modulation for T₁ and T₂ sensitivity. The sequence was compared to an MRF sequence using a 2D spiral rotated in 3D k-space using the tiny golden angle acquisition with similar resolution/readout duration. Both sequences were evaluated using simulations, phantom validation, and in vivo imaging.

RESULTS

In all experiments, the Seiffert spiral MRF sequence performed similar to if not better than the multi-axis 2D spiral MRF sequence. Strong intraclass correlation coefficients (> 0.9) were found between conventional and MRF sequences in phantoms, whereas the in vivo results showed slightly less aliasing artifact with the Seiffert trajectory.

CONCLUSION

In this study, Seiffert spirals were used within the MRF framework to acquire high-resolution T₁ and T₂ relaxation time maps in less than 2.5 min. The reduced aliasing artifacts seen with the Seiffert sequence suggests that sampling over 3D k-space evenly each TR can improve quantification or shorten scan times.

Data-Consistent non-Cartesian deep subspace learning for efficient dynamic MR image reconstruction

Non-Cartesian sampling with subspace-constrained image reconstruction is a popular approach to dynamic MRI, but slow iterative reconstruction limits its clinical application. Data-consistent (DC) deep learning can accelerate reconstruction with good image quality, but has not been formulated for non-Cartesian subspace imaging. In this study, we propose a DC non-Cartesian deep subspace learning framework for fast, accurate dynamic MR image reconstruction. Four novel DC formulations are developed and evaluated: two gradient decent approaches, a directly solved approach, and a conjugate gradient approach. We applied a U-Net model with and without DC layers to reconstruct T1-weighted images for cardiac MR Multitasking (an advanced multidimensional imaging method), comparing our results to the iteratively reconstructed reference. Experimental results show that the proposed framework significantly improves reconstruction accuracy over the U-Net model without DC, while significantly accelerating the reconstruction over conventional iterative reconstruction.

Optimized multi-axis spiral projection MR fingerprinting with subspace reconstruction for rapid whole-brain high-isotropic-resolution quantitative imaging

PURPOSE

To improve image quality and accelerate the acquisition of 3D MR fingerprinting (MRF).

METHODS

Building on the multi-axis spiral-projection MRF technique, a subspace reconstruction with locally low-rank constraint and a modified spiral-projection spatiotemporal encoding scheme called tiny golden-angle shuffling were implemented for rapid whole-brain high-resolution quantitative mapping. Reconstruction parameters such as the locally low-rank regularization parameter and the subspace rank were tuned using retrospective in vivo data and simulated examinations. B₀ inhomogeneity correction using multifrequency interpolation was incorporated into the subspace reconstruction to further improve the image quality by mitigating blurring caused by off-resonance effect.

RESULTS

The proposed MRF acquisition and reconstruction framework yields high-quality 1-mm isotropic whole-brain quantitative maps in 2 min at better quality compared with 6-min acquisitions of prior approaches. The proposed method was validated to not induce bias in T₁ and T₂ mapping. High-quality whole-brain MRF data were also obtained at 0.66-mm isotropic resolution in 4 min using the proposed technique, where the increased resolution was shown to improve visualization of subtle brain structures.

CONCLUSIONS

The proposed tiny golden-angle shuffling, MRF with optimized spiral-projection trajectory and subspace reconstruction enables high-resolution quantitative mapping in ultrafast acquisition time.

MR fingerprinting for semisolid magnetization transfer and chemical exchange saturation transfer quantification

Chemical exchange saturation transfer (CEST) MRI has positioned itself as a promising contrast mechanism, capable of providing molecular information at sufficient resolution and amplified sensitivity. However, it has not yet become a routinely employed clinical technique, due to a variety of confounding factors affecting its contrast-weighted image interpretation and the inherently long scan time. CEST MR fingerprinting (MRF) is a novel approach for addressing these challenges, allowing simultaneous quantitation of several proton exchange parameters using rapid acquisition schemes. Recently, a number of deep-learning algorithms have been developed to further boost the performance and speed of CEST and semi-solid macromolecule magnetization transfer (MT) MRF. This review article describes the fundamental theory behind semisolid MT/CEST-MRF and its main applications. It then details supervised and unsupervised learning approaches for MRF image reconstruction and describes artificial intelligence (AI)-based pipelines for protocol optimization. Finally, practical considerations are discussed, and future perspectives are given, accompanied by basic demonstration code and data.

Three-dimensional high-resolution T1 and T2 mapping of whole macaque brain at 9.4 T using magnetic resonance fingerprinting

PURPOSE

Quantitative T₁ and T₂ mapping in non-human primates with whole-brain coverage is challenged by the requirement of sub-millimeter resolution and the inhomogeneity of the transmit magnetic field (B₁ ⁺ ) covering a large field of view. The goal of the current study is to develop a magnetic resonance fingerprinting (MRF) method for simultaneous T₁ and T₂ mapping of the entire macaque brain within feasible scan time.

METHODS

A three-dimensional (3D) MRF sequence with both inversion- and T₂ -preparation modules was developed and evaluated on a 9.4 T preclinical scanner. Data acquisition used a 3D stack-of-spirals trajectory, with undersampling along both the in-plane and the through-plane directions. The effect of B₁ ⁺ inhomogeneity was accounted for by matching the acquired fingerprint to a dictionary simulated with the B₁ ⁺ factors measured from a separate scan. In vitro and ex vivo studies were performed to evaluate the accuracy and the undersampling capacity of the MRF method. The application of the MRF method for in vivo, brain-wide T₁ and T₂ mapping was demonstrated on macaques at 4, 6, and 12 years of age.

RESULTS

The MRF method enabled highly repeatable T₁ and T₂ mapping at high spatial resolution (0.35 × 0.35 × 1 mm³ ) with an acceleration factor of 24. In vivo studies showed significant age-related T₂ reduction in deep gray nuclei including the globus pallidus, the putamen, and the caudate nucleus.

CONCLUSIONS

This study demonstrates the first MRF study for brain-wide, multi-parametric quantification in non-human primates with sub-millimeter resolution.

3D Echo Planar Time-resolved Imaging (3D-EPTI) for ultrafast multi-parametric quantitative MRI

Multi-parametric quantitative MRI has shown great potential to improve the sensitivity and specificity of clinical diagnosis and to enhance our understanding of complex brain processes, but suffers from long scan time especially at high spatial resolution. To address this longstanding challenge, we introduce a novel approach termed 3D Echo Planar Time-resolved Imaging (3D-EPTI), which significantly increases the acceleration capacity of MRI sampling, and provides high acquisition efficiency for multi-parametric MRI. This is achieved by exploiting the spatiotemporal correlation of MRI data at multiple timescales through new encoding strategies within and between efficient continuous readouts. Specifically, an optimized spatiotemporal CAIPI encoding within the readouts combined with a radial-block sampling strategy across the readouts enables an acceleration rate of 800 fold in the k-t space. A subspace reconstruction was employed to resolve thousands of high-quality multi-contrast images. We have demonstrated the ability of 3D-EPTI to provide robust and repeatable whole-brain simultaneous T₁, T₂, T₂*, PD and B₁⁺ mapping at high isotropic resolution within minutes (e.g., 1-mm isotropic resolution in 3 minutes), and to enable submillimeter multi-parametric imaging to study detailed brain structures.

Tailored Magnetic Resonance Fingerprinting for simultaneous non-synthetic and quantitative imaging: A repeatability study

PURPOSE

The goals of this study include: (a) generating tailored magnetic resonance fingerprinting (TMRF) based non-synthetic imaging; (b) assessing the repeatability of TMRF and deep learning-based mapping of in vitro ISMRM/NIST phantom and in vivo brain data of healthy human subjects.

METHODS

We have acquired qualitative images obtained from the vendor-supplied gold standard, MRF (synthetic), and TMRF (non-synthetic) on one representative healthy human brain. We also acquired thirty datasets on the ISMRM/NIST phantom for the in vitro repeatability study on a GE Discovery 3T MR750w scanner using the TMRF sequence. We compared T₁ and T₂ maps generated from thirty ISMRM/NIST phantom datasets to the spin-echo (SE) based gold standard (GS) method as part of the in vitro repeatability study. R-squared coefficient of determination in a simple linear regression and Bland-Altman analysis were computed for thirty datasets of ISMRM/NIST phantom to assess the accuracy of in vitro quantitative TMRF data. The repeatability of T₁ and T₂ estimates by TMRF was evaluated by calculating the standard deviation (SD) divided by the average of thirty datasets for each sphere, respectively. We acquired ten volunteers for the in vivo repeatability study on the same scanner using the same TMRF sequence. These volunteers were imaged five times with two runs per repetition, resulting in one hundred in vivo datasets. Five contrasts, T₁ and T₂ maps of ten human volunteers acquired over five repetitions, were evaluated in the in vivo repeatability study. We computed the intraclass correlation coefficient (ICC) of the signal-to-noise ratio (SNR), signal intensities, T₁ and T₂ relaxation times in white matter (WM) and gray matter (GM).

RESULTS

The synthetic images generated from MRF show partial volume and flow artifacts compared to non-synthetic images obtained from TMRF images and the gold standard. in vitro studies show that TMRF estimates have less than 5% variations except sphere 14 in the T₂ array (6.36%). TMRF and spin-echo relaxometry measurements were strongly correlated; R² values were 0.9958 and 0.9789 for T₁ and T₂ estimates, respectively. Based on the ICC values, SNR, mean intensity values, and relaxation times of WM and GM for the in vivo studies were consistent. T₁ and T₂ values of WM and GM were similar to previously published values. The mean ± SD of T₁ and T₂ for WM for ten subjects and five repeats are 992 ± 41 ms and 99 ± 6 ms, while the corresponding values for T₁ and T₂ for GM are 1598 ± 73 ms and 152 ± 14 ms.

CONCLUSION

TMRF and deep learning-based reconstruction produce repeatable, non-synthetic multi-contrast images and parametric maps simultaneously. This article is protected by copyright. All rights reserved.

Myocardial T1, T2, T2*, and fat fraction quantification via low-rank motion-corrected cardiac MR fingerprinting

Purpose

Develop a novel 2D cardiac MR fingerprinting (MRF) approach to enable simultaneous T1, T2, T2*, and fat fraction (FF) myocardial tissue characterization in a single breath‐hold scan.

Methods

Simultaneous, co‐registered, multi‐parametric mapping of T1, T2, and FF has been recently achieved with cardiac MRF. Here, we further incorporate T2* quantification within this approach, enabling simultaneous T1, T2, T2*, and FF myocardial tissue characterization in a single breath‐hold scan. T2* quantification is achieved with an eight‐echo readout that requires a long cardiac acquisition window. A novel low‐rank motion‐corrected (LRMC) reconstruction is exploited to correct for cardiac motion within the long acquisition window. The proposed T1/T2/T2*/FF cardiac MRF was evaluated in phantom and in 10 healthy subjects in comparison to conventional mapping techniques.

Results

The proposed approach achieved high quality parametric mapping of T1, T2, T2*, and FF with corresponding normalized RMS error (RMSE) T1 = 5.9%, T2 = 9.6% (T2 values <100 ms), T2* = 3.3% (T2* values <100 ms), and FF = 0.8% observed in phantom scans. In vivo, the proposed approach produced higher left‐ventricular myocardial T1 values than MOLLI (1148 vs 1056 ms), lower T2 values than T2‐GraSE (42.8 vs 50.6 ms), lower T2* values than eight‐echo gradient echo (GRE) (35.0 vs 39.4 ms), and higher FF values than six‐echo GRE (0.8 vs 0.3 %) reference techniques. The proposed approach achieved considerable reduction in motion artifacts compared to cardiac MRF without motion correction, improved spatial uniformity, and statistically higher apparent precision relative to conventional mapping for all parameters.

Conclusion

The proposed cardiac MRF approach enables simultaneous, co‐registered mapping of T1, T2, T2*, and FF in a single breath‐hold for comprehensive myocardial tissue characterization, achieving higher apparent precision than conventional methods.

An MR fingerprinting approach for quantitative inhomogeneous magnetization transfer imaging

PURPOSE

Magnetization transfer (MT) and inhomogeneous MT (ihMT) contrasts are used in MRI to provide information about macromolecular tissue content. In particular, MT is sensitive to macromolecules, and ihMT appears to be specific to myelinated tissue. This study proposes a technique to characterize MT and ihMT properties from a single acquisition, producing both semiquantitative contrast ratios and quantitative parameter maps.

THEORY AND METHODS

Building on previous work that uses multiband RF pulses to efficiently generate ihMT contrast, we propose a cyclic steady-state approach that cycles between multiband and single-band pulses to boost the achieved contrast. Resultant time-variable signals are reminiscent of an MR fingerprinting acquisition, except that the signal fluctuations are entirely mediated by MT effects. A dictionary-based low-rank inversion method is used to reconstruct the resulting images and to produce both semiquantitative MT ratio and ihMT ratio maps, as well as quantitative parameter estimates corresponding to an ihMT tissue model.

RESULTS

Phantom and in vivo brain data acquired at 1.5 Tesla demonstrate the expected contrast trends, with ihMT ratio maps showing contrast more specific to white matter, as has been reported by others. Quantitative estimation of semisolid fraction and dipolar T₁ was also possible and yielded measurements consistent with literature values in the brain.

CONCLUSION

By cycling between multiband and single-band pulses, an entirely MT-mediated fingerprinting method was demonstrated. This proof-of-concept approach can be used to generate semiquantitative maps and quantitatively estimate some macromolecular-specific tissue parameters.

A Spatial Off-Resonance Correction in Spirals for Magnetic Resonance Fingerprinting

In MR Fingerprinting (MRF), balanced Steady-State Free Precession (bSSFP) has advantages over unbalanced SSFP because it retains the spin history achieving a higher signal-to-noise ratio (SNR) and scan efficiency. However, bSSFP-MRF is not frequently used because it is sensitive to off-resonance, producing artifacts and blurring, and affecting the parametric map quality. Here we propose a novel Spatial Off-resonance Correction (SOC) approach for reducing these artifacts in bSSFP-MRF with spiral trajectories. SOC-MRF uses each pixel's Point Spread Function to create system matrices that encode both off-resonance and gridding effects. We iteratively compute the inverse of these matrices to reduce the artifacts. We evaluated the proposed method using brain simulations and actual MRF acquisitions of a standardized T1/T2 phantom and five healthy subjects. The results show that the off-resonance distortions in T1/T2 maps were considerably reduced using SOC-MRF. For T2, the Normalized Root Mean Square Error (NRMSE) was reduced from 17.3 to 8.3% (simulations) and from 35.1 to 14.9% (phantom). For T1, the NRMS was reduced from 14.7 to 7.7% (simulations) and from 17.7 to 6.7% (phantom). For in-vivo, the mean and standard deviation in different ROI in white and gray matter were significantly improved. For example, SOC-MRF estimated an average T2 for white matter of 77ms (the ground truth was 74ms) versus 50 ms of MRF. For the same example the standard deviation was reduced from 18 ms to 6ms. The corrections achieved with the proposed SOC-MRF may expand the potential applications of bSSFP-MRF, taking advantage of its better SNR property.

Automated design of pulse sequences for magnetic resonance fingerprinting using physics-inspired optimization

Magnetic resonance fingerprinting (MRF) is a method to extract quantitative tissue properties such as [Formula: see text] and [Formula: see text] relaxation rates from arbitrary pulse sequences using conventional MRI hardware. MRF pulse sequences have thousands of tunable parameters, which can be chosen to maximize precision and minimize scan time. Here, we perform de novo automated design of MRF pulse sequences by applying physics-inspired optimization heuristics. Our experimental data suggest that systematic errors dominate over random errors in MRF scans under clinically relevant conditions of high undersampling. Thus, in contrast to prior optimization efforts, which focused on statistical error models, we use a cost function based on explicit first-principles simulation of systematic errors arising from Fourier undersampling and phase variation. The resulting pulse sequences display features qualitatively different from previously used MRF pulse sequences and achieve fourfold shorter scan time than prior human-designed sequences of equivalent precision in [Formula: see text] and [Formula: see text] Furthermore, the optimization algorithm has discovered the existence of MRF pulse sequences with intrinsic robustness against shading artifacts due to phase variation.

Generalized low-rank nonrigid motion-corrected reconstruction for MR fingerprinting

PURPOSE

Develop a novel low-rank motion-corrected (LRMC) reconstruction for nonrigid motion-corrected MR fingerprinting (MRF).

METHODS

Generalized motion-corrected (MC) reconstructions have been developed for steady-state imaging. Here we extend this framework to enable nonrigid MC for transient imaging applications with varying contrast, such as MRF. This is achieved by integrating low-rank dictionary-based compression into the generalized MC model to reconstruct MC singular images, reducing motion artifacts in the resulting parametric maps. The proposed LRMC reconstruction was applied for cardiac motion correction in 2D myocardial MRF (T₁ and T₂ ) with extended cardiac acquisition window (~450 ms) and for respiratory MC in free-breathing 3D myocardial and 3D liver MRF. Experiments were performed in phantom and 22 healthy subjects. The proposed approach was compared with reference spin echo (phantom) and with 2D electrocardiogram-triggered/breath-hold MOLLI and T₂ gradient-and-spin echo conventional maps (in vivo 2D and 3D myocardial MRF).

RESULTS

Phantom results were in general agreement with reference spin-echo measurements, presenting relative errors of approximately 5.4% and 5.5% for T₁ and short T₂ (<100 ms), respectively. The proposed LRMC MRF reduced residual blurring artifacts with respect to no MC for cardiac or respiratory motion in all cases (2D and 3D myocardial, 3D abdominal). In 2D myocardial MRF, left-ventricle T₁ values were 1150 ± 41 ms for LRMC MRF and 1010 ± 56 ms for MOLLI; T₂ values were 43.8 ± 2.3 ms for LRMC MRF and 49.5 ± 4.5 ms for T₂ gradient and spin echo. Corresponding measurements for 3D myocardial MRF were 1085 ± 30 ms and 1062 ± 29 ms for T₁ , and 43.5 ± 1.9 ms and 51.7 ± 1.7 ms for T₂ . For 3D liver, LRMC MRF measured liver T₁ at 565 ± 44 ms and liver T₂ at 35.4 ± 2.4 ms.

CONCLUSION

The proposed LRMC reconstruction enabled generalized (nonrigid) MC for 2D and 3D MRF, both for cardiac and respiratory motion. The proposed approach reduced motion artifacts in the MRF maps with respect to no motion compensation and achieved good agreement with reference measurements.

Current Applications and Future Development of Magnetic Resonance Fingerprinting in Diagnosis, Characterization, and Response Monitoring in Cancer

Magnetic resonance imaging (MRI) has enabled non-invasive cancer diagnosis, monitoring, and management in common clinical settings. However, inadequate quantitative analyses in MRI continue to limit its full potential and these often have an impact on clinicians' judgments. Magnetic resonance fingerprinting (MRF) has recently been introduced to acquire multiple quantitative parameters simultaneously in a reasonable timeframe. Initial retrospective studies have demonstrated the feasibility of using MRF for various cancer characterizations. Further trials with larger cohorts are still needed to explore the repeatability and reproducibility of the data acquired by MRF. At the moment, technical difficulties such as undesirable processing time or lack of motion robustness are limiting further implementations of MRF in clinical oncology. This review summarises the latest findings and technology developments for the use of MRF in cancer management and suggests possible future implications of MRF in characterizing tumour heterogeneity and response assessment.