Robust Empirical Bayesian Reconstruction of Distributed Sources for Electromagnetic Brain Imaging

Structured noise champagne: an empirical Bayesian algorithm for electromagnetic brain imaging with structured noise

Introduction

Electromagnetic brain imaging is the reconstruction of brain activity from non-invasive recordings of electroencephalography (EEG), magnetoencephalography (MEG), and also from invasive ones such as the intracranial recording of electrocorticography (ECoG), intracranial electroencephalography (iEEG), and stereo electroencephalography EEG (sEEG). These modalities are widely used techniques to study the function of the human brain. Efficient reconstruction of electrophysiological activity of neurons in the brain from EEG/MEG measurements is important for neuroscience research and clinical applications. An enduring challenge in this field is the accurate inference of brain signals of interest while accounting for all sources of noise that contribute to the sensor measurements. The statistical characteristic of the noise plays a crucial role in the success of the brain source recovery process, which can be formulated as a sparse regression problem.

Method

In this study, we assume that the dominant environment and biological sources of noise that have high spatial correlations in the sensors can be expressed as a structured noise model based on the variational Bayesian factor analysis. To the best of our knowledge, no existing algorithm has addressed the brain source estimation problem with such structured noise. We propose to apply a robust empirical Bayesian framework for iteratively estimating the brain source activity and the statistics of the structured noise. In particular, we perform inference of the variational Bayesian factor analysis (VBFA) noise model iteratively in conjunction with source reconstruction.

Results

To demonstrate the effectiveness of the proposed algorithm, we perform experiments on both simulated and real datasets. Our algorithm achieves superior performance as compared to several existing benchmark algorithms.

Discussion

A key aspect of our algorithm is that we do not require any additional baseline measurements to estimate the noise covariance from the sensor data under scenarios such as resting state analysis, and other use cases wherein a noise or artifactual source occurs only in the active period but not in the baseline period (e.g., neuro-modulatory stimulation artifacts and speech movements).

Flexible Alternating Projection for Spatially Extended Brain Source Localization

OBJECTIVE

Understanding neural sources behind MEG and EEG signals is significant for basic and clinical neuroscience. Existing techniques, such as Recursively Applied and Projected Multiple Signal Classification (RAP-MUSIC) and Alternating Projection (AP), rely on limited current dipoles, representing focal sources with zero spatial extent. However, this oversimplifies realistic neural activity, which exhibits varying spatial extents.

METHODS

To address this, we enhanced the AP approach, creating FLEX-AP, capable of localizing discrete and extended sources. FLEX-AP simultaneously optimizes location and spatial extent of candidate sources.

RESULTS

FLEX-AP demonstrated superior performance, reducing localization errors compared to AP, RAP-MUSIC, and FLEX-RAP-MUSIC, particularly with extended sources. Moreover, FLEX-AP exhibited enhanced robustness against modeling errors in realistic scenarios. Applying FLEX-AP to MEG recordings of auditory responses validated its effectiveness, underscoring potential in advancing neuroscientific investigations.

CONCLUSION

FLEX-AP offers a robust, flexible framework for M/EEG source localization, overcoming limitations of simplistic zero-extent dipole models. By accurately estimating position and spatial extent of neural sources, FLEX-AP bridges the gap between theoretical models and realistic activity, demonstrating utility in simulated and real-world scenarios.

SIGNIFICANCE

FLEX-AP advances source imaging techniques with implications for basic neuroscience and clinical applications. Modeling extended sources precisely enables more accurate investigations of brain dynamics, potentially improving diagnostic and therapeutic approaches for neurological and psychiatric disorders.

Robust interpolation of EEG/MEG sensor time-series via electromagnetic source imaging

Objective.electroencephalography (EEG) and magnetoencephalography (MEG) are widely used non-invasive techniques in clinical and cognitive neuroscience. However, low spatial resolution measurements, partial brain coverage by some sensor arrays, as well as noisy sensors could result in distorted sensor topographies resulting in inaccurate reconstructions of underlying brain dynamics. Solving these problems has been a challenging task. This paper proposes a robust framework based on electromagnetic source imaging for interpolation of unknown or poor quality EEG/MEG measurements.Approach.This framework consists of two steps: (1) estimating brain source activity using a robust inverse algorithm along with the leadfield matrix of available good sensors, and (2) interpolating unknown or poor quality EEG/MEG measurements using the reconstructed brain sources using the leadfield matrices of unknown or poor quality sensors. We evaluate the proposed framework through simulations and several real datasets, comparing its performance to two popular benchmarks-neighborhood interpolation and spherical spline interpolation algorithms.Results.In both simulations and real EEG/MEG measurements, we demonstrate several advantages compared to benchmarks, which are robust to highly correlated brain activity, low signal-to-noise ratio data and accurately estimates cortical dynamics.Significance.These results demonstrate a rigorous platform to enhance the spatial resolution of EEG and MEG, to overcome limitations of partial coverage of EEG/MEG sensor arrays that is particularly relevant to low-channel count optically pumped magnetometer arrays, and to estimate activity in poor/noisy sensors to a certain extent based on the available measurements from other good sensors. Implementation of this framework will enhance the quality of EEG and MEG, thereby expanding the potential applications of these modalities.

XDL-ESI: Electrophysiological Sources Imaging via explainable deep learning framework with validation on simultaneous EEG and iEEG

Electroencephalography (EEG) or Magnetoencephalography (MEG) source imaging aims to estimate the underlying activated brain sources to explain the observed EEG/MEG recordings. Solving the inverse problem of EEG/MEG Source Imaging (ESI) is challenging due to its ill-posed nature. To achieve a unique solution, it is essential to apply sophisticated regularization constraints to restrict the solution space. Traditionally, the design of regularization terms is based on assumptions about the spatiotemporal structure of the underlying source dynamics. In this paper, we propose a novel paradigm for ESI via an Explainable Deep Learning framework, termed as XDL-ESI, which connects the iterative optimization algorithm with deep learning architecture by unfolding the iterative updates with neural network modules. The proposed framework has the advantages of (1) establishing a data-driven approach to model the source solution structure instead of using hand-crafted regularization terms; (2) improving the robustness of source solutions by introducing a topological loss that leverages the geometric spatial information applying varying penalties on distinct localization errors; (3) improving the reconstruction efficiency and interpretability as it inherits the advantages from both the iterative optimization algorithms (interpretability) and deep learning approaches (function approximation). The proposed XDL-ESI framework provides an efficient, accurate, and interpretable paradigm to solve the ESI inverse problem with satisfactory performance in both simulated data and real clinical data. Specially, this approach is further validated using simultaneous EEG and intracranial EEG (iEEG).

Source imaging method based on diagonal covariance bases and its applications to OPM-MEG

Magnetoencephalography (MEG) is a noninvasive imaging technique used in neuroscience and clinical research. The source estimation of MEG involves solving a highly underdetermined inverse problem, which requires additional constraints to restrict the solution space. Traditional methods tend to obscure the extent of the sources. However, an accurate estimation of the source extent is important for studying brain activity or preoperatively estimating pathogenic regions. To improve the estimation accuracy of the extended source extent, the spatial constraint of sources is employed in the Bayesian framework. For example, the source is decomposed into a linear combination of validated spatial basis functions, which is proved to improve the source imaging accuracy. In this work, we further construct the spatial properties of the source using the diagonal covariance bases (DCB), which we summarize as the source imaging method SI-DCB. In this approach, specifically, the covariance matrix of the spatial coefficients is modeled as a weighted combination of diagonal covariance basis functions. The convex analysis is used to estimate noise and model parameters under the Bayesian framework. Extensive numerical simulations showed that SI-DCB outperformed five benchmark methods in accurately estimating the location and extent of patch sources. The effectiveness of SI-DCB was verified through somatosensory stimulation experiments performed on a 31-channel OPM-MEG system. The SI-DCB correctly identified the source area where each brain response occurred. The superior performance of SI-DCB suggests that it can provide a template approach for improving the accuracy of source extent estimations under a sparse Bayesian framework.

Rejuvenating classical brain electrophysiology source localization methods with spatial graph Fourier filters for source extents estimation

EEG/MEG source imaging (ESI) aims to find the underlying brain sources to explain the observed EEG or MEG measurement. Multiple classical approaches have been proposed to solve the ESI problem based on different neurophysiological assumptions. To support clinical decision-making, it is important to estimate not only the exact location of the source signal but also the extended source activation regions. Existing methods may render over-diffuse or sparse solutions, which limit the source extent estimation accuracy. In this work, we leverage the graph structures defined in the 3D mesh of the brain and the spatial graph Fourier transform (GFT) to decompose the spatial graph structure into sub-spaces of low-, medium-, and high-frequency basis. We propose to use the low-frequency basis of spatial graph filters to approximate the extended areas of brain activation and embed the GFT into the classical ESI methods. We validated the classical source localization methods with the corresponding improved version using GFT in both synthetic data and real data. We found the proposed method can effectively reconstruct focal source patterns and significantly improve the performance compared to the classical algorithms.

Joint Learning of Full-Structure Noise in Hierarchical Bayesian Regression Models

We consider the reconstruction of brain activity from electroencephalography (EEG). This inverse problem can be formulated as a linear regression with independent Gaussian scale mixture priors for both the source and noise components. Crucial factors influencing the accuracy of the source estimation are not only the noise level but also its correlation structure, but existing approaches have not addressed the estimation of noise covariance matrices with full structure. To address this shortcoming, we develop hierarchical Bayesian (type-II maximum likelihood) models for observations with latent variables for source and noise, which are estimated jointly from data. As an extension to classical sparse Bayesian learning (SBL), where across-sensor observations are assumed to be independent and identically distributed, we consider Gaussian noise with full covariance structure. Using the majorization-maximization framework and Riemannian geometry, we derive an efficient algorithm for updating the noise covariance along the manifold of positive definite matrices. We demonstrate that our algorithm has guaranteed and fast convergence and validate it in simulations and with real MEG data. Our results demonstrate that the novel framework significantly improves upon state-of-the-art techniques in the real-world scenario where the noise is indeed non-diagonal and full-structured. Our method has applications in many domains beyond biomagnetic inverse problems.

Deep learning based source imaging provides strong sublobar localization of epileptogenic zone from MEG interictal spikes

Electromagnetic source imaging (ESI) offers unique capability of imaging brain dynamics for studying brain functions and aiding the clinical management of brain disorders. Challenges exist in ESI due to the ill-posedness of the inverse problem and thus the need of modeling the underlying brain dynamics for regularizations. Advances in generative models provide opportunities for more accurate and realistic source modeling that could offer an alternative approach to ESI for modeling the underlying brain dynamics beyond equivalent physical source models. However, it is not straightforward to explicitly formulate the knowledge arising from these generative models within the conventional ESI framework. Here we investigate a novel source imaging framework based on mesoscale neuronal modeling and deep learning (DL) that can learn the sensor-source mapping relationship directly from MEG data for ESI. Two DL-based ESI models were trained based on data generated by neural mass models and either generic or personalized head models. The robustness of the two DL models was evaluated by systematic computer simulations and clinical validation in a cohort of 29 drug-resistant focal epilepsy patients who underwent intracranial EEG (iEEG) evaluation or surgical resection. Results estimated from pre-operative MEG interictal spikes were quantified using the overlap with resection regions and the distance to the seizure-onset zone (SOZ) defined by iEEG recordings. The DL-based ESI provided robust results when no personalized head geometry is considered, reaching a spatial dispersion of 21.90 ± 19.03 mm, sublobar concordance of 83 %, and sublobar sensitivity and specificity of 66 and 97 % respectively. When using personalized head geometry derived from individual patients' MRI in the training data, personalized DL-based ESI model can further improve the performance and reached a spatial dispersion of 8.19 ± 8.14 mm, sublobar concordance of 93 %, and sublobar sensitivity and specificity of 77 and 99 % respectively. When compared to the SOZ, the localization error of the personalized approach is 15.78 ± 5.54 mm, outperforming the conventional benchmarks. This work demonstrates that combining generative models and deep learning enables an accurate and robust imaging of epileptogenic zone from MEG recordings with strong sublobar precision, suggesting its added value to enhancing MEG source localization and imaging, and to epilepsy source localization and other clinical applications.

Bayesian Adaptive Beamformer for Robust Electromagnetic Brain Imaging of Correlated Sources in High Spatial Resolution

Reconstructing complex brain source activity at a high spatiotemporal resolution from magnetoencephalography (MEG) or electroencephalography (EEG) remains a challenging problem. Adaptive beamformers are routinely deployed for this imaging domain using the sample data covariance. However adaptive beamformers have long been hindered by 1) high degree of correlation between multiple brain sources, and 2) interference and noise embedded in sensor measurements. This study develops a novel framework for minimum variance adaptive beamformers that uses a model data covariance learned from data using a sparse Bayesian learning algorithm (SBL-BF). The learned model data covariance effectively removes influence from correlated brain sources and is robust to noise and interference without the need for baseline measurements. A multiresolution framework for model data covariance computation and parallelization of the beamformer implementation enables efficient high-resolution reconstruction images. Results with both simulations and real datasets indicate that multiple highly correlated sources can be accurately reconstructed, and that interference and noise can be sufficiently suppressed. Reconstructions at 2-2.5mm resolution (  ∼  150K voxels) are possible with efficient run times of 1-3 minutes. This novel adaptive beamforming algorithm significantly outperforms the state-of-the-art benchmarks. Therefore, SBL-BF provides an effective framework for efficiently reconstructing multiple correlated brain sources with high resolution and robustness to interference and noise.

Bayesian Algorithms for Joint Estimation of Brain Activity and Noise in Electromagnetic Imaging

Simultaneously estimating brain source activity and noise has long been a challenging task in electromagnetic brain imaging using magneto- and electroencephalography. The problem is challenging not only in terms of solving the NP-hard inverse problem of reconstructing unknown brain activity across thousands of voxels from a limited number of sensors, but also for the need to simultaneously estimate the noise and interference. We present a generative model with an augmented leadfield matrix to simultaneously estimate brain source activity and sensor noise statistics in electromagnetic brain imaging (EBI). We then derive three Bayesian inference algorithms for this generative model (expectation-maximization (EBI-EM), convex bounding (EBI-Convex) and fixed-point (EBI-Mackay)) to simultaneously estimate the hyperparameters of the prior distribution for brain source activity and sensor noise. A comprehensive performance evaluation for these three algorithms is performed. Simulations consistently show that the performance of EBI-Convex and EBI-Mackay updates is superior to that of EBI-EM. In contrast to the EBI-EM algorithm, both EBI-Convex and EBI-Mackay updates are quite robust to initialization, and are computationally efficient with fast convergence in the presence of both Gaussian and real brain noise. We also demonstrate that EBI-Convex and EBI-Mackay update algorithms can reconstruct complex brain activity with only a few trials of sensor data, and for resting-state data, achieving significant improvement in source reconstruction and noise learning for electromagnetic brain imaging.

Empirical Bayesian localization of event-related time-frequency neural activity dynamics

Accurate reconstruction of the spatio-temporal dynamics of event-related cortical oscillations across human brain regions is an important problem in functional brain imaging and human cognitive neuroscience with magnetoencephalography (MEG) and electroencephalography (EEG). The problem is challenging not only in terms of localization of complex source configurations from sensor measurements with unknown noise and interference but also for reconstruction of transient event-related time-frequency dynamics of cortical oscillations. We recently proposed a robust empirical Bayesian algorithm for simultaneous reconstruction of complex brain source activity and noise covariance, in the context of evoked and resting-state data. In this paper, we expand upon this empirical Bayesian framework for optimal reconstruction of event-related time-frequency dynamics of regional cortical oscillations, referred to as time-frequency Champagne (TFC). This framework enables imaging of five-dimensional (space, time, and frequency) event-related brain activity from M/EEG data, and can be viewed as a time-frequency optimized adaptive Bayesian beamformer. We evaluate TFC in both simulations and several real datasets, with comparisons to benchmark standards - variants of time-frequency optimized adaptive beamformers (TFBF) as well as the sLORETA algorithm. In simulations, we demonstrate several advantages in estimating time-frequency cortical oscillatory dynamics compared to benchmarks. With real MEG data, we demonstrate across many datasets that the proposed approach is robust to highly correlated brain activity and low SNR data, and is able to accurately reconstruct cortical dynamics with data from just a few epochs.

Electromagnetic Source Imaging via Bayesian Modeling with Smoothness in Spatial and Temporal Domains

Accurate reconstruction of cortical activation from electroencephalography and magnetoencephalography (E/MEG) is a long-standing challenge because of the inherently ill-posed inverse problem. In this paper, a novel algorithm under the empirical Bayesian framework, source imaging with smoothness in spatial and temporal domains (SI-SST), is proposed to address this issue. In SI-SST, current sources are decomposed into the product of spatial smoothing kernel, sparseness encoding coefficients, and temporal basis functions (TBFs). Further smoothness is integrated in the temporal domain with the employment of an underlying autoregressive model. Because sparseness encoding coefficients are constructed depending on overlapped clusters over cortex in this model, we derived a novel update rule based on fixed-point criterion instead of the convexity based approach which becomes invalid in this scenario. Entire variables and hyper parameters are updated alternatively in the variational inference procedure. SI-SST was assessed by multiple metrics with both simulated and experimental datasets. In practice, SI-SST had the superior reconstruction performance in both spatial extents and temporal profiles compared to the benchmarks.

Imaging somatosensory cortex responses measured by OPM-MEG: Variational free energy-based spatial smoothing estimation approach

In recent years, optically pumped magnetometer (OPM)-based magnetoencephalography (MEG) has shown potential for analyzing brain activity. It has a flexible sensor configuration and comparable sensitivity to conventional SQUID-MEG. We constructed a 32-channel OPM-MEG system and used it to measure cortical responses to median and ulnar nerve stimulations. Traditional magnetic source imaging methods tend to blur the spatial extent of sources. Accurate estimation of the spatial size of the source is important for studying the organization of brain somatotopy and for pre-surgical functional mapping. We proposed a new method called variational free energy-based spatial smoothing estimation (FESSE) to enhance the accuracy of mapping somatosensory cortex responses. A series of computer simulations based on the OPM-MEG showed better performance than the three types of competing methods under different levels of signal-to-noise ratios, source patch sizes, and co-registration errors. FESSE was then applied to the source imaging of the OPM-MEG experimental data.

Clinical Validation of the Champagne Algorithm for Epilepsy Spike Localization

Magnetoencephalography (MEG) is increasingly used for presurgical planning in people with medically refractory focal epilepsy. Localization of interictal epileptiform activity, a surrogate for the seizure onset zone whose removal may prevent seizures, is challenging and depends on the use of multiple complementary techniques. Accurate and reliable localization of epileptiform activity from spontaneous MEG data has been an elusive goal. One approach toward this goal is to use a novel Bayesian inference algorithm-the Champagne algorithm with noise learning-which has shown tremendous success in source reconstruction, especially for focal brain sources. In this study, we localized sources of manually identified MEG spikes using the Champagne algorithm in a cohort of 16 patients with medically refractory epilepsy collected in two consecutive series. To evaluate the reliability of this approach, we compared the performance to equivalent current dipole (ECD) modeling, a conventional source localization technique that is commonly used in clinical practice. Results suggest that Champagne may be a robust, automated, alternative to manual parametric dipole fitting methods for localization of interictal MEG spikes, in addition to its previously described clinical and research applications.

Interhemispheric Auditory Cortical Synchronization in Asymmetric Hearing Loss

Objectives:

Auditory cortical activation of the two hemispheres to monaurally presented tonal stimuli has been shown to be asynchronous in normal hearing (NH) but synchronous in the extreme case of adult-onset asymmetric hearing loss (AHL) with single-sided deafness. We addressed the wide knowledge gap between these two anchoring states of interhemispheric temporal organization. The objectives of this study were as follows: (1) to map the trajectory of interhemispheric temporal reorganization from asynchrony to synchrony using magnitude of interaural threshold difference as the independent variable in a cross-sectional study and (2) to evaluate reversibility of interhemispheric synchrony in association with hearing in noise performance by amplifying the aidable poorer ear in a repeated measures, longitudinal study.

Design:

The cross-sectional and longitudinal cohorts were comprised of 49 subjects (AHL; N = 21; 11 male, 10 female; mean age = 48 years) and NH (N = 28; 16 male, 12 female; mean age = 45 years). The maximum interaural threshold difference of the two cohorts spanned from 0 to 65 dB. Magnetoencephalography analyses focused on latency of the M100 peak response from auditory cortex in both hemispheres between 50 msec and 150 msec following monaural tonal stimulation at the frequency (0.5, 1, 2, 3, or 4 kHz) corresponding to the maximum and minimum interaural threshold difference for better and poorer ears separately. The longitudinal AHL cohort was drawn from three subjects in the cross-sectional AHL cohort (all male; ages 49 to 60 years; varied AHL etiologies; no amplification for at least 2 years). All longitudinal study subjects were treated by monaural amplification of the poorer ear and underwent repeated measures examination of the M100 response latency and quick speech in noise hearing in noise performance at baseline, and postamplification months 3, 6, and 12.

Results:

The M100 response peak latency values in the ipsilateral hemisphere lagged those in the contralateral hemisphere for all stimulation conditions. The mean (SD) interhemispheric latency difference values (ipsilateral less contralateral) to better ear stimulation for three categories of maximum interaural threshold difference were as follows: NH (≤ 10 dB)—8.6 (3.0) msec; AHL (15 to 40 dB)—3.0 (1.2) msec; AHL (≥ 45 dB)—1.4 (1.3) msec. In turn, the magnitude of difference values were used to define interhemispheric temporal organization states of asynchrony, mixed asynchrony and synchrony, and synchrony, respectively. Amplification of the poorer ear in longitudinal subjects drove interhemispheric organization change from baseline synchrony to postamplification asynchrony and hearing in noise performance improvement in those with baseline impairment over a 12-month period.

Conclusions:

Interhemispheric temporal organization in AHL was anchored between states of asynchrony in NH and synchrony in single-sided deafness. For asymmetry magnitudes between 15 and 40 dB, the intermediate mixed state of asynchrony and synchrony was continuous and reversible. Amplification of the poorer ear in AHL improved hearing in noise performance and restored normal temporal organization of auditory cortices in the two hemispheres. The return to normal interhemispheric asynchrony from baseline synchrony and improvement in hearing following monoaural amplification of the poorer ear evolved progressively over a 12-month period.

Probabilistic Structure Learning for EEG/MEG Source Imaging With Hierarchical Graph Priors

Brain source imaging is an important method for noninvasively characterizing brain activity using Electroencephalogram (EEG) or Magnetoencephalography (MEG) recordings. Traditional EEG/MEG Source Imaging (ESI) methods usually assume the source activities at different time points are unrelated, and do not utilize the temporal structure in the source activation, making the ESI analysis sensitive to noise. Some methods may encourage very similar activation patterns across the entire time course and may be incapable of accounting the variation along the time course. To effectively deal with noise while maintaining flexibility and continuity among brain activation patterns, we propose a novel probabilistic ESI model based on a hierarchical graph prior. Under our method, a spanning tree constraint ensures that activity patterns have spatiotemporal continuity. An efficient algorithm based on an alternating convex search is presented to solve the resulting problem of the proposed model with guaranteed convergence. Comprehensive numerical studies using synthetic data on a realistic brain model are conducted under different levels of signal-to-noise ratio (SNR) from both sensor and source spaces. We also examine the EEG/MEG datasets in two real applications, in which our ESI reconstructions are neurologically plausible. All the results demonstrate significant improvements of the proposed method over benchmark methods in terms of source localization performance, especially at high noise levels.

A Novel Bayesian Approach for EEG Source Localization

We propose a new method for EEG source localization. An efficient solution to this problem requires choosing an appropriate regularization term in order to constraint the original problem. In our work, we adopt the Bayesian framework to place constraints; hence, the regularization term is closely connected to the prior distribution. More specifically, we propose a new sparse prior for the localization of EEG sources. The proposed prior distribution has sparse properties favoring focal EEG sources. In order to obtain an efficient algorithm, we use the variational Bayesian (VB) framework which provides us with a tractable iterative algorithm of closed-form equations. Additionally, we provide extensions of our method in cases where we observe group structures and spatially extended EEG sources. We have performed experiments using synthetic EEG data and real EEG data from three publicly available datasets. The real EEG data are produced due to the presentation of auditory and visual stimulus. We compare the proposed method with well-known approaches of EEG source localization and the results have shown that our method presents state-of-the-art performance, especially in cases where we expect few activated brain regions. The proposed method can effectively detect EEG sources in various circumstances. Overall, the proposed sparse prior for EEG source localization results in more accurate localization of EEG sources than state-of-the-art approaches.

Robust estimation of noise for electromagnetic brain imaging with the champagne algorithm

Robust estimation of the number, location, and activity of multiple correlated brain sources has long been a challenging task in electromagnetic brain imaging from M/EEG data, one that is significantly impacted by interference from spontaneous brain activity, sensor noise, and other sources of artifacts. Recently, we introduced the Champagne algorithm, a novel Bayesian inference algorithm that has shown tremendous success in M/EEG source reconstruction. Inherent to Champagne and most other related Bayesian reconstruction algorithms is the assumption that the noise covariance in sensor data can be estimated from “baseline” or “control” measurements. However, in many scenarios, such baseline data is not available, or is unreliable, and it is unclear how best to estimate the noise covariance. In this technical note, we propose several robust methods to estimate the contributions to sensors from noise arising from outside the brain without the need for additional baseline measurements. The incorporation of these methods for diagonal noise covariance estimation improves the robust reconstruction of complex brain source activity under high levels of noise and interference, while maintaining the performance features of Champagne. Specifically, we show that the resulting algorithm, Champagne with noise learning, is quite robust to initialization and is computationally efficient. In simulations, performance of the proposed noise learning algorithm is consistently superior to Champagne without noise learning. We also demonstrate that, even without the use of any baseline data, Champagne with noise learning is able to reconstruct complex brain activity with just a few trials or even a single trial, demonstrating significant improvements in source reconstruction for electromagnetic brain imaging.

Tinnitus Neuroimaging

This article reviews the use of human neuroimaging for chronic subjective tinnitus. Evidence-based guidance on the clinical use of imaging to identify relevant auditory lesions when evaluating tinnitus patients is given. After introducing the anatomy and imaging modalities most pertinent to the neuroscience of tinnitus, the article reviews tinnitus-associated alterations in key auditory and nonauditory networks in the central nervous system. Emphasis is placed on how these findings support proposed models of tinnitus and how this line of investigation is relevant to practicing clinicians.