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.

Patient-specific electromechanical models of the heart for the prediction of pacing acute effects in CRT: a preliminary clinical validation.

Cardiac resynchronisation therapy (CRT) is an effective treatment for patients with congestive heart failure and a wide QRS complex. However, up to 30% of patients are non-responders to therapy in terms of exercise capacity or left ventricular reverse remodelling. A number of controversies still remain surrounding patient selection, targeted lead implantation and optimisation of this important treatment. The development of biophysical models to predict the response to CRT represents a potential strategy to address these issues. In this article, we present how the personalisation of an electromechanical model of the myocardium can predict the acute haemodynamic changes associated with CRT. In order to introduce such an approach as a clinical application, we needed to design models that can be individualised from images and electrophysiological mapping of the left ventricle. In this paper the personalisation of the anatomy, the electrophysiology, the kinematics and the mechanics are described. The acute effects of pacing on pressure development were predicted with the in silico model for several pacing conditions on two patients, achieving good agreement with invasive haemodynamic measurements: the mean error on dP/dt(max) is 47.5±35mmHgs(-1), less than 5% error. These promising results demonstrate the potential of physiological models personalised from images and electrophysiology signals to improve patient selection and plan CRT.

LBM-EP: Lattice-Boltzmann method for fast cardiac electrophysiology simulation from 3D images.

Current treatments of heart rhythm troubles require careful planning and guidance for optimal outcomes. Computational models of cardiac electrophysiology are being proposed for therapy planning but current approaches are either too simplified or too computationally intensive for patient-specific simulations in clinical practice. This paper presents a novel approach, LBM-EP, to solve any type of mono-domain cardiac electrophysiology models at near real-time that is especially tailored for patient-specific simulations. The domain is discretized on a Cartesian grid with a level-set representation of patient's heart geometry, previously estimated from images automatically. The cell model is calculated node-wise, while the transmembrane potential is diffused using Lattice-Boltzmann method within the domain defined by the level-set. Experiments on synthetic cases, on a data set from CESC'10 and on one patient with myocardium scar showed that LBM-EP provides results comparable to an FEM implementation, while being 10 - 45 times faster. Fast, accurate, scalable and requiring no specific meshing, LBM-EP paves the way to efficient and detailed models of cardiac electrophysiology for therapy planning.

A personalized biomechanical model for respiratory motion prediction.

Time-resolved imaging of the thorax or abdominal area is affected by respiratory motion. Nowadays, one-dimensional respiratory surrogates are used to estimate the current state of the lung during its cycle, but with rather poor results. This paper presents a framework to predict the 3D lung motion based on a patient-specific finite element model of respiratory mechanics estimated from two CT images at end of inspiration (EI) and end of expiration (EE). We first segment the lung, thorax and sub-diaphragm organs automatically using a machine-learning algorithm. Then, a biomechanical model of the lung, thorax and sub-diaphragm is employed to compute the 3D respiratory motion. Our model is driven by thoracic pressures, estimated automatically from the EE and EI images using a trust-region approach. Finally, lung motion is predicted by modulating the thoracic pressures. The effectiveness of our approach is evaluated by predicting lung deformation during exhale on five DIR-Lab datasets. Several personalization strategies are tested, showing that an average error of 3.88 +/- 1.54 mm in predicted landmark positions can be achieved. Since our approach is generative, it may constitute a 3D surrogate information for more accurate medical image reconstruction and patient respiratory analysis.

Towards patient-specific finite-element simulation of MitralClip procedure.

MitralClip is a novel minimally invasive procedure to treat mitral valve (MV) regurgitation. It consists in clipping the mitral leaflets together to close the regurgitant hole. A careful preoperative planning is necessary to select respondent patients and to determine the clipping sites. Although preliminary indications criteria are established, they lack prediction power with respect to complications and effectiveness of the therapy in specific patients. We propose an integrated framework for personalized simulation of MV function and apply it to simulate MitralClip procedure. A patient-specific dynamic model of the MV apparatus is computed automatically from 4D TEE images. A biomechanical model of the MV, constrained by the observed motion of the mitral annulus and papillary muscles, is employed to simulate valve closure and MitralClip intervention. The proposed integrated framework enables, for the first time, to quantitatively evaluate an MV finite-element model in-vivo, on eleven patients, and to predict the outcome of MitralClip intervention in one of these patients. The simulations are compared to ground truth and to postoperative images, resulting in promising accuracy (average point-to-mesh distance: 1.47 +/- 0.24 mm). Our framework may constitute a tool for MV therapy planning and patient management.

A statistical model for quantification and prediction of cardiac remodelling: application to tetralogy of Fallot.

Cardiac remodelling plays a crucial role in heart diseases. Analyzing how the heart grows and remodels over time can provide precious insights into pathological mechanisms, eventually resulting in quantitative metrics for disease evaluation and therapy planning. This study aims to quantify the regional impacts of valve regurgitation and heart growth upon the end-diastolic right ventricle (RV) in patients with tetralogy of Fallot, a severe congenital heart defect. The ultimate goal is to determine, among clinical variables, predictors for the RV shape from which a statistical model that predicts RV remodelling is built. Our approach relies on a forward model based on currents and a diffeomorphic surface registration algorithm to estimate an unbiased template. Local effects of RV regurgitation upon the RV shape were assessed with Principal Component Analysis (PCA) and cross-sectional multivariate design. A generative 3-D model of RV growth was then estimated using partial least squares (PLS) and canonical correlation analysis (CCA). Applied on a retrospective population of 49 patients, cross-effects between growth and pathology could be identified. Qualitatively, the statistical findings were found realistic by cardiologists. 10-fold cross-validation demonstrated a promising generalization and stability of the growth model. Compared to PCA regression, PLS was more compact, more precise and provided better predictions.

LogDemons revisited: consistent regularisation and incompressibility constraint for soft tissue tracking in medical images.

Non-linear image registration is a standard approach to track soft tissues in medical images. By estimating spatial transformations between images, visible structures can be followed over time. For clinical applications the model of transformation must be consistent with the properties of the biological tissue, such as incompressibility. LogDemons is a fast non-linear registration algorithm that provides diffusion-like diffeomorphic transformations parameterised by stationary velocity fields. Yet, its use for tissue tracking has been limited because of the ad-hoc Gaussian regularisation that prevents implementing other transformation models. In this paper, we propose a mathematical formulation of demons regularisation that fits into LogDemons framework. This formulation enables to ensure volume-preserving deformations by minimising the energy functional directly under the linear divergence-free constraint, yielding little computational overhead. Tests on synthetic incompressible fields showed that our approach outperforms the original logDemons in terms of incompressible deformation recovery. The algorithm showed promising results on one patient for the automatic recovery of myocardium strain from cardiac anatomical and 3D tagged MRI.

Segmentation of focal cortical dysplasia lesions on MRI using level set evolution.

Focal cortical dysplasia (FCD) is the most frequent malformation of cortical development in patients with medically intractable epilepsy. On MRI, FCD lesions are not easily differentiable from the normal cortex and defining their spatial extent is challenging. In this paper, we introduce a method to segment FCD lesions on T1-weighted MRI. It relies on two successive three-dimensional deformable models, whose evolutions are based on the level set framework. The first deformable model is driven by probability maps obtained from three MRI features: cortical thickness, relative intensity and gradient. These features correspond to the visual characteristics of FCD and allow discriminating lesions and normal tissues. In a second stage, the previous result is expanded towards the underlying and overlying cortical boundaries, throughout the whole cortical section. The method was quantitatively evaluated by comparison with manually traced labels in 18 patients with FCD. The automated segmentations achieved a strong agreement with the manuals labels, demonstrating the applicability of the method to assist the delineation of FCD lesions on MRI. This new approach may become a useful tool for the presurgical evaluation of patients with intractable epilepsy related to cortical dysplasia.

Segmentation of focal cortical dysplasia lesions using a feature-based level set.

Focal cortical dysplasia (FCD), a malformation of cortical development, is an important cause of medically intractable epilepsy. FCD lesions are difficult to distinguish from non-lesional cortex and their delineation on MRI is a challenging task. This paper presents a method to segment FCD lesions on T1-weighted MRI, based on a 3D deformable model, implemented using the level set framework. The deformable model is driven by three MRI features: cortical thickness, relative intensity and gradient. These features correspond to the visual characteristics of FCD and allow to differentiate lesions from normal tissues. The proposed method was tested on 18 patients with FCD and its performance was quantitatively evaluated by comparison with the manual tracings of two trained raters. The validation showed that the similarity between the level set segmentation and the manual labels is similar to the agreement between the two human raters. This new approach may become a useful tool for the presurgical evaluation of patients with intractable epilepsy.