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

3D deuterium metabolic imaging (DMI) of the human liver at 7 T using low-rank and subspace model-based reconstruction

PURPOSE

To implement a low-rank and subspace model-based reconstruction for 3D deuterium metabolic imaging (DMI) and compare its performance against Fourier transform-based (FFT) reconstruction in terms of spectral fitting reliability.

METHODS

Both reconstruction methods were applied on simulated and experimental DMI data. Numerical simulations were performed to evaluate the effect of increasing acceleration factors. The impact on spectral fitting results, SNR, and the overall normalized root mean square error (NRMSE) compared to ground-truth data were calculated. A comparative analysis was performed on DMI data acquired from the human liver, including both natural abundance and post-deuterated glucose intake data at 7 T.

RESULTS

Simulation showed the Cramer-Rao lower bound [%] of water, glucose, sum of glutamate and glutamine (Glx), and lipid signals for the low-rank and subspace model-based reconstruction at R = 1.0 was 12.4, 14.7, 17.3, and 11.0 times lower than FFT. At R = 1.1, NRMSE was 1.4%, 1.3%, 0.8%, and 4.2% lower for the water, glucose, Glx, and lipid, respectively, compared to FFT. However, the NRMSE of the Glx and lipid increased by 0.4% and 3.2% at R = 1.3. For the in vivo DMI experiment, SNR was 2.5-3.0 times higher compared to FFT. The fitted amplitude of water and glucose peaks showed Cramer-Rao lower bound [%] values that were approximately 2.3 times lower than FFT.

CONCLUSION

Simulations and in vivo experiments on the human liver demonstrate that low-rank and subspace model-based reconstruction with undersampled data mitigates noise and enhances spectral fitting quality.

Efficient 3D FISP-MRF at 0.55 T using long spiral readouts and concomitant field effect mitigation

PURPOSE

To demonstrate the feasibility of SNR-efficient 3D fast imaging with steady state precession MR fingerprinting (FISP-MRF) using long spiral readouts with mitigation of concomitant field effects at 0.55 T.

METHODS

Fourteen FISP-MRF sequences with different spiral readout lengths (2.9 ms to 22.0 ms) were implemented with the open-source Pulseq framework. Datasets were reconstructed using a low-rank and subspace model-based reconstruction combined with MaxGIRF spatial encoding model. Concomitant field-induced blurring and MRF precision were evaluated using reconstructed images and quantitative maps acquired from an ACR phantom, an ISMRM/NIST system phantom, and 2 healthy volunteers.

RESULTS

A simulation study shows that the SNR of time-series images would increase by ∼2× in white matter as a spiral readout increased from 2.9 ms to 22.0 ms. Empirically, MRF T1 and T2 standard deviations of in-vivo white matter were reduced by ∼50 %, and concomitant field mitigation was necessary. Residual spatial blurring was non-negligible for readouts ≥16.5 ms, suggesting an operating regime (2.9-16.5 ms) for 3D FISP-MRF at 0.55 T.

CONCLUSION

We demonstrate SNR-efficient 3D FISP-MRF using long spiral readouts in conjunction with concomitant field-induced blurring mitigation. A wide operating regime with improved precision is feasible only after concomitant field mitigation.

Time-resolved MR fingerprinting for T2* signal extraction: MR fingerprinting meets echo planar time-resolved imaging

PURPOSE

This study leverages the echo planar time-resolved imaging (EPTI) concept in MR fingerprinting (MRF) framework for a new time-resolved MRF (TRMRF) approach, and explores its capability for fast simultaneous quantification of multiple MR parameters including T1, T2, T2*, proton density, off resonance, and B1 +.

METHODS

The proposed TRMRF method uses the concept of EPTI to track the signal change along the EPI echo train for T2* weighting with a k-t Poisson-based sampling order designed for acquisition. A two-dimensional decomposition algorithm was designed for the image reconstruction, enabling fast and precise subspace modeling. The accuracy of proposed method was evaluated by a T1/T2 phantom. The feasibility was demonstrated through 5 healthy volunteer brain studies.

RESULTS

In the phantom studies, T1, T2, and T2* maps of TRMRF correlated strongly with gold-standard methods. The concordance correlation coefficients are 0.9999, 0.9984 and 0.9978, and R2s are 0.9998, 0.9971, and 0.9983. In the in vivo studies, quantitative maps were acquired with 5 healthy volunteers. TRMRF was demonstrated to have comparable results with spiral MRF and gradient-echo EPTI. TRMRF scans using 16, 10, and 6 s per slice were also evaluated to demonstrate the capability of shorter scan times.

CONCLUSION

A new approach is proposed to exploit the advantage of EPTI in the MRF framework. We demonstrate in phantom and in vivo experiments that T1, T2, T2*, proton density, off resonance, and B1 + can be simultaneously quantified within 6 s/slice by TRMRF.

SuperMRF: deep robust reconstruction for highly accelerated magnetic resonance fingerprinting

Background

Magnetic resonance fingerprinting (MRF) is a rapid imaging technique for simultaneous mapping of multiple tissue properties such as T1 and T2 relaxation times. However, conventional pattern matching reconstruction and iterative low rank reconstruction methods may not take full advantage of the spatiotemporal content of MRF data and can require significant computational resources with long reconstruction times. Deep learning reconstruction using a three-dimensional (3D) convolutional neural network (CNN)-based method may enable high-quality, rapid MRF reconstruction. Evaluation of such proposed deep learning reconstruction methods for MRF is needed to clarify whether deep learning techniques adapted from other MR image reconstruction problems will yield benefits when employed in MRF applications. The objective of this study is to design and evaluate a novel deep learning framework (SuperMRF) that directly transforms undersampled parameter-weighted 3D Cartesian MRF data into quantitative T1 and T2 maps, bypassing traditional pattern-matching in MRF.

Methods

In contrast to conventional MRF where only the temporal evolution of each voxel is used for quantification, SuperMRF exploits both two-dimensional spatial and one-dimensional temporal information with a 3D CNN for reconstruction. Controlled simulation experiments were performed using reference parameter maps from in vivo knee scans of healthy volunteers. To evaluate the robustness to noise, we trained our network using clean data and tested it on simulated noisy data. Conventional inner product-based pattern matching and state-of-the-art iterative low rank reconstruction techniques were used for comparison. The performance of all methods was evaluated with respect to structural similarity index (SSIM), peak signal-to-noise ratio (PSNR), and normalized mean squared error (NMSE). Prospective real-world MRF scans were performed in four volunteer subjects using the trained network from simulations and cartilage and muscle T1 and T2 values were compared between conventional pattern matching, low rank reconstruction, and SuperMRF.

Results

SuperMRF estimated accurate T1 and T2 mapping in a highly accelerated scan (15× undersampling in k-space with a 20-fold reduction in the number of acquired MRF frames) with low error (NMSE of 5%) and high resemblance (SSIM of 94%) to reference quantitative maps. SuperMRF was observed to be superior to the conventional and low rank MRF reconstruction methods in terms of NMSE, SSIM, and robustness to noise. In prospective real-world data, SuperMRF provided comparable T1 and T2 maps as compared to low rank MRF. The only significantly different cartilage and muscle values in prospective data across the three reconstruction methods were those from conventional MRF T2.

Conclusions

Our results demonstrate that the proposed SuperMRF can achieve rapid, robust reconstruction with reduced frames in addition to k-space undersampling, outperforming the conventional and state-of-the-art reconstruction methods in simulation and providing comparable results to low rank reconstruction in prospective real-world subjects.

Deep graph embedding based on Laplacian eigenmaps for MR fingerprinting reconstruction

Magnetic resonance fingerprinting (MRF) is a promising technique for fast quantitative imaging of multiple tissue parameters. However, the highly undersampled schemes utilized in MRF typically lead to noticeable aliasing artifacts in reconstructed images. Existing model-based methods can mitigate aliasing artifacts and enhance reconstruction quality but suffer from long reconstruction times. In addition, data priors used in these methods, such as low-rank and total variation, make it challenging to incorporate non-local and non-linear redundancies in MRF data. Furthermore, existing deep learning-based methods for MRF often lack interpretability and struggle with the high computational overhead caused by the high dimensionality of MRF data. To address these issues, we introduce a novel deep graph embedding framework based on the Laplacian eigenmaps for improved MRF reconstruction. Our work first models the acquired high-dimensional MRF data and the corresponding parameter maps as graph data nodes. Then, we propose an MRF reconstruction framework based on the graph embedding framework, retaining intrinsic graph structures between parameter maps and MRF data. To improve the accuracy of the estimated graph structure and the computational efficiency of the proposed framework, we unroll the iterative optimization process into a deep neural network, incorporating a learned graph embedding module to adaptively learn the Laplacian eigenmaps. By introducing the graph embedding framework into the MRF reconstruction, the proposed method can effectively exploit non-local and non-linear correlations in MRF data. Numerical experiments demonstrate that our approach can reconstruct high-quality MRF data and multiple parameter maps within a significantly reduced computational cost.

Dynamic MRI with Locally Low-Rank Subspace Constraint: Towards 1-Second Temporal Resolution Aided by Deep Learning

MRI is the most effective method for screening high-risk breast cancer patients. While current exams primarily rely on the qualitative evaluation of morphological features before and after contrast administration and less on contrast kinetic information, the latest developments in acquisition protocols aim to combine both. However, balancing between spatial and temporal resolution poses a significant challenge in dynamic MRI. Here, we propose a radial MRI reconstruction framework for Dynamic Contrast Enhanced (DCE) imaging, which offers a joint solution to existing spatial and temporal MRI limitations. It leverages a locally low-rank (LLR) subspace model to represent spatially localized dynamics based on tissue information. Our framework demonstrated substantial improvement in CNR, noise reduction and enables a flexible temporal resolution, ranging from a few seconds to 1-second, aided by a neural network, resulting in images with reduced undersampling penalties. Finally, our reconstruction framework also shows potential benefits for head and neck, and brain MRI applications, making it a viable alternative for a range of DCE-MRI exams.

Generalizable, sequence-invariant deep learning image reconstruction for subspace-constrained quantitative MRI

PURPOSE

To develop a deep subspace learning network that can function across different pulse sequences.

METHODS

A contrast-invariant component-by-component (CBC) network structure was developed and compared against previously reported spatiotemporal multicomponent (MC) structure for reconstructing MR Multitasking images. A total of 130, 167, and 16 subjects were imaged using T1, T1-T2, and T1-T2-T 2 * $$ {\mathrm{T}}_2^{\ast } $$ -fat fraction (FF) mapping sequences, respectively. We compared CBC and MC networks in matched-sequence experiments (same sequence for training and testing), then examined their cross-sequence performance and generalizability by unmatched-sequence experiments (different sequences for training and testing). A "universal" CBC network was also evaluated using mixed-sequence training (combining data from all three sequences). Evaluation metrics included image normalized root mean squared error and Bland-Altman analyses of end-diastolic maps, both versus iteratively reconstructed references.

RESULTS

The proposed CBC showed significantly better normalized root mean squared error than MC in both matched-sequence and unmatched-sequence experiments (p < 0.001), fewer structural details in quantitative error maps, and tighter limits of agreement. CBC was more generalizable than MC (smaller performance loss; p = 0.006 in T1 and p < 0.001 in T1-T2 from matched-sequence testing to unmatched-sequence testing) and additionally allowed training of a single universal network to reconstruct images from any of the three pulse sequences. The mixed-sequence CBC network performed similarly to matched-sequence CBC in T1 (p = 0.178) and T1-T2 (p = 0121), where training data were plentiful, and performed better in T1-T2-T 2 * $$ {\mathrm{T}}_2^{\ast } $$ -FF (p < 0.001) where training data were scarce.

CONCLUSION

Contrast-invariant learning of spatial features rather than spatiotemporal features improves performance and generalizability, addresses data scarcity, and offers a pathway to universal supervised deep subspace learning.

Relaxometry and contrast-free cerebral microvascular quantification using balanced steady-state free precession MR fingerprinting

PURPOSE

This study proposes a novel, contrast-free Magnetic Resonance Fingerprinting (MRF) method using balanced Steady-State Free Precession (bSSFP) sequences for the quantification of cerebral blood volume (CBV), vessel radius (R), and relaxometry parameters (T  1 $$ {}_1 $$ , T  2 $$ {}_2 $$ , T  2 $$ {}_2 $$ *) in the brain.

METHODS

The technique leverages the sensitivity of bSSFP sequences to intra-voxel frequency distributions in both transient and steady-state regimes. A dictionary-matching process is employed, using simulations of realistic mouse microvascular networks to generate the MRF dictionary. The method is validated through in silico and in vivo experiments on six healthy subjects, comparing results with standard MRF methods and literature values.

RESULTS

The proposed method shows strong correlation and agreement with standard MRF methods for T  1 $$ {}_1 $$ and T  2 $$ {}_2 $$ values. High-resolution maps provide detailed visualizations of CBV and microvascular structures, highlighting differences in white matter (WM) and gray matter (GM) regions. The measured GM/WM ratio for CBV is 1.91, consistent with literature values.

CONCLUSION

This contrast-free bSSFP-based MRF method offers an new approach for quantifying CBV, vessel radius, and relaxometry parameters. Further validation against DSC imaging and clinical studies in pathological conditions is warranted to confirm its clinical utility.

Acceleration of Simultaneous Multislice Magnetic Resonance Fingerprinting With Spatiotemporal Convolutional Neural Network

Magnetic Resonance Fingerprinting (MRF) can be accelerated with simultaneous multislice (SMS) imaging for joint T1 and T2 quantification. However, the high inter-slice and in-plane acceleration in SMS-MRF causes severe aliasing artifacts, limiting the multiband (MB) factors to typically 2 or 3. Deep learning has demonstrated superior performance compared to the conventional dictionary matching approach for single-slice MRF, but its effectiveness in SMS-MRF remains unexplored. In this paper, we introduced a new deep learning approach with decoupled spatiotemporal feature learning for SMS-MRF to achieve high MB factors for accurate and volumetric T1 and T2 quantification in neuroimaging. The proposed method leverages information from both spatial and temporal domains to mitigate the significant aliasing in SMS-MRF. Neural networks, trained using either acquired SMS-MRF data or simulated data generated from single-slice MRF acquisitions, were evaluated. The performance was further compared with both dictionary matching and a deep learning approach based on residual channel attention U-Net. Experimental results demonstrated that the proposed method, trained with acquired SMS-MRF data, achieves the best performance in brain T1 and T2 quantification, outperforming dictionary matching and residual channel attention U-Net. With a MB factor of 4, rapid T1 and T2 mapping was achieved with 1.5 s per slice for quantitative brain imaging.

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

OBJECTIVE

We extend the traditional framework for estimating subspace bases in quantitative MRI 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.

METHODS

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.

RESULTS

Compared to the traditional SVD approach, the proposed method better preserves the CRB across all biophysical parameters with minimal 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.

CONCLUSION

The proposed method permits subspace reconstruction with a more compact basis, thereby offering significant computational savings.

SIGNIFICANCE

Efficient subspace reconstruction facilitates the validation and translation of advanced quantitative MRI techniques, e.g., magnetization transfer and diffusion.

3D MR fingerprinting for dynamic contrast-enhanced imaging of whole mouse brain

PURPOSE

Quantitative MRI enables direct quantification of contrast agent concentrations in contrast-enhanced scans. However, the lengthy scan times required by conventional methods are inadequate for tracking contrast agent transport dynamically in mouse brain. We developed a 3D MR fingerprinting (MRF) method for simultaneous T1 and T2 mapping across the whole mouse brain with 4.3-min temporal resolution.

METHOD

We designed a 3D MRF sequence with variable acquisition segment lengths and magnetization preparations on a 9.4T preclinical MRI scanner. Model-based reconstruction approaches were employed to improve the accuracy and speed of MRF acquisition. The method's accuracy for T1 and T2 measurements was validated in vitro, while its repeatability of T1 and T2 measurements was evaluated in vivo (n = 3). The utility of the 3D MRF sequence for dynamic tracking of intracisternally infused gadolinium-diethylenetriamine pentaacetic acid (Gd-DTPA) in the whole mouse brain was demonstrated (n = 5).

RESULTS

Phantom studies confirmed accurate T1 and T2 measurements by 3D MRF with an undersampling factor of up to 48. Dynamic contrast-enhanced MRF scans achieved a spatial resolution of 192 × 192 × 500 μm3 and a temporal resolution of 4.3 min, allowing for the analysis and comparison of dynamic changes in concentration and transport kinetics of intracisternally infused Gd-DTPA across brain regions. The sequence also enabled highly repeatable, high-resolution T1 and T2 mapping of the whole mouse brain (192 × 192 × 250 μm3) in 30 min.

CONCLUSION

We present the first dynamic and multi-parametric approach for quantitatively tracking contrast agent transport in the mouse brain using 3D MRF.

Contrast-Optimized Basis Functions for Self-Navigated Motion Correction in Quantitative MRI

Purpose

The long scan times of quantitative MRI techniques make motion artifacts more likely. For MR-Fingerprinting-like approaches, this problem can be addressed with self-navigated retrospective motion correction based on reconstructions in a singular value decomposition (SVD) subspace. However, the SVD promotes high signal intensity in all tissues, which limits the contrast between tissue types and ultimately reduces the accuracy of registration. The purpose of this paper is to rotate the subspace for maximum contrast between two types of tissue and improve the accuracy of motion estimates.

Methods

A subspace is derived that promotes contrasts between brain parenchyma and CSF, achieved through the generalized eigendecomposition of mean autocorrelation matrices, followed by a Gram-Schmidt process to maintain orthogonality.We tested our motion correction method on 85 scans with varying motion levels, acquired with a 3D hybrid-state sequence optimized for quantitative magnetization transfer imaging.

Results

A comparative analysis shows that the contrast-optimized basis significantly improve the parenchyma-CSF contrast, leading to smoother motion estimates and reduced artifacts in the quantitative maps.

Conclusion

The proposed contrast-optimized subspace improves the accuracy of the motion estimation.

Simultaneous Bilateral T1, T2, and T1ρ Relaxation Mapping of Hip Joint With 3D-MRI Fingerprinting

BACKGROUND

Three-dimensional MR fingerprinting (3D-MRF) has been increasingly used to assess cartilage degeneration, particularly in the knee joint, by looking into multiple relaxation parameters. A comparable 3D-MRF approach can be adapted to assess cartilage degeneration for the hip joint, with changes to accommodate specific challenges of hip joint imaging.

PURPOSE

To demonstrate the feasibility and repeatability of 3D-MRF in the bilateral hip jointly we map proton density (PD), T1, T2, T1ρ, and ∆B1+ in clinically feasible scan times.

STUDY TYPE

Prospective.

SUBJECTS

Eight healthy subjects, three patients with mild osteoarthritis (OA), and one of the OA patients had femoral acetabular impingement (FAI). A National Institute of Standards and Technology/International Society for Magnetic Resonance in Medicine (NIST/ISMRM) system phantom was also used.

FIELD STRENGTH/SEQUENCE

3 T, 3D-MRF sequence for bilateral hip joint mapping. Reference sequences include Volume Interpolated Breath-hold Examination (VIBE) for T1 mapping, and magnetization-prepared fast low-angle shot (TFL) for T2 and T1ρ mapping.

ASSESSMENT

The signal-to-noise ratio (SNR), repeatability, scan time, and accuracy of T1, T2, and T1ρ maps of 3D-MRF sequence were evaluated on a NIST/ISMRM phantom and human subjects. Differences in the parametric maps between OA and healthy subjects were assessed.

STATISTICAL TESTS

Regression, Bland-Altman, Kruskal-Wallis, and Wilcoxon tests were used to assess for accuracy, repeatability, and subregional variation. The P-value <0.05 indicated statistically significant.

RESULTS

A 3D-MRF sequence sensitive to PD, T1, T2, T1ρ, and ∆B1+ within 15 minutes, achieving high SNR and low test-retest coefficient of variance (T1: 3.36%, T2: 3.99%, T1ρ: 5.93%). Mild hip OA patients, including one with mild OA and FAI, showed elevation of 29.4 ± 9% (T2) and 32.4 ± 4.4% (T1ρ) in femoral lateral compartment of the hip joint compared to healthy controls.

DATA CONCLUSION

3D-MRF may be a feasible approach for simultaneous, quantitative mapping of bilateral hip joint cartilage in healthy and mild OA patients.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 1.

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.

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.

3D MR Fingerprinting for Dynamic Contrast-Enhanced Imaging of Whole Mouse Brain

Purpose

Quantitative MRI enables direct quantification of contrast agent concentrations in contrast-enhanced scans. However, the lengthy scan times required by conventional methods are inadequate for tracking contrast agent transport dynamically in mouse brain. We developed a 3D MR fingerprinting (MRF) method for simultaneous T1 and T2 mapping across the whole mouse brain with 4.3-min temporal resolution.

Method

We designed a 3D MRF sequence with variable acquisition segment lengths and magnetization preparations on a 9.4T preclinical MRI scanner. Model-based reconstruction approaches were employed to improve the accuracy and speed of MRF acquisition. The method's accuracy for T1 and T2 measurements was validated in vitro, while its repeatability of T1 and T2 measurements was evaluated in vivo (n=3). The utility of the 3D MRF sequence for dynamic tracking of intracisternally infused Gd-DTPA in the whole mouse brain was demonstrated (n=5).

Results

Phantom studies confirmed accurate T1 and T2 measurements by 3D MRF with an undersampling factor up to 48. Dynamic contrast-enhanced (DCE) MRF scans achieved a spatial resolution of 192 ✕ 192 ✕ 500 μm3 and a temporal resolution of 4.3 min, allowing for the analysis and comparison of dynamic changes in concentration and transport kinetics of intracisternally infused Gd-DTPA across brain regions. The sequence also enabled highly repeatable, high-resolution T1 and T2 mapping of the whole mouse brain (192 ✕ 192 ✕ 250 μm3) in 30 min.

Conclusion

We present the first dynamic and multi-parametric approach for quantitatively tracking contrast agent transport in the mouse brain using 3D MRF.

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.

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.

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 T2 -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 T2 -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.

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.

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.

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.

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 T1 and T2 values obtained with MRF have been validated with phantom studies as well as in normal volunteers and patients. Studies have shown that MRF-derived T1 and T2 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.

An efficient approach to optimal experimental design for magnetic resonance fingerprinting with B-splines

PURPOSE

To introduce a computationally efficient approach to optimizing the data acquisition parameters of MR Fingerprinting experiments with the Cramér-Rao bound.

METHODS

This paper presents a new approach to the optimal experimental design (OED) problem for MR Fingerprinting, which leverages an early observation that the optimized data acquisition parameters of MR Fingerprinting experiments are highly structured. Specifically, the proposed approach captures the desired structure by representing the sequences of data acquisition parameters with a special class of piecewise polynomials known as B-splines. This incorporates low-dimensional spline subspace constraints into the OED problem, which significantly reduces the search space of the problem, thereby improving the computational efficiency. With the rich B-spline representations, the proposed approach also allows for incorporating prior knowledge on the structure of different acquisition parameters, which facilitates the experimental design.

RESULTS

The effectiveness of the proposed approach was evaluated using numerical simulations, phantom experiments, and in vivo experiments. The proposed approach achieves a two-order-of-magnitude improvement of the computational efficiency over the state-of-the-art approaches, while providing a comparable signal-to-noise ratio efficiency benefit. It enables an optimal experimental design problem for MR Fingerprinting with a typical acquisition length to be solved in approximately 1 min.

CONCLUSIONS

The proposed approach significantly improves the computational efficiency of the optimal experimental design for MR Fingerprinting, which enhances its practical utility for a variety of quantitative MRI applications.

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.

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 T1 and T2 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 T1 and T2 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 T1 and T2 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.

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.