Cardiac anisotropy in boundary-element models for the electrocardiogram

The boundary-element method (BEM) is widely used for electrocardiogram (ECG) simulation. Its major disadvantage is its perceived inability to deal with the anisotropic electric conductivity of the myocardial interstitium, which led researchers to represent only intracellular anisotropy or neglect anisotropy altogether. We computed ECGs with a BEM model based on dipole sources that accounted for a “compound” anisotropy ratio. The ECGs were compared with those computed by a finite-difference model, in which intracellular and interstitial anisotropy could be represented without compromise. For a given set of conductivities, we always found a compound anisotropy value that led to acceptable differences between BEM and finite-difference results. In contrast, a fully isotropic model produced unacceptably large differences. A model that accounted only for intracellular anisotropy showed intermediate performance. We conclude that using a compound anisotropy ratio allows BEM-based ECG models to more accurately represent both anisotropies.

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

[2]  A. Tveito,et al.  Modeling the electrical activity of the heart: A Bidomain Model of the ventricles embedded in a torso , 2002 .

[3]  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.

[4]  Martin Buist,et al.  Torso Coupling Techniques for the Forward Problem of Electrocardiography , 2002, Annals of Biomedical Engineering.

[5]  J. Nenonen,et al.  Simulated Epicardial Potential Maps During Paced Activation Reflect Myocardial Fibrous Structure , 1998, Annals of Biomedical Engineering.

[6]  Jens Haueisen,et al.  Boundary Element Computations in the Forward and Inverse Problems of Electrocardiography: Comparison of Collocation and Galerkin Weightings , 2008, IEEE Transactions on Biomedical Engineering.

[7]  R. Macleod,et al.  The Forward Problem of Electrocardiography , 2010 .

[8]  G. Huiskamp,et al.  Simulation of depolarization in a membrane-equations-based model of the anisotropic ventricle , 1998, IEEE Transactions on Biomedical Engineering.

[9]  R. Gulrajani,et al.  Effects of rotational myocardial anisotropy in forward potential computations with equivalent heart dipoles , 1997, Annals of Biomedical Engineering.

[10]  M Tysler,et al.  Noninvasive assessment of local myocardium repolarization changes using high resolution surface ECG mapping. , 2007, Physiological research.

[11]  J Malmivuo,et al.  Computer model analysis of the relationship of ST-segment and ST-segment/heart rate slope response to the constituents of the ischemic injury source. , 1997, Journal of electrocardiology.

[12]  A. van Oosterom,et al.  ECGSIM: an interactive tool for studying the genesis of QRST waveforms , 2004, Heart.

[13]  A. J. Pullan,et al.  The Forward Problem of Electrocardiography , 2010 .

[14]  B. Taccardi,et al.  Potential Fields Generated by Oblique Dipole Layers Modeling Excitation Wavefronts in the Anisotropic Myocardium: Comparison with Potential Fields Elicited by Paced Dog Hearts in a Volume Conductor , 1982, Circulation research.

[15]  David B. Geselowitz,et al.  On the theory of the electrocardiogram , 1989, Proc. IEEE.

[16]  Thom F. Oostendorp,et al.  Application of the fastest route algorithm in the interactive simulation of the effect of local ischemia on the ECG , 2008, Medical & Biological Engineering & Computing.

[17]  R. Gulrajani,et al.  A Simulation Study of the Effects of Torso Inhomogeneities on Electrocardiographic Potentials, using Realistic Heart and Torso Models , 1983, Circulation research.

[18]  S. Rush,et al.  Qualitative effects of thoracic resistivity variations on the interpretation of electrocardiograms: the low resistance surface layer. , 1968, American Heart Journal.

[19]  M. N. Morrow,et al.  An assessment of variable thickness and fiber orientation of the skeletal muscle layer on electrocardiographic calculations , 1991, IEEE Transactions on Biomedical Engineering.

[20]  R M Gulrajani,et al.  A computer heart model incorporating anisotropic propagation. I. Model construction and simulation of normal activation. , 1993, Journal of electrocardiology.

[21]  O. Dossel,et al.  Computer based modeling of the congenital long-QT 2 syndrome in the Visible Man torso: From genes to ECG , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[22]  T. Musha,et al.  Three-dimensional simulation of the ventricular depolarization and repolarization processes and body surface potentials: normal heart and bundle branch block. , 1987, IEEE transactions on bio-medical engineering.

[23]  R M Gulrajani,et al.  Models of the electrical activity of the heart and computer simulation of the electrocardiogram. , 1988, Critical reviews in biomedical engineering.

[24]  Andrew J. Pullan,et al.  A high order coupled finite element/boundary element torso model , 1993, Proceedings of Computers in Cardiology Conference.

[25]  Gerald M. Saidel,et al.  Role of O2 in Regulation of Lactate Dynamics during Hypoxia: Mathematical Model and Analysis , 2004, Annals of Biomedical Engineering.

[26]  Daming Wei,et al.  Comparative simulation of excitation and body surface electrocardiogram with isotropic and anisotropic computer heart models , 1995, IEEE Transactions on Biomedical Engineering.

[27]  A. van Oosterom,et al.  Non-invasive imaging of activation times during drug inducedconduction changes , 2007 .

[28]  Ling Xia,et al.  Truncated Total Least Squares: A New Regularization Method for the Solution of ECG Inverse Problems , 2008, IEEE Transactions on Biomedical Engineering.

[29]  R. Barr,et al.  Determining surface potentials from current dipoles, with application to electrocardiography. , 1966, IEEE transactions on bio-medical engineering.

[30]  R. Gulrajani The forward and inverse problems of electrocardiography. , 1998, IEEE engineering in medicine and biology magazine : the quarterly magazine of the Engineering in Medicine & Biology Society.

[31]  B. Horáček,et al.  Digital model for studies in magnetocardiography , 1973 .

[32]  R M Gulrajani,et al.  A computer heart model incorporating anisotropic propagation. III. Simulation of ectopic beats. , 1993, Journal of electrocardiology.

[33]  L. Guerri,et al.  Oblique dipole layer potentials applied to electrocardiology , 1983, Journal of mathematical biology.

[34]  D. Geselowitz,et al.  Simulation Studies of the Electrocardiogram: I. The Normal Heart , 1978, Circulation research.

[35]  Joakim Sundnes,et al.  Simulation of ST segment changes during subendocardial ischemia using a realistic 3-D cardiac geometry , 2005, IEEE Transactions on Biomedical Engineering.

[36]  Ling Xia,et al.  Effect of cardiac motion on body surface electrocardiographic potentials: an MRI-based simulation study. , 2006, Physics in medicine and biology.

[37]  Daming Wei,et al.  Computer simulation of intracardiac potential with whole-heart model , 2007, Int. J. Bioinform. Res. Appl..

[38]  R. Gulrajani Bioelectricity and biomagnetism , 1998 .

[39]  Mary N. Morrow,et al.  A Comparison of Finite Element and Integral Equation Formulations for the Calculation of Electrocardiographic Potentials , 1985, IEEE Transactions on Biomedical Engineering.

[40]  C. Henriquez Simulating the electrical behavior of cardiac tissue using the bidomain model. , 1993, Critical reviews in biomedical engineering.

[41]  S. Rush,et al.  Quantitative comparison of pre-mortem ECG's with those reconstructed from activation data of a revived heart. , 1980, Journal of electrocardiology.

[42]  A. M. Scher,et al.  Effect of Tissue Anisotropy on Extracellular Potential Fields in Canine Myocardium in Situ , 1982, Circulation research.

[43]  D. Westbury,et al.  Electrical activity of the heart. , 1971, British medical journal.

[44]  M. Lynn,et al.  The application of electromagnetic theory to electrocardiology. I. Derivation of the integral equations. , 1967, Biophysical journal.

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

[46]  Yoram Rudy,et al.  Critical aspects of the forward and inverse problems in electrocardiography , 1985 .

[47]  A. Oosterom Genesis of the T wave as based on an equivalent surface source model , 2001 .

[48]  Bernhard Pfeifer,et al.  Cardiac anisotropy: is it negligible regarding noninvasive activation time imaging? , 2006, IEEE Transactions on Biomedical Engineering.

[49]  Mary N. Morrow,et al.  A Comparison of Finite Element and Integral Equation Formulations for the Calculation of Electrocardiographic Potentials-II , 1987, IEEE Transactions on Biomedical Engineering.

[50]  M. Lynn,et al.  The application of electromagnetic theory to electrocardiology. II. Numerical solution of the integral equations. , 1967, Biophysical journal.

[51]  B. Taccardi,et al.  Anisotropic Mechanisms for Multiphasic Unipolar Electrograms: Simulation Studies and Experimental Recordings , 2000, Annals of Biomedical Engineering.

[52]  Jane Sands Robb,et al.  The normal heart , 1942 .

[53]  Guy E. Mailloux,et al.  Theoretical Evaluation of the McFee and Frank Vectorcardiographic Lead Systems Using a Numerical Inhomogeneous Torso Model , 1982, IEEE Transactions on Biomedical Engineering.

[54]  Mark Potse,et al.  A Comparison of Monodomain and Bidomain Reaction-Diffusion Models for Action Potential Propagation in the Human Heart , 2006, IEEE Transactions on Biomedical Engineering.

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

[56]  A. Pullan A high-order coupled finite element/boundary element torso model , 1996, IEEE Transactions on Biomedical Engineering.

[57]  Robert Modre,et al.  Lead field computation for the electrocardiographic inverse problem - finite elements versus boundary elements , 2005, Comput. Methods Programs Biomed..

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

[59]  Marie-Claude Trudel,et al.  Simulation of QRST integral maps with a membrane-based computer heart model employing parallel processing , 2004, IEEE Transactions on Biomedical Engineering.

[60]  Clive A. Croxton Liquid State Physics–A Statistical Mechanical Introduction: Numerical solution of the integral equations , 1974 .

[61]  Vincent Jacquemet,et al.  Genesis of the P wave: atrial signals as generated by the equivalent double layer source model. , 2005, Europace : European pacing, arrhythmias, and cardiac electrophysiology : journal of the working groups on cardiac pacing, arrhythmias, and cardiac cellular electrophysiology of the European Society of Cardiology.

[62]  A. J. Pullan,et al.  Mathematical models and numerical methods for the forward problem in cardiac electrophysiology , 2002 .

[63]  T. Musha,et al.  Three-Dimensional Simulation of the Ventricular Depolarization and Repolarization Processes and Body Surface Potentials: Nornal Heart and Bundle Branch Block , 1987, IEEE Transactions on Biomedical Engineering.