Source reconstruction of brain electromagnetic fields — Source iteration of minimum norm (SIMN)

A recursive scheme aiming at obtaining sparse and focal brain electromagnetic source distribution is proposed based on the interpretation that the weighted minimum norm is the minimum norm estimates of amplitudes on grid points for the source distribution specified by the diagonal elements of the weight matrix. The source distribution is updated so that, at each grid point, the number of current dipoles equals the total source strength estimate of the pre-specified current dipoles. The source strength of a pre-specified current dipole is estimated by projecting the vector of minimum norm estimate to the space spanned by the three column vectors, corresponding to the three amplitudes of the current dipole, of the resolution matrix. The norm of the projected vector yields the source strength estimate of the current dipole. Exact inverse solutions are obtained by this source iteration of minimum norm (SIMN) algorithm for noiseless MEG signals from multi-point sources provided the sources are sufficiently sparse and there are no substantial cancellations among the signals of the sources. For noisy data, a set of "noise sources" is introduced. The diagonal matrix formed by the "noise source numbers" plays the role of regularization matrix and Tikhonov regularization is applied to initialize the "noise source numbers". Application to the source localization of real EEG data is also presented.

[1]  R. Leahy,et al.  Equivalence of linear approaches in bioelectromagnetic inverse solutions , 2004, IEEE Workshop on Statistical Signal Processing, 2003.

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

[3]  Fusheng Yang,et al.  Standardized shrinking LORETA-FOCUSS (SSLOFO): a new algorithm for spatio-temporal EEG source reconstruction , 2005, IEEE Transactions on Biomedical Engineering.

[4]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[5]  M. Fuchs,et al.  Linear and nonlinear current density reconstructions. , 1999, Journal of clinical neurophysiology : official publication of the American Electroencephalographic Society.

[6]  C. Michel,et al.  Linear inverse solutions with optimal resolution kernels applied to electromagnetic tomography , 1997, Human brain mapping.

[7]  K. Matsuura,et al.  Selective minimum-norm solution of the biomagnetic inverse problem , 1995, IEEE Transactions on Biomedical Engineering.

[8]  C. Aine.  A conceptual overview and critique of functional neuroimaging techniques in humans: I. MRI/FMRI and PET. , 1995, Critical reviews in neurobiology.

[9]  R. Ilmoniemi,et al.  Interpreting magnetic fields of the brain: minimum norm estimates , 2006, Medical and Biological Engineering and Computing.

[10]  S. Gonzalez-Andino,et al.  A critical analysis of linear inverse solutions to the neuroelectromagnetic inverse problem , 1998, IEEE Transactions on Biomedical Engineering.

[11]  S. Thorpe,et al.  A Limit to the Speed of Processing in Ultra-Rapid Visual Categorization of Novel Natural Scenes , 2001, Journal of Cognitive Neuroscience.

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

[13]  Terrence J. Sejnowski,et al.  From single-trial EEG to brain area dynamics , 2002, Neurocomputing.

[14]  Aine Cj,et al.  A conceptual overview and critique of functional neuroimaging techniques in humans: I. MRI/FMRI and PET. , 1995 .

[15]  H O'Beirne,et al.  A light source for pupil stimulation. , 1969, IEEE transactions on bio-medical engineering.

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

[17]  Gene H. Golub,et al.  Matrix computations , 1983 .

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

[19]  I F Gorodnitsky,et al.  Neuromagnetic source imaging with FOCUSS: a recursive weighted minimum norm algorithm. , 1995, Electroencephalography and clinical neurophysiology.

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

[21]  Richard M. Leahy,et al.  A nonlinear simplex algorithm for minimum order solutions , 1988, ICASSP-88., International Conference on Acoustics, Speech, and Signal Processing.

[22]  Richard M. Leahy,et al.  MEG-based imaging of focal neuronal current sources , 1995, 1995 IEEE Nuclear Science Symposium and Medical Imaging Conference Record.

[23]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[24]  A. Ioannides,et al.  Continuous probabilistic solutions to the biomagnetic inverse problem , 1990 .

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