Search

We study the effect of the head shape variations on the EEG/magnetoencephalography (MEG) forward and inverse problems. We build a random head model such that each sample represents the head shape of a different individual and solve the forward problem assuming this random head model, using a polynomial chaos expansion. The random solution of the forward problem is then used to quantify the effect of the geometry when the inverse problem is solved with a standard head model. The results derived with this approach are valid for a continuous family of head models, rather than just for a set of cases. The random model consists of three random surfaces that define layers of different electric conductivity, and we built an example based on a set of 30 deterministic models from adults. Our results show that for a dipolar source model, the effect of the head shape variations on the EEG/MEG inverse problem due to the random head model is slightly larger than the effect of the electronic noise present in the sensors. The variations in the EEG inverse problem solutions are due to the variations in the shape of the volume conductor, while the variations in the MEG inverse problem solutions, larger than the EEG ones, are caused mainly by the variations of the absolute position of the sources in a coordinate system based on anatomical landmarks, in which the magnetometers have a fixed position. Index Terms--EEG/magnetoencephalography (MEG) average head model, polynomial chaos expansion (PCE), sparse grids, stochastic modeling.

Electroencephalography (EEG) is an important tool for studying the brain functions and is becoming popular in clinical practice. In this paper, we develop four parametric EEG models to estimate current sources that are spatially distributed on a surface. Our models approximate the source shape and extent explicitly and can be applied to localize extended sources which are often encountered, e.g., in epilepsy diagnosis. We assume a realistic head model and solve the EEG forward problem using the boundary element method. We present the source models with increasing degrees of freedom, provide the forward solutions, and derive the maximum-likelihood estimates as well as Cramer-Rao bounds of the unknown source parameters. In order to evaluate the applicability of the proposed models, we first compare their estimation performances with the dipole model's using several known source distributions. We then discuss the conditions under which we can distinguish between the proposed extended sources and the focal dipole using the generalized likelihood ratio test. We also apply our models to the electric measurements obtained from a phantom body in which an extended electric source is imbedded. We observe that the proposed model can capture the source extent information satisfactorily and the localization accuracy is better than the dipole model. Index Terms--Cramer-Rao bounds, EEG, extended source modeling, likelihood ratio test.

We develop three parametric models for electroencephalography (EEG) to estimate current sources that are spatially distributed on a line. We assume a realistic head model and solve the EEG forward problem using the boundary element method (BEM). We present the models with increasing degrees of freedom, provide the forward solutions, and derive the maximum-likelihood estimates as well as Cramer-Rao bounds of the unknown source parameters. A series of experiments are conducted to evaluate the applicability of the proposed models. We use numerical examples to demonstrate the usefulness of our line-source models in estimating extended sources. We also apply our models to the real EEG data of N20 response that is known to have an extended source. We observe that the line-source models explain the N20 measurements better than the dipole model. Index Terms--Cramer-Rao bounds, EEG, extended source modeling.

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 Cramer-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. Index Terms--Biomagnetic measurements, magnetoencephalograph, source localization, spatially extended sources, surface-source models.

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 Cramer-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. Index Terms--Cramer-Rao bound, EEG, inverse problem, perturbations analysis, stochastic modeling.

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. Index Terms--BEM, current density, forward problem, MEG.

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 Cramer-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 N20 response after electric stimulation of the median nerve known to be an extended source. Index Terms--Line-source models, magnetoencephalography, N20 responses, source localization.

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. Index Terms--Electroencephalography, epilepsy, magnetoencephalography, model selection.

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. Index Terms--EEG, EEG forward problem, layered media, meshless method, moving least squares approximation, numerical solution, volume conductor.

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 Cramer-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. Index Terms--Brain conductivities, Cramer-Rao bound, electroencephalography, maximum-likelihood estimation, parameter estimation, sensor array processing.

A clamped-mode resonant converter (CMRC) is proposed to be used as a transmitter which broadcasts correction signals in a differential Global Positioning System (DGPS). The digital information of the DGPS modulates the carrier with different modulation methods minimum shift keying (MSK) (recommended for the marine radiobeacon band), stepped binary phase shift keying (S-BPSK) and BPSK. The first two methods allow keeping the converter's operation in soft switching mode. Extensive trials under impedance mismatch and modulation conditions confirm that the soft switching mode is preserved.

An innovative frequency-domain approach to multi-input, multi-output (MIMO) linear filter design for sample-data systems is described and discussed.

A mean-square error lower bound for the discrete-time nonlinear filtering problem was obtained based on the Van Trees version of the Cramer-Rao inequality. The case of singular conditional distribution of the one-step-ahead state vector given the present state was examined. The bound was analyzed for three examples: the recursive estimation of slowly changing parameters of an autoregressive process, tracking a slowly varying frequency of a single cisoid in noise, and tracking parameters of a sinusoidal frequency with sinusoidal phase modulation.

In this paper, a design procedure for clamped-mode resonant converters working with variable load is proposed and analyzed. The operation of these converters with general parallel or series loads is reviewed first. The operation of the transistors as zero-voltage switches is identified and characterized as mode A. Mode A is preferred to simplify the implementation of the power switches and to increment the converter efficiency. As a case study, the design of an RF power amplifier is presented. The amplifier is connected to a variable load with reactive characteristics, through a low-pass filter of three elements. The validity of the proposed design approach is verified with an experimental setup. Index Terms--Power amplifier, resonant power conversion.