Reproducibility and Repeatability of MR Fingerprinting Relaxometry in the Human Brain.

Background Only sparse literature investigates the reproducibility and repeatability of relaxometry methods in MRI. However, statistical data on reproducibility and repeatability of any quantitative method is essential for clinical application. Purpose To evaluate the reproducibility and repeatability of two-dimensional fast imaging with steady-state free precession MR fingerprinting in vivo in human brains. Materials and Methods Two-dimensional section-selective MR fingerprinting based on a steady-state free precession sequence with an external radiofrequency transmit field, or B1+, correction was used to generate T1 and T2 maps. This prospective study was conducted between July 2017 and January 2018 with 10 scanners from a single manufacturer, including different models, at four different sites. T1 and T2 relaxation times and their variation across scanners (reproducibility) as well as across repetitions on a scanner (repeatability) were analyzed. The relative deviations of T1 and T2 to the average (95% confidence interval) were calculated for several brain compartments. Results Ten healthy volunteers (mean age ± standard deviation, 28.5 years ± 6.9; eight men, two women) participated in this study. Reproducibility and repeatability of T1 and T2 measures in the human brain varied across brain compartments (1.8%-20.9%) and were higher in solid tissues than in the cerebrospinal fluid. T1 measures in solid tissue brain compartments were more stable compared with T2 measures. The half-widths of the confidence intervals for relative deviations were 3.4% for mean T1 and 8.0% for mean T2 values across scanners. Intrascanner repeatability half-widths of the confidence intervals for relative deviations were in the range of 2.0%-3.1% for T1 and 3.1%-7.9% for T2. Conclusion This study provides values on reproducibility and repeatability of T1 and T2 relaxometry measured with fast imaging with steady-state free precession MR fingerprinting in brain tissues of healthy volunteers. Reproducibility and repeatability are considerably higher in solid brain compartments than in cerebrospinal fluid and are higher for T1 than for T2. © RSNA, 2019 Online supplemental material is available for this article. See also the editorial by Barkhof and Parker in this issue.

Magnetic resonance field fingerprinting.

PURPOSE: To develop and evaluate the magnetic resonance field fingerprinting method that simultaneously generates T1 , T2 , B0 , and B 1 + maps from a single continuous measurement. METHODS: An encoding pattern was designed to integrate true fast imaging with steady-state precession (TrueFISP), fast imaging with steady-state precession (FISP), and fast low-angle shot (FLASH) sequence segments with varying flip angles, radio frequency (RF) phases, TEs, and gradient moments in a continuous acquisition. A multistep matching process was introduced that includes steps for integrated spiral deblurring and the correction of intravoxel phase dispersion. The method was evaluated in phantoms as well as in vivo studies in brain and lower abdomen. RESULTS: Simultaneous measurement of T1 , T2 , B0 , and B 1 + is achieved with T1 and T2 subsequently being less afflicted by B0 and B 1 + variations. Phantom results demonstrate the stability of generated parameter maps. Higher undersampling factors and spatial resolution can be achieved with the proposed method as compared with solely FISP-based magnetic resonance fingerprinting. High-resolution B0 maps can potentially be further used as diagnostic information. CONCLUSION: The proposed magnetic resonance field fingerprinting method can estimate T1 , T2 , B0 , and B 1 + maps accurately in phantoms, in the brain, and in the lower abdomen.

Physics-Informed Discretization for Reproducible and Robust Radiomic Feature Extraction Using Quantitative MRI.

OBJECTIVE: Given the limited repeatability and reproducibility of radiomic features derived from weighted magnetic resonance imaging (MRI), there may be significant advantages to using radiomics in conjunction with quantitative MRI. This study introduces a novel physics-informed discretization (PID) method for reproducible radiomic feature extraction and evaluates its performance using quantitative MRI sequences including magnetic resonance fingerprinting (MRF) and apparent diffusion coefficient (ADC) mapping. MATERIALS AND METHODS: A multiscanner, scan-rescan dataset comprising whole-brain 3D quantitative (MRF T1, MRF T2, and ADC) and weighted MRI (T1w MPRAGE, T2w SPACE, and T2w FLAIR) from 5 healthy subjects was prospectively acquired. Subjects underwent 2 repeated acquisitions on 3 distinct 3 T scanners each, for a total of 6 scans per subject (30 total scans). First-order statistical (n = 23) and second-order texture (n = 74) radiomic features were extracted from 56 brain tissue regions of interest using the proposed PID method (for quantitative MRI) and conventional fixed bin number (FBN) discretization (for quantitative MRI and weighted MRI). Interscanner radiomic feature reproducibility was measured using the intraclass correlation coefficient (ICC), and the effect of image sequence (eg, MRF T1 vs T1w MPRAGE), as well as image discretization method (ie, PID vs FBN), on radiomic feature reproducibility was assessed using repeated measures analysis of variance. The robustness of PID and FBN discretization to segmentation error was evaluated by simulating segmentation differences in brainstem regions of interest. Radiomic features with ICCs greater than 0.75 following simulated segmentation were determined to be robust to segmentation. RESULTS: First-order features demonstrated higher reproducibility in quantitative MRI than weighted MRI sequences, with 30% (n = 7/23) features being more reproducible in MRF T1 and MRF T2 than weighted MRI. Gray level co-occurrence matrix (GLCM) texture features extracted from MRF T1 and MRF T2 were significantly more reproducible using PID compared with FBN discretization; for all quantitative MRI sequences, PID yielded the highest number of texture features with excellent reproducibility (ICC > 0.9). Comparing texture reproducibility of quantitative and weighted MRI, a greater proportion of MRF T1 (n = 225/370, 61%) and MRF T2 (n = 150/370, 41%) texture features had excellent reproducibility (ICC > 0.9) compared with T1w MPRAGE (n = 148/370, 40%), ADC (n = 115/370, 32%), T2w SPACE (n = 98/370, 27%), and FLAIR (n = 102/370, 28%). Physics-informed discretization was also more robust than FBN discretization to segmentation error, as 46% (n = 103/222, 46%) of texture features extracted from quantitative MRI using PID were robust to simulated 6 mm segmentation shift compared with 19% (n = 42/222, 19%) of weighted MRI texture features extracted using FBN discretization. CONCLUSIONS: The proposed PID method yields radiomic features extracted from quantitative MRI sequences that are more reproducible and robust than radiomic features extracted from weighted MRI using conventional (FBN) discretization approaches. Quantitative MRI sequences also demonstrated greater scan-rescan robustness and first-order feature reproducibility than weighted MRI.

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.

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 (R2 > 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.

A parallel spatial and Bloch manifold regularized iterative reconstruction method for MR Fingerprinting.

Magnetic Resonance Fingerprinting (MRF) reconstructs tissue maps based on a sequence of very highly undersampled images. In order to be able to perform MRF reconstruction, state-of-the-art MRF methods rely on priors such as the MR physics (Bloch equations) and might also use some additional low-rank or spatial regularization. However to our knowledge these three regularizations are not applied together in a joint reconstruction. The reason is that it is indeed challenging to incorporate effectively multiple regularizations in a single MRF optimization algorithm. As a result most of these methods are not robust to noise especially when the sequence length is short. In this paper, we propose a family of new methods where spatial and low-rank regularizations, in addition to the Bloch manifold regularization, are applied on the images. We show on digital phantom and NIST phantom scans, as well as volunteer scans that the proposed methods bring significant improvement in the quality of the estimated tissue maps.

Intracranial Distribution of Intravenously Administered Gadolinium-based Contrast Agent over a Period of 24 Hours: Evaluation with 3D-real IR Imaging and MR Fingerprinting.

PURPOSE: To evaluate the feasibility for the detection of slight contrast effects after intravenous administration of single dose gadolinium-based contrast agent (IV-SD-GBCA), the time course of the GBCA distribution up to 24 h was examined in various fluid spaces and brain parenchyma using 3D-real IR imaging and MR fingerprinting (MRF). METHODS: Twenty-four patients with a suspicion of endolymphatic hydrops were scanned at pre-administration and at 10 min, 4 and 24 h post-IV-SD-GBCA. 3D-real IR images and MRF at the level of the internal auditory canal were obtained. The signal intensity on the 3D-real IR image of the cerebrospinal fluid (CSF) in the cerebellopontine angle cistern (CPA), Sylvian fissure (Syl), lateral ventricle (LV), and cochlear perilymph (CPL) was measured. The T1 and T2 values of cerebellar gray (GM) and white matter (WM) were measured using MRF. Each averaged value at the various time points was compared using an analysis of variance. RESULTS: The signal intensity on the 3D-real IR image in each CSF region peaked at 4 h, and was decreased significantly by 24 h (P< 0.05). All patients had a maximum signal intensity at 4 h in the CPA, and Syl. The mean CPL signal intensity peaked at 4 h and decreased significantly by 24 h (P < 0.05). All patients but two had a maximum signal intensity at 4 h. Regarding the T1 value in the cerebellar WM and GM, the T1 value at 10 min post-IV-GBCA was significantly decreased compared to the pre-contrast scan, but no significant difference was observed at the other time points. There was no significant change in T2 in the gray or white matter at any of the time points. CONCLUSION: Time course of GBCA after IV-SD-GBCA could be evaluated by 3D-real IR imaging in CSF spaces and in the brain by MRF.

Detection of IV-gadolinium Leakage from the Cortical Veins into the CSF Using MR Fingerprinting.

PURPOSE: It has been reported that leakage of intravenously administered gadolinium-based contrast agents (IV-GBCAs) into the cerebrospinal fluid (CSF) from the cortical veins even in healthy subjects can be detected using a highly sensitive pulse sequence such as heavily T2-weighted 3D fluid-attenuated inversion recovery and 3D-real inversion recovery (IR). The purpose of this study was to evaluate the feasibility of MR fingerprinting to detect GBCA leakage from the cortical veins after IV-GBCA. MATERIALS: Fourteen patients with suspected endolymphatic hydrops (EH) who received a single dose of IV-GBCA (39-79 years old) were included. The real IR images as well as MR fingerprinting images were obtained at 4 h after IV-GBCA. T1 and T2 values were obtained using MR fingerprinting and analyzed in ROIs covering intense GBCA leakage, and non-leakage areas of the CSF as determined on real IR images. The scan time for real IR imaging was 10 min and that for MR fingerprinting was 41 s. RESULTS: The mean T1 value of the ROI in the area of GBCA leakage was 2422 ± 261 ms and that in the non-leakage area was 3851 ± 235 ms (P < 0.01). There was no overlap between the T1 values in the area of GBCA leakage and those in the non-leakage area.The mean T2 value in the area of GBCA leakage was 319 ± 90 ms and that in the non-leakage area was 670 ± 166 ms (P < 0.01). There was some overlap between the T2 values in the area of GBCA leakage and those in the non-leakage area. CONCLUSION: Leaked GBCA from the cortical veins into the surrounding CSF can be detected using MR fingerprinting obtained in <1 min.

Effect of spiral undersampling patterns on FISP MRF parameter maps.

PURPOSE: Artifacts arising from undersampling are not always treatable as incoherent noise for the pattern matching process in Magnetic Resonance Fingerprinting (MRF). To estimate the effect of undersampling artifacts on MRF quantitative results, spiral sampling trajectories and their temporal variation is examined. METHODS: The effect of sampling trajectories and their variation during the MRF experiment was assessed by characterizing aliasing artifacts. Temporal rearrangements of sampling trajectories were tested and evaluated in simulations and scans of phantoms and in a volunteer brain. RESULTS: Results show that some temporal variations of sampling patterns can lead to spatial biases in MRF parameter maps. Observed effects are consistent with derived performance indicators for different interleaving schemes, leading to substantially improved MRF sampling patterns. CONCLUSION: With the help of the presented simulation framework, MRF implementations can be investigated and improved. This was demonstrated for a spiral FISP (Fast imaging with steady-state free precession) MRF implementation, where a significantly improved interleaving scheme was identified, and confirmed by experiment.

Magnetic Resonance Fingerprinting Reconstruction Using Recurrent Neural Networks.

Magnetic Resonance Fingerprinting (MRF) is an imaging technique acquiring unique time signals for different tissues. Although the acquisition is highly accelerated, the reconstruction time remains a problem, as the state-of-the-art template matching compares every signal with a set of possible signals. To overcome this limitation, deep learning based approaches, e.g. Convolutional Neural Networks (CNNs) have been proposed. In this work, we investigate the applicability of Recurrent Neural Networks (RNNs) for this reconstruction problem, as the signals are correlated in time. Compared to previous methods based on CNNs, RNN models yield significantly improved results using in-vivo data.

Deep Learning for Magnetic Resonance Fingerprinting: A New Approach for Predicting Quantitative Parameter Values from Time Series.

The purpose of this work is to evaluate methods from deep learning for application to Magnetic Resonance Fingerprinting (MRF). MRF is a recently proposed measurement technique for generating quantitative parameter maps. In MRF a non-steady state signal is generated by a pseudo-random excitation pattern. A comparison of the measured signal in each voxel with the physical model yields quantitative parameter maps. Currently, the comparison is done by matching a dictionary of simulated signals to the acquired signals. To accelerate the computation of quantitative maps we train a Convolutional Neural Network (CNN) on simulated dictionary data. As a proof of principle we show that the neural network implicitly encodes the dictionary and can replace the matching process.