https://hal-imt.archives-ouvertes.fr/hal-01183552Strohmeier, DanielDanielStrohmeierTU - Technische Universität IlmenauGramfort, AlexandreAlexandreGramfortLTCI - Laboratoire Traitement et Communication de l'Information - Télécom ParisTech - IMT - Institut Mines-Télécom [Paris] - CNRS - Centre National de la Recherche ScientifiqueHaueisen, JensJensHaueisenTU - Technische Universität IlmenauMEG/EEG source imaging with a non-convex penalty in the time-frequency domainHAL CCSD2015MEGEEGinverse problemstructured sparsityGabor transformiterative reweighted optimization algorithm[INFO] Computer Science [cs][INFO.INFO-IM] Computer Science [cs]/Medical Imaging[MATH.MATH-OC] Mathematics [math]/Optimization and Control [math.OC]Gramfort, AlexandreRecherche collaborative en neurosciences computationnelles - Revealing thalamocortical interactions in humans with integrated MEG/EEG, human intracranial recordings and computational neural modeling - - THALAMEEG2014 - ANR-14-NEUC-0002 - CRCNS - VALID - 2015-08-11 17:28:002020-12-08 09:48:202015-08-12 14:09:33enConference papersapplication/pdf1Due to the excellent temporal resolution, MEG/EEG source imaging is an important measurement modality to study dynamic processes in the brain. As the bioelectromagnetic inverse problem is ill-posed, constraints have to be imposed on the source estimates to find a unique solution. These constraints can be applied either in the standard or a transformed domain. The Time-Frequency Mixed Norm Estimate applies a composite convex regularization functional promoting structured sparsity in the time-frequency domain by combining an L2,1-mixed-norm and an 1-norm penalty on the coefficients of the Gabor TF decomposition of the source signals, to improve the reconstruction of spatially sparse neural activations with non-stationary and transient signals. Due to the L1-norm based constraints, the resulting source estimates are however biased in amplitude and often suboptimal in terms of source selection. In this work, we present the iterative reweighted Time-Frequency Mixed Norm Estimate, which employs a composite non-convex penalty formed by the sum of an L2,0.5-quasinorm and an 0.5-quasinorm penalty. The resulting non-convex problem is solved with a reweighted convex optimization scheme, in which each iteration is equivalent to a weighted Time-Frequency Mixed-Norm Estimate solved efficiently using a block coordinate descent scheme and an active set strategy. We compare our approach to alternative solvers using simulations and analysis of MEG data and demonstrate the benefit of the iterative reweighted Time-Frequency Mixed Norm Estimate with regard to active source identification, amplitude bias correction, and temporal unmixing of activations.