Accelerated cardiac T1 mapping with recurrent networks and cyclic, model-based loss

MyoFold: Joint myocardial tissue composition and wall motion quantification via a highly folded sequence

PURPOSE

This study aims to develop and evaluate a novel cardiovascular MR sequence, MyoFold, designed for the simultaneous quantifications of myocardial tissue composition and wall motion.

METHODS

MyoFold is designed as a 2D single breathing-holding sequence, integrating joint T1/T2 mapping and cine imaging. The sequence uses a 2-fold accelerated balanced SSFP (bSSFP) for data readout and incorporates electrocardiogram synchronization to align with the cardiac cycle. MyoFold initially acquires six single-shot inversion-recovery images, completed during the diastole of six successive heartbeats. T2 preparation (T2-prep) is applied to introduce T2 weightings for the last three images. Subsequently, over the following six heartbeats, segmented bSSFP is performed for the movie of the entire cardiac cycle, synchronized with an electrocardiogram. A neural network trained using numerical simulations of MyoFold is used for T1 and T2 calculations. MyoFold was validated through phantom and in vivo experiments, with comparisons made against MOLLI, SASHA, T2-prep bSSFP, and the conventional cine.

RESULTS

In phantom studies, MyoFold exhibited a 10% overestimation in T1 measurements, whereas T2 measurements demonstrated high accuracy. In vivo experiments revealed that MyoFold T1 had comparable accuracy to SASHA and precision similar to MOLLI. MyoFold demonstrated good agreement with T2-prep bSSFP in myocardial T2 measurements. No significant differences were observed in the quantification of left-ventricle wall thickness and function between MyoFold and the conventional cine.

CONCLUSION

MyoFold presents as a rapid, simple, and multitasking approach for quantitative cardiovascular MR examinations, offering simultaneous assessment of tissue composition and wall motion. The sequence's multitasking capabilities make it a promising tool for comprehensive cardiac evaluations in clinical settings.

Recurrent neural network-based simultaneous cardiac T1, T2, and T1ρ mapping

The purpose of the current study was to explore the feasibility of training a deep neural network to accelerate the process of generating T1, T2, and T1ρ maps for a recently proposed free-breathing cardiac multiparametric mapping technique, where a recurrent neural network (RNN) was utilized to exploit the temporal correlation among the multicontrast images. The RNN-based model was developed for rapid and accurate T1, T2, and T1ρ estimation. Bloch simulation was performed to simulate a dataset of more than 10 million signals and time correspondences with different noise levels for network training. The proposed RNN-based method was compared with a dictionary-matching method and a conventional mapping method to evaluate the model's effectiveness in phantom and in vivo studies at 3 T, respectively. In phantom studies, the RNN-based method and the dictionary-matching method achieved similar accuracy and precision in T1, T2, and T1ρ estimations. In in vivo studies, the estimated T1, T2, and T1ρ values obtained by the two methods achieved similar accuracy and precision for 10 healthy volunteers (T1: 1228.70 ± 53.80 vs. 1228.34 ± 52.91 ms, p > 0.1; T2: 40.70 ± 2.89 vs. 41.19 ± 2.91 ms, p > 0.1; T1ρ: 45.09 ± 4.47 vs. 45.23 ± 4.65 ms, p > 0.1). The RNN-based method can generate cardiac multiparameter quantitative maps simultaneously in just 2 s, achieving 60-fold acceleration compared with the dictionary-matching method. The RNN-accelerated method offers an almost instantaneous approach for reconstructing accurate T1, T2, and T1ρ maps, being much more efficient than the dictionary-matching method for the free-breathing multiparametric cardiac mapping technique, which may pave the way for inline mapping in clinical applications.

Arterial Input Function (AIF) Correction Using AIF Plus Tissue Inputs with a Bi-LSTM Network

Background: The arterial input function (AIF) is vital for myocardial blood flow quantification in cardiac MRI to indicate the input time-concentration curve of a contrast agent. Inaccurate AIFs can significantly affect perfusion quantification. Purpose: When only saturated and biased AIFs are measured, this work investigates multiple ways of leveraging tissue curve information, including using AIF + tissue curves as inputs and optimizing the loss function for deep neural network training. Methods: Simulated data were generated using a 12-parameter AIF mathematical model for the AIF. Tissue curves were created from true AIFs combined with compartment-model parameters from a random distribution. Using Bloch simulations, a dictionary was constructed for a saturation-recovery 3D radial stack-of-stars sequence, accounting for deviations such as flip angle, T2* effects, and residual longitudinal magnetization after the saturation. A preliminary simulation study established the optimal tissue curve number using a bidirectional long short-term memory (Bi-LSTM) network with just AIF loss. Further optimization of the loss function involves comparing just AIF loss, AIF with compartment-model-based parameter loss, and AIF with compartment-model tissue loss. The optimized network was examined with both simulation and hybrid data, which included in vivo 3D stack-of-star datasets for testing. The AIF peak value accuracy and ktrans results were assessed. Results: Increasing the number of tissue curves can be beneficial when added tissue curves can provide extra information. Using just the AIF loss outperforms the other two proposed losses, including adding either a compartment-model-based tissue loss or a compartment-model parameter loss to the AIF loss. With the simulated data, the Bi-LSTM network reduced the AIF peak error from -23.6 ± 24.4% of the AIF using the dictionary method to 0.2 ± 7.2% (AIF input only) and 0.3 ± 2.5% (AIF + ten tissue curve inputs) of the network AIF. The corresponding ktrans error was reduced from -13.5 ± 8.8% to -0.6 ± 6.6% and 0.3 ± 2.1%. With the hybrid data (simulated data for training; in vivo data for testing), the AIF peak error was 15.0 ± 5.3% and the corresponding ktrans error was 20.7 ± 11.6% for the AIF using the dictionary method. The hybrid data revealed that using the AIF + tissue inputs reduced errors, with peak error (1.3 ± 11.1%) and ktrans error (-2.4 ± 6.7%). Conclusions: Integrating tissue curves with AIF curves into network inputs improves the precision of AI-driven AIF corrections. This result was seen both with simulated data and with applying the network trained only on simulated data to a limited in vivo test dataset.

Virtual MOLLI Target: Generative Adversarial Networks Toward Improved Motion Correction in MRI Myocardial T1 Mapping

BACKGROUND

The modified Look-Locker inversion recovery (MOLLI) sequence is commonly used for myocardial T1 mapping. However, it acquires images with different inversion times, which causes difficulty in motion correction for respiratory-induced misregistration to a given target image.

HYPOTHESIS

Using a generative adversarial network (GAN) to produce virtual MOLLI images with consistent heart positions can reduce respiratory-induced misregistration of MOLLI datasets.

STUDY TYPE

Retrospective.

POPULATION

1071 MOLLI datasets from 392 human participants.

FIELD STRENGTH/SEQUENCE

Modified Look-Locker inversion recovery sequence at 3 T.

ASSESSMENT

A GAN model with a single inversion time image as input was trained to generate virtual MOLLI target (VMT) images at different inversion times which were subsequently used in an image registration algorithm. Four VMT models were investigated and the best performing model compared with the standard vendor-provided motion correction (MOCO) technique.

STATISTICAL TESTS

The effectiveness of the motion correction technique was assessed using the fitting quality index (FQI), mutual information (MI), and Dice coefficients of motion-corrected images, plus subjective quality evaluation of T1 maps by three independent readers using Likert score. Wilcoxon signed-rank test with Bonferroni correction for multiple comparison. Significance levels were defined as P < 0.01 for highly significant differences and P < 0.05 for significant differences.

RESULTS

The best performing VMT model with iterative registration demonstrated significantly better performance (FQI 0.88 ± 0.03, MI 1.78 ± 0.20, Dice 0.84 ± 0.23, quality score 2.26 ± 0.95) compared to other approaches, including the vendor-provided MOCO method (FQI 0.86 ± 0.04, MI 1.69 ± 0.25, Dice 0.80 ± 0.27, quality score 2.16 ± 1.01).

DATA CONCLUSION

Our GAN model generating VMT images improved motion correction, which may assist reliable T1 mapping in the presence of respiratory motion. Its robust performance, even with considerable respiratory-induced heart displacements, may be beneficial for patients with difficulties in breath-holding.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 1.

The influence of left bundle branch block on myocardial T1 mapping

Tissue characterisation using T1 mapping has become an established magnetic resonance imaging (MRI) technique to detect myocardial diseases. This retrospective study aimed to determine the influence of left bundle branch block (LBBB) on T1 mapping at 1.5 T. Datasets of 36 patients with LBBB and 27 healthy controls with T1 mapping (Modified Look-Locker inversion-recovery (MOLLI), 5(3)3 sampling) were included. T1 relaxation times were determined on mid-cavity short-axis images. R2 maps were generated as a pixel-wise indicator for the goodness of the fit of T1 maps. R2 values were significantly lower in patients with LBBB than in healthy controls (whole myocardium/septum, 0.997, IQR, 0.00 vs. 0.998, IQR, 0.00; p = 0.008/0.998, IQR, 0.00 vs. 0.999, IQR, 0.00; p = 0.027). Manual correction of semi-automated evaluation tended to improve R2 values but not significantly. Strain analysis was performed and the systolic dyssynchrony index (SDIglobal) was calculated as a measure for left ventricular dyssynchrony. While MRI is generally prone to artefacts, lower goodness of the fit in LBBB may be mainly attributable to asynchronous contraction. Therefore, careful checking of the source data and, if necessary, manual post-processing is important. New techniques might improve the goodness of the fit of T1 mapping by reducing sampling in the motion prone diastole of LBBB patients.

Present and Future Innovations in AI and Cardiac MRI

Cardiac MRI is used to diagnose and treat patients with a multitude of cardiovascular diseases. Despite the growth of clinical cardiac MRI, complicated image prescriptions and long acquisition protocols limit the specialty and restrain its impact on the practice of medicine. Artificial intelligence (AI)-the ability to mimic human intelligence in learning and performing tasks-will impact nearly all aspects of MRI. Deep learning (DL) primarily uses an artificial neural network to learn a specific task from example data sets. Self-driving scanners are increasingly available, where AI automatically controls cardiac image prescriptions. These scanners offer faster image collection with higher spatial and temporal resolution, eliminating the need for cardiac triggering or breath holding. In the future, fully automated inline image analysis will most likely provide all contour drawings and initial measurements to the reader. Advanced analysis using radiomic or DL features may provide new insights and information not typically extracted in the current analysis workflow. AI may further help integrate these features with clinical, genetic, wearable-device, and "omics" data to improve patient outcomes. This article presents an overview of AI and its application in cardiac MRI, including in image acquisition, reconstruction, and processing, and opportunities for more personalized cardiovascular care through extraction of novel imaging markers.

Scanner-Independent MyoMapNet for Accelerated Cardiac MRI T1 Mapping Across Vendors and Field Strengths

BACKGROUND

In cardiac T1 mapping, a series of T1 -weighted (T1 w) images are collected and numerically fitted to a two or three-parameter model of the signal recovery to estimate voxel-wise T1 values. To reduce the scan time, one can collect fewer T1 w images, albeit at the cost of precision or/and accuracy. Recently, the feasibility of using a neural network instead of conventional two- or three-parameter fit modeling has been demonstrated. However, prior studies used data from a single vendor and field strength; therefore, the generalizability of the models has not been established.

PURPOSE

To develop and evaluate an accelerated cardiac T1 mapping approach based on MyoMapNet, a convolution neural network T1 estimator that can be used across different vendors and field strengths by incorporating the relevant scanner information as additional inputs to the model.

STUDY TYPE

Retrospective, multicenter.

POPULATION

A total of 1423 patients with known or suspected cardiac disease (808 male, 57 ± 16 years), from three centers, two vendors (Siemens, Philips), and two field strengths (1.5 T, 3 T). The data were randomly split into 60% training, 20% validation, and 20% testing.

FIELD STRENGTH/SEQUENCE

A 1.5 T and 3 T, Modified Look-Locker inversion recovery (MOLLI) for native and postcontrast T1 .

ASSESSMENT

Scanner-independent MyoMapNet (SI-MyoMapNet) was developed by altering the deep learning (DL) architecture of MyoMapNet to incorporate scanner vendor and field strength as inputs. Epicardial and endocardial contours and blood pool (by manually drawing a large region of interest in the blood pool) of the left ventricle were manually delineated by three readers, with 2, 8, and 9 years of experience, and SI-MyoMapNet myocardial and blood pool T1 values (calculated from four T1 w images) were compared with conventional MOLLI T1 values (calculated from 8 to 11 T1 w images).

STATISTICAL TESTS

Equivalency test with 95% confidence interval (CI), linear regression slope, Pearson correlation coefficient (r), Bland-Altman analysis.

RESULTS

The proposed SI-MyoMapNet successfully created T1 maps. Native and postcontrast T1 values measured from SI-MyoMapNet were strongly correlated with MOLLI, despite using only four T1 w images, at both field-strengths and vendors (all r > 0.86). For native T1 , SI-MyoMapNet and MOLLI were in good agreement for myocardial and blood T1 values in institution 1 (myocardium: 5 msec, 95% CI [3, 8]; blood: -10 msec, 95%CI [-16, -4]), in institution 2 (myocardium: 6 msec, 95% CI [0, 11]; blood: 0 msec, [-18, 17]), and in institution 3 (myocardium: 7 msec, 95% CI [-8, 22]; blood: 8 msec, [-14, 30]). Similar results were observed for postcontrast T1 .

DATA CONCLUSION

Inclusion of field strength and vendor as additional inputs to the DL architecture allows generalizability of MyoMapNet across different vendors or field strength.

EVIDENCE LEVEL

2.

TECHNICAL EFFICACY

Stage 2.

DeepFittingNet: A deep neural network-based approach for simplifying cardiac T1 and T2 estimation with improved robustness

PURPOSE

To develop and evaluate a deep neural network (DeepFittingNet) for T1 /T2 estimation of the most commonly used cardiovascular MR mapping sequences to simplify data processing and improve robustness.

THEORY AND METHODS

DeepFittingNet is a 1D neural network composed of a recurrent neural network (RNN) and a fully connected (FCNN) neural network, in which RNN adapts to the different number of input signals from various sequences and FCNN subsequently predicts A, B, and Tx of a three-parameter model. DeepFittingNet was trained using Bloch-equation simulations of MOLLI and saturation-recovery single-shot acquisition (SASHA) T1 mapping sequences, and T2 -prepared balanced SSFP (T2 -prep bSSFP) T2 mapping sequence, with reference values from the curve-fitting method. Several imaging confounders were simulated to improve robustness. The trained DeepFittingNet was tested using phantom and in-vivo signals, and compared to the curve-fitting algorithm.

RESULTS

In testing, DeepFittingNet performed T1 /T2 estimation of four sequences with improved robustness in inversion-recovery T1 estimation. The mean bias in phantom T1 and T2 between the curve-fitting and DeepFittingNet was smaller than 30 and 1 ms, respectively. Excellent agreements between both methods was found in the left ventricle and septum T1 /T2 with a mean bias <6 ms. There was no significant difference in the SD of both the left ventricle and septum T1 /T2 between the two methods.

CONCLUSION

DeepFittingNet trained with simulations of MOLLI, SASHA, and T2 -prep bSSFP performed T1 /T2 estimation tasks for all these most used sequences. Compared with the curve-fitting algorithm, DeepFittingNet improved the robustness for inversion-recovery T1 estimation and had comparable performance in terms of accuracy and precision.