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

DGEDDGAN: A dual-domain generator and edge-enhanced dual discriminator generative adversarial network for MRI reconstruction

Magnetic resonance imaging (MRI) as a critical clinical tool in medical imaging, requires a long scan time for producing high-quality MRI images. To accelerate the speed of MRI while reconstructing high-quality images with sharper edges and fewer aliases, a novel dual-domain generator and edge-enhancement dual discriminator generative adversarial network structure named DGEDDGAN for MRI reconstruction is proposed, in which one discriminator is responsible for holistic image reconstruction, whereas the other is adopted to enhance the edge preservation. A dual-domain U-Net structure that cascades the frequency domain and image domain is designed for the generator. The densely connected residual block is used to replace the traditional U-Net convolution block to improve the feature reuse capability while overcoming the gradient vanishing problem. The coordinate attention mechanism in each skip connection is employed to effectively reduce the loss of spatial information and enforce the feature selection capability. Extensive experiments on two publicly available datasets i.e., IXI dataset and CC-359, demonstrate that the proposed method can reconstruct the high-quality MRI images with more edge details and fewer artifacts, outperforming several state-of-the-art methods under various sampling rates and masks. The time of single-image reconstruction is below 13 ms, which meets the demand of faster processing.

Deep graph embedding based on Laplacian eigenmaps for MR fingerprinting reconstruction

Magnetic resonance fingerprinting (MRF) is a promising technique for fast quantitative imaging of multiple tissue parameters. However, the highly undersampled schemes utilized in MRF typically lead to noticeable aliasing artifacts in reconstructed images. Existing model-based methods can mitigate aliasing artifacts and enhance reconstruction quality but suffer from long reconstruction times. In addition, data priors used in these methods, such as low-rank and total variation, make it challenging to incorporate non-local and non-linear redundancies in MRF data. Furthermore, existing deep learning-based methods for MRF often lack interpretability and struggle with the high computational overhead caused by the high dimensionality of MRF data. To address these issues, we introduce a novel deep graph embedding framework based on the Laplacian eigenmaps for improved MRF reconstruction. Our work first models the acquired high-dimensional MRF data and the corresponding parameter maps as graph data nodes. Then, we propose an MRF reconstruction framework based on the graph embedding framework, retaining intrinsic graph structures between parameter maps and MRF data. To improve the accuracy of the estimated graph structure and the computational efficiency of the proposed framework, we unroll the iterative optimization process into a deep neural network, incorporating a learned graph embedding module to adaptively learn the Laplacian eigenmaps. By introducing the graph embedding framework into the MRF reconstruction, the proposed method can effectively exploit non-local and non-linear correlations in MRF data. Numerical experiments demonstrate that our approach can reconstruct high-quality MRF data and multiple parameter maps within a significantly reduced computational cost.

Feasibility of 3D MRI fingerprinting for rapid knee cartilage T1, T2, and T1ρ mapping at 0.55T: Comparison with 3T

Low-field strength scanners present an opportunity for more inclusive imaging exams and bring several challenges including lower signal-to-noise ratio (SNR) and longer scan times. Magnetic resonance fingerprinting (MRF) is a rapid quantitative multiparametric method that can enable multiple quantitative maps simultaneously. To demonstrate the feasibility of an MRF sequence for knee cartilage evaluation in a 0.55T system we performed repeatability and accuracy experiments with agar-gel phantoms. Additionally, five healthy volunteers (age 32 ± 4 years old, 2 females) were scanned at 3T and 0.55T. The MRI acquisition protocols include a stack-of-stars T1ρ-enabled MRF sequence, a VIBE sequence with variable flip angles (VFA) for T1 mapping, and fat-suppressed turbo flash (TFL) sequences for T2 and T1ρ mappings. Double-Echo steady-state (DESS) sequence was also used for cartilage segmentation. Acquisitions were performed at two different field strengths, 0.55T and 3T, with the same sequences but protocols were slightly different to accommodate differences in signal-to-noise ratio and relaxation times. Cartilage segmentation was done using five compartments. T1, T2, and T1ρ values were measured in the knee cartilage using both MRF and conventional relaxometry sequences. The MRF sequence demonstrated excellent repeatability in a test-retest experiment with model agar-gel phantoms, as demonstrated with correlation and Bland-Altman plots. Underestimation of T1 values was observed on both field strengths, with the average global difference between reference values and the MRF being 151 ms at 0.55T and 337 ms at 3T. At 0.55T, MRF measurements presented significant biases but strong correlations with the reference measurements. Although a larger error was present in T1 measurements, MRF measurements trended similarly to the conventional measurements for human subjects and model agar-gel phantoms.

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.

Deep magnetic resonance fingerprinting based on Local and Global Vision Transformer

To mitigate systematic errors in magnetic resonance fingerprinting (MRF), the precomputed dictionary is usually computed with minimal granularity across the entire range of tissue parameters. However, the dictionary grows exponentially with the number of parameters increase, posing significant challenges to the computational efficiency and matching accuracy of pattern-matching algorithms. Existing works, primarily based on convolutional neural networks (CNN), focus solely on local information to reconstruct multiple parameter maps, lacking in-depth investigations on the MRF mechanism. These methods may not exploit long-distance redundancies and the contextual information within voxel fingerprints introduced by the Bloch equation dynamics, leading to limited reconstruction speed and accuracy. To overcome these limitations, we propose a novel end-to-end neural network called the Local and Global Vision Transformer (LG-ViT) for MRF parameter reconstruction. Our proposed LG-ViT employs a multi-stage architecture that effectively reduces the computational overhead associated with the high-dimensional MRF data and the transformer model. Specifically, a local Transformer encoder is proposed to capture contextual information embedded within voxel fingerprints and local correlations introduced by the interconnected human tissues. Additionally, a global Transformer encoder is proposed to leverage long-distance dependencies arising from shared characteristics among different tissues across various spatial regions. By incorporating MRF physics-based data priors and effectively capturing local and global correlations, our proposed LG-ViT can achieve fast and accurate MRF parameter reconstruction. Experiments on both simulation and in vivo data demonstrate that the proposed method enables faster and more accurate MRF parameter reconstruction compared to state-of-the-art deep learning-based methods.

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.

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.

Learned Tensor Low-CP-Rank and Bloch Response Manifold Priors for Non-Cartesian MRF Reconstruction

Magnetic resonance fingerprinting (MRF) can rapidly perform simultaneous imaging of multiple tissue parameters. However, the rapid acquisition schemes used in MRF inevitably introduce aliasing artifacts in the recovered tissue fingerprints, reducing the accuracy of the predicted parameter maps. Current regularized reconstruction methods are based on iterative procedures which are usually time-consuming. In addition, most of the current deep learning-based methods for MRF often lack interpretability owing to the black-box nature, and most deep learning-based methods are not applicable for non-Cartesian scenarios, which limits the practical applications. In this paper, we propose a joint reconstruction model incorporating MRF-physics prior and the data correlation constraint for non-Cartesian MRF reconstruction. To avoid time-consuming iterative procedures, we unroll the reconstruction model into a deep neural network. Specifically, we propose a learned CANDECOMP/PARAFAC (CP) decomposition module to exploit the tensor low-rank priors of high-dimensional MRF data, which avoids computationally burdensome singular value decomposition. Inspired by the MRF-physics, we also propose a Bloch response manifold module to learn the mapping between reconstructed MRF data and the multiple parameter maps. Numerical experiments show that the proposed network can reconstruct high-quality MRF data and multiple parameter maps within significantly reduced computational time.

Improving motion robustness of 3D MR fingerprinting with a fat navigator

PURPOSE

To develop a 3D MR fingerprinting (MRF) method in combination with fat navigators to improve its motion robustness for neuroimaging.

METHODS

A rapid fat navigator was developed using the stack-of-spirals acquisition and non-Cartesian spiral GRAPPA. The fat navigator module was implemented in the 3D MRF sequence with high scan efficiency. The developed method was first validated in phantoms and five healthy subjects with intentional head motion. The method was further applied to infants with neonatal opioid withdrawal symptoms. The 3D MRF scans with fat navigators acquired with and without acceleration along the partition-encoding direction were both examined in the study.

RESULTS

Both phantom and in vivo results demonstrated that the added fat navigator modules did not influence the quantification accuracy in MRF. In combination with non-Cartesian spiral GRAPPA, a rapid fat navigator sampling with whole-brain coverage was achieved in ˜0.5 s at 3T, reducing its sensitivity to potential motion. Based on the motion waveforms extracted from fat navigators, the motion robustness of the 3D MRF was largely improved. With the proposed method, the motion-corrupted MRF datasets yielded T1 and T2 maps with significantly reduced artifacts and high correlations with measurements from the reference motion-free MRF scans.

CONCLUSION

We developed a 3D MRF method coupled with rapid fat navigators to improve its motion robustness for quantitative neuroimaging. Our results demonstrate that (1) accurate tissue quantification was preserved with the fat navigator modules and (2) the motion robustness for quantitative tissue mapping was largely improved with the developed method.

Updates on Compositional MRI Mapping of the Cartilage: Emerging Techniques and Applications

Osteoarthritis (OA) is a widely occurring degenerative joint disease that is severely debilitating and causes significant socioeconomic burdens to society. Magnetic resonance imaging (MRI) is the preferred imaging modality for the morphological evaluation of cartilage due to its excellent soft tissue contrast and high spatial resolution. However, its utilization typically involves subjective qualitative assessment of cartilage. Compositional MRI, which refers to the quantitative characterization of cartilage using a variety of MRI methods, can provide important information regarding underlying compositional and ultrastructural changes that occur during early OA. Cartilage compositional MRI could serve as early imaging biomarkers for the objective evaluation of cartilage and help drive diagnostics, disease characterization, and response to novel therapies. This review will summarize current and ongoing state-of-the-art cartilage compositional MRI techniques and highlight emerging methods for cartilage compositional MRI including MR fingerprinting, compressed sensing, multiexponential relaxometry, improved and robust radio-frequency pulse sequences, and deep learning-based acquisition, reconstruction, and segmentation. The review will also briefly discuss the current challenges and future directions for adopting these emerging cartilage compositional MRI techniques for use in clinical practice and translational OA research studies. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 2.

Optimized three-dimensional ultrashort echo time: Magnetic resonance fingerprinting for myelin tissue fraction mapping

Quantitative assessment of brain myelination has gained attention for both research and diagnosis of neurological diseases. However, conventional pulse sequences cannot directly acquire the myelin-proton signals due to its extremely short T2 and T2* values. To obtain the myelin-proton signals, dedicated short T2 acquisition techniques, such as ultrashort echo time (UTE) imaging, have been introduced. However, it remains challenging to isolate the myelin-proton signals from tissues with longer T2. In this article, we extended our previous two-dimensional ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) with dual-echo acquisition to three dimensional (3D). Given a relatively low proton density (PD) of myelin-proton, we utilized Cramér-Rao Lower Bound to encode myelin-proton with the maximal SNR efficiency for optimizing the MR fingerprinting design, in order to improve the sensitivity of the sequence to myelin-proton. In addition, with a second echo of approximately 3 ms, myelin-water component can be also captured. A myelin-tissue (myelin-proton and myelin-water) fraction mapping can be thus calculated. The optimized 3D UTE-MRF with dual-echo acquisition is tested in simulations, physical phantom and in vivo studies of both healthy subjects and multiple sclerosis patients. The results suggest that the rapidly decayed myelin-proton and myelin-water signal can be depicted with UTE signals of our method at clinically relevant resolution (1.8 mm isotropic) in 15 min. With its good sensitivity to myelin loss in multiple sclerosis patients demonstrated, our method for the whole brain myelin-tissue fraction mapping in clinical friendly scan time has the potential for routine clinical imaging.

Accelerating MR Parameter Mapping Using Nonlinear Compressive Manifold Learning and Regularized Pre-Imaging

In this study, we presented a novel method to reconstruct the MR parametric maps from highly undersampled k-space data. Specifically, we utilized a nonlinear model to sparsely represent the unknown MR parameter-weighted images in high-dimensional feature space. Each image at a specific time point is assumed to belong to a low-dimensional manifold which is learned from training images created based on the parametric model. The final reconstruction is carried out by venturing the sparse representation of the images in the feature space back to the input space, using the pre-imaging technique. Particularly, among an infinite number of solutions that satisfy the data consistency, the one that is closest to the manifold is selected as the desired solution. The underlying optimization problem is solved using kernel trick, sparse coding, and split Bregman iteration algorithm. In addition, both spatial and temporal regularizations were utilized to further improve the reconstruction quality. The proposed method was validated on both phantom and in vivo human brain T2 mapping data. Results showed the proposed method was superior to the conventional linear model-based reconstruction methods, in terms of artifact removal and quantitative estimate accuracy. The proposed method could be potentially beneficial for quantitative MR applications.