Magnetization transfer in magnetic resonance fingerprinting

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.

Reducing femoral flow artefacts in radial magnetic resonance fingerprinting of the prostate using region-optimised virtual coils

High acceleration factors in radial magnetic resonance fingerprinting (MRF) of the prostate lead to strong streak-like artefacts from flow in the femoral blood vessels, possibly concealing important anatomical information. Region-optimised virtual (ROVir) coils is a beamforming-based framework to create virtual coils that maximise signal in a region of interest while minimising signal in a region of interference. In this study, the potential of removing femoral flow streak artefacts in prostate MRF using ROVir coils is demonstrated in silico and in vivo. The ROVir framework was applied to radial MRF k-space data in an automated pipeline designed to maximise prostate signal while minimising signal from the femoral vessels. The method was tested in 15 asymptomatic volunteers at 3 T. The presence of streaks was visually assessed and measurements of whole prostate T1, T2 and signal-to-noise ratio (SNR) with and without streak correction were examined. In addition, a purpose-built simulation framework in which blood flow through the femoral vessels can be turned on and off was used to quantitatively evaluate ROVir's ability to suppress streaks in radial prostate MRF. In vivo it was shown that removing selected ROVir coils visibly reduces streak-like artefacts from the femoral blood flow, without increasing the reconstruction time. On average, 80% of the prostate SNR was retained. A similar reduction of streaks was also observed in silico, while the quantitative accuracy of T1 and T2 mapping was retained. In conclusion, ROVir coils efficiently suppress streaking artefacts from blood flow in radial MRF of the prostate, thereby improving the visual clarity of the images, without significant sacrifices to acquisition time, reconstruction time and accuracy of quantitative values. This is expected to help enable T1 and T2 mapping of prostate cancer in clinically viable times, aiding differentiation between prostate cancer from noncancer and healthy prostate tissue.

High SNR full brain relaxometry at 7T by accelerated MR-STAT

PURPOSE

To demonstrate the feasibility and robustness of the Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) framework for fast, high SNR relaxometry at 7T.

METHODS

To deploy MR-STAT on 7T-systems, we designed optimized flip-angles using the BLAKJac-framework that incorporates the SAR-constraints. Transmit RF-inhomogeneities were mitigated by including a measuredB 1 + $$ {B}_1^{+} $$ -map in the reconstruction. Experiments were performed on a gel-phantom and on five volunteers to explore the robustness of the sequence and its sensitivity toB 1 + $$ {B}_1^{+} $$ inhomogeneities. The SNR-gain at 7T was explored by comparing phantom and in vivo results to MR-STAT at 3T in terms of SNR-efficiency.

RESULTS

The higher SNR at 7T enabled two-fold acceleration with respect to current 2D MR-STAT protocols at lower field strengths. The resulting scan had whole-brain coverage, with 1 x 1 x 3 mm3 resolution (1.5 mm slice-gap) and was acquired within 3 min including theB 1 + $$ {B}_1^{+} $$ -mapping. AfterB 1 + $$ {B}_1^{+} $$ -correction, the estimated T1 and T2 in a phantom showed a mean relative error of, respectively, 1.7% and 4.4%. In vivo, the estimated T1 and T2 in gray and white matter corresponded to the range of values reported in literature with a variation over the subjects of 1.0%-2.1% (WM-GM) for T1 and 4.3%-5.3% (WM-GM) for T2. We measured a higher SNR-efficiency at 7T (R = 2) than at 3T for both T1 and T2 with, respectively, a 4.1 and 2.3 times increase in SNR-efficiency.

CONCLUSION

We presented an accelerated version of MR-STAT tailored to high field (7T) MRI using a low-SAR flip-angle train and showed high quality parameter maps with an increased SNR-efficiency compared to MR-STAT at 3T.

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.

Robust T2 estimation with balanced steady state free precession

PURPOSE

To develop a novel signal representation for balanced steady state free precession (bSSFP) displaying its T2 independence on B1 and on magnetization transfer (MT) effects.

METHODS

A signal model for bSSFP is developed that shows only an explicit dependence (up to a scaling factor) on E2 (and, therefore, T2) and a novel parameter c (with implicit dependence on the flip angle and E1). Moreover, it is shown that MT effects, entering the bSSFP signal via a binary spin bath model, can be captured by a redefinition of T1 and, therefore, leading to modification of E1, resulting in the same signal model. Various sets of phase-cycled bSSFP brain scans (different flip angles, different TR, different RF pulse durations, and different number of phase cycles) were recorded at 3 T. The parameters T2 (E2) and c were estimated using a variable projection (VARPRO) method and Monte-Carlo simulations were performed to assess T2 estimation precision.

RESULTS

Initial experiments confirmed the expected independence of T2 on various protocol settings, such as TR, the flip angle, B1 field inhomogeneity, and the RF pulse duration. Any variation (within the explored range) appears to directly affect the estimation of the parameter c only-in agreement with theory.

CONCLUSION

BSSFP theory predicts an extraordinary feature that all MT and B1-related variational aspects do not enter T2 estimation, making it a potentially robust methodology for T2 quantification, pending validation against existing standards.

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.

DTI-MR fingerprinting for rapid high-resolution whole-brain T1 , T2 , proton density, ADC, and fractional anisotropy mapping

PURPOSE

This study aims to develop a high-efficiency and high-resolution 3D imaging approach for simultaneous mapping of multiple key tissue parameters for routine brain imaging, including T1 , T2 , proton density (PD), ADC, and fractional anisotropy (FA). The proposed method is intended for pushing routine clinical brain imaging from weighted imaging to quantitative imaging and can also be particularly useful for diffusion-relaxometry studies, which typically suffer from lengthy acquisition time.

METHODS

To address challenges associated with diffusion weighting, such as shot-to-shot phase variation and low SNR, we integrated several innovative data acquisition and reconstruction techniques. Specifically, we used M1-compensated diffusion gradients, cardiac gating, and navigators to mitigate phase variations caused by cardiac motion. We also introduced a data-driven pre-pulse gradient to cancel out eddy currents induced by diffusion gradients. Additionally, to enhance image quality within a limited acquisition time, we proposed a data-sharing joint reconstruction approach coupled with a corresponding sequence design.

RESULTS

The phantom and in vivo studies indicated that the T1 and T2 values measured by the proposed method are consistent with a conventional MR fingerprinting sequence and the diffusion results (including diffusivity, ADC, and FA) are consistent with the spin-echo EPI DWI sequence.

CONCLUSION

The proposed method can achieve whole-brain T1 , T2 , diffusivity, ADC, and FA maps at 1-mm isotropic resolution within 10 min, providing a powerful tool for investigating the microstructural properties of brain tissue, with potential applications in clinical and research settings.

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.

Dynamic and rapid deep synthesis of chemical exchange saturation transfer and semisolid magnetization transfer MRI signals

Model-driven analysis of biophysical phenomena is gaining increased attention and utility for medical imaging applications. In magnetic resonance imaging (MRI), the availability of well-established models for describing the relations between the nuclear magnetization, tissue properties, and the externally applied magnetic fields has enabled the prediction of image contrast and served as a powerful tool for designing the imaging protocols that are now routinely used in the clinic. Recently, various advanced imaging techniques have relied on these models for image reconstruction, quantitative tissue parameter extraction, and automatic optimization of acquisition protocols. In molecular MRI, however, the increased complexity of the imaging scenario, where the signals from various chemical compounds and multiple proton pools must be accounted for, results in exceedingly long model simulation times, severely hindering the progress of this approach and its dissemination for various clinical applications. Here, we show that a deep-learning-based system can capture the nonlinear relations embedded in the molecular MRI Bloch-McConnell model, enabling a rapid and accurate generation of biologically realistic synthetic data. The applicability of this simulated data for in-silico, in-vitro, and in-vivo imaging applications is then demonstrated for chemical exchange saturation transfer (CEST) and semisolid macromolecule magnetization transfer (MT) analysis and quantification. The proposed approach yielded 63-99% acceleration in data synthesis time while retaining excellent agreement with the ground truth (Pearson's r > 0.99, p < 0.0001, normalized root mean square error < 3%).

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.

Deep-Learning-Based Contrast Synthesis From MRF Parameter Maps in the Knee Joint

BACKGROUND

Magnetic resonance fingerprinting (MRF) is a method to speed up acquisition of quantitative MRI data. However, MRF does not usually produce contrast-weighted images that are required by radiologists, limiting reachable total scan time improvement. Contrast synthesis from MRF could significantly decrease the imaging time.

PURPOSE

To improve clinical utility of MRF by synthesizing contrast-weighted MR images from the quantitative data provided by MRF, using U-nets that were trained for the synthesis task utilizing L1- and perceptual loss functions, and their combinations.

STUDY TYPE

Retrospective.

POPULATION

Knee joint MRI data from 184 subjects from Northern Finland 1986 Birth Cohort (ages 33-35, gender distribution not available).

FIELD STRENGTH AND SEQUENCE

A 3 T, multislice-MRF, proton density (PD)-weighted 3D-SPACE (sampling perfection with application optimized contrasts using different flip angle evolution), fat-saturated T2-weighted 3D-space, water-excited double echo steady state (DESS).

ASSESSMENT

Data were divided into training, validation, test, and radiologist's assessment sets in the following way: 136 subjects to training, 3 for validation, 3 for testing, and 42 for radiologist's assessment. The synthetic and target images were evaluated using 5-point Likert scale by two musculoskeletal radiologists blinded and with quantitative error metrics.

STATISTICAL TESTS

Friedman's test accompanied with post hoc Wilcoxon signed-rank test and intraclass correlation coefficient. The statistical cutoff P <0.05 adjusted by Bonferroni correction as necessary was utilized.

RESULTS

The networks trained in the study could synthesize conventional images with high image quality (Likert scores 3-4 on a 5-point scale). Qualitatively, the best synthetic images were produced with combination of L1- and perceptual loss functions and perceptual loss alone, while L1-loss alone led to significantly poorer image quality (Likert scores below 3). The interreader and intrareader agreement were high (0.80 and 0.92, respectively) and significant. However, quantitative image quality metrics indicated best performance for the pure L1-loss.

DATA CONCLUSION

Synthesizing high-quality contrast-weighted images from MRF data using deep learning is feasible. However, more studies are needed to validate the diagnostic accuracy of these synthetic images.

EVIDENCE LEVEL

4.

TECHNICAL EFFICACY

Stage 1.

Assessment of MRF for simultaneous T1 and T2 quantification and water-fat separation in the liver at 0.55 T

OBJECTIVE

The goal of this work was to assess the feasibility of performing MRF in the liver on a 0.55 T scanner and to examine the feasibility of water-fat separation using rosette MRF at 0.55 T.

MATERIALS AND METHODS

Spiral and rosette MRF sequences were implemented on a commercial 0.55 T scanner. The accuracy of both sequences in T1 and T2 quantification was validated in the ISMRM/NIST system phantom. The efficacy of rosette MRF in water-fat separation was evaluated in simulations and water/oil phantoms. Both spiral and rosette MRF were performed in the liver of healthy subjects.

RESULTS

In the ISMRM/NIST phantom, both spiral and rosette MRF achieved good agreement with reference values in T1 and T2 measurements. In addition, rosette MRF enables water-fat separation and can generate water- and fat- specific T1 maps, T2 maps, and proton density images from the same dataset for a spatial resolution of 1.56 × 1.56 × 5mm3 within the acquisition time of 15 s.

CONCLUSION

It is feasible to measure T1 and T2 simultaneously in the liver using MRF on a 0.55 T system with lower performance gradients compared to state-of-the-art 1.5 T and 3 T systems within an acquisition time of 15 s. In addition, rosette MRF enables water-fat separation along with T1 and T2 quantification with no time penalty.

Age-Dependent Changes in Knee Cartilage T1 , T2 , and T1p Simultaneously Measured Using MRI Fingerprinting

BACKGROUND

Magnetic resonance fingerprinting (MRF) techniques have been recently described for simultaneous multiparameter cartilage mapping of the knee although investigation of their ability to detect early cartilage degeneration remains limited.

PURPOSE

To investigate age-dependent changes in knee cartilage T1 , T2 , and T1p relaxation times measured using a three-dimensional (3D) MRF sequence in healthy volunteers.

STUDY TYPE

Prospective.

SUBJECTS

The study group consisted of 24 healthy asymptomatic human volunteers (15 males with mean age 34.9 ± 14.4 years and 9 females with mean age 44.5 ± 13.1 years).

FIELD STRENGTH/SEQUENCE

A 3.0 T gradient-echo-based 3D-MRF sequence was used to simultaneously create proton density-weighted images and T1 , T2 , and T1p maps of knee cartilage.

ASSESSMENT

Mean global cartilage and regional cartilage (lateral femur, lateral tibia, medial femur, medial tibia, and patella) T1 , T2 , and T1ρ relaxation times of the knee were measured.

STATISTICAL TESTS

Kruskal-Wallis tests were used to compared cartilage T1 , T2 , and T1ρ relaxation times between different age groups, while Spearman correlation coefficients was used to determine the association between age and cartilage T1 , T2 , and T1ρ relaxation times. The value of P < 0.05 was considered statistically significant.

RESULTS

Higher age groups showed higher global and regional cartilage T1 , T2 , and T1ρ . There was a significant difference between age groups in global cartilage T2 and T1ρ but no significant difference (P = 0.13) in global cartilage T1. Significant difference was also present between age groups in cartilage T2 and T1ρ for medial femur cartilage and medial tibia cartilage. There were significant moderate correlations between age and T2 and T1ρ for global cartilage (R2  = 0.63-0.64), medial femur cartilage (R2  = 0.50-0.56), and medial tibia cartilage (R2  = 0.54-0.66).

CONCLUSION

Cartilage T2 and T1p relaxation times simultaneously measured using a 3D-MRF sequence in healthy volunteers showed age-dependent changes in knee cartilage, primarily within the medial compartment.

Gestational age-related changes in relaxation times of neonatal brain by quantitative synthetic magnetic resonance imaging

OBJECTIVE

This study aimed to explore the correlation between T1 and T2 relaxation times of synthetic MRI (SyMRI) and gestational age (GA) in each hemisphere of preterm and term newborns at the initial 28 days of birth.

METHODS

Seventy preterm and full-term infants were prospectively included in this study. All subjects completed 3.0 T routine MRI and SyMRI (MAGiC) one-stop scanning within 28 days of birth (aged 34-42 W at examination). The SyMRI postprocessing software (v8.0.4) was used to measure the T1 and T2 relaxation values of each brain region. The linear regression equations of quantitative relaxation values with GA were established to compare the variation speed in each brain region.

RESULTS

A significant linear and negative correlation was found between relaxation times and GA in the neonate cerebral cortex and subcortical gray and white matter regions (All p<.05). The relaxation time of the left centrum semiovale decreased with maximum variance with increasing GA among all white matter regions (T1: b = -51.45, β = -0.65, p < .0001; T2: b = -8.77, β = -0.71, p < .0001), whereas the right posterior limb of internal capsule showed minimal variance (T1: b = -27.94, β = -0.60, p < .0001; T2: b = -3.25, β = -0.68, p < .0001). Among all gray matter regions, the right globus pallidus and thalamus indicated the most significant decreasing degree of T1 and T2 relaxation values with GA (right globus pallidus T1: b = -33.14, β = -0.64, p < .0001; right thalamus T2: b = -3.94, β = -0.81, p < .0001), and the right and left occipital lobes indicated the least significant decreasing degree of T1 and T2 relaxation values with GA, respectively (right occipital lobes T1: b = -11.18, β = -0.26, p = .028; left occipital lobes T2: b = -1.22, β = -0.27, p = .024).

CONCLUSIONS

SyMRI could quantitatively evaluate the linear changes of T1 and T2 relaxation values with GA in brain gray and white matter of preterm and term neonates.

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.

Nonaqueous magnetization following adiabatic and selective pulses in brain: T1 and cross-relaxation dynamics

Inversion pulses are commonly employed in MRI for T 1 $$ {T}_1 $$ -weighted contrast and relaxation measurements. In the brain, it is often assumed that adiabatic pulses saturate the nonaqueous magnetization. We investigated this assumption using solid-state NMR to monitor the nonaqueous signal directly following adiabatic inversion and compared this with signals following hard and soft inversion pulses. The effects of the different preparations on relaxation dynamics were explored. Inversion recovery experiments were performed on ex vivo bovine and porcine brains using 360-MHz (8.4 T) and 200-MHz (4.7 T) NMR spectrometers, respectively, using broadband rectangular, adiabatic, and sinc inversion pulses as well as a long rectangular saturation pulse. Analogous human brain MRI experiments were performed at 3 T using single-slice echo-planar imaging. Relaxation data were fitted by mono- and biexponential decay models. Further fitting analysis was performed using only two inversion delay times. Adiabatic and sinc inversion left much of the nonaqueous magnetization along B 0 $$ {B}_0 $$ and resulted in biexponential relaxation. Saturation of both aqueous and nonaqueous magnetization components led to effectively monoexponential T 1 $$ {T}_1 $$ relaxation. Typical adiabatic inversion pulses do not, as has been widely assumed, saturate the nonaqueous proton magnetization in white matter. Unequal magnetization states in aqueous and nonaqueous ¹ H reservoirs prepared by soft and adiabatic pulses result in biexponential T 1 $$ {T}_1 $$ relaxation. Both pools must be prepared in the same magnetization state (e.g., saturated or inverted) in order to observe consistent monoexponential relaxation.

Normative quantitative relaxation atlases for characterization of cortical regions using magnetic resonance fingerprinting

Quantitative magnetic resonance (MR) has been used to study cyto- and myelo-architecture of the human brain non-invasively. However, analyzing brain cortex using high-resolution quantitative MR acquisition can be challenging to perform using 3T clinical scanners. MR fingerprinting (MRF) is a highly efficient and clinically feasible quantitative MR technique that simultaneously provides T1 and T2 relaxation maps. Using 3D MRF from 40 healthy subjects (mean age = 25.6 ± 4.3 years) scanned on 3T magnetic resonance imaging, we generated whole-brain gyral-based normative MR relaxation atlases and investigated cortical-region-based T1 and T2 variations. Gender and age dependency of T1 and T2 variations were additionally analyzed. The coefficient of variation of T1 and T2 for each cortical-region was 3.5% and 7.3%, respectively, supporting low variability of MRF measurements across subjects. Significant differences in T1 and T2 were identified among 34 brain regions (P < 0.001), lower in the precentral, postcentral, paracentral lobule, transverse temporal, lateral occipital, and cingulate areas, which contain sensorimotor, auditory, visual, and limbic functions. Significant correlations were identified between age and T1 and T2 values. This study established whole-brain MRF T1 and T2 atlases of healthy subjects using a clinical 3T scanner, which can provide a quantitative and region-specific baseline for future brain studies and pathology detection.

Quantitative MRI in cardiometabolic disease: From conventional cardiac and liver tissue mapping techniques to multi-parametric approaches

Cardiometabolic disease refers to the spectrum of chronic conditions that include diabetes, hypertension, atheromatosis, non-alcoholic fatty liver disease, and their long-term impact on cardiovascular health. Histological studies have confirmed several modifications at the tissue level in cardiometabolic disease. Recently, quantitative MR methods have enabled non-invasive myocardial and liver tissue characterization. MR relaxation mapping techniques such as T₁, T_1ρ, T₂ and T₂* provide a pixel-by-pixel representation of the corresponding tissue specific relaxation times, which have been shown to correlate with fibrosis, altered tissue perfusion, oedema and iron levels. Proton density fat fraction mapping approaches allow measurement of lipid tissue in the organ of interest. Several studies have demonstrated their utility as early diagnostic biomarkers and their potential to bear prognostic implications. Conventionally, the quantification of these parameters by MRI relies on the acquisition of sequential scans, encoding and mapping only one parameter per scan. However, this methodology is time inefficient and suffers from the confounding effects of the relaxation parameters in each single map, limiting wider clinical and research applications. To address these limitations, several novel approaches have been proposed that encode multiple tissue parameters simultaneously, providing co-registered multiparametric information of the tissues of interest. This review aims to describe the multi-faceted myocardial and hepatic tissue alterations in cardiometabolic disease and to motivate the application of relaxometry and proton-density cardiac and liver tissue mapping techniques. Current approaches in myocardial and liver tissue characterization as well as latest technical developments in multiparametric quantitative MRI are included. Limitations and challenges of these novel approaches, and recommendations to facilitate clinical validation are also discussed.

MR Fingerprinting for Liver Tissue Characterization: A Histopathologic Correlation Study

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

Acceleration Strategies for MR-STAT: Achieving High-Resolution Reconstructions on a Desktop PC Within 3 Minutes

MR-STAT is an emerging quantitative magnetic resonance imaging technique which aims at obtaining multi-parametric tissue parameter maps from single short scans. It describes the relationship between the spatial-domain tissue parameters and the time-domain measured signal by using a comprehensive, volumetric forward model. The MR-STAT reconstruction solves a large-scale nonlinear problem, thus is very computationally challenging. In previous work, MR-STAT reconstruction using Cartesian readout data was accelerated by approximating the Hessian matrix with sparse, banded blocks, and can be done on high performance CPU clusters with tens of minutes. In the current work, we propose an accelerated Cartesian MR-STAT algorithm incorporating two different strategies: firstly, a neural network is trained as a fast surrogate to learn the magnetization signal not only in the full time-domain but also in the compressed low-rank domain; secondly, based on the surrogate model, the Cartesian MR-STAT problem is re-formulated and split into smaller sub-problems by the alternating direction method of multipliers. The proposed method substantially reduces the computational requirements for runtime and memory. Simulated and in-vivo balanced MR-STAT experiments show similar reconstruction results using the proposed algorithm compared to the previous sparse Hessian method, and the reconstruction times are at least 40 times shorter. Incorporating sensitivity encoding and regularization terms is straightforward, and allows for better image quality with a negligible increase in reconstruction time. The proposed algorithm could reconstruct both balanced and gradient-spoiled in-vivo data within 3 minutes on a desktop PC, and could thereby facilitate the translation of MR-STAT in clinical settings.

Repeatability of simultaneous 3D 1H MRF/23Na MRI in brain at 7 T

Proton MRI can provide detailed morphological images, but it reveals little information about cell homeostasis. On the other hand, sodium MRI can provide metabolic information but cannot resolve fine structures. The complementary nature of proton and sodium MRI raises the prospect of their combined use in a single experiment. In this work, we assessed the repeatability of normalized proton density (PD), T₁, T₂, and normalized sodium density-weighted quantification measured with simultaneous 3D ¹H MRF/²³Na MRI in the brain at 7 T, from ten healthy volunteers who were scanned three times each. The coefficients of variation (CV) and the intra-class correlation (ICC) were calculated for the mean and standard deviation (SD) of these 4 parameters in grey matter, white matter, and cerebrospinal fluid. As result, the CVs were lower than 3.3% for the mean values and lower than 6.9% for the SD values. The ICCs were higher than 0.61 in all 24 measurements. We conclude that the measurements of normalized PD, T₁, T₂, and normalized sodium density-weighted from simultaneous 3D ¹H MRF/²³Na MRI in the brain at 7 T showed high repeatability. We estimate that changes > 6.6% (> 2 CVs) in mean values of both ¹H and ²³Na metrics could be detectable with this method.

B0 Correction for 3T Amide Proton Transfer (APT) MRI Using a Simplified Two-Pool Lorentzian Model of Symmetric Water and Asymmetric Solutes

Amide proton transfer (APT)-weighted MRI is a promising molecular imaging technique that has been employed in clinic for detection and grading of brain tumors. MTRasym, the quantification method of APT, is easily influenced by B0 inhomogeneity and causes artifacts. Current model-free interpolation methods have enabled moderate B0 correction for middle offsets, but have performed poorly at limbic offsets. To address this shortcoming, we proposed a practical B0 correction approach that is suitable under time-limited sparse acquisition scenarios and for B1 ≥ 1 μT under 3T. In this study, this approach employed a simplified Lorentzian model containing only two pools of symmetric water and asymmetric solutes, to describe the Z-spectral shape with wide and ‘invisible’ CEST peaks. The B0 correction was then performed on the basis of the fitted two-pool Lorentzian lines, instead of using conventional model-free interpolation. The approach was firstly evaluated on densely sampled Z-spectra data by using the spline interpolation of all acquired 16 offsets as the gold standard. When only six offsets were available for B0 correction, our method outperformed conventional methods. In particular, the errors at limbic offsets were significantly reduced (n = 8, p < 0.01). Secondly, our method was assessed on the six-offset APT data of nine brain tumor patients. Our MTRasym (3.5 ppm), using the two-pool model, displayed a similar contrast to the vendor-provided B0-orrected MTRasym (3.5 ppm). While the vendor failed in correcting B0 at 4.3 and 2.7 ppm for a large portion of voxels, our method enabled well differentiation of B0 artifacts from tumors. In conclusion, the proposed approach could alleviate analysis errors caused by B0 inhomogeneity, which is useful for facilitating the comprehensive metabolic analysis of brain tumors.

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Circumventing the curse of dimensionality in magnetic resonance fingerprinting through a deep learning approach

Magnetic resonance fingerprinting (MRF) is a rapidly developing approach for fast quantitative MRI. A typical drawback of dictionary-based MRF is an explosion of the dictionary size as a function of the number of reconstructed parameters, according to the "curse of dimensionality", which determines an explosion of resource requirements. Neural networks (NNs) have been proposed as a feasible alternative, but this approach is still in its infancy. In this work, we design a deep learning approach to MRF using a fully connected network (FCN). In the first part we investigate, by means of simulations, how the NN performance scales with the number of parameters to be retrieved in comparison with the standard dictionary approach. Four MRF sequences were considered: IR-FISP, bSSFP, IR-FISP-B₁ , and IR-bSSFP-B₁ , the latter two designed to be more specific for B1+ parameter encoding. Estimation accuracy, memory usage, and computational time required to perform the estimation task were considered to compare the scalability capabilities of the dictionary-based and the NN approaches. In the second part we study optimal training procedures by including different data augmentation and preprocessing strategies during training to achieve better accuracy and robustness to noise and undersampling artifacts. The study is conducted using the IR-FISP MRF sequence exploiting both simulations and in vivo acquisitions. Results demonstrate that the NN approach outperforms the dictionary-based approach in terms of scalability capabilities. Results also allow us to heuristically determine the optimal training strategy to make an FCN able to predict T₁ ,  T₂ , and M₀ maps that are in good agreement with those obtained with the original dictionary approach. k-SVD denoising is proposed and found to be critical as a preprocessing step to handle undersampled data.

Using magnetic resonance fingerprinting to characterize periventricular nodular heterotopias in pharmacoresistant epilepsy

OBJECTIVE

We aimed to use a novel magnetic resonance fingerprinting (MRF) technique to examine in vivo tissue property characteristics of periventricular nodular heterotopia (PVNH). These characteristics were further correlated with stereotactic-electroencephalographic (SEEG) ictal onset findings.

METHODS

We included five patients with PVNH who had SEEG-guided surgery and at least 1 year of seizure freedom or substantial seizure reduction. High-resolution MRF scans were acquired at 3 T, generating three-dimensional quantitative T₁ and T₂  maps. We assessed the differences between T₁ and T₂  values from the voxels in the nodules located in the SEEG-defined seizure onset zone (SOZ) and non-SOZ, on -individual and group levels. Receiver operating characteristic analyses were performed to obtain the optimal classification performance. Quantification of SEEG ictal onset signals from the nodules was performed by calculating power spectrum density (PSD). The association between PSD and T₁ /T₂  values was further assessed at different frequency bands.

RESULTS

Individual-level analysis showed T₁ was significantly higher in SOZ voxels than non-SOZ voxels (p < .05), with an average 73% classification accuracy. Group-level analysis also showed higher T₁ was significantly associated with SOZ voxels (p < .001). At the optimal cutoff (normalized T₁ of 1.1), a 76% accuracy for classifying SOZ nodules from non-SOZ nodules was achieved. T₁  values were significantly associated with ictal onset PSD at the ultraslow, θ, β, γ, and ripple bands (p < .05). T₂  values were significantly associated with PSD only at the ultraslow band (p < .05).

SIGNIFICANCE

Quantitative MRF measures, especially T₁ , can provide additional noninvasive information to separate nodules in SOZ and non-SOZ. The T₁ and T₂ tissue property changes carry electrophysiological underpinnings relevant to the epilepsy, as shown by their significant positive associations with power changes during the SEEG seizure onset. The use of MRF as a supplementary noninvasive tool may improve presurgical evaluation for patients with PVNH and pharmacoresistant epilepsy.

Characterization of cardiac amyloidosis using cardiac magnetic resonance fingerprinting

BACKGROUND

Cardiac amyloidosis (CA) is an infiltrative cardiomyopathy with poor prognosis absent appropriate treatment. Elevated native myocardial T₁ and T₂ have been reported for CA, and tissue characterization by cardiac MRI may expedite diagnosis and treatment. Cardiac Magnetic Resonance Fingerprinting (cMRF) has the potential to enable tissue characterization for CA through rapid, simultaneous T₁ and T₂ mapping. Furthermore, cMRF signal timecourses may provide additional information beyond myocardial T₁ and T₂.

METHODS

Nine CA patients and five controls were scanned at 3 T using a prospectively gated cMRF acquisition. Two cMRF-based analysis approaches were examined: (1) relaxometric-based linear discriminant analysis (LDA) using native T₁ and T₂, and (2) signal timecourse-based LDA. The Fisher coefficient was used to compare the separability of patient and control groups from both approaches. Leave-two-out cross-validation was employed to evaluate the classification error rates of both approaches.

RESULTS

Elevated myocardial T₁ and T₂ was observed in patients vs controls (T₁: 1395 ± 121 vs 1240 ± 36.4 ms, p < 0.05; T₂: 36.8 ± 3.3 vs 31.8 ± 2.6 ms, p < 0.05). LDA scores were elevated in patients for relaxometric-based LDA (0.56 ± 0.28 vs 0.18 ± 0.13, p < 0.05) and timecourse-based LDA (0.97 ± 0.02 vs 0.02 ± 0.02, p < 0.05). The Fisher coefficient was greater for timecourse-based LDA (60.8) vs relaxometric-based LDA (1.6). Classification error rates were lower for timecourse-based LDA vs relaxometric-based LDA (12.6 ± 24.3 vs 22.5 ± 30.1%, p < 0.05).

CONCLUSIONS

These findings suggest that cMRF may be a valuable technique for the detection and characterization of CA. Analysis of cMRF signal timecourse data may improve tissue characterization as compared to analysis of native T₁ and T₂ alone.

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

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

An MR fingerprinting approach for quantitative inhomogeneous magnetization transfer imaging

PURPOSE

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

THEORY AND METHODS

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

RESULTS

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

CONCLUSION

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

Feasibility of Magnetic Resonance Fingerprinting on Aging MRI Hardware

The purpose of this work is to evaluate the feasibility of performing magnetic resonance fingerprinting (MRF) on older and lower-performance MRI hardware as a means to bring advanced imaging to the aging MRI install base. Phantom and in vivo experiments were performed on a 1.5T Siemens Aera (installed 2015) and 1.5T Siemens Symphony (installed 2002). A 2D spiral MRF sequence for simultaneous T₁/T₂/M₀ mapping was implemented on both scanners with different gradient trajectories to accommodate system specifications. In phantom, for T₁/T₂ values in a physiologically relevant range (T₁: 195-1539 ms; T₂: 20-267 ms), scanners had strong correlation (R² > 0.999) with average absolute percent difference of 8.1% and 10.1%, respectively. Comparison of the two trajectories on the newer scanner showed differences of 2.6% (T₁) and 10.9% (T₂), suggesting a partial explanation of the observed inter-scanner bias. Inter-scanner agreement was better when the same trajectory was used, with differences of 6.0% (T₁) and 4.0% (T₂). Intra-scanner coefficient of variation (CV) of T₁ and T₂ estimates in phantom were <2.0% and in vivo were ≤3.5%. In vivo inter-scanner white matter CV was 4.8% (T₁) and 5.1% (T₂). White matter measurements on the aging scanner after two months were consistent, with differences of 1.9% (T₁) and 3.9% (T₂). In conclusion, MRF is feasible on an aging MRI scanner and required only changes to the gradient trajectory.

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

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

Magnetic resonance fingerprinting residual signals can disassociate human grey matter regions

The importance of accurate structural discrimination of the human grey matter regions has motivated the development of observer-independent reproducible methods that account for inter-individual architectonic variations. We introduce a non-invasive statistical residual analysis framework, employing unique tissue-specific magnetic resonance fingerprinting (MRF) signals after adjusting for the effect of T₁ and T₂^* MR relaxometry parameters (here termed MRF residuals). A 7 T Siemens MR scanner was used to acquire MRF signals, quantitative transmit magnetic field (B₁⁺) maps and T₁-weighted anatomical images of eleven cortical areas (5L, 5M, 5Ci, 7A, 7P, 7PC, hIP3, BA2, BA4a, BA4p and BA6) from six female participants. MRF residual signal for each voxel was calculated as the difference between the actual and best matching MRF signal evolutions from a precomputed MRF dictionary covering a range of T₁, T₂^* and B₁⁺ values. To compare MRF residuals between regions of interest, normalised autocorrelation was used as a shape-based statistical signal characterisation method and the Euclidean distance between autocorrelation profiles of residuals was used to measure the interareal dissimilarity. In the eleven cortical areas in both cerebral hemispheres of six participants, the proposed MRF residual analysis consistently showed interareal dissimilarity profiles that concorded with histological studies, indicating that MRF residuals potentially contain tissue microstructural information. MRF residual signals provide additional area-specific information that is complementary to the MR relaxometry-based (T₁, T₂^*) information used previously for distinguishing microstructural differences between human cerebral cortex regions in vivo. The proposed approach led to more accurate identification of structural variations across cortical areas of interest.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

A deep learning approach for magnetic resonance fingerprinting: Scaling capabilities and good training practices investigated by simulations

MR fingerprinting (MRF) is an innovative approach to quantitative MRI. A typical disadvantage of dictionary-based MRF is the explosive growth of the dictionary as a function of the number of reconstructed parameters, an instance of the curse of dimensionality, which determines an explosion of resource requirements. In this work, we describe a deep learning approach for MRF parameter map reconstruction using a fully connected architecture. Employing simulations, we have investigated how the performance of the Neural Networks (NN) approach scales with the number of parameters to be retrieved, compared to the standard dictionary approach. We have also studied optimal training procedures by comparing different strategies for noise addition and parameter space sampling, to achieve better accuracy and robustness to noise. Four MRF sequences were considered: IR-FISP, bSSFP, IR-FISP-B₁, and IR-bSSFP-B₁. A comparison between NN and the dictionary approaches in reconstructing parameter maps as a function of the number of parameters to be retrieved was performed using a numerical brain phantom. Results demonstrated that training with random sampling and different levels of noise variance yielded the best performance. NN performance was at least as good as the dictionary-based approach in reconstructing parameter maps using Gaussian noise as a source of artifacts: the difference in performance increased with the number of estimated parameters because the dictionary method suffers from the coarse resolution of the parameter space sampling. The NN proved to be more efficient in memory usage and computational burden, and has great potential for solving large-scale MRF problems.

Magnetization-prepared GRASP MRI for rapid 3D T1 mapping and fat/water-separated T1 mapping

PURPOSE

This study aimed to (i) develop Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI using a stack-of-stars trajectory for rapid free-breathing T1 mapping and (ii) extend MP-GRASP to multi-echo acquisition (MP-Dixon-GRASP) for fat/water-separated (water-specific) T1 mapping.

METHODS

An adiabatic non-selective 180° inversion-recovery pulse was added to a gradient-echo-based golden-angle stack-of-stars sequence for magnetization-prepared 3D single-echo or 3D multi-echo acquisition. In combination with subspace-based GRASP-Pro reconstruction, the sequence allows for standard T1 mapping (MP-GRASP) or fat/water-separated T1 mapping (MP-Dixon-GRASP), respectively. The accuracy of T1 mapping using MP-GRASP was evaluated in a phantom and volunteers (brain and liver) against clinically accepted reference methods. The repeatability of T1 estimation was also assessed in the phantom and volunteers. The performance of MP-Dixon-GRASP for water-specific T1 mapping was evaluated in a fat/water phantom and volunteers (brain and liver).

RESULTS

ROI-based mean T1 values are correlated between the references and MP-GRASP in the phantom (R2 = 1.0), brain (R2 = 0.96), and liver (R2 = 0.73). MP-GRASP achieved good repeatability of T1 estimation in the phantom (R2 = 1.0), brain (R2 = 0.99), and liver (R2 = 0.82). Water-specific T1 is different from in-phase and out-of-phase composite T1 (composite T1 when fat and water signal are mixed in phase or out of phase) both in the phantom and volunteers.

CONCLUSION

This work demonstrated the initial performance of MP-GRASP and MP-Dixon-GRASP MRI for rapid 3D T1 mapping and 3D fat/water-separated T1 mapping in the brain (without motion) and in the liver (during free breathing). With fat/water-separated T1 estimation, MP-Dixon-GRASP could be potentially useful for imaging patients with fatty-liver diseases.

Fast and accurate modeling of transient-state, gradient-spoiled sequences by recurrent neural networks

Fast and accurate modeling of MR signal responses are typically required for various quantitative MRI applications, such as MR fingerprinting. This work uses a new extended phase graph (EPG)-Bloch model for accurate simulation of transient-state, gradient-spoiled MR sequences, and proposes a recurrent neural network (RNN) as a fast surrogate of the EPG-Bloch model for computing large-scale MR signals and derivatives. The computational efficiency of the RNN model is demonstrated by comparisons with other existing models, showing one to three orders of acceleration compared with the latest GPU-accelerated, open-source EPG package. By using numerical and in vivo brain data, two used cases, namely, MRF dictionary generation and optimal experimental design, are also provided. Results show that the RNN surrogate model can be efficiently used for computing large-scale dictionaries of transient-state signals and derivatives within tens of seconds, resulting in several orders of magnitude acceleration with respect to state-of-the-art implementations. The practical application of transient-state quantitative techniques can therefore be substantially facilitated.

Feasibility of MR fingerprinting using a high-performance 0.55 T MRI system

BACKGROUND

MR fingerprinting (MRF) is a versatile method for rapid multi-parametric quantification. The application of MRF for lower MRI field could enable multi-contrast imaging and improve exam efficiency on these systems. The purpose of this work is to demonstrate the feasibility of 3D whole-brain T1 and T2 mapping using MR fingerprinting on a contemporary 0.55 T MRI system.

MATERIALS AND METHODS

A 3D whole brain stack-of-spirals FISP MRF sequence was implemented for 0.55 T. Quantification was validated using the NIST/ISMRM Quantitative MRI phantom, and T1 and T2 values of white matter, gray matter, and cerebrospinal fluid were measured in 19 healthy subjects. To assess MRF performance in the lower SNR regime of 0.55 T, measurement precision was calculated from 100 simulated pseudo-replicas of in vivo data and within-session measurement repeatability was evaluated.

RESULTS

T1 and T2 values calculated by MRF were strongly correlated to standard measurements in the ISMRM/NIST MRI system phantom (R² > 0.99), with a small constant bias of approximately 5 ms in T2 values. 3D stack-of-spirals MRF was successfully applied for whole brain quantitative T1 and T2 at 0.55 T, with spatial resolution of 1.2 mm × 1.2 mm × 5 mm, and acquisition time of 8.5 min. Moreover, the T1 and T2 quantifications had precision <5%, despite the lower SNR of 0.55 T.

CONCLUSION

A 3D whole-brain stack-of-spirals FISP MRF sequence is feasible for T1 and T2 mapping at 0.55 T.

Simultaneous T1 , T2 , and T1ρ relaxation mapping of the lower leg muscle with MR fingerprinting

PURPOSE

To develop a novel MR-fingerprinting (MRF) pulse sequence that is insensitive to B 1 + and B₀ imperfections for simultaneous T₁ , T₂ , and T_1ρ relaxation mapping.

METHODS

We implemented a totally balanced spin-lock (TB-SL) module to encode T_1ρ relaxation into an existing MRF framework that encoded T₁ and T₂ . The spin-lock module used two 180° pulses with compensatory phases to reduce T_1ρ sensitivity to B₁ and B₀ inhomogeneities. We compared T_1ρ measured using TB-SL MRF in Bloch simulations, model agar phantoms, and in vivo experiments to those with a self-compensated spin-lock preparation module (SC-SL). The TB-SL MRF repeatability was evaluated in maps acquired in the lower leg skeletal muscle of 12 diabetic peripheral neuropathy patients, scanned two times each during visits separated by about 30 days.

RESULTS

The phantom relaxation times measured with TB-SL and SC-SL MRF were in good agreement with reference values in regions with low B₁ inhomogeneities. Compared with SC-SL, TB-SL MRF showed in experiments greater robustness against severe B₁ inhomogeneities and in Bloch simulations greater robustness against B₁ and B₀ . We measured with TB-SL MRF an average T₁ = 950.1 ± 28.7 ms, T₂ = 26.0 ± 1.2 ms, and T_1ρ = 31.7 ± 3.2 ms in skeletal muscle across patients. Bland-Altman analysis demonstrated low bias between TB-SL and SC-SL MRF and between TB-SL MRF maps acquired in two visits. The coefficient of variation was less than 3% for all measurements.

CONCLUSION

The proposed TB-SL MRF sequence is fast and insensitive to B 1 + and B₀ imperfections. It can simultaneously map T₁ , T₂ , T_1ρ , and B 1 + in a single scan and can potentially be used to study muscle composition.

Cardiac magnetic resonance fingerprinting: Trends in technical development and potential clinical applications

Quantitative cardiac magnetic resonance has emerged in recent years as an approach for evaluating a range of cardiovascular conditions, with T1 and T2 mapping at the forefront of these developments. Cardiac Magnetic Resonance Fingerprinting (cMRF) provides a rapid and robust framework for simultaneous quantification of myocardial T1 and T2 in addition to other tissue properties. Since the advent of cMRF, a number of technical developments and clinical validation studies have been reported. This review provides an overview of cMRF, recent technical developments, healthy subject and patient studies, anticipated technical improvements, and potential clinical applications. Recent technical developments include slice profile and pulse efficiency corrections, improvements in image reconstruction, simultaneous multislice imaging, 3D whole-ventricle imaging, motion-resolved imaging, fat-water separation, and machine learning for rapid dictionary generation. Future technical developments in cMRF, such as B0 and B1 field mapping, acceleration of acquisition and reconstruction, imaging of patients with implanted devices, and quantification of additional tissue properties are also described. Potential clinical applications include characterization of infiltrative, inflammatory, and ischemic cardiomyopathies, tissue characterization in the left atrium and right ventricle, post-cardiac transplantation assessment, reduction of contrast material, pre-procedural planning for electrophysiology interventions, and imaging of patients with implanted devices.

Toward Quantification: Microstructure and Magnetic Resonance Fingerprinting

Quantitative magnetic resonance imaging (MRI) is a long-standing challenge. We advocate that the origin of the problem is the simplification applied in commonly used models of the MRI signal relation to the target parameters of biological tissues. Two research fields are briefly reviewed as ways to respond to the challenge of quantitative MRI, both experiencing an exponential growth right now. Microstructure MRI strives to build physiology-based models from cells to signal and, given the signal, back to the cells again. Magnetic resonance fingerprinting aims at efficient simultaneous determination of multiple signal parameters. The synergy of these yet disjoined approaches promises truly quantitative MRI with specific target-oriented diagnostic tools rather than universal imaging methods.

3D magnetic resonance fingerprinting with quadratic RF phase

PURPOSE

To implement 3D magnetic resonance fingerprinting (MRF) with quadratic RF phase (qRF-MRF) for simultaneous quantification of T₁ , T₂ , ΔB₀ , and T 2 ∗ .

METHODS

3D MRF data with effective undersampling factor of 3 in the slice direction were acquired with quadratic RF phase patterns for T₁ , T₂ , and T 2 ∗ sensitivity. Quadratic RF phase encodes the off-resonance by modulating the on-resonance frequency linearly in time. Transition to 3D brings practical limitations for reconstruction and dictionary matching because of increased data and dictionary sizes. Randomized singular value decomposition (rSVD)-based compression in time and reduction in dictionary size with a quadratic interpolation method are combined to be able to process prohibitively large data sets in feasible reconstruction and matching times.

RESULTS

Accuracy of 3D qRF-MRF maps in various resolutions and orientations are compared to 3D fast imaging with steady-state precession (FISP) for T₁ and T₂ contrast and to 2D qRF-MRF for T 2 ∗ contrast and ΔB₀ . The precision of 3D qRF-MRF was 1.5-2 times higher than routine clinical scans. 3D qRF-MRF ΔB₀ maps were further processed to highlight the susceptibility contrast.

CONCLUSION

Natively co-registered 3D whole brain T₁ , T₂ , T 2 ∗ , ΔB₀ , and QSM maps can be acquired in as short as 5 min with 3D qRF-MRF.

MR fingerprinting for rapid simultaneous T1 , T2 , and T1ρ relaxation mapping of the human articular cartilage at 3T

PURPOSE

To implement a novel technique for simultaneous, quantitative multiparametric mapping of the knee articular cartilage.

METHODS

A novel MR fingerprinting pulse sequence is proposed and implemented for simultaneous measurements of proton density, T1 , T2, and T1ρ relaxation times at 3T. The repeatability and reproducibility of the proposed technique were assessed in model phantoms. Institutional review board-approved MR fingerprinting imaging sequence was performed on healthy volunteers and patients with mild knee osteoarthritis. The Wilcoxon test was used to compare healthy controls and patients. The intra- and intersubject repeatability were assessed with coefficient of variation and the RMS coefficient of variation, respectively RESULTS: The Bland-Altman plots demonstrated an average difference of 4.67 ms, -0.09 ms, and 0.05 ms between 2 scans in the same scanner; and 9.68 ms, 0.29 ms, and -0.72 ms between the scans acquired on 2 different scanners for T1 , T2 , and T1ρ , respectively. The in vivo knee study showed excellent repeatability with RMS coefficient of variation less than 3%, 6%, and 5% for T1 , T2 , and T1ρ , respectively. The Wilcoxon test showed a significant difference between control and mild osteoarthritis patients for T1 (P = .04), T2 (P = .01), and T1ρ (P = .02) relaxation time in medial tibial cartilage, as well as for T2 relaxation time (P = .02) in medial femoral cartilage.

CONCLUSION

The proposed MRF sequence is fast and can simultaneously measure the T1 , T2 , T1ρ , and B 1 + maps in a single scan. It is able to discriminate between mild osteoarthritis patients and healthy volunteers.

Confounding factors in breast magnetic resonance fingerprinting: B 1 + , slice profile, and diffusion effects

PURPOSE

Magnetic resonance fingerprinting (MRF) offers rapid quantitative imaging but may be subject to confounding effects (CE) if these are not included in the model-based reconstruction. This study characterizes the influence of in-plane B 1 + , slice profile and diffusion effects on T₁ and T₂ estimation in the female breast at 1.5T.

METHODS

Simulations were used to predict the influence of each CE on the accuracy of MRF and to investigate the influence of electronic noise and spiral aliasing artefacts. The experimentally observed bias in regions of fibroglandular tissue (FGT) and fatty tissue (FT) was analyzed for undersampled spiral breast MRF data of 6 healthy volunteers by performing MRF reconstruction with and without a CE.

RESULTS

Theoretic analysis predicts T₁ under-/T₂ overestimation if the nominal flip angles are underestimated and inversely, T₁ under-/T₂ overestimation if omitting slice profile correction, and T₁ under-/T₂ underestimation if omitting diffusion in the signal model. Averaged over repeated signal simulations, including spiral aliasing artefacts affected precision more than accuracy. Strong in-plane B 1 + effects occurred in vivo, causing T₂ left-right inhomogeneity between both breasts. Their correction decreased the T₂ difference from 29 to 5 ms in FGT and from 29 to 9 ms in FT. Slice profile correction affected FGT T₂ most strongly, resulting in -22% smaller values. For the employed spoiler gradient strengths, diffusion did not affect the parameter maps, corresponding well with theoretic predictions.

CONCLUSION

Understanding CEs and their relative significance for an MRF sequence is important when defining an MRF signal model for accurate parameter mapping.

Magnetic resonance fingerprinting of the pancreas at 1.5 T and 3.0 T

Magnetic resonance imaging of the pancreas is increasingly used as an important diagnostic modality for characterisation of pancreatic lesions. Pancreatic MRI protocols are mostly qualitative due to time constraints and motion sensitivity. MR Fingerprinting is an innovative acquisition technique that provides qualitative data and quantitative parameter maps from a single free-breathing acquisition with the potential to reduce exam times. This work investigates the feasibility of MRF parameter mapping for pancreatic imaging in the presence of free-breathing exam. Sixteen healthy participants were prospectively imaged using MRF framework. Regions-of-interest were drawn in multiple solid organs including the pancreas and T1 and T2 values determined. MRF T1 and T2 mapping was performed successfully in all participants (acquisition time:2.4-3.6 min). Mean pancreatic T1 values were 37-43% lower than those of the muscle, spleen, and kidney at both 1.5 and 3.0 T. For these organs, the mean pancreatic T2 values were nearly 40% at 1.5 T and < 12% at 3.0 T. The feasibility of MRF at 1.5 T and 3 T was demonstrated in the pancreas. By enabling fast and free-breathing quantitation, MRF has the potential to add value during the clinical characterisation and grading of pathological conditions, such as pancreatitis or cancer.

Probing myelin content of the human brain with MRI: A review

Rapid and efficient transmission of electric signals among neurons of vertebrates is ensured by myelin-insulating sheaths surrounding axons. Human cognition, sensation, and motor functions rely on the integrity of these layers, and demyelinating diseases often entail serious cognitive and physical impairments. Magnetic resonance imaging radically transformed the way these disorders are monitored, offering an irreplaceable tool to noninvasively examine the brain structure. Several advanced techniques based on MRI have been developed to provide myelin-specific contrasts and a quantitative estimation of myelin density in vivo. Here, the vast offer of acquisition strategies developed to date for this task is reviewed. Advantages and pitfalls of the different approaches are compared and discussed.

Normal volumetric and T1 relaxation time values at 1.5 T in segmented pediatric brain MRI using a MP2RAGE acquisition

OBJECTIVES

This study introduced a tailored MP2RAGE-based brain acquisition for a comprehensive assessment of the normal maturing brain.

METHODS

Seventy normal patients (35 girls and 35 boys) from 1 to 16 years of age were recruited within a prospective monocentric study conducted from a single University Hospital. Brain MRI examinations were performed at 1.5 T using a 20-channel head coil and an optimized 3D MP2RAGE sequence with a total acquisition time of 6:36 min. Automated 38 region segmentation was performed using the MorphoBox (template registration, bias field correction, brain extraction, and tissue classification) which underwent a major adaptation of three age-group T1-weighted templates. Volumetry and T1 relaxometry reference ranges were established using a logarithmic model and a modified Gompertz growth respectively.

RESULTS

Detailed automated brain segmentation and T1 mapping were successful in all patients. Using these data, an age-dependent model of normal brain maturation with respect to changes in volume and T1 relaxometry was established. After an initial rapid increase until 24 months of life, the total intracranial volume was found to converge towards 1400 mL during adolescence. The expected volumes of white matter (WM) and cortical gray matter (GM) showed a similar trend with age. After an initial major decrease, T1 relaxation times were observed to decrease progressively in all brain structures. The T1 drop in the first year of life was more pronounced in WM (from 1000-1100 to 650-700 ms) than in GM structures.

CONCLUSION

The 3D MP2RAGE sequence allowed to establish brain volume and T1 relaxation time normative ranges in pediatrics.

KEY POINTS

• The 3D MP2RAGE sequence provided a reliable quantitative assessment of brain volumes and T1 relaxation times during childhood. • An age-dependent model of normal brain maturation was established. • The normative ranges enable an objective comparison to a normal cohort, which can be useful to further understand, describe, and identify neurodevelopmental disorders in children.

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.