High-dimensionality undersampled patch-based reconstruction (HD-PROST) for accelerated multi-contrast MRI

Accelerated model-based T1, T2* and proton density mapping using a Bayesian approach with automatic hyperparameter estimation

PURPOSE

To achieve automatic hyperparameter estimation for the model-based recovery of quantitative MR maps from undersampled data, we propose a Bayesian formulation that incorporates the signal model and sparse priors among multiple image contrasts.

THEORY

We introduce a novel approximate message passing framework "AMP-PE" that enables the automatic and simultaneous recovery of hyperparameters and quantitative maps.

METHODS

We employed the variable-flip-angle method to acquire multi-echo measurements using gradient echo sequence. We explored undersampling schemes to incorporate complementary sampling patterns across different flip angles and echo times. We further compared AMP-PE with conventional compressed sensing approaches such as thel 1 $$ {l}_1 $$ -norm minimization, PICS and other model-based approaches such as GraSP, MOBA.

RESULTS

Compared to conventional compressed sensing approaches such as thel 1 $$ {l}_1 $$ -norm minimization and PICS, AMP-PE achieved superior reconstruction performance with lower errors inT 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ mapping and comparable performance inT 1 $$ {\mathrm{T}}_1 $$ and proton density mappings. When compared to other model-based approaches including GraSP and MOBA, AMP-PE exhibited greater robustness and outperformed GraSP in reconstruction error. AMP-PE offers faster speed than MOBA. AMP-PE performed better than MOBA at higher sampling rates and worse than MOBA at a lower sampling rate. Notably, AMP-PE eliminates the need for hyperparameter tuning, which is a requisite for all the other approaches.

CONCLUSION

AMP-PE offers the benefits of model-based recovery with the additional key advantage of automatic hyperparameter estimation. It works adeptly in situations where ground-truth is difficult to obtain and in clinical environments where it is desirable to automatically adapt hyperparameters to individual protocol, scanner and patient.

3D B1+ corrected simultaneous myocardial T1 and T1ρ mapping with subject-specific respiratory motion correction and water-fat separation

PURPOSE

To develop a 3D free-breathing cardiac multi-parametric mapping framework that is robust to confounders of respiratory motion, fat, and B1+ inhomogeneities and validate it for joint myocardial T1 and T1ρ mapping at 3T.

METHODS

An electrocardiogram-triggered sequence with dual-echo Dixon readout was developed, where nine cardiac cycles were repeatedly acquired with inversion recovery and T1ρ preparation pulses for T1 and T1ρ sensitization. A subject-specific respiratory motion model relating the 1D diaphragmatic navigator to the respiration-induced 3D translational motion of the heart was constructed followed by respiratory motion binning and intra-bin 3D translational and inter-bin non-rigid motion correction. Spin history B1+ inhomogeneities were corrected with optimized dual flip angle strategy. After water-fat separation, the water images were matched to the simulated dictionary for T1 and T1ρ quantification. Phantoms and 10 heathy subjects were imaged to validate the proposed technique.

RESULTS

The proposed technique achieved strong correlation (T1: R2 = 0.99; T1ρ: R2 = 0.98) with the reference measurements in phantoms. 3D cardiac T1 and T1ρ maps with spatial resolution of 2 × 2 × 4 mm were obtained with scan time of 5.4 ± 0.5 min, demonstrating comparable T1 (1236 ± 59 ms) and T1ρ (50.2 ± 2.4 ms) measurements to 2D separate breath-hold mapping techniques. The estimated B1+ maps showed spatial variations across the left ventricle with the septal and inferior regions being 10%-25% lower than the anterior and septal regions.

CONCLUSION

The proposed technique achieved efficient 3D joint myocardial T1 and T1ρ mapping at 3T with respiratory motion correction, spin history B1+ correction and water-fat separation.

Partition-based k-space synthesis for multi-contrast parallel imaging

PURPOSE

Multi-contrast magnetic resonance imaging is a significant and essential medical imaging technique. However, multi-contrast imaging has longer acquisition time and is easy to cause motion artifacts. In particular, the acquisition time for a T2-weighted image is prolonged due to its longer repetition time (TR). On the contrary, T1-weighted image has a shorter TR. Therefore, utilizing complementary information across T1 and T2-weighted image is a way to decrease the overall imaging time. Previous T1-assisted T2 reconstruction methods have mostly focused on image domain using whole-based image fusion approaches. The image domain reconstruction method has the defects of high computational complexity and limited flexibility. To address this issue, we propose a novel multi-contrast imaging method called partition-based k-space synthesis (PKS) which can achieve better reconstruction quality of T2-weighted image by feature fusion.

METHODS

Concretely, we first decompose fully-sampled T1 k-space data and under-sampled T2 k-space data into two sub-data, separately. Then two new objects are constructed by combining the two sub-T1/T2 data. After that, the two new objects as the whole data to realize the reconstruction of T2-weighted image.

RESULTS

Experimental results showed that the developed PKS scheme can achieve comparable or better results than using traditional k-space parallel imaging (SAKE) that processes each contrast independently. At the same time, our method showed good adaptability and robustness under different contrast-assisted and T1-T2 ratios. Efficient target modal image reconstruction under various conditions were realized and had excellent performance in restoring image quality and preserving details.

CONCLUSIONS

This work proposed a PKS multi-contrast method to assist in target mode image reconstruction. We have conducted extensive experiments on different multi-contrast, diverse ratios of T1 to T2 and different sampling masks to demonstrate the generalization and robustness of our proposed model.

ACGRHA-Net: Accelerated multi-contrast MR imaging with adjacency complementary graph assisted residual hybrid attention network

Multi-contrast magnetic resonance (MR) imaging is an advanced technology used in medical diagnosis, but the long acquisition process can lead to patient discomfort and limit its broader application. Shortening acquisition time by undersampling k-space data introduces noticeable aliasing artifacts. To address this, we propose a method that reconstructs multi-contrast MR images from zero-filled data by utilizing a fully-sampled auxiliary contrast MR image as a prior to learn an adjacency complementary graph. This graph is then combined with a residual hybrid attention network, forming the adjacency complementary graph assisted residual hybrid attention network (ACGRHA-Net) for multi-contrast MR image reconstruction. Specifically, the optimal structural similarity is represented by a graph learned from the fully sampled auxiliary image, where the node features and adjacency matrices are designed to precisely capture structural information among different contrast images. This structural similarity enables effective fusion with the target image, improving the detail reconstruction. Additionally, a residual hybrid attention module is designed in parallel with the graph convolution network, allowing it to effectively capture key features and adaptively emphasize these important features in target contrast MR images. This strategy prioritizes crucial information while preserving shallow features, thereby achieving comprehensive feature fusion at deeper levels to enhance multi-contrast MR image reconstruction. Extensive experiments on the different datasets, using various sampling patterns and accelerated factors demonstrate that the proposed method outperforms the current state-of-the-art reconstruction methods.

A Feasibility Study of AI-Assisted Compressed Sensing in Prostate T2-Weighted Imaging

RATIONALE AND OBJECTIVES

To evaluate the image quality and PI-RADS scoring performance of prostate T2-weighted imaging (T2WI) based on AI-assisted compressed sensing (ACS).

MATERIALS AND METHODS

In this prospective study, adult male urological outpatients or inpatients underwent prostate MRI, including T2WI, diffusion-weighted imaging and apparent diffusion coefficient maps. Three accelerated scanning protocols using parallel imaging (PI) and ACS: T2WIPI, T2WIACS1 and T2WIACS2 were evaluated through comparative analysis. Quantitative analysis included signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), slope profile, and edge rise distance (ERD). Image quality was qualitatively assessed using a five-point Likert scale (ranging from 1 = non-diagnostic to 5 = excellent). PI-RADS scores were determined for the largest or most suspicious lesions in each patient. The Friedman test and one-way ANOVA with post hoc tests were utilized for group comparisons, with statistical significance set at P < 0.05.

RESULTS

This study included 40 participants. Compared to PI, ACS reduced acquisition time by over 50%, significantly enhancing the CNR of sagittal and axial T2WI (P < 0.05), significantly improving the image quality of sagittal and axial T2WI (P < 0.05). No significant differences were observed in slope profile, ERD, and PI-RADS scores between groups (P > 0.05).

CONCLUSION

ACS reduced prostate T2WI acquisition time by half while improving image quality without affecting PI-RADS scores.

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.

Accelerated Cartesian cardiac T2 mapping based on a calibrationless locally low-rank tensor constraint

Background

Cardiac T2 mapping is a valuable tool for diagnosing myocardial edema, inflammation, and infiltration, yet its spatial resolution is limited by the single-shot balanced steady-state free precession acquisition and duration of the cardiac quiescent period, which may reduce sensitivity in detecting focal lesions in the myocardium. To improve spatial resolution without extending the acquisition window, this study examined a novel accelerated Cartesian cardiac T2 mapping technique.

Methods

We introduce a novel improved-resolution cardiac T2 mapping approach leveraging a calibrationless space-contrast-coil locally low-rank tensor (SCC-LLRT)-constrained reconstruction algorithm in conjunction with Cartesian undersampling trajectory. The method was validated with phantom imaging and in vivo imaging that involved 13 healthy participants and 20 patients. The SCC-LLRT algorithm was compared with a conventional locally low-rank (LLR)-constrained algorithm and a nonlinear inversion (NLINV) reconstruction algorithm. The improved-resolution T2 mapping (1.4 mm × 1.4 mm) was compared globally and regionally with the regular-resolution T2 mapping (2.3 mm × 1.9 mm) according to the 16-segment model of the American Heart Association. The agreement between the improved-resolution and regular-resolution T2 mappings was evaluated by linear regression and Bland-Altman analyses. Image quality was scored by two experienced reviewers on a five-point scale (1, worst; 5, best).

Results

In healthy participants, SCC-LLRT significantly reduced artifacts (4.50±0.39) compared with LLR (2.31±0.60; P<0.001) and NLINV (3.65±0.56; P<0.01), suppressed noise (4.12±0.35) compared with NLINV (2.65±0.50; P<0.001), and improved the overall image quality (4.38±0.40) compared with LLR (2.54±0.41; P<0.001) and NLINV (3.04±0.50; P<0.001). Compared with the regular-resolution T2 mapping, the proposed method significantly improved the sharpness of myocardial boundaries (4.46±0.60 vs. 3.04±0.50; P<0.001) and the conspicuity of papillary muscles and fine structures (4.46±0.63 vs. 2.65±0.30; P<0.001). Myocardial T2 values obtained with the proposed method correlated significantly with those from regular-resolution T2 mapping in both healthy participants (r=0.79; P<0.01) and patients (r=0.94; P<0.001).

Conclusions

The proposed SCC-LLRT-constrained reconstruction algorithm in conjunction with Cartesian undersampling pattern achieved improved-resolution cardiac T2 mapping of comparable accuracy, precision, and scan-rescan reproducibility compared with the regular-resolution T2 mapping. The higher resolution improved the sharpness of myocardial borders and the conspicuity of image fine details, which may increase diagnostic confidence in cardiac T2 mapping for detecting small lesions.

Alternating low-rank tensor reconstruction for improved multiparametric mapping with cardiovascular MR Multitasking

PURPOSE

To develop a novel low-rank tensor reconstruction approach leveraging the complete acquired data set to improve precision and repeatability of multiparametric mapping within the cardiovascular MR Multitasking framework.

METHODS

A novel approach that alternated between estimation of temporal components and spatial components using the entire data set acquired (i.e., including navigator data and imaging data) was developed to improve reconstruction. The precision and repeatability of the proposed approach were evaluated on numerical simulations, 10 healthy subjects, and 10 cardiomyopathy patients at multiple scan times for 2D myocardial T1/T2 mapping with MR Multitasking and were compared with those of the previous navigator-derived fixed-basis approach.

RESULTS

In numerical simulations, the proposed approach outperformed the previous fixed-basis approach with lower T1 and T2 error against the ground truth at all scan times studied and showed better motion fidelity. In human subjects, the proposed approach showed no significantly different sharpness or T1/T2 measurement and significantly improved T1 precision by 20%-25%, T2 precision by 10%-15%, T1 repeatability by about 30%, and T2 repeatability by 25%-35% at 90-s and 50-s scan times The proposed approach at the 50-s scan time also showed comparable results with that of the previous fixed-basis approach at the 90-s scan time.

CONCLUSION

The proposed approach improved precision and repeatability for quantitative imaging with MR Multitasking while maintaining comparable motion fidelity, T1/T2 measurement, and septum sharpness and had the potential for further reducing scan time from 90 s to 50 s.

Free-breathing 3D whole-heart joint T1/T2 mapping and water/fat imaging at 0.55 T

PURPOSE

To develop and validate a highly efficient motion compensated free-breathing isotropic resolution 3D whole-heart joint T1/T2 mapping sequence with anatomical water/fat imaging at 0.55 T.

METHODS

The proposed sequence takes advantage of shorter T1 at 0.55 T to acquire three interleaved water/fat volumes with inversion-recovery preparation, no preparation, and T2 preparation, respectively. Image navigators were used to facilitate nonrigid motion-compensated image reconstruction. T1 and T2 maps were jointly calculated by a dictionary matching method. Validations were performed with simulation, phantom, and in vivo experiments on 10 healthy volunteers and 1 patient. The performance of the proposed sequence was compared with conventional 2D mapping sequences including modified Look-Locker inversion recovery and T2-prepared balanced steady-SSFP sequence.

RESULTS

The proposed sequence has a good T1 and T2 encoding sensitivity in simulation, and excellent agreement with spin-echo reference T1 and T2 values was observed in a standardized T1/T2 phantom (R2 = 0.99). In vivo experiments provided good-quality co-registered 3D whole-heart T1 and T2 maps with 2-mm isotropic resolution in a short scan time of about 7 min. For healthy volunteers, left-ventricle T1 mean and SD measured by the proposed sequence were both comparable with those of modified Look-Locker inversion recovery (640 ± 35 vs. 630 ± 25 ms [p = 0.44] and 49.9 ± 9.3 vs. 54.4 ± 20.5 ms [p = 0.42]), whereas left-ventricle T2 mean and SD measured by the proposed sequence were both slightly lower than those of T2-prepared balanced SSFP (53.8 ± 5.5 vs. 58.6 ± 3.3 ms [p < 0.01] and 5.2 ± 0.9 vs. 6.1 ± 0.8 ms [p = 0.03]). Myocardial T1 and T2 in the patient measured by the proposed sequence were in good agreement with conventional 2D sequences and late gadolinium enhancement.

CONCLUSION

The proposed sequence simultaneously acquires 3D whole-heart T1 and T2 mapping with anatomical water/fat imaging at 0.55 T in a fast and efficient 7-min scan. Further investigation in patients with cardiovascular disease is now warranted.

Evaluation of myocarditis with a free-breathing three-dimensional isotropic whole-heart joint T1 and T2 mapping sequence

BACKGROUND

The diagnosis of myocarditis by cardiovascular magnetic resonance (CMR) requires the use of T2 and T1 weighted imaging, ideally incorporating parametric mapping. Current two-dimensional (2D) mapping sequences are acquired sequentially and involve multiple breath-holds resulting in prolonged scan times and anisotropic image resolution. We developed an isotropic free-breathing three-dimensional (3D) whole-heart sequence that allows simultaneous T1 and T2 mapping and validated it in patients with suspected myocarditis.

METHODS

Eighteen healthy volunteers and 28 patients with suspected myocarditis underwent conventional 2D T1 and T2 mapping with whole-heart coverage and 3D joint T1/T2 mapping on a 1.5T scanner. Acquisition time, image quality, and diagnostic performance were compared. Qualitative analysis was performed using a 4-point Likert scale. Bland-Altman plots were used to assess the quantitative agreement between 2D and 3D sequences.

RESULTS

The 3D T1/T2 sequence was acquired in 8 min 26 s under free breathing, whereas 2D T1 and T2 sequences were acquired with breath-holds in 11 min 44 s (p = 0.0001). All 2D images were diagnostic. For 3D images, 89% (25/28) of T1 and 96% (27/28) of T2 images were diagnostic with no significant difference in the proportion of diagnostic images for the 3D and 2D T1 (p = 0.2482) and T2 maps (p = 1.0000). Systematic bias in T1 was noted with biases of 102, 115, and 152 ms for basal-apical segments, with a larger bias for higher T1 values. Good agreement between T2 values for 3D and 2D techniques was found (bias of 1.8, 3.9, and 3.6 ms for basal-apical segments). The sensitivity and specificity of the 3D sequence for diagnosing acute myocarditis were 74% (95% confidence interval [CI] 49%-91%) and 83% (36%-100%), respectively, with a c-statistic (95% CI) of 0.85 (0.79-0.91) and no statistically significant difference between the 2D and 3D sequences for the detection of acute myocarditis for T1 (p = 0.2207) or T2 (p = 1.0000).

CONCLUSION

Free-breathing whole-heart 3D joint T1/T2 mapping was comparable to 2D mapping sequences with respect to diagnostic performance, but with the added advantages of free breathing and shorter scan times. Further work is required to address the bias noted at high T1 values, but this did not significantly impact diagnostic accuracy.

A hybrid deep image prior and compressed sensing reconstruction method for highly accelerated 3D coronary magnetic resonance angiography

Introduction

High-resolution whole-heart coronary magnetic resonance angiography (CMRA) often suffers from unreasonably long scan times, rendering imaging acceleration highly desirable. Traditional reconstruction methods used in CMRA rely on either hand-crafted priors or supervised learning models. Although the latter often yield superior reconstruction quality, they require a large amount of training data and memory resources, and may encounter generalization issues when dealing with out-of-distribution datasets.

Methods

To address these challenges, we introduce an unsupervised reconstruction method that combines deep image prior (DIP) with compressed sensing (CS) to accelerate 3D CMRA. This method incorporates a slice-by-slice DIP reconstruction and 3D total variation (TV) regularization, enabling high-quality reconstruction under a significant acceleration while enforcing continuity in the slice direction. We evaluated our method by comparing it to iterative SENSE, CS-TV, CS-wavelet, and other DIP-based variants, using both retrospectively and prospectively undersampled datasets.

Results

The results demonstrate the superiority of our 3D DIP-CS approach, which improved the reconstruction accuracy relative to the other approaches across both datasets. Ablation studies further reveal the benefits of combining DIP with 3D TV regularization, which leads to significant improvements of image quality over pure DIP-based methods. Evaluation of vessel sharpness and image quality scores shows that DIP-CS improves the quality of reformatted coronary arteries.

Discussion

The proposed method enables scan-specific reconstruction of high-quality 3D CMRA from a five-minute acquisition, without relying on fully-sampled training data or placing a heavy burden on memory resources.

High Spatial-Resolution and Acquisition-Efficiency Cardiac MR T1 Mapping Based on Radial bSSFP and a Low-Rank Tensor Constraint

BACKGROUND

Cardiac T1 mapping is valuable for evaluating myocardial fibrosis, yet its resolution and acquisition efficiency are limited, potentially obscuring visualization of small pathologies.

PURPOSE

To develop a technique for high-resolution cardiac T1 mapping with a less-than-100-millisecond acquisition window based on radial MOdified Look-Locker Inversion recovery (MOLLI) and a calibrationless space-contrast-coil locally low-rank tensor (SCC-LLRT) constrained reconstruction.

STUDY TYPE

Prospective.

SUBJECTS/PHANTOM

Sixteen healthy subjects (age 25 ± 3 years, 44% females) and 12 patients with suspected cardiomyopathy (age 57 ± 15 years, 42% females), NiCl2-agar phantom.

FIELD STRENGTH/SEQUENCE

3-T, standard MOLLI, radial MOLLI, inversion-recovery spin-echo, late gadolinium enhancement.

ASSESSMENT

SCC-LLRT was compared to a conventional locally low-rank (LLR) method through simulations using Normalized Root-Mean-Square Error (NRMSE) and Structural Similarity Index Measure (SSIM). Radial MOLLI was compared to standard MOLLI across phantom, healthy subjects, and patients. Three independent readers subjectively evaluated the quality of T1 maps using a 5-point scale (5 = best).

STATISTICAL TESTS

Paired t-test, Wilcoxon signed-rank test, intraclass correlation coefficient analysis, linear regression, Bland-Altman analysis. P < 0.05 was considered statistically significant.

RESULTS

In simulations, SCC-LLRT demonstrated a significant improvement in NRMSE and SSIM compared to LLR. In phantom, both radial MOLLI and standard MOLLI provided consistent T1 estimates across different heart rates. In healthy subjects, radial MOLLI exhibited a significantly lower mean T1 (1115 ± 39 msec vs. 1155 ± 36 msec), similar T1 SD (74 ± 14 msec vs. 67 ± 23 msec, P = 0.20), and similar T1 reproducibility (28 ± 18 msec vs. 22 ± 15 msec, P = 0.34) compared to standard MOLLI. In patients, the proposed method significantly improved the sharpness of myocardial boundaries (4.50 ± 0.65 vs. 3.25 ± 0.43), the conspicuity of papillary muscles and fine structures (4.33 ± 0.74 vs. 3.33 ± 0.47), and artifacts (4.75 ± 0.43 vs. 3.83 ± 0.55). The reconstruction time for a single slice was 5.2 hours.

DATA CONCLUSION

The proposed method enables high-resolution cardiac T1 mapping with a short acquisition window and improved image quality.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 1.

Tensor Methods in Biomedical Image Analysis

In the past decade, tensors have become increasingly attractive in different aspects of signal and image processing areas. The main reason is the inefficiency of matrices in representing and analyzing multimodal and multidimensional datasets. Matrices cannot preserve the multidimensional correlation of elements in higher-order datasets and this highly reduces the effectiveness of matrix-based approaches in analyzing multidimensional datasets. Besides this, tensor-based approaches have demonstrated promising performances. These together, encouraged researchers to move from matrices to tensors. Among different signal and image processing applications, analyzing biomedical signals and images is of particular importance. This is due to the need for extracting accurate information from biomedical datasets which directly affects patient's health. In addition, in many cases, several datasets have been recorded simultaneously from a patient. A common example is recording electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) of a patient with schizophrenia. In such a situation, tensors seem to be among the most effective methods for the simultaneous exploitation of two (or more) datasets. Therefore, several tensor-based methods have been developed for analyzing biomedical datasets. Considering this reality, in this paper, we aim to have a comprehensive review on tensor-based methods in biomedical image analysis. The presented study and classification between different methods and applications can show the importance of tensors in biomedical image enhancement and open new ways for future studies.

Non-rigid motion-compensated 3D whole-heart T2 mapping in a hybrid 3T PET-MR system

PURPOSE

Simultaneous PET-MRI improves inflammatory cardiac disease diagnosis. However, challenges persist in respiratory motion and mis-registration between free-breathing 3D PET and 2D breath-held MR images. We propose a free-breathing non-rigid motion-compensated 3D T2 -mapping sequence enabling whole-heart myocardial tissue characterization in a hybrid 3T PET-MR system and provides non-rigid respiratory motion fields to correct also simultaneously acquired PET data.

METHODS

Free-breathing 3D whole-heart T2 -mapping was implemented on a hybrid 3T PET-MRI system. Three datasets were acquired with different T2 -preparation modules (0, 28, 55 ms) using 3-fold undersampled variable-density Cartesian trajectory. Respiratory motion was estimated via virtual 3D image navigators, enabling multi-contrast non-rigid motion-corrected MR reconstruction. T2 -maps were computed using dictionary-matching. Approach was tested in phantom, 8 healthy subjects, 14 MR only and 2 PET-MR patients with suspected cardiac disease and compared with spin echo reference (phantom) and clinical 2D T2 -mapping (in-vivo).

RESULTS

Phantom results show a high correlation (R2  = 0.996) between proposed approach and gold standard 2D T2 mapping. In-vivo 3D T2 -mapping average values in healthy subjects (39.0 ± 1.4 ms) and patients (healthy tissue) (39.1 ± 1.4 ms) agree with conventional 2D T2 -mapping (healthy = 38.6 ± 1.2 ms, patients = 40.3 ± 1.7 ms). Bland-Altman analysis reveals bias of 1.8 ms and 95% limits of agreement (LOA) of -2.4-6 ms for healthy subjects, and bias of 1.3 ms and 95% LOA of -1.9 to 4.6 ms for patients.

CONCLUSION

Validated efficient 3D whole-heart T2 -mapping at hybrid 3T PET-MRI provides myocardial inflammation characterization and non-rigid respiratory motion fields for simultaneous PET data correction. Comparable T2 values were achieved with both 3D and 2D methods. Improved image quality was observed in the PET images after MR-based motion correction.

Accelerated liver water T1 mapping using single-shot continuous inversion-recovery spiral imaging

PURPOSE

Liver T1 mapping techniques typically require long breath holds or long scan time in free-breathing, need correction for B 1 + inhomogeneities and process composite (water and fat) signals. The purpose of this work is to accelerate the multi-slice acquisition of liver water selective T1 (wT1) mapping in a single breath hold, improving the k-space sampling efficiency.

METHODS

The proposed continuous inversion-recovery (IR) Look-Locker methodology combines a single-shot gradient echo spiral readout, Dixon processing and a dictionary-based analysis for liver wT1 mapping at 3 T. The sequence parameters were adapted to obtain short scan times. The influence of fat, B 1 + inhomogeneities and TE on the estimation of T1 was first assessed using simulations. The proposed method was then validated in a phantom and in 10 volunteers, comparing it with MRS and the modified Look-Locker inversion-recovery (MOLLI) method. Finally, the clinical feasibility was investigated by comparing wT1 maps with clinical scans in nine patients.

RESULTS

The phantom results are in good agreement with MRS. The proposed method encodes the IR-curve for the liver wT1 estimation, is minimally sensitive to B 1 + inhomogeneities and acquires one slice in 1.2 s. The volunteer results confirmed the multi-slice capability of the proposed method, acquiring nine slices in a breath hold of 11 s. The present work shows robustness to B 1 + inhomogeneities (wT 1 , No B 1 + = 1.07 wT 1 , B 1 + - 45.63 , R 2 = 0.99 ) , good repeatability (wT 1 , 2 ° = 1 . 0 wT 1 , 1 ° - 2.14 , R 2 = 0.96 ) and is in better agreement with MRS (wT 1 = 0.92 wT 1 MRS + 103.28 , R 2 = 0.38 ) than is MOLLI (wT 1 MOLLI = 0.76 wT 1 MRS + 254.43 , R 2 = 0.44 ) . The wT1 maps in patients captured diverse lesions, thus showing their clinical feasibility.

CONCLUSION

A single-shot spiral acquisition can be combined with a continuous IR Look-Locker method to perform rapid repeatable multi-slice liver water T1 mapping at a rate of 1.2 s per slice without a B 1 + map. The proposed method is suitable for nine-slice liver clinical applications acquired in a single breath hold of 11 s.

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.

mcLARO: Multi-contrast learned acquisition and reconstruction optimization for simultaneous quantitative multi-parametric mapping

PURPOSE

To develop a method for rapid sub-millimeter T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and QSM mapping in a single scan using multi-contrast learned acquisition and reconstruction optimization (mcLARO).

METHODS

A pulse sequence was developed by interleaving inversion recovery and T2 magnetization preparations and single-echo and multi-echo gradient echo acquisitions, which sensitized k-space data to T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and magnetic susceptibility. The proposed mcLARO optimized both the multi-contrast k-space under-sampling pattern and image reconstruction based on image feature fusion using a deep learning framework. The proposed mcLARO method withR = 8 $$ R=8 $$ under-sampling was validated in a retrospective ablation study and compared with other deep learning reconstruction methods, including MoDL and Wave-MoDL, using fully sampled data as reference. Various under-sampling ratios in mcLARO were investigated. mcLARO was also evaluated in a prospective study using separately acquired conventionally sampled quantitative maps as reference standard.

RESULTS

The retrospective ablation study showed improved image sharpness of mcLARO compared to the baseline network without the multi-contrast sampling pattern optimization or image feature fusion module. The under-sampling ratio experiment showed that higher under-sampling ratios resulted in blurrier images and lower precision of quantitative values. The prospective study showed that small or negligible bias and narrow 95% limits of agreement on regional T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and QSM values by mcLARO (5:39 mins) compared to reference scans (40:03 mins in total). mcLARO outperformed MoDL and Wave-MoDL in terms of image sharpness, noise suppression, and artifact removal.

CONCLUSION

mcLARO enabled fast sub-millimeter T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and QSM mapping in a single scan.

Highly efficient free-breathing 3D whole-heart imaging in 3-min: single center study in adults with congenital heart disease

BACKGROUND

Three dimensional, whole-heart (3DWH) MRI is an established non-invasive imaging modality in patients with congenital heart disease (CHD) for the diagnosis of cardiovascular morphology and for clinical decision making. Current techniques utilise diaphragmatic navigation (dNAV) for respiratory motion correction and gating and are frequently limited by long acquisition times. This study proposes and evaluates the diagnostic performance of a respiratory gating-free framework, which considers respiratory image-based navigation (iNAV), and highly accelerated variable density Cartesian sampling in concert with non-rigid motion correction and low-rank patch-based denoising (iNAV-3DWH-PROST). The method is compared to the clinical dNAV-3DWH sequence in adult patients with CHD.

METHODS

In this prospective single center study, adult patients with CHD who underwent the clinical dNAV-3DWH MRI were also scanned with the iNAV-3DWH-PROST. Diagnostic confidence (4-point Likert scale) and diagnostic accuracy for common cardiovascular lesions was assessed by three readers. Scan times and diagnostic confidence were compared using the Wilcoxon-signed rank test. Co-axial vascular dimensions at three anatomic landmarks were measured, and agreement between the research and the corresponding clinical sequence was assessed with Bland-Altman analysis.

RESULTS

The study included 60 participants (mean age ± [SD]: 33 ± 14 years; 36 men). The mean acquisition time of iNAV-3DWH-PROST was significantly lower compared with the conventional clinical sequence (3.1 ± 0.9 min vs 13.9 ± 3.9 min, p < 0.0001). Diagnostic confidence was higher for the iNAV-3DWH-PROST sequence compared with the clinical sequence (3.9 ± 0.2 vs 3.4 ± 0.8, p < 0.001), however there was no significant difference in diagnostic accuracy. Narrow limits of agreement and mean bias less than 0.08 cm were found between the research and the clinical vascular measurements.

CONCLUSIONS

The iNAV-3DWH-PROST framework provides efficient, high quality and robust 3D whole-heart imaging in significantly shorter scan time compared to the standard clinical sequence.

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.

AI-assisted compressed sensing and parallel imaging sequences for MRI of patients with nasopharyngeal carcinoma: comparison of their capabilities in terms of examination time and image quality

OBJECTIVE

To compare examination time and image quality between artificial intelligence (AI)-assisted compressed sensing (ACS) technique and parallel imaging (PI) technique in MRI of patients with nasopharyngeal carcinoma (NPC).

METHODS

Sixty-six patients with pathologically confirmed NPC underwent nasopharynx and neck examination using a 3.0-T MRI system. Transverse T2-weighted fast spin-echo (FSE) sequence, transverse T1-weighted FSE sequence, post-contrast transverse T1-weighted FSE sequence, and post-contrast coronal T1-weighted FSE were obtained by both ACS and PI techniques, respectively. The signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and duration of scanning of both sets of images analyzed by ACS and PI techniques were compared. The images from the ACS and PI techniques were scored for lesion detection, margin sharpness of lesions, artifacts, and overall image quality using the 5-point Likert scale.

RESULTS

The examination time with ACS technique was significantly shorter than that with PI technique (p < 0.0001). The comparison of SNR and CNR showed that ACS technique was significantly superior with PI technique (p < 0.005). Qualitative image analysis showed that the scores of lesion detection, margin sharpness of lesions, artifacts, and overall image quality were higher in the ACS sequences than those in the PI sequences (p < 0.0001). Inter-observer agreement was evaluated for all qualitative indicators for each method, in which the results showed satisfactory-to-excellent agreement (p < 0.0001).

CONCLUSION

Compared with the PI technique, the ACS technique for MR examination of NPC can not only shorten scanning time but also improve image quality.

CLINICAL RELEVANCE STATEMENT

The artificial intelligence (AI)-assisted compressed sensing (ACS) technique shortens examination time for patients with nasopharyngeal carcinoma, while improving the image quality and examination success rate, which will benefit more patients.

KEY POINTS

• Compared with the parallel imaging (PI) technique, the artificial intelligence (AI)-assisted compressed sensing (ACS) technique not only reduced examination time, but also improved image quality. • Artificial intelligence (AI)-assisted compressed sensing (ACS) pulls the state-of-the-art deep learning technique into the reconstruction procedure and helps find an optimal balance of imaging speed and image quality.

Simultaneous Highly Efficient Contrast-Free Lumen and Vessel Wall MR Imaging for Anatomical Assessment of Aortic Disease

BACKGROUND

Bright-blood lumen and black-blood vessel wall imaging are required for the comprehensive assessment of aortic disease. These images are usually acquired separately, resulting in long examinations and potential misregistration between images.

PURPOSE

To characterize the performance of an accelerated and respiratory motion-compensated three-dimensional (3D) cardiac MRI technique for simultaneous contrast-free aortic lumen and vessel wall imaging with an interleaved T2 and inversion recovery prepared sequence (iT2Prep-BOOST).

STUDY TYPE

Prospective.

POPULATION

A total of 30 consecutive patients with aortopathy referred for a clinically indicated cardiac MRI examination (9 females, mean age ± standard deviation: 32 ± 12 years).

FIELD STRENGTH/SEQUENCE

1.5-T; bright-blood MR angiography (diaphragmatic navigator-gated T2-prepared 3D balanced steady-state free precession [bSSFP], T2Prep-bSSFP), breath-held black-blood two-dimensional (2D) half acquisition single-shot turbo spin echo (HASTE), and 3D bSSFP iT2Prep-BOOST.

ASSESSMENT

iT2Prep-BOOST bright-blood images were compared to T2prep-bSSFP images in terms of aortic vessel dimensions, lumen-to-myocardium contrast ratio (CR), and image quality (diagnostic confidence, vessel sharpness and presence of artifacts, assessed by three cardiologists on a 4-point scale, 1: nondiagnostic to 4: excellent). The iT2Prep-BOOST black-blood images were compared to 2D HASTE images for quantification of wall thickness. A visual comparison between computed tomography (CT) and iT2Prep-BOOST was performed in a patient with chronic aortic dissection.

STATISTICAL TESTS

Paired t-tests, Wilcoxon signed-rank tests, intraclass correlation coefficient (ICC), Bland-Altman analysis. A P value < 0.05 was considered statistically significant.

RESULTS

Bright-blood iT2Prep-BOOST resulted in significantly improved image quality (mean ± standard deviation 3.8 ± 0.5 vs. 3.3 ± 0.8) and CR (2.9 ± 0.8 vs. 1.8 ± 0.5) compared with T2Prep-bSSFP, with a shorter scan time (7.8 ± 1.7 minutes vs. 12.9 ± 3.4 minutes) while providing a complementary 3D black-blood image. Aortic lumen diameter and vessel wall thickness measurements in bright-blood and black-blood images were in good agreement with T2Prep-bSSFP and HASTE images (<0.02 cm and <0.005 cm bias, respectively) and good intrareader (ICC > 0.96) and interreader (ICC > 0.94) agreement was observed for all measurements.

DATA CONCLUSION

iT2Prep-BOOST might enable time-efficient simultaneous bright- and black-blood aortic imaging, with improved image quality compared to T2Prep-bSSFP and HASTE imaging, and comparable measurements for aortic wall and lumen dimensions.

EVIDENCE LEVEL

2.

TECHNICAL EFFICACY

Stage 2.

Motion-compensated T1 mapping in cardiovascular magnetic resonance imaging: a technical review

T 1 mapping is becoming a staple magnetic resonance imaging method for diagnosing myocardial diseases such as ischemic cardiomyopathy, hypertrophic cardiomyopathy, myocarditis, and more. Clinically, most T 1 mapping sequences acquire a single slice at a single cardiac phase across a 10 to 15-heartbeat breath-hold, with one to three slices acquired in total. This leaves opportunities for improving patient comfort and information density by acquiring data across multiple cardiac phases in free-running acquisitions and across multiple respiratory phases in free-breathing acquisitions. Scanning in the presence of cardiac and respiratory motion requires more complex motion characterization and compensation. Most clinical mapping sequences use 2D single-slice acquisitions; however newer techniques allow for motion-compensated reconstructions in three dimensions and beyond. To further address confounding factors and improve measurement accuracy, T 1 maps can be acquired jointly with other quantitative parameters such as T 2 , T 2 ∗ , fat fraction, and more. These multiparametric acquisitions allow for constrained reconstruction approaches that isolate contributions to T 1 from other motion and relaxation mechanisms. In this review, we examine the state of the literature in motion-corrected and motion-resolved T 1 mapping, with potential future directions for further technical development and clinical translation.

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.

Improved accuracy and precision of fat-suppressed isotropic 3D T2 mapping MRI of the knee with dictionary fitting and patch-based denoising

PURPOSE

To develop an isotropic three-dimensional (3D) T2 mapping technique for the quantitative assessment of the composition of knee cartilage with high accuracy and precision.

METHODS

A T2-prepared water-selective isotropic 3D gradient-echo pulse sequence was used to generate four images at 3 T. These were used for three T2 map reconstructions: standard images with an analytical T2 fit (AnT2Fit); standard images with a dictionary-based T2 fit (DictT2Fit); and patch-based-denoised images with a dictionary-based T2 fit (DenDictT2Fit). The accuracy of the three techniques was first optimized in a phantom study against spin-echo imaging, after which knee cartilage T2 values and coefficients of variation (CoV) were assessed in ten subjects in order to establish accuracy and precision in vivo. Data given as mean ± standard deviation.

RESULTS

After optimization in the phantom, whole-knee cartilage T2 values of the healthy volunteers were 26.6 ± 1.6 ms (AnT2Fit), 42.8 ± 1.8 ms (DictT2Fit, p < 0.001 versus AnT2Fit), and 40.4 ± 1.7 ms (DenDictT2Fit, p = 0.009 versus DictT2Fit). The whole-knee T2 CoV reduced from 51.5% ± 5.6% to 30.5 ± 2.4 and finally to 13.1 ± 1.3%, respectively (p < 0.001 between all). The DictT2Fit improved the data reconstruction time: 48.7 ± 11.3 min (AnT2Fit) versus 7.3 ± 0.7 min (DictT2Fit, p < 0.001). Very small focal lesions were observed in maps generated with DenDictT2Fit.

CONCLUSIONS

Improved accuracy and precision for isotropic 3D T2 mapping of knee cartilage were demonstrated by using patch-based image denoising and dictionary-based reconstruction.

KEY POINTS

• Dictionary T2 fitting improves the accuracy of three-dimensional (3D) knee T2 mapping. • Patch-based denoising results in high precision in 3D knee T2 mapping. • Isotropic 3D knee T2 mapping enables the visualization of small anatomical details.

Single-heartbeat cardiac cine imaging via jointly regularized nonrigid motion-corrected reconstruction

The aim of the current study was to develop a novel approach for 2D breath-hold cardiac cine imaging from a single heartbeat, by combining cardiac motion-corrected reconstructions and nonrigidly aligned patch-based regularization. Conventional cardiac cine imaging is obtained via motion-resolved reconstructions of data acquired over multiple heartbeats. Here, we achieve single-heartbeat cine imaging by incorporating nonrigid cardiac motion correction into the reconstruction of each cardiac phase, in conjunction with a motion-aligned patch-based regularization. The proposed Motion-Corrected CINE (MC-CINE) incorporates all acquired data into the reconstruction of each (motion-corrected) cardiac phase, resulting in a better posed problem than motion-resolved approaches. MC-CINE was compared with iterative sensitivity encoding (itSENSE) and Extra-Dimensional Golden Angle Radial Sparse Parallel (XD-GRASP) in 14 healthy subjects in terms of image sharpness, reader scoring (range: 1-5) and reader ranking (range: 1-9) of image quality, and single-slice left ventricular assessment. MC-CINE was significantly superior to both itSENSE and XD-GRASP using 20 heartbeats, two heartbeats, and one heartbeat. Iterative SENSE, XD-GRASP, and MC-CINE achieved a sharpness of 74%, 74%, and 82% using 20 heartbeats, and 53%, 66%, and 82% with one heartbeat, respectively. The corresponding results for reader scoring were 4.0, 4.7, and 4.9 with 20 heartbeats, and 1.1, 3.0, and 3.9 with one heartbeat. The corresponding results for reader ranking were 5.3, 7.3, and 8.6 with 20 heartbeats, and 1.0, 3.2, and 5.4 with one heartbeat. MC-CINE using a single heartbeat presented nonsignificant differences in image quality to itSENSE with 20 heartbeats. MC-CINE and XD-GRASP at one heartbeat both presented a nonsignificant negative bias of less than 2% in ejection fraction relative to the reference itSENSE. It was concluded that the proposed MC-CINE significantly improves image quality relative to itSENSE and XD-GRASP, enabling 2D cine from a single heartbeat.

Latest Advances in Image Acceleration: All Dimensions are Fair Game

Magnetic resonance imaging (MRI) is a versatile modality that can generate high-resolution images with a variety of tissue contrasts. However, MRI is a slow technique and requires long acquisition times, which increase with higher temporal and spatial resolution and/or when multiple contrasts and large volumetric coverage is required. In order to speedup MR data acquisition, several approaches have been introduced in the literature. Most of these techniques acquire less data than required and exploit intrinsic redundancies in the MR images to recover the information that was not sampled. This article presents a review of MR acquisition and reconstruction methods that have exploited redundancies in the temporal, spatial, and contrast/parametric dimensions to accelerate image data acquisition, focusing on cardiac and abdominal MR imaging applications. The review describes how each of these dimensions has been separately exploited for speeding up MR acquisition to then discuss more advanced techniques where multiple dimensions are exploited together for further reducing scan times. Finally, future directions for multidimensional image acceleration and remaining technical challenges are discussed. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: 1.

Parallel non-Cartesian spatial-temporal dictionary learning neural networks (stDLNN) for accelerating 4D-MRI

Dynamic magnetic resonance imaging (MRI) acquisitions are relatively slow due to physical and physiological limitations. The spatial-temporal dictionary learning (DL) approach accelerates dynamic MRI by learning spatial-temporal correlations, but the regularization parameters need to be manually adjusted, the performance at high acceleration rate is limited, and the reconstruction can be time-consuming. Deep learning techniques have shown good performance in accelerating MRI due to the powerful representational capabilities of neural networks. In this work, we propose a parallel non-Cartesian spatial-temporal dictionary learning neural networks (stDLNN) framework that combines dictionary learning with deep learning algorithms and utilizes the spatial-temporal prior information of dynamic MRI data to achieve better reconstruction quality and efficiency. The coefficient estimation modules (CEM) are designed in the framework to adaptively adjust the regularization coefficients. Experimental results show that combining dictionary learning with deep neural networks and using spatial-temporal dictionaries can obviously improve the image quality and computational efficiency compared with the state-of-the-art non-Cartesian imaging methods for accelerating the 4D-MRI especially at high acceleration rate.

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.

MR multitasking-based dynamic imaging for cerebrovascular evaluation (MT-DICE): Simultaneous quantification of permeability and leakage-insensitive perfusion by dynamic T 1 / T 2 * mapping

PURPOSE

To develop an MR multitasking-based dynamic imaging for cerebrovascular evaluation (MT-DICE) technique for simultaneous quantification of permeability and leakage-insensitive perfusion with a single-dose contrast injection.

METHODS

MT-DICE builds on a saturation-recovery prepared multi-echo fast low-angle shot sequence. The k-space is randomly sampled for 7.6 min, with single-dose contrast agent injected 1.5 min into the scan. MR multitasking is used to model the data into six dimensions, including three spatial dimensions for whole-brain coverage, a saturation-recovery time dimension, and a TE dimension for dynamicT 1 $$ {\mathrm{T}}_1 $$ andT 2 * $$ {\mathrm{T}}_2^{\ast } $$ quantification, respectively, and a contrast dynamics dimension for capturing contrast kinetics. The derived pixel-wiseT 1 / T 2 * $$ {\mathrm{T}}_1/{\mathrm{T}}_2^{\ast } $$ time series are converted into contrast concentration-time curves for calculation of kinetic metrics. The technique was assessed for its agreement with reference methods inT 1 $$ {\mathrm{T}}_1 $$ andT 2 * $$ {\mathrm{T}}_2^{\ast } $$ measurements in eight healthy subjects and, in three of them, inter-session repeatability of permeability and leakage-insensitive perfusion parameters. Its feasibility was also demonstrated in four patients with brain tumors.

RESULTS

MT-DICET 1 / T 2 * $$ {\mathrm{T}}_1/{\mathrm{T}}_2^{\ast } $$ values of normal gray matter and white matter were in excellent agreement with reference values (intraclass correlation coefficients = 0.860/0.962 for gray matter and 0.925/0.975 for white matter ). Both permeability and perfusion parameters demonstrated good to excellent intersession agreement with the lowest intraclass correlation coefficients at 0.694. Contrast kinetic parameters in all healthy subjects and patients were within the literature range.

CONCLUSION

Based on dynamicT 1 / T 2 * $$ {\mathrm{T}}_1/{\mathrm{T}}_2^{\ast } $$ mapping, MT-DICE allows for simultaneous quantification of permeability and leakage-insensitive perfusion metrics with a single-dose contrast injection.

patchIT: A Multipurpose Patch Creation Tool for Image Processing Applications

Patch-based approaches in image processing are often preferable to working with the entire image. They provide an alternative representation of the image as a set of partial local sub-images (patches) which is a vital preprocessing step in many image processing applications. In this paper, a new software tool called patchIT is presented, providing an integrated framework suitable for the systematic and automatized extraction of patches from images based on user-defined geometrical and spatial criteria. Patches can be extracted in both a sliding and random manner and can be exported either as images, MATLAB .mat files, or raw text files. The proposed tool offers further functionality, including masking operations that act as spatial filters, identifying candidate patch areas, as well as geometric transformations by applying patch value indexing. It also efficiently handles issues that arise in large-scale patch processing scenarios in terms of memory and time requirements. In addition, a use case in cartographic research is presented that utilizes patchIT for map evaluation purposes based on a visual heterogeneity indicator. The tool supports all common image file formats and efficiently processes bitonal, grayscale, color, and multispectral images. PatchIT is freely available to the scientific community under the third version of GNU General Public License (GPL v3) on the GitHub platform.

Simultaneous multi-parametric acquisition and reconstruction techniques in cardiac magnetic resonance imaging: Basic concepts and status of clinical development

Simultaneous multi-parametric acquisition and reconstruction techniques (SMART) are gaining attention for their potential to overcome some of cardiovascular magnetic resonance imaging's (CMR) clinical limitations. The major advantages of SMART lie within their ability to simultaneously capture multiple "features" such as cardiac motion, respiratory motion, T1/T2 relaxation. This review aims to summarize the overarching theory of SMART, describing key concepts that many of these techniques share to produce co-registered, high quality CMR images in less time and with less requirements for specialized personnel. Further, this review provides an overview of the recent developments in the field of SMART by describing how they work, the parameters they can acquire, their status of clinical testing and validation, and by providing examples for how their use can improve the current state of clinical CMR workflows. Many of the SMART are in early phases of development and testing, thus larger scale, controlled trials are needed to evaluate their use in clinical setting and with different cardiac pathologies.

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.

Repeatability of simultaneous 3D 1H MRF/23Na MRI in brain at 7 T

Proton MRI can provide detailed morphological images, but it reveals little information about cell homeostasis. On the other hand, sodium MRI can provide metabolic information but cannot resolve fine structures. The complementary nature of proton and sodium MRI raises the prospect of their combined use in a single experiment. In this work, we assessed the repeatability of normalized proton density (PD), T₁, T₂, and normalized sodium density-weighted quantification measured with simultaneous 3D ¹H MRF/²³Na MRI in the brain at 7 T, from ten healthy volunteers who were scanned three times each. The coefficients of variation (CV) and the intra-class correlation (ICC) were calculated for the mean and standard deviation (SD) of these 4 parameters in grey matter, white matter, and cerebrospinal fluid. As result, the CVs were lower than 3.3% for the mean values and lower than 6.9% for the SD values. The ICCs were higher than 0.61 in all 24 measurements. We conclude that the measurements of normalized PD, T₁, T₂, and normalized sodium density-weighted from simultaneous 3D ¹H MRF/²³Na MRI in the brain at 7 T showed high repeatability. We estimate that changes > 6.6% (> 2 CVs) in mean values of both ¹H and ²³Na metrics could be detectable with this method.

Simultaneous 3D acquisition of 1 H MRF and 23 Na MRI

PURPOSE

To develop a 3D MR technique to simultaneously acquire proton multiparametric maps (T₁ , T₂ , and proton density) and sodium density weighted images over the whole brain.

METHODS

We implemented a 3D stack-of-stars MR pulse sequence which consists of interleaved proton (¹ H) and sodium (²³ Na) excitations, tailored slice encoding gradients that can encode the same slice for both nuclei, and simultaneous readout with different radial trajectories (¹ H, full-radial; ²³ Na, center-out radial). The receive chain of our 7T scanner was modified to enable simultaneous acquisition of ¹ H and ²³ Na signal. A heuristically optimized flip angle train was implemented for proton MR fingerprinting (MRF). The SNR and the accuracy of proton T₁ and T₂ were evaluated in phantoms. Finally, in vivo application of the method was demonstrated in five healthy subjects.

RESULTS

The SNR for the simultaneous measurement was almost identical to that for the single-nucleus measurements (<2% change). The proton T₁ and T₂ maps remained similar to the results from a reference 2D MRF technique (normalized RMS error in T₁ ≈ 4.2% and T₂ ≈ 11.3%). Measurements in healthy subjects corroborated these results and demonstrated the feasibility of our method for in vivo application. The in vivo T₁ values measured using our method were lower than the results measured by other conventional techniques.

CONCLUSIONS

With the 3D simultaneous implementation, we were able to acquire sodium and proton density weighted images in addition to proton T₁ , T₂ , and B1+ from ¹ H MRF that covers the whole brain volume within 21 min.

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.

Cardiac MR: From Theory to Practice

Cardiovascular disease (CVD) is the leading single cause of morbidity and mortality, causing over 17. 9 million deaths worldwide per year with associated costs of over $800 billion. Improving prevention, diagnosis, and treatment of CVD is therefore a global priority. Cardiovascular magnetic resonance (CMR) has emerged as a clinically important technique for the assessment of cardiovascular anatomy, function, perfusion, and viability. However, diversity and complexity of imaging, reconstruction and analysis methods pose some limitations to the widespread use of CMR. Especially in view of recent developments in the field of machine learning that provide novel solutions to address existing problems, it is necessary to bridge the gap between the clinical and scientific communities. This review covers five essential aspects of CMR to provide a comprehensive overview ranging from CVDs to CMR pulse sequence design, acquisition protocols, motion handling, image reconstruction and quantitative analysis of the obtained data. (1) The basic MR physics of CMR is introduced. Basic pulse sequence building blocks that are commonly used in CMR imaging are presented. Sequences containing these building blocks are formed for parametric mapping and functional imaging techniques. Commonly perceived artifacts and potential countermeasures are discussed for these methods. (2) CMR methods for identifying CVDs are illustrated. Basic anatomy and functional processes are described to understand the cardiac pathologies and how they can be captured by CMR imaging. (3) The planning and conduct of a complete CMR exam which is targeted for the respective pathology is shown. Building blocks are illustrated to create an efficient and patient-centered workflow. Further strategies to cope with challenging patients are discussed. (4) Imaging acceleration and reconstruction techniques are presented that enable acquisition of spatial, temporal, and parametric dynamics of the cardiac cycle. The handling of respiratory and cardiac motion strategies as well as their integration into the reconstruction processes is showcased. (5) Recent advances on deep learning-based reconstructions for this purpose are summarized. Furthermore, an overview of novel deep learning image segmentation and analysis methods is provided with a focus on automatic, fast and reliable extraction of biomarkers and parameters of clinical relevance.

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.

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.

Simultaneous T1 , T2 , and T1ρ cardiac magnetic resonance fingerprinting for contrast agent-free myocardial tissue characterization

PURPOSE

To develop a simultaneous T₁ , T₂ , and T_1ρ cardiac magnetic resonance fingerprinting (MRF) approach to enable comprehensive contrast agent-free myocardial tissue characterization in a single breath-hold scan.

METHODS

A 2D gradient-echo electrocardiogram-triggered cardiac MRF sequence with low flip angles, varying magnetization preparation, and spiral trajectory was acquired at 1.5 T to encode T₁ , T₂ , and T_1⍴ simultaneously. The MRF images were reconstructed using low-rank inversion, regularized with a multicontrast patch-based higher-order reconstruction. Parametric maps were generated and matched in the singular value domain to extended phase graph-based dictionaries. The proposed approach was tested in phantoms and 10 healthy subjects and compared against conventional methods in terms of coefficients of determination and best fits for the phantom study, and in terms of Bland-Altman agreement, average values and coefficient of variation of T₁ , T₂ , and T_1⍴ for the healthy subjects study.

RESULTS

The T₁ , T₂ , and T_1⍴ MRF values showed excellent correlation with conventional spin-echo and clinical mapping methods in phantom studies (r² > 0.97). Measured MRF values in myocardial tissue (mean ± SD) were 1133 ± 33 ms, 38.8 ± 3.5 ms, and 52.0 ± 4.0 ms for T₁ , T₂ and T1_⍴ , respectively, against 1053 ± 47 ms, 50.4 ± 3.9 ms, and 55.9 ± 3.3 ms for T₁ modified Look-Locker inversion imaging, T₂ gradient and spin echo, and T_1⍴ turbo field echo, respectively.

CONCLUSION

A cardiac MRF approach for simultaneous quantification of myocardial T₁ , T₂ , and T_1ρ in a single breath-hold MR scan of about 16 seconds has been proposed. The approach has been investigated in phantoms and healthy subjects showing good agreement with reference spin echo measurements and conventional clinical maps.

Simultaneous comprehensive liver T1 , T2 , T 2 ∗ , T1ρ , and fat fraction characterization with MR fingerprinting

PURPOSE

To develop a novel simultaneous co-registered T₁ , T₂ , T 2 ∗ , T_1ρ , and fat fraction abdominal MR fingerprinting (MRF) approach for fully comprehensive liver-tissue characterization in a single breath-hold scan.

METHODS

A gradient-echo liver MRF sequence with low fixed flip angle, multi-echo radial readout, and varying magnetization preparation pulses for multiparametric encoding is performed at 1.5 T. The T 2 ∗ and fat fraction are estimated from a graph/cut water/fat separation method using a six-peak fat model. Water/fat singular images obtained are then matched to an MRF dictionary, estimating water-specific T₁ , T₂ , and T_1ρ . The proposed approach was tested in phantoms and 10 healthy subjects and compared against conventional sequences.

RESULTS

For the phantom studies, linear fits show excellent coefficients of determination (r² > 0.9) for every parametric map. For in vivo studies, the average values measured within regions of interest drawn on liver, spleen, muscle, and fat are statistically different from the reference scans (p < 0.05) for T₁ , T₂ , and T_1⍴ but not for T 2 ∗ and fat fraction, whereas correlation between MRF and reference scans is excellent for each parameter (r² > 0.92 for every parameter).

CONCLUSION

The proposed multi-echo inversion-recovery, T₂ , and T_1⍴ prepared liver MRF sequence presented in this work allows for quantitative T₁ , T₂ , T 2 ∗ , T_1⍴ , and fat fraction liver-tissue characterization in a single breath-hold scan of 18 seconds. The approach showed good agreement and correlation with respect to reference clinical maps.

Respiratory Motion-Registered Isotropic Whole-Heart T2 Mapping in Patients With Acute Non-ischemic Myocardial Injury

Background: T₂ mapping is a magnetic resonance imaging technique that can be used to detect myocardial edema and inflammation. However, the focal nature of myocardial inflammation may render conventional 2D approaches suboptimal and make whole-heart isotropic 3D mapping desirable. While self-navigated 3D radial T₂ mapping has been demonstrated to work well at a magnetic field strength of 3T, it results in too noisy maps at 1.5T. We therefore implemented a novel respiratory motion-resolved compressed-sensing reconstruction in order to improve the 3D T₂ mapping precision and accuracy at 1.5T, and tested this in a heterogeneous patient cohort. Materials and Methods: Nine healthy volunteers and 25 consecutive patients with suspected acute non-ischemic myocardial injury (sarcoidosis, n = 19; systemic sclerosis, n = 2; acute graft rejection, n = 2, and myocarditis, n = 2) were included. The free-breathing T₂ maps were acquired as three ECG-triggered T₂-prepared 3D radial volumes. A respiratory motion-resolved reconstruction was followed by image registration of the respiratory states and pixel-wise T₂ mapping. The resulting 3D maps were compared to routine 2D T₂ maps. The T₂ values of segments with and without late gadolinium enhancement (LGE) were compared in patients. Results: In the healthy volunteers, the myocardial T₂ values obtained with the 2D and 3D techniques were similar (45.8 ± 1.8 vs. 46.8 ± 2.9 ms, respectively; P = 0.33). Conversely, in patients, T₂ values did differ between 2D (46.7 ± 3.6 ms) and 3D techniques (50.1 ± 4.2 ms, P = 0.004). Moreover, with the 2D technique, T₂ values of the LGE-positive segments were similar to those of the LGE-negative segments (T_2LGE-= 46.2 ± 3.7 vs. T_2LGE+ = 47.6 ± 4.1 ms; P = 0.49), whereas the 3D technique did show a significant difference (T_2LGE- = 49.3 ± 6.7 vs. T_2LGE+ = 52.6 ± 8.7 ms, P = 0.006). Conclusion: Respiratory motion-registered 3D radial imaging at 1.5T led to accurate isotropic 3D whole-heart T₂ maps, both in the healthy volunteers and in a small patient cohort with suspected non-ischemic myocardial injury. Significantly higher T₂ values were found in patients as compared to controls in 3D but not in 2D, suggestive of the technique's potential to increase the sensitivity of CMR at earlier stages of disease. Further study will be needed to demonstrate its accuracy.

Generalized low-rank nonrigid motion-corrected reconstruction for MR fingerprinting

PURPOSE

Develop a novel low-rank motion-corrected (LRMC) reconstruction for nonrigid motion-corrected MR fingerprinting (MRF).

METHODS

Generalized motion-corrected (MC) reconstructions have been developed for steady-state imaging. Here we extend this framework to enable nonrigid MC for transient imaging applications with varying contrast, such as MRF. This is achieved by integrating low-rank dictionary-based compression into the generalized MC model to reconstruct MC singular images, reducing motion artifacts in the resulting parametric maps. The proposed LRMC reconstruction was applied for cardiac motion correction in 2D myocardial MRF (T₁ and T₂ ) with extended cardiac acquisition window (~450 ms) and for respiratory MC in free-breathing 3D myocardial and 3D liver MRF. Experiments were performed in phantom and 22 healthy subjects. The proposed approach was compared with reference spin echo (phantom) and with 2D electrocardiogram-triggered/breath-hold MOLLI and T₂ gradient-and-spin echo conventional maps (in vivo 2D and 3D myocardial MRF).

RESULTS

Phantom results were in general agreement with reference spin-echo measurements, presenting relative errors of approximately 5.4% and 5.5% for T₁ and short T₂ (<100 ms), respectively. The proposed LRMC MRF reduced residual blurring artifacts with respect to no MC for cardiac or respiratory motion in all cases (2D and 3D myocardial, 3D abdominal). In 2D myocardial MRF, left-ventricle T₁ values were 1150 ± 41 ms for LRMC MRF and 1010 ± 56 ms for MOLLI; T₂ values were 43.8 ± 2.3 ms for LRMC MRF and 49.5 ± 4.5 ms for T₂ gradient and spin echo. Corresponding measurements for 3D myocardial MRF were 1085 ± 30 ms and 1062 ± 29 ms for T₁ , and 43.5 ± 1.9 ms and 51.7 ± 1.7 ms for T₂ . For 3D liver, LRMC MRF measured liver T₁ at 565 ± 44 ms and liver T₂ at 35.4 ± 2.4 ms.

CONCLUSION

The proposed LRMC reconstruction enabled generalized (nonrigid) MC for 2D and 3D MRF, both for cardiac and respiratory motion. The proposed approach reduced motion artifacts in the MRF maps with respect to no motion compensation and achieved good agreement with reference measurements.

Channel Attention Networks for Robust MR Fingerprint Matching

Magnetic Resonance Fingerprinting (MRF) enables simultaneous mapping of multiple tissue parameters such as T1 and T2 relaxation times. The working principle of MRF relies on varying acquisition parameters pseudo-randomly, so that each tissue generates its unique signal evolution during scanning. Even though MRF provides faster scanning, it has disadvantages such as erroneous and slow generation of the corresponding parametric maps, which needs to be improved. Moreover, there is a need for explainable architectures for understanding the guiding signals to generate accurate parametric maps. In this paper, we addressed both of these shortcomings by proposing a novel neural network architecture (CONV-ICA) consisting of a channel-wise attention module and a fully convolutional network. The proposed approach, evaluated over 3 simulated MRF signals, reduces error in the reconstruction of tissue parameters by 8.88% for T1 and 75.44% for T2 with respect to state-of-the-art methods. Another contribution of this study is a new channel selection method: attention-based channel selection. Furthermore, the effect of patch size and temporal frames of MRF signal on channel reduction are analyzed by employing a channel-wise attention.

Systematic Analysis and Review of Magnetic Resonance Imaging (MRI) Reconstruction Techniques

BACKGROUND

Magnetic Resonance Imaging (MRI) plays an important role in the field of medical diagnostic imaging as it poses non-invasive acquisition and high soft-tissue contrast. However, a huge time is needed for the MRI scanning process that results in motion artifacts, degrades image quality, misinterprets the data, and may cause discomfort to the patient. Thus, the main goal of MRI research is to accelerate data acquisition processing without affecting the quality of the image.

INTRODUCTION

This paper presents a survey based on distinct conventional MRI reconstruction methodologies. In addition, a novel MRI reconstruction strategy is proposed based on weighted Compressive Sensing (CS), Penalty-aided minimization function, and Meta-heuristic optimization technique.

METHODS

An illustrative analysis is done concerning adapted methods, datasets used, execution tools, performance measures, and values of evaluation metrics. Moreover, the issues of existing methods and the research gaps considering conventional MRI reconstruction schemes are elaborated to obtain improved contribution for devising significant MRI reconstruction techniques.

RESULTS

The proposed method will reduce conventional aliasing artifact problems, may attain lower Mean Square Error (MSE), higher Peak Signal-to-Noise Ratio (PSNR), and Structural SIMilarity (SSIM) index.

CONCLUSION

The issues of existing methods and the research gaps considering conventional MRI reconstruction schemes are elaborated to devising an improved significant MRI reconstruction technique.

Optimizing multicontrast MRI reconstruction with shareable feature aggregation and selection

This paper proposes a new method for optimizing feature sharing in deep neural network-based, rapid, multicontrast magnetic resonance imaging (MC-MRI). Using the shareable information of MC images for accelerated MC-MRI reconstruction, current algorithms stack the MC images or features without optimizing the sharing protocols, leading to suboptimal reconstruction results. In this paper, we propose a novel feature aggregation and selection scheme in a deep neural network to better leverage the MC features and improve the reconstruction results. First, we propose to extract and use the shareable information by mapping the MC images into multiresolution feature maps with multilevel layers of the neural network. In this way, the extracted features capture complementary image properties, including local patterns from the shallow layers and semantic information from the deep layers. Then, an explicit selection module is designed to compile the extracted features optimally. That is, larger weights are learned to incorporate the constructive, shareable features; and smaller weights are assigned to the unshareable information. We conduct comparative studies on publicly available T2-weighted and T2-weighted fluid attenuated inversion recovery brain images, and the results show that the proposed network consistently outperforms existing algorithms. In addition, the proposed method can recover the images with high fidelity under 16 times acceleration. The ablation studies are conducted to evaluate the effectiveness of the proposed feature aggregation and selection mechanism. The results and the visualization of the weighted features show that the proposed method does effectively improve the usage of the useful features and suppress useless information, leading to overall enhanced reconstruction results. Additionally, the selection module can zero-out repeated and redundant features and improve network efficiency.

Synergistic multi-contrast cardiac magnetic resonance image reconstruction

Cardiac magnetic resonance imaging (CMR) is an important tool for the non-invasive diagnosis of a variety of cardiovascular diseases. Parametric mapping with multi-contrast CMR is able to quantify tissue alterations in myocardial disease and promises to improve patient care. However, magnetic resonance imaging is an inherently slow imaging modality, resulting in long acquisition times for parametric mapping which acquires a series of cardiac images with different contrasts for signal fitting or dictionary matching. Furthermore, extra efforts to deal with respiratory and cardiac motion by triggering and gating further increase the scan time. Several techniques have been developed to speed up CMR acquisitions, which usually acquire less data than that required by the Nyquist-Shannon sampling theorem, followed by regularized reconstruction to mitigate undersampling artefacts. Recent advances in CMR parametric mapping speed up CMR by synergistically exploiting spatial-temporal and contrast redundancies. In this article, we will review the recent developments in multi-contrast CMR image reconstruction for parametric mapping with special focus on low-rank and model-based reconstructions. Deep learning-based multi-contrast reconstruction has recently been proposed in other magnetic resonance applications. These developments will be covered to introduce the general methodology. Current technical limitations and potential future directions are discussed. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.