Optimizing MR Scan Design for Model-Based ${T}_{1}$ , ${T}_{2}$ Estimation From Steady-State Sequences

Neuroimaging of Parkinson's disease by quantitative susceptibility mapping

Parkinson's disease (PD) is a common neurodegenerative disease, and apart from a few rare genetic causes, its pathogenesis remains largely unclear. Recent scientific interest has been captured by the involvement of iron biochemistry and the disruption of iron homeostasis, particularly within the brain regions specifically affected in PD. The advent of Quantitative Susceptibility Mapping (QSM) has enabled non-invasive quantification of brain iron in vivo by MRI, which has contributed to the understanding of iron-associated pathogenesis and has the potential for the development of iron-based biomarkers in PD. This review elucidates the biochemical underpinnings of brain iron accumulation, details advancements in iron-sensitive MRI technologies, and discusses the role of QSM as a biomarker of iron deposition in PD. Despite considerable progress, several challenges impede its clinical application after a decade of QSM studies. The initiation of multi-site research is warranted for developing robust, interpretable, and disease-specific biomarkers for monitoring PD disease progression.

Cramér-Rao Bound Optimized Subspace Reconstruction in Quantitative MRI

We extend the traditional framework for estimating subspace bases that maximize the preserved signal energy to additionally preserve the Cramér-Rao bound (CRB) of the biophysical parameters and, ultimately, improve accuracy and precision in the quantitative maps. To this end, we introduce an approximate compressed CRB based on orthogonalized versions of the signal's derivatives with respect to the model parameters. This approximation permits singular value decomposition (SVD)-based minimization of both the CRB and signal losses during compression. Compared to the traditional SVD approach, the proposed method better preserves the CRB across all biophysical parameters with negligible cost to the preserved signal energy, leading to reduced bias and variance of the parameter estimates in simulation. In vivo, improved accuracy and precision are observed in two quantitative neuroimaging applications, permitting the use of smaller basis sizes in subspace reconstruction and offering significant computational savings.

Optimal flip angles for in vivo liver 3D T1 mapping and B1+ mapping at 3T

PURPOSE

The spoiled gradient recalled echo (SPGR) sequence with variable flip angles (FAs) enables whole liverT 1 $$ {T}_1 $$ mapping at high spatial resolutions but is strongly affected byB 1 + $$ {B}_1^{+} $$ inhomogeneities. The aim of this work was to study how the precision of acquiredT 1 $$ {T}_1 $$ maps is affected by theT 1 $$ {T}_1 $$ andB 1 + $$ {B}_1^{+} $$ ranges observed in the liver at 3T, as well as how noise propagates from the acquired signals into the resultingT 1 $$ {T}_1 $$ map.

THEORY

TheT 1 $$ {T}_1 $$ variance was estimated through the Fisher information matrix with a total noise variance including, for the first time, theB 1 + $$ {B}_1^{+} $$ map noise as well as contributions from the SPGR noise.

METHODS

Simulations were used to find the optimal FAs for both theB 1 + $$ {B}_1^{+} $$ mapping andT 1 $$ {T}_1 $$ mapping. The simulations results were validated in 10 volunteers.

RESULTS

Four optimized SPGR FAs of 2°, 2°, 15°, and 15° (TR = 4.1 ms) andB 1 + $$ {B}_1^{+} $$ map FAs of 65° and 130° achieved aT 1 $$ {T}_1 $$ coefficient of variation of 6.2 ± 1.7% across 10 volunteers and validated our theoretical model. Four optimal FAs outperformed five uniformly spaced FAs, saving the patient one breath-hold. For the liverB 1 + $$ {B}_1^{+} $$ andT 1 $$ {T}_1 $$ parameter space at 3T, a higher return inT 1 $$ {T}_1 $$ precision was obtained by investing FAs in the SPGR acquisition rather than in theB 1 + $$ {B}_1^{+} $$ map.

CONCLUSION

A novel framework was developed and validated to calculate the SPGRT 1 $$ {T}_1 $$ variance. This framework efficiently identifies optimal FA values and determines the total number of SPGR andB 1 + $$ {B}_1^{+} $$ measurements needed to achieve a desiredT 1 $$ {T}_1 $$ precision.

Optimization of spin-lock times in T1ρ mapping of knee cartilage: Cramér-Rao bounds versus matched sampling-fitting

PURPOSE

To compare different optimization approaches for choosing the spin-lock times (TSLs), in spin-lattice relaxation time in the rotating frame (T_1ρ ) mapping.

METHODS

Optimization criteria for TSLs based on Cramér-Rao lower bounds (CRLB) are compared with matched sampling-fitting (MSF) approaches for T_1ρ mapping on synthetic data, model phantoms, and knee cartilage. The MSF approaches are optimized using robust methods for noisy cost functions. The MSF approaches assume that optimal TSLs depend on the chosen fitting method. An iterative non-linear least squares (NLS) and artificial neural networks (ANN) are tested as two possible T_1ρ fitting methods for MSF approaches.

RESULTS

All optimized criteria were better than non-optimized ones. However, we observe that a modified CRLB and an MSF based on the mean of the normalized absolute error (MNAE) were more robust optimization approaches, performing well in all tested cases. The optimized TSLs obtained the best performance with synthetic data (3.5-8.0% error), model phantoms (1.5-2.8% error), and healthy volunteers (7.7-21.1% error), showing stable and improved quality results, comparing to non-optimized approaches (4.2-13.3% error on synthetic data, 2.1-6.2% error on model phantoms, 9.8-27.8% error on healthy volunteers).

CONCLUSION

A modified CRLB and the MSF based on MNAE are robust optimization approaches for choosing TSLs in T_1ρ mapping. All optimized criteria allowed good results even using rapid scans with two TSLs when a complex-valued fitting is done with iterative NLS or ANN.

Efficiency analysis for quantitative MRI of T1 and T2 relaxometry methods

This study presents a comparison of quantitative MRI methods based on an efficiency metric that quantifies their intrinsic ability to extract information about tissue parameters. Under a regime of unbiased parameter estimates, an intrinsic efficiency metricηwas derived for fully-sampled experiments which can be used to both optimize and compare sequences. Here we optimize and compare several steady-state and transient gradient-echo based qMRI methods, such as magnetic resonance fingerprinting (MRF), for jointT1andT2mapping. The impact of undersampling was also evaluated, assuming incoherent aliasing that is treated as noise by parameter estimation.In vivovalidation of the efficiency metric was also performed. Transient methods such as MRF can be up to 3.5 times more efficient than steady-state methods, when spatial undersampling is ignored. If incoherent aliasing is treated as noise during least-squares parameter estimation, the efficiency is reduced in proportion to the SNR of the data, with reduction factors of 5 often seen for practical SNR levels.In vivovalidation showed a very good agreement between the theoretical and experimentally predicted efficiency. This work presents and validates an efficiency metric to optimize and compare the performance of qMRI methods. Transient methods were found to be intrinsically more efficient than steady-state methods, however the effect of spatial undersampling can significantly erode this advantage.

An information theory model for optimizing quantitative magnetic resonance imaging acquisitions

Acquisition parameter selection is currently performed empirically for many quantitative MRI (qMRI) acquisitions. Tuning parameters for different scan times, tissues, and resolutions requires some amount of trial and error. There is an opportunity to quantitatively optimize these acquisition parameters in order to minimize variability of quantitative maps and post-processing techniques such as synthetic image generation. The objective of this work is to introduce and evaluate a quantitative method for selecting parameters that minimize image variability. An information theory framework was developed for this purpose and applied to a 3D-quantification using an interleaved Look-Locker acquisition sequence with T2 preparation pulse (3D-QALAS) signal model for qMRI. In this framework, mutual information is used to measure the information gained by a measurement as a function of acquisition parameters, quantifying the information content of potential acquisitions and allowing informed parameter selection. The information theory framework was tested on artificial data generated from a representative mathematical phantom, measurements acquired on a qMRI multiparametric imaging standard phantom, and in vivo measurements in a human brain. The phantom measurements showed that higher mutual information calculated by the model correlated with smaller coefficient of variation in the reconstructed parametric maps, and in vivo measurements demonstrated that information-based calibration of acquisition parameters resulted in a decrease in parametric map variability consistent with model predictions.

Myelin water fraction estimation using small-tip fast recovery MRI

PURPOSE

To demonstrate the feasibility of an optimized set of small-tip fast recovery (STFR) MRI scans for rapidly estimating myelin water fraction (MWF) in the brain.

METHODS

We optimized a set of STFR scans to minimize the Cramér-Rao Lower Bound of MWF estimates. We evaluated the RMSE of MWF estimates from the optimized scans in simulation. We compared STFR-based MWF estimates (both modeling exchange and not modeling exchange) to multi-echo spin echo (MESE)-based estimates. We used the optimized scans to acquire in vivo data from which a MWF map was estimated. We computed the STFR-based MWF estimates using PERK, a recently developed kernel regression technique, and the MESE-based MWF estimates using both regularized non-negative least squares (NNLS) and PERK.

RESULTS

In simulation, the optimized STFR scans led to estimates of MWF with low RMSE across a range of tissue parameters and across white matter and gray matter. The STFR-based MWF estimates that modeled exchange compared well to MESE-based MWF estimates in simulation. When the optimized scans were tested in vivo, the MWF map that was estimated using a 3-compartment model with exchange was closer to the MESE-based MWF map.

CONCLUSIONS

The optimized STFR scans appear to be well suited for estimating MWF in simulation and in vivo when we model exchange in training. In this case, the STFR-based MWF estimates are close to the MESE-based estimates.

Computational MRI with Physics-based Constraints: Application to Multi-contrast and Quantitative Imaging

Compressed sensing takes advantage of low-dimensional signal structure to reduce sampling requirements far below the Nyquist rate. In magnetic resonance imaging (MRI), this often takes the form of sparsity through wavelet transform, finite differences, and low rank extensions. Though powerful, these image priors are phenomenological in nature and do not account for the mechanism behind the image formation. On the other hand, MRI signal dynamics are governed by physical laws, which can be explicitly modeled and used as priors for reconstruction. These explicit and implicit signal priors can be synergistically combined in an inverse problem framework to recover sharp, multi-contrast images from highly accelerated scans. Furthermore, the physics-based constraints provide a recipe for recovering quantitative, bio-physical parameters from the data. This article introduces physics-based modeling constraints in MRI and shows how they can be used in conjunction with compressed sensing for image reconstruction and quantitative imaging. We describe model-based quantitative MRI, as well as its linear subspace approximation. We also discuss approaches to selecting user-controllable scan parameters given knowledge of the physical model. We present several MRI applications that take advantage of this framework for the purpose of multi-contrast imaging and quantitative 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.

NOVIFAST: A Fast Algorithm for Accurate and Precise VFA MRI Mapping

In quantitative magnetic resonance mapping, the variable flip angle (VFA) steady state spoiled gradient recalled echo (SPGR) imaging technique is popular as it provides a series of high resolution weighted images in a clinically feasible time. Fast, linear methods that estimate maps from these weighted images have been proposed, such as DESPOT1 and iterative re-weighted linear least squares. More accurate, non-linear least squares (NLLS) estimators are in play, but these are generally much slower and require careful initialization. In this paper, we present NOVIFAST, a novel NLLS-based algorithm specifically tailored to VFA SPGR mapping. By exploiting the particular structure of the SPGR model, a computationally efficient, yet accurate and precise map estimator is derived. Simulation and in vivo human brain experiments demonstrate a twenty-fold speed gain of NOVIFAST compared with conventional gradient-based NLLS estimators while maintaining a high precision and accuracy. Moreover, NOVIFAST is eight times faster than the efficient implementations of the variable projection (VARPRO) method. Furthermore, NOVIFAST is shown to be robust against initialization.

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.

Dictionary-Free MRI PERK: Parameter Estimation via Regression with Kernels

This paper introduces a fast, general method for dictionary-free parameter estimation in quantitative magnetic resonance imaging (QMRI) parameter estimation via regression with kernels (PERK). PERK first uses prior distributions and the nonlinear MR signal model to simulate many parameter-measurement pairs. Inspired by machine learning, PERK then takes these parameter-measurement pairs as labeled training points and learns from them a nonlinear regression function using kernel functions and convex optimization. PERK admits a simple implementation as per-voxel nonlinear lifting of MRI measurements followed by linear minimum mean-squared error regression. We demonstrate PERK for $ {\textit {T}_{1}}, {\textit {T}_{2}}$ estimation, a well-studied application where it is simple to compare PERK estimates against dictionary-based grid search estimates and iterative optimization estimates. Numerical simulations as well as single-slice phantom and in vivo experiments demonstrate that PERK and other tested methods produce comparable $ {\textit {T}_{1}}, {\textit {T}_{2}}$ estimates in white and gray matter, but PERK is consistently at least $140\times $ faster. This acceleration factor may increase by several orders of magnitude for full-volume QMRI estimation problems involving more latent parameters per voxel.