Anisotropic filtering for model-based segmentation of 4D cylindrical echocardiographic images

Atlas-Based Shared-Boundary Deformable Multi-Surface Models through Multi-Material and Two-Manifold Dual Contouring

This paper presents a multi-material dual “contouring” method used to convert a digital 3D voxel-based atlas of basal ganglia to a deformable discrete multi-surface model that supports surgical navigation for an intraoperative MRI-compatible surgical robot, featuring fast intraoperative deformation computation. It is vital that the final surface model maintain shared boundaries where appropriate so that even as the deep-brain model deforms to reflect intraoperative changes encoded in ioMRI, the subthalamic nucleus stays in contact with the substantia nigra, for example, while still providing a significantly sparser representation than the original volumetric atlas consisting of hundreds of millions of voxels. The dual contouring (DC) algorithm is a grid-based process used to generate surface meshes from volumetric data. The DC method enables the insertion of vertices anywhere inside the grid cube, as opposed to the marching cubes (MC) algorithm, which can insert vertices only on the grid edges. This multi-material DC method is then applied to initialize, by duality, a deformable multi-surface simplex model, which can be used for multi-surface atlas-based segmentation. We demonstrate our proposed method on synthetic and deep-brain atlas data, and a comparison of our method’s results with those of traditional DC demonstrates its effectiveness.

Model-Agnostic Temporal Regularizer for Object Localization Using Motion Fields

Video analysis often requires locating and tracking target objects. In some applications, the localization system has access to the full video, which allows fine-grain motion information to be estimated. This paper proposes capturing this information through motion fields and using it to improve the localization results. The learned motion fields act as a model-agnostic temporal regularizer that can be used with any localization system based on keypoints. Unlike optical flow-based strategies, our motion fields are estimated from the model domain, based on the trajectories described by the object keypoints. Therefore, they are not affected by poor imaging conditions. The benefits of the proposed strategy are shown on three applications: 1) segmentation of cardiac magnetic resonance; 2) facial model alignment; and 3) vehicle tracking. In each case, combining popular localization methods with the proposed regularizer leads to improvement in overall accuracies and reduces gross errors.

Distance regularized two level sets for segmentation of left and right ventricles from cine-MRI

This paper presents a new level set method for segmentation of cardiac left and right ventricles. We extend the edge based distance regularized level set evolution (DRLSE) model in Li et al. (2010) to a two-level-set formulation, with the 0-level set and k-level set representing the endocardium and epicardium, respectively. The extraction of endocardium and epicardium is obtained as a result of the interactive curve evolution of the 0 and k level sets derived from the proposed variational level set formulation. The initialization of the level set function in the proposed two-level-set DRLSE model is generated from roughly located endocardium, which can be performed by applying the original DRLSE model. Experimental results have demonstrated the effectiveness of the proposed two-level-set DRLSE model.

New developments in paediatric cardiac functional ultrasound imaging

Ultrasound imaging can be used to estimate the morphology as well as the motion and deformation of tissues. If the interrogated tissue is actively deforming, this deformation is directly related to its function and quantification of this deformation is normally referred as 'strain imaging'. Tissue can also be deformed by applying an internal or external force and the resulting, induced deformation is a function of the mechanical tissue characteristics. In combination with the load applied, these strain maps can be used to estimate or reconstruct the mechanical properties of tissue. This technique was named 'elastography' by Ophir et al. in 1991. Elastography can be used for atherosclerotic plaque characterisation, while the contractility of the heart or skeletal muscles can be assessed with strain imaging. Rather than using the conventional video format (DICOM) image information, radio frequency (RF)-based ultrasound methods enable estimation of the deformation at higher resolution and with higher precision than commercial methods using Doppler (tissue Doppler imaging) or video image data (2D speckle tracking methods). However, the improvement in accuracy is mainly achieved when measuring strain along the ultrasound beam direction, so it has to be considered a 1D technique. Recently, this method has been extended to multiple directions and precision further improved by using spatial compounding of data acquired at multiple beam steered angles. Using similar techniques, the blood velocity and flow can be determined. RF-based techniques are also beneficial for automated segmentation of the ventricular cavities. In this paper, new developments in different techniques of quantifying cardiac function by strain imaging, automated segmentation, and methods of performing blood flow imaging are reviewed and their application in paediatric cardiology is discussed.

Combining multiple dynamic models and deep learning architectures for tracking the left ventricle endocardium in ultrasound data

We present a new statistical pattern recognition approach for the problem of left ventricle endocardium tracking in ultrasound data. The problem is formulated as a sequential importance resampling algorithm such that the expected segmentation of the current time step is estimated based on the appearance, shape, and motion models that take into account all previous and current images and previous segmentation contours produced by the method. The new appearance and shape models decouple the affine and nonrigid segmentations of the left ventricle to reduce the running time complexity. The proposed motion model combines the systole and diastole motion patterns and an observation distribution built by a deep neural network. The functionality of our approach is evaluated using a dataset of diseased cases containing 16 sequences and another dataset of normal cases comprised of four sequences, where both sets present long axis views of the left ventricle. Using a training set comprised of diseased and healthy cases, we show that our approach produces more accurate results than current state-of-the-art endocardium tracking methods in two test sequences from healthy subjects. Using three test sequences containing different types of cardiopathies, we show that our method correlates well with interuser statistics produced by four cardiologists.

A Comparison of Hyperelastic Warping of PET Images with Tagged MRI for the Analysis of Cardiac Deformation

The objectives of the following research were to evaluate the utility of a deformable image registration technique known as hyperelastic warping for the measurement of local strains in the left ventricle through the analysis of clinical, gated PET image datasets. Two normal human male subjects were sequentially imaged with PET and tagged MRI imaging. Strain predictions were made for systolic contraction using warping analyses of the PET images and HARP based strain analyses of the MRI images. Coefficient of determination R (2) values were computed for the comparison of circumferential and radial strain predictions produced by each methodology. There was good correspondence between the methodologies, with R (2) values of 0.78 for the radial strains of both hearts and from an R (2) = 0.81 and R (2) = 0.83 for the circumferential strains. The strain predictions were not statistically different (P ≤ 0.01). A series of sensitivity results indicated that the methodology was relatively insensitive to alterations in image intensity, random image noise, and alterations in fiber structure. This study demonstrated that warping was able to provide strain predictions of systolic contraction of the LV consistent with those provided by tagged MRI Warping.

Automatic 4-D Registration in Dynamic MR Renography

Dynamic contrast-enhanced 4-D MR renography has the potential for broad clinical applications, but suffers from respiratory motion that limits analysis and interpretation. Since each examination yields at least over 10 - 20 serial 3-D images of the abdomen, manual registration is prohibitively labor-intensive. Besides in-plane motion and translation, out-of-plane motion and rotation are observed in the image series. In this paper, a novel robust and automated technique for removing out-of-plane translation and rotation with sub-voxel accuracy in 4-D dynamic MR images is presented. The method was evaluated on simulated motion data derived directly from a clinical patients data. The method was also tested on 24 clinical patient kidney data sets. Registration results were compared with a mutual information method, in which differences between manually co-registered time-intensity curves and tested time-intensity curves were compared. Evaluation results showed that our method agreed well with these ground truth data.

True 4D Image Denoising on the GPU

The use of image denoising techniques is an important part of many medical imaging applications. One common application is to improve the image quality of low-dose (noisy) computed tomography (CT) data. While 3D image denoising previously has been applied to several volumes independently, there has not been much work done on true 4D image denoising, where the algorithm considers several volumes at the same time. The problem with 4D image denoising, compared to 2D and 3D denoising, is that the computational complexity increases exponentially. In this paper we describe a novel algorithm for true 4D image denoising, based on local adaptive filtering, and how to implement it on the graphics processing unit (GPU). The algorithm was applied to a 4D CT heart dataset of the resolution 512  × 512  × 445  × 20. The result is that the GPU can complete the denoising in about 25 minutes if spatial filtering is used and in about 8 minutes if FFT-based filtering is used. The CPU implementation requires several days of processing time for spatial filtering and about 50 minutes for FFT-based filtering. The short processing time increases the clinical value of true 4D image denoising significantly.

Left ventricular border tracking using cardiac motion models and optical flow

The use of automated methods is becoming increasingly important for assessing cardiac function quantitatively and objectively. In this study, we propose a method for tracking three-dimensional (3-D) left ventricular contours. The method consists of a local optical flow tracker and a global tracker, which uses a statistical model of cardiac motion in an optical-flow formulation. We propose a combination of local and global trackers using gradient-based weights. The algorithm was tested on 35 echocardiographic sequences, with good results (surface error: 1.35 ± 0.46 mm, absolute volume error: 5.4 ± 4.8 mL). This demonstrates the method's potential in automated tracking in clinical quality echocardiograms, facilitating the quantitative and objective assessment of cardiac function.

Parameterization of real-time 3D speckle tracking framework for cardiac strain assessment

Cross-correlation based 3D speckle tracking algorithm can be used to automatically track myocardial motion on three dimensional real-time (RT3D) echocardiography. The goal of this study was to experimentally investigate the effects of different parameters associated with such algorithm to ensure accurate cardiac strain measurements. The investigation was performed on 10 chronic obstructive pulmonary disease RT3DE cardiac ultrasound images. The following two parameters were investigated: 1) the gradient threshold of the anisotropic diffusion pre-filtering and 2) the window size of the cross correlation template matching in the speckle tracking. Results suggest that the optimal gradient threshold of the anisotropic filter depends on the average gradient of the background speckle noise, and that an optimal pair of template size and search window size can be identified determines the cross-correlation level and computational cost.

Automated border detection in three-dimensional echocardiography: principles and promises

Several automated border detection approaches for three-dimensional echocardiography have been developed in recent years, allowing quantification of a range of clinically important parameters. In this review, the background and principles of these approaches and the different classes of methods are described from a practical perspective, as well as the research trends to achieve a robust method.

Online 3-D reconstruction of the right atrium from echocardiography data via a topographic cellular contour extraction algorithm

A computational method providing online, automated 3-D reconstruction of the right atrium of the human heart is presented in this paper. Endocardial boundaries were extracted from transesophageal ultrasound data by a topographic cellular contour extraction (TCCE) algorithm. The TCCE method was implemented on a continuous-time, analog, massively parallel processor, and on a digital serial processor. Processing speeds were 500 or 40 frames per second, depending on the hardware used. Extracted boundary point sets were rendered into a 3-D mesh and the volume of the cavity was quantified. Accuracy of volume quantification was validated on 60 in vitro static phantoms and 12 clinical subjects. For the clinical recordings, reference volumes were estimated from manually delineated endocardial boundaries. The average error of volume quantification by the TCCE method was 8% +/-5% ( n = 12), the average of the interobserver variability between two independent human experts was 5% +/-2% ( n = 6). Interactive planning of atrial septal defect closure in pediatric cardiology is presented as an example that demonstrates a potential clinical application of the method.

Computer-assisted venous thrombosis volume quantification

Venous thrombosis (VT) volume assessment, by verifying its risk of progression when anticoagulant or thrombolytic therapies are prescribed, is often necessary to screen life-threatening complications. Commonly, VT volume estimation is done by manual delineation of few contours in the ultrasound (US) image sequence, assuming that the VT has a regular shape and constant radius, thus producing significant errors. This paper presents and evaluates a comprehensive functional approach based on the combination of robust anisotropic diffusion and deformable contours to calculate VT volume in a more accurate manner when applied to freehand 2-D US image sequences. Robust anisotropic filtering reduces image speckle noise without generating incoherent edge discontinuities. Prior knowledge of the VT shape allows initializing the deformable contour, which is then guided by the noise-filtering outcome. Segmented contours are subsequently used to calculate VT volume. The proposed approach is integrated into a system prototype compatible with existing clinical US machines that additionally tracks the acquired images 3-D position and provides a dense Delaunay triangulation required for volume calculation. A predefined robust anisotropic diffusion and deformable contour parameter set enhances the system usability. Experimental results pertinence is assessed by comparison with manual and tetrahedron-based volume computations, using images acquired by two medical experts of eight plastic phantoms and eight in vitro VTs, whose independently measured volume is the reference ground truth. Results show a mean difference between 16 and 35 mm(3) for volumes that vary from 655 to 2826 mm(3). Two in vivo VT volumes are also calculated to illustrate how this approach could be applied in clinical conditions when the real value is unknown. Comparative results for the two experts differ from 1.2% to 10.08% of the smallest estimated value when the image acquisition cadences are similar.

Region-based endocardium tracking on real-time three-dimensional ultrasound

Matrix-phased array transducers for real-time 3-D ultrasound enable fast, noninvasive visualization of cardiac ventricles. Typically, 3-D ultrasound images are semiautomatically segmented to extract the left ventricular endocardial surface at end-diastole and end-systole. Automatic segmentation and propagation of this surface throughout the entire cardiac cycle is a challenging and cumbersome task. If the position of the endocardial surface is provided at one or two time frames during the cardiac cycle, automated tracking of the surface over the remaining time frames could reduce the workload of cardiologists and optimize analysis of 3-D ultrasound data. In this paper, we applied a region-based tracking algorithm to track the endocardial surface between two reference frames that were manually segmented. To evaluate the tracking of the endocardium, the method was applied to 40 open-chest dog datasets with 484 frames in total. Ventricular geometry and volumes derived from region-based endocardial surfaces and manual tracing were quantitatively compared, showing strong correlation between the two approaches. Statistical analysis showed that the errors from tracking were within the range of interobserver variability of manual tracing. Moreover, our algorithm performed well on ischemia datasets, suggesting that the method is robust-to-abnormal wall motion. In conclusion, the proposed optical flow-based surface tracking method is very efficient and accurate, providing dynamic "interpolation" of segmented endocardial surfaces.

Sparse registration for three-dimensional stress echocardiography

Three-dimensional (3-D) stress echocardiography is a novel technique for diagnosing cardiac dysfunction. It involves evaluating wall motion of the left ventricle, by visually analyzing ultrasound images obtained in rest and in different stages of stress. Since the acquisitions are performed minutes apart, variabilities may exist in the visualized cross-sections. To improve anatomical correspondence between rest and stress, aligning the images is essential. We developed a new intensity-based, sparse registration method to retrieve standard anatomical views from 3-D stress images that were equivalent to the manually selected views in the rest images. Using sparse image planes, the influence of common image artifacts could be reduced. We investigated different similarity measures and different levels of sparsity. The registration was tested using data of 20 patients and quantitatively evaluated based on manually defined anatomical landmarks. Alignment was best using sparse registration with two long-axis and two short-axis views; registration errors were reduced significantly, to the range of interobserver variabilities. In 91% of the cases, the registration result was qualitatively assessed as better than or equal to the manual alignment. In conclusion, sparse registration improves the alignment of rest and stress images, with a performance similar to manual alignment. This is an important step towards objective quantification in 3-D stress echocardiography.

Enhancing obstetric and gynecology ultrasound images by adaptation of the speckle reducing anisotropic diffusion filter

OBJECTIVE

So far there is no ideal speckle reduction filtering technique that is capable of enhancing and reducing the level of noise in medical ultrasound (US) images, while efficiently responding to medical experts' validation criteria which quite often include a subjective component. This paper presents an interactive tool called evolutionary speckle reducing anisotropic diffusion filter (EVOSRAD) that performs adaptive speckle filtering on ultrasound B-mode still images. The medical expert runs the algorithm interactively, having a permanent control over the output, and guiding the filtering process towards obtaining enhanced images that agree to his/her subjective quality criteria.

METHODS AND MATERIAL

We employ an interactive evolutionary algorithm (IGA) to adapt on-line the parameters of a speckle reducing anisotropic diffusion (SRAD) filter. For a given input US image, the algorithm evolves the parameters of the SRAD filter according to subjective criteria of the medical expert who runs the interactive algorithm. The method and its validation are applied to a test bed comprising both real and simulated obstetrics and gynecology (OB/GYN) ultrasound images.

RESULTS

The potential of the method is analyzed in comparison to other speckle reduction filters: the original SRAD filter, the anisotropic diffusion, offset and median filters. Results obtained show the good potential of the method on several classes of OB/GYN ultrasound images, as well as on a synthetic image simulating a real fetal US image. Quality criteria for the evaluation and validation of the method include subjective scoring given by the medical expert who runs the interactive method, as well as objective global and local quality criteria.

CONCLUSIONS

The method presented allows the medical expert to design its own filters according to the degree of medical expertise as well as to particular and often subjective assessment criteria. A filter is designed for a given class of ultrasound images and for a given medical expert who will later use the respective filter in clinical practice. The process of designing a filter is simple and employs an interactive visualization and scoring stage that does not require image processing knowledge. Results show that filters tailored using the presented method achieve better quality scores than other more generic speckle filtering techniques.

Image filtering via generalized scale

In medical imaging, low signal-to-noise ratio (SNR) and/or contrast-to-noise ratio (CNR) often cause many image processing algorithms to perform poorly. Postacquisition image filtering is an important off-line image processing approach widely employed to enhance the SNR and CNR. A major drawback of many filtering techniques is image degradation by diffusing/blurring edges and/or fine structures. In this paper, we introduce a scale-based filtering method that employs scale-dependent diffusion conductance to perform filtering. This approach utilizes novel object scale information via a concept called generalized scale, which imposes no shape, size, or anisotropic constraints unlike previously published ball scale-based filtering strategies. The object scale allows us to better control the filtering process by constraining smoothing in regions with fine details and in the vicinity of boundaries while permitting effective smoothing in the interior of homogeneous regions. A new quantitative evaluation strategy that captures the SNR to CNR trade-off behavior of filtering methods is presented. The evaluations based on the Brainweb data sets show superior performance of generalized scale-based diffusive filtering over two existing methods, namely, ball scale-based and nonlinear complex diffusion processes. Qualitative experiments based on both phantom and patient magnetic resonance images demonstrate that the generalized scale-based approach leads to better preservation of fine details and edges.

Real-time tracking and shape analysis of atrial septal defects in 3D echocardiography

RATIONALE AND OBJECTIVES

Real-time cardiac ultrasound (US) allows monitoring the heart motion during intracardiac beating heart procedures. Our application assists pediatric atrial septal defect (ASD) closure techniques using real-time 3D US guidance and rigid instruments. ASD tracking is also an important tool for facilitating systematic clinical studies of the dynamic behavior of the intra-atrial communication. One major image processing challenge is associated with the required processing of information at high frame rate, especially given the low image quality.

MATERIALS AND METHODS

We present an optimization scheme for a block flow technique, which combines the probability-based velocity computation for an entire block (a 3D volume centered on the ASD) with cyclic template matching. The adapted similarity imposes constraints both locally (from frame to frame) to conserve energy, and globally (from a reference template) to minimize cumulative errors. The algorithm is optimized for fast and reliable results. For tests, we use three intra-operational 4D ultrasound sequences of clinical infant beating hearts with ASD.

RESULTS

Computing velocity at the block level with an optimized scheme, our technique tracks ASD motion at a frequency of 60 frames/s on clinical 4D datasets. Results are stable and accurate for changes in resolution and block size. In particular, we show robust real-time tracking and preliminary segmentation results of the ASD shape, size and orientation as a function of time.

CONCLUSIONS

We present an optimized block flow technique for real-time tracking of ASD to assist in minimally invasive beating heart surgery. Our method proposes the standard use of references for processing repetitive data. This paper represents, to our knowledge, the first study on the dynamic morphology of ASD that takes into account the angular effect introduced by the slanted position of the intra-atrial communication with respect to the US probe.

Application of a new segmentation tool based on interactive simplex meshes to cardiac images and pulmonary MRI data

RATIONALE AND OBJECTIVES

Medical image segmentation is still very time consuming and is therefore seldom integrated into clinical routine. Various three-dimensional (3D) segmentation approaches could facilitate the work, but they are rarely used in clinical setups because of complex initialization and parametrization of such models.

MATERIALS AND METHODS

We developed a new semiautomatic 3D-segmentation tool based on deformable simplex meshes. The user can define attracting points in the original image data. The new deformation algorithm guarantees that the surface model will pass through these interactively set points. The user can directly influence the evolution of the deformable model and gets direct feedback during the segmentation process.

RESULTS

The segmentation tool was evaluated for cardiac image data and magnetic resonance imaging lung images. Comparison with manual segmentation showed high accuracy. Time needed for delineation of the various structures could be reduced in some cases. The model was not sensitive to noise in the input data and model initialization.

CONCLUSIONS

The tool is suitable for fast interactive segmentation of any kind of 3D or 3D time-resolved medical image data. It enables the clinician to influence a complex 3D-segmentation algorithm and makes this algorithm controllable. The better the quality of the data, the less interaction is required. The tool still works when the processed images have low quality.

A global approach for solving evolutive heat transfer for image denoising and inpainting

This paper proposes an alternative to partial differential equations (PDEs) for solving problems in computer vision based on evolutive heat transfer. Traditionally, the method for solving such physics-based problems is to discretize and solve a PDE by a purely mathematical process. Instead of using the PDE, we propose to use the global heat principle and to decompose it into basic laws. We show that some of these laws admit an exact global version since they arise from conservative principles. We also show that the assumptions made about the other basic Iaws can be made wisely, taking into account knowledge about the problem and the domain. The numerical scheme is derived in a straightforward way from the modeled problem, thus providing a physical explanation for each step in the solution. The advantage of such an approach is that it minimizes the approximations made during the whole process and it modularizes it, allowing changing the application to a great number of problems. We apply the scheme to two applications: image denoising and inpainting which are modeled with heat transfer. For denoising, we propose a new approximation for the conductivity coefficient and we add thin lines to the features in order to block diffusion.

Ultrasound image segmentation: a survey

This paper reviews ultrasound segmentation paper methods, in a broad sense, focusing on techniques developed for medical B-mode ultrasound images. First, we present a review of articles by clinical application to highlight the approaches that have been investigated and degree of validation that has been done in different clinical domains. Then, we present a classification of methodology in terms of use of prior information. We conclude by selecting ten papers which have presented original ideas that have demonstrated particular clinical usefulness or potential specific to the ultrasound segmentation problem.

An electromechanical model of the heart for image analysis and simulation

This paper presents a new three-dimensional electromechanical model of the two cardiac ventricles designed both for the simulation of their electrical and mechanical activity, and for the segmentation of time series of medical images. First, we present the volumetric biomechanical models built. Then the transmembrane potential propagation is simulated, based on FitzHugh-Nagumo reaction-diffusion equations. The myocardium contraction is modeled through a constitutive law including an electromechanical coupling. Simulation of a cardiac cycle, with boundary conditions representing blood pressure and volume constraints, leads to the correct estimation of global and local parameters of the cardiac function. This model enables the introduction of pathologies and the simulation of electrophysiology interventions. Moreover, it can be used for cardiac image analysis. A new proactive deformable model of the heart is introduced to segment the two ventricles in time series of cardiac images. Preliminary results indicate that this proactive model, which integrates a priori knowledge on the cardiac anatomy and on its dynamical behavior, can improve the accuracy and robustness of the extraction of functional parameters from cardiac images even in the presence of noisy or sparse data. Such a model also allows the simulation of cardiovascular pathologies in order to test therapy strategies and to plan interventions.

SPASM: A 3D-ASM for segmentation of sparse and arbitrarily oriented cardiac MRI data

A new technique (SPASM) based on a 3D-ASM is presented for automatic segmentation of cardiac MRI image data sets consisting of multiple planes with arbitrary orientations, and with large undersampled regions. Model landmark positions are updated in a two-stage iterative process. First, landmark positions close to intersections with images are updated. Second, the update information is propagated to the regions without image information, such that new locations for the whole set of the model landmarks are obtained. Feature point detection is performed by a fuzzy inference system, based on fuzzy C-means clustering. Model parameters were optimized on a computer cluster and the computational load distributed by grid computing. SPASM was applied to image data sets with an increasing sparsity (from 2 to 11 slices) comprising images with different orientations and stemming from different MRI acquisition protocols. Segmentation outcomes and calculated volumes were compared to manual segmentation on a dense short-axis data configuration in a 3D manner. For all data configurations, (sub-)pixel accuracy was achieved. Performance differences between data configurations were significantly different (p<0.05) for SA data sets with less than 6 slices, but not clinically relevant (volume differences<4 ml). Comparison to results from other 3D model-based methods showed that SPASM performs comparable to or better than these other methods, but SPASM uses considerably less image data. Sensitivity to initial model placement proved to be limited within a range of position perturbations of approximately 20 mm in all directions.

An implementation of image sharpening based on morphological operations for ubiquitous echo

Ubiquitous echo is a portable ultrasound imaging equipment. We discuss an image sharpening method based on geometrical information by mathematical morphology with double structuring element (DSE) for on-line processing on ubiquitous echo. The sharpening method improves the contrast of tissue boundaries without speckle emphasis. The computational complexity of the morphological operations is reduced by chain rule of the operations and decomposition of the DSE not to delay refreshing the ultrasound moving image.

A method for the dynamic analysis of the heart using a Lyapounov based denoising algorithm

Heart tracking in ultrasound sequences is a difficult task due to speckle noise, low SNR and lack of contrast. Therefore it is usually difficult to obtain robust estimates of the heart cavities since feature detectors produce a large number of outliers. This paper presents an algorithm which combines two main operations: i) a novel denoising algorithm based on the Lyapounov equation and ii) a robust tracker, recently proposed by the authors, based on a model of the outlier features. Experimental results are provided, showing that the proposed algorithm is computationally efficient and leads to accurate estimates of the left ventricle during the cardiac cycle.

Ultrasound image segmentation using spectral clustering

Segmentation of ultrasound images is necessary in a variety of clinical applications, but the development of automatic techniques is still an open problem. Spectral clustering techniques have recently become popular for data and image analysis. In particular, image segmentation has been proposed via the normalized cut (NCut) criterion. This article describes an initial investigation to determine the suitability of such segmentation techniques for ultrasound images. The adaptation of the NCut technique to ultrasound is described first. Segmentation is then performed on simulated ultrasound images. Tests are also performed on abdominal and fetal images with the segmentation results compared to manual segmentation. The success of the segmentation on these test cases warrants further research into NCut-based segmentation of ultrasound images.

MRI diffusion-based filtering: a note on performance characterisation

Frequently MRI data is characterised by a relatively low signal to noise ratio (SNR) or contrast to noise ratio (CNR). When developing automated Computer Assisted Diagnostic (CAD) techniques the errors introduced by the image noise are not acceptable. Thus, to limit these errors, a solution is to filter the data in order to increase the SNR. More importantly, the image filtering technique should be able to reduce the level of noise, but not at the expense of feature preservation. In this paper we detail the implementation of a number of 3D diffusion-based filtering techniques and we analyse their performance when they are applied to a large collection of MR datasets of varying type and quality.

Combinative multi-scale level set framework for echocardiographic image segmentation

In the automatic segmentation of echocardiographic images, a priori shape knowledge has been used to compensate for poor features in ultrasound images. This shape knowledge is often learned via an off-line training process, which requires tedious human effort and is highly expertise-dependent. More importantly, a learned shape template can only be used to segment a specific class of images with similar boundary shape. In this paper, we present a multi-scale level set framework for segmentation of endocardial boundaries at each frame in a multiframe echocardiographic image sequence. We point out that the intensity distribution of an ultrasound image at a very coarse scale can be approximately modeled by Gaussian. Then we combine region homogeneity and edge features in a level set approach to extract boundaries automatically at this coarse scale. At finer scale levels, these coarse boundaries are used to both initialize boundary detection and serve as an external constraint to guide contour evolution. This constraint functions similar to a traditional shape prior. Experimental results validate this combinative framework.

Deformable biomechanical models: Application to 4D cardiac image analysis

This article describes a methodology for creating a generic volumetric biomechanical model from different image modalities and segmenting time series of medical images using this model. The construction of such a generic model consists of three stages: geometric meshing, non-rigid deformation of the mesh in images of various modalities, and image-to-mesh information mapping through rasterization. The non-rigid deformation stage, which relies on a combination of global and local deformations, can then be used to segment time series of images, e.g. cine MRI or gated SPECT cardiac images. We believe that this type of deformable biomechanical model can play an important role in the extraction of useful quantitative local parameters of cardiac function. The biomechanical model of the heart will be coupled with an electrical model of another collaborative project in order to simulate and analyze a larger class of pathologies.

3-D active appearance models: segmentation of cardiac MR and ultrasound images

A model-based method for three-dimensional image segmentation was developed and its performance assessed in segmentation of volumetric cardiac magnetic resonance (MR) images and echocardiographic temporal image sequences. Comprehensive design of a three-dimensional (3-D) active appearance model (AAM) is reported for the first time as an involved extension of the AAM framework introduced by Cootes et al. The model's behavior is learned from manually traced segmentation examples during an automated training stage. Information about shape and image appearance of the cardiac structures is contained in a single model. This ensures a spatially and/or temporally consistent segmentation of three-dimensional cardiac images. The clinical potential of the 3-D AAM is demonstrated in short-axis cardiac MR images and four-chamber echocardiographic sequences. The method's performance was assessed by comparison with manually identified independent standards in 56 clinical MR and 64 clinical echo image sequences. The AAM method showed good agreement with the independent standard using quantitative indexes of border positioning errors, endo- and epicardial volumes, and left ventricular mass. In MR, the endocardial volumes, epicardial volumes, and left ventricular wall mass correlation coefficients between manual and AAM were R2 = 0.94, 0.97, 0.82, respectively. For echocardiographic analysis, the area correlation was R2 = 0.79. The AAM method shows high promise for successful application to MR and echocardiographic image analysis in a clinical setting.