M/EEG source reconstruction based on Gabor thresholding in the source space

Daniel Strohmeier (1) , Alexandre Gramfort (2, 3, 4) , Jens Haueisen (1) , Matti Hämäläinen (3) , Matthieu Kowalski (5)

1 TU - Technische Universität Ilmenau
2 PARIETAL - Modelling brain structure, function and variability based on high-field MRI data
3 Athinoula A. Martinos Center for Biomedical Imaging
4 LNAO - Laboratoire de Neuroimagerie Assistée par Ordinateur
5 L2S - Laboratoire des signaux et systèmes

Thanks to their high temporal resolution, source reconstruction based on Magnetoencephalography (MEG) and/or Electroencephalography (EEG) is an important tool for noninvasive functional brain imaging. Since the MEG/EEG inverse problem is ill-posed, inverse solvers employ priors on the sources. While priors are generally applied in the time domain, the time-frequency (TF) characteristics of brain signals are rarely employed as a spatio-temporal prior. In this work, we present an inverse solver which employs a structured sparse prior formed by the sum of $\ell_{21}$ and $\ell_{1}$ norms on the coefficients of the Gabor TF decomposition of the source activations. The resulting convex optimization problem is solved using a first-order scheme based on proximal operators. We provide empirical evidence based on EEG simulations that the proposed method is able to recover neural activations that are spatially sparse, temporally smooth and non-stationary. We compare our approach to alternative solvers based also on convex sparse priors, and demonstrate the benefit of promoting sparse Gabor decompositions via a mathematically principled iterative thresholding procedure.

Domaines: Traitement du signal et de l'image [eess.SP]; Imagerie médicale; Traitement du signal et de l'image [eess.SP]

Daniel Strohmeier, Alexandre Gramfort, Jens Haueisen, Matti Hämäläinen, Matthieu Kowalski. M/EEG source reconstruction based on Gabor thresholding in the source space. NFSI 2011, May 2011, Banff, Canada. ⟨hal-00621044⟩