Multicompartment Magnetic Resonance Fingerprinting

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.

White matter changes measured by multi-component MR Fingerprinting in multiple sclerosis

T2-hyperintense lesions are the key imaging marker of multiple sclerosis (MS). Previous studies have shown that the white matter surrounding such lesions is often also affected by MS. Our aim was to develop a new method to visualize and quantify the extent of white matter tissue changes in MS based on relaxometry properties. We applied a fast, multi-parametric quantitative MRI approach and used a multi-component MR Fingerprinting (MC-MRF) analysis. We assessed the differences in the MRF component representing prolongedrelaxation time between patients with MS and controls and studied the relation between this component's volume and structural white matter damage identified on FLAIR MRI scans in patients with MS. A total of 48 MS patients at two different sites and 12 healthy controls were scanned with FLAIR and MRF-EPI MRI scans. MRF scans were analyzed with a joint-sparsity multi-component analysis to obtain magnetization fraction maps of different components, representing tissues such as myelin water, white matter, gray matter and cerebrospinal fluid. In the MS patients, an additional component was identified with increased transverse relaxation times compared to the white matter, likely representing changes in free water content. Patients with MS had a higher volume of the long- component in the white matter of the brain compared to healthy controls (B (95%-CI) = 0.004 (0.0006-0.008), p = 0.02). Furthermore, this MRF component had a moderate correlation (correlation coefficient R 0.47) with visible structural white matter changes on the FLAIR scans. Also, the component was found to be more extensive compared to structural white matter changes in 73% of MS patients. In conclusion, our MRF acquisition and analysis captured white matter tissue changes in MS patients compared to controls. In patients these tissue changes were more extensive compared to visually detectable white matter changes on FLAIR scans. Our method provides a novel way to quantify the extent of white matter changes in MS patients, which is underestimated using only conventional clinical MRI scans.

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.

Accuracy and repeatability of joint sparsity multi-component estimation in MR Fingerprinting

MR fingerprinting (MRF) is a promising method for quantitative characterization of tissues. Often, voxel-wise measurements are made, assuming a single tissue-type per voxel. Alternatively, the Sparsity Promoting Iterative Joint Non-negative least squares Multi-Component MRF method (SPIJN-MRF) facilitates tissue parameter estimation for identified components as well as partial volume segmentations. The aim of this paper was to evaluate the accuracy and repeatability of the SPIJN-MRF parameter estimations and partial volume segmentations. This was done (1) through numerical simulations based on the BrainWeb phantoms and (2) using in vivo acquired MRF data from 5 subjects that were scanned on the same week-day for 8 consecutive weeks. The partial volume segmentations of the SPIJN-MRF method were compared to those obtained by two conventional methods: SPM12 and FSL. SPIJN-MRF showed higher accuracy in simulations in comparison to FSL- and SPM12-based segmentations: Fuzzy Tanimoto Coefficients (FTC) comparing these segmentations and Brainweb references were higher than 0.95 for SPIJN-MRF in all the tissues and between 0.6 and 0.7 for SPM12 and FSL in white and gray matter and between 0.5 and 0.6 in CSF. For the in vivo MRF data, the estimated relaxation times were in line with literature and minimal variation was observed. Furthermore, the coefficient of variation (CoV) for estimated tissue volumes with SPIJN-MRF were 10.5% for the myelin water, 6.0% for the white matter, 5.6% for the gray matter, 4.6% for the CSF and 1.1% for the total brain volume. CoVs for CSF and total brain volume measured on the scanned data for SPIJN-MRF were in line with those obtained with SPM12 and FSL. The CoVs for white and gray matter volumes were distinctively higher for SPIJN-MRF than those measured with SPM12 and FSL. In conclusion, the use of SPIJN-MRF provides accurate and precise tissue relaxation parameter estimations taking into account intrinsic partial volume effects. It facilitates obtaining tissue fraction maps of prevalent tissues including myelin water which can be relevant for evaluating diseases affecting the white matter.

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.

Data-driven myelin water imaging based on T1 and T2 relaxometry

Long acquisition times preclude the application of multiecho spin echo (MESE) sequences for myelin water fraction (MWF) mapping in daily clinical practice. In search of alternative methods, previous studies of interest explored the biophysical modeling of MWF from measurements of different tissue properties that can be obtained in scan times shorter than those required for the MESE. In this work, a novel data-driven estimation of MWF maps from fast relaxometry measurements is proposed and investigated. T₁ and T₂ relaxometry maps were acquired in a cohort of 20 healthy subjects along with a conventional MESE sequence. Whole-brain quantitative mapping was achieved with a fast protocol in 6 min 24 s. Reference MWF maps were derived from the MESE sequence (TA = 11 min 17 s) and their data-driven estimation from relaxometry measurements was investigated using three different modeling strategies: two general linear models (GLMs) with linear and quadratic regressors, respectively; a random forest regression model; and two deep neural network architectures, a U-Net and a conditional generative adversarial network (cGAN). Models were validated using a 10-fold crossvalidation. The resulting maps were visually and quantitatively compared by computing the root mean squared error (RMSE) between the estimated and reference MWF maps, the intraclass correlation coefficients (ICCs) between corresponding MWF values in different brain regions, and by performing Bland-Altman analysis. Qualitatively, the estimated maps appear to generally provide a similar, yet more blurred MWF contrast in comparison with the reference, with the cGAN model best capturing MWF variabilities in small structures. By estimating the average adjusted coefficient of determination of the GLM with quadratic regressors, we showed that 87% of the variability in the MWF values can be explained by relaxation times alone. Further quantitative analysis showed an average RMSE smaller than 0.1% for all methods. The ICC was greater than 0.81 for all methods, and the bias smaller than 2.19%. It was concluded that this work confirms the notion that relaxometry parameters contain a large part of the information on myelin water and that MWF maps can be generated from T₁ /T₂ data with minimal error. Among the investigated modeling approaches, the cGAN provided maps with the best trade-off between accuracy and blurriness. Fast relaxometry, like the 6 min 24 s whole-brain protocol used in this work in conjunction with machine learning, may thus have the potential to replace time-consuming MESE acquisitions.

Correcting for imaging gradients-related bias of T2 relaxation times at high-resolution MRI

Purpose

High‐resolution animal imaging is an integral part of preclinical drug development and the investigation of diseases' pathophysiology. Quantitative mapping of T₂ relaxation times (qT₂) is a valuable tool for both preclinical and research applications, providing high sensitivity to subtle tissue pathologies. High‐resolution T₂ mapping, however, suffers from severe underestimation of T₂ values due to molecular diffusion. This affects both single‐echo and multi‐echo spin echo (SSE and MESE), on top of the well‐known contamination of MESE signals by stimulated echoes, and especially on high‐field and preclinical scanners in which high imaging gradients are used in comparison to clinical scanners.

Methods

Diffusion bias due to imaging gradients was analyzed by quantifying the effective b‐value for each coherence pathway in SSE and MESE protocols, and incorporating this information in a joint T₂‐diffusion reconstruction algorithm. Validation was done on phantoms and in vivo mouse brain using a 9.4T and a 7T MRI scanner.

Results

Underestimation of T₂ values due to strong imaging gradients can reach up to 70%, depending on scan parameters and on the sample's diffusion coefficient. The algorithm presented here produced T₂ values that agreed with reference spectroscopic measurements, were reproducible across scan settings, and reduced the average bias of T₂ values from −33.5 ± 20.5% to −0.1 ± 3.6%.

Conclusions

A new joint T₂‐diffusion reconstruction algorithm is able to negate imaging gradient–related underestimation of T₂ values, leading to reliable mapping of T₂ values at high resolutions.

Click here for author‐reader discussions

The Road Toward Reproducibility of Parametric Mapping of the Heart: A Technical Review

Parametric mapping of the heart has become an essential part of many cardiovascular magnetic resonance imaging exams, and is used for tissue characterization and diagnosis in a broad range of cardiovascular diseases. These pulse sequences are used to quantify the myocardial T₁, T₂, T2*, and T_1ρ relaxation times, which are unique surrogate indices of fibrosis, edema and iron deposition that can be used to monitor a disease over time or to compare patients to one another. Parametric mapping is now well-accepted in the clinical setting, but its wider dissemination is hindered by limited inter-center reproducibility and relatively long acquisition times. Recently, several new parametric mapping techniques have appeared that address both of these problems, but substantial hurdles remain for widespread clinical adoption. This review serves both as a primer for newcomers to the field of parametric mapping and as a technical update for those already well at home in it. It aims to establish what is currently needed to improve the reproducibility of parametric mapping of the heart. To this end, we first give an overview of the metrics by which a mapping technique can be assessed, such as bias and variability, as well as the basic physics behind the relaxation times themselves and what their relevance is in the prospect of myocardial tissue characterization. This is followed by a summary of routine mapping techniques and their variations. The problems in reproducibility and the sources of bias and variability of these techniques are reviewed. Subsequently, novel fast, whole-heart, and multi-parametric techniques and their merits are treated in the light of their reproducibility. This includes state of the art segmentation techniques applied to parametric maps, and how artificial intelligence is being harnessed to solve this long-standing conundrum. We finish up by sketching an outlook on the road toward inter-center reproducibility, and what to expect in the future.

(An overview of) Synergistic reconstruction for multimodality/multichannel imaging methods

Imaging is omnipresent in modern society with imaging devices based on a zoo of physical principles, probing a specimen across different wavelengths, energies and time. Recent years have seen a change in the imaging landscape with more and more imaging devices combining that which previously was used separately. Motivated by these hardware developments, an ever increasing set of mathematical ideas is appearing regarding how data from different imaging modalities or channels can be synergistically combined in the image reconstruction process, exploiting structural and/or functional correlations between the multiple images. Here we review these developments, give pointers to important challenges and provide an outlook as to how the field may develop in the forthcoming years. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.

Whole-Brain Imaging of Subvoxel T1-Diffusion Correlation Spectra in Human Subjects

T1 relaxation and water mobility generate eloquent MRI tissue contrasts with great diagnostic value in many neuroradiological applications. However, conventional methods do not adequately quantify the microscopic heterogeneity of these important biophysical properties within a voxel, and therefore have limited biological specificity. We describe a new correlation spectroscopic (CS) MRI method for measuring how T1 and mean diffusivity (MD) co-vary in microscopic tissue environments. We develop a clinical pulse sequence that combines inversion recovery (IR) with single-shot isotropic diffusion encoding (IDE) to efficiently acquire whole-brain MRIs with a wide range of joint T1-MD weightings. Unlike conventional diffusion encoding, the IDE preparation ensures that all subvoxel water pools are weighted by their MDs regardless of the sizes, shapes, and orientations of their corresponding microscopic diffusion tensors. Accordingly, IR-IDE measurements are well-suited for model-free, quantitative spectroscopic analysis of microscopic water pools. Using numerical simulations, phantom experiments, and data from healthy volunteers we demonstrate how IR-IDE MRIs can be processed to reconstruct maps of two-dimensional joint probability density functions, i.e., correlation spectra, of subvoxel T1-MD values. In vivo T1-MD spectra show distinct cerebrospinal fluid and parenchymal tissue components specific to white matter, cortical gray matter, basal ganglia, and myelinated fiber pathways, suggesting the potential for improved biological specificity. The one-dimensional marginal distributions derived from the T1-MD correlation spectra agree well with results from other relaxation spectroscopic and quantitative MRI studies, validating the T1-MD contrast encoding and the spectral reconstruction. Mapping subvoxel T1-diffusion correlations in patient populations may provide a more nuanced, comprehensive, sensitive, and specific neuroradiological assessment of the non-specific changes seen on fluid-attenuated inversion recovery (FLAIR) and diffusion-weighted MRIs (DWIs) in cancer, ischemic stroke, or brain injury.

3D sodium (23 Na) magnetic resonance fingerprinting for time-efficient relaxometric mapping

PURPOSE

To develop a framework for 3D sodium (²³ Na) MR fingerprinting (MRF), based on irreducible spherical tensor operators with tailored flip angle (FA) pattern and time-efficient data acquisition for simultaneous quantification of T₁ , T 2 l ∗ , T 2 s ∗ , and T 2 ∗ in addition to ΔB₀ .

METHODS

²³ Na-MRF was implemented in a 3D sequence and irreducible spherical tensor operators were exploited in the simulations. Furthermore, the Cramér Rao lower bound was used to optimize the flip angle pattern. A combination of single and double echo readouts was implemented to increase the readout efficiency. A study was conducted to compare results in a multicompartment phantom acquired with MRF and reference methods. Finally, the relaxation times in the human brain were measured in four healthy volunteers.

RESULTS

Phantom experiments revealed a mean difference of 1.0% between relaxation times acquired with MRF and results determined with the reference methods. Simultaneous quantification of the longitudinal and transverse relaxation times in the human brain was possible within 32 min using 3D ²³ Na-MRF with a nominal resolution of (5 mm)³ . In vivo measurements in four volunteers yielded average relaxation times of: T_1,brain = (35.0 ± 3.2) ms, T 2 l , brain ∗ = (29.3 ± 3.8) ms and T 2 s , brain ∗ = (5.5 ± 1.3) ms in brain tissue, whereas T_1,CSF = (61.9 ± 2.8) ms and T 2 , CSF ∗ = (46.3 ± 4.5) ms was found in cerebrospinal fluid.

CONCLUSION

The feasibility of in vivo 3D relaxometric sodium mapping within roughly ½ h is demonstrated using MRF in the human brain, moving sodium relaxometric mapping toward clinically relevant measurement times.

Magnetic resonance fingerprinting: an overview

Magnetic resonance fingerprinting (MRF) is an evolving quantitative MRI framework consisting of unique data acquisition, processing, visualization, and interpretation steps. MRF is capable of simultaneously producing multiple high-resolution property maps including T1, T2, M0, ADC, and T2* measurements. While a relatively new technology, MRF has undergone rapid development for a variety of clinical applications from brain tumor characterization and epilepsy imaging to characterization of prostate cancer, cardiac imaging, among others. This paper will provide a brief overview of current state of MRF technology including highlights of technical and clinical advances. We will conclude with a brief discussion of the challenges that need to be overcome to establish MRF as a quantitative imaging biomarker.

Magnetic resonance fingerprinting: from evolution to clinical applications

In 2013, Magnetic Resonance Fingerprinting (MRF) emerged as a method for fast, quantitative Magnetic Resonance Imaging. This paper reviews the current status of MRF up to early 2020 and aims to highlight the advantages MRF can offer medical imaging professionals. By acquiring scan data as pseudorandom samples, MRF elicits a unique signal evolution, or 'fingerprint', from each tissue type. It matches 'randomised' free induction decay acquisitions against pre-computed simulated tissue responses to generate a set of quantitative images of T1 , T2 and proton density (PD) with co-registered voxels, rather than as traditional relative T1 - and T2 -weighted images. MRF numeric pixel values retain accuracy and reproducibility between 2% and 8%. MRF acquisition is robust to strong undersampling of k-space. Scan sequences have been optimised to suppress sub-sampling artefacts, while artificial intelligence and machine learning techniques have been employed to increase matching speed and precision. MRF promises improved patient comfort with reduced scan times and fewer image artefacts. Quantitative MRF data could be used to define population-wide numeric biomarkers that classify normal versus diseased tissue. Certification of clinical centres for MRF scan repeatability would permit numeric comparison of sequential images for any individual patient and the pooling of multiple patient images across large, cross-site imaging studies. MRF has to date shown promising results in early clinical trials, demonstrating reliable differentiation between malignant and benign prostate conditions, and normal and sclerotic hippocampal tissue. MRF is now undergoing small-scale trials at several sites across the world; moving it closer to routine clinical application.

Multidimensional correlation spectroscopic imaging of exponential decays: From theoretical principles to in vivo human applications

Multiexponential modeling of relaxation or diffusion MR signal decays is a popular approach for estimating and spatially mapping different microstructural tissue compartments. While this approach can be quite powerful, it is also limited by the fact that one-dimensional multiexponential modeling is an ill-posed inverse problem with substantial ambiguities. In this article, we present an overview of a recent multidimensional correlation spectroscopic imaging approach to this problem. This approach helps to alleviate ill-posedness by making advantageous use of multidimensional contrast encoding (e.g., 2D diffusion-relaxation encoding or 2D relaxation-relaxation encoding) combined with a regularized spatial-spectral estimation procedure. Theoretical calculations, simulations, and experimental results are used to illustrate the benefits of this approach relative to classical methods. In addition, we demonstrate an initial proof-of-principle application of this kind of approach to in vivo human MRI experiments.

Sparse Recovery Beyond Compressed Sensing: Separable Nonlinear Inverse Problems

Extracting information from nonlinear measurements is a fundamental challenge in data analysis. In this work, we consider separable inverse problems, where the data are modeled as a linear combination of functions that depend nonlinearly on certain parameters of interest. These parameters may represent neuronal activity in a human brain, frequencies of electromagnetic waves, fluorescent probes in a cell, or magnetic relaxation times of biological tissues. Separable nonlinear inverse problems can be reformulated as underdetermined sparse-recovery problems, and solved using convex programming. This approach has had empirical success in a variety of domains, from geophysics to medical imaging, but lacks a theoretical justification. In particular, compressed-sensing theory does not apply, because the measurement operators are deterministic and violate incoherence conditions such as the restricted-isometry property. Our main contribution is a theory for sparse recovery adapted to deterministic settings. We show that convex programming succeeds in recovering the parameters of interest, as long as their values are sufficiently distinct with respect to the correlation structure of the measurement operator. The theoretical results are illustrated through numerical experiments for two applications: heat-source localization and estimation of brain activity from electroencephalography data.

Magnetic resonance fingerprinting review part 2: Technique and directions

Magnetic resonance fingerprinting (MRF) is a general framework to quantify multiple MR-sensitive tissue properties with a single acquisition. There have been numerous advances in MRF in the years since its inception. In this work we highlight some of the recent technical developments in MRF, focusing on sequence optimization, modifications for reconstruction and pattern matching, new methods for partial volume analysis, and applications of machine and deep learning. Level of Evidence: 2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:993-1007.

Free-running cardiac magnetic resonance fingerprinting: Joint T1/T2 map and Cine imaging

PURPOSE

To develop and evaluate a novel non-ECG triggered 2D magnetic resonance fingerprinting (MRF) sequence allowing for simultaneous myocardial T₁ and T₂ mapping and cardiac Cine imaging.

METHODS

Cardiac MRF (cMRF) has been recently proposed to provide joint T₁/T₂ myocardial mapping by triggering the acquisition to mid-diastole and relying on a subject-dependent dictionary of MR signal evolutions to generate the maps. In this work, we propose a novel "free-running" (non-ECG triggered) cMRF framework for simultaneous myocardial T₁ and T₂ mapping and cardiac Cine imaging in a single scan. Free-running cMRF is based on a transient state bSSFP acquisition with tiny golden angle radial readouts, varying flip angle and multiple adiabatic inversion pulses. The acquired data is retrospectively gated into several cardiac phases, which are reconstructed with an approach that combines parallel imaging, low rank modelling and patch-based high-order tensor regularization. Free-running cMRF was evaluated in a standardized phantom and ten healthy subjects. Comparison with reference spin-echo, MOLLI, SASHA, T₂-GRASE and Cine was performed.

RESULTS

T₁ and T₂ values obtained with the proposed approach were in good agreement with reference phantom values (ICC(A,1) > 0.99). Reported values for myocardium septum T₁ were 1043 ± 48 ms, 1150 ± 100 ms and 1160 ± 79 ms for MOLLI, SASHA and free-running cMRF respectively and for T₂ of 51.7 ± 4.1 ms and 44.6 ± 4.1 ms for T₂-GRASE and free-running cMRF respectively. Good agreement was observed between free-running cMRF and conventional Cine 2D ejection fraction (bias = -0.83%).

CONCLUSION

The proposed free-running cardiac MRF approach allows for simultaneous assessment of myocardial T₁ and T₂ and Cine imaging in a single scan.

Fast multi-component analysis using a joint sparsity constraint for MR fingerprinting

PURPOSE

To develop an efficient algorithm for multi-component analysis of magnetic resonance fingerprinting (MRF) data without making a priori assumptions about the exact number of tissues or their relaxation properties.

METHODS

Different tissues or components within a voxel are potentially separable in MRF because of their distinct signal evolutions. The observed signal evolution in each voxel can be described as a linear combination of the signals for each component with a non-negative weight. An assumption that only a small number of components are present in the measured field of view is usually imposed in the interpretation of multi-component data. In this work, a joint sparsity constraint is introduced to utilize this additional prior knowledge in the multi-component analysis of MRF data. A new algorithm combining joint sparsity and non-negativity constraints is proposed and compared to state-of-the-art multi-component MRF approaches in simulations and brain MRF scans of 11 healthy volunteers.

RESULTS

Simulations and in vivo measurements show reduced noise in the estimated tissue fraction maps compared to previously proposed methods. Applying the proposed algorithm to the brain data resulted in 4 or 5 components, which could be attributed to different brain structures, consistent with previous multi-component MRF publications.

CONCLUSIONS

The proposed algorithm is faster than previously proposed methods for multi-component MRF and the simulations suggest improved accuracy and precision of the estimated weights. The results are easier to interpret compared to voxel-wise methods, which combined with the improved speed is an important step toward clinical evaluation of multi-component MRF.

Water-fat Dixon cardiac magnetic resonance fingerprinting

Purpose

Cardiac magnetic resonance fingerprinting (cMRF) has been recently introduced to simultaneously provide T₁, T₂, and M₀ maps. Here, we develop a 3‐point Dixon‐cMRF approach to enable simultaneous water specific T₁, T₂, and M₀ mapping of the heart and fat fraction (FF) estimation in a single breath‐hold scan.

Methods

Dixon‐cMRF is achieved by combining cMRF with several innovations that were previously introduced for other applications, including a 3‐echo GRE acquisition with golden angle radial readout and a high‐dimensional low‐rank tensor constrained reconstruction to recover the highly undersampled time series images for each echo. Water–fat separation of the Dixon‐cMRF time series is performed to allow for water‐ and fat‐specific T₁, T₂, and M₀ estimation, whereas FF estimation is extracted from the M₀ maps. Dixon‐cMRF was evaluated in a standardized T₁–T₂ phantom, in a water–fat phantom, and in healthy subjects in comparison to current clinical standards: MOLLI, SASHA, T₂‐GRASE, and 6‐point Dixon proton density FF (PDFF) mapping.

Results

Dixon‐cMRF water T₁ and T₂ maps showed good agreement with reference T₁ and T₂ mapping techniques (R² > 0.99 and maximum normalized RMSE ~5%) in a standardized phantom. Good agreement was also observed between Dixon‐cMRF FF and reference PDFF (R² > 0.99) and between Dixon‐cMRF water T₁ and T₂ and water selective T₁ and T₂ maps (R² > 0.99) in a water–fat phantom. In vivo Dixon‐cMRF water T₁ values were in good agreement with MOLLI and water T₂ values were slightly underestimated when compared to T₂‐GRASE. Average myocardium septal T₁ values were 1129 ± 38 ms, 1026 ± 28 ms, and 1045 ± 32 ms for SASHA, MOLLI, and the proposed water Dixon‐cMRF. Average T₂ values were 51.7 ± 2.2 ms and 42.8 ± 2.6 ms for T₂‐GRASE and water Dixon‐cMRF, respectively. Dixon‐cMRF FF maps showed good agreement with in vivo PDFF measurements (R² > 0.98) and average FF in the septum was measured at 1.3%.

Conclusion

The proposed Dixon‐cMRF allows to simultaneously quantify myocardial water T₁, water T₂, and FF in a single breath‐hold scan, enabling multi‐parametric T₁, T₂, and fat characterization. Moreover, reduced T₁ and T₂ quantification bias caused by water–fat partial volume was demonstrated in phantom experiments.

Cardiac Magnetic Resonance Fingerprinting: Technical Developments and Initial Clinical Validation

PURPOSE OF REVIEW

Magnetic resonance imaging (MRI) has enabled non-invasive myocardial tissue characterization in a wide range of cardiovascular diseases by quantifying several tissue specific parameters such as T1, T2, and T2* relaxation times. Simultaneous assessment of these parameters has recently gained interest to potentially improve diagnostic accuracy and enable further understanding of the underlying disease. However, these quantitative maps are usually acquired sequentially and are not necessarily co-registered, making multi-parametric analysis challenging. Magnetic resonance fingerprinting (MRF) has been recently introduced to unify and streamline parametric mapping into a single simultaneous, multi-parametric, fully co-registered, and efficient scan. Feasibility of cardiac MRF has been demonstrated and initial clinical validation studies are ongoing. Provide an overview of the cardiac MRF framework, recent technical developments and initial undergoing clinical validation.

RECENT FINDINGS

Cardiac MRF has enabled the acquisition of co-registered T1 and T2 maps in a single, efficient scan. Initial results demonstrate feasibility of cardiac MRF in healthy subjects and small patient cohorts. Current in vivo results show a small bias and comparable precision in T1 and T2 with respect to conventional clinical parametric mapping approaches. This bias may be explained by several confounding factors such as magnetization transfer and field inhomogeneities, which are currently not included in the cardiac MRF model. Initial clinical validation for cardiac MRF has demonstrated good reproducibility in healthy subjects and heart transplant patients, reduced artifacts in inflammatory cardiomyopathy patients and good differentiation between hypertrophic cardiomyopathy and healthy controls. Cardiac MRF has emerged as a novel technique for simultaneous, multi-parametric, and co-registered mapping of different tissue parameters. Initial efforts have focused on enabling T1, T2, and fat quantification; however this approach has the potential of enabling quantification of several other parameters (such as T2*, diffusion, perfusion, and flow) from a single scan. Initial results in healthy subjects and patients are promising, thus further clinical validation is now warranted.

Sparsity and locally low rank regularization for MR fingerprinting

PURPOSE

Develop a sparse and locally low rank (LLR) regularized reconstruction to accelerate MR fingerprinting (MRF).

METHODS

Recent works have introduced low rank reconstructions to MRF, based on temporal compression operators learned from the MRF dictionary. In other MR applications, LLR regularization has been introduced to exploit temporal redundancy in local regions of the image. Here, we propose to include spatial sparsity and LLR regularization terms in the MRF reconstruction. This approach, so called SLLR-MRF, further reduces aliasing in the time-point images and enables higher acceleration factors. The proposed approach was evaluated in simulations, T1 /T2 phantom acquisition, and in vivo brain acquisitions in 5 healthy subjects with different undersampling factors. Acceleration was also used in vivo to enable acquisitions with higher in-plane spatial resolution in comparable scan time.

RESULTS

Simulations, phantom, and in vivo results show that low rank MRF reconstructions with high acceleration factors (<875 time-point images, 1 radial spoke per time-point) have residual aliasing artifacts that propagate into the parametric maps. The artifacts are reduced with the proposed SLLR-MRF resulting in considerable improvements in precision, without changes in accuracy. In vivo results show improved parametric maps for the proposed SLLR-MRF, potentially enabling MRF acquisitions with 1 radial spoke per time-point in approximately 2.6 s (~600 time-point images) for 2 × 2 mm and 9.6 s (1750 time-point images) for 1 × 1 mm in-plane resolution.

CONCLUSION

The proposed SLLR-MRF reconstruction further improves parametric map quality compared with low rank MRF, enabling shorter scan times and/or increased spatial resolution.