Analysis of temporal encoding efficiency of MR fingerprinting sequences.

PURPOSE: MR Fingerprinting (MRF) relies on highly-undersampled images to simultaneously estimate multiple tissue parameters of interest. While a good understanding of the encoding principle behind MRF exists, we want to shed light on the question of when parameters are encoded during an MRF acquisition. THEORY AND METHODS: We analyze the importance of each time point by leaving it out during matching (leave-one-out, LOO) assuming linear reconstruction is applied and study its influence on the reconstructed parameter map. To accelerate the analysis, we treat LOO as a small perturbation (LOOP) to the full matching problem and derive an analytical formula by leveraging Stolk and Sbrizzi's analysis on the interplay of k-space sampling and transient state dynamics. To study the influence of geometry and parameter distribution, we deploy LOOP on randomly sliced 3D brain geometries with randomized T1/T2 values to identify primary encoding regions independent of geometry. RESULTS: LOOP can be evaluated orders of magnitude faster than conventional matching and visualizes temporal encoding efficiency. We use the findings to accelerate an MRF sequence by truncation as well as to remove undersampling artifacts through an iterative omission scheme in an ill-working MRF sequence in both in-silico and in-vivo experiments. CONCLUSION: LOOP is an analytical approach to quantify and visualize parameter encoding in MRF.

Personalization of biomechanical simulations of the left ventricle by in-vivo cardiac DTI data: Impact of fiber interpolation methods.

Simulations of cardiac electrophysiology and mechanics have been reported to be sensitive to the microstructural anisotropy of the myocardium. Consequently, a personalized representation of cardiac microstructure is a crucial component of accurate, personalized cardiac biomechanical models. In-vivo cardiac Diffusion Tensor Imaging (cDTI) is a non-invasive magnetic resonance imaging technique capable of probing the heart's microstructure. Being a rather novel technique, issues such as low resolution, signal-to noise ratio, and spatial coverage are currently limiting factors. We outline four interpolation techniques with varying degrees of data fidelity, different amounts of smoothing strength, and varying representation error to bridge the gap between the sparse in-vivo data and the model, requiring a 3D representation of microstructure across the myocardium. We provide a workflow to incorporate in-vivo myofiber orientation into a left ventricular model and demonstrate that personalized modelling based on fiber orientations from in-vivo cDTI data is feasible. The interpolation error is correlated with a trend in personalized parameters and simulated physiological parameters, strains, and ventricular twist. This trend in simulation results is consistent across material parameter settings and therefore corresponds to a bias introduced by the interpolation method. This study suggests that using a tensor interpolation approach to personalize microstructure with in-vivo cDTI data, reduces the fiber uncertainty and thereby the bias in the simulation results.

Comparison of interpolation methods of predominant cardiomyocyte orientation from in vivo and ex vivo cardiac diffusion tensor imaging data.

Cardiac electrophysiology and cardiac mechanics both depend on the average cardiomyocyte long-axis orientation. In the realm of personalized medicine, knowledge of the patient-specific changes in cardiac microstructure plays a crucial role. Patient-specific computational modelling has emerged as a tool to better understand disease progression. In vivo cardiac diffusion tensor imaging (cDTI) is a vital tool to non-destructively measure the average cardiomyocyte long-axis orientation in the heart. However, cDTI suffers from long scan times, rendering volumetric, high-resolution acquisitions challenging. Consequently, interpolation techniques are needed to populate bio-mechanical models with patient-specific average cardiomyocyte long-axis orientations. In this work, we compare five interpolation techniques applied to in vivo and ex vivo porcine input data. We compare two tensor interpolation approaches, one rule-based approximation, and two data-driven, low-rank models. We demonstrate the advantage of tensor interpolation techniques, resulting in lower interpolation errors than do low-rank models and rule-based methods adapted to cDTI data. In an ex vivo comparison, we study the influence of three imaging parameters that can be traded off against acquisition time: in-plane resolution, signal to noise ratio, and number of acquired short-axis imaging slices.

A 3D personalized cardiac myocyte aggregate orientation model using MRI data-driven low-rank basis functions.

Cardiac myocyte aggregate orientation has a strong impact on cardiac electrophysiology and mechanics. Studying the link between structural characteristics, strain, and stresses over the cardiac cycle and cardiac function requires a full volumetric representation of the microstructure. In this work, we exploit the structural similarity across hearts to extract a low-rank representation of predominant myocyte orientation in the left ventricle from high-resolution magnetic resonance ex-vivo cardiac diffusion tensor imaging (cDTI) in porcine hearts. We compared two reduction methods, Proper Generalized Decomposition combined with Singular Value Decomposition and Proper Orthogonal Decomposition. We demonstrate the existence of a general set of basis functions of aggregated myocyte orientation which defines a data-driven, personalizable, parametric model featuring higher flexibility than existing atlas and rule-based approaches. A more detailed representation of microstructure matching the available patient data can improve the accuracy of personalized computational models. Additionally, we approximate the myocyte orientation of one ex-vivo human heart and demonstrate the feasibility of transferring the basis functions to humans.