822 research outputs found
A two-way regularization method for MEG source reconstruction
The MEG inverse problem refers to the reconstruction of the neural activity
of the brain from magnetoencephalography (MEG) measurements. We propose a
two-way regularization (TWR) method to solve the MEG inverse problem under the
assumptions that only a small number of locations in space are responsible for
the measured signals (focality), and each source time course is smooth in time
(smoothness). The focality and smoothness of the reconstructed signals are
ensured respectively by imposing a sparsity-inducing penalty and a roughness
penalty in the data fitting criterion. A two-stage algorithm is developed for
fast computation, where a raw estimate of the source time course is obtained in
the first stage and then refined in the second stage by the two-way
regularization. The proposed method is shown to be effective on both synthetic
and real-world examples.Comment: Published in at http://dx.doi.org/10.1214/11-AOAS531 the Annals of
Applied Statistics (http://www.imstat.org/aoas/) by the Institute of
Mathematical Statistics (http://www.imstat.org
Inverse Modeling for MEG/EEG data
We provide an overview of the state-of-the-art for mathematical methods that
are used to reconstruct brain activity from neurophysiological data. After a
brief introduction on the mathematics of the forward problem, we discuss
standard and recently proposed regularization methods, as well as Monte Carlo
techniques for Bayesian inference. We classify the inverse methods based on the
underlying source model, and discuss advantages and disadvantages. Finally we
describe an application to the pre-surgical evaluation of epileptic patients.Comment: 15 pages, 1 figur
State-space solutions to the dynamic magnetoencephalography inverse problem using high performance computing
Determining the magnitude and location of neural sources within the brain
that are responsible for generating magnetoencephalography (MEG) signals
measured on the surface of the head is a challenging problem in functional
neuroimaging. The number of potential sources within the brain exceeds by an
order of magnitude the number of recording sites. As a consequence, the
estimates for the magnitude and location of the neural sources will be
ill-conditioned because of the underdetermined nature of the problem. One
well-known technique designed to address this imbalance is the minimum norm
estimator (MNE). This approach imposes an regularization constraint that
serves to stabilize and condition the source parameter estimates. However,
these classes of regularizer are static in time and do not consider the
temporal constraints inherent to the biophysics of the MEG experiment. In this
paper we propose a dynamic state-space model that accounts for both spatial and
temporal correlations within and across candidate intracortical sources. In our
model, the observation model is derived from the steady-state solution to
Maxwell's equations while the latent model representing neural dynamics is
given by a random walk process.Comment: Published in at http://dx.doi.org/10.1214/11-AOAS483 the Annals of
Applied Statistics (http://www.imstat.org/aoas/) by the Institute of
Mathematical Statistics (http://www.imstat.org
Frontal and superior temporal auditory processing abnormalities in schizophrenia
AbstractBackgroundAlthough magnetoencephalography (MEG) studies show superior temporal gyrus (STG) auditory processing abnormalities in schizophrenia at 50 and 100ms, EEG and corticography studies suggest involvement of additional brain areas (e.g., frontal areas) during this interval. Study goals were to identify 30 to 130ms auditory encoding processes in schizophrenia (SZ) and healthy controls (HC) and group differences throughout the cortex.MethodsThe standard paired-click task was administered to 19 SZ and 21 HC subjects during MEG recording. Vector-based Spatial–temporal Analysis using L1-minimum-norm (VESTAL) provided 4D maps of activity from 30 to 130ms. Within-group t-tests compared post-stimulus 50ms and 100ms activity to baseline. Between-group t-tests examined 50 and 100ms group differences.ResultsBilateral 50 and 100ms STG activity was observed in both groups. HC had stronger bilateral 50 and 100ms STG activity than SZ. In addition to the STG group difference, non-STG activity was also observed in both groups. For example, whereas HC had stronger left and right inferior frontal gyrus activity than SZ, SZ had stronger right superior frontal gyrus and left supramarginal gyrus activity than HC.ConclusionsLess STG activity was observed in SZ than HC, indicating encoding problems in SZ. Yet auditory encoding abnormalities are not specific to STG, as group differences were observed in frontal and SMG areas. Thus, present findings indicate that individuals with SZ show abnormalities in multiple nodes of a concurrently activated auditory network
Fast Optimal Transport Averaging of Neuroimaging Data
Knowing how the Human brain is anatomically and functionally organized at the
level of a group of healthy individuals or patients is the primary goal of
neuroimaging research. Yet computing an average of brain imaging data defined
over a voxel grid or a triangulation remains a challenge. Data are large, the
geometry of the brain is complex and the between subjects variability leads to
spatially or temporally non-overlapping effects of interest. To address the
problem of variability, data are commonly smoothed before group linear
averaging. In this work we build on ideas originally introduced by Kantorovich
to propose a new algorithm that can average efficiently non-normalized data
defined over arbitrary discrete domains using transportation metrics. We show
how Kantorovich means can be linked to Wasserstein barycenters in order to take
advantage of an entropic smoothing approach. It leads to a smooth convex
optimization problem and an algorithm with strong convergence guarantees. We
illustrate the versatility of this tool and its empirical behavior on
functional neuroimaging data, functional MRI and magnetoencephalography (MEG)
source estimates, defined on voxel grids and triangulations of the folded
cortical surface.Comment: Information Processing in Medical Imaging (IPMI), Jun 2015, Isle of
Skye, United Kingdom. Springer, 201
Applying neural networks for improving the MEG inverse solution
Magnetoencephalography (MEG) and electroencephalography (EEG) are appealing non-invasive methods for recording brain activity with high temporal resolution. However, locating the brain source currents from recordings picked up by the sensors on the scalp introduces an ill-posed inverse problem. The MEG inverse problem one of the most difficult inverse problems in medical imaging. The current standard in approximating the MEG inverse problem is to use multiple distributed inverse solutions – namely dSPM, sLORETA and L2 MNE – to estimate the source current distribution in the brain. This thesis investigates if these inverse solutions can be "post-processed" by a neural network to provide improved accuracy on source locations.
Recently, deep neural networks have been used to approximate other ill-posed inverse medical imaging problems with accuracy comparable to current state-of- the-art inverse reconstruction algorithms. Neural networks are powerful tools for approximating problems with limited prior knowledge or problems that require high levels of abstraction. In this thesis a special case of a deep convolutional network, the U-Net, is applied to approximate the MEG inverse problem using the standard inverse solutions (dSPM, sLORETA and L2 MNE) as inputs.
The U-Net is capable of learning non-linear relationships between the inputs and producing predictions about the site of single-dipole activation with higher accuracy than the L2 minimum-norm based inverse solutions with the following resolution metrics: dipole localization error (DLE), spatial dispersion (SD) and overall amplitude (OA). The U-Net model is stable and performs better in aforesaid resolution metrics than the inverse solutions with multi-dipole data previously unseen by the U-Net
- …