Sensitivity of the Projected Subtraction Approach to Mesh Degeneracies and Its Impact on the Forward Problem in EEG

Objective: Subtraction-based techniques are known for being theoretically rigorous and accurate methods for solving the forward problem in electroencephalography (EEG-FP) by means of the finite-element method. Within them, the projected subtraction (PS) approach is generally adopted because of its computational efficiency. Although this technique received the attention of the community, its sensitivity to degenerated elements is still poorly understood. In this paper, we investigate the impact of low-quality tetrahedra on the results computed with the PS approach. Methods: We derived upper bounds on the relative error of the element source vector as a function of geometrical features describing the tetrahedral discretization of the domain. These error bounds were then utilized for showing the instability of the PS method with regards to the mesh quality. To overcome this issue, we proposed an alternative technique, coined projected gradient subtraction (PGS) approach, that exploits the stability of the corresponding bounds. Results: Computer simulations showed that the PS method is extremely sensitive to the mesh shape and size, leading to unacceptable solutions of the EEG-FP in case of using suboptimal tessellations. This was not the case of the PGS approach, which led to stable and accurate results in a comparable amount of time. Conclusion: Solutions of the EEG-FP computed with the PS method are highly sensitive to degenerated elements. Such errors can be mitigated by the PGS approach, which showed better performance than the PS technique. Significance: The PGS is an efficient method for computing high-quality lead field matrices even in the presence of degenerated elements.

[1]  S. Kochen,et al.  Approximate average head models for EEG source imaging , 2009, Journal of Neuroscience Methods.

[2]  C H Wolters,et al.  Accuracy and run-time comparison for different potential approaches and iterative solvers in finite element method based EEG source analysis. , 2009, Applied numerical mathematics : transactions of IMACS.

[3]  Jean Gotman,et al.  Extent of cortical generators visible on the scalp: Effect of a subdural grid , 2014, NeuroImage.

[4]  Leandro Beltrachini,et al.  Neural Processing of Emotional Facial and Semantic Expressions in Euthymic Bipolar Disorder (BD) and Its Association with Theory of Mind (ToM) , 2012, PloS one.

[5]  R D Pascual-Marqui,et al.  Standardized low-resolution brain electromagnetic tomography (sLORETA): technical details. , 2002, Methods and findings in experimental and clinical pharmacology.

[6]  Per-Olof Persson,et al.  A Simple Mesh Generator in MATLAB , 2004, SIAM Rev..

[7]  Moritz Dannhauer,et al.  Modeling of the human skull in EEG source analysis , 2011, Human brain mapping.

[8]  M. Sigman,et al.  Cortical deficits of emotional face processing in adults with ADHD: Its relation to social cognition and executive function , 2011, Social neuroscience.

[9]  J. Haueisen,et al.  Influence of head models on EEG simulations and inverse source localizations , 2006, Biomedical engineering online.

[10]  J. D. Munck,et al.  A fast method to compute the potential in the multisphere model (EEG application) , 1993, IEEE Transactions on Biomedical Engineering.

[11]  Carsten H. Wolters,et al.  A full subtraction approach for finite element method based source analysis using constrained Delaunay tetrahedralisation , 2009, NeuroImage.

[12]  C. H. Woltersa,et al.  Numerical approaches for dipole modeling in finite element method based source analysis , 2007 .

[13]  H. Si,et al.  Adaptive tetrahedral mesh generation by constrained Delaunay refinement , 2008 .

[14]  M. Clerc,et al.  Comparison of Boundary Element and Finite Element Approaches to the EEG Forward Problem , 2012, Biomedizinische Technik. Biomedical engineering.

[15]  L. Formaggia,et al.  Degenerate tetrahedra removal , 2016 .

[16]  Cheng Guan Koay,et al.  Analytically exact spiral scheme for generating uniformly distributed points on the unit sphere , 2011, J. Comput. Sci..

[17]  D. Louis Collins,et al.  Unbiased average age-appropriate atlases for pediatric studies , 2011, NeuroImage.

[18]  Thomas R. Knösche,et al.  A guideline for head volume conductor modeling in EEG and MEG , 2014, NeuroImage.

[19]  Joachim Schöberl,et al.  NETGEN An advancing front 2D/3D-mesh generator based on abstract rules , 1997 .

[20]  Jonathan Richard Shewchuk,et al.  Aggressive Tetrahedral Mesh Improvement , 2007, IMR.

[21]  David A. Boas,et al.  Tetrahedral mesh generation from volumetric binary and grayscale images , 2009, 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[22]  J. Shewchuk What Is a Good Linear Finite Element? Interpolation, Conditioning, Anisotropy, and Quality Measures , 2002 .

[23]  Scott Makeig,et al.  Simultaneous head tissue conductivity and EEG source location estimation , 2016, NeuroImage.

[24]  Carlos H. Muravchik,et al.  On the EEG/MEG forward problem solution for distributed cortical sources , 2009, Medical & Biological Engineering & Computing.

[25]  Andreas Schulze-Bonhage,et al.  The role of blood vessels in high-resolution volume conductor head modeling of EEG , 2015, NeuroImage.

[26]  T. Preston,et al.  Finite Elements for Electrical Engineers , 1984 .

[27]  Carlos H. Muravchik,et al.  General bounds for electrode mislocation on the EEG inverse problem , 2011, Comput. Methods Programs Biomed..

[28]  David R. Wozny,et al.  The electrical conductivity of human cerebrospinal fluid at body temperature , 1997, IEEE Transactions on Biomedical Engineering.

[29]  Carlos H. Muravchik,et al.  Analysis of parametric estimation of head tissue conductivities using Electrical Impedance Tomography , 2013, Biomed. Signal Process. Control..