Space–time event sparse penalization for magneto-/electroencephalography

This article presents a new spatio-temporal method for M/EEG source reconstruction based on the assumption that only a small number of events, localized in space and/or time, are responsible for the measured signal. Each space-time event is represented using a basis function expansion which reflects the most relevant (or measurable) features of the signal. This model of neural activity leads naturally to a Bayesian likelihood function which balances the model fit to the data with the complexity of the model, where the complexity is related to the number of included events. A novel Expectation-Maximization algorithm which maximizes the likelihood function is presented. The new method is shown to be effective on several MEG simulations of neurological activity as well as data from a self-paced finger tapping experiment.

[1]  Polina Golland,et al.  A distributed spatio-temporal EEG/MEG inverse solver , 2009, NeuroImage.

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

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

[4]  Manbir Singh,et al.  An Evaluation of Methods for Neuromagnetic Image Reconstruction , 1987, IEEE Transactions on Biomedical Engineering.

[5]  Gabriel Curio,et al.  Current multipole expansion to estimate lateral extent of neuronal activity: a theoretical analysis , 2000, IEEE Trans. Biomed. Eng..

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

[7]  Tom E. Bishop,et al.  Blind Image Restoration Using a Block-Stationary Signal Model , 2006, 2006 IEEE International Conference on Acoustics Speech and Signal Processing Proceedings.

[8]  A. Brix.  Bayesian Data Analysis, 2nd edn , 2005 .

[9]  J. Tropp.  Algorithms for simultaneous sparse approximation. Part II: Convex relaxation , 2006, Signal Process..

[10]  Dmitry M. Malioutov,et al.  A sparse signal reconstruction perspective for source localization with sensor arrays , 2005, IEEE Transactions on Signal Processing.

[11]  J. Culham,et al.  The role of parietal cortex in visuomotor control: What have we learned from neuroimaging? , 2006, Neuropsychologia.

[12]  David L. Donoho,et al.  De-noising by soft-thresholding , 1995, IEEE Trans. Inf. Theory.

[13]  M. Paradiso,et al.  Neuroscience: Exploring the Brain , 1996 .

[14]  Leonard J. Trejo,et al.  Feature Extraction of Event-Related Potentials Using Wavelets: An Application to Human Performance Monitoring , 1999, Brain and Language.

[15]  John K Kruschke,et al.  Bayesian data analysis. , 2010, Wiley interdisciplinary reviews. Cognitive science.

[16]  New York Dover,et al.  ON THE CONVERGENCE PROPERTIES OF THE EM ALGORITHM , 1983 .

[17]  Stephen J. Wright,et al.  Simultaneous Variable Selection , 2005, Technometrics.

[18]  Barry D. Van Veen,et al.  Cortical patch basis model for spatially extended neural activity , 2006, IEEE Transactions on Biomedical Engineering.

[19]  Karl J. Friston,et al.  A mesostate-space model for EEG and MEG , 2007, NeuroImage.

[20]  Anders M. Dale,et al.  Vector-based spatial–temporal minimum L1-norm solution for MEG , 2006, NeuroImage.

[21]  D. Lehmann,et al.  Low resolution electromagnetic tomography: a new method for localizing electrical activity in the brain. , 1994, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.

[22]  Emmanuel J. Candès,et al.  Decoding by linear programming , 2005, IEEE Transactions on Information Theory.

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

[24]  Bhaskar D. Rao,et al.  Sparse signal reconstruction from limited data using FOCUSS: a re-weighted minimum norm algorithm , 1997, IEEE Trans. Signal Process..

[25]  Aleksandar Dogandzic,et al.  Estimating evoked dipole responses in unknown spatially correlated noise with EEG/MEG arrays , 2000, IEEE Trans. Signal Process..

[26]  Benoit Cottereau,et al.  Multiresolution imaging of MEG cortical sources using an explicit piecewise model , 2007, NeuroImage.

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

[28]  K. Lehnertz,et al.  Single trial analysis of event related potentials: non-linear de-noising with wavelets , 2000, Clinical Neurophysiology.

[29]  Joe Brewer,et al.  Kronecker products and matrix calculus in system theory , 1978 .

[30]  D. Donoho.  For most large underdetermined systems of linear equations the minimal 𝓁1‐norm solution is also the sparsest solution , 2006 .

[31]  K. Matsuura,et al.  A robust reconstruction of sparse biomagnetic sources , 1997, IEEE Transactions on Biomedical Engineering.

[32]  Stephen J. Wright,et al.  Sparse Reconstruction by Separable Approximation , 2008, IEEE Transactions on Signal Processing.

[33]  Ahmet Ademoglu,et al.  Analysis of event-related potentials (ERP) by damped sinusoids , 1998, Biological Cybernetics.

[34]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[35]  Tapabrata Maiti,et al.  Bayesian Data Analysis (2nd ed.) (Book) , 2004 .

[36]  Bin He,et al.  Sparse Source Imaging in EEG , 2007, 2007 Joint Meeting of the 6th International Symposium on Noninvasive Functional Source Imaging of the Brain and Heart and the International Conference on Functional Biomedical Imaging.

[37]  丸山 徹.  Convex Analysisの二,三の進展について , 1977 .

[38]  Richard M. Leahy,et al.  Electromagnetic brain mapping , 2001, IEEE Signal Process. Mag..

[39]  Polina Golland,et al.  A Distributed Spatio-temporal EEG/MEG Inverse Solver , 2008, MICCAI.

[40]  Joel A. Tropp,et al.  ALGORITHMS FOR SIMULTANEOUS SPARSE APPROXIMATION , 2006 .

[41]  Nelson J. Trujillo-Barreto,et al.  Bayesian M/EEG source reconstruction with spatio-temporal priors , 2008, NeuroImage.

[42]  Diego Clonda,et al.  Bayesian spatio-temporal approach for EEG source reconstruction: conciliating ECD and distributed models , 2006, IEEE Transactions on Biomedical Engineering.

[43]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[44]  Sylvain Baillet,et al.  A Bayesian approach to introducing anatomo-functional priors in the EEG/MEG inverse problem , 1997, IEEE Transactions on Biomedical Engineering.

[45]  Richard M. Leahy,et al.  MEG Source Imaging Using Multipolar Expansions , 1999, IPMI.

[46]  M. Yuan,et al.  Model selection and estimation in regression with grouped variables , 2006 .

[47]  Robert D. Nowak,et al.  An EM algorithm for wavelet-based image restoration , 2003, IEEE Trans. Image Process..

[48]  Robert T. Knight,et al.  Five-dimensional neuroimaging: Localization of the time–frequency dynamics of cortical activity , 2008, NeuroImage.

[49]  W. Drongelen,et al.  Localization of brain electrical activity via linearly constrained minimum variance spatial filtering , 1997, IEEE Transactions on Biomedical Engineering.

[50]  Andrew K. Bolstad,et al.  Magneto-/Electroencephalography with Space-Time Sparse Priors , 2007, 2007 IEEE/SP 14th Workshop on Statistical Signal Processing.