Multiple dipole modeling and localization from spatio-temporal MEG data
- 1 June 1992
- journal article
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Biomedical Engineering
- Vol. 39 (6) , 541-557
- https://doi.org/10.1109/10.141192
Abstract
An array of biomagnetometers may be used to measure the spatio-temporal neuromagnetic field or magnetoencephalogram (MEG) produced by neural activity in the brain. A popular model for the neural activity produced in response to a given sensory stimulus is a set of current dipoles, where each dipole represents the primary current associated with the combined activation of a large number of neurons located in a small volume of the brain. An important problem in the interpretation of MEG data from evoked response experiments is the localization of these neural current dipoles. We present here a linear algebraic framework for three common spatio-temporal dipole models: i) unconstrained dipoles, ii) dipoles with a fixed location, and iii) dipoles with a fixed orientation and location. In all cases, we assume that the location, orientation, and magnitude of the dipoles are unknown. With a common model, we show how the parameter estimation problem may be decomposed into the estimation of the time invariant parameters using nonlinear least-squares minimization, followed by linear estimation of the associated time varying parameters. A subspace formulation is presented and used to derive a suboptimal least-squares subspace scanning method. The resulting algorithm is a special case of the well-known MUltiple SIgnal Classification (MUSIC) method, in which the solution (multiple dipole locations) is found by scanning potential locations using a simple one dipole model. Principal components analysis (PCA) dipole fitting has also been used to individually fit single dipoles in a multiple dipole problem. Analysis is presented here to show why PCA dipole fitting will fail in general, whereas the subspace method presented here will generally succeed. Numerically efficient means of calculating the cost functions are presented, and problems of model order selection and missing moments are discussed. Results from a simulation and a somatosensory experiment are presented.Keywords
This publication has 30 references indexed in Scilit:
- MUSIC, maximum likelihood and Cramer-Rao bound: further results and comparisonsPublished by Institute of Electrical and Electronics Engineers (IEEE) ,2003
- Sensor array processing based on subspace fittingIEEE Transactions on Signal Processing, 1991
- Detection of the number of signals: a predicted eigen-threshold approachIEEE Transactions on Signal Processing, 1991
- Instrumentation and techniques for high-resolution magnetic imagingPublished by SPIE-Intl Soc Optical Eng ,1990
- Maximum likelihood methods for direction-of-arrival estimationIEEE Transactions on Acoustics, Speech, and Signal Processing, 1990
- Multiple source location-the signal subspace approachIEEE Transactions on Acoustics, Speech, and Signal Processing, 1990
- Principal component analysis of event-related potentials: A note on misallocation of varianceElectroencephalography and Clinical Neurophysiology/Evoked Potentials Section, 1986
- A new interpretation of the generators of BAEP waves I–V: Results of a spatio-temporal dipole modelElectroencephalography and Clinical Neurophysiology/Evoked Potentials Section, 1985
- Principal component analysis of event-related potentials: Simulation studies demonstrate misallocation of variance across componentsElectroencephalography and Clinical Neurophysiology/Evoked Potentials Section, 1984
- APPLICATION OF DIPOLE LOCALIZATION METHODS TO SOURCE IDENTIFICATION OF HUMAN EVOKED POTENTIALS*Annals of the New York Academy of Sciences, 1982