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.