A Markov random field approach for topology-preserving registration: application to object-based tomographic image interpolation

Time-efficient three-dimensional transmural scar assessment provides relevant substrate characterization for ventricular tachycardia features and long-term recurrences in ischemic cardiomyopathy

Delayed gadolinium-enhanced cardiac magnetic resonance (LGE-CMR) imaging requires novel and time-efficient approaches to characterize the myocardial substrate associated with ventricular arrhythmia in patients with ischemic cardiomyopathy. Using a translational approach in pigs and patients with established myocardial infarction, we tested and validated a novel 3D methodology to assess ventricular scar using custom transmural criteria and a semiautomatic approach to obtain transmural scar maps in ventricular models reconstructed from both 3D-acquired and 3D-upsampled-2D-acquired LGE-CMR images. The results showed that 3D-upsampled models from 2D LGE-CMR images provided a time-efficient alternative to 3D-acquired sequences to assess the myocardial substrate associated with ischemic cardiomyopathy. Scar assessment from 2D-LGE-CMR sequences using 3D-upsampled models was superior to conventional 2D assessment to identify scar sizes associated with the cycle length of spontaneous ventricular tachycardia episodes and long-term ventricular tachycardia recurrences after catheter ablation. This novel methodology may represent an efficient approach in clinical practice after manual or automatic segmentation of myocardial borders in a small number of conventional 2D LGE-CMR slices and automatic scar detection.

Lung 4D CT Image Registration Based on High-Order Markov Random Field

To solve the problem that traditional image registration methods based on continuous optimization for large motion lung 4D CT image sequences are easy to fall into local optimal solutions and lead to serious misregistration, a novel image registration method based on high-order Markov Random Field (MRF) is proposed. By analyzing the effect of the deformation field constraint of the potential functions with different order cliques in MRF model, energy functions with high-order cliques form are designed separately for 2D and 3D images to preserve topology of the deformation field. In order to preserve the topology of the deformation field more effectively, it is necessary to apply a smooth term and a topology preservation term simultaneously in the energy function and use logarithmic function to impose a penalty on the Jacobian matrix with high-order cliques in the topology preservation term. For the complexity of the designed energy function with high-order cliques form, Markov Chain Monte Carlo (MCMC) algorithm is used to solve the optimization problem of the designed energy function. To address the high computational requirements in lung 4D CT image registration, a multi-level processing strategy is adopted to reduce the space complexity of the proposed registration method and promotes the computational efficiency. In the DIR-lab dataset with 4D CT images and the COPD (Chronic Obstructive Pulmonary Disease) dataset with 3D CT images, the average target registration error (TRE) of our proposed method can reach 0.95 mm respectively.

Registration and Modeling From Spaced and Misaligned Image Volumes

We address the problem of object modeling from 3D and 3D+T data made up of images, which contain different parts of an object of interest, are separated by large spaces, and are misaligned with respect to each other. These images have only a limited number of intersections, hence making their registration particularly challenging. Furthermore, such data may result from various medical imaging modalities and can, therefore, present very diverse spatial configurations. Previous methods perform registration and object modeling (segmentation and interpolation) sequentially. However, sequential registration is ill-suited for the case of images with few intersections. We propose a new methodology, which, regardless of the spatial configuration of the data, performs the three stages of registration, segmentation, and shape interpolation from spaced and misaligned images simultaneously. We integrate these three processes in a level set framework, in order to benefit from their synergistic interactions. We also propose a new registration method that exploits segmentation information rather than pixel intensities, and that accounts for the global shape of the object of interest, for increased robustness and accuracy. The accuracy of registration is compared against traditional mutual information based methods, and the total modeling framework is assessed against traditional sequential processing and validated on artificial, CT, and MRI data.

Topology preserving non-rigid image registration using time-varying elasticity model for MRI brain volumes

In this paper, we present a new non-rigid image registration method that imposes a topology preservation constraint on the deformation. We propose to incorporate the time varying elasticity model into the deformable image matching procedure and constrain the Jacobian determinant of the transformation over the entire image domain. The motion of elastic bodies is governed by a hyperbolic partial differential equation, generally termed as elastodynamics wave equation, which we propose to use as a deformation model. We carried out clinical image registration experiments on 3D magnetic resonance brain scans from IBSR database. The results of the proposed registration approach in terms of Kappa index and relative overlap computed over the subcortical structures were compared against the existing topology preserving non-rigid image registration methods and non topology preserving variant of our proposed registration scheme. The Jacobian determinant maps obtained with our proposed registration method were qualitatively and quantitatively analyzed. The results demonstrated that the proposed scheme provides good registration accuracy with smooth transformations, thereby guaranteeing the preservation of topology.

Integrated segmentation and interpolation of sparse data

We address the two inherently related problems of segmentation and interpolation of 3D and 4D sparse data and propose a new method to integrate these stages in a level set framework. The interpolation process uses segmentation information rather than pixel intensities for increased robustness and accuracy. The method supports any spatial configurations of sets of 2D slices having arbitrary positions and orientations. We achieve this by introducing a new level set scheme based on the interpolation of the level set function by radial basis functions. The proposed method is validated quantitatively and/or subjectively on artificial data and MRI and CT scans and is compared against the traditional sequential approach, which interpolates the images first, using a state-of-the-art image interpolation method, and then segments the interpolated volume in 3D or 4D. In our experiments, the proposed framework yielded similar segmentation results to the sequential approach but provided a more robust and accurate interpolation. In particular, the interpolation was more satisfactory in cases of large gaps, due to the method taking into account the global shape of the object, and it recovered better topologies at the extremities of the shapes where the objects disappear from the image slices. As a result, the complete integrated framework provided more satisfactory shape reconstructions than the sequential approach.

Multi-Stencil Streamline Fast Marching: A General 3-D Framework to Determine Myocardial Thickness and Transmurality in Late Enhancement Images

We propose a fully 3-D methodology for the computation of myocardial nonviable tissue transmurality in contrast enhanced magnetic resonance images. The outcome is a continuous map defined within the myocardium where not only current state-of-the-art measures of transmurality can be calculated, but also information on the location of nonviable tissue is preserved. The computation is done by means of a partial differential equation framework we have called multi-stencil streamline fast marching. Using it, the myocardial and scarred tissue thickness is simultaneously computed. Experimental results show that the proposed 3-D method allows for the computation of transmurality in myocardial regions where current 2-D methods are not able to as conceived, and it also provides more robust and accurate results in situations where the assumptions on which current 2-D methods are based-i.e., there is a visible endocardial contour and its corresponding epicardial points lie on the same slice-, are not met.

Estimating the granularity coefficient of a Potts-Markov random field within a Markov chain Monte Carlo algorithm

This paper addresses the problem of estimating the Potts parameter β jointly with the unknown parameters of a Bayesian model within a Markov chain Monte Carlo (MCMC) algorithm. Standard MCMC methods cannot be applied to this problem because performing inference on β requires computing the intractable normalizing constant of the Potts model. In the proposed MCMC method, the estimation of β is conducted using a likelihood-free Metropolis-Hastings algorithm. Experimental results obtained for synthetic data show that estimating β jointly with the other unknown parameters leads to estimation results that are as good as those obtained with the actual value of β. On the other hand, choosing an incorrect value of β can degrade estimation performance significantly. To illustrate the interest of this method, the proposed algorithm is successfully applied to real bidimensional SAR and tridimensional ultrasound images.

A magnetic resonance software simulator for the evaluation of myocardial deformation estimation

This paper proposes a methodology to design a physiologically realistic computer simulator of images of the left ventricle myocardium based on a patient-specific biomechanical model. The simulator takes a magnetic resonance image of a given patient at end diastole, uses a manual segmentation of that image to model the geometry of the myocardium and sets the parameters of the constitutive model used for biomechanical simulation according to a regional labeling of the contractility of the myocardium for that patient. The simulated deformations are used to warp the magnetic resonance dataset throughout the cardiac cycle to generate different image modalities. The simulator is validated by quantifying its ability to model actual deformations in a set of patients affected by an acute myocardial infarction. Specifically a high correlation has been encountered between the ejection fraction derived from the simulated end systolic deformation of the myocardium and the myocardium segmented from actual data. Additionally, most of the parameters that describe the simulated deformation compare well with reported values. Overall, the simulator is intended as a testbed for extensive comparisons of myocardial motion tracking methods due to its ability to relate the impaired myocardial function with the associated ventricular remodeling, a novel contribution in the literature of cardiac image simulators.