A cortical potential imaging study from simultaneous extra- and intracranial electrical recordings by means of the finite element method

In the present study, we have validated the cortical potential imaging (CPI) technique for estimating cortical potentials from scalp EEG using simultaneously recorded electrocorticogram (ECoG) in the presence of strong local inhomogeneity, i.e., Silastic ECoG grid(s). The finite element method (FEM) was used to model the realistic postoperative head volume conductor, which includes the scalp, skull, cerebrospinal fluid (CSF) and brain, as well as the Silastic ECoG grid(s) implanted during the surgical evaluation in epilepsy patients, from the co-registered magnetic resonance (MR) and computer tomography (CT) images. A series of computer simulations were conducted to evaluate the present FEM-based CPI technique and to assess the effect of the Silastic ECoG grid on the scalp EEG forward solutions. The present simulation results show that the Silastic ECoG grid has substantial influence on the scalp potential forward solution due to the distortion of current pathways in the presence of the extremely low conductive materials. On the other hand, its influence on the estimated cortical potential distribution is much less than that on the scalp potential distribution. With appropriate numerical modeling and inverse estimation techniques, we have demonstrated the feasibility of estimating the cortical potentials from the scalp EEG with the implanted Silastic ECoG gird(s), in both computer simulations and in human experimentation. In an epilepsy patient undergoing surgical evaluation, the cortical potentials were reconstructed from the simultaneously recorded scalp EEG, in which main features of spatial patterns during interictal spike were preserved and over 0.75 correlation coefficient value was obtained between the recorded and estimated cortical potentials. The FEM-based CPI technique provides a means of connecting the simultaneous recorded ECoG and the scalp EEG and promises to become an effective tool to evaluate and validate CPI techniques using clinic data.     

Movement related activity in the high gamma range of the human EEG

Electrocorticographic (ECoG) recordings obtained using intracranially implanted electrodes in epilepsy patients indicate that high gamma band (HGB) activity of sensorimotor cortex is focally increased during voluntary movement. These movement related HGB modulations may play an important role in sensorimotor cortex function. It is however currently not clear to what extent this type of neural activity can be detected using non-invasive electroencephalography (EEG) and how similar HGB responses in healthy human subjects are to those observed in epilepsy patients. Using the same arm reaching task, we have investigated spectral power changes both in intracranial ECoG recordings in epilepsy patients and in non-invasive EEG recordings optimized for detecting HGB activity in healthy subjects. Our results show a common HGB response pattern both in ECoG and EEG recorded above the sensorimotor cortex contralateral to the side of arm movement. In both cases, HGB activity increased around movement onset in the 60-90 Hz range and became most pronounced at reaching movement end. Additionally, we found EEG HGB activity above the frontal midline possibly generated by the anterior supplementary motor area (SMA), a region that was however not covered by the intracranial electrodes used in the present study. In summary, our findings show that HGB activity from human sensorimotor cortex can be non-invasively detected in healthy subjects using EEG, opening a new perspective for investigating the role of high gamma range neuronal activity both in function and dysfunction of the human cortical sensorimotor network.     

颅内电极监测癫痫发作初始期脑电定位意义

目的探讨颅内电极监测癫痫发作初始期脑电对致痫灶的定位价值。方法14例难治性癫痫患者,其临床、影像学检查及头皮脑电记录等不能精确定位致痫灶,经颅骨钻孔、埋置颅内深部和/或皮层电极,长程脑电监测并分析发作初始期脑电的异常放电节律和范围,对比术后病理学检查和随访结果,分析发作初期脑电定位致痫灶的准确性。结果记录到了47次临床发作。发作初期异常放电可分为4种形式:低幅快节律、尖波节律、高幅棘波节律和棘慢波节律。按异常放电范围可分为:局限性放电、区域性放电和广泛性放电。11例(78.5%)患者准确定位了致痫灶,术后癫痫发作消失(64.3%)或极少发作(14.2%),2例(14.2%)发作减少90%以上,1例(7.1%)术后无改善。结论颅内埋置电极脑电监测癫痫发作初始期异常放电形式、部位和异常放电范围是定位致痫灶可靠的方法。     

Human intracranially-recorded cortical responses evoked by painful electrical stimulation of the sural nerve

Intracranial recordings were obtained from 5 epilepsy patients to help identify the generators of the scalp somatosensory evoked potential (SEP) components that appear to be involved in orienting attention towards a potentially threatening, painful sural nerve electrical stimulus. The intracranial recording data support, for the most part, the generators suggested by our scalp SEP studies. The generators of the central negativity at 70-110 ms post-stimulus and the contralateral temporal negativity at 100-180 ms are located in the somatosensory association areas in the medial wall of the parietal cortex and in the parietal operculum and insula, respectively. The negative potential at 130-200 ms recorded from over the fronto-central scalp appears to be generated in the medial prefrontal cortex and primary somatosensory cortex foot area. The intracranial recording data suggest that the positive scalp potential at 280-320 ms, which corresponds to the pain-related P2, has multiple generators, including the anterior cingulate cortex, inferior parietal cortex, and possibly the somatosensory association areas in the medial wall of the parietal cortex. Finally, the positive scalp potential at 320-400 ms has generators in brain areas that others have shown to generate the P3a, including the dorsolateral and medial prefrontal cortices, temporal parietal junction, and the posterior hippocampus, which supports our hypothesis that this potential is a pain-evoked P3a. The putative functional roles of the brain areas generating these components and the response properties of the intracranial peaks recorded from these brain areas are in most cases consistent with the attention- and pain-related properties of their corresponding scalp SEP components. The intracranial recordings also demonstrate that the source configuration underlying the SEP components are often more complex than was suggested from the scalp studies. This complexity implies that the dipole source localization analysis of these components will at best provide only a very crude estimate of source location and activity, and that caution must be used when interpreting a change in the scalp component amplitude.     

Magnetoencephalographic yield of interictal spikes in temporal lobe epilepsy. Comparison with scalp EEG recordings

To compare magnetoencephalography (MEG) with scalp electroencephalography (EEG) in the detection of interictal spikes in temporal lobe epilepsy (TLE), we simultaneously recorded MEG and scalp EEG with a whole-scalp neuromagnetometer in 46 TLE patients. We visually searched interictal spikes on MEG and EEG channels and classified them into three types according to their presentation on MEG alone (M-spikes), EEG alone (E-spikes), or concomitantly on both modalities (M/E-spikes). The M-spikes and M/E-spikes were localized with MEG equivalent current dipole modeling. We analyzed the relative contribution of MEG and EEG in the overall yield of spike detection and also compared M-spikes with M/E-spikes in terms of dipole locations and strengths. During the 30- to 40-min MEG recordings, interictal spikes were obtained in 36 (78.3%) of the 46 patients. Among the 36 patients, most spikes were M/E-spikes (68.3%), some were M-spikes (22.1%), and some were E-spikes (9.7%). In comparison with EEG, MEG gave better spike yield in patients with lateral TLE. Sources of M/E- and M-spikes were situated in the same anatomical regions, whereas the average dipole strength was larger for M/E- than M-spikes. In conclusion, some interictal spikes appeared selectively on either MEG or EEG channels in TLE patients although more spikes were simultaneously identified on both modalities. Thus, simultaneous MEG and EEG recordings help to enhance spike detection. Identification of M-spikes would offer important localization of irritative foci, especially in patients with lateral TLE.     

Localization of abnormal discharges causing insular epilepsy by magnetoencephalography

The insula, one of the five cerebral lobes of the brain, is located deep within the brain and lies mainly beneath the temporal lobe. Insular epilepsy can be easily confused and misdiagnosed as temporal lobe epilepsy (TLE) because of the similar clinical symptoms and scalp electroencephalography (EEG) findings due to the insula location and neuronal connections with the temporal lobe. Magnetoencephalography (MEG) has higher sensitivity and spatial resolution than scalp EEG, and thus can often identify epileptic discharges not revealed by scalp EEG. Simultaneous scalp EEG and MEG were performed to detect and localize epileptic discharges in two patients known to have insular epilepsy associated with cavernous angioma in the insula. Epileptic discharges were detected as abnormal spikes in the EEG and MEG findings. In Patient 1, the sources of all MEG spikes detected simultaneously by EEG and MEG (E/M-spikes) were localized in the anterior temporal lobe, similar to TLE. In contrast, the sources of all MEG spikes detected only by MEG (M-spikes) were adjacent to the insular lesion. In Patient 2, the sources of all MEG spikes detected simultaneously by EEG and MEG (E/M-spikes) were localized in the anterior temporal lobe. These findings indicate that MEG allows us to detect insular activity that is undetectable by scalp EEG. In conclusion, simultaneous EEG and MEG are helpful for detecting spikes and obtaining additional information about the epileptic origin and propagation in patients with insular epilepsy.     

Towards source volume estimation of interictal spikes in focal epilepsy using magnetoencephalography

Interictal spikes are a hallmark of cortical epileptogenicity; their spatial distribution in the cortex defines the so-called ‘irritative’ zone or spiking volume (SV). Delineating the SV precisely is a challenge during the presurgical evaluation of patients with epilepsy. Magnetoencephalography (MEG) recordings enable determination of the brain sources of epileptic spikes using source localization procedures. Most previous clinical MEG studies have relied on dipole modeling of epileptic spikes, which does not permit a volumetric estimation of the spiking cortex.In the present study, we propose a new source modeling procedure, Volumetric Imaging of Epileptic Spikes (VIES). In VIES, the SV is identified as the 3D region where sources of the high frequency activities (> 20 Hz) associated with epileptic spikes are distributed. We localized these sources using a beamforming approach (DICS, Dynamic Imaging of Coherent Neural Sources). To determine the optimal parameters and accuracy of the method, we compared the SV obtained by VIES with the SV defined by the invasive gold standard, intracranial stereotactic EEG recordings (SEEG), in 21 patients with focal epilepsy. Using rigorous validation criteria based on the exact anatomical location of SEEG contacts, we found that the overall sensitivity of VIES for detecting spiking SEEG contacts was 76% and its specificity for correctly identifying non-spiking SEEG contacts was 67%, indicating a good agreement between VIES and SEEG. Moreover, we found that classical dipole clustering was not informative in 9/21 patients, while VIES enable to delineate the SV in all patients. For the 12 patients having a SV delineated both with VIES and dipole clustering, VIES method had higher sensitivity and lower specificity. This proof-of-concept study shows that VIES is a promising approach to non-invasive estimation of the SV in focal epilepsy.     

Dynamic imaging of ictal oscillations using non-invasive high-resolution EEG

Scalp electroencephalography (EEG) has been established as a major component of the pre-surgical evaluation for epilepsy surgery. However, its ability to localize seizure onset zones (SOZ) has been significantly restricted by its low spatial resolution and indirect correlation with underlying brain activities. Here we report a novel non-invasive dynamic seizure imaging (DSI) approach based upon high-density EEG recordings. This novel approach was particularly designed to image the dynamic changes of ictal rhythmic discharges that evolve through time, space and frequency. This method was evaluated in a group of 8 epilepsy patients and results were rigorously validated using intracranial EEG (iEEG) (n=3) and surgical outcome (n=7). The DSI localized the ictal activity in concordance with surgically resected zones and ictal iEEG recordings in the cohort of patients. The present promising results support the ability to precisely and accurately image dynamic seizure activity from non-invasive measurements. The successful establishment of such a non-invasive seizure imaging modality for surgical evaluation will have a significant impact in the management of medically intractable epilepsy.     

Epilepsy surgery

If the first two or three antiepileptic drugs used do not control epilepsy, there is little chance that subsequent medications will be effective. In individuals with refractory focal epilepsy, neurosurgery can have a 60-70% chance of bringing long-term remission and these cases should be referred to a specialised centre for evaluation. The standard evaluation includes clinical review, brain imaging with magnetic resonance imaging, recording of seizures with prolonged scalp electroencephalography (EEG) and video, neuropsychological and psychiatric assessments. The aim is to establish converging evidence that there is a single epileptic focus and that the rest of the brain is functioning normally. In some individuals further evaluation with functional imaging and intracranial EEG recordings may be necessary. The most commonly performed resective operation is an anterior temporal lobe resection to remove a sclerotic hippocampus, followed by lesionectomies and neocortical resections. Palliative manoeuvres, to reduce seizure frequency and severity include corpus callosotomy, subpial transection and vagal nerve stimulation.     

Brain electrical tomography (BET) analysis of induced gamma band responses during a simple object recognition task

The formation of cortical object representations requires the activation of cell assemblies, correlated by induced oscillatory bursts above 20 Hz (gamma band), which are characterized by trial-by-trial latency fluctuations around a mean of approximately 300 ms after stimulus onset. The present electroencephalogram (EEG) study was intended to uncover to the generators of induced gamma band responses (GBRs) and to analyze phase-synchronization between these sources. A standard object recognition task was used to elicit gamma activity. At the scalp surface (electrode space), we found an augmentation of induced GBRs after the presentation of meaningful (familiar) as opposed to meaningless (unfamiliar) stimuli, which was accompanied by a dense pattern of significant phase-locking values between distant recording sites. Subsequently, intracranial current density distributions compatible with the observed scalp voltage topographies were estimated by means of VARETA (Variable Resolution Electromagnetic Tomography). In source space brain electrical tomographies (BETs) revealed widespread generators of induced GBRs at temporal, parietal, posterior, and frontal areas. Phase-locking analysis was calculated between re-constructed electrode signals based on separate forward solutions of the observed generators, thereby eliminating the possibly confounding influence of activity from areas not under observation. The results support the view that induced GBRs signify synchronous neuronal activity in a broadly distributed network during object recognition. The localization of the generators of event-related potentials (ERPs), evoked gamma activity, and induced alpha activity revealed different sources as compared to the induced GBR and, thus, seem to mirror complementary functions during the present task as compared to induced high-frequency brain dynamics.     

Frequency-specific electrocorticographic correlates of working memory delay period fMRI activity

Electrocorticography (ECoG) and functional MRI (BOLD-fMRI) have been used previously to measure brain activity during working memory delay periods. These studies have separately reported oscillation changes in the theta (4-8 Hz) band and BOLD-fMRI increases during delay periods when information is maintained in memory. However, it is not known how intracranial cortical field potential (CFP) changes relate to BOLD-fMRI responses during delay periods. To answer this question, fMRI was obtained from six epilepsy patients during a visual working memory task. Then, following subdural macroelectrode implant, continuous ECoG was used to record CFPs during the same task. Time-frequency analyses showed delay period gamma band oscillation amplitude increases on electrodes located near fMRI activity, while in the theta band changes were higher for electrodes located away from fMRI activation. The amplitude of the ECoG gamma band response was significantly positively correlated with the fMRI response, while a negative correlation was found for the theta band. The findings are consistent with previous reports of local field potential (LFP) coupling in the gamma band with BOLD-fMRI responses during visual stimulation in monkeys, but are novel in that the relationship reported here persists after the disappearance of visual stimuli while information is being maintained in memory. We conclude that there is a relationship between BOLD-fMRI increases and human working memory delay period gamma oscillation increases and theta decreases. The spectral profile change provides a basis for comparison of working memory delay period BOLD-fMRI with field potential recordings in animals and other human intracranial EEG studies.     

脑磁图在癫(癎)诊断和治疗中的应用

目的应用脑磁图(MEG)技术对致灶进行定位,比较其与头皮脑电图(EEG)在神经影像学方法的应用价值,探讨MEG技术对癫灶定位的应用前景。方法对113例癫患者进行手术治疗,术前均通过临床症候学、头皮EEG、MRI、MEG检查,进行MEG与其他检查方法和临床症候学在定侧、定叶诊断的对比研究。手术在皮层EEG及脑深部EEG监测下进行,手术治疗结果以Engel疗效分级评价。所有手术标本常规行光镜检查。结果113例患者中MEG定位局限于单个叶的为91例,头皮EEG仅为30例。术前致灶定位依据多种检查结果和临床症候学综合定位。MEG与MRI、临床症候学在定侧诊断准确性方面比较差异有显著性。MEG与MRI在定叶诊断准确性方面的比较差异有显著性。36例术前头皮EEG表现为双侧或全导癫波的癫患者,其中有34例MEG表现为单侧癫波,具有定侧诊断的意义。结论头皮EEG、MRI、正电子发射计算机断层扫描(PET)、临床症候学均不足以做出独立精确的致灶定位诊断,综合比较MEG比上述方法定侧和定叶的准确性高。MEG空间分辨率、时间分辨率高,有助于区分致灶和镜灶。     

Five-dimensional neuroimaging: localization of the time-frequency dynamics of cortical activity

The spatiotemporal dynamics of cortical oscillations across human brain regions remain poorly understood because of a lack of adequately validated methods for reconstructing such activity from noninvasive electrophysiological data. In this paper, we present a novel adaptive spatial filtering algorithm optimized for robust source time-frequency reconstruction from magnetoencephalography (MEG) and electroencephalography (EEG) data. The efficacy of the method is demonstrated with simulated sources and is also applied to real MEG data from a self-paced finger movement task. The algorithm reliably reveals modulations both in the beta band (12-30 Hz) and high gamma band (65-90 Hz) in sensorimotor cortex. The performance is validated by both across-subjects statistical comparisons and by intracranial electrocorticography (ECoG) data from two epilepsy patients. Interestingly, we also reliably observed high frequency activity (30-300 Hz) in the cerebellum, although with variable locations and frequencies across subjects. The proposed algorithm is highly parallelizable and runs efficiently on modern high-performance computing clusters. This method enables the ultimate promise of MEG and EEG for five-dimensional imaging of space, time, and frequency activity in the brain and renders it applicable for widespread studies of human cortical dynamics during cognition.     

Localization of human supratemporal auditory areas from intracerebral auditory evoked potentials using distributed source models

While source localization methods are increasingly developed to identify brain areas underlying scalp electro/magnetoencephalographic data (EEG/MEG), these methods have not yet been used to identify the sources of intracerebral signals which offer highly detailed information. Here, we adapted the minimum current estimates method to intracranial data in order to localize supratemporal sources of intracerebral auditory 1-kHz-tone-evoked potentials occurring within 100 ms after stimulus onset. After an evaluation of localization method and despite inter-subject variability, we found a common spatiotemporal pattern of activities, which involved the first Heschl's gyrus (H1) and sulcus (HS), the Planum Temporale (PT), H2/H3 when present, and the superior temporal gyrus (STG). Four time periods of activity were distinguished, corresponding to the time range of the scalp components P0, Na, Pa/Pb, and N100. The sources of the earliest components P0 (16-19 ms) and Na (20-25 ms) could be identified in the postero-medial portion of HS or H1. Then, several areas became simultaneously active after 25 ms. The Pa/Pb time range (30-50 ms) was characterized by a medio-lateral and postero-anterior propagation of activity over the supratemporal plane involving successively H1/HS, the Planum Temporale, H2/H3 when present, and the STG. Finally, we found to a large extent that the N100 (55-100 ms) involved almost the same areas as those active during the Pa/Pb complex, with a similar propagation of activities. Reconstructing scalp data from these sources on fictive EEG/MEG channels reproduced classical auditory evoked waveforms and topographies. In conclusion, the spatiotemporal pattern of activation of supratemporal auditory areas could be identified on the individual anatomy using current estimates from intracerebral data. Such detailed localization approach could also be used prior to epilepsy surgery to help identify epileptogenic foci and preserve functional cortical areas.     

Recording of fast activity at the onset of partial seizures: depth EEG vs. scalp EEG

Rapid discharges (25-80 Hz), a characteristic EEG pattern often recorded at seizure onset in partial epilepsies, are often considered as electrophysiological signatures of the epileptogenic zone. While the recording of rapid discharges from intracranial electrodes has long been established, their observation from the scalp is challenging. The prevailing view is that rapid discharges cannot be seen clearly (or at all) in scalp EEG because they have low signal-to-noise ratio. To date, however, no studies have investigated the ‘observability’ of rapid discharges, i.e. under what conditions and to what extent they can be visible in recorded EEG signals. Here, we used a model-based approach to examine the impact of several factors (distance to sources, skull conductivity, source area, source synchrony, and background activity) on the observability of rapid discharges in simultaneously simulated depth EEG and scalp EEG signals. In our simulations, the rapid discharge was clearly present in depth EEG signals but mostly almost not visible in scalp EEG signals. We identified some of the factors that may limit the observability of the rapid discharge on the scalp. Notably, surrounding background activity was found to be the most critical factor. The findings are discussed in relation to the presurgical evaluation of epilepsy.     

Comparison of magnetoencephalographic spikes with and without concurrent electroencephalographic spikes in extratemporal epilepsy

Interictal spikes in patients with epilepsy may be detected by either electroencephalography (EEG) (E-spikes) or magnetoencephalography (MEG) (M-spikes), or both MEG and EEG (E/M-spikes). Localization and amplitude were compared between E/M-spikes and M-spikes in 7 adult patients with extratemporal epilepsy to evaluate the clinical significance of MEG spikes. MEG and EEG were simultaneously measured using a helmet-shaped MEG system with planar-type gradiometers and scalp electrodes of the international 10-20 system. Sources of E/M-spikes and M-spikes were estimated by an equivalent current dipole (ECD) model for MEG at peak latency. Each subject showed 9 to 20 (mean 13.4) E/M-spikes and 9 to 31 (mean 16.3) M-spikes. No subjects showed significant differences in the ECD locations between E/M- and M-spikes. ECD moments of the E/M-spikes were significantly larger in 2 patients and not significantly different in the other 5 patients. The similar localizations of E/M-spikes and M-spikes suggest that combination of MEG and EEG is useful to detect more interictal spikes in patients with extratemporal epilepsy. The smaller tendency of ECD amplitude of the M-spikes than E/M-spikes suggests that scalp EEG may overlook small tangential spikes due to background brain noise. Localization value of M-spikes is clinically equivalent to that of E/M-spikes.     

Localizing epileptogenic regions in partial epilepsy using three-dimensional statistical parametric maps of background EEG source spectra

This preliminary study sought to localize epileptogenic regions in patients with partial epilepsy by analysis of interictal EEG activity utilizing variable resolution electromagnetic tomography (VARETA), a three-dimensional quantitative electroencephalographic (QEEG) frequency-domain distributed source modeling technique. The very narrow band (VNB) spectra spanned the frequency range 0.39 Hz to 19.1 Hz, in 0.39 Hz steps. These VNB spectra were compared to normative data and transformed to provide Z-scores for every scalp derivation, and the spatial distributions of the probable EEG generators of the most abnormal values were displayed on slices from a probabilistic MRI atlas. Each voxel was color-coded to represent the significance of the deviation relative to age appropriate normative values. We compared the resulting three-dimensional images to the localization of epileptogenic regions based on invasive intracranial EEG recordings of seizure onsets. The VARETA image indicated abnormal interictal spectral power values in regions of seizure onset identified by invasive monitoring, mainly in delta and theta range (1.5 to 8.0 Hz). The VARETA localization of the most abnormal voxel was congruent with the epileptogenic regions identified by intracranial recordings with regard to hemisphere in all 6 cases, and with regard to lobe in 5 cases. In contrast, abnormal findings with routine EEG agreed with invasive monitoring with regard to hemisphere in 3 cases and with regard to lobe in 2 cases. These results suggest that analysis of background interictal EEG utilizing distributed source models should be investigated further in clinical epilepsy.     

Hierarchical Bayesian inference for the EEG inverse problem using realistic FE head models: depth localization and source separation for focal primary currents

The estimation of the activity-related ion currents by measuring the induced electromagnetic fields at the head surface is a challenging and severely ill-posed inverse problem. This is especially true in the recovery of brain networks involving deep-lying sources by means of EEG/MEG recordings which is still a challenging task for any inverse method. Recently, hierarchical Bayesian modeling (HBM) emerged as a unifying framework for current density reconstruction (CDR) approaches comprising most established methods as well as offering promising new methods. Our work examines the performance of fully-Bayesian inference methods for HBM for source configurations consisting of few, focal sources when used with realistic, high-resolution finite element (FE) head models. The main foci of interest are the correct depth localization, a well-known source of systematic error of many CDR methods, and the separation of single sources in multiple-source scenarios. Both aspects are very important in the analysis of neurophysiological data and in clinical applications. For these tasks, HBM provides a promising framework and is able to improve upon established CDR methods such as minimum norm estimation (MNE) or sLORETA in many aspects. For challenging multiple-source scenarios where the established methods show crucial errors, promising results are attained. Additionally, we introduce Wasserstein distances as performance measures for the validation of inverse methods in complex source scenarios.     

EEG source analysis and fMRI reveal two electrical sources in the fronto-parietal operculum during subepidermal finger stimulation

Using functional magnetic resonance imaging (fMRI) and electroencephalographic (EEG) source dipole analysis in 10 normal subjects, two electrical source dipoles in the contralateral fronto-parietal operculum were identified during repetitive painful subepidermal stimulation of the right index finger. The anterior source dipole peaking at 79 +/- 8 ms (mean +/- SD) was located in the frontal operculum, and oriented tangentially toward the cortical surface. The posterior source dipole peaking at 118 +/- 12 ms was located in the upper bank of the Sylvian fissure corresponding to the second somatosensory cortex (S2). The orientations of the posterior source dipoles displayed large variability, but differed significantly (P < 0.05) from the orientations of the anterior source dipoles. Electrical sources and fMRI clusters were also observed in ipsilateral fronto-parietal operculum. However, due to low signal-to-noise ratio of ipsilateral EEG sources in individual recordings, separation of sources into anterior and posterior clusters was not performed. Combined fMRI and source dipole EEG analysis of individual data suggests the presence of two distinct electrical sources in the fronto-parietal operculum participating in processing of somatosensory stimuli. The anterior region of the fronto-parietal operculum shows earlier peak activation than the posterior region.     

A realistic three dimensional FEM of the human head

A realistic three-dimensional finite element model (FEM) of the human head has been developed. Separate layers for the scalp, skull, cerebrospinal fluid (CSF) and brain were modelled. Hexahedral elements, a special master matrix assembly technique and an iterative successive over-relaxation (SOR) solution scheme were employed. This approach enabled rapid modelling with minimal memory requirements, which makes this method practical if used for electrical impedance tomography (EIT) or source localization inverse problems. Compared to scalp electrodes, subdural voltage sensing electrodes were three to four times more sensitive close to an oedema or source region, if it was peripheral, but this decreased to 30%-40% for central oedema or source regions. Scalp current injecting electrodes are preferable, since the maximum allowable current is 10 times larger than that of the subdural ones. The distance of voltage sensing electrodes from a region to be imaged highly affects sensitivity, so depth electrodes will be more sensitive, provided that they are close to the region of interest. Finally, the electrode size has significant effects on the input or transfer impedance.