A Weighted Approach for Sparse Signal Support Estimation with Application to EEG Source Localization

In sparse signal recovery problems, <inline-formula><tex-math notation="LaTeX">$\ell _1$</tex-math></inline-formula> -norm minimization is typically used as an alternative to more complex <inline-formula><tex-math notation="LaTeX">$\ell _0$</tex-math></inline-formula>-norm minimization. The range space property (RSP) provides the conditions under which the least <inline-formula><tex-math notation="LaTeX">$\ell _1$</tex-math></inline-formula>-norm solution is equal to at most one of the least <inline-formula><tex-math notation="LaTeX">$\ell _0$</tex-math></inline-formula>-norm solutions. These conditions depend on the sensing matrix and the support of the underlying sparse solution. In this paper, we address the problem of recovering sparse signals by weighting the corresponding sensing matrix with a diagonal matrix. We show that by appropriately choosing the weights, we can formulate an <inline-formula> <tex-math notation="LaTeX">$\ell _1$</tex-math></inline-formula>-norm minimization problem that satisfies the RSP, even if the original problem does not. By solving the weighted problem we can obtain the support of the original problem. We provide the conditions which the weights must satisfy, for both noise free and noisy cases. Although the precise conditions involve information about the support of the sparse vector, the class of good weights is very wide, and in most cases encompasses an estimate of the underlying vector obtained via a conventional method, i.e., a method that does not encourage sparsity. The proposed approach is a good candidate for Electroencephalography (EEG) sparse source localization, where the corresponding sensing matrix has high coherence. The performance of the proposed approach is evaluated via simulations and also via experiments on localizing active sources in the brain corresponding to an auditory task from EEG recordings of a human subject.

[1]  EEG signatures of auditory activity correlate with simultaneously recorded fMRI responses in humans , 2010, NeuroImage.

[2]  Kensuke Sekihara,et al.  A novel adaptive beamformer for MEG source reconstruction effective when large background brain activities exist , 2006, IEEE Transactions on Biomedical Engineering.

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

[4]  M. Fuchs,et al.  An improved boundary element method for realistic volume-conductor modeling , 1998, IEEE Transactions on Biomedical Engineering.

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

[6]  O. L. Frost,et al.  An algorithm for linearly constrained adaptive array processing , 1972 .

[7]  Andreas Daffertshofer,et al.  The relationship between structural and functional connectivity: Graph theoretical analysis of an EEG neural mass model , 2010, NeuroImage.

[8]  David Poeppel,et al.  Reconstructing spatio-temporal activities of neural sources using an MEG vector beamformer technique , 2001, IEEE Transactions on Biomedical Engineering.

[9]  Arnaud Delorme,et al.  EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis , 2004, Journal of Neuroscience Methods.

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

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

[12]  Bart Vanrumste,et al.  Review on solving the forward problem in EEG source analysis , 2007, Journal of NeuroEngineering and Rehabilitation.

[13]  Michael Elad,et al.  Optimally sparse representation in general (nonorthogonal) dictionaries via ℓ1 minimization , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[14]  D. Mathalon,et al.  Event-related EEG time-frequency analysis: an overview of measures and an analysis of early gamma band phase locking in schizophrenia. , 2008, Schizophrenia bulletin.

[15]  M Congedo,et al.  A review of classification algorithms for EEG-based brain–computer interfaces , 2007, Journal of neural engineering.

[16]  B. Rao,et al.  Forward sequential algorithms for best basis selection , 1999 .

[17]  Dezhong Yao,et al.  Lp Norm Iterative Sparse Solution for EEG Source Localization , 2007, IEEE Transactions on Biomedical Engineering.

[18]  D. Tucker,et al.  EEG coherency. I: Statistics, reference electrode, volume conduction, Laplacians, cortical imaging, and interpretation at multiple scales. , 1997, Electroencephalography and clinical neurophysiology.

[19]  B.D. Rao,et al.  Comparison of basis selection methods , 1996, Conference Record of The Thirtieth Asilomar Conference on Signals, Systems and Computers.

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

[21]  K.-R. Muller,et al.  The Berlin brain-computer interface: EEG-based communication without subject training , 2006, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[22]  Bin He,et al.  Electrophysiological Imaging of Brain Activity and Connectivity—Challenges and Opportunities , 2011, IEEE Transactions on Biomedical Engineering.

[23]  J. Capon High-resolution frequency-wavenumber spectrum analysis , 1969 .

[24]  Ingo Fründ,et al.  Stimulus intensity affects early sensory processing: sound intensity modulates auditory evoked gamma-band activity in human EEG. , 2007, International Journal of Psychophysiology.

[25]  A. Lee Swindlehurst,et al.  Matching Pursuit and Source Deflation for Sparse EEG/MEG Dipole Moment Estimation , 2013, IEEE Transactions on Biomedical Engineering.

[26]  Giovanni B. Frisoni,et al.  EEG markers are associated to gray matter changes in thalamus and basal ganglia in subjects with mild cognitive impairment , 2012, NeuroImage.

[27]  Athina P. Petropulu,et al.  EEG sparse source localization via Range Space Rotation , 2015, 2015 IEEE 6th International Workshop on Computational Advances in Multi-Sensor Adaptive Processing (CAMSAP).

[28]  R. O. Schmidt,et al.  Multiple emitter location and signal Parameter estimation , 1986 .

[29]  Joel A. Tropp,et al.  Signal Recovery From Random Measurements Via Orthogonal Matching Pursuit , 2007, IEEE Transactions on Information Theory.

[30]  G. Jacobson,et al.  Auditory evoked gamma band potential in normal subjects. , 1997, Journal of the American Academy of Audiology.

[31]  Stephen P. Boyd,et al.  Enhancing Sparsity by Reweighted ℓ1 Minimization , 2007, 0711.1612.

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

[33]  Jing Zhang,et al.  Auditory cortical responses evoked by pure tones in healthy and sensorineural hearing loss subjects: functional MRI and magnetoencephalography. , 2006, Chinese medical journal.

[34]  R. DeVore,et al.  Compressed sensing and best k-term approximation , 2008 .

[35]  Ruimin Wang,et al.  Classification of Four-Class Motor Imagery Employing Single-Channel Electroencephalography , 2014, PloS one.

[36]  Barry D. Van Veen,et al.  MEG and EEG source localization in beamspace , 2006, IEEE Transactions on Biomedical Engineering.

[37]  K. Whittingstall,et al.  Dipole localization accuracy using grand-average EEG data sets , 2004, Clinical Neurophysiology.

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

[39]  Lei Ding,et al.  Motor imagery classification by means of source analysis for brain–computer interface applications , 2004, Journal of neural engineering.

[40]  Julius P. A. Dewald,et al.  Evaluation of different cortical source localization methods using simulated and experimental EEG data , 2005, NeuroImage.

[41]  Richard M. Leahy,et al.  Brainstorm: A User-Friendly Application for MEG/EEG Analysis , 2011, Comput. Intell. Neurosci..

[42]  E.J. Candes Compressive Sampling , 2022 .

[43]  Michael A. Saunders,et al.  Atomic Decomposition by Basis Pursuit , 1998, SIAM J. Sci. Comput..

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

[45]  Bart Vanrumste,et al.  Journal of Neuroengineering and Rehabilitation Open Access Review on Solving the Inverse Problem in Eeg Source Analysis , 2022 .

[46]  Bin He,et al.  EEG Source Imaging Enhances the Decoding of Complex Right-Hand Motor Imagery Tasks , 2016, IEEE Transactions on Biomedical Engineering.

[47]  Ander Ramos-Murguialday,et al.  Classification of different reaching movements from the same limb using EEG , 2017, Journal of neural engineering.

[48]  Georgios B. Giannakis,et al.  Weighted sparse signal decomposition , 2012, 2012 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[49]  Ning Jiang,et al.  Enhanced Low-Latency Detection of Motor Intention From EEG for Closed-Loop Brain-Computer Interface Applications , 2014, IEEE Transactions on Biomedical Engineering.

[50]  S. Posse,et al.  Intensity coding of auditory stimuli: an fMRI study , 1998, Neuropsychologia.

[51]  Laleh Najafizadeh,et al.  Capturing dynamic patterns of task-based functional connectivity with EEG , 2013, NeuroImage.

[52]  A. Engel,et al.  Cognitive functions of gamma-band activity: memory match and utilization , 2004, Trends in Cognitive Sciences.

[53]  Bhaskar D. Rao,et al.  Sparse Bayesian learning for basis selection , 2004, IEEE Transactions on Signal Processing.

[54]  R. Pascual-Marqui Review of methods for solving the EEG inverse problem , 1999 .

[55]  Yun-Bin Zhao,et al.  RSP-Based Analysis for Sparsest and Least $\ell_1$-Norm Solutions to Underdetermined Linear Systems , 2013, IEEE Transactions on Signal Processing.

[56]  Clemens Brunner,et al.  Mu rhythm (de)synchronization and EEG single-trial classification of different motor imagery tasks , 2006, NeuroImage.

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

[58]  Bin He,et al.  Dynamic imaging of ictal oscillations using non-invasive high-resolution EEG , 2011, NeuroImage.