- Technical advance
- Open Access
- Open Peer Review
Automatic measurement of aortic annulus diameter in 3-dimensional Transoesophageal echocardiography
- Jørn Bersvendsen^{1, 2, 3}Email author,
- Jan O Beitnes^{4},
- Stig Urheim^{4},
- Svend Aakhus^{4} and
- Eigil Samset^{1, 2, 3}
https://doi.org/10.1186/1471-2342-14-31
© Bersvendsen et al.; licensee BioMed Central Ltd. 2014
- Received: 6 August 2013
- Accepted: 6 June 2014
- Published: 8 September 2014
Abstract
Background
Transcatheter aortic valve implantation involves percutaneously implanting a biomechanical aortic valve to treat severe aortic stenosis. In order to select a proper device, precise sizing of the aortic valve annulus must be completed.
Methods
In this paper, we describe a fully automatic segmentation method to measure the aortic annulus diameter in patients with aortic calcification, operating on 3-dimensional transesophageal echocardiographic images. The method is based on state estimation of a subdivision surface representation of the left ventricular outflow tract and aortic root. The state estimation is solved by an extended Kalman filter driven by edge detections normal to the subdivision surface.
Results
The method was validated on echocardiographic recordings of 16 patients. Comparison against two manual measurements showed agreements (mean ±SD) of -0.3±1.6 and -0.2±2.3 mm for perimeter-derived diameters, compared to an interobserver agreement of -0.1±2.1 mm.
Conclusions
With this study, we demonstrated the feasibility of an efficient and fully automatic measurement of the aortic annulus in patients with aortic disease. The algorithm robustly measured the aortic annulus diameter, providing measurements indistinguishable from those done by cardiologists.
Keywords
- Segmentation
- Subdivision surface
- 3-Dimensional Echocardiography
- Aortic valve
- Transcatheter aortic valve implantation
Background
Transcatheter Aortic Valve Implantation (TAVI) involves percutaneously implanting a biomechanical aortic valve to treat severe aortic stenosis. Because of its minimally invasive nature, TAVI is a viable alternative for patients who are at too high risk to undergo conventional surgical aortic valve replacement.
Precise sizing of the aortic annulus prior to TAVI is required for determining procedure eligibility and for selecting the correct implant size and type. Errors in prosthesis sizing may lead to complications during or after the procedure, such as Paravalvular Aortic Regurgitation (PAR) [1].
In current clinical practice, measurement of the annulus diameter before TAVI is typically done by 2D transthoracic echocardiography, 2D Transoesophageal Echocardiography (TEE) or Multi-Slice Computed Tomography (MSCT).
It has been shown that sizing based on MSCT, as opposed to 2D TEE, results in fewer instances of post operational PAR [2], as 2D modalities can fail to accurately describe the 3D structure of the aortic valve [3–5]. Strong correlations between 3D TEE and MSCT measurements of the annulus diameter [4] indicate feasibility of similar results for a method based on 3D TEE.
We propose an algorithm for automatic annulus measurements operating on 3D TEE images, using a real-time volumetric tracking and segmentation framework presented by Orderud et al. [6, 7]. The framework uses an extended Kalman filter to solve a state space estimation formulation of the segmentation problem, and has been applied on the left ventricle.
In this paper, we apply the same framework to model the left ventricular outflow tract and aortic root. We propose a two-stage approach by performing segmentation based on a stiff and deformable surface sequentially. Combined with assimilation of forward and backward tracking, we obtain a fully automatic measurement of the aortic annulus diameter in 3D TEE images.
We validated our results by comparing automatic measurements of 16 recordings to manual measurements made by two cardiologists.
Methods
Segmentation
The method presented here is an application of a previously presented real-time volumetric segmentation framework, operating on deformable subdivision surfaces [6, 7]. The segmentation is represented as a state estimation problem and solved with an extended Kalman filter.
Surface model
The surface is deformed locally by displacing each control point in the direction normal to the cylinder long axis, to maximize deformation per degree of freedom. The local deformation transform is denoted T _{ l }(p;x _{ l }), where x _{ l } is the state vector of local deformations.
where x_{ g }=[t _{ x }, t _{ y }, t _{ z }, s, θ _{ x },θ _{ y },θ _{ z }]^{⊤} is the global transform state vector, allows for translation, scaling and rotation of the model. The composite transform is given by T(p;x)=T _{ g }(T _{ l }(p;x _{ l });x _{ g }) where ${\mathbf{\text{x}}}^{\top}=\left[{\mathbf{\text{x}}}_{g}^{\top},\phantom{\rule{0.3em}{0ex}}{\mathbf{\text{x}}}_{l}^{\top}\right]$ is the state vector.
The aortic annulus plane is represented by a disc placed in the middle of the surface model. The disc shares the same global transform T _{ g } but is not deformable.
Motion model
where A is a diagonal matrix specifying the regularization strength for each state. Note that 0≤a _{ ij }≤1 ensures stability. The diagonal elements of A were chosen separately for translation, scaling, rotation and deformation states.
where Q _{0} is the process noise covariance matrix. A low noise value will decrease the prediction covariance relative to the measurement covariance, which in turn will make the Kalman filter rely more on the previous estimate than the detected edges. Q _{0} therefore functions as a fairness parameter.
Edge detection
300 evenly distributed points on the model surface are defined. For the valve disc, 40 edge points are defined.
where $\widehat{\mathbf{\text{x}}}={\widehat{\mathbf{\text{x}}}}_{k|k-1}$.
Edge displacements are detected by searching along n _{ g } around p _{ g } using the least mean squares fit to an intensity step or peak function. Outlier edges are rejected based on the intensity step function height and differences between neighboring edges. The capture range is determined by the search length along n _{ g }.
Each measured edge displacement v _{ i } has an associated measurement noise with estimated variance r _{ i } which is the sum of squared deviations in the intensity fit. The variance estimates are normalized such that $\sum _{i}{r}_{i}={r}_{\mathit{\text{edge}}}$.
Measurement update
This computation is efficient as it does not require inversion of matrices with size dependent on the number of measurements.
Forward and backward tracking
A common problem with segmentation of time-series is that the segmentation lags behind the recording. We solve this by tracking forward and backward in time.
This bidirectional tracking makes the segmentation robust to significant movement of the LVOT and aortic root during the cardiac cycle.
Two-phase segmentation
The described tracking algorithm is run in two passes; stiff segmentation and deformable segmentation.
Stiff segmentation
In the first pass, the subdivision surface is made stiff by removing all deformation states. The initial mesh is then oriented along an estimated LVOT long axis, derived from the ultrasound recorded roll angle.
The Kalman filter iterates over each consecutive frame once to ensure rough convergence. This is then repeated for the actual segmentation. A simplified motion model is used where x _{0,k } = x _{0} and Q _{0,k } = Q _{0} are constant during the cardiac cycle. Q _{0} was chosen to be a diagonal matrix of process noise standard deviations.
The resulting global pose states ${\widehat{\text{x}}}_{\text{stiff}},k$, with estimated covariance P _{stiff},k, aligns the subdivision surface to the recording for each frame. This captures the global movement of the LVOT and aortic root during the cardiac cycle.
Deformable segmentation
After stiff segmentation, the deformation states are reintroduced. The state vectors ${\widehat{x}}_{\text{stiff}},k$ and covariance matrices P _{stiff},k from the stiff segmentation are used for x _{0,k } and Q _{0,k } in the motion model. Bidirectional tracking is performed over a single heart cycle.
Different prediction parameters A ar e used in the stiff and deformable phases to reflect the increased confidence of the global transform states after stiff segmentation. For deformable segmentation, the regularization strength A is increased for the these states, ensuring that the movement of the aortic structure is tracked.
Automatic annulus measurements
The aortic annulus is extracted by the intersection of the deformed surface model and the aortic valve disk. An ellipse is fitted to the intersection points by least mean squares optimization and the major and minor axes, area and perimeter are extracted.
The mid systolic frame was defined as the frame with maximum detected aortic annulus area.
Transoesophageal echocardiography
Acquisition
Patient characteristics
Age ^{∗}, yr | 69±16 |
Sex ^{∗}, male/female | 7/6 |
Disease, n | |
Aortic stenosis | 12 |
Aortic insufficiency | 2 |
Normal | 2 |
Aortic morphology, n | |
Tricuspid | 14 |
Bicuspid | 2 |
LV EF < 50 % ^{∗}, n | 3 |
Manual measurements
Sensitivity analysis
To assess the robustness of the algorithm with respect to the roll-angle derived initial LVOT long axis estimate, the following test was carried out. For a single recording, the initial state x _{0} prior to stiff segmentation was randomly perturbed 500 times, and the resulting perimeter-derived annulus diameter distributions were recorded. This was done separately for perturbations to translation and rotation.
Translation The perturbations were on the form d n/||n|| where the elements of $\mathbf{\text{n}}\in {\mathbb{R}}^{3}$ were uniformely distributed. The analysis was repeated for each d=2,4,…,12 mm, which is within the capture range of ± 14 mm.Rotation For each iteration, a vector n/||n|| was generated, where the elements of $\mathbf{\text{n}}\in {\mathbb{R}}^{3}$ were uniformely distributed. The initial model was then rotated an angle ϕ about the line v(t)=v _{0}+n t where v _{0} is the initial annulus center. The analysis was repeated for each ϕ=1,2,…,15 deg.
To assess robustness with respect to scaling, the segmentation was repeated for 100 linearly spaced initial annulus diameters D _{0}∈[ 15,35] mm.
Statistical analysis
Agreement between the automatic method and the two manual observers was analyzed using Blant and Altmans method and two-way absolute agreement intraclass correlation coefficients. All statistical analysis was performed using IBM SPSS Statistics for Windows, Version 20.0 (IBM Corp., Armonk, NY, USA).
Results
The described algorithm was executed on all 16 3D TEE datasets. The segmentation time was (mean ±SD) 9.9±7.3 s on a standard laptop.
Sensitivity analysis
Figure 3 shows the results of the sensitivity analysis for translational and rotational perturbations in the initial state. The sample standard deviations of the resulting diameter distributions were ≤0.29, ≤0.28 and 0.10 mm for the translation, rotation and scaling sensitivity tests respectively.
Comparison of manual and automatic measurements
Automatic and manual measurements of arotic annulus diameters
D _{auto} | D _{1} | D _{2} | |
---|---|---|---|
Minimum | 24.3 ± 3.0 | ||
Maximum | 26.8 ± 3.5 | ||
Sagittal | 23.9 ± 2.7 | 23.9 ± 3.1 | |
Coronal | 25.4 ± 3.3 | 25.3 ± 3.0 | |
Area-derived | 25.5 ± 3.2 | 24.9 ± 3.0 | 25.0 ± 3.1 |
Perimeter-derived | 25.6 ± 3.2 | 25.9 ± 3.2 | 25.8 ± 3.1 |
Comparison of automatic and manual measurements of arotic annulus diameter
Bias | Intraclass correlation | |
---|---|---|
Sagittal | ||
D _{1} versus D _{2} | 0.063 ± 1.5 | 0.88 |
Coronal | ||
D _{1} versus D _{2} | 0.063 ± 2.4 | 0.74 |
Perimeter-derived | ||
D _{ a u t o } versus D _{1} | -0.35 ± 1.6 | 0.87 |
D _{ a u t o } versus D _{2} | -0.23 ± 2.3 | 0.75 |
D _{1} versus D _{2} | 0.12 ± 2.1 | 0.78 |
Area-derived | ||
D _{ a u t o } versus D _{1} | 0.62 ± 1.7 | 0.85 |
D _{ a u t o } versus D _{2} | 0.46 ± 2.3 | 0.74 |
D _{1} versus D _{2} | -0.16 ± 2.1 | 0.77 |
The interobserver Intraclass Correlation Coefficients (ICC) were 0.78 and 0.77 for perimeter and area derived diameters respectively. The intraclass correlation coefficients between the automatic method and each of the manual observers were 0.87 and 0.75 for the perimeter derived diameters and 0.85 and 0.74 for the area derived diameters.
Discussion
Sensitivity analysis
Even for significant perturbations (± 12 mm translation, ± 15 deg rotation or ± 10 mm initial annulus diameter), the standard deviation of the automatic measurement was significantly lower than the interobserver variation. This indicates that the algorithm is robust with respect to errors in the initial LVOT long axis derived from the recorded roll angle, as well as the assumed initial annulus diameter.
Comparison of manual and automatic measurements
The algorithm performance was indistinguishable from human observers’ performance. The automatic algorithm successfully segmented the LVOT and aortic root and measured the aortic annulus diameter in all 16 images, with mean computation time 9.9 s. Interobserver correlation coefficient for the manual measurements was comparable to that reported by others [4]. The correlation and deviation between the automatic and each manual measurements were comparable to the interobserver reliability.
Perimeter-derived measurements showed the closest agreement with the manual observers. Since the perimeter-derived diameter changes the least during the cardiac cycle [8], this measurement should not be affected by errors in mid-systole detection. This is therefore a more robust measurement compared to area, major and minor diameters, and was chosen as the algorithm’s main output diameter.
The annulus plane is normally defined as the plane spanned by the hinge points of the three valve cusps. However, the hinge points are not explicitly defined in the described model. Because the disk representing the aortic valve is non-deformable and shares the same pose transform as the surface model, the detected annulus will align perpendicular to the long axis of the LVOT.
Since the hinge point plane and the perpendicular plane are closely aligned, we propose that a perpendicular plane is a good estimation of the anatomical annulus plane. In the rare cases where these planes are not aligned, we submit that a perpendicular plane is of clinical relevance since a prosthetic valve is more likely to align with the LVOT long axis than the native valve’s hinge points.
The largest absolute deviation between the two manual observers was 5.4 mm. Poor image quality, low frame rate (7.7 vps) and a wide sinus of Valsalva lead to the significant interobserver deviation. The largest absolute deviation between automatic and manual measurements was 5.5 mm. In this case, the automatic method grossly overestimated the annulus diameter, resulting from a very wide sinus of Valsalva combined with dropouts close to the annulus. Although these discrepancies would indicate a difference in device selection, it is extremely unlikely that these specific images would be used as the basis of device selection.
Two recordings had visible stitching artifacts. In both recordings all absolute deviations between the the automatic and manual measurements were ≤0.6 mm, indicating robustness against stitching artifacts.
Recently, a validation study of the first description of an automated aortic root modeling and quantification algorithm for 3D TEE images was published [9]. The study reported annulus diameter agreement (mean ±SD) of 1.1±1.3 and 3.6±2.3 mm for sagittal and coronal diameters respectively, which is comparable to our results. However, manual identification of peak systole and end diastole was required, and manual segmentation adjustments were needed in 23 of 69 TEE recordings. The reported interobserver variability was 0.2±0.56 and 0.0±0.61 mm. Although our presented algorithm is fully automatic and therefore has no interobserver variability, these values are comparable to our sensitivity analysis results. The reported computation and adjustment time was 2.3±0.6 minutes, which is significantly longer than our results.
This method is based on machine learning and statistical shape models [10]. However, these algorithms require a large database of recordings annotated with manual ground truth segmentations. The presented method is simpler and does not rely on a history of previous segmentations.
Within the presented framework, there are several quality measures available that can potentially be used to automatically identify poor segmentations, e.g. number of discarded edge profiles, deviation from segmentation surface to detected edges or the state covariance estimates. These should be further investigated with an available gold standard to create criteria that can automatically judge the segmentation quality.
This study used a limited sample size of 16 patients. Further studies with a larger number of patients should be performed.
Comparison of measurements in prospective 3D TEE images with ECG gated Multi-Slice Computed Tomography (MSCT) gold standard should be carried out to investigate if the algorithm can render MSCT superfluous for a significant portion of TAVI candidates.
Conclusions
With this study, we demonstrated the feasibility of an efficient and fully automatic measurement of the aortic annulus in patients with aortic disease. The algorithm robustly measured the aortic annulus diameter, providing measurements indistinguishable from those done by cardiologists.
Declarations
Acknowledgements
The authors gratefully acknowledge the Research Council of Norway for their research grant, as well as the support of Center for Cardiological Innovation, funded by the Research Council of Norway.
Authors’ Affiliations
References
- Zamorano JL, Badano LP, Bruce C, Chan K-L, Gonçalves A, Hahn RT, Keane MG, La Canna G, Monaghan MJ, Nihoyannopoulos P, Silvestry FE, Vanoverschelde J-L, Gillam LD, Vahanian A, Di Bello V, Buck T: EAE/ASE recommendations for the use of echocardiography in new transcatheter interventions for valvular heart disease. Eur J Echocardiogr : J Working Group Echocardiogr Eur Soc Cardiol. 2011, 12 (8): 557-584. 10.1093/ejechocard/jer086. doi:10.1093/ejechocard/jer086,View ArticleGoogle Scholar
- Jilaihawi H, Kashif M, Fontana G, Furugen A, Shiota T, Friede G, Makhija R, Doctor N, Leon MB, Makkar RR: Cross-sectional computed tomographic assessment improves accuracy of aortic annular sizing for transcatheter aortic valve replacement and reduces the incidence of paravalvular aortic regurgitation. J Am Coll Cardiol. 2012, 59 (14): 1275-1286. 10.1016/j.jacc.2011.11.045. doi:10.1016/j.jacc.2011.11.045,View ArticlePubMedGoogle Scholar
- Ng ACT, Delgado V, van der Kley F, Shanks M, van de Veire NRL, Bertini M, Nucifora G, van Bommel RJ, Tops LF, de Weger A, Tavilla G, de Roos A, Kroft LJ, Leung DY, Schuijf J, Schalij MJ, Bax JJ: Comparison of aortic root dimensions and geometries before and after transcatheter aortic valve implantation by 2- and 3-dimensional transesophageal echocardiography and multislice computed tomography. Circ Cardiovasc Imaging. 2010, 3 (1): 94-102. 10.1161/CIRCIMAGING.109.885152. doi:10.1161/CIRCIMAGING.109.885152,View ArticlePubMedGoogle Scholar
- Altiok E, Koos R, Schröder J, Brehmer K, Hamada S, Becker M, Mahnken AH, Almalla M, Dohmen G, Autschbach R, Marx N, Hoffmann R: Comparison of two-dimensional and three-dimensional imaging techniques for measurement of aortic annulus diameters before transcatheter aortic valve implantation. Heart (British Cardiac Soc). 2011, 97 (19): 1578-1584. 10.1136/hrt.2011.223974. doi:10.1136/hrt.2011.223974,View ArticleGoogle Scholar
- Muraru D, Badano LP, Vannan M, Iliceto S: Assessment of aortic valve complex by three-dimensional echocardiography: a framework for its effective application in clinical practice. Eur Heart J Cardiovascular Imaging. 2012, 13 (7): 541-555. 10.1093/ehjci/jes075. doi:10.1093/ehjci/jes075,View ArticleGoogle Scholar
- Orderud F: A Framework for real-time left ventricular tracking in 3D+T echocardiography, using nonlinear deformable contours and kalman filter based tracking. Computers in Cardiology, 2006. 2006, [http://ieeexplore.ieee.org/xpl/articleDetails.jsp?arnumber=4511804],Google Scholar
- Orderud F, Rabben SI: Real-time 3d segmentation of the left ventricle using deformable subdivision surfaces. Computer Vision and Pattern Recognition, 2008.CVPR 2008. IEEE Conference On. 2008, 1-8. doi:10.1109/CVPR.2008.4587442,View ArticleGoogle Scholar
- Hamdan A, Guetta V, Konen E, Goitein O, Segev A, Raanani E, Spiegelstein D, Hay I, Di Segni E, Eldar M, Schwammenthal E: Deformation dynamics and mechanical properties of the aortic annulus by 4-dimensional computed tomography: insights into the functional anatomy of the aortic valve complex and implications for transcatheter aortic valve therapy. J Am Coll Cardiol. 2012, 59 (2): 119-27. 10.1016/j.jacc.2011.09.045. doi:10.1016/j.jacc.2011.09.045,View ArticlePubMedGoogle Scholar
- Calleja A, Thavendiranathan P, Ionasec RI, Houle H, Liu S, Voigt I, Sai Sudhakar C, Crestanello J, Ryan T, Vannan MA: Automated quantitative 3-d modeling of the aortic valve and root by 3-d transesophageal echocardiography in normals, aortic regurgitation and, aortic stenosis: comparison to computed tomography in normals and clinical implications. Circ: Cardiovasc Imaging. 2012, doi:10.1161/CIRCIMAGING.112.976993,Google Scholar
- Ionasec RI, Voigt I, Georgescu B, Wang Y, Houle H, Vega-Higuera F, Navab N, Comaniciu D: Patient-specific modeling and quantification of the aortic and mitral valves from 4-d cardiac ct and tee. Med Imaging IEEE Trans on. 2010, 29 (9): 1636-1651. doi:10.1109/TMI.2010.2048756,View ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/14/31/prepub
Pre-publication history
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.