Daniel Strohmeier

MEG and EEG data analysis with MNE-Python

Magnetoencephalography and electroencephalography (M/EEG) measure the weak electromagnetic signals generated by neuronal activity in the brain. Using these signals to characterize and locate neural activation in the brain is a challenge that requires expertise in physics, signal processing, statistics, and numerical methods. As part of the MNE software suite, MNE-Python is an open-source software package that addresses this challenge by providing state-of-the-art algorithms implemented in Python that cover multiple methods of data preprocessing, source localization, statistical analysis, and estimation of functional connectivity between distributed brain regions. All algorithms and utility functions are implemented in a consistent manner with well-documented interfaces, enabling users to create M/EEG data analysis pipelines by writing Python scripts. Moreover, MNE-Python is tightly integrated with the core Python libraries for scientific comptutation (NumPy, SciPy) and visualization (matplotlib and Mayavi), as well as the greater neuroimaging ecosystem in Python via the Nibabel package. The code is provided under the new BSD license allowing code reuse, even in commercial products. Although MNE-Python has only been under heavy development for a couple of years, it has rapidly evolved with expanded analysis capabilities and pedagogical tutorials because multiple labs have collaborated during code development to help share best practices. MNE-Python also gives easy access to preprocessed datasets, helping users to get started quickly and facilitating reproducibility of methods by other researchers. Full documentation, including dozens of examples, is available at http://martinos.org/mne.

MNE software for processing MEG and EEG data

Magnetoencephalography and electroencephalography (M/EEG) measure the weak electromagnetic signals originating from neural currents in the brain. Using these signals to characterize and locate brain activity is a challenging task, as evidenced by several decades of methodological contributions. MNE, whose name stems from its capability to compute cortically-constrained minimum-norm current estimates from M/EEG data, is a software package that provides comprehensive analysis tools and workflows including preprocessing, source estimation, time-frequency analysis, statistical analysis, and several methods to estimate functional connectivity between distributed brain regions. The present paper gives detailed information about the MNE package and describes typical use cases while also warning about potential caveats in analysis. The MNE package is a collaborative effort of multiple institutes striving to implement and share best methods and to facilitate distribution of analysis pipelines to advance reproducibility of research. Full documentation is available at http://martinos.org/mne.

Time-Frequency Mixed-Norm Estimates: Sparse M/EEG imaging with non-stationary source activations

Magnetoencephalography (MEG) and electroencephalography (EEG) allow functional brain imaging with high temporal resolution. While solving the inverse problem independently at every time point can give an image of the active brain at every millisecond, such a procedure does not capitalize on the temporal dynamics of the signal. Linear inverse methods (minimum-norm, dSPM, sLORETA, beamformers) typically assume that the signal is stationary: regularization parameter and data covariance are independent of time and the time varying signal-to-noise ratio (SNR). Other recently proposed non-linear inverse solvers promoting focal activations estimate the sources in both space and time while also assuming stationary sources during a time interval. However such a hypothesis holds only for short time intervals. To overcome this limitation, we propose time-frequency mixed-norm estimates (TF-MxNE), which use time-frequency analysis to regularize the ill-posed inverse problem. This method makes use of structured sparse priors defined in the time-frequency domain, offering more accurate estimates by capturing the non-stationary and transient nature of brain signals. State-of-the-art convex optimization procedures based on proximal operators are employed, allowing the derivation of a fast estimation algorithm. The accuracy of the TF-MxNE is compared with recently proposed inverse solvers with help of simulations and by analyzing publicly available MEG datasets.

Functional brain imaging with M/EEG using structured sparsity in time-frequency dictionaries

Magnetoencephalography (MEG) and electroencephalography (EEG) allow functional brain imaging with high temporal resolution. While time-frequency analysis is often used in the field, it is not commonly employed in the context of the ill-posed inverse problem that maps the MEG and EEG measurements to the source space in the brain. In this work, we detail how convex structured sparsity can be exploited to achieve a principled and more accurate functional imaging approach. Importantly, time-frequency dictionaries can capture the non-stationary nature of brain signals and state-of-the-art convex optimization procedures based on proximal operators allow the derivation of a fast estimation algorithm. We compare the accuracy of our new method to recently proposed inverse solvers with help of simulations and analysis of real MEG data.

Real-Time MEG Source Localization Using Regional Clustering

With its millisecond temporal resolution, Magnetoencephalography (MEG) is well suited for real-time monitoring of brain activity. Real-time feedback allows the adaption of the experiment to the subject’s reaction and increases time efficiency by shortening acquisition and off-line analysis. Two formidable challenges exist in real-time analysis: the low signal-to-noise ratio (SNR) and the limited time available for computations. Since the low SNR reduces the number of distinguishable sources, we propose an approach which downsizes the source space based on a cortical atlas and allows to discern the sources in the presence of noise. Each cortical region is represented by a small set of dipoles, which is obtained by a clustering algorithm. Using this approach, we adapted dynamic statistical parametric mapping for real-time source localization. In terms of point spread and crosstalk between regions the proposed clustering technique performs better than selecting spatially evenly distributed dipoles. We conducted real-time source localization on MEG data from an auditory experiment. The results demonstrate that the proposed real-time method localizes sources reliably in the superior temporal gyrus. We conclude that real-time source estimation based on MEG is a feasible, useful addition to the standard on-line processing methods, and enables feedback based on neural activity during the measurements.

Improved MEG/EEG source localization with reweighted mixed-norms

MEG/EEG source imaging allows for the noninvasive analysis of brain activity with high temporal and good spatial resolution. As the bioelectromagnetic inverse problem is ill-posed, a priori information is required to find a unique source estimate. For the analysis of evoked brain activity, spatial sparsity of the neuronal activation can be assumed. Due to the convexity, ℓ-norm based constraints are often used for this, which however lead to source estimates biased in amplitude and often suboptimal in terms of source selection. As an alternative, non-convex regularization functionals such as ℓ p-quasinorms with 0 <; p <; 1 can be used. In this work, we present a MEG/EEG inverse solver based on a ℓ 2,0.5-quasinorm penalty promoting spatial sparsity as well as temporal stationarity of the brain activity. For solving the resulting non-convex optimization problem, we propose the iterative reweighted Mixed Norm Estimate, which is based on reweighted convex optimization and combines a block coordinate descent scheme and an active set strategy to solve each surrogate problem efficiently. We provide empirical evidence based on simulations and analysis of MEG data that the proposed method outperforms the standard Mixed Norm Estimate in terms of active source identification and amplitude bias.

The Iterative Reweighted Mixed-Norm Estimate for Spatio-Temporal MEG/EEG Source Reconstruction

Source imaging based on magnetoencephalography (MEG) and electroencephalography (EEG) allows for the non-invasive analysis of brain activity with high temporal and good spatial resolution. As the bioelectromagnetic inverse problem is ill-posed, constraints are required. For the analysis of evoked brain activity, spatial sparsity of the neuronal activation is a common assumption. It is often taken into account using convex constraints based on the l1-norm. The resulting source estimates are however biased in amplitude and often suboptimal in terms of source selection due to high correlations in the forward model. In this work, we demonstrate that an inverse solver based on a block-separable penalty with a Frobenius norm per block and a l0.5-quasinorm over blocks addresses both of these issues. For solving the resulting non-convex optimization problem, we propose the iterative reweighted Mixed Norm Estimate (irMxNE), an optimization scheme based on iterative reweighted convex surrogate optimization problems, which are solved efficiently using a block coordinate descent scheme and an active set strategy. We compare the proposed sparse imaging method to the dSPM and the RAP-MUSIC approach based on two MEG data sets. We provide empirical evidence based on simulations and analysis of MEG data that the proposed method improves on the standard Mixed Norm Estimate (MxNE) in terms of amplitude bias, support recovery, and stability.

SPHARA - A Generalized Spatial Fourier Analysis for Multi-Sensor Systems with Non-Uniformly Arranged Sensors: Application to EEG

Important requirements for the analysis of multichannel EEG data are efficient techniques for signal enhancement, signal decomposition, feature extraction, and dimensionality reduction. We propose a new approach for spatial harmonic analysis (SPHARA) that extends the classical spatial Fourier analysis to EEG sensors positioned non-uniformly on the surface of the head. The proposed method is based on the eigenanalysis of the discrete Laplace-Beltrami operator defined on a triangular mesh. We present several ways to discretize the continuous Laplace-Beltrami operator and compare the properties of the resulting basis functions computed using these discretization methods. We apply SPHARA to somatosensory evoked potential data from eleven volunteers and demonstrate the ability of the method for spatial data decomposition, dimensionality reduction and noise suppression. When employing SPHARA for dimensionality reduction, a significantly more compact representation can be achieved using the FEM approach, compared to the other discretization methods. Using FEM, to recover 95% and 99% of the total energy of the EEG data, on average only 35% and 58% of the coefficients are necessary. The capability of SPHARA for noise suppression is shown using artificial data. We conclude that SPHARA can be used for spatial harmonic analysis of multi-sensor data at arbitrary positions and can be utilized in a variety of other applications.

Reconstruction of quasi-radial dipolar activity using three-component magnetic field measurements

OBJECTIVE While standard magnetoencephalographic systems record only one component of the biomagnetic field, novel vector-biomagnetometers enable measurement of all three components of the field at each sensing point. Because information content in standard one-component magnetoencephalography (MEG) is often not adequate to reconstruct quasi-radial dipolar activity, we tested the hypothesis that quasi-radial activity can be estimated using three-component MEG. METHODS We stimulated the right median nerve in 11 healthy volunteers and recorded the somatosensory evoked fields over the contralateral hemisphere using a novel vector-biomagnetometer system comprised of SQUID-based magnetometer triplets. Source reconstruction for the early cortical components N20m and P25m was subsequently performed. RESULTS Both tangential and quasi-radial dipolar activity could be reconstructed in 10 of the 11 participants. Dipole locations were found in the vicinity of the central sulcus, and dipole orientations were predominantly tangential for N20m and quasi-radial for P25m. The mean location difference between the tangential and quasi-radial dipoles was 11.9 mm and the mean orientation difference was 97.5°. CONCLUSIONS Quasi-radial dipolar activity can be reconstructed from three-component magnetoencephalographic measurements. SIGNIFICANCE Three-component MEG provides higher information content than does standard MEG.

Rod Driven Frequency Entrainment and Resonance Phenomena

A controversy exists on photic driving in the human visual cortex evoked by intermittent photic stimulation. Frequency entrainment and resonance phenomena are reported for frequencies higher than 12 Hz in some studies while missing in others. We hypothesized that this might be due to different experimental conditions, since both high and low intensity light stimulation were used. However, most studies do not report radiometric measurements, which makes it impossible to categorize the stimulation according to photopic, mesopic, and scotopic vision. Low intensity light stimulation might lead to scotopic vision, where rod perception dominates. In this study, we investigated photic driving for rod-dominated visual input under scotopic conditions. Twelve healthy volunteers were stimulated with low intensity light flashes at 20 stimulation frequencies, leading to rod activation only. The frequencies were multiples of the individual alpha frequency (α) of each volunteer in the range from 0.40 to 2.30∗α. Three hundred and six-channel whole head magnetoencephalography recordings were analyzed in time, frequency, and spatiotemporal domains with the Topographic Matching Pursuit algorithm. We found resonance phenomena and frequency entrainment for stimulations at or close to the individual alpha frequency (0.90–1.10∗α) and half of the alpha frequency (0.40–0.55∗α). No signs of resonance and frequency entrainment phenomena were revealed around 2.00∗α. Instead, on-responses at the beginning and off-responses at the end of each stimulation train were observed for the first time in a photic driving experiment at frequencies of 1.30–2.30∗α, indicating that the flicker fusion threshold was reached. All results, the resonance and entrainment as well as the fusion effects, provide evidence for rod-dominated photic driving in the visual cortex.