Over the years, medical image tracking has gained considerable attention from both medical and research communities due to its widespread utility in a multitude of clinical applications, from functional assessment during diagnosis and therapy planning to structure tracking or image fusion during image-guided interventions. Despite the ever-increasing number of image tracking methods available, most still consist of independent implementations with specific target applications, lacking the versatility to deal with distinct end-goals without the need for methodological tailoring and/or exhaustive tuning of numerous parameters. With this in mind, we have developed the medical image tracking toolbox (MITT)-a software package designed to ease customization of image tracking solutions in the medical field. While its workflow principles make it suitable to work with 2-D or 3-D image sequences, its modules offer versatility to set up computationally efficient tracking solutions, even for users with limited programming skills. MITT is implemented in both C/C++ and MATLAB, including several variants of an object-based image tracking algorithm and allowing to track multiple types of objects (i.e., contours, multi-contours, surfaces, and multi-surfaces) with several customization features. In this paper, the toolbox is presented, its features discussed, and illustrative examples of its usage in the cardiology field provided, demonstrating its versatility, simplicity, and time efficiency.

A plethora of techniques for cardiac deformation imaging with 3D ultrasound, typically referred to as 3D speckle tracking techniques, are available from academia and industry. Although the benefits of single methods over alternative ones have been reported in separate publications, the intrinsic differences in the data and definitions used makes it hard to compare the relative performance of different solutions. To address this issue, we have recently proposed a framework to simulate realistic 3D echocardiographic recordings and used it to generate a common set of ground-truth data for 3D speckle tracking algorithms, which was made available online. The aim of this study was therefore to use the newly developed database to contrast non-commercial speckle tracking solutions from research groups with leading expertise in the field. The five techniques involved cover the most representative families of existing approaches, namely block-matching, radio-frequency tracking, optical flow and elastic image registration. The techniques were contrasted in terms of tracking and strain accuracy. The feasibility of the obtained strain measurements to diagnose pathology was also tested for ischemia and dyssynchrony.

Real-time 3D Echocardiography (RT3DE) has been proven to be an accurate tool for left ventricular (LV) volume assessment. However, identification of the LV endocardium remains a challenging task, mainly because of the low tissue/blood contrast of the images combined with typical artifacts. Several semi and fully automatic algorithms have been proposed for segmenting the endocardium in RT3DE data in order to extract relevant clinical indices, but a systematic and fair comparison between such methods has so far been impossible due to the lack of a publicly available common database. Here, we introduce a standardized evaluation framework to reliably evaluate and compare the performance of the algorithms developed to segment the LV border in RT3DE. A database consisting of 45 multivendor cardiac ultrasound recordings acquired at different centers with corresponding reference measurements from three experts are made available. The algorithms from nine research groups were quantitatively evaluated and compared using the proposed online platform. The results showed that the best methods produce promising results with respect to the experts' measurements for the extraction of clinical indices, and that they offer good segmentation precision in terms of mean distance error in the context of the experts' variability range. The platform remains open for new submissions.

Myocardial deformation imaging can provide valuable insights in myocardial mechanics and help in the diagnosis, prognosis and follow-up of cardiac diseases. However, extracting these indices in 3D is challenging due to the limitations in spatial and temporal resolution of the current volumetric ultrasound systems. For this purpose, we developed an anatomical free-form deformation image registration framework which is locally adapted to the anatomy of the heart. In this work we explored whether incorporating a myocardial volume conservation regularizer would improve strain estimates. We evaluated our technique on in silico echo sequences featuring realistic speckle textures and showed the volume conservation regularizer to be beneficial in reducing strain errors further when used in combination with a smoothness penalty. This combination led to more physiological boundary conditions. It also made distinguishing ischemic from normal segments easier in clinical images.

A novel semi-automatic algorithm for aortic valve (AV) wall segmentation is presented for 3D transesophageal echocardiography (TEE) datasets. The proposed methodology uses a 3D cylindrical formulation of the B-spline Explicit Active Surfaces (BEAS) framework in a dual-stage energy evolution process, comprising a threshold-based and a localized region-based stage. Hereto, intensity and shape-based features are combined to accurately delineate the AV wall from the ascending aorta (AA) to the left ventricular outflow tract (LVOT). Shape-prior information is included using a profile-based statistical shape model (SSM), and embedded in BEAS through two novel regularization terms: one confining the segmented AV profiles to shapes seen in the SSM (hard regularization) and another penalizing according to the profile's degree of likelihood (soft regularization). The proposed energy functional takes thus advantage of the intensity data in regions with strong image content, while complementing it with shape knowledge in regions with nearly absent image data. The proposed algorithm has been validated in 20 3D-TEE datasets with both stenotic and non-stenotic valves. It was shown to be accurate, robust and computationally efficient, taking less than 1 second to segment the AV wall from the AA to the LVOT with an average accuracy of 0.78 mm. Semi-automatically extracted measurements at four relevant anatomical levels (LVOT, aortic annulus, sinuses of Valsalva and sinotubular junction) showed an excellent agreement with experts' ones, with a higher reproducibility than manually-extracted measures.

Tracking of moving objects in an image sequence is an important task in many application (e.g. medical imaging, robotics). However, this task is usually difficult due to inherent problems that could happen in sequences (i.e. possible occlusion of the object, large interframe motion). In this paper, we describe a new approach to integrate a priori motion information into a level-set-based active contour approach. Specifically, we introduce a new constraint that enforces the conservation of the levels of the implicit function along the image sequence. This constraint is formulated as a motion prior energy and is used in a tracking algorithm. The method is validated quantitatively on a clinical application based on 10 echocar-diographic and 5 cine-MRI sequences (≈ 700 images).

Despite the availability of multiple solutions for assessing myocardial strain by ultrasound, little is currently known about the relative performance of the different methods. In this study, we sought to contrast two strain estimation techniques directly (speckle tracking and elastic registration) in an in vivo setting by comparing both to a gold standard reference measurement. In five open-chest sheep instrumented with ultrasonic microcrystals, 2-D images were acquired with a GE Vivid7 ultrasound system. Radial (ε(RR)), longitudinal (ε(LL)), and circumferential strain (ε(CC)) were estimated during four inotropic stages: at rest, during esmolol and dobutamine infusion, and during acute ischemia. The correlation of the end-systolic strain values of a well-validated speckle tracking approach and an elastic registration method against sonomicrometry were comparable for ε(LL) ( r=0.70 versus r=0.61, respectively; p=0.32) and ε(CC) ( r=0.73 versus r=0.80 respectively; p=0.31). However, the elastic registration method performed considerably better for ε(RR) ( r=0.64 versus r=0.85 respectively; p=0.09). Moreover, the bias and limits of agreement with respect to the reference strain estimates were statistically significantly smaller in this direction . This could be related to regularization which is imposed during the motion estimation process as opposed to an a posteriori regularization step in the speckle tracking method. Whether one method outperforms the other in detecting dysfunctional regions remains the topic of future research.

A new formulation of active contours based on explicit functions has been recently suggested. This novel framework allows real-time 3-D segmentation since it reduces the dimensionality of the segmentation problem. In this paper, we propose a B-spline formulation of this approach, which further improves the computational efficiency of the algorithm. We also show that this framework allows evolving the active contour using local region-based terms, thereby overcoming the limitations of the original method while preserving computational speed. The feasibility of real-time 3-D segmentation is demonstrated using simulated and medical data such as liver computer tomography and cardiac ultrasound images.

We introduce in this paper a formulation of region-based active contours using explicit surfaces. By formally relating the explicit formulation to the usual implicit formulation, we show that it allows to take advantage of the local or global region-based functionals designed in the level-set framework, while providing a fast algorithm yielding close to real-time segmentation of 3D objects. We compare the segmentation performance of the proposed method with the fast level-set based method of Shi [8]. In similar conditions, the proposed method is 26x to 46x faster. Furthermore, we show that it can handle challenging, inhomogeneous data, while keeping a low computational burden.