Lead field computation for the electrocardiographic inverse problem - finite elements versus boundary elements

In order to be able to solve the inverse problem of electrocardiography, the lead field matrix (transfer matrix) has to be calculated. The two methods applied for computing this matrix, which are compared in this study, are the boundary element method (BEM) and the finite element method (FEM). The performance of both methods using a spherical model was investigated. For a comparable discretization level, the BEM yields smaller relative errors compared to analytical solutions. The BEM needs less computation time, but a larger amount of memory. Inversely calculated myocardial activation times using either the FEM or BEM computed lead field matrices give similar activation time patterns. The FEM, however, is also capable of considering anisotropic conductivities. This property might have an impact for future development, when also individual myocardial fiber architecture can be considered in the inverse formulation.

[1]  R. Leahy,et al.  EEG and MEG: forward solutions for inverse methods , 1999, IEEE Transactions on Biomedical Engineering.

[2]  Robert Modre,et al.  Noninvasive myocardial activation time imaging: a novel inverse algorithm applied to clinical ECG mapping data , 2002, IEEE Transactions on Biomedical Engineering.

[3]  F. Hanser,et al.  Ventricular surface activation time imaging from electrocardiogram mapping data , 2004, Medical and Biological Engineering and Computing.

[4]  P. Hunter,et al.  Mathematical model of geometry and fibrous structure of the heart. , 1991, The American journal of physiology.

[5]  G. Huiskamp,et al.  A new method for myocardial activation imaging , 1997, IEEE Transactions on Biomedical Engineering.

[6]  C. Henriquez,et al.  A finite volume model of cardiac propagation , 1997, Annals of Biomedical Engineering.

[7]  B. Vanrumste,et al.  Comparing iterative solvers for linear systems associated with the finite difference discretisation of the forward problem in electro-encephalographic source analysis , 2006, Medical and Biological Engineering and Computing.

[8]  Robert Modre,et al.  Model-based imaging of cardiac electrical excitation in humans , 2002, IEEE Transactions on Medical Imaging.

[9]  Robert Modre,et al.  An iterative algorithm for myocardial activation time imaging , 2001, Comput. Methods Programs Biomed..

[10]  C. Henriquez,et al.  Magnetic resonance myocardial fiber-orientation mapping with direct histological correlation. , 1998, American journal of physiology. Heart and circulatory physiology.

[11]  P. Hunter,et al.  A Deformable Finite Element Derived Finite Difference Method for Cardiac Activation Problems , 2003, Annals of Biomedical Engineering.

[12]  F. Hanser,et al.  Clinical ECG mapping and imaging of cardiac electrical excitation. , 2002, Journal of electrocardiology.

[13]  G. Huiskamp,et al.  A Bidomain Model Based BEM-FEM Coupling Formulation for Anisotropic Cardiac Tissue , 2004, Annals of Biomedical Engineering.

[14]  David B. Geselowitz,et al.  A bidomain model for anisotropic cardiac muscle , 2006, Annals of Biomedical Engineering.

[15]  P. Wach,et al.  Stability of Activation Time Imaging From Single Beat Data Under Clinical Conditions , 2001 .

[16]  J Malmivuo,et al.  A software implementation for detailed volume conductor modelling in electrophysiology using finite difference method. , 1999, Computer methods and programs in biomedicine.

[17]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[18]  D. H. Campen,et al.  Optimization of Cardiac Fiber Orientation for Homogeneous Fiber Strain During Ejection , 1999, Annals of Biomedical Engineering.

[19]  Gulrajani Rm,et al.  Models of the electrical activity of the heart and computer simulation of the electrocardiogram. , 1988 .

[20]  C R Johnson,et al.  Computational and numerical methods for bioelectric field problems. , 1997, Critical reviews in biomedical engineering.

[21]  Robert Modre,et al.  On modeling the Wilson terminal in the boundary and finite element method , 2002, IEEE Transactions on Biomedical Engineering.

[22]  H. Nowak,et al.  Application of high-order boundary elements to the electrocardiographic inverse problem. , 1999, Computer methods and programs in biomedicine.

[23]  Jens Haueisen,et al.  Neuromagnetic field computation using the multiple multipole method , 1996 .

[24]  B. Roth.  Electrical conductivity values used with the bidomain model of cardiac tissue , 1997, IEEE Transactions on Biomedical Engineering.

[25]  R. Taylor.  The Finite Element Method, the Basis , 2000 .

[26]  M. Schocke,et al.  Atrial Noninvasive Activation Mapping of Paced Rhythm Data , 2003, Journal of cardiovascular electrophysiology.