Raviart–Thomas-type sources adapted to applied EEG and MEG: implementation and results

This paper studies numerically electroencephalography and magnetoencephalography (EEG and MEG), two non-invasive imaging modalities in which external measurements of the electric potential and the magnetic field are, respectively, utilized to reconstruct the primary current density (neuronal activity) of the human brain. The focus is on adapting a Raviart–Thomas-type source model to meet the needs of EEG and MEG applications. The goal is to construct a model that provides an accurate approximation of dipole source currents and can be flexibly applied to different reconstruction strategies as well as to realistic computation geometries. The finite element method is applied in the simulation of the data. Least-squares fit interpolation is used to establish Cartesian source directions, which guarantee that the recovered current field is minimally dependent on the underlying finite element mesh. Implementation is explained in detail and made accessible, e.g., by using quadrature-free formulae and the Gaussian one-point rule in numerical integration. Numerical results are presented concerning, for example, the iterative alternating sequential inverse algorithm as well as resolution, smoothness and local refinement of the finite element mesh. Both spherical and pseudo-realistic head models, as well as real MEG data, are utilized in the numerical experiments.

[1]  Harri Hakula,et al.  Conditionally Gaussian Hypermodels for Cerebral Source Localization , 2008, SIAM J. Imaging Sci..

[2]  M. Bertero,et al.  Projected Landweber method and preconditioning , 1997 .

[3]  Katrina Wendel,et al.  The Influence of CSF on EEG Sensitivity Distributions of Multilayered Head Models , 2008, IEEE Transactions on Biomedical Engineering.

[4]  Nevzat G. Gencer,et al.  Parallel implementation of the accelerated BEM approach for EMSI of the human brain , 2008, Medical & Biological Engineering & Computing.

[5]  J.B. Allen,et al.  A unified approach to short-time Fourier analysis and synthesis , 1977, Proceedings of the IEEE.

[6]  Y. D'Asseler,et al.  A finite difference method with reciprocity used to incorporate anisotropy in electroencephalogram dipole source localization , 2005 .

[7]  J.C. Mosher,et al.  Multiple dipole modeling and localization from spatio-temporal MEG data , 1992, IEEE Transactions on Biomedical Engineering.

[8]  R. Ilmoniemi,et al.  Magnetoencephalography-theory, instrumentation, and applications to noninvasive studies of the working human brain , 1993 .

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

[10]  M. Scherg,et al.  Two bilateral sources of the late AEP as identified by a spatio-temporal dipole model. , 1985, Electroencephalography and clinical neurophysiology.

[11]  Yutaka Ohtake,et al.  Mesh regularization and adaptive smoothing , 2001, Comput. Aided Des..

[12]  E. Somersalo,et al.  Visualization of Magnetoencephalographic Data Using Minimum Current Estimates , 1999, NeuroImage.

[13]  Peter C. Hansen,et al.  MEG. An introduction to methods , 2010 .

[14]  Sampsa Pursiainen.  EEG/MEG forward simulation through h- and p-type finite elements , 2008 .

[15]  A. Schnitzler,et al.  The neural basis of intermittent motor control in humans , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Théodore Papadopoulo,et al.  A Trilinear Immersed Finite Element Method for Solving the Electroencephalography Forward Problem , 2010, SIAM J. Sci. Comput..

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

[18]  H. Spekreijse,et al.  Mathematical dipoles are adequate to describe realistic generators of human brain activity , 1988, IEEE Transactions on Biomedical Engineering.

[19]  Sampsa Pursiainen.  Coarse-to-fine reconstruction in linear inverse problems with application to limited-angle computerized tomography , 2008 .

[20]  Sampsa Pursiainen,et al.  Computational methods in electromagnetic biomedical inverse problems , 2008 .

[21]  J. P. Ary,et al.  Location of Sources of Evoked Scalp Potentials: Corrections for Skull and Scalp Thicknesses , 1981, IEEE Transactions on Biomedical Engineering.

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

[23]  Jukka Nenonen,et al.  Representation of bioelectric current sources using Whitney elements in the finite element method , 2005, Physics in medicine and biology.

[24]  Gernot Beer,et al.  The Boundary Element Method with Programming: For Engineers and Scientists , 2008 .

[25]  Thom F. Oostendorp,et al.  The conductivity of the human skull: results of in vivo and in vitro measurements , 2000, IEEE Transactions on Biomedical Engineering.

[26]  R Van Uitert,et al.  Finite Elenent EEG and MEG Simulations for Realistic Head Models: Quadratic vs. Linear Approximations , 2001 .

[27]  Chris R. Johnson,et al.  Volume Currents in Forward and Inverse Magnetoencephalographic Simulations Using Realistic Head Models , 2004, Annals of Biomedical Engineering.

[28]  J. J. Ermer,et al.  Rapidly recomputable EEG forward models for realistic head shapes. , 2001, Physics in medicine and biology.

[29]  Bart Vanrumste,et al.  Dipole estimation errors due to differences in modeling anisotropic conductivities in realistic head models for EEG source analysis , 2008, Physics in medicine and biology.

[30]  Olivier D. Faugeras,et al.  Symmetric BEM Formulation for the M/EEG Forward Problem , 2003, IPMI.

[31]  Simon Wotherspoon.  Numerical Integration 2 , 1994 .

[32]  Lauri Parkkonen,et al.  Dynamical MEG source modeling with multi‐target Bayesian filtering , 2009, Human brain mapping.

[33]  Olivier D. Faugeras,et al.  A common formalism for the Integral formulations of the forward EEG problem , 2005, IEEE Transactions on Medical Imaging.

[34]  Jaakko Malmivuo,et al.  Effect of skull resistivity on the spatial resolutions of EEG and MEG , 2004, IEEE Transactions on Biomedical Engineering.

[35]  Sampsa Pursiainen,et al.  Forward simulation and inverse dipole localization with the lowest order Raviart—Thomas elements for electroencephalography , 2011 .

[36]  Isaac Amidror,et al.  Scattered data interpolation methods for electronic imaging systems: a survey , 2002, J. Electronic Imaging.

[37]  J. Sarvas.  Basic mathematical and electromagnetic concepts of the biomagnetic inverse problem. , 1987, Physics in medicine and biology.

[38]  Z. Zhang,et al.  A fast method to compute surface potentials generated by dipoles within multilayer anisotropic spheres. , 1995, Physics in medicine and biology.