Personalization of a cardiac electromechanical model using reduced order unscented Kalman filtering from regional volumes.

Patient-specific cardiac modeling can help in understanding pathophysiology and therapy planning. However it requires to combine functional and anatomical data in order to build accurate models and to personalize the model geometry, kinematics, electrophysiology and mechanics. Personalizing the electromechanical coupling from medical images is a challenging task. We use the Bestel-Clément-Sorine (BCS) electromechanical model of the heart, which provides reasonable accuracy with a reasonable number of parameters (14 for each ventricle) compared to the available clinical data at the organ level. We propose a personalization strategy from cine MRI data in two steps. We first estimate global parameters with an automatic calibration algorithm based on the Unscented Transform which allows to initialize the parameters while matching the volume and pressure curves. In a second step we locally personalize the contractilities of all AHA (American Heart Association) zones of the left ventricle using the reduced order unscented Kalman filtering on Regional Volumes. This personalization strategy was validated synthetically and tested successfully on eight healthy and three pathological cases.

Benchmarking framework for myocardial tracking and deformation algorithms: an open access database.

In this paper we present a benchmarking framework for the validation of cardiac motion analysis algorithms. The reported methods are the response to an open challenge that was issued to the medical imaging community through a MICCAI workshop. The database included magnetic resonance (MR) and 3D ultrasound (3DUS) datasets from a dynamic phantom and 15 healthy volunteers. Participants processed 3D tagged MR datasets (3DTAG), cine steady state free precession MR datasets (SSFP) and 3DUS datasets, amounting to 1158 image volumes. Ground-truth for motion tracking was based on 12 landmarks (4 walls at 3 ventricular levels). They were manually tracked by two observers in the 3DTAG data over the whole cardiac cycle, using an in-house application with 4D visualization capabilities. The median of the inter-observer variability was computed for the phantom dataset (0.77 mm) and for the volunteer datasets (0.84 mm). The ground-truth was registered to 3DUS coordinates using a point based similarity transform. Four institutions responded to the challenge by providing motion estimates for the data: Fraunhofer MEVIS (MEVIS), Bremen, Germany; Imperial College London - University College London (IUCL), UK; Universitat Pompeu Fabra (UPF), Barcelona, Spain; Inria-Asclepios project (INRIA), France. Details on the implementation and evaluation of the four methodologies are presented in this manuscript. The manually tracked landmarks were used to evaluate tracking accuracy of all methodologies. For 3DTAG, median values were computed over all time frames for the phantom dataset (MEVIS=1.20mm, IUCL=0.73 mm, UPF=1.10mm, INRIA=1.09 mm) and for the volunteer datasets (MEVIS=1.33 mm, IUCL=1.52 mm, UPF=1.09 mm, INRIA=1.32 mm). For 3DUS, median values were computed at end diastole and end systole for the phantom dataset (MEVIS=4.40 mm, UPF=3.48 mm, INRIA=4.78 mm) and for the volunteer datasets (MEVIS=3.51 mm, UPF=3.71 mm, INRIA=4.07 mm). For SSFP, median values were computed at end diastole and end systole for the phantom dataset(UPF=6.18 mm, INRIA=3.93 mm) and for the volunteer datasets (UPF=3.09 mm, INRIA=4.78 mm). Finally, strain curves were generated and qualitatively compared. Good agreement was found between the different modalities and methodologies, except for radial strain that showed a high variability in cases of lower image quality.

Understanding the mechanisms amenable to CRT response: from pre-operative multimodal image data to patient-specific computational models.

This manuscript describes our recent developments towards better understanding of the mechanisms amenable to cardiac resynchronization therapy response. We report the results from a full multimodal dataset corresponding to eight patients from the euHeart project. The datasets include echocardiography, MRI and electrophysiological studies. We investigate two aspects. The first one focuses on pre-operative multimodal image data. From 2D echocardiography and 3D tagged MRI images, we compute atlas based dyssynchrony indices. We complement these indices with presence and extent of scar tissue and correlate them with CRT response. The second one focuses on computational models. We use pre-operative imaging to generate a patient-specific computational model. We show results of a fully automatic personalized electromechanical simulation. By case-per-case discussion of the results, we highlight the potential and key issues of this multimodal pipeline for the understanding of the mechanisms of CRT response and a better patient selection.

Generating anatomical models of the heart and the aorta from medical images for personalized physiological simulations.

The anatomy and motion of the heart and the aorta are essential for patient-specific simulations of cardiac electrophysiology, wall mechanics and hemodynamics. Within the European integrated project euHeart, algorithms have been developed that allow to efficiently generate patient-specific anatomical models from medical images from multiple imaging modalities. These models, for instance, account for myocardial deformation, cardiac wall motion, and patient-specific tissue information like myocardial scar location. Furthermore, integration of algorithms for anatomy extraction and physiological simulations has been brought forward. Physiological simulations are linked closer to anatomical models by encoding tissue properties, like the muscle fibers, into segmentation meshes. Biophysical constraints are also utilized in combination with image analysis to assess tissue properties. Both examples show directions of how physiological simulations could provide new challenges and stimuli for image analysis research in the future.

Realistic simulation of cardiac magnetic resonance studies modeling anatomical variability, trabeculae, and papillary muscles.

Simulated magnetic resonance imaging brain studies have been generated for over a decade. Despite their useful potential, simulated cardiac studies are only emerging. This article focuses on the realistic simulation of cardiac magnetic resonance imaging datasets. The methodology is based on the XCAT phantom, which is modified to increase realism of the simulated images. Modifications include the modeling of trabeculae and papillary muscles based on clinical measurements and published data. To develop and evaluate our approach, the clinical database included 40 patients for anatomical measurements, 10 patients for papillary muscle modeling, and 10 patients for local gray value statistics. The virtual database consisted of 40 digital voxel phantoms. Histograms from different tissues were obtained from the real datasets and compared with histograms of the simulated datasets with the Chi-square dissimilarity metric (χ(2)) and Kullback-Leibler divergence. For the original phantom, χ(2) values averaged 0.65 ± 0.06 and Kullboek-Leibler values averaged 0.69 ± 0.38. For the modified phantom, χ(2) values averaged 0.34 ± 0.12 and Kullboek-Leibler values averaged 0.32 ± 0.15. The proposed approach demonstrated a noticeable improvement of the local appearance of the simulated images with respect to the ones obtained originally.

Simulation of late gadolinium enhancement cardiac magnetic resonance studies.

In this study we propose a pipeline for simulation of late gadolinium enhancement images. We used a modified version of the XCAT phantom to improve simulation realism. Modifications included the modeling of trabeculae and papillary muscles, and the increase of sublabels to resemble tissue intensity variability. Magnetic properties for each body tissue were sampled in three settings: from Gaussian distributions, combining Rayleigh-Gaussian distributions, and from Rayleigh distributions. Thirty-two simulated datasets were compared with 32 clinical datasets from infarcted patients. Histograms were obtained for five tissues: lung, pericardium, myocardium, blood and hyper-enhanced area. Real and simulated histograms were compared with the Chi-square dissimilarity metric (χ(2)) and Kullback-Leibler divergence (KL). The generated simulated images look similar to real images according to both metrics. Rayleigh and the Rayleigh-Gaussian models obtained comparable average results (respectively: χ(2)= 0.16 ± 0.12 and 0.18 ± 0.11; KL=0.15 ± 0.17 and 0.16 ± 0.18).

3D mesh based wall thickness measurement: identification of left ventricular hypertrophy phenotypes.

Left ventricular hypertrophy (LVH) is a complex cardiac condition mainly identified by the thickening of the myocardial wall. Although most of the contemporary cardiac imaging modalities provide high resolution 3D images, the wall thickness (WT) is still measured within the acquired planes. This way of measurement may introduce an error as cardiac wall is not necessarily orthogonal to the plane. In this study we analyze how different approaches to measure WT can affect an automatic identification of hypertrophy. The compared approaches are: WT measured along surface normal and the one provided by a medial surface. For both approaches we evaluated their ability to identify LVH phenotypes by testing with two classifiers: Transductive Confidence Machine-k Nearest Neighbor (TCM-kNN) and Linear Discriminant Analysis (LDA). Fifty three subjects were included in this study: 18 patients with hypertrophic cardiomyopathy (HCM), 13 patients with hypertensive heart disease (HDD) and 22 sedentary subjects (CG). Medial surface based approach allowed obtaining higher classification accuracy in HDD patients, while normal based approach allowed for higher classification accuracy in HCM patients.

Myocardial motion estimation in tagged MR sequences by using alphaMI-based non rigid registration.

Tagged Magnetic Resonance Imaging (MRI) is currently the reference MR modality for myocardial motion and strain analysis. NMI-based non rigid registration has proven to be an accurate method to retrieve cardiac deformation fields. The use of alphaMI permits higher dimensional features to be implemented in myocardial deformation estimation through image registration. This paper demonstrates that this is feasible with a set of Haar wavelet features of high dimension. While we do not demonstrate performance improvement for this set of features, there is no significant degradation as compared to implementing the registration method with the traditional NMI metric. We use Entropic Spanning Graphs (ESGs) to estimate the alphaMI of the wavelet feature vectors WFVs since this is not possible with histograms. To the best of our knowledge, this is the first time that ESGs are used for non rigid registration.