Robust sliding-window reconstruction for Accelerating the acquisition of MR fingerprinting

Reconstruction of multi-phase parametric maps in 4D-magnetic resonance fingerprinting (4D-MRF) by optimization of local T1 and T2 sensitivities

BACKGROUND

Time-resolved magnetic resonance fingerprinting (MRF), or 4D-MRF, has been demonstrated its feasibility in motion management in radiotherapy (RT). However, the prohibitive long acquisition time is one of challenges of the clinical implementation of 4D-MRF. The shortening of acquisition time causes data insufficiency in each respiratory phase, leading to poor accuracies and consistencies of the predicted tissues' properties of each phase.

PURPOSE

To develop a technique for the reconstruction of multi-phase parametric maps in four-dimensional magnetic resonance fingerprinting (4D-MRF) through the optimization of local T1 and T2 sensitivities.

METHODS

The proposed technique employed an iterative optimization to tailor the data arrangement of each phase by manipulation of inter-phase frames, such that the T1 and T2 sensitivities, which were quantified by the modified Minkowski distance, of the truncated signal evolution curve was maximized. The multi-phase signal evolution curves were modified by sliding window reconstruction and inter-phase frame sharing (SWIFS). Motion correction (MC) and dot product matching were sequentially performed on the modified signal evolution and dictionary to reconstruct the multi-parametric maps. The proposed technique was evaluated by numerical simulations using the extended cardiac-torso (XCAT) phantom with regular and irregular breathing patterns, and by in vivo MRF data of three health volunteers and six liver cancer patients acquired at a 3.0 T scanner.

RESULTS

In simulation study, the proposed SWIFS approach achieved the overall mean absolute percentage error (MAPE) of 8.62% ± 1.59% and 16.2% ± 3.88% for the eight-phases T1 and T2 maps, respectively, in the sagittal view with irregular breathing patterns. In contrast, the overall MAPE of T1 and T2 maps generated by the conventional approach with multiple MRF repetitions were 22.1% ± 11.0% and 30.8% ± 14.9%, respectively. For in-vivo study, the predicted mean T1 and T2 of liver by the proposed SWIFS approach were 795 ms ± 38.9 ms and 58.3 ms ± 11.7 ms, respectively.

CONCLUSIONS

Both simulation and in vivo results showed that the approach empowered by T1 and T2 sensitivities optimization and sliding window under the shortened acquisition of MRF had superior performance in the estimation of multi-phase T1 and T2 maps as compared to the conventional approach with oversampling of MRF data.

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Nonuniform sliding-window reconstruction for accelerated dual contrast agent quantification with MR fingerprinting

OBJECTIVE

MR fingerprinting (MRF) can enable preclinical studies of cell tracking by quantifying multiple contrast agents simultaneously, but faster scan times are required for in vivo applications. Sliding window (SW)-MRF is one option for accelerating MRF, but standard implementations are not sufficient to preserve the accuracy of T2*, which is critical for tracking iron-labelled cells in vivo.

PURPOSE

To develop a SW approach to MRF which preserves the T2* accuracy required for accelerated concentration mapping of iron-labelled cells on single-channel preclinical systems.

METHODS

A nonuniform SW was applied to the MRF sequence and dictionary. Segments of the sequence most sensitive to T2* were subject to a shorter window length, preserving the T2* sensitivity. Phantoms containing iron-labelled CD8+ T cells and gadolinium were used to compare 24× undersampled uniform and nonuniform SW-MRF parameter maps. Dual concentration maps were generated for both uniform and nonuniform MRF and compared.

RESULTS

Lin's concordance correlation coefficient, compared to gold standard parameter values, was much greater for nonuniform SW-MRF than for uniform SW-MRF. A Wilcoxon signed-rank test showed no significant difference between nonuniform SW-MRF and gold standards. Nonuniform SW-MRF outperformed the uniform SW-MRF concentration maps for all parameters, providing a balance between T2* sensitivity of short window lengths, and SNR of longer window lengths.

CONCLUSIONS

Nonuniform SW-MRF improves the accuracy of matching compared to uniform SW-MRF, allowing higher accelerated concentration mapping for preclinical systems.

Respiratory-Correlated 4-Dimensional Magnetic Resonance Fingerprinting for Liver Cancer Radiation Therapy Motion Management

PURPOSE

The objective of this study was to develop a respiratory-correlated (RC) 4-dimensional (4D) imaging technique based on magnetic resonance fingerprinting (MRF) (RC-4DMRF) for liver tumor motion management in radiation therapy.

METHODS AND MATERIALS

Thirteen patients with liver cancer were prospectively enrolled in this study. k-space MRF signals of the liver were acquired during free-breathing using the fast acquisition with steady-state precession sequence on a 3T scanner. The signals were binned into 8 respiratory phases based on respiratory surrogates, and interphase displacement vector fields were estimated using a phase-specific low-rank optimization method. Hereafter, the tissue property maps, including T1 and T2 relaxation times, and proton density, were reconstructed using a pyramid motion-compensated method that alternatively optimized interphase displacement vector fields and subspace images. To evaluate the efficacy of RC-4DMRF, amplitude motion differences and Pearson correlation coefficients were determined to assess measurement agreement in tumor motion between RC-4DMRF and cine magnetic resonance imaging (MRI); mean absolute percentage errors of the RC-4DMRF-derived tissue maps were calculated to reveal tissue quantification accuracy using digital human phantom; and tumor-to-liver contrast-to-noise ratio of RC-4DMRF images was compared with that of planning CT and contrast-enhanced MRI (CE-MRI) images. A paired Student t test was used for statistical significance analysis with a P value threshold of .05.

RESULTS

RC-4DMRF achieved excellent agreement in motion measurement with cine MRI, yielding the mean (± standard deviation) Pearson correlation coefficients of 0.95 ± 0.05 and 0.93 ± 0.09 and amplitude motion differences of 1.48 ± 1.06 mm and 0.81 ± 0.64 mm in the superior-inferior and anterior-posterior directions, respectively. Moreover, RC-4DMRF achieved high accuracy in tissue property quantification, with mean absolute percentage errors of 8.8%, 9.6%, and 5.0% for T1, T2, and proton density, respectively. Notably, the tumor contrast-to-noise ratio in RC-4DMRI-derived T1 maps (6.41 ± 3.37) was found to be the highest among all tissue property maps, approximately equal to that of CE-MRI (6.96 ± 1.01, P = .862), and substantially higher than that of planning CT (2.91 ± 1.97, P = .048).

CONCLUSIONS

RC-4DMRF demonstrated high accuracy in respiratory motion measurement and tissue properties quantification, potentially facilitating tumor motion management in liver radiation therapy.

Tailored magnetic resonance fingerprinting of post-operative pediatric brain tumor patients

PURPOSE

Brain and spinal cord tumors are the second most common cancer in children and account for one out of four cancers diagnosed. However, the long acquisition times associated with acquiring both data types prohibit using quantitative MR (qMR) in pediatric imaging protocols. This study aims to demonstrate the tailored magnetic resonance fingerprinting's (TMRF) ability to simultaneously provide quantitative maps (T1, T2) and multi-contrast qualitative images (T1 weighted, T1 FLAIR, T2 weighted) rapidly in pediatric brain tumor patients.

METHODS

In this work, we imaged five pediatric patients with brain tumors (resected/residual) using TMRF at 3 T. We compared the TMRF-derived T2 weighted images with those from the vendor-supplied sequence (as the gold standard, GS) for healthy and pathological tissue signal intensities. The relaxometric maps from TMRF were subjected to a region of interest (ROI) analysis to differentiate between healthy and pathological tissues. We performed the Wilcoxon rank sum test to check for significant differences between the two tissue types.

RESULTS

We found significant differences (p < 0.05) in both T1 and T2 ROI values between the two tissue types. A strong correlation was found between the TMRF-based T2 weighted and GS signal intensities for the healthy (correlation coefficient, r = 0.99) and pathological tissues (r = 0.88).

CONCLUSION

The TMRF implementation provides the two relaxometric maps and can potentially save ~2 min if it replaces the T2-weighted imaging in the current protocol.

Emerging Trends in Fast MRI Using Deep-Learning Reconstruction on Undersampled k-Space Data: A Systematic Review

Magnetic Resonance Imaging (MRI) is an essential medical imaging modality that provides excellent soft-tissue contrast and high-resolution images of the human body, allowing us to understand detailed information on morphology, structural integrity, and physiologic processes. However, MRI exams usually require lengthy acquisition times. Methods such as parallel MRI and Compressive Sensing (CS) have significantly reduced the MRI acquisition time by acquiring less data through undersampling k-space. The state-of-the-art of fast MRI has recently been redefined by integrating Deep Learning (DL) models with these undersampled approaches. This Systematic Literature Review (SLR) comprehensively analyzes deep MRI reconstruction models, emphasizing the key elements of recently proposed methods and highlighting their strengths and weaknesses. This SLR involves searching and selecting relevant studies from various databases, including Web of Science and Scopus, followed by a rigorous screening and data extraction process using the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. It focuses on various techniques, such as residual learning, image representation using encoders and decoders, data-consistency layers, unrolled networks, learned activations, attention modules, plug-and-play priors, diffusion models, and Bayesian methods. This SLR also discusses the use of loss functions and training with adversarial networks to enhance deep MRI reconstruction methods. Moreover, we explore various MRI reconstruction applications, including non-Cartesian reconstruction, super-resolution, dynamic MRI, joint learning of reconstruction with coil sensitivity and sampling, quantitative mapping, and MR fingerprinting. This paper also addresses research questions, provides insights for future directions, and emphasizes robust generalization and artifact handling. Therefore, this SLR serves as a valuable resource for advancing fast MRI, guiding research and development efforts of MRI reconstruction for better image quality and faster data acquisition.

Tailored magnetic resonance fingerprinting

Neuroimaging of certain pathologies requires both multi-parametric qualitative and quantitative imaging. The role of the quantitative MRI (qMRI) is well accepted but suffers from long acquisition times leading to patient discomfort, especially in geriatric and pediatric patients. Previous studies show that synthetic MRI can be used in order to reduce the scan time and provide qMRI as well as multi-contrast data. However, this approach suffers from artifacts such as partial volume and flow. In order to increase the scan efficiency (the number of contrasts and quantitative maps acquired per unit time), we designed, simulated, and demonstrated rapid, simultaneous, multi-contrast qualitative (T₁ weighted, T₁ fluid attenuated inversion recovery (FLAIR), T₂ weighted, water, and fat), and quantitative imaging (T₁ and T₂ maps) through the approach of tailored MR fingerprinting (TMRF) to cover whole-brain in approximately four minutes. We performed TMRF on in vivo four healthy human brains and in vitro ISMRM/NIST phantom and compared with vendor supplied gold standard (GS) and MRF sequences. All scans were performed on a 3 T GE Premier system and images were reconstructed offline using MATLAB. The reconstructed qualitative images were then subjected to custom DL denoising and gradient anisotropic diffusion denoising. The quantitative tissue parametric maps were reconstructed using a dense neural network to gain computational speed compared to dictionary matching. The grey matter and white matter tissues in qualitative and quantitative data for the in vivo datasets were segmented semi-automatically. The SNR and mean contrasts were plotted and compared across all three methods. The GS images show better SNR in all four subjects compared to MRF and TMRF (GS > TMRF>MRF). The T₁ and T₂ values of MRF are relatively overestimated as compared to GS and TMRF. The scan efficiency for TMRF is 1.72 min^-1 which is higher compared to GS (0.32 min^-1) and MRF (0.90 min^-1).

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.

Optimisation of data acquisition towards continuous cardiac Magnetic Resonance Fingerprinting applications

PURPOSE

Assess and optimise acquisition parameters for continuous cardiac Magnetic Resonance Fingerprinting (MRF).

METHODS

Different acquisition schemes (flip angle amplitude, lobe size, T2-preparation pulses) for cardiac MRF were assessed in simulations and phantom and demonstrated in one healthy volunteer. Three different experimental designs were evaluated using central composite and fractional factorial designs. Relative errors for T1 and T2 were calculated for a wide range of realistic T1 and T2 value combinations. The effect of different designs on the accuracy of T1 and T2 was assessed using response surface modelling and Cohen's f calculations.

RESULTS

Larger flip angle amplitudes lead to an improvement of T2 accuracy and precision for simulations and phantom experiments. Similar effects could also be shown qualitatively in in-vivo scans. Accuracy and precision of T1 were robust to different design parameters with improved values for faster flip angle variation. Cohen's f showed that T2-preparation pulses influence the accuracy of T2. The number of pulses used is the most important parameter. Without T2-preparation pulses, RMSE were 3.0 ± 8.09 % for T1 and 16.24 ± 14.47 % for T2. Using those pulses reduced the RMSE to 2.3 ± 8.4 % for T1 and 14.11 ± 13.46 % for T2. Nonetheless, even if the improvement is significant, RMSE are still too high for reliable quantification.

CONCLUSION

In contrast to previous study using triggered MRF sequences using < 30° flip angles, large flip angle amplitudes led to better results for continuous cardiac MRF sequences. T2-preparation pulse can improve the accuracy of T2 estimation but lead to longer scan times.

Artificial intelligence in cardiac magnetic resonance fingerprinting

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

MR fingerprinting of the prostate

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

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Technical overview of magnetic resonance fingerprinting and its applications in radiation therapy

Magnetic resonance fingerprinting (MRF) is an emerging imaging technique for rapid and simultaneous quantification of multiple tissue properties. The technique has been developed for quantitative imaging of different organs. The obtained quantitative measures have the potential to improve multiple steps of a typical radiotherapy workflow and potentially further improve integration of magnetic resonance imaging guided clinical decision making. In this review paper, we first provide a technical overview of the MRF method from data acquisition to postprocessing, along with recent development in advanced reconstruction methods. We further discuss critical aspects that could influence its usage in radiation therapy, such as accuracy and precision, repeatability and reproducibility, geometric distortion, and motion robustness. Finally, future directions for MRF application in radiation therapy are discussed.

Quantitative T1 and T2 mapping by magnetic resonance fingerprinting (MRF) of the placenta before and after maternal hyperoxia

INTRODUCTION

MR relaxometry has been used to assess placental exchange function, but methods to date are not sufficiently fast to be robust to placental motion. Magnetic resonance fingerprinting (MRF) permits rapid, voxel-wise, intrinsically co-registered T₁ and T₂ mapping. After characterizing measurement error, we scanned pregnant women during air and oxygen breathing to demonstrate MRF's ability to detect placental oxygenation changes.

METHODS

The accuracy of FISP-based, sliding-window reconstructed MRF was tested on phantoms. MRF scans in 9-s breath holds were acquired at 3T in 31 pregnant women during air and oxygen breathing. A mixed effects model was used to test for changes in placenta relaxation times between physiological states, to assess the dependency on gestational age (GA), and the impact of placental motion.

RESULTS

MRF estimates of known phantom relaxation times resulted in mean absolute errors for T₁ of 92 ms (4.8%), but T₂ was less accurate at 16 ms (13.6%). During normoxia, placental T₁ = 1825 ± 141 ms (avg ± standard deviation) and T₂ = 60 ± 16 ms (gestational age range 24.3-36.7, median 32.6 weeks). In the statistical model, placental T₂ rose and T₁ remained contant after hyperoxia, and no GA dependency was observed for T₁ or T₂.

DISCUSSION

Well-characterized, motion-robust MRF was used to acquire T₁ and T₂ maps of the placenta. Changes with hyperoxia are consistent with a net increase in oxygen saturation. Toward the goal of whole-placenta quantitative oxygenation imaging over time, we aim to implement 3D MRF with integrated motion correction to improve T₂ accuracy.

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.

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.

Rapid three-dimensional multiparametric MRI with quantitative transient-state imaging

Novel methods for quantitative, transient-state multiparametric imaging are increasingly being demonstrated for assessment of disease and treatment efficacy. Here, we build on these by assessing the most common Non-Cartesian readout trajectories (2D/3D radials and spirals), demonstrating efficient anti-aliasing with a k-space view-sharing technique, and proposing novel methods for parameter inference with neural networks that incorporate the estimation of proton density. Our results show good agreement with gold standard and phantom references for all readout trajectories at 1.5 T and 3 T. Parameters inferred with the neural network were within 6.58% difference from the parameters inferred with a high-resolution dictionary. Concordance correlation coefficients were above 0.92 and the normalized root mean squared error ranged between 4.2 and 12.7% with respect to gold-standard phantom references for T1 and T2. In vivo acquisitions demonstrate sub-millimetric isotropic resolution in under five minutes with reconstruction and inference times < 7 min. Our 3D quantitative transient-state imaging approach could enable high-resolution multiparametric tissue quantification within clinically acceptable acquisition and reconstruction times.

An Accurate Dictionary Creation Method for MR Fingerprinting Using a Fast Bloch Simulator

This study proposes an accurate method for creating a dictionary for magnetic resonance fingerprinting (MRF) using a fast Bloch image simulator. An MRF sequence based on a fast imaging with steady precession sequence and a numerical phantom were used for dictionary generation. Cartesian and spiral readout gradients were used for the Bloch image simulation. The validity and usefulness of the method for accurate dictionary creation were demonstrated by MRF parameter maps obtained by pattern matching with the dictionaries generated by the proposed method.

Retrospective rigid motion correction of three-dimensional magnetic resonance fingerprinting of the human brain

PURPOSE

To obtain three-dimensional (3D), quantitative and motion-robust imaging with magnetic resonance fingerprinting (MRF).

METHODS

Our acquisition is based on a 3D spiral projection k-space scheme. We compared different orderings of trajectory interleaves in terms of rigid motion-correction robustness. In all tested orderings, we considered the whole dataset as a sum of 56 segments of 7-s duration, acquired sequentially with the same flip angle schedule. We performed a separate image reconstruction for each segment, producing whole-brain navigators that were aligned to the first segment using normalized correlation. The estimated rigid motion was used to correct the k-space data, and the aligned data were matched with the dictionary to obtain motion-corrected maps.

RESULTS

A significant improvement on the motion-affected maps after motion correction is evident with the suppression of motion artifacts. Correlation with the motionless baseline improved by 20% on average for both T₁ and T₂ estimations after motion correction. In addition, the average motion-induced quantification bias of 70 ms for T₁ and 18 ms for T₂ values was reduced to 12 ms and 6 ms, respectively, improving the reliability of quantitative estimations.

CONCLUSION

We established a method that allows correcting 3D rigid motion on a 7-s timescale during the reconstruction of MRF data using self-navigators, improving the image quality and the quantification robustness.

Magnetic resonance fingerprinting review part 2: Technique and directions

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

Coupled Dictionary Learning for Multi-Contrast MRI Reconstruction

Magnetic resonance (MR) imaging tasks often involve multiple contrasts, such as T1-weighted, T2-weighted and fluid-attenuated inversion recovery (FLAIR) data. These contrasts capture information associated with the same underlying anatomy and thus exhibit similarities in either structure level or gray level. In this paper, we propose a coupled dictionary learning based multi-contrast MRI reconstruction (CDLMRI) approach to leverage the dependency correlation between different contrasts for guided or joint reconstruction from their under-sampled k -space data. Our approach iterates between three stages: coupled dictionary learning, coupled sparse denoising, and enforcing k -space consistency. The first stage learns a set of dictionaries that not only are adaptive to the contrasts, but also capture correlations among multiple contrasts in a sparse transform domain. By capitalizing on the learned dictionaries, the second stage performs coupled sparse coding to remove the aliasing and noise in the corrupted contrasts. The third stage enforces consistency between the denoised contrasts and the measurements in the k -space domain. Numerical experiments, consisting of retrospective under-sampling of various MRI contrasts with a variety of sampling schemes, demonstrate that CDLMRI is capable of capturing structural dependencies between different contrasts. The learned priors indicate notable advantages in multi-contrast MR imaging and promising applications in quantitative MR imaging such as MR fingerprinting.

Machine Learning for Rapid Magnetic Resonance Fingerprinting Tissue Property Quantification

Magnetic Resonance Fingerprinting (MRF) is an MRI-based method that can provide quantitative maps of multiple tissue properties simultaneously from a single rapid acquisition. Tissue property maps are generated by matching the complex signal evolutions collected at the scanner to a dictionary of signals derived using Bloch equation simulations. However, in some circumstances, the process of dictionary generation and signal matching can be time-consuming, reducing the utility of this technique. Recently, several groups have proposed using machine learning to accelerate the extraction of quantitative maps from MRF data. This article will provide an overview of current research that combines MRF and machine learning, as well as present original research demonstrating how machine learning can speed up dictionary generation for cardiac MRF.

Initial assessment of 3D magnetic resonance fingerprinting (MRF) towards quantitative brain imaging for radiation therapy

PURPOSE

Magnetic resonance fingerprinting (MRF) provides quantitative T1/T2 maps, enabling applications in clinical radiotherapy such as large-scale, multi-center clinical trials for longitudinal assessment of therapy response. We evaluated the feasibility of a quantitative three-dimensional-MRF (3D-MRF) towards its radiotherapy applications of primary brain tumors.

METHODS

A fast whole-brain 3D-MRF sequence initially developed for diagnostic radiology was optimized using flexible body coils, which is the typical MR imaging setup for radiotherapy treatment planning and for MR imaging (MRI)-guided treatment delivery. Optimization criteria included the accuracy and the precision of T1/T2 quantifications of polyvinylpyrrolidone (PVP) solutions, compared to those from the 3D-MRF using a 32-channel head coil. The accuracy of T1/T2 quantifications from the optimized MRF was first examined in healthy volunteers with two different coil setups. The intra- and inter-scanner variations of image intensity from the optimized sequence were quantified by longitudinal scans of the PVP solutions on two 3T scanners. Using a 3D-printed MRI geometry phantom, susceptibility-induced distortion with the optimized 3D-MRF was quantified as the Dice coefficient of phantom contours, compared to those from CT images. By introducing intentional head motion during 10% of the scan, the robustness of the optimized 3D-MRF towards motion was evaluated through visual inspection of motion artifacts and through quantitative analysis of image sharpness in brain MRF maps.

RESULTS

The optimized sequence acquired whole-brain T1, T2 and proton density maps and with a resolution of 1.2 × 1.2 × 3 mm³ in 10 min, similar to the total acquisition time of 3D T1- and T2-weighted images of the same resolution. In vivo T1 and T2 values of the white and gray matter were consistent with literature. The intra- and inter-scanner variability of the intensity-normalized MRF T1 was 1.0% ± 0.7% and 2.3% ± 1.0% respectively, in contrast to 5.3% ± 3.8% and 3.2% ± 1.6% from the normalized T1-weighted MRI. Repeatability and reproducibility of MRF T1 were independent of intensity normalization. Both phantom and human data demonstrated that the optimized 3D-MRF is more robust to subject motion and artifacts from subject-specific susceptibility difference. Compared to CT contours, the Dice coefficient of phantom contours from 3D-MRF was 0.93, improved from 0.87 from the T1-weighted MRI.

CONCLUSION

Compared to conventional MRI, the optimized 3D-MRF demonstrated improved repeatability across time points and reproducibility across scanners for better tissue quantification, as well as improved robustness to subject-specific susceptibility and motion artifacts under a typical MR imaging setup for radiotherapy. More importantly, quantitative MRF T1/T2 measurements lead to promising potentials towards longitudinal quantitative assessment of treatment response for better adaptive therapy and for large-scale, multi-center clinical trials.

Submillimeter MR fingerprinting using deep learning-based tissue quantification

PURPOSE

To develop a rapid 2D MR fingerprinting technique with a submillimeter in-plane resolution using a deep learning-based tissue quantification approach.

METHODS

A rapid and high-resolution MR fingerprinting technique was developed for brain T₁ and T₂ quantification. The 2D acquisition was performed using a FISP-based MR fingerprinting sequence and a spiral trajectory with 0.8-mm in-plane resolution. A deep learning-based method was used to replace the standard template matching method for improved tissue characterization. A novel network architecture (i.e., residual channel attention U-Net) was proposed to improve high-resolution details in the estimated tissue maps. Quantitative brain imaging was performed with 5 adults and 2 pediatric subjects, and the performance of the proposed approach was compared with several existing methods in the literature.

RESULTS

In vivo measurements with both adult and pediatric subjects show that high-quality T₁ and T₂ mapping with 0.8-mm in-plane resolution can be achieved in 7.5 seconds per slice. The proposed deep learning method outperformed existing algorithms in tissue quantification with improved accuracy. Compared with the standard U-Net, high-resolution details in brain tissues were better preserved by the proposed residual channel attention U-Net. Experiments on pediatric subjects further demonstrated the potential of the proposed technique for fast pediatric neuroimaging. Alongside reduced data acquisition time, a 5-fold acceleration in tissue property mapping was also achieved with the proposed method.

CONCLUSION

A rapid and high-resolution MR fingerprinting technique was developed, which enables high-quality T₁ and T₂ quantification with 0.8-mm in-plane resolution in 7.5 seconds per slice.

Placental MRI: Developing Accurate Quantitative Measures of Oxygenation

The Human Placenta Project has focused attention on the need for noninvasive magnetic resonance imaging (MRI)-based techniques to diagnose and monitor placental function throughout pregnancy. The hope is that the management of placenta-related pathologies would be improved if physicians had more direct, real-time measures of placental health to guide clinical decision making. As oxygen alters signal intensity on MRI and oxygen transport is a key function of the placenta, many of the MRI methods under development are focused on quantifying oxygen transport or oxygen content of the placenta. For example, measurements from blood oxygen level-dependent imaging of the placenta during maternal hyperoxia correspond to outcomes in twin pregnancies, suggesting that some aspects of placental oxygen transport can be monitored by MRI. Additional methods are being developed to accurately quantify baseline placental oxygenation by MRI relaxometry. However, direct validation of placental MRI methods is challenging and therefore animal studies and ex vivo studies of human placentas are needed. Here we provide an overview of the current state of the art of oxygen transport and quantification with MRI. We suggest that as these techniques are being developed, increased focus be placed on ensuring they are robust and reliable across individuals and standardized to enable predictive diagnostic models to be generated from the data. The field is still several years away from establishing the clinical benefit of monitoring placental function in real time with MRI, but the promise of individual personalized diagnosis and monitoring of placental disease in real time continues to motivate this effort.

Magnetic resonance fingerprinting of temporal lobe white matter in mesial temporal lobe epilepsy

Objective

Mesial temporal lobe epilepsy (MTLE) is a network disorder. We aimed to quantify the white matter alterations in the temporal lobe of MTLE patients with hippocampal sclerosis (MTLE‐HS) by using magnetic resonance fingerprinting (MRF), a novel imaging technique, which allows simultaneous measurements of multiple parameters with a single acquisition.

Methods

We consecutively recruited 27 unilateral MTLE‐HS patients and 22 healthy controls. Measurements including T1, T2, and PD values in the temporopolar white matter and temporal stem were recorded and analyzed.

Results

We found increased T2 value in both sides, and increased T1 value in the ipsilateral temporopolar white matter of MTLE‐HS patients, as compared with healthy controls. The T1 and T2 values were higher in the ipsilateral than the contralateral side. In the temporal stem, increased T1 and T2 values in the ipsilateral side of the MTLE‐HS patients were also observed. Only increased T2 values were observed in the contralateral temporal stem. No significant differences in PD values were observed in either the temporopolar white matter or temporal stem of the MTLE‐HS patients. Correlation analysis revealed that T1 and T2 values in the ipsilateral temporopolar white matter were negatively correlated with the age at epilepsy onset.

Interpretation

By using MRF, we were able to assess the alterations of T1 and T2 in the temporal lobe white matter of MTLE‐HS patients. MRF could be a promising imaging technique in identifying mild changes in MTLE patients, which might optimize the pre‐surgical evaluation and therapeutic interventions in these patients.

Magnetic resonance fingerprinting Part 1: Potential uses, current challenges, and recommendations

Magnetic resonance fingerprinting (MRF) is a powerful quantitative MRI technique capable of acquiring multiple property maps simultaneously in a short timeframe. The MRF framework has been adapted to a wide variety of clinical applications, but faces challenges in technical development, and to date has only demonstrated repeatability and reproducibility in small studies. In this review, we discuss the current implementations of MRF and their use in a clinical setting. Based on this analysis, we highlight areas of need that must be addressed before MRF can be fully adopted into the clinic and make recommendations to the MRF community on standardization and validation strategies of MRF techniques. Level of Evidence: 2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:675-692.

Feasibility of Quantitative Magnetic Resonance Fingerprinting in Ovarian Tumors for T1 and T2 Mapping in a PET/MR Setting

Multiparametric magnetic resonance imaging (MRI) can be used to characterize many cancer subtypes including ovarian cancer. Quantitative mapping of MRI relaxation values, such as T 1 and T 2 mapping, is promising for improving tumor assessment beyond conventional qualitative T 1- and T 2-weighted images. However, quantitative MRI relaxation mapping methods often involve long scan times due to sequentially measuring many parameters. Magnetic resonance fingerprinting (MRF) is a new method that enables fast quantitative MRI by exploiting the transient signals caused by the variation of pseudorandom sequence parameters. These transient signals are then matched to a simulated dictionary of T 1 and T 2 values to create quantitative maps. The ability of MRF to simultaneously measure multiple parameters, could represent a new approach to characterizing cancer and assessing treatment response. This feasibility study investigates MRF for simultaneous T 1, T 2, and relative proton density (rPD) mapping using ovarian cancer as a model system.

Sparsity and locally low rank regularization for MR fingerprinting

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) for simultaneous quantification of long and ultrashort T2 tissues

PURPOSE

To demonstrate an ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) method that allows quantifying relaxation times for muscle and bone in the musculoskeletal system and generating bone enhanced images that mimic CT scans.

METHODS

A fast imaging steady-state free precession MRF sequence with half pulse excitation and half projection readout was designed to sample fast T₂ decay signals. Varying echo time (TE) of a sinusoidal pattern was applied to enhance sensitivity for tissues with short and ultrashort T₂ values. The performance of UTE-MRF was evaluated via simulations, phantom, and in vivo experiments.

RESULTS

A minimal TE of 0.05 ms was achieved. Simulations indicated the sinusoidal TE sampling increased T₂ quantification accuracy in the cortical bone and tendon but had little impact on long T₂ muscle quantifications. For the rubber phantom, the averaged relaxometries from UTE-MRF (T₁ = 162 ms and T₂ = 1.07 ms) compared well with the gold standard (T₁ = 190 ms and T 2 ∗ = 1.03 ms). For the long T₂ agarose phantom, the linear regression slope between UTE-MRF and gold standard was 1.07 (R² = 0.991) for T₁ and 1.04 (R² = 0.994) for T₂ . In vivo experiments showed the detection of the cortical bone (averaged T₂ = 1.0 ms) and Achilles tendon (averaged T₂ = 15 ms). Scalp structures from the bone enhanced image show high similarity with CT.

CONCLUSION

The UTE-MRF with sinusoidal TEs can simultaneously quantify T₁ , T₂ , proton density, and B₀ in long, short, even ultrashort T₂ musculoskeletal structures. Bone enhanced images can be achieved in the brain with UTE-MRF.

Partial volume mapping using magnetic resonance fingerprinting

Magnetic resonance fingerprinting (MRF) is a quantitative imaging technique that maps multiple tissue properties through pseudorandom signal excitation and dictionary-based reconstruction. The aim of this study is to estimate and validate partial volumes from MRF signal evolutions (PV-MRF), and to characterize possible sources of error. Partial volume model inversion (pseudoinverse) and dictionary-matching approaches to calculate brain tissue fractions (cerebrospinal fluid, gray matter, white matter) were compared in a numerical phantom and seven healthy subjects scanned at 3 T. Results were validated by comparison with ground truth in simulations and ROI analysis in vivo. Simulations investigated tissue fraction errors arising from noise, undersampling artifacts, and model errors. An expanded partial volume model was investigated in a brain tumor patient. PV-MRF with dictionary matching is robust to noise, and estimated tissue fractions are sensitive to model errors. A 6% error in pure tissue T₁ resulted in average absolute tissue fraction error of 4% or less. A partial volume model within these accuracy limits could be semi-automatically constructed in vivo using k-means clustering of MRF-mapped relaxation times. Dictionary-based PV-MRF robustly identifies pure white matter, gray matter and cerebrospinal fluid, and partial volumes in subcortical structures. PV-MRF could also estimate partial volumes of solid tumor and peritumoral edema. We conclude that PV-MRF can attribute subtle changes in relaxation times to altered tissue composition, allowing for quantification of specific tissues which occupy a fraction of a voxel.

Fast 3D brain MR fingerprinting based on multi-axis spiral projection trajectory

PURPOSE

To develop a fast, sub-millimeter 3D magnetic resonance fingerprinting (MRF) technique for whole-brain quantitative scans.

METHODS

An acquisition trajectory based on multi-axis spiral projection imaging (maSPI) was implemented for 3D MRF with steady-state precession and slab excitation. By appropriately assigning the in-plane and through-plane rotations of spiral interleaves in a novel acquisition scheme, an maSPI-based acquisition was implemented, and the total acquisition time was reduced by up to a factor of 8 compared to stack-of-spiral (SOS)-based acquisition. A sliding-window method was also used to further reduce the required number of time points for a faster acquisition. The experiments were conducted both on a phantom and in vivo.

RESULTS

The results from the phantom measurements with the proposed and gold standard methods were consistent with a good linear correlation and an R² value approaching 0.99. The in vivo experiments achieved whole-brain parametric maps with isotropic resolutions of 1 mm and 0.8 mm in 5.0 and 6.0 min, respectively, with potential for further acceleration. An in vivo experiment with intentionally moving subjects demonstrated that the maSPI scheme largely outperforms the SOS scheme in terms of robustness to head motion.

CONCLUSION

3D MRF with an maSPI acquisition scheme enables fast and robust scans for high-resolution parametric mapping.

Optimal Experiment Design for Magnetic Resonance Fingerprinting: Cramér-Rao Bound Meets Spin Dynamics

Magnetic resonance (MR) fingerprinting is a new quantitative imaging paradigm, which simultaneously acquires multiple MR tissue parameter maps in a single experiment. In this paper, we present an estimation-theoretic framework to perform experiment design for MR fingerprinting. Specifically, we describe a discrete-time dynamic system to model spin dynamics, and derive an estimation-theoretic bound, i.e., the Cramér-Rao bound, to characterize the signal-to-noise ratio (SNR) efficiency of an MR fingerprinting experiment. We then formulate an optimal experiment design problem, which determines a sequence of acquisition parameters to encode MR tissue parameters with the maximal SNR efficiency, while respecting the physical constraints and other constraints from the image decoding/reconstruction process. We evaluate the performance of the proposed approach with numerical simulations, phantom experiments, and in vivo experiments. We demonstrate that the optimized experiments substantially reduce data acquisition time and/or improve parameter estimation. For example, the optimized experiments achieve about a factor of two improvement in the accuracy of T2 maps, while keeping similar or slightly better accuracy of T1 maps. Finally, as a remarkable observation, we find that the sequence of optimized acquisition parameters appears to be highly structured rather than randomly/pseudo-randomly varying as is prescribed in the conventional MR fingerprinting experiments.

Magnetic Resonance Fingerprinting: Implications and Opportunities for PET/MR

Magnetic Resonance Imaging (MRI) can be used to assess anatomical structure, and its sensitivity to a variety of tissue properties enables superb contrast between tissues as well as the ability to characterize these tissues. However, despite vast potential for quantitative and functional evaluation, MRI is typically used qualitatively, in which the underlying tissue properties are not measured, and thus the brightness of each pixel is not quantitatively meaningful. Positron Emission Tomography (PET) is an inherently quantitative imaging modality that interrogates functional activity within a tissue, probed by a molecule of interest coupled with an appropriate tracer. These modalities can complement one another to provide clinical information regarding both structure and function, but there are still technical and practical hurdles in the way of the integrated use of both modalities. Recent advances in MRI have moved the field in an increasingly quantitative direction, which is complementary to PET, and could also potentially help solve some of the challenges in PET/MR. Magnetic Resonance Fingerprinting (MRF) is a recently described MRI-based technique which can efficiently and simultaneously quantitatively map several tissue properties in a single exam. Here, the basic principles behind the quantitative approach of MRF are laid out, and the potential implications for combined PET/MR are discussed.

Rigid motion-corrected magnetic resonance fingerprinting

PURPOSE

Develop a method for rigid body motion-corrected magnetic resonance fingerprinting (MRF).

METHODS

MRF has shown some robustness to abrupt motion toward the end of the acquisition. Here, we study the effects of different types of rigid body motion during the acquisition on MRF and propose a novel approach to correct for this motion. The proposed method (MC-MRF) follows 4 steps: (1) sliding window reconstruction is performed to produce high-quality auxiliary dynamic images; (2) rotation and translation motion is estimated from the dynamic images by image registration; (3) estimated motion is used to correct acquired k-space data with corresponding rotations and phase shifts; and (4) motion-corrected data are reconstructed with low-rank inversion. MC-MRF was validated in a standard T1 /T2 phantom and 2D in vivo brain acquisitions in 7 healthy subjects. Additionally, the effect of through-plane motion in 2D MC-MRF was investigated.

RESULTS

Simulation results show that motion in MRF can introduce artifacts in T1 and T2 maps, depending when it occurs. MC-MRF improved parametric map quality in all phantom and in vivo experiments with in-plane motion, comparable to the no-motion ground truth. Reduced parametric map quality, even after motion correction, was observed for acquisitions with through-plane motion, particularly for smaller structures in T2 maps.

CONCLUSION

Here, a novel method for motion correction in MRF (MC-MRF) is proposed, which improves parametric map quality and accuracy in comparison to no-motion correction approaches. Future work will include validation of 3D MC-MRF to enable also through-plane motion correction.

Quantifying amide proton exchange rate and concentration in chemical exchange saturation transfer imaging of the human brain

Current chemical exchange saturation transfer (CEST) neuroimaging protocols typically acquire CEST-weighted images, and, as such, do not essentially provide quantitative proton-specific exchange rates (or brain pH) and concentrations. We developed a dictionary-free MR fingerprinting (MRF) technique to allow CEST parameter quantification with a reduced data set. This was accomplished by subgrouping proton exchange models (SPEM), taking amide proton transfer (APT) as an example, into two-pool (water and semisolid macromolecules) and three-pool (water, semisolid macromolecules, and amide protons) models. A variable radiofrequency saturation scheme was used to generate unique signal evolutions for different tissues, reflecting their CEST parameters. The proposed MRF-SPEM method was validated using Bloch-McConnell equation-based digital phantoms with known ground-truth, which showed that MRF-SPEM can achieve a high degree of accuracy and precision for absolute CEST parameter quantification and CEST phantoms. For in-vivo studies at 3 T, using the same model as in the simulations, synthetic Z-spectra were generated using rates and concentrations estimated from the MRF-SPEM reconstruction and compared with experimentally measured Z-spectra as the standard for optimization. The MRF-SPEM technique can provide rapid and quantitative human brain CEST mapping.

Rigid motion correction for magnetic resonance fingerprinting with sliding-window reconstruction and image registration

Magnetic resonance fingerprinting (MRF) can be used to simultaneously obtain multiple parameter maps from a single pulse sequence. However, patient motion during MRF acquisition may result in blurring and artifacts in estimated parameter maps. In this work, a novel motion correction method was proposed to correct for rigid motion in MRF. The proposed method involved sliding-window reconstruction to obtain intermediate images followed by image registration to estimate rigid motion information between these images. Finally, the motion-corrupted k-space data were corrected with the estimated motion parameters and then reconstructed to obtain the parameter maps via the conventional MRF processing pipeline. The proposed method was evaluated using both simulations and in vivo MRF experiments with intently different types of motion. For motion-corrupted data, the proposed method yielded brain T₁, T₂ and proton density maps with obviously reduced blurring and artifacts and lower normalized root-mean-square error, compared to MRF without motion correction. In conclusion, motion-corrected MRF using the proposed method has the potential to produce accurate parameter maps in the presence of in-plane rigid motion.

MR fingerprinting Deep RecOnstruction NEtwork (DRONE)

PURPOSE

Demonstrate a novel fast method for reconstruction of multi-dimensional MR Fingerprinting (MRF) data using Deep Learning methods.

METHODS

A neural network (NN) is defined using the TensorFlow framework and trained on simulated MRF data computed with the Extended Phase Graph formalism. The NN reconstruction accuracy for noiseless and noisy data is compared to conventional MRF template matching as a function of training data size, and quantified in simulated numerical brain phantom data and ISMRM/NIST phantom data measured on 1.5T and 3T scanners with an optimized MRF EPI and MRF FISP sequences with spiral readout. The utility of the method is demonstrated in a healthy subject in vivo at 1.5 T.

RESULTS

Network training required 10 to 74 minutes and once trained, data reconstruction required approximately 10 ms for the MRF EPI and 76 ms for the MRF FISP sequence. Reconstruction of simulated, noiseless brain data using the NN resulted in a root-mean-square error (RMSE) of 2.6 ms for T1 and 1.9 ms for T2. The reconstruction error in the presence of noise was less than 10% for both T1 and T2 for signal-to-noise greater than 25 dB. Phantom measurements yielded good agreement (R2=0.99/0.99 for MRF EPI T1/T2 and 0.94/0.98 for MRF FISP T1/T2) between the T1 and T2 estimated by the NN and reference values from the ISMRM/NIST phantom.

CONCLUSION

Reconstruction of MRF data with a NN is accurate, 300–5000 fold faster and more robust to noise and undersampling than conventional MRF dictionary matching.

Accelerated magnetic resonance fingerprinting using soft-weighted key-hole (MRF-SOHO)

Object

To develop a novel approach for highly accelerated Magnetic Resonance Fingerprinting (MRF) acquisition.

Materials and methods

The proposed method combines parallel imaging, soft-gating and key-hole approaches to highly accelerate MRF acquisition. Slowly varying flip angles (FA), commonly used during MRF acquisition, lead to a smooth change in the signal contrast of consecutive time-point images. This assumption enables sharing of high frequency data between different time-points, similar to what is done in some dynamic MR imaging methods such as key-hole. The proposed approach exploits this information using a SOft-weighted key-HOle (MRF-SOHO) reconstruction to achieve high acceleration factors and/or increased resolution without compromising image quality or increasing scan time. MRF-SOHO was validated on a standard T₁/T₂ phantom and in in-vivo brain acquisitions reconstructing T₁, T₂ and proton density parametric maps.

Results

Accelerated MRF-SOHO using less data per time-point and less time-point images enabled a considerable reduction in scan time (up to 4.6x), while obtaining similar T₁ and T₂ accuracy and precision when compared to zero-filled MRF reconstruction. For the same number of spokes and time-points, the proposed method yielded an enhanced performance in quantifying parameters than the zero-filled MRF reconstruction, which was verified with 2, 1 and 0.7 (sub-millimetre) resolutions.

Conclusion

The proposed MRF-SOHO enabled a 4.6x scan time reduction for an in-plane spatial resolution of 2x2 mm² when compared to zero-filled MRF and enabled sub-millimetric (0.7x0.7 mm²) resolution MRF.

Detection of Lesions in Mesial Temporal Lobe Epilepsy by Using MR Fingerprinting

Purpose To improve diagnosis of hippocampal sclerosis (HS) in patients with mesial temporal lobe epilepsy (MTLE) by using MR fingerprinting and compare with visual assessment of T1- and T2-weighted MR images. Materials and Methods For this prospective study performed between April and November 2016, T1 and T2 maps were obtained and tissue segmentation performed in consecutive patients with drug-resistant MTLE with unilateral or bilateral HS. T1 and T2 maps were compared between 33 patients with MTLE (23 women and 10 men; mean age, 32.6 years; age range, 16-60 years) and 30 healthy participants (20 women and 10 men; mean age, 28.8 years; age range, 18-40 years). Differences in individual bilateral hippocampi were compared by using a Wilcoxon signed rank test, whereas the Wilcoxon rank-sum test was used for difference analysis between healthy control participants and patients with MTLE. Results The diagnosis rate (ie, ratio of HS diagnosed on the basis of a 2.5-minute MR fingerprinting examination compared with standard methods: MRI, electroencephalography, and PET) was 32 of 33 (96.9%; 95% confidence interval: 84.9%, 100%), reflecting improved accuracy of diagnosis (P = 1.92 × 10^-12) over routine MR examinations that had a diagnostic rate of 23 of 33 (69.7%; 95% confidence interval: 51.5%, 81.6%). The comparison between atrophic and normal-appearing hippocampus in 33 patients with MTLE and healthy control participants demonstrated that both T1 and T2 values in HS lesions were higher than those of normal hippocampal tissue of healthy participants (T1: 1361 msec ± 85 vs 1249 msec ± 59, respectively; T2: 135 msec ± 15 vs 104 msec ± 9, respectively; P < .0001). Conclusion MR fingerprinting allowed for multiparametric mapping of temporal lobe within 2.5 minutes and helped to identify lesions suspicious for HS in patients with MTLE with improved accuracy.

3D MR fingerprinting with accelerated stack-of-spirals and hybrid sliding-window and GRAPPA reconstruction

PURPOSE

Whole-brain high-resolution quantitative imaging is extremely encoding intensive, and its rapid and robust acquisition remains a challenge. Here we present a 3D MR fingerprinting (MRF) acquisition with a hybrid sliding-window (SW) and GRAPPA reconstruction strategy to obtain high-resolution T1, T2 and proton density (PD) maps with whole brain coverage in a clinically feasible timeframe.

METHODS

3D MRF data were acquired using a highly under-sampled stack-of-spirals trajectory with a steady-state precession (FISP) sequence. For data reconstruction, kx-ky under-sampling was mitigated using SW combination along the temporal axis. Non-uniform fast Fourier transform (NUFFT) was then applied to create Cartesian k-space data that are fully-sampled in the in-plane direction, and Cartesian GRAPPA was performed to resolve kz under-sampling to create an alias-free SW dataset. T1, T2 and PD maps were then obtained using dictionary matching.

RESULTS

Phantom study demonstrated that the proposed 3D-MRF acquisition/reconstruction method is able to produce quantitative maps that are consistent with conventional quantification techniques. Retrospectively under-sampled in vivo acquisition revealed that SW + GRAPPA substantially improves quantification accuracy over the current state-of-the-art accelerated 3D MRF. Prospectively under-sampled in vivo study showed that whole brain T1, T2 and PD maps with 1 mm3 resolution could be obtained in 7.5 min.

CONCLUSIONS

3D MRF stack-of-spirals acquisition with hybrid SW + GRAPPA reconstruction may provide a feasible approach for rapid, high-resolution quantitative whole-brain imaging.