A new finite element formulation for the forward problem of electro-magnetic source imaging

The voltage measurements on the scalp and magnetic field measurements near the scalp can be used to compute spatio-temporal evolution of electric currents across the cerebral cortex. Such a representation can be called Electro-Magnetic Source Image. In order to interpret the measured fields accurately, a realistic numerical model of the head should be prepared that includes physical properties and geometry of the individual's head. For that purpose, in this study, Finite Element Method (FEM) is used to solve the scalar potential function. A new FEM formulation is derived that assumes constant current density in isoparametric volume elements. A software is prepared by optimizing the code to run in a PC with 128 MB memory. It is shown that, for a current dipole at 2.5 cm from the origin, the potential distribution can be solved within an accuracy of 0.5% in comparison with the exact solutions, with a mesh of 29001 nodes and 7000 elements. In order to solve the magnetic held, the surface integration tools developed in the same FEM code are used. For the same dipole location the radial magnetic fields are calculated within an accuracy of 2.2%. Obviously, the accuracy in results can be increased by increasing the number of nodes.

[1]  A. Dale,et al.  Improved Localizadon of Cortical Activity by Combining EEG and MEG with MRI Cortical Surface Reconstruction: A Linear Approach , 1993, Journal of Cognitive Neuroscience.

[2]  N. G. Gencer,et al.  Differential characterization of neural sources with the bimodal truncated SVD pseudo-inverse for EEG and MEG measurements , 1998, IEEE Transactions on Biomedical Engineering.

[3]  J. Le,et al.  Method to reduce blur distortion from EEG's using a realistic head model , 1993, IEEE Transactions on Biomedical Engineering.

[4]  John S. George,et al.  Three-Dimensional Volumetric Reconstruction for Neuromagnetic Source Localization , 1989 .

[5]  David R. Owen,et al.  FINITE ELEMENT PROGRAMMING , 1980, The Finite Element Method Using MATLAB.

[6]  D. A. Driscoll,et al.  EEG electrode sensitivity--an application of reciprocity. , 1969, IEEE transactions on bio-medical engineering.

[7]  H Hosaka,et al.  Part II: magnetic field produced by a current dipole. , 1976, Journal of electrocardiology.

[8]  Toshimitsu Musha,et al.  Electric Dipole Tracing in the Brain by Means of the Boundary Element Method and Its Accuracy , 1987, IEEE Transactions on Biomedical Engineering.

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

[10]  M. Hämäläinen,et al.  Realistic conductivity geometry model of the human head for interpretation of neuromagnetic data , 1989, IEEE Transactions on Biomedical Engineering.

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

[12]  D. Geselowitz On bioelectric potentials in an inhomogeneous volume conductor. , 1967, Biophysical journal.

[13]  F. Perrin,et al.  The finite element method for a realistic head model of electrical brain activities: preliminary results. , 1991, Clinical physics and physiological measurement : an official journal of the Hospital Physicists' Association, Deutsche Gesellschaft fur Medizinische Physik and the European Federation of Organisations for Medical Physics.