Low rank alternating direction method of multipliers reconstruction for MR fingerprinting

Cardiac MR Fingerprinting: Overview, Technical Developments, and Applications

Cardiovascular magnetic resonance (CMR) is an established imaging modality with proven utility in assessing cardiovascular diseases. The ability of CMR to characterize myocardial tissue using T1- and T2-weighted imaging, parametric mapping, and late gadolinium enhancement has allowed for the non-invasive identification of specific pathologies not previously possible with modalities like echocardiography. However, CMR examinations are lengthy and technically complex, requiring multiple pulse sequences and different anatomical planes to comprehensively assess myocardial structure, function, and tissue composition. To increase the overall impact of this modality, there is a need to simplify and shorten CMR exams to improve access and efficiency, while also providing reproducible quantitative measurements. Multiparametric MRI techniques that measure multiple tissue properties offer one potential solution to this problem. This review provides an in-depth look at one such multiparametric approach, cardiac magnetic resonance fingerprinting (MRF). The article is structured as follows. First, a brief review of single-parametric and (non-Fingerprinting) multiparametric CMR mapping techniques is presented. Second, a general overview of cardiac MRF is provided covering pulse sequence implementation, dictionary generation, fast k-space sampling methods, and pattern recognition. Third, recent technical advances in cardiac MRF are covered spanning a variety of topics, including simultaneous multislice and 3D sampling, motion correction algorithms, cine MRF, synthetic multicontrast imaging, extensions to measure additional clinically important tissue properties (proton density fat fraction, T2*, and T1ρ), and deep learning methods for image reconstruction and parameter estimation. The last section will discuss potential clinical applications, concluding with a perspective on how multiparametric techniques like MRF may enable streamlined CMR protocols. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 1.

Deep plug-and-play MRI reconstruction based on multiple complementary priors

Magnetic resonance imaging (MRI) is widely used in clinical diagnosis as a safe, non-invasive, high-resolution medical imaging technology, but long scanning time has been a major challenge for this technology. The undersampling reconstruction method has become an important technical means to accelerate MRI by reducing the data sampling rate while maintaining high-quality imaging. However, traditional undersampling reconstruction techniques such as compressed sensing mainly rely on relatively single sparse or low-rank prior information to reconstruct the image, which has limitations in capturing the comprehensive features of images, resulting in the insufficient performance of the reconstructed image in terms of details and key information. In this paper, we propose a deep plug-and-play multiple complementary priors MRI reconstruction model, which combines traditional low-rank matrix recovery model methods and deep learning methods, and integrates global, local and nonlocal priors to improve reconstruction quality. Specifically, we capture the global features of the image through the matrix nuclear norm, and use the deep convolutional neural network denoiser Swin-Conv-UNet (SCUNet) and block-matching and 3-D filtering (BM3D) algorithm to preserve the local details and structural texture of the image, respectively. In addition, we utilize an efficient half-quadratic splitting (HQS) algorithm to solve the proposed model. The experimental results show that our proposed method has better reconstruction ability than the existing popular methods in terms of visual effects and numerical results.

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.

Bias-reduced neural networks for parameter estimation in quantitative MRI

PURPOSE

To develop neural network (NN)-based quantitative MRI parameter estimators with minimal bias and a variance close to the Cramér-Rao bound.

THEORY AND METHODS

We generalize the mean squared error loss to control the bias and variance of the NN's estimates, which involves averaging over multiple noise realizations of the same measurements during training. Bias and variance properties of the resulting NNs are studied for two neuroimaging applications.

RESULTS

In simulations, the proposed strategy reduces the estimates' bias throughout parameter space and achieves a variance close to the Cramér-Rao bound. In vivo, we observe good concordance between parameter maps estimated with the proposed NNs and traditional estimators, such as nonlinear least-squares fitting, while state-of-the-art NNs show larger deviations.

CONCLUSION

The proposed NNs have greatly reduced bias compared to those trained using the mean squared error and offer significantly improved computational efficiency over traditional estimators with comparable or better accuracy.

Manifold Regularizer for High-Resolution fMRI Joint Reconstruction and Dynamic Quantification

Oscillating Steady-State Imaging (OSSI) is a recently developed fMRI acquisition method that can provide 2 to 3 times higher SNR than standard fMRI approaches. However, because the OSSI signal exhibits a nonlinear oscillation pattern, one must acquire and combine nc (e.g., 10) OSSI images to get an image that is free of oscillation for fMRI, and fully sampled acquisitions would compromise temporal resolution. To improve temporal resolution and accurately model the nonlinearity of OSSI signals, instead of using subspace models that are not well suited for the data, we build the MR physics for OSSI signal generation as a regularizer for the undersampled reconstruction. Our proposed physics-based manifold model turns the disadvantages of OSSI acquisition into advantages and enables joint reconstruction and quantification. OSSI manifold model (OSSIMM) outperforms subspace models and reconstructs high-resolution fMRI images with a factor of 12 acceleration and without spatial or temporal smoothing. Furthermore, OSSIMM can dynamically quantify important physics parameters, including R2∗ maps, with a temporal resolution of 150 ms.

Sensitivity of unconstrained quantitative magnetization transfer MRI to Amyloid burden in preclinical Alzheimer's disease

Magnetization transfer MRI is sensitive to semi-solid macromolecules, including amyloid beta, and has previously been used to discriminate Alzheimer's disease (AD) patients from controls. Here, we fit an unconstrained 2-pool quantitative MT (qMT) model, i.e., without constraints on the longitudinal relaxation rateR 1 s of semi-solids, and investigate the sensitivity of the estimated parameters to amyloid accumulation in preclinical subjects. We scanned 15 cognitively normal volunteers, of which 9 were amyloid positive by [18F]Florbetaben PET. A 12 min hybrid-state qMT scan with an effective resolution of 1.24 mm isotropic and whole-brain coverage was acquired to estimate the unconstrained 2-pool qMT parameters. Group comparisons and correlations with Florbetaben PET standardized uptake value ratios were analyzed at the lobar level. We find that the exchange rate and semi-solid pool'sR 1 s were sensitive to the amyloid concentration, while morphometric measures of cortical thickness derived from structural MRI were not. Changes in the exchange rate are consistent with previous reports in clinical AD, while changes inR 1 s have not been reported previously as its value is typically constrained in the literature. Our results demonstrate that qMT MRI may be a promising surrogate marker of amyloid beta without the need for contrast agents or radiotracers.

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.

Fourier Domain Robust Denoising Decomposition and Adaptive Patch MRI Reconstruction

The sparsity of the Fourier transform domain has been applied to magnetic resonance imaging (MRI) reconstruction in k -space. Although unsupervised adaptive patch optimization methods have shown promise compared to data-driven-based supervised methods, the following challenges exist in MRI reconstruction: 1) in previous k -space MRI reconstruction tasks, MRI with noise interference in the acquisition process is rarely considered. 2) Differences in transform domains should be resolved to achieve the high-quality reconstruction of low undersampled MRI data. 3) Robust patch dictionary learning problems are usually nonconvex and NP-hard, and alternate minimization methods are often computationally expensive. In this article, we propose a method for Fourier domain robust denoising decomposition and adaptive patch MRI reconstruction (DDAPR). DDAPR is a two-step optimization method for MRI reconstruction in the presence of noise and low undersampled data. It includes the low-rank and sparse denoising reconstruction model (LSDRM) and the robust dictionary learning reconstruction model (RDLRM). In the first step, we propose LSDRM for different domains. For the optimization solution, the proximal gradient method is used to optimize LSDRM by singular value decomposition and soft threshold algorithms. In the second step, we propose RDLRM, which is an effective adaptive patch method by introducing a low-rank and sparse penalty adaptive patch dictionary and using a sparse rank-one matrix to approximate the undersampled data. Then, the block coordinate descent (BCD) method is used to optimize the variables. The BCD optimization process involves valid closed-form solutions. Extensive numerical experiments show that the proposed method has a better performance than previous methods in image reconstruction based on compressed sensing or deep learning.

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.

Bias-Reduced Neural Networks for Parameter Estimation in Quantitative MRI

Purpose

To develop neural network (NN)-based quantitative MRI parameter estimators with minimal bias and a variance close to the Cramér-Rao bound.

Theory and Methods

We generalize the mean squared error loss to control the bias and variance of the NN's estimates, which involves averaging over multiple noise realizations of the same measurements during training. Bias and variance properties of the resulting NNs are studied for two neuroimaging applications.

Results

In simulations, the proposed strategy reduces the estimates' bias throughout parameter space and achieves a variance close to the Cramér-Rao bound. In vivo, we observe good concordance between parameter maps estimated with the proposed NNs and traditional estimators, such as non-linear least-squares fitting, while state-of-the-art NNs show larger deviations.

Conclusion

The proposed NNs have greatly reduced bias compared to those trained using the mean squared error and offer significantly improved computational efficiency over traditional estimators with comparable or better accuracy.

Rapid quantitative magnetization transfer imaging: Utilizing the hybrid state and the generalized Bloch model

PURPOSE

To explore efficient encoding schemes for quantitative magnetization transfer (qMT) imaging with few constraints on model parameters.

THEORY AND METHODS

We combine two recently proposed models in a Bloch-McConnell equation: the dynamics of the free spin pool are confined to the hybrid state, and the dynamics of the semi-solid spin pool are described by the generalized Bloch model. We numerically optimize the flip angles and durations of a train of radio frequency pulses to enhance the encoding of three qMT parameters while accounting for all eight parameters of the two-pool model. We sparsely sample each time frame along this spin dynamics with a three-dimensional radial koosh-ball trajectory, reconstruct the data with subspace modeling, and fit the qMT model with a neural network for computational efficiency.

RESULTS

We extracted qMT parameter maps of the whole brain with an effective resolution of 1.24 mm from a 12.6-min scan. In lesions of multiple sclerosis subjects, we observe a decreased size of the semi-solid spin pool and longer relaxation times, consistent with previous reports.

CONCLUSION

The encoding power of the hybrid state, combined with regularized image reconstruction, and the accuracy of the generalized Bloch model provide an excellent basis for efficient quantitative magnetization transfer imaging with few constraints on model parameters.

Unconstrained quantitative magnetization transfer imaging: disentangling T1 of the free and semi-solid spin pools

Since the inception of magnetization transfer (MT) imaging, it has been widely assumed that Henkelman's two spin pools have similar longitudinal relaxation times, which motivated many researchers to constrain them to each other. However, several recent publications reported a T 1 s of the semi-solid spin pool that is much shorter than T 1 f of the free pool. While these studies tailored experiments for robust proofs-of-concept, we here aim to quantify the disentangled relaxation processes on a voxel-by-voxel basis in a clinical imaging setting, i.e., with an effective resolution of 1.24mm isotropic and full brain coverage in 12min. To this end, we optimized a hybrid-state pulse sequence for mapping the parameters of an unconstrained MT model. We scanned four people with relapsing-remitting multiple sclerosis (MS) and four healthy controls with this pulse sequence and estimated T 1 f ≈ 1.84 s and T 1 s ≈ 0.34 s in healthy white matter. Our results confirm the reports that T 1 s ≪ T 1 f and we argue that this finding identifies MT as an inherent driver of longitudinal relaxation in brain tissue. Moreover, we estimated a fractional size of the semi-solid spin pool of m 0 s ≈ 0.212 , which is larger than previously assumed. An analysis of T 1 f in normal-appearing white matter revealed statistically significant differences between individuals with MS and controls.

Minimization of eddy current artifacts in sequences with periodic dynamics

PURPOSE

To minimize eddy current artifacts in periodic pulse sequences with balanced gradient moments as, for example, used for quantitative MRI.

THEORY AND METHODS

Eddy current artifacts in balanced sequences result from large jumps in k-space. In quantitative MRI, one often samples some spin dynamics repeatedly while acquiring different parts of k-space. We swap individual k-space lines between different repetitions in order to minimize jumps in temporal succession without changing the overall trajectory. This reordering can be formulated as a traveling salesman problem and we tackle the discrete optimization with a simulated annealing algorithm.

RESULTS

Compared to the default ordering, we observe a substantial reduction of artifacts in the reconstructed images and the derived quantitative parameter maps. Comparing two variants of our algorithm, one that resembles the pairing approach originally proposed by Bieri et al., and one that minimizes all k-space jumps equally, we observe slightly lower artifact levels in the latter.

CONCLUSION

The proposed reordering scheme effectively reduces eddy current artifacts in sequences with balanced gradient moments. In contrast to previous approaches, we capitalize on the periodicity of the sampled signal dynamics, enabling both efficient k-space sampling and minimizing artifacts caused by eddy currents.

Emerging Trends in Magnetic Resonance Fingerprinting for Quantitative Biomedical Imaging Applications: A Review

Magnetic resonance imaging (MRI) stands as a vital medical imaging technique, renowned for its ability to offer high-resolution images of the human body with remarkable soft-tissue contrast. This enables healthcare professionals to gain valuable insights into various aspects of the human body, including morphology, structural integrity, and physiological processes. Quantitative imaging provides compositional measurements of the human body, but, currently, either it takes a long scan time or is limited to low spatial resolutions. Undersampled k-space data acquisitions have significantly helped to reduce MRI scan time, while compressed sensing (CS) and deep learning (DL) reconstructions have mitigated the associated undersampling artifacts. Alternatively, magnetic resonance fingerprinting (MRF) provides an efficient and versatile framework to acquire and quantify multiple tissue properties simultaneously from a single fast MRI scan. The MRF framework involves four key aspects: (1) pulse sequence design; (2) rapid (undersampled) data acquisition; (3) encoding of tissue properties in MR signal evolutions or fingerprints; and (4) simultaneous recovery of multiple quantitative spatial maps. This paper provides an extensive literature review of the MRF framework, addressing the trends associated with these four key aspects. There are specific challenges in MRF for all ranges of magnetic field strengths and all body parts, which can present opportunities for further investigation. We aim to review the best practices in each key aspect of MRF, as well as for different applications, such as cardiac, brain, and musculoskeletal imaging, among others. A comprehensive review of these applications will enable us to assess future trends and their implications for the translation of MRF into these biomedical imaging applications.

High-resolution myelin-water fraction and quantitative relaxation mapping using 3D ViSTa-MR fingerprinting

Purpose

This study aims to develop a high-resolution whole-brain multi-parametric quantitative MRI approach for simultaneous mapping of myelin-water fraction (MWF), T1, T2, and proton-density (PD), all within a clinically feasible scan time.

Methods

We developed 3D ViSTa-MRF, which combined Visualization of Short Transverse relaxation time component (ViSTa) technique with MR Fingerprinting (MRF), to achieve high-fidelity whole-brain MWF and T1/T2/PD mapping on a clinical 3T scanner. To achieve fast acquisition and memory-efficient reconstruction, the ViSTa-MRF sequence leverages an optimized 3D tiny-golden-angle-shuffling spiral-projection acquisition and joint spatial-temporal subspace reconstruction with optimized preconditioning algorithm. With the proposed ViSTa-MRF approach, high-fidelity direct MWF mapping was achieved without a need for multi-compartment fitting that could introduce bias and/or noise from additional assumptions or priors.

Results

The in-vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide fast multi-parametric mapping with high SNR and good quality. The in-vivo results of 1mm- and 0.66mm-iso datasets indicate that the MWF values measured by the proposed method are consistent with standard ViSTa results that are 30x slower with lower SNR. Furthermore, we applied the proposed method to enable 5-minute whole-brain 1mm-iso assessment of MWF and T1/T2/PD mappings for infant brain development and for post-mortem brain samples.

Conclusions

In this work, we have developed a 3D ViSTa-MRF technique that enables the acquisition of whole-brain MWF, quantitative T1, T2, and PD maps at 1mm and 0.66mm isotropic resolution in 5 and 15 minutes, respectively. This advancement allows for quantitative investigations of myelination changes in the brain.

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.

Global deep learning optimization of chemical exchange saturation transfer magnetic resonance fingerprinting acquisition schedule

Chemical exchange saturation transfer (CEST) MRI is a promising molecular imaging technique but suffers from long scan times and complicated processing. CEST was recently combined with magnetic resonance fingerprinting (MRF) to address these shortcomings. However, the CEST-MRF signal depends on multiple acquisition and tissue parameters so selecting an optimal acquisition schedule is challenging. In this work, we propose a novel dual-network deep learning framework to optimize the CEST-MRF acquisition schedule. The quality of the optimized schedule was assessed in a digital brain phantom and compared with alternate deep learning optimization approaches. The effect of schedule length on the reconstruction error was also investigated. A healthy subject was scanned with optimized and random schedules and with a conventional CEST sequence for comparison. The optimized schedule was also tested in a subject with metastatic renal cell carcinoma. Reproducibility was assessed via test-retest experiments and the concordance correlation coefficient calculated for white matter (WM) and grey matter (GM). The optimized schedule was 12% shorter but yielded equal or lower normalized root mean square error for all parameters. The proposed optimization also provided a lower error compared with alternate methodologies. Longer schedules generally yielded lower error. In vivo maps obtained with the optimized schedule showed reduced noise and improved delineation of GM and WM. CEST curves synthesized from the optimized parameters were highly correlated (r = 0.99) with measured conventional CEST. The mean concordance correlation coefficient in WM/GM for all tissue parameters was 0.990/0.978 for the optimized schedule but only 0.979/0.975 for the random schedule. The proposed schedule optimization is widely applicable to MRF pulse sequences and provides accurate and reproducible tissue maps with reduced noise at a shorter scan time than a randomly generated schedule.

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.

Quantitative longitudinal mapping of radiation-treated prostate cancer using MR fingerprinting with radial acquisition and subspace reconstruction

MR fingerprinting (MRF) enables fast multiparametric quantitative imaging with a single acquisition and has been shown to improve diagnosis of prostate cancer. However, most prostate MRF studies were performed with spiral acquisitions that are sensitive to B0 inhomogeneities and consequent blurring. In this work, a radial MRF acquisition with a novel subspace reconstruction technique was developed to enable fast T1/T2 mapping in the prostate in under 4 min. The subspace reconstruction exploits the extensive temporal correlations in the MRF dictionary to pre-compute a low dimensional space for the solution and thus reduce the number of radial spokes to accelerate the acquisition. Iterative reconstruction with the subspace model and additional regularization of the signal representation in the subspace is performed to minimize the number of spokes and maintain matching quality and SNR. Reconstruction accuracy was assessed using the ISMRM NIST phantom. In-vivo validation was performed on two healthy subjects and two prostate cancer patients undergoing radiation therapy. The longitudinal repeatability was quantified using the concordance correlation coefficient (CCC) in one of the healthy subjects by repeated scans over 1 year. One prostate cancer patient was scanned at three time points, before initiating therapy and following brachytherapy and external beam radiation. Changes in the T1/T2 maps obtained with the proposed method were quantified. The prostate, peripheral and transitional zones, and visible dominant lesion were delineated for each study, and the statistics and distribution of the quantitative mapping values were analyzed. Significant image quality improvements compared with standard reconstruction methods were obtained with the proposed subspace reconstruction method. A notable decrease in the spread of the T1/T2 values without biasing the estimated mean values was observed with the subspace reconstruction and agreed with reported literature values. The subspace reconstruction enabled visualization of small differences in T1/T2 values in the tumor region within the peripheral zone. Longitudinal imaging of a volunteer subject yielded CCC of 0.89 for MRF T1, and 0.81 for MRF T2 in the prostate gland. Longitudinal imaging of the prostate patient confirmed the feasibility of capturing radiation treatment related changes. This work is a proof-of-concept for a high resolution and fast quantitative mapping using golden-angle radial MRF combined with a subspace reconstruction technique for longitudinal treatment response assessment in subjects undergoing radiation treatment.

RIFLE: Imputation and Robust Inference from Low Order Marginals

The ubiquity of missing values in real-world datasets poses a challenge for statistical inference and can prevent similar datasets from being analyzed in the same study, precluding many existing datasets from being used for new analyses. While an extensive collection of packages and algorithms have been developed for data imputation, the overwhelming majority perform poorly if there are many missing values and low sample sizes, which are unfortunately common characteristics in empirical data. Such low-accuracy estimations adversely affect the performance of downstream statistical models. We develop a statistical inference framework for regression and classification in the presence of missing data without imputation. Our framework, RIFLE (Robust InFerence via Low-order moment Estimations), estimates low-order moments of the underlying data distribution with corresponding confidence intervals to learn a distributionally robust model. We specialize our framework to linear regression and normal discriminant analysis, and we provide convergence and performance guarantees. This framework can also be adapted to impute missing data. In numerical experiments, we compare RIFLE to several state-of-the-art approaches (including MICE, Amelia, MissForest, KNN-imputer, MIDA, and Mean Imputer) for imputation and inference in the presence of missing values. Our experiments demonstrate that RIFLE outperforms other benchmark algorithms when the percentage of missing values is high and/or when the number of data points is relatively small. RIFLE is publicly available at https://github.com/optimization-for-data-driven-science/RIFLE.

Multiparametric MRI: From Simultaneous Rapid Acquisition Methods and Analysis Techniques Using Scoring, Machine Learning, Radiomics, and Deep Learning to the Generation of Novel Metrics

ABSTRACT

With the recent advancements in rapid imaging methods, higher numbers of contrasts and quantitative parameters can be acquired in less and less time. Some acquisition models simultaneously obtain multiparametric images and quantitative maps to reduce scan times and avoid potential issues associated with the registration of different images. Multiparametric magnetic resonance imaging (MRI) has the potential to provide complementary information on a target lesion and thus overcome the limitations of individual techniques. In this review, we introduce methods to acquire multiparametric MRI data in a clinically feasible scan time with a particular focus on simultaneous acquisition techniques, and we discuss how multiparametric MRI data can be analyzed as a whole rather than each parameter separately. Such data analysis approaches include clinical scoring systems, machine learning, radiomics, and deep learning. Other techniques combine multiple images to create new quantitative maps associated with meaningful aspects of human biology. They include the magnetic resonance g-ratio, the inner to the outer diameter of a nerve fiber, and the aerobic glycolytic index, which captures the metabolic status of tumor tissues.

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

OBJECTIVE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

KomaMRI.jl: An open-source framework for general MRI simulations with GPU acceleration

PURPOSE

To develop an open-source, high-performance, easy-to-use, extensible, cross-platform, and general MRI simulation framework (Koma).

METHODS

Koma was developed using the Julia programming language. Like other MRI simulators, it solves the Bloch equations with CPU and GPU parallelization. The inputs are the scanner parameters, the phantom, and the pulse sequence that is Pulseq-compatible. The raw data is stored in the ISMRMRD format. For the reconstruction, MRIReco.jl is used. A graphical user interface utilizing web technologies was also designed. Two types of experiments were performed: one to compare the quality of the results and the execution speed, and the second to compare its usability. Finally, the use of Koma in quantitative imaging was demonstrated by simulating Magnetic Resonance Fingerprinting (MRF) acquisitions.

RESULTS

Koma was compared to two well-known open-source MRI simulators, JEMRIS and MRiLab. Highly accurate results (with mean absolute differences below 0.1% compared to JEMRIS) and better GPU performance than MRiLab were demonstrated. In an experiment with students, Koma was proved to be easy to use, eight times faster on personal computers than JEMRIS, and 65% of test subjects recommended it. The potential for designing acquisition and reconstruction techniques was also shown through the simulation of MRF acquisitions, with conclusions that agree with the literature.

CONCLUSIONS

Koma's speed and flexibility have the potential to make simulations more accessible for education and research. Koma is expected to be used for designing and testing novel pulse sequences before implementing them in the scanner with Pulseq files, and for creating synthetic data to train machine learning models.

Low-rank motion correction for accelerated free-breathing first-pass myocardial perfusion imaging

PURPOSE

Develop a novel approach for accelerated 2D free-breathing myocardial perfusion via low-rank motion-corrected (LRMC) reconstructions.

METHODS

Myocardial perfusion imaging requires high spatial and temporal resolution, despite scan time constraints. Here, we incorporate LRMC models into the reconstruction-encoding operator, together with high-dimensionality patch-based regularization, to produce high quality, motion-corrected myocardial perfusion series from free-breathing acquisitions. The proposed framework estimates beat-to-beat nonrigid respiratory (and any other incidental) motion and the dynamic contrast subspace from the actual acquired data, which are then incorporated into the proposed LRMC reconstruction. LRMC was compared with iterative SENSitivity Encoding (SENSE) (itSENSE) and low-rank plus sparse (LpS) reconstruction in 10 patients based on image-quality scoring and ranking by two clinical expert readers.

RESULTS

LRMC achieved significantly improved results relative to itSENSE and LpS in terms of image sharpness, temporal coefficient of variation, and expert reader evaluation. Left ventricle image sharpness was approximately 75%, 79%, and 86% for itSENSE, LpS and LRMC, respectively, indicating improved image sharpness for the proposed approach. Corresponding temporal coefficient of variation results were 23%, 11% and 7%, demonstrating improved temporal fidelity of the perfusion signal with the proposed LRMC. Corresponding clinical expert reader scores (1-5, from poor to excellent image quality) were 3.3, 3.9 and 4.9, demonstrating improved image quality with the proposed LRMC, in agreement with the automated metrics.

CONCLUSION

LRMC produces motion-corrected myocardial perfusion in free-breathing acquisitions with substantially improved image quality when compared with iterative SENSE and LpS reconstructions.

Mitigating undersampling errors in MR fingerprinting by sequence optimization

PURPOSE

To develop a method for MR Fingerprinting (MRF) sequence optimization that takes both the applied undersampling pattern and a realistic reference map into account.

METHODS

A predictive model for the undersampling error leveraging on perturbation theory was exploited to optimize the MRF flip angle sequence for improved robustness against undersampling artifacts. In this framework parameter maps from a previously acquired MRF scan were used as reference. Sequences were optimized for different sequence lengths, smoothness constraints and undersampling factors. Numerical simulations and in vivo measurements in eight healthy subjects were performed to assess the effect of the performed optimization. The optimized MRF sequences were compared to a conventionally shaped flip angle pattern and an optimized pattern based on the Cramér-Rao lower bound (CRB).

RESULTS

Numerical simulations and in vivo results demonstrate that the undersampling errors can be suppressed by flip angle optimization. Analysis of the in vivo results show that a sequence optimized for improved robustness against undersampling with a flip angle train of length 400 yielded significantly lower median absolute errors in T 1 : 5 . 6 % ± 2 . 9 % and T 2 : 7 . 9 % ± 2 . 3 % compared to the conventional ( T 1 : 8 . 0 % ± 1 . 9 % , T 2 : 14 . 5 % ± 2 . 6 % ) and CRB-based ( T 1 : 21 . 6 % ± 4 . 1 % , T 2 : 31 . 4 % ± 4 . 4 % ) sequences.

CONCLUSION

The proposed method is able to optimize the MRF flip angle pattern such that significant mitigation of the artifacts from strong k-space undersampling in MRF is achieved.

Quantitative MRI by nonlinear inversion of the Bloch equations

PURPOSE

Development of a generic model-based reconstruction framework for multiparametric quantitative MRI that can be used with data from different pulse sequences.

METHODS

Generic nonlinear model-based reconstruction for quantitative MRI estimates parametric maps directly from the acquired k-space by numerical optimization. This requires numerically accurate and efficient methods to solve the Bloch equations and their partial derivatives. In this work, we combine direct sensitivity analysis and pre-computed state-transition matrices into a generic framework for calibrationless model-based reconstruction that can be applied to different pulse sequences. As a proof-of-concept, the method is implemented and validated for quantitative T 1 $$ {\mathrm{T}}_1 $$ and T 2 $$ {\mathrm{T}}_2 $$ mapping with single-shot inversion-recovery (IR) FLASH and IR bSSFP sequences in simulations, phantoms, and the human brain.

RESULTS

The direct sensitivity analysis enables a highly accurate and numerically stable calculation of the derivatives. The state-transition matrices efficiently exploit repeating patterns in pulse sequences, speeding up the calculation by a factor of 10 for the examples considered in this work, while preserving the accuracy of native ordinary differential equations solvers. The generic model-based method reproduces quantitative results of previous model-based reconstructions based on the known analytical solutions for radial IR FLASH. For IR bSFFP it produces accurate T 1 $$ {\mathrm{T}}_1 $$ and T 2 $$ {\mathrm{T}}_2 $$ maps for the National Insitute of Standards and Technology (NIST) phantom in numerical simulations and experiments. Feasibility is also shown for human brain, although results are affected by magnetization transfer effects.

CONCLUSION

By developing efficient tools for numerical optimizations using the Bloch equations as forward model, this work enables generic model-based reconstruction for quantitative MRI.

Efficient performance analysis and optimization of transient-state sequences for multiparametric magnetic resonance imaging

In transient-state multiparametric MRI sequences such as Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT), MR fingerprinting, or hybrid-state imaging, the flip angle pattern of the RF excitation varies over the sequence. This gives considerable freedom to choose an optimal pattern of flip angles. For pragmatic reasons, most optimization methodologies choose for a single-voxel approach (i.e., without taking the spatial encoding scheme into account). Particularly in MR-STAT, the context of spatial encoding is important. In the current study, we present a methodology, called BLock Analysis of a K-space-domain Jacobian (BLAKJac), which is sufficiently fast to optimize a sequence in the context of a predetermined phase-encoding pattern. Based on MR-STAT acquisitions and reconstructions, we show that sequences optimized using BLAKJac are more reliable in terms of actually achieved precision than conventional single-voxel-optimized sequences. In addition, BLAKJac provides analytical tools that give insights into the performance of the sequence in a very limited computation time. Our experiments are based on MR-STAT, but the theory is equally valid for other transient-state multiparametric methods.

Free-running 3D whole-heart T1 and T2 mapping and cine MRI using low-rank reconstruction with non-rigid cardiac motion correction

PURPOSE

To introduce non-rigid cardiac motion correction into a novel free-running framework for the simultaneous acquisition of 3D whole-heart myocardial <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> maps and cine images, enabling a <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>∼</mml:mo></mml:mrow> <mml:annotation>$$ \sim $$</mml:annotation></mml:semantics> </mml:math> 3-min scan.

METHODS

Data were acquired using a free-running 3D golden-angle radial readout interleaved with inversion recovery and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> -preparation pulses. After correction for translational respiratory motion, non-rigid cardiac-motion-corrected reconstruction with dictionary-based low-rank compression and patch-based regularization enabled 3D whole-heart <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> mapping at any given cardiac phase as well as whole-heart cardiac cine imaging. The framework was validated and compared with established methods in 11 healthy subjects.

RESULTS

Good quality 3D <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> maps and cine images were reconstructed for all subjects. Septal <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> values using the proposed approach ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>1200</mml:mn> <mml:mo>±</mml:mo> <mml:mn>50</mml:mn></mml:mrow> <mml:annotation>$$ 1200\pm 50 $$</mml:annotation></mml:semantics> </mml:math> ms) were higher than those from a 2D MOLLI sequence ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>1063</mml:mn> <mml:mo>±</mml:mo> <mml:mn>33</mml:mn></mml:mrow> <mml:annotation>$$ 1063\pm 33 $$</mml:annotation></mml:semantics> </mml:math> ms), which is known to underestimate <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> , while <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> values from the proposed approach ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>51</mml:mn> <mml:mo>±</mml:mo> <mml:mn>4</mml:mn></mml:mrow> <mml:annotation>$$ 51\pm 4 $$</mml:annotation></mml:semantics> </mml:math> ms) were in good agreement with those from a 2D GraSE sequence ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>51</mml:mn> <mml:mo>±</mml:mo> <mml:mn>2</mml:mn></mml:mrow> <mml:annotation>$$ 51\pm 2 $$</mml:annotation></mml:semantics> </mml:math> ms).

CONCLUSION

The proposed technique provides 3D <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> maps and cine images with isotropic spatial resolution in a single <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>∼</mml:mo></mml:mrow> <mml:annotation>$$ \sim $$</mml:annotation></mml:semantics> </mml:math> 3.3-min scan.

Multicomponent MR fingerprinting reconstruction using joint-sparsity and low-rank constraints

PURPOSE

To develop an efficient algorithm for multicomponent MR fingerprinting (MC-MRF) reconstructions directly from highly undersampled data without making prior assumptions about tissue relaxation times and expected number of tissues.

METHODS

The proposed method reconstructs MC-MRF maps from highly undersampled data by iteratively applying a joint-sparsity constraint to the estimated tissue components. Intermediate component maps are obtained by a low-rank multicomponent alternating direction method of multipliers (MC-ADMM) including the non-negativity of tissue weights as an extra regularization term. Over iterations, the used dictionary compression is adjusted. The proposed method (k-SPIJN) is compared with a two-step approach in which image reconstruction and multicomponent estimations are performed sequentially and tested in numerical simulations and in vivo by applying different undersampling factors in eight healthy volunteers. In the latter case, fully sampled data serves as the reference.

RESULTS

The proposed method shows improved precision and accuracy in simulations compared with a state-of-art sequential approach. Obtained in vivo magnetization fraction maps for different tissue types show reduced systematic errors and reduced noise-like effects. Root mean square errors in estimated magnetization fraction maps significantly reduce from 13.0% <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>±</mml:mo></mml:mrow> <mml:annotation>$$ \pm $$</mml:annotation></mml:semantics> </mml:math> 5.8% with the conventional, two-step approach to 9.6% <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>±</mml:mo></mml:mrow> <mml:annotation>$$ \pm $$</mml:annotation></mml:semantics> </mml:math> 3.9% and 9.6% <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>±</mml:mo></mml:mrow> <mml:annotation>$$ \pm $$</mml:annotation></mml:semantics> </mml:math> 3.2% with the proposed MC-ADMM and k-SPIJN methods, respectively. Mean standard deviation in homogeneous white matter regions reduced significantly from 8.6% to 2.9% (two step vs. k-SPIJN).

CONCLUSION

The proposed MC-ADMM and k-SPIJN reconstruction methods estimate MC-MRF maps from highly undersampled data resulting in improved image quality compared with the existing method.

3D magnetic resonance fingerprinting for rapid simultaneous T1, T2, and T1ρ volumetric mapping of human articular cartilage at 3 T

Quantitative MRI can detect early biochemical changes in cartilage; however, the conventional techniques only measure one parameter (e.g., T1 , T2 , and T1ρ ) at a time while also being comparatively slow. We implemented a 3D magnetic resonance fingerprinting (3D-MRF) technique for simultaneous, volumetric mapping of T1 , T2 , and T1ρ in knee articular cartilage in under 9 min. It is evaluated on 11 healthy volunteers (mean age: 53 ± 9 years), five mild knee osteoarthritis (OA) patients (Kellgren-Lawrence (KL) score: 2, mean age: 60 ± 4 years), and the National Institute of Standards and Technology (NIST)/International Society for Magnetic Resonance in Medicine (ISMRM) system phantom. Proton density image, and T1 , T2, T1ρ relaxation times, and B1 + were estimated in the NIST/ISMRM system phantom as well as in the human knee medial and lateral femur, medial and lateral tibia, and patellar cartilage. The repeatability and reproducibility of the proposed technique were assessed in the phantom using analysis of the Bland-Altman plots. The intrasubject repeatability was assessed with the coefficient of variation (CV) and root mean square CV (rmsCV). The Mann-Whitney U test was used to assess the difference between healthy subjects and mild knee OA patients. The Bland-Altman plots in the NIST/ISMRM phantom demonstrated an average difference of 0.001% ± 015%, 1.2% ± 7.1%, and 0.47% ± 3% between two scans from the same 3-T scanner (repeatability), and 0.002% ± 015%, 0.62% ± 10.5%, and 0.97% ± 14% between the scans acquired on two different 3-T scanners (reproducibility) for T1 , T2 , and T1ρ , respectively. The in vivo knee study showed excellent repeatability with rmsCV less than 1%, 2%, and 1% for T1 , T2 , and T1ρ , respectively. T1ρ relaxation time in the mild knee OA patients was significantly higher (p < 0.05) than in healthy subjects. The proposed 3D-MRF sequence is fast, reproducible, robust to B1 + inhomogeneity, and can simultaneously measure the T1 , T2 , T1ρ , and B1 + volumetric maps of the knee joint in a single scan within a clinically feasible scan time.

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

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

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.

Cramér-Rao bound-informed training of neural networks for quantitative MRI

PURPOSE

To improve the performance of neural networks for parameter estimation in quantitative MRI, in particular when the noise propagation varies throughout the space of biophysical parameters.

THEORY AND METHODS

A theoretically well-founded loss function is proposed that normalizes the squared error of each estimate with respective Cramér-Rao bound (CRB)-a theoretical lower bound for the variance of an unbiased estimator. This avoids a dominance of hard-to-estimate parameters and areas in parameter space, which are often of little interest. The normalization with corresponding CRB balances the large errors of fundamentally more noisy estimates and the small errors of fundamentally less noisy estimates, allowing the network to better learn to estimate the latter. Further, proposed loss function provides an absolute evaluation metric for performance: A network has an average loss of 1 if it is a maximally efficient unbiased estimator, which can be considered the ideal performance. The performance gain with proposed loss function is demonstrated at the example of an eight-parameter magnetization transfer model that is fitted to phantom and in vivo data.

RESULTS

Networks trained with proposed loss function perform close to optimal, that is, their loss converges to approximately 1, and their performance is superior to networks trained with the standard mean-squared error (MSE). The proposed loss function reduces the bias of the estimates compared to the MSE loss, and improves the match of the noise variance to the CRB. This performance gain translates to in vivo maps that align better with the literature.

CONCLUSION

Normalizing the squared error with the CRB during the training of neural networks improves their performance in estimating biophysical parameters.

A Self-Supervised Deep Learning Reconstruction for Shortening the Breathhold and Acquisition Window in Cardiac Magnetic Resonance Fingerprinting

The aim of this study is to shorten the breathhold and diastolic acquisition window in cardiac magnetic resonance fingerprinting (MRF) for simultaneous T1, T2, and proton spin density (M0) mapping to improve scan efficiency and reduce motion artifacts. To this end, a novel reconstruction was developed that combines low-rank subspace modeling with a deep image prior, termed DIP-MRF. A system of neural networks is used to generate spatial basis images and quantitative tissue property maps, with training performed using only the undersampled k-space measurements from the current scan. This approach avoids difficulties with obtaining in vivo MRF training data, as training is performed de novo for each acquisition. Calculation of the forward model during training is accelerated by using GRAPPA operator gridding to shift spiral k-space data to Cartesian grid points, and by using a neural network to rapidly generate fingerprints in place of a Bloch equation simulation. DIP-MRF was evaluated in simulations and at 1.5 T in a standardized phantom, 18 healthy subjects, and 10 patients with suspected cardiomyopathy. In addition to conventional mapping, two cardiac MRF sequences were acquired, one with a 15-heartbeat(HB) breathhold and 254 ms acquisition window, and one with a 5HB breathhold and 150 ms acquisition window. In simulations, DIP-MRF yielded decreased nRMSE compared to dictionary matching and a sparse and locally low rank (SLLR-MRF) reconstruction. Strong correlation (R2 &gt; 0.999) with T1 and T2 reference values was observed in the phantom using the 5HB/150 ms scan with DIP-MRF. DIP-MRF provided better suppression of noise and aliasing artifacts in vivo, especially for the 5HB/150 ms scan, and lower intersubject and intrasubject variability compared to dictionary matching and SLLR-MRF. Furthermore, it yielded a better agreement between myocardial T1 and T2 from 15HB/254 ms and 5HB/150 ms MRF scans, with a bias of −9 ms for T1 and 2 ms for T2. In summary, this study introduces an extension of the deep image prior framework for cardiac MRF tissue property mapping, which does not require pre-training with in vivo scans, and has the potential to reduce motion artifacts by enabling a shortened breathhold and acquisition window.

Vendor-neutral sequences and fully transparent workflows improve inter-vendor reproducibility of quantitative MRI

PURPOSE

We developed an end-to-end workflow that starts with a vendor-neutral acquisition and tested the hypothesis that vendor-neutral sequences decrease inter-vendor variability of T1, magnetization transfer ratio (MTR), and magnetization transfer saturation-index (MTsat) measurements.

METHODS

We developed and deployed a vendor-neutral 3D spoiled gradient-echo (SPGR) sequence on three clinical scanners by two MRI vendors. We then acquired T1 maps on the ISMRM-NIST system phantom, as well as T1, MTR, and MTsat maps in three healthy participants. We performed hierarchical shift function analysis in vivo to characterize the differences between scanners when the vendor-neutral sequence is used instead of commercial vendor implementations. Inter-vendor deviations were compared for statistical significance to test the hypothesis.

RESULTS

In the phantom, the vendor-neutral sequence reduced inter-vendor differences from 8% to 19.4% to 0.2% to 5% with an overall accuracy improvement, reducing ground truth T1 deviations from 7% to 11% to 0.2% to 4%. In vivo, we found that the variability between vendors is significantly reduced (p = 0.015) for all maps (T1, MTR, and MTsat) using the vendor-neutral sequence.

CONCLUSION

We conclude that vendor-neutral workflows are feasible and compatible with clinical MRI scanners. The significant reduction of inter-vendor variability using vendor-neutral sequences has important implications for qMRI research and for the reliability of multicenter clinical trials.

Simultaneous bilateral T1 , T2 , and T1ρ relaxation mapping of the hip joint with magnetic resonance fingerprinting

Quantitative MRI can detect early biochemical changes in cartilage, but its bilateral use in clinical routines is challenging. The aim of this prospective study was to demonstrate the feasibility of magnetic resonance fingerprinting for bilateral simultaneous T₁ , T₂ , and T_1ρ mapping of the hip joint. The study population consisted of six healthy volunteers with no known trauma or pain in the hip. Monoexponential T₁ , T₂ , and T_1ρ relaxation components were assessed in femoral lateral, superolateral, and superomedial, and inferior, as well as acetabular, superolateral, and superomedial subregions in left and right hip cartilage. Aligned ranked nonparametric factorial analysis was used to assess the side's impact on the subregions. Kruskal-Wallis and Wilcoxon tests were used to compare subregions, and coefficient of variation to assess repeatability. Global averages of T₁ (676.0 ± 45.4 and 687.6 ± 44.5 ms), T₂ (22.5 ± 2.6 and 22.1 ± 2.5 ms), and T_1ρ (38.2 ± 5.5 and 38.2 ± 5.5 ms) were measured in the left and right hip, and articular cartilage, respectively. The Kruskal-Wallis test showed a significant difference between different subregions' relaxation times regardless of the hip side (p < 0.001 for T₁ , p = 0.012 for T₂ , and p < 0.001 for T_1ρ ). The Wilcoxon test showed that T₁ of femoral layers was significantly (p < 0.003) higher than that for acetabular cartilage. The experiments showed excellent repeatability with CV_rms of 1%, 2%, and 4% for T₁ , T₂ , and T_1ρ, respectively. It was concluded that bilateral T₁ , T₂ , and T_1ρ relaxation times, as well as B₁⁺ maps, can be acquired simultaneously from hip joints using the proposed MRF sequence.

High-Quality MR Fingerprinting Reconstruction Using Structured Low-Rank Matrix Completion and Subspace Projection

Due to the capability of fast multiparametric quantitative imaging, magnetic resonance fingerprinting (MRF) is becoming a promising quantitative magnetic resonance imaging approach. However, the artifacts caused by the highly undersampled data acquisition lead to inaccurate estimation of the tissue parameter maps. Based on the assumption that the 3-D MRF data can be modeled as a piecewise smooth signal, with the discontinuities localized to the zero sets of a bandlimited function, we exploit the low-rank property of the structured Toeplitz matrix constructed from the Fourier measurements. In addition, we adopt the subspace projection scheme to improve the accuracy of parameter estimation. In order to efficiently solve the regularized problem, we propose an iterative two-stage algorithm, which alternately updates the k -space data and projects the space-time matrix into the dictionary space. Numerical experiments demonstrate that the proposed algorithm shows significant improvement in MRF time-series images reconstruction and can provide more accurate parameter maps over the state-of-the-art algorithms.

Low‐rank inversion reconstruction for through‐plane accelerated radial MR fingerprinting applied to relaxometry at 0. 35 T

Purpose

To reduce scan time, methods to accelerate phase‐encoded/non‐Cartesian MR fingerprinting (MRF) acquisitions for variable density spiral acquisitions have recently been developed. These methods are not applicable to MRF acquisitions, wherein a single k‐space spoke is acquired per frame. Therefore, we propose a low‐rank inversion method to resolve MRF contrast dynamics from through‐plane accelerated Cartesian/radial measurements applied to quantitative relaxation‐time mapping on a 0.35T system.

Methods

An algorithm was implemented to reconstruct through‐plane aliased low‐rank images describing the contrast dynamics occurring because of the transient‐state MRF acquisition. T₁ and T₂ times from accelerated acquisitions were compared with those from unaccelerated linear reconstructions in a standardized system phantom and within in vivo brain and prostate experiments on a hybrid 0.35T MRI/linear accelerator.

Results

No significant differences between T₁ and T₂ times for the accelerated reconstructions were observed compared to fully sampled acquisitions (p = 0.41 and p = 0.36, respectively). The mean absolute errors in T₁ and T₂ were 5.6% and 2.9%, respectively, between the full and accelerated acquisitions. The SDs in T₁ and T₂ decreased with the advanced accelerated reconstruction compared with the unaccelerated reconstruction (p = 0.02 and p = 0.03, respectively). The quality of the T₁ and T₂ maps generated with the proposed approach are comparable to those obtained using the unaccelerated data sets.

Conclusions

Through‐plane accelerated MRF with radial k‐space coverage was demonstrated at a low field strength of 0.35 T. This method enabled 3D T₁ and T₂ mapping at 0.35 T with a 3‐min scan.

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.

Rapid MR relaxometry using deep learning: An overview of current techniques and emerging trends

Quantitative mapping of MR tissue parameters such as the spin-lattice relaxation time (T1 ), the spin-spin relaxation time (T2 ), and the spin-lattice relaxation in the rotating frame (T1ρ ), referred to as MR relaxometry in general, has demonstrated improved assessment in a wide range of clinical applications. Compared with conventional contrast-weighted (eg T1 -, T2 -, or T1ρ -weighted) MRI, MR relaxometry provides increased sensitivity to pathologies and delivers important information that can be more specific to tissue composition and microenvironment. The rise of deep learning in the past several years has been revolutionizing many aspects of MRI research, including image reconstruction, image analysis, and disease diagnosis and prognosis. Although deep learning has also shown great potential for MR relaxometry and quantitative MRI in general, this research direction has been much less explored to date. The goal of this paper is to discuss the applications of deep learning for rapid MR relaxometry and to review emerging deep-learning-based techniques that can be applied to improve MR relaxometry in terms of imaging speed, image quality, and quantification robustness. The paper is comprised of an introduction and four more sections. Section 2 describes a summary of the imaging models of quantitative MR relaxometry. In Section 3, we review existing "classical" methods for accelerating MR relaxometry, including state-of-the-art spatiotemporal acceleration techniques, model-based reconstruction methods, and efficient parameter generation approaches. Section 4 then presents how deep learning can be used to improve MR relaxometry and how it is linked to conventional techniques. The final section concludes the review by discussing the promise and existing challenges of deep learning for rapid MR relaxometry and potential solutions to address these challenges.

3D MR fingerprinting using Seiffert spirals

PURPOSE

Seiffert spirals were recently explored as an efficient way to traverse 3D k-space compared to traditional 3D techniques. Several studies have shown the ability of 3D MR fingerprinting (MRF) techniques to acquire T₁ and T₂ relaxation maps in a short period of time. However, these sequences do not sample across a large region of 3D k-space every TR, especially in the way that Seiffert trajectories can.

METHODS

A 3D MRF sequence was designed using 8 Seiffert spirals rotated in 3D k-space, with flip angle modulation for T₁ and T₂ sensitivity. The sequence was compared to an MRF sequence using a 2D spiral rotated in 3D k-space using the tiny golden angle acquisition with similar resolution/readout duration. Both sequences were evaluated using simulations, phantom validation, and in vivo imaging.

RESULTS

In all experiments, the Seiffert spiral MRF sequence performed similar to if not better than the multi-axis 2D spiral MRF sequence. Strong intraclass correlation coefficients (> 0.9) were found between conventional and MRF sequences in phantoms, whereas the in vivo results showed slightly less aliasing artifact with the Seiffert trajectory.

CONCLUSION

In this study, Seiffert spirals were used within the MRF framework to acquire high-resolution T₁ and T₂ relaxation time maps in less than 2.5 min. The reduced aliasing artifacts seen with the Seiffert sequence suggests that sampling over 3D k-space evenly each TR can improve quantification or shorten scan times.

Data-Consistent non-Cartesian deep subspace learning for efficient dynamic MR image reconstruction

Non-Cartesian sampling with subspace-constrained image reconstruction is a popular approach to dynamic MRI, but slow iterative reconstruction limits its clinical application. Data-consistent (DC) deep learning can accelerate reconstruction with good image quality, but has not been formulated for non-Cartesian subspace imaging. In this study, we propose a DC non-Cartesian deep subspace learning framework for fast, accurate dynamic MR image reconstruction. Four novel DC formulations are developed and evaluated: two gradient decent approaches, a directly solved approach, and a conjugate gradient approach. We applied a U-Net model with and without DC layers to reconstruct T1-weighted images for cardiac MR Multitasking (an advanced multidimensional imaging method), comparing our results to the iteratively reconstructed reference. Experimental results show that the proposed framework significantly improves reconstruction accuracy over the U-Net model without DC, while significantly accelerating the reconstruction over conventional iterative reconstruction.

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.

Three-dimensional high-resolution T1 and T2 mapping of whole macaque brain at 9.4 T using magnetic resonance fingerprinting

PURPOSE

Quantitative T₁ and T₂ mapping in non-human primates with whole-brain coverage is challenged by the requirement of sub-millimeter resolution and the inhomogeneity of the transmit magnetic field (B₁ ⁺ ) covering a large field of view. The goal of the current study is to develop a magnetic resonance fingerprinting (MRF) method for simultaneous T₁ and T₂ mapping of the entire macaque brain within feasible scan time.

METHODS

A three-dimensional (3D) MRF sequence with both inversion- and T₂ -preparation modules was developed and evaluated on a 9.4 T preclinical scanner. Data acquisition used a 3D stack-of-spirals trajectory, with undersampling along both the in-plane and the through-plane directions. The effect of B₁ ⁺ inhomogeneity was accounted for by matching the acquired fingerprint to a dictionary simulated with the B₁ ⁺ factors measured from a separate scan. In vitro and ex vivo studies were performed to evaluate the accuracy and the undersampling capacity of the MRF method. The application of the MRF method for in vivo, brain-wide T₁ and T₂ mapping was demonstrated on macaques at 4, 6, and 12 years of age.

RESULTS

The MRF method enabled highly repeatable T₁ and T₂ mapping at high spatial resolution (0.35 × 0.35 × 1 mm³ ) with an acceleration factor of 24. In vivo studies showed significant age-related T₂ reduction in deep gray nuclei including the globus pallidus, the putamen, and the caudate nucleus.

CONCLUSIONS

This study demonstrates the first MRF study for brain-wide, multi-parametric quantification in non-human primates with sub-millimeter resolution.

MR Fingerprinting-A Radiogenomic Marker for Diffuse Gliomas

(1) Background: Advanced MR imaging (MRI) of brain tumors is mainly based on qualitative contrast images. MR Fingerprinting (MRF) offers a novel approach. The purpose of this study was to use MRF-derived T1 and T2 relaxation maps to differentiate diffuse gliomas according to isocitrate dehydrogenase (IDH) mutation. (2) Methods: Twenty-four patients with histologically verified diffuse gliomas (14 IDH-mutant, four 1p/19q-codeleted, 10 IDH-wildtype) were enrolled. MRF T1 and T2 relaxation times were compared to apparent diffusion coefficient (ADC), relative cerebral blood volume (rCBV) within solid tumor, peritumoral edema, and normal-appearing white matter (NAWM), using contrast-enhanced MRI, diffusion-, perfusion-, and susceptibility-weighted imaging. For perfusion imaging, a T2* weighted perfusion sequence with leakage correction was used. Correlations of MRF T1 and T2 times with two established conventional sequences for T1 and T2 mapping were assessed (a fast double inversion recovery-based MR sequence ('MP2RAGE') for T1 quantification and a multi-contrast spin echo-based sequence for T2 quantification). (3) Results: MRF T1 and T2 relaxation times were significantly higher in the IDH-mutant than in IDH-wildtype gliomas within the solid part of the tumor (p = 0.024 for MRF T1, p = 0.041 for MRF T2). MRF T1 and T2 relaxation times were significantly higher in the IDH-wildtype than in IDH-mutant gliomas within peritumoral edema less than or equal to 1cm adjacent to the tumor (p = 0.038 for MRF T1 mean, p = 0.010 for MRF T2 mean). In the solid part of the tumor, there was a high correlation between MRF and conventionally measured T1 and T2 values (r = 0.913, p < 0.001 for T1, r = 0.775, p < 0.001 for T2), as well as between MRF and ADC values (r = 0.813, p < 0.001 for T2, r = 0.697, p < 0.001 for T1). The correlation was weak between the MRF and rCBV values (r = -0.374, p = 0.005 for T2, r = -0.181, p = 0.181 for T1). (4) Conclusions: MRF enables fast, single-sequence based, multi-parametric, quantitative tissue characterization of diffuse gliomas and may have the potential to differentiate IDH-mutant from IDH-wildtype gliomas.

Myocardial T1, T2, T2*, and fat fraction quantification via low-rank motion-corrected cardiac MR fingerprinting

Purpose

Develop a novel 2D cardiac MR fingerprinting (MRF) approach to enable simultaneous T1, T2, T2*, and fat fraction (FF) myocardial tissue characterization in a single breath‐hold scan.

Methods

Simultaneous, co‐registered, multi‐parametric mapping of T1, T2, and FF has been recently achieved with cardiac MRF. Here, we further incorporate T2* quantification within this approach, enabling simultaneous T1, T2, T2*, and FF myocardial tissue characterization in a single breath‐hold scan. T2* quantification is achieved with an eight‐echo readout that requires a long cardiac acquisition window. A novel low‐rank motion‐corrected (LRMC) reconstruction is exploited to correct for cardiac motion within the long acquisition window. The proposed T1/T2/T2*/FF cardiac MRF was evaluated in phantom and in 10 healthy subjects in comparison to conventional mapping techniques.

Results

The proposed approach achieved high quality parametric mapping of T1, T2, T2*, and FF with corresponding normalized RMS error (RMSE) T1 = 5.9%, T2 = 9.6% (T2 values <100 ms), T2* = 3.3% (T2* values <100 ms), and FF = 0.8% observed in phantom scans. In vivo, the proposed approach produced higher left‐ventricular myocardial T1 values than MOLLI (1148 vs 1056 ms), lower T2 values than T2‐GraSE (42.8 vs 50.6 ms), lower T2* values than eight‐echo gradient echo (GRE) (35.0 vs 39.4 ms), and higher FF values than six‐echo GRE (0.8 vs 0.3 %) reference techniques. The proposed approach achieved considerable reduction in motion artifacts compared to cardiac MRF without motion correction, improved spatial uniformity, and statistically higher apparent precision relative to conventional mapping for all parameters.

Conclusion

The proposed cardiac MRF approach enables simultaneous, co‐registered mapping of T1, T2, T2*, and FF in a single breath‐hold for comprehensive myocardial tissue characterization, achieving higher apparent precision than conventional methods.

An MR fingerprinting approach for quantitative inhomogeneous magnetization transfer imaging

PURPOSE

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

THEORY AND METHODS

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

RESULTS

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

CONCLUSION

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