Biomedical subjects

Arye Nehorai

Publications and source records attributed to Arye Nehorai.

10 recordsLinked to original sources

PubMed2006

Performance analysis of reduced-rank beamformers for estimating dipole source signals using EEG/MEG.

We study the performance of various beamformers for estimating a current dipole source at a known location using electroencephalography (EEG) and magnetoencephalography (MEG). We present our beamformers in the form of the generalized sidelobe canceler (GSC). Under this structure, the beamformer can be solved by finding a filter that achieves the minimum mean-squared error (MMSE) between the mainbeam response and filtered observed signal. We express the MMSE as a function of the filter's rank and use it as a criterion to evaluate the performance of the beamformers. We do not make any assumptions on the rank of the interference-plus-noise covariance matrix. Instead, we treat it as low-rank and derive a general expression for the MMSE. We present numerical examples to compare the MSE performance of beamformers commonly studied in the literature: principal components (PCs), cross-spectral metrics (CSMs), and eigencanceler (EIG) beamformers. Our results show that good estimates of the dipole source signals can be achieved using reduced-rank beamformers even for low signal-to-noise ratio (SNR) values.

Action Potentials↗

PubMed2006

Effects of geometric head model perturbations on the EEG forward and inverse problems.

We study the effect of geometric head model perturbations on the electroencephalography (EEG) forward and inverse problems. Small magnitude perturbations of the shape of the head could represent uncertainties in the head model due to errors on images or techniques used to construct the model. They could also represent small scale details of the shape of the surfaces not described in a deterministic model, such as the sulci and fissures of the cortical layer. We perform a first-order perturbation analysis, using a meshless method for computing the sensitivity of the solution of the forward problem to the geometry of the head model. The effect on the forward problem solution is treated as noise in the EEG measurements and the Cramér-Rao bound is computed to quantify the effect on the inverse problem performance. Our results show that, for a dipolar source, the effect of the perturbations on the inverse problem performance is under the level of the uncertainties due to the spontaneous brain activity. Thus, the results suggest that an extremely detailed model of the head may be unnecessary when solving the EEG inverse problem.

Brain↗

PubMed2006

Surface-source modeling and estimation using biomagnetic measurements.

We propose a number of electric source models that are spatially distributed on an unknown surface for biomagnetism. These can be useful to model, e.g., patches of electrical activity on the cortex. We use a realistic head (or another organ) model and discuss the special case of a spherical head model with radial sensors resulting in more efficient computations of the estimates for magnetoencephalography. We derive forward solutions, maximum likelihood (ML) estimates, and Cramér-Rao bound (CRB) expressions for the unknown source parameters. A model selection method is applied to decide on the most appropriate model. We also present numerical examples to compare the performances and computational costs of the different models and illustrate when it is possible to distinguish between surface and focal sources or line sources. Finally, we apply our methods to real biomagnetic data of phantom human torso and demonstrate the applicability of them.

Action Potentials↗

PubMed2005

Ellipsoidal head model for fetal magnetoencephalography: forward and inverse solutions.

Fetal magnetoencephalography (fMEG) is a non-invasive technique where measurements of the magnetic field outside the maternal abdomen are used to infer the source location and signals of the fetus' neural activity. There are a number of aspects related to fMEG modelling that must be addressed, such as the conductor volume, fetal position and orientation, gestation period, etc. We propose a solution to the forward problem of fMEG based on an ellipsoidal head geometry. This model has the advantage of highlighting special characteristics of the field that are inherent to the anisotropy of the human head, such as the spread and orientation of the field in relationship with the localization and position of the fetal head. Our forward solution is presented in the form of a kernel matrix that facilitates the solution of the inverse problem through decoupling of the dipole localization parameters from the source signals. Then, we use this model and the maximum likelihood technique to solve the inverse problem assuming the availability of measurements from multiple trials. The applicability and performance of our methods are illustrated through numerical examples based on a real 151-channel SQUID fMEG measurement system (SARA). SARA is an MEG system especially designed for fetal assessment and is currently used for heart and brain studies. Finally, since our model requires knowledge of the best-fitting ellipsoid's centre location and semiaxes lengths, we propose a method for estimating these parameters through a least-squares fit on anatomical information obtained from three-dimensional ultrasound images.

Algorithms↗

PubMed2005

A meshless method for solving the EEG forward problem.

We present a numerical method to solve the quasistatic Maxwell equations and compute the electroencephalography (EEG) forward problem solution. More generally, we develop a computationally efficient method to obtain the electric potential distribution generated by a source of electric activity inside a three-dimensional body of arbitrary shape and layers of different electric conductivities. The method needs only a set of nodes on the surface and inside the head, but not a mesh connecting the nodes. This represents an advantage over traditional methods like boundary elements or finite elements since the generation of the mesh is typically computationally intensive. The performance of the proposed method is compared with the boundary element method (BEM) by numerically solving some EEG forward problems examples. For a large number of nodes and the same precision, our method has lower computational load than BEM due to a faster convergence rate and to the sparsity of the linear system to be solved.

Algorithms↗

PubMed2005

Model selection in spatio-temporal electromagnetic source analysis.

Several methods [model selection procedures (MSPs)] to determine the number of sources in electroencephalogram (EEG) and magnetoencphalogram (MEG) data have previously been investigated in an instantaneous analysis. In this paper, these MSPs are extended to a spatio-temporal analysis if possible. It is seen that the residual variance (RV) tends to overestimate the number of sources. The Akaike information criterion (AIC) and the Wald test on amplitudes (WA) and the Wald test on locations (WL) have the highest probabilities of selecting the correct number of sources. The WA has the advantage that it offers the opportunity to test which source is active at which time sample.

Algorithms↗

PubMed2005

Distinguishing between moving and stationary sources using EEG/MEG measurements with an application to epilepsy.

Performances of electroencephalography (EEG) and magnetoencephalography (MEG) source estimation methods depend on the validity of the assumed model. In many cases, the model structure is related to physical information. We discuss a number of statistical selection methods to distinguish between two possible models using least-squares estimation and assuming a spherical head model. The first model has a single moving source whereas the second has two stationary sources; these may result in similar EEG/MEG measurements. The need to decide between such models occurs for example in Jacksonian seizures (e.g., epilepsy) or in intralobular activities, where a model with either two stationary dipole sources or a single moving dipole source may be possible. We also show that all of the selection methods discussed choose the correct model with probability one when the number of trials goes to infinity. Finally we present numerical examples and compare the performances of the methods by varying parameters such as the signal-to-noise ratio, source depth, and separation of sources, and also apply the methods to real MEG data for epilepsy.

Algorithms↗

PubMed2005

Line-source modeling and estimation with magnetoencephalography.

We propose a number of source models that are spatially distributed on a line for magnetoencephalography (MEG) using both a spherical head with radial sensors for more efficient computation and a realistic head model for more accurate results. We develop these models with increasing degrees of freedom, derive forward solutions, maximum-likelihood (ML) estimates, and Cramér-Rao bound (CRB) expressions for the unknown source parameters. A model selection method is applied to select the most appropriate model. We also present numerical examples to compare the performances and computational costs of the different models, to determine the regions where better estimates are possible and when it is possible to distinguish between line and focal sources. We demonstrate the usefulness of the proposed line-source models over the previously available focal source model in certain distributed source cases. Finally, we apply our methods to real MEG data, the N2O response after electric stimulation of the median nerve known to be an extended source.

Algorithms↗

PubMed2005

MEG forward problem formulation using equivalent surface current densities.

We present a formulation for the magnetoencephalography (MEG) forward problem with a layered head model. Traditionally the magnetic field is computed based on the electric potential on the interfaces between the layers. We propose to express the effect of the volumetric currents in terms of an equivalent surface current density on each interface, and obtain the magnetic field based on them. The boundary elements method is used to compute the equivalent current density and the magnetic field for a realistic head geometry. We present numerical results showing that the MEG forward problem is solved correctly with this formulation, and compare it with the performance of the traditional formulation. We conclude that the traditional formulation generally performs better, but still the new formulation is useful in certain situations.

Animals↗

PubMed2004

Estimating brain conductivities and dipole source signals with EEG arrays.

Techniques based on electroencephalography (EEG) measure the electric potentials on the scalp and process them to infer the location, distribution, and intensity of underlying neural activity. Accuracy in estimating these parameters is highly sensitive to uncertainty in the conductivities of the head tissues. Furthermore, dissimilarities among individuals are ignored when standarized values are used. In this paper, we apply the maximum-likelihood and maximum a posteriori (MAP) techniques to simultaneously estimate the layer conductivity ratios and source signal using EEG data. We use the classical 4-sphere model to approximate the head geometry, and assume a known dipole source position. The accuracy of our estimates is evaluated by comparing their standard deviations with the Cramér-Rao bound (CRB). The applicability of these techniques is illustrated with numerical examples on simulated EEG data. Our results show that the estimates have low bias and attain the CRB for sufficiently large number of experiments. We also present numerical examples evaluating the sensitivity to imprecise assumptions on the source position and skull thickness. Finally, we propose extensions to the case of unknown source position and present examples for real data.

Brain↗