Finite-element-method (FEM) model generation of time-resolved 3D echocardiographic geometry data for mitral-valve volumetry
© Verhey et al; licensee BioMed Central Ltd. 2006
Received: 16 December 2005
Accepted: 03 March 2006
Published: 03 March 2006
Mitral Valve (MV) 3D structural data can be easily obtained using standard transesophageal echocardiography (TEE) devices but quantitative pre- and intraoperative volume analysis of the MV is presently not feasible in the cardiac operation room (OR). Finite element method (FEM) modelling is necessary to carry out precise and individual volume analysis and in the future will form the basis for simulation of cardiac interventions.
With the present retrospective pilot study we describe a method to transfer MV geometric data to 3D Slicer 2 software, an open-source medical visualization and analysis software package. A newly developed software program (ROIExtract) allowed selection of a region-of-interest (ROI) from the TEE data and data transformation for use in 3D Slicer. FEM models for quantitative volumetric studies were generated.
ROI selection permitted the visualization and calculations required to create a sequence of volume rendered models of the MV allowing time-based visualization of regional deformation. Quantitation of tissue volume, especially important in myxomatous degeneration can be carried out. Rendered volumes are shown in 3D as well as in time-resolved 4D animations.
The visualization of the segmented MV may significantly enhance clinical interpretation. This method provides an infrastructure for the study of image guided assessment of clinical findings and surgical planning. For complete pre- and intraoperative 3D MV FEM analysis, three input elements are necessary: 1. time-gated, reality-based structural information, 2. continuous MV pressure and 3. instantaneous tissue elastance. The present process makes the first of these elements available.
Volume defect analysis is essential to fully understand functional and geometrical dysfunction of but not limited to the valve. 3D Slicer was used for semi-automatic valve border detection and volume-rendering of clinical 3D echocardiographic data. FEM based models were also calculated.
A Philips/HP Sonos 5500 ultrasound device stores volume data as time-resolved 4D volume data sets. Data sets for three subjects were used. Since 3D Slicer does not process time-resolved data sets, we employed a standard movie maker to animate the individual time-based models and visualizations.
Calculation time and model size were minimized.
Pressures were also easily available. We speculate that calculation of instantaneous elastance may be possible using instantaneous pressure values and tissue deformation data derived from the animated FEM.
Mitral valve defects are a common source of severe mitral regurgitation. 3D echocardiography, especially the intraoperative TEE, allows a detailed assessment of the valve leaflets and commissures . It is currently available in most cardiac surgical operating rooms. In some centers, intraoperative 3D echocardiography is used to evaluate geometry and to plan surgical interventions prior to mitral valve remodeling surgery. However, quantitation of mitral valve geometry is limited to rather imprecise measures. Thus, the cardiac surgeon has no sophisticated, immediate, quantitative analysis of the preoperative 3D mitral valve geometry. Intraoperative quantitative analysis of the dynamic geometry of the mitral valve might provide new information upon which to base more precise, patient-specific planning of the surgical intervention. After the operation, assessment of the repair in a precise and quantitative fashion while the valve is fully operational may allow a better understanding of long term outcome predictors.
Because the mitral valve cannot be realistically described by a symmetric mathematical model , the modern approach consists of using a FEM mesh which can approximate the cardiac geometry [3–5]. A very good overview on this theme is given by J. Mackele .
Initial attempts at FEM in the heart have been carried out with 3D segmentation and tracking using sophisticated and expensive cardiac MRI [7–9]. MRI is impractical in the cardiac surgical operating room and is complicated by the fact that the mitral valve and the papillary muscles are active materials behaving differently during systole and diastole. An ideal model would provide material properties specific to each patient as first mentioned by McCulloch , but until now patient-specific modeling in the operating room has not been possible.
FEM modeling of 3D intraoperative echo data provides an excellent tool for incorporating material properties, volumetric data and boundary pressures to more accurately record, and then to simulate mitral valve dynamic performance. Accurate simulation will be the foundation of surgical planning. The limitation until now in applying FEM intraoperatively has been the technical complexity of this technique. The purpose of this study is to take the first step towards introducing FEM for mitral valve volumetry into the operating room environment. The goal is to facilitate transfer of geometric data from the 3D ultrasound data set into FEM models and there from to further visualization and simulation processing.
Materials and methods
General data set parameters
reference number in the article
number of time steps
resolution time [s]
Images were gated for both beat-to-beat variability and respiratory motion. All images were stripped of patient identifiers. The time dependent 3D data sets acquired from the Philips Sonos were arranged and stored on a Windows XP based standard PC with TomTec Echo-View© software  in a propeller-like geometry, in which the images intersect along the scanning axis . For mitral valve geometry reconstruction, a newly developed software program preprocessed the data by reformatting the data into parallel slices.
Time-resolved images of the rendered mitral valve geometry resulting from the 3D Slicer were transferred to a standard movie maker (Dreamweaver Fireworks 4)  running on a Fujitsu-Siemens laptop with Windows XP operation system in order to obtain animations for visualization purposes.
The quantification of distances and volumes can be carried out with 3D Slicer.
Data acquisition and processing times.
ultrasound acquisition [s]
pre-processing: reformatting, ROI-application [s]
FEM modelling with 3D Slicer [s]
post-processing moviemaker [s]
total processing time [s]
The Philips Sonos 5500 ultrasound system required around 11 minutes acquisition time per patient. The application of the reformatting software, the ROI algorithm with manual placement of the necessary landmarks took approximately 7 minutes per patient. In three patient data sets, conversion from TomTec Echo-View data to the 3D Slicer FEM model was carried out in 2 minutes for each time step of the heart cycle. On average, the processing time took 54 minutes per patient. Image acquisition and movie production on the above computer was 30 seconds per sequence or approximately 6 minutes per patient. Total time for the procedure was approximately 78 minutes per patient (see Table 2).
The primary intention of this study was to demonstrate the feasibility of transporting individual patient's mitral valve geometry data into a FEM model. Standard PC/workstation computer technology was utilized to accomplish the transfer from a common TEE-machine (Philips Sonos 5500).
The data acquisition software were commercially available TomTec Echo View© package and open-source 3D Slicer 2 software, plus recently developed components for data set reformatting. Accomplishing this transfer forms a foundation for quantitative approaches to intraoperative surgical planning and outcomes assessment in valvular reconstructive surgery.
This method would also be applicable to next generation "live 3D" systems, 3D ultrasound data sets obtained from matrix array transducers.
The total time required for acquisition to a completed FEM model was approximately 1hr 15 minutes and can be accomplished during the time period when the patient is being prepared for cardiopulmonary bypass (generally 1 to 1.5 hr). Thus the feasibility in terms of duration is demonstrated. The total time of processing could be surely decreased in the future if a single computer platform would be used. All used programs run e.g. on Windows XP systems. For this feasibility study we used a combination of three platforms just for our convenience and in order not to disturb the intervention itself.
In terms of procedure accuracy, reproducibility and duration, the primary limitation is the dependence of the newly developed software on manual entries of the four landmarks. Manual entries must be done for multiple frames within the TEE data sets and as a consequence the method involves some degree of inter- and intraobserver variability which is a general problem for ultrasonic imaging. The reproducibility of these measurements will require further study.
In Figure 3 the left of the three orthogonal planes at the bottom demonstrates an uncontinuous heart wall. This shows that the quality of the calibration of the transducer might be enhanced from the medical point of view. These data sets basically were taken during an intervention and not in an in-vitro experiment. From the technical point of view this doesn't limit the usability of the method described in the study.
A limitation of the present study is that it is focused only on the deployment of the transfer method. The tool for modelling is not within the scope of this article. Nevertheless, quantification the method makes it immediately available for mitral valve leaflet motion, shape and volume analysis. Instantaneous analysis has a number of potential applications for mitral valve function assessment and surgical planning. These applications require both comprehensive automated valve leaflet motion analysis as well as quantification of dynamic mechanical properties through a biomechanical FEM of the mitral valve region.
The scope of this study was to produce a prototype in which the feasibility of the method could be assessed. In a fully operational system, we postulate clinical applications such as enhanced/automated leaflet motion abnormality detection, assessment of regional relaxation which encompasses the entire ventricle, assessment and guidance of ventricular remodeling operations, and serial assessment of recovery of regional wall function post myocardial stunning
FEM meshes have been used for approximately 30 years [6, 15] in the analysis of many anatomical structures and organs such as major vessels [16, 17], heart valves  and ventricles [4, 19], lung , corneoscleral shell , plastic and reconstructive craniofacial surgery  and the femur . A FEM model can be created to determine the deformation of the mitral valve loaded by intraventricular pressure. Steady-state fluid dynamics and structural analyses can be carried out using commercial codes based on FEM . At a sequence of time-steps in the cardiac cycle, the valve leaflet can be considered to be a quasi-incompressible, transversely isotropic hyperelastic material based on the analysis of Feng . Until now, biomechanical cardiac FEM models have been based on a simplified ellipsoidal and cylindrical geometries  or asymmetric modelling . A FEM created in this way is not patient-specific and does not accurately represent precise regional deformations in a cardiac structure such as the mitral valve loaded by intraventricular pressure. The method described here will allow patient specificity and the precise representation of leaflet deformation. Precise deformation data is required to allow elastance determinations critical for simulation of the response of the tissue to virtual interventions.
For complete intraoperative 3D mitral valve finite element analysis, three input elements are necessary: 1. time-gated, reality-based structural information, 2. continuous left ventricle pressure, and 3. instantaneous tissue elastance. The first of these elements is now available using the methods presented herein. The later two parameters will be required for further robust modeling and analysis. FEM analysis has not been feasible for mitral valve in the intraoperative setting. The major roadblock was the complexity and the practicality of transfer of structural 3D data to a FEM analysis program. This study describes a method to rapidly transfer 3D structural data from the TEE device into a FEM model which can be loaded easily by standard FEM analysis programs. Once measured pressure and calculated elastance are added to the model, near real-time dynamic stress-strain information in the operating room will be achievable.
- Delabays A, Jeanrenaud X, Chassot PG, Von Segesser LK, Kappenberger L: Localization and quantification of mitral valve prolapse using three-dimensional echocardiography. European Journal of Echocardiography 2004, 5: 422–429. 10.1016/j.euje.2004.03.007View ArticleGoogle Scholar
- Lim KH, Yeo JH, Duran CM: Three-dimensional asymmetrical modeling of the mitral valve: a finite element study with dynamic boundaries. J Heart Valve Dis 2005, 14: 386–392.Google Scholar
- Liu F, Lu WX, Xia L, Wu GH: The construction of three-dimensional composite finite element mechanical model of human left ventricle. JSME International Journal Series C Mechanical Systems Machine Elements & Manufacturing 2001, 44: 125–133.View ArticleGoogle Scholar
- Verhey JF, Nathan NS: Feasibility of rapid and automated importation of 3D echocardiographic left ventricular (LV) geometry into a finite element (FEM) analysis model. BioMedical Engineering OnLine 2004, 32: 1–15.Google Scholar
- Usyk TP, Belik ME, Michailova A, McCulloch AD: Three-dimensional model of cardiac electromechanics: cell to organ. Conference Proceedings Second Joint EMBS BMES Conference 2002, 1220–1221.Google Scholar
- Mackerle J: Finite element modelling and simulations in cardiovascular mechanics and cardiology: A bibliography 1993–2004. Computer Methods in Biomechanics & Biomedical Engineering 2005, 8: 59 -581. 10.1080/10255840500141486View ArticleGoogle Scholar
- Crowley JJ, Huang CL, Gates AR, Basu A, Shapiro LM, Carpenter TA, Hall LD: A quantitative description of dynamic left ventricular geometry in anaesthetized rats using magnetic resonance imaging. Exp Physiol 1997, 82: 887–904.View ArticleGoogle Scholar
- Pham QC, Vincent F, Clarysse P, Croisille P, Magnin IE: A FEM-based deformable model for the 3D segmentation and tracking of the heart in cardiac MRI. Ispa 2001, 250–254.Google Scholar
- Lembcke A, Wiese TH, Enzweiler CN, Kivelitz DE, Dushe S, Dohmen PM, Borges AC, Rogalla P, Hamm B: Quantification of mitral valve regurgitation by left ventricular volume and flow measurements using electron beam computed tomography: comparison with magnetic resonance imaging. J Comput Assist Tomogr 2003, 27: 385–391. 10.1097/00004728-200305000-00015View ArticleGoogle Scholar
- McCulloch A, Waldman L, Rogers J, Guccione J: Large-scale finite element analysis of the beating heart. Critical Reviews in Biomedical Engineering 1992, 20: 427–449.Google Scholar
- TomTec Imaging Systems GH: 4D LV-Analysis TEE. Unterschleissheim, Germany, ; 2002. [http://www.tomtec.de]Google Scholar
- Nelson TR, Downey DB, Pretorius DH, Fenster A: Three-Dimensional Ultrasound. , Lippincott Williams & Wilkins, Philadelphia, New York, Baltimore; 1999.Google Scholar
- Gering DT, Nabavi A, Kikinis R, Hata N, O'Donnell LJ, Grimson WEL, Jolesz FA, Black PM, Wells WM: An integrated visualization system for surgical planning and guidance using image fusion and an open MR. Journal of Magnetic Resonance Imaging 2001, 13: 967–975. 10.1002/jmri.1139View ArticleGoogle Scholar
- Macromedia: Fireworks 4. 2000.Google Scholar
- Ghista DN, Kobayashi AS, Davis N, Ray G: Finite element analysis in biomechanics. International Conference on Variational Methods in Engineering Southampton Univ Press 1973, 5: 50–81.Google Scholar
- Aamo OM, Fossen TI: Finite element modelling of moored vessels. Mathematical & Computer Modelling of Dynamical Systems 2001, 7: 47–75. 10.1076/mcmd.188.8.131.5232View ArticleGoogle Scholar
- Veress AI, Vince DG, Anderson PM, Cornhill JF, Herderick EE, Kuban BD, Greenberg NL, Thomas JD: Vascular mechanics of the coronary artery. Zeitschrift fur Kardiologie 2000, 89: 92–100. 10.1007/s003920070106View ArticleGoogle Scholar
- Beck A, Thubrikar MJ, Robicsek F: Stress analysis of the aortic valve with and without the sinuses of valsalva. Journal of Heart Valve Disease 2001, 10: 1–11.Google Scholar
- Young AA: Model tags: direct three-dimensional tracking of heart wall motion from tagged magnetic resonance images. Medical Image Analysis 1999, 3: 361–372. 10.1016/S1361-8415(99)80029-2View ArticleGoogle Scholar
- Klepfer RN, Johnson CR, Macleod RS: The effects of inhomogeneities and anisotropies on electrocardiographic fields: a 3-D finite-element study. IEEE Transactions on Biomedical Engineering 1997, 44: 706–719. 10.1109/10.605427View ArticleGoogle Scholar
- Coquart L, Depeursinge C, Curnier A, Ohayon R: A fluid-structure interaction problem in biomechanics: prestressed vibrations of the eye by the finite element method. Journal of Biomechanics 1992, 25: 1105–1118. 10.1016/0021-9290(92)90067-BView ArticleGoogle Scholar
- Girod S, Teschner M, Schrell U, Kevekordes B, Girod B: Computer-aided 3-D simulation and prediction of craniofacial surgery: a new approach. Journal of Cranio Maxillo Facial Surgery 2001, 29: 156–158. 10.1054/jcms.2000.0203View ArticleGoogle Scholar
- Higa M, Nishimura I, Tanino H, Itoh H, Matsuno T, Mitamura Y: Shape optimization of artificial hip prosthesis with 3D FEM. Journal of the Japan Society of Precision Engineering/Seimitsu Kogaku Kaishi 2002, 68: 948–952.View ArticleGoogle Scholar
- Migliavacca F, Pennati G, Di Martino E, Dubini G, Pietrabissa R: Pressure drops in a distensible model of end-to-side anastomosis in systemic-to-pulmonary shunts. Computer Methods in Biomechanics & Biomedical Engineering 2002, 5: 243–248. 10.1080/10255840290010689View ArticleGoogle Scholar
- Feng B, Veress AI, Sitek A, Gullberg GT, Roy DG: Estimation of mechanical properties from gated SPECT and cine MRI data using a finite-element mechanical model of the left ventricle. IEEE Transactions on Nuclear Science 2001, 48: 725–733. 10.1109/23.940154View ArticleGoogle Scholar
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 cited.