EEG/MEG forward simulation through h- and p-type finite elements

Electro/Magnetoencephalography (EEG/MEG) is a non-invasive imaging modality, in which a primary current density generated by the neural activity in the brain is to be reconstructed from external electric potential/magnetic field measurements. This work focuses on effective and accurate simulation of the EEG/MEG forward model through the h- and p-versions of the finite element method (h- and p-FEM). The goal is to compare the effectiveness of these two versions in forward simulation. Both h- and p-type forward simulations are described and implemented, and the technical solutions found are discussed. These include, for example, suitable ways to generate a finite element mesh for a real head geometry through the use of different element types. Performances of the two implemented forward simulation types are compared by measuring directly the forward modeling error, as well as by computing reconstructions through a regularized FOCUSS (FOCal Underdetermined System Solver) algorithm. The results obtained suggest that the p-type performs better in terms of the forward modeling error. However, both types perform well in regularized FOCUSS reconstruction.

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

[2]  I. Doležel,et al.  Higher-Order Finite Element Methods , 2003 .

[3]  Bhaskar D. Rao,et al.  Subset selection in noise based on diversity measure minimization , 2003, IEEE Trans. Signal Process..

[4]  Jouko Lampinen,et al.  Bayesian analysis of the neuromagnetic inverse problem with ℓ p -norm priors , 2005, NeuroImage.

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

[6]  J P Kaipio,et al.  Effects of electrode properties on EEG measurements and a related inverse problem. , 2000, Medical engineering & physics.

[7]  E. Somersalo,et al.  Existence and uniqueness for electrode models for electric current computed tomography , 1992 .

[8]  R. Leahy,et al.  On MEG forward modelling using multipolar expansions. , 2002, Physics in medicine and biology.

[9]  Mark Ainsworth,et al.  Hierarchic finite element bases on unstructured tetrahedral meshes , 2003 .

[10]  E. Somersalo,et al.  Statistical and computational inverse problems , 2004 .

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

[12]  Karl J. Friston,et al.  Systematic Regularization of Linear Inverse Solutions of the EEG Source Localization Problem , 2002, NeuroImage.

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

[14]  R. Ilmoniemi,et al.  EEG minimum-norm estimation compared with MEG dipole fitting in the localization of somatosensory sources at S1 , 2004, Clinical Neurophysiology.

[15]  Carlos H. Muravchik,et al.  MEG forward problem formulation using equivalent surface current densities , 2005, IEEE Transactions on Biomedical Engineering.

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