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.

A fast MR fingerprinting simulator for direct error estimation and sequence optimization.

Magnetic resonance fingerprinting (MRF) is a novel quantitative MR technique that simultaneously provides multiple tissue property maps. When optimizing MRF scans, modeling undersampling errors and field imperfections in cost functions for direct measurement of quantitative errors will make the optimization results more practical and robust. However, optimizing such cost function is computationally expensive and impractical for MRF optimization with tens of thousands of iterations. Here, we introduce a fast MRF simulator to simulate aliased images from actual scan scenarios including undersampling and system imperfections, which substantially reduces computational time and allows for direct error estimation of the quantitative maps and efficient sequence optimization. We evaluate the performance and computational speed of the proposed approach by simulations and in vivo experiments. The simulations from the proposed method closely approximate the signals and MRF maps from in vivo scans, with 158 times shorter processing time than the conventional simulation method using Non-uniform Fourier transform. We also demonstrate the power of applying the fast MRF simulator in MRF sequence optimization. The optimized sequences are validated with in vivo scans to assess the image quality and accuracy. The optimized sequences produce artifact-free T1 and T2 maps in 2D and 3D scans with equivalent mapping accuracy as the human-designed sequence but at shorter scan times. Incorporating the proposed simulator in the MRF optimization framework makes direct estimation of undersampling errors during the optimization process feasible, and provide optimized MRF sequences that are robust against undersampling artifacts and field inhomogeneity.

Automated design of pulse sequences for magnetic resonance fingerprinting using physics-inspired optimization.

Magnetic resonance fingerprinting (MRF) is a method to extract quantitative tissue properties such as [Formula: see text] and [Formula: see text] relaxation rates from arbitrary pulse sequences using conventional MRI hardware. MRF pulse sequences have thousands of tunable parameters, which can be chosen to maximize precision and minimize scan time. Here, we perform de novo automated design of MRF pulse sequences by applying physics-inspired optimization heuristics. Our experimental data suggest that systematic errors dominate over random errors in MRF scans under clinically relevant conditions of high undersampling. Thus, in contrast to prior optimization efforts, which focused on statistical error models, we use a cost function based on explicit first-principles simulation of systematic errors arising from Fourier undersampling and phase variation. The resulting pulse sequences display features qualitatively different from previously used MRF pulse sequences and achieve fourfold shorter scan time than prior human-designed sequences of equivalent precision in [Formula: see text] and [Formula: see text] Furthermore, the optimization algorithm has discovered the existence of MRF pulse sequences with intrinsic robustness against shading artifacts due to phase variation.

Radiomic analysis of magnetic resonance fingerprinting in adult brain tumors.

PURPOSE: This is a radiomics study investigating the ability of texture analysis of MRF maps to improve differentiation between intra-axial adult brain tumors and to predict survival in the glioblastoma cohort. METHODS: Magnetic resonance fingerprinting (MRF) acquisition was performed on 31 patients across 3 groups: 17 glioblastomas, 6 low-grade gliomas, and 8 metastases. Using regions of interest for the solid tumor and peritumoral white matter on T1 and T2 maps, second-order texture features were calculated from gray-level co-occurrence matrices and gray-level run length matrices. Selected features were compared across the three tumor groups using Wilcoxon rank-sum test. Receiver operating characteristic curve analysis was performed for each feature. Kaplan-Meier method was used for survival analysis with log rank tests. RESULTS: Low-grade gliomas and glioblastomas had significantly higher run percentage, run entropy, and information measure of correlation 1 on T1 than metastases (p < 0.017). The best separation of all three tumor types was seen utilizing inverse difference normalized and homogeneity values for peritumoral white matter in both T1 and T2 maps (p < 0.017). In solid tumor T2 maps, lower values in entropy and higher values of maximum probability and high-gray run emphasis were associated with longer survival in glioblastoma patients (p < 0.05). Several texture features were associated with longer survival in glioblastoma patients on peritumoral white matter T1 maps (p < 0.05). CONCLUSION: Texture analysis of MRF-derived maps can improve our ability to differentiate common adult brain tumors by characterizing tumor heterogeneity, and may have a role in predicting outcomes in patients with glioblastoma.

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 T1 , T2 , ΔB0 , and T2∗ . METHODS: 3D MRF data with effective undersampling factor of 3 in the slice direction were acquired with quadratic RF phase patterns for T1 , T2 , and T2∗ 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 T1 and T2 contrast and to 2D qRF-MRF for T2∗ contrast and ΔB0 . The precision of 3D qRF-MRF was 1.5-2 times higher than routine clinical scans. 3D qRF-MRF ΔB0 maps were further processed to highlight the susceptibility contrast. CONCLUSION: Natively co-registered 3D whole brain T1 , T2 , T2∗ , ΔB0 , and QSM maps can be acquired in as short as 5 min with 3D qRF-MRF.

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.

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.

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

Magnetic resonance fingerprinting: a technical review.

Multiparametric quantitative imaging is gaining increasing interest due to its widespread advantages in clinical applications. Magnetic resonance fingerprinting is a recently introduced approach of fast multiparametric quantitative imaging. In this article, magnetic resonance fingerprinting acquisition, dictionary generation, reconstruction, and validation are reviewed.

Three-dimensional MR Fingerprinting for Quantitative Breast Imaging.

Purpose To develop a fast three-dimensional method for simultaneous T1 and T2 quantification for breast imaging by using MR fingerprinting. Materials and Methods In this prospective study, variable flip angles and magnetization preparation modules were applied to acquire MR fingerprinting data for each partition of a three-dimensional data set. A fast postprocessing method was implemented by using singular value decomposition. The proposed technique was first validated in phantoms and then applied to 15 healthy female participants (mean age, 24.2 years ± 5.1 [standard deviation]; range, 18-35 years) and 14 female participants with breast cancer (mean age, 55.4 years ± 8.8; range, 39-66 years) between March 2016 and April 2018. The sensitivity of the method to B1 field inhomogeneity was also evaluated by using the Bloch-Siegert method. Results Phantom results showed that accurate and volumetric T1 and T2 quantification was achieved by using the proposed technique. The acquisition time for three-dimensional quantitative maps with a spatial resolution of 1.6 × 1.6 × 3 mm3 was approximately 6 minutes. For healthy participants, averaged T1 and T2 relaxation times for fibroglandular tissues at 3.0 T were 1256 msec ± 171 and 46 msec ± 7, respectively. Compared with normal breast tissues, higher T2 relaxation time (68 msec ± 13) was observed in invasive ductal carcinoma (P < .001), whereas no statistical difference was found in T1 relaxation time (1183 msec ± 256; P = .37). Conclusion A method was developed for breast imaging by using the MR fingerprinting technique, which allows simultaneous and volumetric quantification of T1 and T2 relaxation times for breast tissues. © RSNA, 2018 Online supplemental material is available for this article.

Development of high-resolution 3D MR fingerprinting for detection and characterization of epileptic lesions.

BACKGROUND: Conventional MRI can be limited in detecting subtle epileptic lesions or identifying active/epileptic lesions among widespread, multifocal lesions. PURPOSE: We developed a high-resolution 3D MR fingerprinting (MRF) protocol to simultaneously provide quantitative T1 , T2 , proton density, and tissue fraction maps for detection and characterization of epileptic lesions. STUDY TYPE: Prospective. POPULATION: National Institute of Standards and Technology (NIST) / International Society for Magnetic Resonance in Medicine (ISMRM) phantom, five healthy volunteers and 15 patients with medically intractable epilepsy undergoing presurgical evaluation with noninvasive or invasive electroclinical data. FIELD STRENGTH/SEQUENCE: 3D MRF scans and routine clinical epilepsy MR protocols were acquired at 3 T. ASSESSMENT: The accuracy of the T1 and T2 values were first evaluated using the NIST/ISMRM phantom. The repeatability was then estimated with both phantom and volunteers based on the coefficient of variance (CV). For epilepsy patients, all the maps were qualitatively reviewed for lesion detection by three independent reviewers (S.E.J., M.L., I.N.) blinded to clinical data. Region of interest (ROI) analysis was performed on T1 and T2 maps to quantify the multiparametric signal differences between lesion and normal tissues. Findings from qualitative review and quantitative ROI analysis were compared with patients' electroclinical data to assess concordance. STATISTICAL TESTS: Phantom results were compared using R-squared, and patient results were compared using linear regression models. RESULTS: The phantom study showed high accuracy with the standard values, with an R2 of 0.99. The volunteer study showed high repeatability, with an average CV of 4.3% for T1 and T2 in various tissue regions. For the 15 patients, MRF showed additional findings in four patients, with the remaining 11 patients showing findings consistent with conventional MRI. The additional MRF findings were highly concordant with patients' electroclinical presentation. DATA CONCLUSION: The 3D MRF protocol showed potential to identify otherwise inconspicuous epileptogenic lesions from the patients with negative conventional MRI diagnosis, as well as to correlate with different levels of epileptogenicity when widespread lesions were present. LEVEL OF EVIDENCE: 3. Technical Efficacy Stage: 3. J. Magn. Reson. Imaging 2019;49:1333-1346.

Fast 3D magnetic resonance fingerprinting for a whole-brain coverage.

PURPOSE: The purpose of this study was to accelerate the acquisition and reconstruction time of 3D magnetic resonance fingerprinting scans. METHODS: A 3D magnetic resonance fingerprinting scan was accelerated by using a single-shot spiral trajectory with an undersampling factor of 48 in the x-y plane, and an interleaved sampling pattern with an undersampling factor of 3 through plane. Further acceleration came from reducing the waiting time between neighboring partitions. The reconstruction time was accelerated by applying singular value decomposition compression in k-space. Finally, a 3D premeasured B1 map was used to correct for the B1 inhomogeneity. RESULTS: The T1 and T2 values of the International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology MRI phantom showed a good agreement with the standard values, with an average concordance correlation coefficient of 0.99, and coefficient of variation of 7% in the repeatability scans. The results from in vivo scans also showed high image quality in both transverse and coronal views. CONCLUSIONS: This study applied a fast acquisition scheme for a fully quantitative 3D magnetic resonance fingerprinting scan with a total acceleration factor of 144 as compared with the Nyquist rate, such that 3D T1 , T2 , and proton density maps can be acquired with whole-brain coverage at clinical resolution in less than 5 min. Magn Reson Med 79:2190-2197, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Bayesian estimation of multicomponent relaxation parameters in magnetic resonance fingerprinting.

PURPOSE: To estimate multiple components within a single voxel in magnetic resonance fingerprinting when the number and types of tissues comprising the voxel are not known a priori. THEORY: Multiple tissue components within a single voxel are potentially separable with magnetic resonance fingerprinting as a result of differences in signal evolutions of each component. The Bayesian framework for inverse problems provides a natural and flexible setting for solving this problem when the tissue composition per voxel is unknown. Assuming that only a few entries from the dictionary contribute to a mixed signal, sparsity-promoting priors can be placed upon the solution. METHODS: An iterative algorithm is applied to compute the maximum a posteriori estimator of the posterior probability density to determine the magnetic resonance fingerprinting dictionary entries that contribute most significantly to mixed or pure voxels. RESULTS: Simulation results show that the algorithm is robust in finding the component tissues of mixed voxels. Preliminary in vivo data confirm this result, and show good agreement in voxels containing pure tissue. CONCLUSIONS: The Bayesian framework and algorithm shown provide accurate solutions for the partial-volume problem in magnetic resonance fingerprinting. The flexibility of the method will allow further study into different priors and hyperpriors that can be applied in the model. Magn Reson Med 80:159-170, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Low rank approximation methods for MR fingerprinting with large scale dictionaries.

PURPOSE: This work proposes new low rank approximation approaches with significant memory savings for large scale MR fingerprinting (MRF) problems. THEORY AND METHODS: We introduce a compressed MRF with randomized singular value decomposition method to significantly reduce the memory requirement for calculating a low rank approximation of large sized MRF dictionaries. We further relax this requirement by exploiting the structures of MRF dictionaries in the randomized singular value decomposition space and fitting them to low-degree polynomials to generate high resolution MRF parameter maps. In vivo 1.5T and 3T brain scan data are used to validate the approaches. RESULTS: T1 , T2 , and off-resonance maps are in good agreement with that of the standard MRF approach. Moreover, the memory savings is up to 1000 times for the MRF-fast imaging with steady-state precession sequence and more than 15 times for the MRF-balanced, steady-state free precession sequence. CONCLUSION: The proposed compressed MRF with randomized singular value decomposition and dictionary fitting methods are memory efficient low rank approximation methods, which can benefit the usage of MRF in clinical settings. They also have great potentials in large scale MRF problems, such as problems considering multi-component MRF parameters or high resolution in the parameter space. Magn Reson Med 79:2392-2400, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Development of a Combined MR Fingerprinting and Diffusion Examination for Prostate Cancer.

Purpose To develop and evaluate an examination consisting of magnetic resonance (MR) fingerprinting-based T1, T2, and standard apparent diffusion coefficient (ADC) mapping for multiparametric characterization of prostate disease. Materials and Methods This institutional review board-approved, HIPAA-compliant retrospective study of prospectively collected data included 140 patients suspected of having prostate cancer. T1 and T2 mapping was performed with fast imaging with steady-state precession-based MR fingerprinting with ADC mapping. Regions of interest were drawn by two independent readers in peripheral zone lesions and normal-appearing peripheral zone (NPZ) tissue identified on clinical images. T1, T2, and ADC were recorded for each region. Histopathologic correlation was based on systematic transrectal biopsy or cognitively targeted biopsy results, if available. Generalized estimating equations logistic regression was used to assess T1, T2, and ADC in the differentiation of (a) cancer versus NPZ, (b) cancer versus prostatitis, (c) prostatitis versus NPZ, and (d) high- or intermediate-grade tumors versus low-grade tumors. Analysis was performed for all lesions and repeated in a targeted biopsy subset. Discriminating ability was evaluated by using the area under the receiver operating characteristic curve (AUC). Results In this study, 109 lesions were analyzed, including 39 with cognitively targeted sampling. T1, T2, and ADC from cancer (mean, 1628 msec ± 344, 73 msec ± 27, and 0.773 × 10-3 mm2/sec ± 0.331, respectively) were significantly lower than those from NPZ (mean, 2247 msec ± 450, 169 msec ± 61, and 1.711 × 10-3 mm2/sec ± 0.269) (P < .0001 for each) and together produced the best separation between these groups (AUC = 0.99). ADC and T2 together produced the highest AUC of 0.83 for separating high- or intermediate-grade tumors from low-grade cancers. T1, T2, and ADC in prostatitis (mean, 1707 msec ± 377, 79 msec ± 37, and 0.911 × 10-3 mm2/sec ± 0.239) were significantly lower than those in NPZ (P < .0005 for each). Interreader agreement was excellent, with an intraclass correlation coefficient greater than 0.75 for both T1 and T2 measurements. Conclusion This study describes the development of a rapid MR fingerprinting- and diffusion-based examination for quantitative characterization of prostatic tissue. © RSNA, 2017 Online supplemental material is available for this article.

Slice profile and B1 corrections in 2D magnetic resonance fingerprinting.

PURPOSE: The goal of this study is to characterize and improve the accuracy of 2D magnetic resonance fingerprinting (MRF) scans in the presence of slice profile (SP) and B1 imperfections, which are two main factors that affect quantitative results in MRF. METHODS: The SP and B1 imperfections are characterized and corrected separately. The SP effect is corrected by simulating the radiofrequency pulse in the dictionary, and the B1 is corrected by acquiring a B1 map using the Bloch-Siegert method before each scan. The accuracy, precision, and repeatability of the proposed method are evaluated in phantom studies. The effects of both SP and B1 imperfections are also illustrated and corrected in the in vivo studies. RESULTS: The SP and B1 corrections improve the accuracy of the T1 and T2 values, independent of the shape of the radiofrequency pulse. The T1 and T2 values obtained from different excitation patterns become more consistent after corrections, which leads to an improvement of the robustness of the MRF design. CONCLUSION: This study demonstrates that MRF is sensitive to both SP and B1 effects, and that corrections can be made to improve the accuracy of MRF with only a 2-s increase in acquisition time. Magn Reson Med 78:1781-1789, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

SVD compression for magnetic resonance fingerprinting in the time domain.

Magnetic resonance (MR) fingerprinting is a technique for acquiring and processing MR data that simultaneously provides quantitative maps of different tissue parameters through a pattern recognition algorithm. A predefined dictionary models the possible signal evolutions simulated using the Bloch equations with different combinations of various MR parameters and pattern recognition is completed by computing the inner product between the observed signal and each of the predicted signals within the dictionary. Though this matching algorithm has been shown to accurately predict the MR parameters of interest, one desires a more efficient method to obtain the quantitative images. We propose to compress the dictionary using the singular value decomposition, which will provide a low-rank approximation. By compressing the size of the dictionary in the time domain, we are able to speed up the pattern recognition algorithm, by a factor of between 3.4-4.8, without sacrificing the high signal-to-noise ratio of the original scheme presented previously.

Quantitative imaging with electrical impedance spectroscopy.

Electrical impedance spectroscopy (EIS) is a noninvasive modality that can be used to determine the electrical admittivity inside a body given a discrete set of current/voltage measurements made on the surface. Of particular interest is the use of EIS in the diagnosis of breast cancer, as the admittivity spectra of malignant and benign tumors differ significantly. Due to the fact that x-ray mammography is the current standard method of breast imaging to detect tumors, it is natural to see if we can use the admittivity distribution along with the mammogram image to improve the diagnosis, with the hopes that the specificity of these two methods combined will be greatly improved from using the mammogram image on its own. EIS is a highly ill-posed inverse problem, but regularization, in the form of structural prior information from the mammogram image as well as modeling error, allows for the problem to be solved for in a computationally efficient manner with improved results. To interpret the solution from the EIS inverse problem, a classification scheme is added, providing a quantitative image which maps out the tissue classification of the inside of the breast. The computational methods for solving the EIS inverse problem and the classification scheme are discussed and computed examples are presented to demonstrate the high simulated sensitivity and specificity of the method.