A systematic evaluation of the spherical model accuracy in EEG dipole localization

High-resolution EEG source localization in personalized segmentation-free head model with multi-dipole fitting

Objective. Electroencephalograms (EEGs) are often used to monitor brain activity. Several source localization methods have been proposed to estimate the location of brain activity corresponding to EEG readings. However, only a few studies evaluated source localization accuracy from measured EEG using personalized head models in a millimeter resolution. In this study, based on a volume conductor analysis of a high-resolution personalized human head model constructed from magnetic resonance images, a finite difference method was used to solve the forward problem and to reconstruct the field distribution.Approach. We used a personalized segmentation-free head model developed using machine learning techniques, in which the abrupt change of electrical conductivity occurred at the tissue interface is suppressed. Using this model, a smooth field distribution was obtained to address the forward problem. Next, multi-dipole fitting was conducted using EEG measurements for each subject (N= 10 male subjects, age: 22.5 ± 0.5), and the source location and electric field distribution were estimated.Main results.For measured somatosensory evoked potential for electrostimulation to the wrist, a multi-dipole model with lead field matrix computed with the volume conductor model was found to be superior than a single dipole model when using personalized segmentation-free models (6/10). The correlation coefficient between measured and estimated scalp potentials was 0.89 for segmentation-free head models and 0.71 for conventional segmented models. The proposed method is straightforward model development and comparable localization difference of the maximum electric field from the target wrist reported using fMR (i.e. 16.4 ± 5.2 mm) in previous study. For comparison, DUNEuro based on sLORETA was (EEG: 17.0 ± 4.0 mm). In addition, somatosensory evoked magnetic fields obtained by Magnetoencephalography was 25.3 ± 8.5 mm using three-layer sphere and sLORETA.Significance. For measured EEG signals, our procedures using personalized head models demonstrated that effective localization of the somatosensory cortex, which is located in a non-shallower cortex region. This method may be potentially applied for imaging brain activity located in other non-shallow regions.

Technical Considerations in EEG Source Imaging

SUMMARY

EEG source imaging is an established technique for identifying the origin of interictal and ictal epileptiform discharges in patients with epilepsy, and it is an important tool in neurophysiology research. Accurate and reliable EEG source imaging requires appropriate choices of how the head, skull, and scalp are modeled, and understanding of the different approaches to modeling is important to guide these choices. Similarly, numerous different approaches to modeling the electrical sources within the brain exist, and appropriate understanding of the strengths and limitations of each are essential to obtaining accurate, reliable, and interpretable solutions. This review aims to describe the essential theoretical basis for these head and source models while also discussing the practical implications of each in clinical or research applications.

From theory to practical fundamentals of electroencephalographic source imaging in localizing the epileptogenic zone

With continued advancement in computational technologies, the analysis of electroencephalography (EEG) has shifted from pure visual analysis to a noninvasive computational technique called EEG source imaging (ESI), which involves mathematical modeling of dipolar and distributed sources of a given scalp EEG pattern. ESI is a noninvasive phase I test for presurgical localization of the seizure onset zone in focal epilepsy. It is a relatively inexpensive modality, as it leverages scalp EEG and magnetic resonance imaging (MRI) data already collected typically during presurgical evaluation. With an adequate number of electrodes and combined with patient-specific MRI-based head models, ESI has proven to be a valuable and accurate clinical diagnostic tool for localizing the epileptogenic zone. Despite its advantages, however, ESI is routinely used at only a minority of epilepsy centers. This paper reviews the current evidence and practical fundamentals for using ESI of interictal and ictal epileptic activity during the presurgical evaluation of drug-resistant patients. We identify common errors in processing and interpreting ESI studies, describe the differences in approach needed for localizing interictal and ictal EEG discharges through practical examples, and describe best practices for optimizing the diagnostic information available from these studies.

Equivalent current dipole sources of neurofeedback training-induced alpha activity through temporal/spectral analytic techniques

Much of the work in alpha NFT has focused on evaluating changes in alpha amplitude. However, the generation mechanism of training-induced alpha activity has not yet been clarified. The present study aimed to identify sources of training-induced alpha activity through four temporal/spectral analytic techniques, i.e., the max peak average (MPA), positive average (PA), negative average (NA) and event-related spectral perturbation average (ERSPA) methods. Thirty-five healthy participants were recruited into an alpha group receiving feedback of 8–12-Hz amplitudes, and twenty-eight healthy participants were recruited into a control group receiving feedback of random 4-Hz amplitudes from the range of 7 to 20 Hz. Twelve sessions were performed within 4 weeks (3 sessions per week). The control group had no change in the amplitude spectrum. In contrast, twenty-nine participants in the alpha group showed significant alpha amplitude increases exclusively and were identified as “responders”. A whole-head EEG was recorded for the “responders” after NFT. The epochs of training-induced alpha activity from whole-head EEG were averaged by four different methods for equivalent current dipole source analysis. High agreement and Cohen’s kappa coefficients on dipole source localization between each method were observed, showing that the dipole clusters of training-induced alpha activity were consistently located in the precuneus, posterior cingulate cortex (PCC) and middle temporal gyrus. The residual variance (goodness of fit) for dipole estimation of the MPA was significantly smaller than that of the others. Our findings indicate that the precuneus, PCC and middle temporal gyrus play important roles in enhancing training-induced alpha activity. The four averaging methods (especially the MPA method) were suitable for investigating sources of brainwaves. Additionally, three dipoles can be used for dipole source analysis of training-induced alpha activity in future research, especially the training sites are around the central regions.

Practical Fundamentals of Clinical MEG Interpretation in Epilepsy

Magnetoencephalography (MEG) is a neurophysiologic test that offers a functional localization of epileptic sources in patients considered for epilepsy surgery. The understanding of clinical MEG concepts, and the interpretation of these clinical studies, are very involving processes that demand both clinical and procedural expertise. One of the major obstacles in acquiring necessary proficiency is the scarcity of fundamental clinical literature. To fill this knowledge gap, this review aims to explain the basic practical concepts of clinical MEG relevant to epilepsy with an emphasis on single equivalent dipole (sECD), which is one the most clinically validated and ubiquitously used source localization method, and illustrate and explain the regional topology and source dynamics relevant for clinical interpretation of MEG-EEG.

Posterior insular activity contributes to the late laser-evoked potential component in EEG recordings

OBJECTIVE

Nociceptive activity in some brain areas has concordantly been reported in EEG source models, such as the anterior/mid-cingulate cortex and the parasylvian area. Whereas the posterior insula has been constantly reported to be active in intracortical and fMRI studies, non-invasive EEG and MEG recordings mostly failed to detect activity in this region. This study aimed to determine an appropriate inverse modeling approach in EEG recordings to model posterior insular activity, assuming the late LEP (laser evoked potential) time window to yield a better separation from other ongoing cortical activity.

METHODS

In 12 healthy volunteers, nociceptive stimuli of three intensities were applied. LEP were recorded using 32-channel EEG recordings. Source analysis was performed in specific time windows defined in the grand-average dataset. Two distinct dipole-pairs located close to the operculo-insular area were compared.

RESULTS

Our results show that posterior insular activity yields a substantial contribution to the latest part (positive component) of the LEP.

CONCLUSIONS

Even though the initial insular activity onset is in the early LEP time window,modelingthe insular activity in the late LEP time window might result in better separation from other ongoing cortical activity.

SIGNIFICANCE

Modeling the late LEP activity might enable to distinguish posterior insular activity.

Effects of head models and dipole source parameters on EEG fields

Head model and an efficient method for computing the forward EEG (electroencephalography)problem are essential to dipole source localization(DSL). In this paper, we use less expensive ovoid geometry to approximate human head, aiming at investigating the effects of head shape and dipole source parameters on EEG fields. The application of point least squares (PLS) based on meshless method was introduced for solving EEG forward problem and numerical simulation is implemented in three kinds of ovoid head models. We present the performances of the surface potential in the face of varying dipole source parameters in detail. The results show that the potential patterns are similar for different dipole position in different head shapes, but the peak value of potential is significantly influenced by the head shape. Dipole position induces a great effect on the peak value of potential and shift of peak potential. The degree of variation between sphere head model and non-sphere head models is seen at the same time. We also show that PLS method with the trigonometric basis is superior to the constant basis, linear basis, and quadratic basis functions in accuracy and efficiency.

Relationship of EEG sources of neonatal seizures to acute perinatal brain lesions seen on MRI: a pilot study

Even though it is known that neonatal seizures are associated with acute brain lesions, the relationship of electroencephalographic (EEG) seizures to acute perinatal brain lesions visible on magnetic resonance imaging (MRI) has not been objectively studied. EEG source localization is successfully used for this purpose in adults, but it has not been sufficiently explored in neonates. Therefore, we developed an integrated method for ictal EEG dipole source localization based on a realistic head model to investigate the utility of EEG source imaging in neonates with postasphyxial seizures. We describe here our method and compare the dipole seizure localization results with acute perinatal lesions seen on brain MRI in 10 full-term infants with neonatal encephalopathy. Through experimental studies, we also explore the sensitivity of our method to the electrode positioning errors and the variations in neonatal skull geometry and conductivity. The localization results of 45 focal seizures from 10 neonates are compared with the visual analysis of EEG and MRI data, scored by expert physicians. In 9 of 10 neonates, dipole locations showed good relationship with MRI lesions and clinical data. Our experimental results also suggest that the variations in the used values for skull conductivity or thickness have little effect on the dipole localization, whereas inaccurate electrode positioning can reduce the accuracy of source estimates. The performance of our fused method indicates that ictal EEG source imaging is feasible in neonates and with further validation studies, this technique can become a useful diagnostic tool.

Effects of forward model errors on EEG source localization

Subject-specific four-layer boundary element method (BEM) electrical forward head models for four participants, generated from magnetic resonance (MR) head images using NFT ( www.sccn.ucsd.edu/wiki/NFT ), were used to simulate electroencephalographic (EEG) scalp potentials at 256 recorded electrode positions produced by single current dipoles of a 3-D grid in brain space. Locations of these dipoles were then estimated using gradient descent within five template head models fit to the electrode positions. These were: a spherical model, three-layer and four-layer BEM head models based on the Montreal Neurological Institute (MNI) template head image, and these BEM models warped to the recorded electrode positions. Smallest localization errors (4.1-6.2 mm, medians) were obtained using the electrode-position warped four-layer BEM models, with largest localization errors (~20 mm) for most basal brain locations. When we increased the brain-to-skull conductivity ratio assumed in the template model scalp projections from the simulated value (25:1) to a higher value (80:1) used in earlier studies, the estimated dipole locations moved outwards (12.4 mm, median). We also investigated the effects of errors in co-registering the electrode positions, of reducing electrode counts, and of adding a fifth, isotropic white matter layer to one individual head model. Results show that when individual subject MR head images are not available to construct subject-specific head models, accurate EEG source localization should employ a four- or five-layer BEM template head model incorporating an accurate skull conductivity estimate and warped to 64 or more accurately 3-D measured and co-registered electrode positions.

A realistic volume conductor model of the neonatal head: methods, challenges and applications

Developing a realistic volume conductor head model is an important step towards a non-invasive investigation of neuro-electrical activity in the brain. For adults, different volume conductor head models have been designed and successfully used for electroencephalography (EEG) source analysis. However, creating appropriate neonatal volume conductor head model for EEG source analysis is a challenging task mainly due to insufficient knowledge of head tissue conductivities and complex anatomy of the developing newborn brain. In this work, we present a pipeline for modeling a realistic volume conductor model of the neonatal head, where we address the modeling challenges and propose our solutions. We also discuss the use of our realistic volume conductor head model for neonatal EEG source analysis.

Cortical oscillatory dynamics in a social interaction model

In this study we sought to investigate cortical oscillatory dynamics accompanying three major kinds of social behavior: aggressive, friendly, and avoidant. Behavioral and EEG data were collected in 48 participants during a computer game modeling social interactions with virtual 'persons'. 3D source reconstruction and independent component analysis were applied to EEG data. Results showed that social behavior was partly reactive and partly proactive with subject's personality playing an important role in shaping this behavior. Most salient differences were found between avoidance and approach behaviors, whereas the two kinds of approach behavior (i.e., aggression and friendship) did not differ from each other. Comparative to avoidance, approach behaviors were associated with higher induced responses in most frequency bands which were mostly observed in cortical areas overlapping with the default mode network. The difference between approach- and avoidance-related oscillatory dynamics was more salient in subjects predisposed to approach behaviors (i.e., in aggressive or sociable subjects) and was less pronounced in subjects predisposed to avoidance behavior (i.e., in high trait anxiety scorers). There was a trend to higher low frequency phase-locking in motor area in approach than in avoid condition. Results are discussed in light of the concept linking induced responses with top-down and evoked responses with bottom-up processes.

The default mode network and EEG α oscillations: an independent component analysis

The default mode network (DMN) has been principally investigated using positron emission tomography (PET) and functional magnetic resonance imaging (fMRI) and has received mixed support in electroencephalographic (EEG) studies. In particular, the existing evidence is too inconsistent to allow formulation of specific hypotheses linking DMN activity to traditional EEG frequency bands. In this study, we aimed to test whether blind decomposition methods are able to identify in EEG data spatial patterns resembling the DMN as it is described in PET and fMRI studies. Further we aimed to test a degree of task-relatedness of DMN patterns identified in the traditional EEG frequency bands. To answer these questions we collected data both in a resting state and during performance of two experimental tasks: an explicit judgment of facial affect and a social game task. Individual differences in amount of self-referential thoughts during the resting state were measured by a short self-report scale. Only alpha band spatial patterns simultaneously showed a considerable overlap with the DMN and high correlations with presumptive DMN function-related outcomes both in the resting state and during the social game task. Spontaneous self-referential thoughts were associated with enhanced alpha activity in the posterior DMN hub, whereas processing of DMN function-related external stimuli disrupted this activity and simultaneously caused partial alpha phase-locking to external events. This evidence implies that synchronization of internal mental processes, as opposed to the processing of external stimuli, might be the primary function of alpha oscillations which is bound to be related to activity of the DMN.

Cortical representation of the human hand assessed by two levels of high-resolution EEG recordings

Increasing interest in cortical plasticity has prompted the growing use of somatosensory evoked potentials (SEPs) to estimate changes in the cortical representation of body regions. Here, we tested the effect of different sites of hand stimulation and of the density of spatial sampling in the quality of estimation of somatosensory sources. Sources of two SEP components from the primary somatosensory cortex (N20/P20 and P45) were estimated using two levels of spatial sampling (64- vs. 128-channel) and stimulation of four distal sites in the upper limbs, including single digits (first vs. fifth) and distal nerves with comparable cortical projection (superficial branch of the radial nerve and distal ulnar nerve). The most robust separation of somatosensory sources was achieved by comparing the cortical representations of the first digit and the distal ulnar nerve territories on the N20/P20 component of SEPs. Although both the 64- and the 128-electrode montages correctly discriminated these two areas, only the 128-electrode montage was able to significantly separate sources in the other cases, notably when using first versus fifth digit stimulation. Trustworthy distinction of cortical representations was not obtainable when using the P45 component, probably because of greater activation volume, radial orientation of sources in areas 1-2 and increased variability with attention and vigilance. Assessment of tangential SEP components to stimulation of first digit versus ulnar nerve appears the best option to assess plastic somatosensory changes, especially when using relatively low-electrode sampling.

Utility of independent component analysis for interpretation of intracranial EEG

Electrode arrays are sometimes implanted in the brains of patients with intractable epilepsy to better localize seizure foci before epilepsy surgery. Analysis of intracranial EEG (iEEG) recordings is typically performed in the electrode channel domain without explicit separation of the sources that generate the signals. However, intracranial EEG signals, like scalp EEG signals, could be linear mixtures of local activity and volume-conducted activity arising in multiple source areas. Independent component analysis (ICA) has recently been applied to scalp EEG data, and shown to separate the signal mixtures into independently generated brain and non-brain source signals. Here, we applied ICA to unmix source signals from intracranial EEG recordings from four epilepsy patients during a visually cued finger movement task in the presence of background pathological brain activity. This ICA decomposition demonstrated that the iEEG recordings were not maximally independent, but rather are linear mixtures of activity from multiple sources. Many of the independent component (IC) projections to the iEEG recording grid were consistent with sources from single brain regions, including components exhibiting classic movement-related dynamics. Notably, the largest IC projection to each channel accounted for no more than 20–80% of the channel signal variance, implying that in general intracranial recordings cannot be accurately interpreted as recordings of independent brain sources. These results suggest that ICA can be used to identify and monitor major field sources of local and distributed functional networks generating iEEG data. ICA decomposition methods are useful for improving the fidelity of source signals of interest, likely including distinguishing the sources of pathological brain activity.

Dipole source analyses of early median nerve SEP components obtained from subdural grid recordings

The median nerve N20 and P22 SEP components constitute the initial response of the primary somatosensory cortex to somatosensory stimulation of the upper extremity. Knowledge of the underlying generators is important both for basic understanding of the initial sequence of cortical activation and to identify landmarks for eloquent areas to spare in resection planning of cortex in epilepsy surgery. We now set out to localize the N20 and P22 using subdural grid recording with special emphasis on the question of the origin of P22: Brodmann area 4 versus area 1. Electroencephalographic dipole source analysis of the N20 and P22 responses obtained from subdural grids over the primary somatosensory cortex after median nerve stimulation was performed in four patients undergoing epilepsy surgery. Based on anatomical landmarks, equivalent current dipoles of N20 and P22 were localized posterior to (n = 2) or on the central sulcus (n = 2). In three patients, the P22 dipole was located posterior to the N20 dipole, whereas in one patient, the P22 dipole was located on the same coordinate in anterior-posterior direction. On average, P22 sources were found to be 6.6 mm posterior [and 1 mm more superficial] compared with the N20 sources. These data strongly suggest a postcentral origin of the P22 SEP component in Brodmann area 1 and render a major precentral contribution to the earliest stages of processing from the primary motor cortex less likely.

A robust current pattern for the detection of intraventricular hemorrhage in neonates using electrical impedance tomography

We compared two 16-electrode electrical impedance tomography (EIT) current patterns on their ability to reconstruct and quantify small amounts of bleeding inside a neonatal human head using both simulated and phantom data. The current patterns used were an adjacent injection RING pattern (with electrodes located equidistantly on the equator of a sphere) and an EEG current pattern based on the 10-20 EEG electrode layout. Structures mimicking electrically important structures in the infant skull were included in a spherical numerical forward model and their effects on reconstructions were determined. The EEG pattern was found to be a better topology to localize and quantify anomalies within lateral ventricular regions. The RING electrode pattern could not reconstruct anomaly location well, as it could not distinguish different axial positions. The quantification accuracy of the RING pattern was as good as the EEG pattern in noise-free environments. However, the EEG pattern showed better quantification ability than the RING pattern when noise was added. The performance of the EEG pattern improved further with respect to the RING pattern when a fontanel was included in forward models. Significantly better resolution and contrast of reconstructed anomalies was achieved when generated from a model containing such an opening and 50 dB added noise. The EEG method was further applied to reconstruct data from a realistic neonatal head model. Overall, acceptable reconstructions and quantification results were obtained using this model and the homogeneous spherical forward model.

Realistic and spherical head modeling for EEG forward problem solution: a comparative cortex-based analysis

The accuracy of forward models for electroencephalography (EEG) partly depends on head tissues geometry and strongly affects the reliability of the source reconstruction process, but it is not yet clear which brain regions are more sensitive to the choice of different model geometry. In this paper we compare different spherical and realistic head modeling techniques in estimating EEG forward solutions from current dipole sources distributed on a standard cortical space reconstructed from Montreal Neurological Institute (MNI) MRI data. Computer simulations are presented for three different four-shell head models, two with realistic geometry, either surface-based (BEM) or volume-based (FDM), and the corresponding sensor-fitted spherical-shaped model. Point Spread Function (PSF) and Lead Field (LF) cross-correlation analyses were performed for 26 symmetric dipole sources to quantitatively assess models' accuracy in EEG source reconstruction. Realistic geometry turns out to be a relevant factor of improvement, particularly important when considering sources placed in the temporal or in the occipital cortex.

Advances in spike localization with EEG dipole modeling

EEG interpretation by visual inspection of waveforms, using the assumption that activity at a given electrode is a representation of only the activity of the cortex immediately beneath it, has been the traditional form of EEG analysis since its inception. The relatively recent advent of digital EEG has allowed more advanced analysis of EEG data and has shown that the simple visual inspection described above is a simplistic form of analysis. This is especially true when one is attempting to localize an epileptogenic focus using EEG spikes or seizure onset data. Spatiotemporal analysis of scalp voltage fields has allowed for improved localization of likely cerebral origins of such waveforms. Equivalent dipole source modeling is one such technique and, although not perfect, provides improved characterization of spike and seizure sources as compared to previous methods when properly interpreted. The use of other modern techniques, such as 3D MRI reconstructions and realistic head models, can further improve accuracy of dipole localization and allow for the synthesis of EEG and imaging data, which may be invaluable, especially in cases of pre-surgical epilepsy evaluation.

Mismatch negativity with visual-only and audiovisual speech

The functional organization of cortical speech processing is thought to be hierarchical, increasing in complexity and proceeding from primary sensory areas centrifugally. The current study used the mismatch negativity (MMN) obtained with electrophysiology (EEG) to investigate the early latency period of visual speech processing under both visual-only (VO) and audiovisual (AV) conditions. Current density reconstruction (CDR) methods were used to model the cortical MMN generator locations. MMNs were obtained with VO and AV speech stimuli at early latencies (approximately 82-87 ms peak in time waveforms relative to the acoustic onset) and in regions of the right lateral temporal and parietal cortices. Latencies were consistent with bottom-up processing of the visible stimuli. We suggest that a visual pathway extracts phonetic cues from visible speech, and that previously reported effects of AV speech in classical early auditory areas, given later reported latencies, could be attributable to modulatory feedback from visual phonetic processing.

Development of a realistic head model for EEG event-detection and source localization in newborn infants

In this work we present an integrated method for electroencephalography (EEG) source localization in newborn infants, based on a realistic head model. To build a realistic head model we propose an interactive hybrid segmentation method for T1 magnetic resonance images (MRI), consisting of active contours, fuzzy c-means (FCM) clustering and mathematical morphology. Subsequently, we solve the localization problem using a spike train detection algorithm and an algorithm that deals with the forward and inverse problem. The performance of this fused method indicates that our realistic head model is suitable for the accurate localization of the EEG activity. We will present both initial qualitative and quantitative results.

Rapid interactions between the ventral visual stream and emotion-related structures rely on a two-pathway architecture

Visual attention can be driven by the affective significance of visual stimuli before full-fledged processing of the stimuli. Two kinds of models have been proposed to explain this phenomenon: models involving sequential processing along the ventral visual stream, with secondary feedback from emotion-related structures ("two-stage models"); and models including additional short-cut pathways directly reaching the emotion-related structures ("two-pathway models"). We tested which type of model would best predict real magnetoencephalographic responses in subjects presented with arousing visual stimuli, using realistic models of large-scale cerebral architecture and neural biophysics. The results strongly support a "two-pathway" hypothesis. Both standard models including the retinotectal pathway and nonstandard models including cortical-cortical long-range fasciculi appear plausible.

Spatial relationship of source localizations in patients with focal epilepsy: Comparison of MEG and EEG with a three spherical shells and a boundary element volume conductor model

Epilepsy surgery is an option for patients with pharmacoresistant focal epilepsies, but it requires a precise focus localization procedure. Magnetoencephalography (MEG) and electroencephalography (EEG) can be used for analysis of interictal activity. The aim of this prospective study was to compare clusters of source localization results with MEG and EEG using a three spherical shells (3SS) and a boundary element method (BEM) volume conductor model. The study was closed when 100 patients met the inclusion criteria. Simultaneous MEG and EEG were recorded during presurgical evaluation. Epileptiform signals were analyzed using an equivalent current dipole model. Centroids of source localizations from MEG, EEG, 3SS, and BEM in their respective combinations were compared. In a 3SS model, MEG source localizations were 5.6 mm inferior to those obtained by EEG, while in a BEM model MEG source localizations were 6.3 mm anterior and 4.8 mm superior. The mean scattering of source localizations between both volume conductor models was 19.5 mm for EEG and 9.6 mm for MEG. For MEG no systematic difference between BEM and 3SS source localizations was found. For EEG, source localizations with BEM were 5.9 mm posterior and 11.7 mm inferior to those determined using 3SS. No differences were found between the 46 temporal and the 54 extratemporal lobe epilepsy patients. The observed systematic differences of source localizations of epileptic spikes due to the applied source signal modality and volume conductor model should be considered in presurgical evaluation when only one source signal and volume conductor model is available.

Clinical Application of Dipole Models in the Localization of Epileptiform Activity

SUMMARY

Routine clinical interpretation of EEG using visual inspection of traces is a time-honored, but simplistic, form of analysis. This is particularly true in attempts to localize an epileptogenic focus by means of EEG spike or seizure waveforms. Improved understanding of the cortical substrates of these potentials has allowed us to identify their likely cerebral origins through spatio-temporal analysis of scalp voltage fields. Equivalent dipole modeling is one such technique. Although an imperfect representation of spike or seizure sources, proper interpretation of dipole models can lead to a far better characterization of their localization and propagation. Modern techniques of 3-D MRI reconstruction and realistic head models have both improved localization accuracy and provided a means of displaying results in an image of the individual's brain.

Computational aspects of the EEG forward problem solution for real head model using finite element method

The real head model has been used for the accurate scalp potential modeling. The realistic shapes of head tissues were derived from a set of 2-D magnetic resonance images (MRI) by extracting surface boundaries for the major tissues such as the scalp, the skull, the cerebrospinal fluid (CSF), the white matter, and the gray matter. From boundary data a 3-D volume generic head model has been constructed and a mesh for an arbitrary complexity head shape has been generated for finite-element method (FEM) modeling. This paper first addresses the use of this realistic finite elements head model to solve the EEG forward problem. The accuracy and computational time of the potential modeling are then examined.

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.

Auditory change detection in schizophrenia: sources of activity, related neuropsychological function and symptoms in patients with a first episode in adolescence, and patients 14 years after an adolescent illness-onset

BACKGROUND

The event-related brain response mismatch negativity (MMN) registers changes in auditory stimulation with temporal lobe sources reflecting short-term echoic memory and frontal sources a deviance-induced switch in processing. Impairment, controversially present at the onset of schizophrenia, develops rapidly and can remain independent of clinical improvement. We examined the characteristics of the scalp-recorded MMN and related these to tests of short-term memory and set-shifting. We assessed whether the equivalent dipole sources are affected already at illness-onset in adolescence and how these features differ after a 14-year course following an adolescent onset. The strength, latency, orientation and location of frontal and temporal lobe sources of MMN activity early and late in the course of adolescent-onset schizophrenia are analysed and illustrated.

METHODS

MMN, a measure of auditory change-detection, was elicited by short deviant tones in a 3-tone oddball-presentation and recorded from 32 scalp electrodes. Four dipole sources were placed following hypothesis-led calculations using brain electrical source analysis on brain atlas and MR-images. A short neuropsychological test battery was administered. We compared 28 adolescent patients with a first episode of schizophrenia and 18 patients 14 years after diagnosis in adolescence with two age-matched control groups from the community (n = 22 and 18, respectively).

RESULTS

MMN peaked earlier in the younger than the older subjects. The amplitude was reduced in patients, especially the younger group, and was here associated with negative symptoms and slow set-shifting. In first-episode patients the temporal lobe sources were more ventral than in controls, while the left cingular and right inferior-mid frontal sources were more caudal. In the older patients the left temporal locus remained ventral (developmental stasis), the right temporal locus extended more antero-laterally (illness progression), and the right frontal source moved antero-laterally (normalised).

CONCLUSION

From the start of the illness there were differences in the dipole-model between healthy and patient groups. Separate characteristics of the sources of the activity differences showed an improvement, stasis or deterioration with illness-duration. The precise nature of the changes in the sources of MMN activity and their relationship to selective information processing and storage depend on the specific psychopathology and heterogeneous course of the illness.

Maturation processes in automatic change detection as revealed by event-related brain potentials and dipole source localization: significance for adult AD/HD

Mismatch negativity (MMN) is an event-related potential reflecting automatic attention-related information processing marking the detection of auditory change. The bilateral scalp distribution develops by 14 years of age, and is elicited with adult latencies by 17 years. But consistent with reports of continued brain maturation after adolescence, we show here that features of the temporal and frontal lobe dipole sources also continue to develop in the third decade of life. This has consequences for studies of the developmental course of MMN anomalies, from childhood into adulthood, in attention-deficit/hyperactivity disorder. Two groups of healthy subjects with mean ages of 17 and 30 years were presented with a 3-tone auditory oddball. The duration-deviant MMN was recorded during attention to a visual discrimination (auditory-passive condition) and an active auditory discrimination. MMN amplitudes were smaller in the older subjects and the MMN lasted longer over the right hemisphere. Latencies and moments of the four dipoles in the temporal and frontal lobes did not distinguish the two subject-groups. But both temporal lobe sources were located significantly more ventrally and further left in the young adult than in the adolescent subjects. The left cingular source moved posteriorly and the right inferior frontal source moved antero-medially in the older subjects. Brain development in the third decade may cause the two frontal sources to move apart on the rostro-caudal axis but the temporal lobe sources to move left on the lateral and down on the dorsoventral axes. Thus special care is necessary in interpreting putative dysfunctional neurobiological changes in developmental attention-deficit disorders where as-yet-unspecified sub-groups may show a late developmental lag, partial lag, or no lag at all, associated with other impairments.

A finite difference method with reciprocity used to incorporate anisotropy in electroencephalogram dipole source localization

Many implementations of electroencephalogram (EEG) dipole source localization neglect the anisotropical conductivities inherent to brain tissues, such as the skull and white matter anisotropy. An examination of dipole localization errors is made in EEG source analysis, due to not incorporating the anisotropic properties of the conductivity of the skull and white matter. First, simulations were performed in a 5 shell spherical head model using the analytical formula. Test dipoles were placed in three orthogonal planes in the spherical head model. Neglecting the skull anisotropy results in a dipole localization error of, on average, 13.73 mm with a maximum of 24.51 mm. For white matter anisotropy these values are 11.21 mm and 26.3 mm, respectively. Next, a finite difference method (FDM), presented by Saleheen and Kwong (1997 IEEE Trans. Biomed. Eng. 44 800-9), is used to incorporate the anisotropy of the skull and white matter. The FDM method has been validated for EEG dipole source localization in head models with all compartments isotropic as well as in a head model with white matter anisotropy. In a head model with skull anisotropy the numerical method could only be validated if the 3D lattice was chosen very fine (grid size < or = 2 mm).

Frontal and temporal sources of mismatch negativity in healthy controls, patients at onset of schizophrenia in adolescence and others at 15 years after onset

Mismatch negativity (MMN) is an event-related potential measure of auditory change detection. It is widely reported to be smaller in patients with schizophrenia and may not improve along with otherwise successful clinical treatment. The main aim of this report is to explore ways of measuring and presenting four features of frequency-deviant MMN dipole sources (dipole moment, peak latency, brain location and orientation) and to relate these to the processes of psychopathology and illness progression. Data from early onset patients (EOS) at the start of the illness in adolescence, and others who had their first break in adolescence 15 years ago (S-15Y) were compared with two groups of age-matched healthy controls (C-EOS, C-15Y). A four-source model fitted the MMN waveform recorded from all four groups, whether MMN amplitude was more (EOS) or less (S-15Y) reduced. The locations were in the left superior temporal and anterior cingulate gyri, right superior temporal and inferior/mid frontal cortices. Dipole latencies confirmed a bottom-up sequence of processing and dipole moments were larger in the temporal lobes and on the left. Patients showed small dipole location changes that were more marked in the S-15Y than the EOS group (more rostral for the left anterior cingulate, more caudal for the right mid-frontal dipole) consistent with illness progression. The modelling of MMN dipole sources on brain atlas and anatomical images suggests that there is a degree of dissociation during illness between small progressive anatomical changes and some functional recovery indexed by scalp recordings from patients with an onset in adolescence 15 years before compared to adolescents in their first episode.

Tridimensional animated brain mapping from conventional paper-ink EEG recordings

One of the most powerful functional evaluations of the electrical activity of the brain is the EEG imaging, but wide clinical use is limited by its costs. It is also of clinical, academic and scientific interest to obtain brain electrical maps from old paper/ink, patient recordings. The aim of the present study was the development of a computer system designed to obtain bidimensional and tridimensional maps, continuous movie map display and mosaic presentation from conventional paper/ink EEG recordings. The wave amplitude was manually measured with a translucent template from conventional 8- to 16-channel EEG paper recordings using 10-20 monopolar montage for one or both hemispheres. The computer system allows the selection of the number and location of electrodes, input of amplitude values, and the map display mode. The electrical brain field was generated from amplitude measurements by a spherical splines interpolation algorithm on a conventional Pentium-based computer. The interpolated surface was represented on a semi-sphere modeled skull. The EEG maps displayed with pseudo-color or gray scales can be rotated, zoomed in or zoomed out and/or printed for clinical reports. Movie animation or mosaic display of space-temporal EEG voltage changes were generated by processing sequential amplitude measurements. This system represents a cost-effective method for EEG mapping from conventional paper/Ink EEG equipments.

Reconstruction of 3-D Head Geometry From Digitized Point Sets: An Evaluation Study

In this paper, we evaluate different methods to estimate patient-specific scalp, skull, and brain surfaces from a set of digitized points from the target's scalp surface. The reconstruction problem is treated as a registration problem: An a priori surface model, consisting of the scalp, skull, and brain surfaces, is registered to the digitized surface points. The surface model is generated from segmented magnetic resonance (MR) volume images. We study both affine and free-form deformation (FFD) registration, the use of average models, the averaging of individual registration results, a model selection procedure, and statistical deformation models. The registration algorithms are mainly previously published, and the objective of this paper is to evaluate these methods in this particular application with sparse data. The main interest of this paper is to generate geometric head models for biomedical applications, such as electroencephalography and magnetoencephalographic. However, the methods can also be applied to other anatomical regions and to other application areas. The methods were validated using 15 MR volume images, from which the scalp, skull, and brain were manually segmented. The best results were achieved by averaging the results of the FFD registrations of the database: the mean distance from the manually segmented target surface to a deformed a priori model surface for the studied anatomical objects was 1.68-2.08 mm, depending on the point set used. The results support the use of the evaluated methods for the reconstruction of geometric models in applications with sparse data.

Electric Source Imaging in Temporal Lobe Epilepsy

The objective of this study was to determine the validity of interictal spike (IIS) source localization in temporal lobe epilepsies (TLE) using stereoelectroencephalography as a validating method. Twenty patients with drug-resistant TLE were studied with high-resolution EEG and stereoelectroencephalography. Sixty-four scalp channels, a realistic head model, and different algorithms were used. For each patient, the intracerebral interictal distribution was studied and classified into one of three groups: L (mainly lateral), ML (mediolateral), and M (medial). In group L (three patients), surface IIS were recorded with a high signal-to-noise ratio. Source localizations designated all or part of the intracerebral interictal distribution. In group ML (11 patients), 8 patients had surface IIS, only 5 of which were localizable. High-resolution EEG permitted localization of the more lateral portion and definition of its rostrocaudal extension. A common pattern was identified in three patients with a predominant role of the temporal pole. In group M (six patients), four patients had rare surface IIS, none of which were localizable. Surface EEG does not record IIS limited to medial temporal lobe structures. In TLE with a mediolateral or a lateral interictal distribution, only the lateral component is detectable on surface EEG and accurately localizable by source localization tools.

Controversies in neurophysiology. MEG is superior to EEG in localization of interictal epileptiform activity: Pro

Both EEG and magnetoencephalography (MEG), with a time resolution of 1 ms or less, provide unique neurophysiologic data not obtainable by other neuroimaging techniques. MEG and EEG have often been compared to each other now although the two are complementary. Now that MEG has emerged as a mature clinical technology, it is worthwhile to compare the relative strengths of each for the localization of interictal epileptiform activity and to describe the strengths of MEG relative to EEG in the localization of interictal epileptiform activity. The sources of MEG and EEG signals will first be reviewed. Issues relevant to solving the forward problem and the inverse problem in MEG and EEG will be addressed followed by a comparison of research concerning the detection and localization of interictal epileptiform activity by MEG and EEG. The emphasis will be upon techniques and software routinely used in clinical applications but some emerging areas of MEG research which are entering clinical practice will also be reviewed.

SIGNIFICANCE

MEG is a new noninvasive neurophysiologic technique which provides unique information for the clinical evaluation of patients with epilepsy, revealing aspects of neuronal function that previously could only be obtained by invasive EEG monitoring, and giving a new window for research of neuronal activity.

Effects of skull thickness, anisotropy, and inhomogeneity on forward EEG/ERP computations using a spherical three-dimensional resistor mesh model

Bone thickness, anisotropy, and inhomogeneity have been reported to induce important variations in electroencephalogram (EEG) scalp potentials. To study this effect, we used an original three-dimensional (3-D) resistor mesh model described in spherical coordinates, consisting of 67,464 elements and 22,105 nodes arranged in 36 different concentric layers. After validation of the model by comparison with the analytic solution, potential variations induced by geometric and electrical skull modifications were investigated at the surface in the dipole plane and along the dipole axis, for several eccentricities and bone thicknesses. The resistor mesh permits one to obtain various configurations, as local modifications are introduced very easily. This has allowed several head models to be designed to study the effects of skull properties (thickness, anisotropy, and heterogeneity) on scalp surface potentials. Results show a decrease of potentials in bone, depending on bone thickness, and a very small decrease through the scalp layer. Nevertheless, similar scalp potentials can be obtained using either a thick scalp layer and a thin skull layer, and vice versa. It is thus important to take into account skull and scalp thicknesses, because the drop of potential in bone depends on both. The use of three different layers for skull instead of one leads to small differences in potential values and patterns. In contrast, the introduction of a hole in the skull highly increases the maximum potential value (by a factor of 11.5 in our case), because of the absence of potential drop in the corresponding volume. The inverse solution without any a priori knowledge indicates that the model with the hole gives the largest errors in both position and dipolar moment. Our results indicate that the resistor mesh model can be used as a robust and user-friendly simulation tool in EEG or event-related potentials. It makes it possible to build up real head models directly from anatomic magnetic resonance imaging without tessellation, and is able to take into account head heterogeneities very simply by changing volume elements conductivity. Hum. Brain Mapping 21:84-95, 2004.

Temporal and spatial congruence of components of motion-onset evoked responses investigated by whole-head magneto-electroencephalography

Motion-onset related components in averaged whole head co-recorded MEG and EEG responses of 24 adults to a low-contrast checkerboard pattern were studied. The aims were to identify these components, to characterize quantitatively their maps and to localize the underlying sources by equivalent-current-dipole (ECD) analyses with a spherical head model.After a weak P1, a large start-elicited negativity arises, comprising the novel N2a (occipital positive and parieto-central negative, peak-latency 141 ms) and the N2 like N2b (bilateral parieto-temporal, 175 ms) component. It is followed by a large positive stop-related component, P2 (156 ms after motion-offset). The corresponding MEG components N2am and N2bm showed bilateral dipole fields with considerable overlap. P1m has a single dipole field around the midline. N2a(m) and N2b(m) can be modelled with two bilateral ECDs with significant different locations. The study shows that accurate mapping and ECD analyses can distinguish two neighbouring areas of the visual cortex, 21+/-4 (SE) mm separated, which activities are reflected in both spatio-temporally closely related N2(m) components. N2a(m) and N2b(m) originate in the extrastriate cortex, possibly close to or in V3/V3A and MT/V5 respectively. Motion-evoked activity in (near) V3/V3A is novel on the basis of EEG data.

The 3D topography of MEG source localization accuracy: effects of conductor model and noise

OBJECTIVE

To evaluate the effect that different head conductor models have on the source estimation accuracy of magnetoencephalography (MEG) under realistic conditions.

METHODS

Magnetic fields evoked by current dipoles were simulated using a highly refined 3-layer realistically shaped conductor model. Noise from a real MEG measurement was added to the simulated fields. Source parameters (location, strength, orientation) were estimated from the noisy signals using 3 spherically symmetric models and several one- and 3-layer realistically shaped boundary-element models. The effect of different measurement sensors (gradiometers, magnetometers) was also tested.

RESULTS

The noise typically present in brain signals masked the errors due to the different conductor models so that in most situations the models gave comparable results. Active cortical areas around the vertex and in the temporal, frontoparietal, and occipital regions were typically found with 2-4 mm accuracy, whereas source localization in several anterior frontal lobe and deep brain structures yielded errors exceeding 2 cm. Localization in anterior frontal regions may benefit most from the use of realistically shaped models.

CONCLUSIONS

The traditionally used sphere model is an adequate model for most research purposes. Any means that increase the signal-to-noise ratio are of highest importance in attempting to improve the source estimation accuracy.

Electrical impedance tomography of human brain function using reconstruction algorithms based on the finite element method

Electrical impedance tomography (EIT) is a recently developed technique which enables the internal conductivity of an object to be imaged using rings of external electrodes. In a recent study, EIT during cortical evoked responses showed encouraging changes in the raw impedance measurements, but reconstructed images were noisy. A simplified reconstruction algorithm was used which modelled the head as a homogeneous sphere. In the current study, the development and validation of an improved reconstruction algorithm are described in which realistic geometry and conductivity distributions have been incorporated using the finite element method. Data from computer simulations and spherical or head-shaped saline-filled tank phantoms, in which the skull was represented by a concentric shell of plaster of Paris or a real human skull, have been reconstructed into images. There were significant improvements in image quality as a result of the incorporation of accurate geometry and extracerebral layers in the reconstruction algorithm. Image quality, assessed by blinded subjective expert observers, also improved significantly when data from the previous evoked response study were reanalysed with the new algorithm. In preliminary images collected during epileptic seizures, the new algorithm generated EIT conductivity changes which were consistent with the electrographic ictal activity. Incorporation of realistic geometry and conductivity into the reconstruction algorithm significantly improves the quality of EIT images and lends encouragement to the belief that EIT may provide a low-cost, portable functional neuroimaging system in the foreseeable future.

Comparison of performance of spherical and realistic head models in dipole localization from noisy EEG

The performance of a three-shell spherical head model versus the performance of a realistic head model is investigated when solving the inverse problem with a single dipole, in the presence of noise. This is evaluated by calculating the average dipole location error for 1000 noisy scalp potential sets, originating from the same test dipole and having the same noise level. The average location errors are obtained utilizing a local linearization, which is validated with a Monte-Carlo simulation. When the difference between the average location error utilizing a spherical and a realistic head model, represented by deltaR, is large for a large number of test dipoles, then it is worth using the more computationally demanding realistic head model. However, if deltaR is small for a large number of test dipoles, then it does not matter which model is used. For 27 electrodes, an electroencephalogram (EEG) epoch of one time sample and spatially white Gaussian noise, we found that the importance of the realistic head model over the spherical head model reduces by increasing the noise level. We further found that increasing the number of scalp electrodes from 27 to 44 has limited impact on the importance of the realistic head model over the spherical head model in EEG dipole source analysis. By increasing the number of time samples to six, the performance of the realistic head model in the inverse calculation gains importance compared with the three-shell spherical head model. Finally, we used spatially and temporally correlated background EEG instead of Gaussian noise. The advantage of the realistic head model over the spherical head model is reduced when applying correlated noise compared to Gaussian noise.