Skip to main content


Swipe om te navigeren naar een ander artikel

01-06-2011 | ICIN | Uitgave 6/2011 Open Access

Netherlands Heart Journal 6/2011

Automated analysis of three-dimensional stress echocardiography

Netherlands Heart Journal > Uitgave 6/2011
K. Y. E. Leung, M. van Stralen, M. G. Danilouchkine, G. van Burken, M. L. Geleijnse, J. H. C. Reiber, N. de Jong, A. F. W. van der Steen, J. G. Bosch


Stress echocardiography (stress echo) is widely applied for evaluating cardiac function [1, 2]. Motion of the left ventricle is assessed by comparing ultrasound images taken in rest and in stress stages. In the past decade, real-time three-dimensional (3D) imaging has become available for quantitative and objective assessment of cardiac function [3, 4]. 3D imaging offers the opportunity to overcome the limitations of traditional two-dimensional stress echo [5]. With 3D imaging, non-foreshortened anatomical views can be defined since different view choices can be made after acquisition. Also, it provides better possibilities for quantifying the true 3D wall motion. However, due to suboptimal image quality, lower spatial and temporal resolution, and the sheer amount of data, quantitative analysis by hand is tedious and subjective.
Automated analysis methods may allow objective and faster assessment of clinical parameters such as left ventricular volume, ejection fraction and wall motion. In practice, the clinician makes a diagnosis using expert knowledge, gathered by analysing many patient images [6, 7]. To develop successful automated methods which can deal with different pathologies and varying image quality, it makes sense to take expert knowledge into account.
In this paper, we present automated methods for 3D stress echo which make use of knowledge gained from databases of patient data. The knowledge is encompassed in so-called statistical models. Using Principal Component Analysis, patient variability can be represented concisely as mathematical descriptors that form a vector of orthogonal parameters [8]. The parameters to be modelled (for example, spatial coordinates of the endocardial border) are compared with the average parameters and typical variations in the patient database. Essentially, the methods are geared towards fitting the image data at hand to the model. In other words, the variation is estimated according to a set of extra criteria (such as image intensity values, ‘smoothness’ of the endocardial border, physically plausible motion, etc.).
Our proposed analysis scheme consists of three steps: identification of correct anatomical views in 3D data, detection of endocardial borders from which motion can be derived, and automated classification of motion. In addition, we present a software package which we have developed dedicated to analysis of 3D stress echo.

Identification of anatomical views

First we describe a method for detecting four-chamber, two-chamber, and short-axis views in a 3D image. These views are modelled via appearance (image intensity values) and pose parameters (translation, rotation, and scaling). These parameters are identified by fitting the database model to the image at hand within an optimisation framework [9].
The views in a rest image can be used to extract the anatomically corresponding views in a stress image of the same patient. This is achieved via image registration [10], so that rest and stress images are aligned via minimisation of an objective criterion (in this case the correlation between image intensity values) [11]. Qualitative and quantitative evaluation in 20 end-diastolic and 20 end-systolic images shows that the views found are in many cases better (18%) or similar (75%) to the manual identification of views. Anatomically aligned views are thus obtained automatically, providing more possibilities for accurate visualisation and quantification.

Quantification of volume and wall motion

Detection of endocardial borders provides important clinical parameters such as volume, ejection fraction, and wall motion. Here, the active appearance model technique [8, 12] is adapted to detect borders in an end-diastolic 3D echocardiogram. The endocardial borders are modelled via spatial coordinates, distributed in a cylindrical and spherical representation which is suited to the shape of the left ventricle (Fig. 1). Image intensity values are mapped to a normal (Gaussian) distribution for appropriate statistical analysis [12]. The model represents variations of ventricular shapes and the typical appearance of the ventricle and myocardium in echocardiograms, including the typical artifacts. The model is trained for image analysis by estimating the relation between model parameters and image change via regression analysis: in the actual matching, the intensity difference between model and image is used to update the model’s parameters and drive the model closer to the image (Fig. 2). Evaluation on 99 patient images shows a successful matching in 91% of cases, with a median surface error of 2.65 mm (average 2.91 mm, standard deviation 1.03 mm)[13].
The borders throughout the cardiac cycle are obtained by tracking the motion of the cardiac wall. A 3D optical flow algorithm is applied, which uses spatial and temporal intensity gradients to track the endocardial border between two time frames. To ensure an overall temporal continuity, the tracking method takes into account physically plausible motion, learned statistically from a patient database. Typical image artifacts, such as echo drop-out and near-field artifacts, are detected and dealt with appropriately using an expectation maximisation based classification approach [14]. The method is capable of tracking the endocardial borders accurately (in 35 patients: surface errors: 1.2 ± 0.5 mm, volume errors: 1.4 ± 6.7 ml, ejection fraction errors: 0.9 ± 4.8%).

Automated motion classification

Statistical models of the endocardial border motion throughout the cardiac cycle can be used to discriminate between normal and pathological motion patterns. The statistical model is transformed using the Orthomax rotation criterion, so that the local wall motion can be represented using fewer descriptors (Fig. 3) [6]. An automated classification method is then applied to the descriptors. Wall motion classification is demonstrated in 129 two-dimensional echocardiographic sequences. Using these models, local motion abnormalities can be detected accurately in 77–91% of the cases. In principle, the method can be applied directly to 3D borders; the evaluation is a subject of further investigation.

Tools for clinical practice

To enable 3D stress echocardiography in routine clinical practice, we have also developed a software package dedicated to 3D stress echo [15]. The essential functionality for analysing 3D stress echo is investigated: side-by-side, synchronised display of anatomically aligned rest and stress images (Fig. 4). Two expert observers analysed 34 noncontrast clinically available rest and stress images visually. Interobserver agreement for segmental myocardial ischaemia is better using the dedicated software (96%) than commercially available software (79%) without these functionalities.


In this paper, we have presented automated methods for identification of anatomical views, detection of endocardial borders, and classification of wall motion abnormalities. The proposed methods make use of statistical modelling of expert knowledge from patient databases. Considerable advances have been made to allow a more objective and accurate analysis. In particular, these methods are better capable of handling typical image artifacts and following expert tracing conventions. We have also demonstrated the advantages of using dedicated software for the analysis of 3D stress echo in the clinical setting. Based on this research, opportunities remain for further development, e.g. in the analysis of contrast echocardiograms, technical enhancement of image quality and spatiotemporal resolution, and improvement of clinical workflow. The developed tools can provide useful quantitative and objective parameters to help the clinical expert in the diagnosis of left ventricular function.


This research is supported by the Dutch Technology Foundation STW (grant 06666), applied science division of NWO and the Technology Program of the Ministry of Economic Affairs, the Netherlands

Open Access

This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Open AccessThis is an open access article distributed under the terms of the Creative Commons Attribution Noncommercial License (https://​creativecommons.​org/​licenses/​by-nc/​2.​0), which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.

Onze productaanbevelingen

Netherlands Heart Journal

Het Netherlands Heart Journal wordt uitgegeven in samenwerking met de Nederlandse Vereniging voor Cardiologie en de Nederlandse Hartstichting. Het tijdschrift is Engelstalig en wordt gratis beschikbaa ...

Over dit artikel

Andere artikelen Uitgave 6/2011

Netherlands Heart Journal 6/2011 Naar de uitgave

Rhythm Puzzle - Answer

Answer to the rhythm puzzle

Rhythm Puzzle - Question

A simple broad complex tachycardia?