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

Computing thickness of irregularly-shaped thin walls using a locally semi-implicit scheme with extrapolation to solve the Laplace equation: Application to the right ventricle

Cardiac Magnetic Resonance (CMR) Imaging is currently considered the gold standard imaging modality in cardiology. However, it is accompanied by a tradeoff between spatial resolution and acquisition time. Providing accurate measures of thin walls relative to the image resolution may prove challenging. One such anatomical structure is the cardiac right ventricle. Methods for measuring thickness of wall-like anatomical structures often rely on the Laplace equation to provide point-to-point correspondences between both boundaries. This work presents limex, a novel method to solve the Laplace equation using ghost nodes and providing extrapolated values, which is tested on three different datasets: a mathematical phantom, a set of biventricular segmentations from CMR images of ten pigs and the database used at the RV Segmentation Challenge held at MICCAI'12. Thickness measurements using the proposed methodology are more accurate than state-of-the-art methods, especially with the coarsest image resolutions, yielding mean L1 norms of the error between 43.28% and 86.52% lower than the second-best methods on the different test datasets. It is also computationally affordable. Limex has outperformed other state-of-the-art methods in classifying RV myocardial segments by their thickness.

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.

Vessel network extraction and analysis of mouse pulmonary vasculature via X-ray micro-computed tomographic imaging

In this work, non-invasive high-spatial resolution three-dimensional (3D) X-ray micro-computed tomography (μCT) of healthy mouse lung vasculature is performed. Methodologies are presented for filtering, segmenting, and skeletonizing the collected 3D images. Novel methods for the removal of spurious branch artefacts from the skeletonized 3D image are introduced, and these novel methods involve a combination of distance transform gradients, diameter-length ratios, and the fast marching method (FMM). These new techniques of spurious branch removal result in the consistent removal of spurious branches without compromising the connectivity of the pulmonary circuit. Analysis of the filtered, skeletonized, and segmented 3D images is performed using a newly developed Vessel Network Extraction algorithm to fully characterize the morphology of the mouse pulmonary circuit. The removal of spurious branches from the skeletonized image results in an accurate representation of the pulmonary circuit with significantly less variability in vessel diameter and vessel length in each generation. The branching morphology of a full pulmonary circuit is characterized by the mean diameter per generation and number of vessels per generation. The methods presented in this paper lead to a significant improvement in the characterization of 3D vasculature imaging, allow for automatic separation of arteries and veins, and for the characterization of generations containing capillaries and intrapulmonary arteriovenous anastomoses (IPAVA).

A Second Order Multi-Stencil Fast Marching Method with a Non-Constant Local Cost Model

The Fast Marching method is widely employed in several fields of image processing. Some years ago a Multi-Stencil version (MSFM) was introduced to improve its accuracy by solving the equation for a set of stencils and choosing the best solution at each considered node. The following work proposes a modified numerical scheme for MSFM to take into account the variation of the local cost, which has proven to be second order. The influence of the stencil set choice on the algorithm outcome with respect to stencil orthogonality and axis swapping is also explored, where stencils are taken from neighborhoods of varying radius. The experimental results show that the proposed schemes improve the accuracy of their original counterparts, and that the use of permutation-invariant stencil sets provides robustness against shifted vector coordinates in the stencil set.