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

Human brain imaging with high-density electroencephalography: Techniques and applications

Electroencephalography (EEG) is a technique for non-invasively measuring neuronal activity in the human brain using electrodes placed on the participant's scalp. With the advancement of digital technologies, EEG analysis has evolved over time from the qualitative analysis of amplitude and frequency modulations to a comprehensive analysis of the complex spatiotemporal characteristics of the recorded signals. EEG is now considered a powerful tool for measuring neural processes in the same time frame in which they happen (i.e. the subsecond range). However, it is commonly argued that EEG suffers from low spatial resolution, which makes it difficult to localize the generators of EEG activity accurately and reliably. Today, the availability of high-density EEG (hdEEG) systems, combined with methods for incorporating information on head anatomy and sophisticated source-localization algorithms, has transformed EEG into an important neuroimaging tool. hdEEG offers researchers and clinicians a rich and varied range of applications. It can be used not only for investigating neural correlates in motor and cognitive neuroscience experiments, but also for clinical diagnosis, particularly in the detection of epilepsy and the characterization of neural impairments in a wide range of neurological disorders. Notably, the integration of hdEEG systems with other physiological recordings, such as kinematic and/or electromyography data, might be especially beneficial to better understand the neuromuscular mechanisms associated with deconditioning in ageing and neuromotor disorders, by mapping the neurokinematic and neuromuscular connectivity patterns directly in the brain.

CutFEM-based MEG forward modeling improves source separability and sensitivity to quasi-radial sources: A somatosensory group study

Source analysis of magnetoencephalography (MEG) data requires the computation of the magnetic fields induced by current sources in the brain. This so-called MEG forward problem includes an accurate estimation of the volume conduction effects in the human head. Here, we introduce the Cut finite element method (CutFEM) for the MEG forward problem. CutFEM's meshing process imposes fewer restrictions on tissue anatomy than tetrahedral meshes while being able to mesh curved geometries contrary to hexahedral meshing. To evaluate the new approach, we compare CutFEM with a boundary element method (BEM) that distinguishes three tissue compartments and a 6-compartment hexahedral FEM in an n = 19 group study of somatosensory evoked fields (SEF). The neural generators of the 20 ms post-stimulus SEF components (M20) are reconstructed using both an unregularized and a regularized inversion approach. Changing the forward model resulted in reconstruction differences of about 1 centimeter in location and considerable differences in orientation. The tested 6-compartment FEM approaches significantly increase the goodness of fit to the measured data compared with the 3-compartment BEM. They also demonstrate higher quasi-radial contributions for sources below the gyral crowns. Furthermore, CutFEM improves source separability compared with both other approaches. We conclude that head models with 6 compartments rather than 3 and the new CutFEM approach are valuable additions to MEG source reconstruction, in particular for sources that are predominantly radial.

Integration of interictal EEG source localization in presurgical epilepsy evaluation - A single-center prospective study

OBJECTIVE

To investigate cost in working hours for initial integration of interictal EEG source localization (ESL) into clinical practice of a tertiary epilepsy center, and to examine concordance of results obtained with three different ESL pipelines.

METHODS

This prospective study covered the first year of using ESL in the Epilepsy-Center Berlin-Brandenburg. Patients aged ≥14 years with drug-resistant focal epilepsy referred for noninvasive presurgical evaluation were included. Interictal ESL was based on low-density EEG and individual head models. Source maxima were obtained from two freely available software packages and one commercial provider. One physician and computer scientist documented their working hours for setting up and processing ESL. Additionally, a survey was conducted among epilepsy centers in Germany to assess the current role of ESL in presurgical evaluation.

RESULTS

Of 40 patients included, 22 (55%) had enough interictal spikes for ESL. The physician's working times decreased from median 4.7 hours [interquartile range 3.9-6.4] in the first third of cases to 2.0 hours [1.9-2.4] in the remaining two thirds; P < 0.01. In addition, computer scientist and physician spent a total of 35.5 and 33.0 working hours on setting up the digital infrastructure, and on training and testing. Sublobar agreement between all three pipelines was 20%, mean measurement of agreement (kappa) 0.13. Finally, the survey revealed that 53% of epilepsy centers in Germany currently use ESL for presurgical evaluation.

SIGNIFICANCE

This study provides information regarding expected effort and costs for integration of ESL into an epilepsy surgery program. Low result agreement across different ESL pipelines calls for further standardization.

Neural stochastic differential equations network as uncertainty quantification method for EEG source localization

EEG source localization remains a challenging problem given the uncertain conductivity values of the volume conductor models (VCMs). As uncertain conductivities vary across people, they may considerably impact the forward and inverse solutions of the EEG, leading to an increase in localization mistakes and misdiagnoses of brain disorders. Calibration of conductivity values using uncertainty quantification (UQ) techniques is a promising approach to reduce localization errors. The widely-known UQ methods involve Bayesian approaches, which utilize prior conductivity values to derive their posterior inference and estimate their optimal calibration. However, these approaches have two significant drawbacks: solving for posterior inference is intractable, and choosing inappropriate priors may lead to increased localization mistakes. This study used the Neural Stochastic Differential equations Network (SDE-Net), a combination of dynamical systems and deep learning techniques that utilizes the Wiener process to minimize conductivity uncertainties in the VCM and improve the inverse problem. Results revealed that SDE-Net generated a lower localization error rate in the inverse problem compared to Bayesian techniques. Future studies may employ new stochastic dynamical systems-based techniques as a UQ technique to address further uncertainties in the EEG Source Localization problem. Our code can be found here:https://github.com/rrwabina/SDENet-UQ-ESL.

Automated interictal source localisation based on high-density EEG

PURPOSE

To study the accuracy of automated interictal EEG source localisation based on high-density EEG, and to compare it to low-density EEG.

METHODS

Thirty patients operated for pharmacoresistant focal epilepsy were retrospectively examined. Twelve months after resective brain surgery, 18 were seizure-free or had 'auras' only, while 12 had persistence of disabling seizures. Presurgical 257-channel EEG lasting 3-20 h was down-sampled to 25, 40, and 204 channels for separate analyses. For each electrode setup, interictal spikes were detected, clustered, and averaged automatically before validation by an expert reviewer. An individual 6-layer finite difference head model and the standardised low-resolution electromagnetic tomography were used to localise the maximum source activity of the most prevalent spike. Sublobar concordance with the resected brain area was visually assessed and related to favourable vs. unfavourable postsurgical outcome.

RESULTS

Depending on the EEG setup, epileptic spikes were detected in 21-24 patients (70-80%). The median number of single spikes per average was 470 (range 17-15,066). Diagnostic sensitivity of EEG source localisation was 58-75%, specificity was 50-67%, and overall accuracy was 55-71%. There were no significant differences between low- and high-density EEG setups with 25 to 257 electrodes.

CONCLUSION

Automated high-density EEG source localisation provides meaningful information in the majority of cases. With hundreds of single spikes averaged, diagnostic accuracy is similar in high- and low-density EEG. Therefore, low-density EEG may be sufficient for interictal EEG source localisation if high numbers of spikes are available.

Automated ictal EEG source imaging: A retrospective, blinded clinical validation study

OBJECTIVE

EEG source imaging (ESI) is a validated tool in the multimodal workup of patients with drug resistant focal epilepsy. However, it requires special expertise and it is underutilized. To circumvent this, automated analysis pipelines have been developed and validated for the interictal discharges. In this study, we present the clinical validation of an automated ESI for ictal EEG signals.

METHODS

We have developed an automated analysis pipeline of ictal EEG activity, based on spectral analysis in source space, using an individual head model of six tissues. The analysis was done blinded to all other data. As reference standard, we used the concordance with the resected area and one-year postoperative outcome.

RESULTS

We analyzed 50 consecutive patients undergoing epilepsy surgery (34 temporal and 16 extra-temporal). Thirty patients (60%) became seizure-free. The accuracy of the automated ESI was 74% (95% confidence interval: 59.66-85.37%).

CONCLUSIONS

Automated ictal ESI has a high accuracy for localizing the seizure onset zone.

SIGNIFICANCE

Automating the ESI of the ictal EEG signals will facilitate implementation of this tool in the presurgical evaluation.

Automated Head Tissue Modelling Based on Structural Magnetic Resonance Images for Electroencephalographic Source Reconstruction

In the last years, technological advancements for the analysis of electroencephalography (EEG) recordings have permitted to investigate neural activity and connectivity in the human brain with unprecedented precision and reliability. A crucial element for accurate EEG source reconstruction is the construction of a realistic head model, incorporating information on electrode positions and head tissue distribution. In this paper, we introduce MR-TIM, a toolbox for head tissue modelling from structural magnetic resonance (MR) images. The toolbox consists of three modules: 1) image pre-processing – the raw MR image is denoised and prepared for further analyses; 2) tissue probability mapping – template tissue probability maps (TPMs) in individual space are generated from the MR image; 3) tissue segmentation – information from all the TPMs is integrated such that each voxel in the MR image is assigned to a specific tissue. MR-TIM generates highly realistic 3D masks, five of which are associated with brain structures (brain and cerebellar grey matter, brain and cerebellar white matter, and brainstem) and the remaining seven with other head tissues (cerebrospinal fluid, spongy and compact bones, eyes, muscle, fat and skin). Our validation, conducted on MR images collected in healthy volunteers and patients as well as an MR template image from an open-source repository, demonstrates that MR-TIM is more accurate than alternative approaches for whole-head tissue segmentation. We hope that MR-TIM, by yielding an increased precision in head modelling, will contribute to a more widespread use of EEG as a brain imaging technique.

Brain Waves: Emergence of Localized, Persistent, Weakly Evanescent Cortical Loops

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of nonevanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full-brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the commonly agreed paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation because of the cortex wave modes in the highly folded areas will exhibit no apparent correlation with the fiber directions. Nonlinear coupling of those linear weakly evanescent wave modes then provides a universal mechanism for the emergence of synchronized brain wave field activity. The resonant and nonresonant terms of nonlinear coupling between multiple modes produce both synchronous spiking-like high-frequency wave activity as well as low-frequency wave rhythms. Numerical simulation of forced multiple-mode dynamics shows that, as forcing increases, there is a transition from damped to oscillatory regime that can then transition quickly to a nonoscillatory state when a critical excitation threshold is reached. The resonant nonlinear coupling results in the emergence of low-frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these weakly evanescent cortical wave modes have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and magnetoencephalography, and for the understanding of brain activity in general, including mechanisms of memory.

Influence of Patient-Specific Head Modeling on EEG Source Imaging

Electromagnetic source imaging (ESI) techniques have become one of the most common alternatives for understanding cognitive processes in the human brain and for guiding possible therapies for neurological diseases. However, ESI accuracy strongly depends on the forward model capabilities to accurately describe the subject's head anatomy from the available structural data. Attempting to improve the ESI performance, we enhance the brain structure model within the individual-defined forward problem formulation, combining the head geometry complexity of the modeled tissue compartments and the prior knowledge of the brain tissue morphology. We validate the proposed methodology using 25 subjects, from which a set of magnetic-resonance imaging scans is acquired, extracting the anatomical priors and an electroencephalography signal set needed for validating the ESI scenarios. Obtained results confirm that incorporating patient-specific head models enhances the performed accuracy and improves the localization of focal and deep sources.

Universal theory of brain waves: from linear loops to nonlinear synchronized spiking and collective brain rhythms

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of non-evanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and to the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the standard paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation due to the cortex wave modes in highly folded areas will exhibit no apparent correlation with the fiber directions. The results are consistent with numerous recent experimental animal and human brain studies demonstrating the existence of electrostatic field activity in the form of traveling waves (including studies where neuronal connections were severed) and with wave loop induced peaks observed in EEG spectra. In addition, we demonstrate that the resonant and non-resonant terms of the nonlinear coupling between multiple modes produce both synchronous spiking-like high frequency wave activity as well as low frequency wave rhythms as a result of their unique dispersion properties. Numerical simulation of forced multiple mode dynamics shows that as forcing increases there is a transition from damped to oscillatory regime that subsequently decays away as over-excitation is reached. The resonant nonlinear coupling results in the emergence of low frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these cortical wave modes, and this new mechanism for understanding the nature of spiking behavior, have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and MEG, and in general for the understanding of brain activity, including mechanisms of memory.

Unification of optimal targeting methods in transcranial electrical stimulation

One of the major questions in high-density transcranial electrical stimulation (TES) is: given a region of interest (ROI) and electric current limits for safety, how much current should be delivered by each electrode for optimal targeting of the ROI? Several solutions, apparently unrelated, have been independently proposed depending on how "optimality" is defined and on how this optimization problem is stated mathematically. The least squares (LS), weighted LS (WLS), or reciprocity-based approaches are the simplest ones and have closed-form solutions. An extended optimization problem can be stated as follows: maximize the directional intensity at the ROI, limit the electric fields at the non-ROI, and constrain total injected current and current per electrode for safety. This problem requires iterative convex or linear optimization solvers. We theoretically prove in this work that the LS, WLS and reciprocity-based closed-form solutions are specific solutions to the extended directional maximization optimization problem. Moreover, the LS/WLS and reciprocity-based solutions are the two extreme cases of the intensity-focality trade-off, emerging under variation of a unique parameter of the extended directional maximization problem, the imposed constraint to the electric fields at the non-ROI. We validate and illustrate these findings with simulations on an atlas head model. The unified approach we present here allows a better understanding of the nature of the TES optimization problem and helps in the development of advanced and more effective targeting strategies.

Multi-Resolution Graph Based Volumetric Cortical Basis Functions From Local Anatomic Features

OBJECTIVE

Modern clinical MRI collects millimeter scale anatomic information, but scalp electroencephalography source localization is ill posed, and cannot resolve individual sources at that resolution. Dimensionality reduction in the space of cortical sources is needed to improve computational and storage complexity, yet volumetric methods still employ simplistic grid coarsening that eliminates fine scale anatomic structure. We present an approach to extend near-arbitrary spatial scaling to volumetric localization.

METHODS

Starting from a voxelwise brain parcellation, sub-parcels are identified from local cortical connectivity with an iterated graph cut approach. Spatial basis functions in each parcel are constructed using either a decomposition of the local leadfield matrix or spectral basis functions of local cortical connectivity graphs.

RESULTS

We present quantitative evaluation with extensive simulations and use multiple sets of real data to highlight how parameter changes impact computed reconstructions. Our results show that volumetric basis functions can improve accuracy by as much as 30%, while reducing computational complexity by over two orders of magnitude. In real data from epilepsy surgical candidates, accurate localization of seizure onset regions is demonstrated.

CONCLUSION

Spatial dimensionality reduction with volumetric basis functions improves reconstruction accuracy while reducing computational complexity.

SIGNIFICANCE

Near-arbitrary spatial dimensionality reduction will enable volumetric reconstruction with modern computationally intensive algorithms and anatomically driven multi-resolution methods.

Variation in Reported Human Head Tissue Electrical Conductivity Values

Electromagnetic source characterisation requires accurate volume conductor models representing head geometry and the electrical conductivity field. Head tissue conductivity is often assumed from previous literature, however, despite extensive research, measurements are inconsistent. A meta-analysis of reported human head electrical conductivity values was therefore conducted to determine significant variation and subsequent influential factors. Of 3121 identified publications spanning three databases, 56 papers were included in data extraction. Conductivity values were categorised according to tissue type, and recorded alongside methodology, measurement condition, current frequency, tissue temperature, participant pathology and age. We found variation in electrical conductivity of the whole-skull, the spongiform layer of the skull, isotropic, perpendicularly- and parallelly-oriented white matter (WM) and the brain-to-skull-conductivity ratio (BSCR) could be significantly attributed to a combination of differences in methodology and demographics. This large variation should be acknowledged, and care should be taken when creating volume conductor models, ideally constructing them on an individual basis, rather than assuming them from the literature. When personalised models are unavailable, it is suggested weighted average means from the current meta-analysis are used. Assigning conductivity as: 0.41 S/m for the scalp, 0.02 S/m for the whole skull, or when better modelled as a three-layer skull 0.048 S/m for the spongiform layer, 0.007 S/m for the inner compact and 0.005 S/m for the outer compact, as well as 1.71 S/m for the CSF, 0.47 S/m for the grey matter, 0.22 S/m for WM and 50.4 for the BSCR.

A Finite-Difference Solution for the EEG Forward Problem in Inhomogeneous Anisotropic Media

Accurate source localization of electroencephalographic (EEG) signals requires detailed information about the geometry and physical properties of head tissues. Indeed, these strongly influence the propagation of neural activity from the brain to the sensors. Finite difference methods (FDMs) are head modelling approaches relying on volumetric data information, which can be directly obtained using magnetic resonance (MR) imaging. The specific goal of this study is to develop a computationally efficient FDM solution that can flexibly integrate voxel-wise conductivity and anisotropy information. Given the high computational complexity of FDMs, we pay particular attention to attain a very low numerical error, as evaluated using exact analytical solutions for spherical volume conductor models. We then demonstrate the computational efficiency of our FDM numerical solver, by comparing it with alternative solutions. Finally, we apply the developed head modelling tool to high-resolution MR images from a real experimental subject, to demonstrate the potential added value of incorporating detailed voxel-wise conductivity and anisotropy information. Our results clearly show that the developed FDM can contribute to a more precise head modelling, and therefore to a more reliable use of EEG as a brain imaging tool.

Seizure Onset Zone Localization from Ictal High-Density EEG in Refractory Focal Epilepsy

Epilepsy surgery is the most efficient treatment option for patients with refractory epilepsy. Before surgery, it is of utmost importance to accurately delineate the seizure onset zone (SOZ). Non-invasive EEG is the most used neuroimaging technique to diagnose epilepsy, but it is hard to localize the SOZ from EEG due to its low spatial resolution and because epilepsy is a network disease, with several brain regions becoming active during a seizure. In this work, we propose and validate an approach based on EEG source imaging (ESI) combined with functional connectivity analysis to overcome these problems. We considered both simulations and real data of patients. Ictal epochs of 204-channel EEG and subsets down to 32 channels were analyzed. ESI was done using realistic head models and LORETA was used as inverse technique. The connectivity pattern between the reconstructed sources was calculated, and the source with the highest number of outgoing connections was selected as SOZ. We compared this algorithm with a more straightforward approach, i.e. selecting the source with the highest power after ESI as the SOZ. We found that functional connectivity analysis estimated the SOZ consistently closer to the simulated EZ/RZ than localization based on maximal power. Performance, however, decreased when 128 electrodes or less were used, especially in the realistic data. The results show the added value of functional connectivity analysis for SOZ localization, when the EEG is obtained with a high-density setup. Next to this, the method can potentially be used as objective tool in clinical settings.

Using reciprocity for relating the simulation of transcranial current stimulation to the EEG forward problem

To explore the relationship between transcranial current stimulation (tCS) and the electroencephalography (EEG) forward problem, we investigate and compare accuracy and efficiency of a reciprocal and a direct EEG forward approach for dipolar primary current sources both based on the finite element method (FEM), namely the adjoint approach (AA) and the partial integration approach in conjunction with a transfer matrix concept (PI). By analyzing numerical results, comparing to analytically derived EEG forward potentials and estimating computational complexity in spherical shell models, AA turns out to be essentially identical to PI. It is then proven that AA and PI are also algebraically identical even for general head models. This relation offers a direct link between the EEG forward problem and tCS. We then demonstrate how the quasi-analytical EEG forward solutions in sphere models can be used to validate the numerical accuracies of FEM-based tCS simulation approaches. These approaches differ with respect to the ease with which they can be employed for realistic head modeling based on MRI-derived segmentations. We show that while the accuracy of the most easy to realize approach based on regular hexahedral elements is already quite high, it can be significantly improved if a geometry-adaptation of the elements is employed in conjunction with an isoparametric FEM approach. While the latter approach does not involve any additional difficulties for the user, it reaches the high accuracies of surface-segmentation based tetrahedral FEM, which is considerably more difficult to implement and topologically less flexible in practice. Finally, in a highly realistic head volume conductor model and when compared to the regular alternative, the geometry-adapted hexahedral FEM is shown to result in significant changes in tCS current flow orientation and magnitude up to 45° and a factor of 1.66, respectively.

The Role of Skull Modeling in EEG Source Imaging for Patients with Refractory Temporal Lobe Epilepsy

We investigated the influence of different skull modeling approaches on EEG source imaging (ESI), using data of six patients with refractory temporal lobe epilepsy who later underwent successful epilepsy surgery. Four realistic head models with different skull compartments, based on finite difference methods, were constructed for each patient: (i) Three models had skulls with compact and spongy bone compartments as well as air-filled cavities, segmented from either computed tomography (CT), magnetic resonance imaging (MRI) or a CT-template and (ii) one model included a MRI-based skull with a single compact bone compartment. In all patients we performed ESI of single and averaged spikes marked in the clinical 27-channel EEG by the epileptologist. To analyze at which time point the dipole estimations were closer to the resected zone, ESI was performed at two time instants: the half-rising phase and peak of the spike. The estimated sources for each model were validated against the resected area, as indicated by the postoperative MRI. Our results showed that single spike analysis was highly influenced by the signal-to-noise ratio (SNR), yielding estimations with smaller distances to the resected volume at the peak of the spike. Although averaging reduced the SNR effects, it did not always result in dipole estimations lying closer to the resection. The proposed skull modeling approaches did not lead to significant differences in the localization of the irritative zone from clinical EEG data with low spatial sampling density. Furthermore, we showed that a simple skull model (MRI-based) resulted in similar accuracy in dipole estimation compared to more complex head models (based on CT- or CT-template). Therefore, all the considered head models can be used in the presurgical evaluation of patients with temporal lobe epilepsy to localize the irritative zone from low-density clinical EEG recordings.

The role of blood vessels in high-resolution volume conductor head modeling of EEG

Reconstruction of the electrical sources of human EEG activity at high spatio-temporal accuracy is an important aim in neuroscience and neurological diagnostics. Over the last decades, numerous studies have demonstrated that realistic modeling of head anatomy improves the accuracy of source reconstruction of EEG signals. For example, including a cerebro-spinal fluid compartment and the anisotropy of white matter electrical conductivity were both shown to significantly reduce modeling errors. Here, we for the first time quantify the role of detailed reconstructions of the cerebral blood vessels in volume conductor head modeling for EEG. To study the role of the highly arborized cerebral blood vessels, we created a submillimeter head model based on ultra-high-field-strength (7T) structural MRI datasets. Blood vessels (arteries and emissary/intraosseous veins) were segmented using Frangi multi-scale vesselness filtering. The final head model consisted of a geometry-adapted cubic mesh with over 17×10(6) nodes. We solved the forward model using a finite-element-method (FEM) transfer matrix approach, which allowed reducing computation times substantially and quantified the importance of the blood vessel compartment by computing forward and inverse errors resulting from ignoring the blood vessels. Our results show that ignoring emissary veins piercing the skull leads to focal localization errors of approx. 5 to 15mm. Large errors (>2cm) were observed due to the carotid arteries and the dense arterial vasculature in areas such as in the insula or in the medial temporal lobe. Thus, in such predisposed areas, errors caused by neglecting blood vessels can reach similar magnitudes as those previously reported for neglecting white matter anisotropy, the CSF or the dura - structures which are generally considered important components of realistic EEG head models. Our findings thus imply that including a realistic blood vessel compartment in EEG head models will be helpful to improve the accuracy of EEG source analyses particularly when high accuracies in brain areas with dense vasculature are required.

Electrical source imaging of interictal spikes using multiple sparse volumetric priors for presurgical epileptogenic focus localization

Electrical source imaging of interictal spikes observed in EEG recordings of patients with refractory epilepsy provides useful information to localize the epileptogenic focus during the presurgical evaluation. However, the selection of the time points or time epochs of the spikes in order to estimate the origin of the activity remains a challenge. In this study, we consider a Bayesian EEG source imaging technique for distributed sources, i.e. the multiple volumetric sparse priors (MSVP) approach. The approach allows to estimate the time courses of the intensity of the sources corresponding with a specific time epoch of the spike. Based on presurgical averaged interictal spikes in six patients who were successfully treated with surgery, we estimated the time courses of the source intensities for three different time epochs: (i) an epoch starting 50 ms before the spike peak and ending at 50% of the spike peak during the rising phase of the spike, (ii) an epoch starting 50 ms before the spike peak and ending at the spike peak and (iii) an epoch containing the full spike time period starting 50 ms before the spike peak and ending 230 ms after the spike peak. To identify the primary source of the spike activity, the source with the maximum energy from 50 ms before the spike peak till 50% of the spike peak was subsequently selected for each of the time windows. For comparison, the activity at the spike peaks and at 50% of the peaks was localized using the LORETA inversion technique and an ECD approach. Both patient-specific spherical forward models and patient-specific 5-layered finite difference models were considered to evaluate the influence of the forward model. Based on the resected zones in each of the patients, extracted from post-operative MR images, we compared the distances to the resection border of the estimated activity. Using the spherical models, the distances to the resection border for the MSVP approach and each of the different time epochs were in the same range as the LORETA and ECD techniques. We found distances smaller than 23 mm, with robust results for all the patients. For the finite difference models, we found that the distances to the resection border for the MSVP inversions of the full spike time epochs were generally smaller compared to the MSVP inversions of the time epochs before the spike peak. The results also suggest that the inversions using the finite difference models resulted in slightly smaller distances to the resection border compared to the spherical models. The results we obtained are promising because the MSVP approach allows to study the network of the estimated source-intensities and allows to characterize the spatial extent of the underlying sources.

•

A Bayesian ESI technique is evaluated to localize interictal spike activity.

•

Averaged spikes in six patients were used that were seizure free after surgery.

•

We compared the technique with the LORETA an ECD technique.

•

We evaluated both spherical and 5-layered finite difference forward models.

•

Our approach is potentially useful to delineate the irritative zone.

Influence of the head model on EEG and MEG source connectivity analyses

The results of brain connectivity analysis using reconstructed source time courses derived from EEG and MEG data depend on a number of algorithmic choices. While previous studies have investigated the influence of the choice of source estimation method or connectivity measure, the effects of the head modeling errors or simplifications have not been studied sufficiently. In the present simulation study, we investigated the influence of particular properties of the head model on the reconstructed source time courses as well as on source connectivity analysis in EEG and MEG. Therefore, we constructed a realistic head model and applied the finite element method to solve the EEG and MEG forward problems. We considered the distinction between white and gray matter, the distinction between compact and spongy bone, the inclusion of a cerebrospinal fluid (CSF) compartment, and the reduction to a simple 3-layer model comprising only the skin, skull, and brain. Source time courses were reconstructed using a beamforming approach and the source connectivity was estimated by the imaginary coherence (ICoh) and the generalized partial directed coherence (GPDC). Our results show that in both EEG and MEG, neglecting the white and gray matter distinction or the CSF causes considerable errors in reconstructed source time courses and connectivity analysis, while the distinction between spongy and compact bone is just of minor relevance, provided that an adequate skull conductivity value is used. Large inverse and connectivity errors are found in the same regions that show large topography errors in the forward solution. Moreover, we demonstrate that the very conservative ICoh is relatively safe from the crosstalk effects caused by imperfect head models, as opposed to the GPDC.

Frequency-pattern functional tomography of magnetoencephalography data allows new approach to the study of human brain organization

A method based on a set of new theorems for the analysis of multichannel time series is described, based on precise Fourier transform and coherence analysis of the restored signals from a detailed set of frequency components. Magnetic field recordings of spontaneous and evoked activity by means of magnetic encephalography demonstrated that multichannel precise Fourier spectrum contains a very large set of harmonics with high coherence. The inverse problem can be solved with great precision based on coherent harmonics, so the technique is a promising platform of general analysis in brain imaging. The analysis method makes it possible to reconstruct sites and timing of electrical activity generated by both spontaneous and evoked brain function at different depths in the brain in the millisecond time range.

Influence of isotropic skull models on EEG source localization

We investigated the influence of using simplified skull models on electroencephalogram (EEG) source localization. The simplified skull models were derived from CT and MR images, with isotropic conductivity modeled as either heterogeneous or homogeneous. A total of four simplified head models were compared against a reference model with a skull accurately segmented with CT images. Our results show that the use of a simplified geometry for the skull, can lead to errors of approximately 1 cm for sources located in the central and temporal regions of the brain.

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.

Influence of skull modeling approaches on EEG source localization

Electroencephalographic source localization (ESL) relies on an accurate model representing the human head for the computation of the forward solution. In this head model, the skull is of utmost importance due to its complex geometry and low conductivity compared to the other tissues inside the head. We investigated the influence of using different skull modeling approaches on ESL. These approaches, consisting in skull conductivity and geometry modeling simplifications, make use of X-ray computed tomography (CT) and magnetic resonance (MR) images to generate seven different head models. A head model with an accurately segmented skull from CT images, including spongy and compact bone compartments as well as some air-filled cavities, was used as the reference model. EEG simulations were performed for a configuration of 32 and 128 electrodes, and for both noiseless and noisy data. The results show that skull geometry simplifications have a larger effect on ESL than those of the conductivity modeling. This suggests that accurate skull modeling is important in order to achieve reliable results for ESL that are useful in a clinical environment. We recommend the following guidelines to be taken into account for skull modeling in the generation of subject-specific head models: (i) If CT images are available, i.e., if the geometry of the skull and its different tissue types can be accurately segmented, the conductivity should be modeled as isotropic heterogeneous. The spongy bone might be segmented as an erosion of the compact bone; (ii) when only MR images are available, the skull base should be represented as accurately as possible and the conductivity can be modeled as isotropic heterogeneous, segmenting the spongy bone directly from the MR image; (iii) a large number of EEG electrodes should be used to obtain high spatial sampling, which reduces the localization errors at realistic noise levels.

Effects of sutures and fontanels on MEG and EEG source analysis in a realistic infant head model

In infants, the fontanels and sutures as well as conductivity of the skull influence the volume currents accompanying primary currents generated by active neurons and thus the associated electroencephalography (EEG) and magnetoencephalography (MEG) signals. We used a finite element method (FEM) to construct a realistic model of the head of an infant based on MRI images. Using this model, we investigated the effects of the fontanels, sutures and skull conductivity on forward and inverse EEG and MEG source analysis. Simulation results show that MEG is better suited than EEG to study early brain development because it is much less sensitive than EEG to distortions of the volume current caused by the fontanels and sutures and to inaccurate estimates of skull conductivity. Best results will be achieved when MEG and EEG are used in combination.

Subspace electrode selection methodology for the reduction of the effect of uncertain conductivity values in the EEG dipole localization: a simulation study using a patient-specific head model

The simulation of the electroencephalogram (EEG) using a realistic head model needs the correct conductivity values of several tissues. However, these values are not precisely known and have an influence on the accuracy of the EEG source analysis problem. This paper presents a novel numerical methodology for the increase of accuracy of the EEG dipole source localization problem. The presented subspace electrode selection (SES) methodology is able to limit the effect of uncertain conductivity values on the solution of the EEG inverse problem, yielding improved source localization accuracy. We redefine the traditional cost function associated with the EEG inverse problem and introduce a selection procedure of EEG potentials. In each iteration of the presented EEG cost function minimization procedure, potentials are selected that are affected as little as possible by the uncertain conductivity value. Using simulation data, the proposed SES methodology is able to enhance the neural source localization accuracy dependent on the dipole location, the assumed versus actual conductivity and the possible noise in measurements.

Evaluation of multiple-sphere head models for MEG source localization

Magnetoencephalography (MEG) source analysis has largely relied on spherical conductor models of the head to simplify forward calculations of the brain's magnetic field. Multiple- (or overlapping, local) sphere models, where an optimal sphere is selected for each sensor, are considered an improvement over single-sphere models and are computationally simpler than realistic models. However, there is limited information available regarding the different methods used to generate these models and their relative accuracy. We describe a variety of single- and multiple-sphere fitting approaches, including a novel method that attempts to minimize the field error. An accurate boundary element method simulation was used to evaluate the relative field measurement error (12% on average) and dipole fit localization bias (3.5 mm) of each model over the entire brain. All spherical models can contribute in the order of 1 cm to the localization bias in regions of the head that depart significantly from a sphere (inferior frontal and temporal). These spherical approximation errors can give rise to larger localization differences when all modeling effects are taken into account and with more complex source configurations or other inverse techniques, as shown with a beamformer example. Results differed noticeably depending on the source location, making it difficult to recommend a fitting method that performs best in general. Given these limitations, it may be advisable to expand the use of realistic head models.

Influence of skull conductivity perturbations on EEG dipole source analysis

PURPOSE

Electroencephalogram (EEG) source analysis is a noninvasive technique used in the presurgical of epilepsy. In this study, the dipole location and orientation errors due to skull conductivity perturbations were investigated in two groups of three-dimensional head models: A spherical head model and a realistic head model.

METHODS

In each group, the head model had a brain-to-skull conductivity ratio (Rsigma) within the range of 10-40. Solving the forward problem in the head model with skull conductivity perturbations along with the inverse problem in the baseline head model with Rsigma=20 permitted the derivation of the dipole estimation errors.

RESULTS

Perturbations in the skull conductivity generated dipole location and orientation errors: The larger the perturbations, the larger the errors and the error ranges. The dipole orientation error due to skull conductivity perturbations was not great (maximal mean of < 5 degrees) in this study, but the dipole location error was notable (maximal mean of > 5 mm).

CONCLUSIONS

Therefore, the influence of skull conductivity perturbations on EEG dipole source analysis cannot be neglected. This study suggests that it is necessary to measure the skull conductivity of the individual patients in order to achieve accurate EEG source analysis.

Uncertainty and sensitivity analysis for anisotropic inhomogeneous head tissue conductivity in human head modelling

The accuracy of an electroencephalography (EEG) forward problem partially depends on the head tissue conductivities. These conductivities are anisotropic and inhomogeneous in nature. This paper investigates the effects of conductivity uncertainty and analyses its sensitivity on an EEG forward problem for a spherical and a realistic head models. We estimate the uncertain conductivities using an efficient constraint based on an optimization method and perturb it by means of the volume and directional constraints. Assigning the uncertain conductivities, we construct spherical and realistic head models by means of a stochastic finite element method for fixed dipolar sources. We also compute EEG based on the constructed head models. We use a probabilistic sensitivity analysis method to determine the sensitivity indexes. These indexes characterize the conductivities with the most or the least effects on the computed outputs. These results demonstrate that conductivity uncertainty has significant effects on EEG. These results also show that the uncertain conductivities of the scalp, the radial direction of the skull and transversal direction in the white matter are more sensible.

Sensitivity of beamformer source analysis to deficiencies in forward modeling

Beamforming approaches have recently been developed for the field of electroencephalography (EEG) and magnetoencephalography (MEG) source analysis and opened up new applications within various fields of neuroscience. While the number of beamformer applications thus increases fast-paced, fundamental methodological considerations, especially the dependence of beamformer performance on leadfield accuracy, is still quite unclear. In this article, we present a systematic study on the influence of improper volume conductor modeling on the source reconstruction performance of an EEG-data based synthetic aperture magnetometry (SAM) beamforming approach. A finite element model of a human head is derived from multimodal MR images and serves as a realistic volume conductor model. By means of a theoretical analysis followed by a series of computer simulations insight is gained into beamformer performance with respect to reconstruction errors in peak location, peak amplitude, and peak width resulting from geometry and anisotropy volume conductor misspecifications, sensor noise, and insufficient sensor coverage. We conclude that depending on source position, sensor coverage, and accuracy of the volume conductor model, localization errors up to several centimeters must be expected. As we could show that the beamformer tries to find the best fitting leadfield (least squares) with respect to its scanning space, this result can be generalized to other localization methods. More specific, amplitude, and width of the beamformer peaks significantly depend on the interaction between noise and accuracy of the volume conductor model. The beamformer can strongly profit from a high signal-to-noise ratio, but this requires a sufficiently realistic volume conductor model.

A systematic study of head tissue inhomogeneity and anisotropy on EEG forward problem computing

In this study, we propose a stochastic method to analyze the effects of inhomogeneous anisotropic tissue conductivity on electroencephalogram (EEG) in forward computation. We apply this method to an inhomogeneous and anisotropic spherical human head model. We apply stochastic finite element method based on Legendre polynomials, Karhunen-Loeve expansion and stochastic Galerkin methods. We apply Volume and Wang's constraints to restrict the anisotropic conductivities for both the white matter (WM) and the skull tissue compartments. The EEGs resulting from deterministic and stochastic FEMs are compared using statistical measurement techniques. Based on these comparisons, we find that EEGs generated by incorporating WM and skull inhomogeneous anisotropic tissue properties individually result in an average of 56.5 and 57.5% relative errors, respectively. Incorporating these tissue properties for both layers together generate 43.5% average relative error. Inhomogeneous scalp tissue causes 27% average relative error and a full inhomogeneous anisotropic model brings in an average of 45.5% relative error. The study results demonstrate that the effects of inhomogeneous anisotropic tissue conductivity are significant on EEG.

Comparison between realistic and spherical approaches in EEG forward modelling

In electroencephalography (EEG) a valid conductor model of the head (forward model) is necessary for predicting measurable scalp voltages from intra-cranial current distributions. All inverse models, capable of inferring the spatial distribution of the neural sources generating measurable electrical and magnetic signals outside the brain are normally formulated in terms of a pre-estimated forward model, which implies considering one (or more) current dipole(s) inside the head and computing the electrical potentials generated at the electrode sites on the scalp surface. Therefore, the accuracy of the forward model strongly affects the reliability of the source reconstruction process independently of the specific inverse model. So far, it is as yet unclear which brain regions are more sensitive to the choice of different model geometry, from both quantitative and qualitative points of view. In this paper, we compare the finite difference method-based realistic model with the four-layers sensor-fitted spherical model using simulated cortical sources in the MNI152 standard space. We focused on the investigation of the spatial variation of the lead fields produced by simulated cortical sources which were placed on the reconstructed mesh of the neocortex along the surface electrodes of a 62-channel configuration. This comparison is carried out by evaluating a point spread function all over the brain cortex, with the aim of finding the lead fields mismatch between realistic and spherical geometry. Realistic geometry turns out to be a relevant factor of improvement which is particularly important when considering sources placed in the temporal or in the occipital cortex. In these situations, using a realistic head model will allow a better spatial discrimination of neural sources when compared to the spherical model.

Current themes in neuroimaging of epilepsy: brain networks, dynamic phenomena, and clinical relevance

Brain scanning methods were first applied in patients with epilepsy more than 30years ago. A very substantial literature now exists in this field, which is exponentially increasing. Contemporary neuroimaging studies in epilepsy reflect new concepts in the epilepsies, as well as current methodological developments. In particular, this area is emphasising the role of networks in epileptogenicity, the existence of dynamic phenomena which can be captured by imaging, and is beginning to validate the implementation of neuroimaging in the clinic. Here, recent studies of the last 5years are reviewed, covering the full range of neuroimaging methods with SPECT, PET and MRI in epilepsy.

Improved EEG source analysis using low-resolution conductivity estimation in a four-compartment finite element head model

Bioelectric source analysis in the human brain from scalp electroencephalography (EEG) signals is sensitive to geometry and conductivity properties of the different head tissues. We propose a low-resolution conductivity estimation (LRCE) method using simulated annealing optimization on high-resolution finite element models that individually optimizes a realistically shaped four-layer volume conductor with regard to the brain and skull compartment conductivities. As input data, the method needs T1- and PD-weighted magnetic resonance images for an improved modeling of the skull and the cerebrospinal fluid compartment and evoked potential data with high signal-to-noise ratio (SNR). Our simulation studies showed that for EEG data with realistic SNR, the LRCE method was able to simultaneously reconstruct both the brain and the skull conductivity together with the underlying dipole source and provided an improved source analysis result. We have also demonstrated the feasibility and applicability of the new method to simultaneously estimate brain and skull conductivity and a somatosensory source from measured tactile somatosensory-evoked potentials of a human subject. Our results show the viability of an approach that computes its own conductivity values and thus reduces the dependence on assigning values from the literature and likely produces a more robust estimate of current sources. Using the LRCE method, the individually optimized four-compartment volume conductor model can, in a second step, be used for the analysis of clinical or cognitive data acquired from the same subject.

Dipole estimation errors due to not incorporating anisotropic conductivities in realistic head models for EEG source analysis

EEG source analysis is a valuable tool for brain functionality research and for diagnosing neurological disorders, such as epilepsy. It requires a geometrical representation of the human head or a head model, which is often modeled as an isotropic conductor. However, it is known that some brain tissues, such as the skull or white matter, have an anisotropic conductivity. Many studies reported that the anisotropic conductivities have an influence on the calculated electrode potentials. However, few studies have assessed the influence of anisotropic conductivities on the dipole estimations. In this study, we want to determine the dipole estimation errors due to not taking into account the anisotropic conductivities of the skull and/or brain tissues. Therefore, head models are constructed with the same geometry, but with an anisotropically conducting skull and/or brain tissue compartment. These head models are used in simulation studies where the dipole location and orientation error is calculated due to neglecting anisotropic conductivities of the skull and brain tissue. Results show that not taking into account the anisotropic conductivities of the skull yields a dipole location error between 2 and 25 mm, with an average of 10 mm. When the anisotropic conductivities of the brain tissues are neglected, the dipole location error ranges between 0 and 5 mm. In this case, the average dipole location error was 2.3 mm. In all simulations, the dipole orientation error was smaller than 10 degrees . We can conclude that the anisotropic conductivities of the skull have to be incorporated to improve the accuracy of EEG source analysis. The results of the simulation, as presented here, also suggest that incorporation of the anisotropic conductivities of brain tissues is not necessary. However, more studies are needed to confirm these suggestions.

Accuracy and run-time comparison for different potential approaches and iterative solvers in finite element method based EEG source analysis

Accuracy and run-time play an important role in medical diagnostics and research as well as in the field of neuroscience. In Electroencephalography (EEG) source reconstruction, a current distribution in the human brain is reconstructed noninvasively from measured potentials at the head surface (the EEG inverse problem). Numerical modeling techniques are used to simulate head surface potentials for dipolar current sources in the human cortex, the so-called EEG forward problem.In this paper, the efficiency of algebraic multigrid (AMG), incomplete Cholesky (IC) and Jacobi preconditioners for the conjugate gradient (CG) method are compared for iteratively solving the finite element (FE) method based EEG forward problem. The interplay of the three solvers with a full subtraction approach and two direct potential approaches, the Venant and the partial integration method for the treatment of the dipole singularity is examined. The examination is performed in a four-compartment sphere model with anisotropic skull layer, where quasi-analytical solutions allow for an exact quantification of computational speed versus numerical error. Specifically-tuned constrained Delaunay tetrahedralization (CDT) FE meshes lead to high accuracies for both the full subtraction and the direct potential approaches. Best accuracies are achieved by the full subtraction approach if the homogeneity condition is fulfilled. It is shown that the AMG-CG achieves an order of magnitude higher computational speed than the CG with the standard preconditioners with an increasing gain factor when decreasing mesh size. Our results should broaden the application of accurate and fast high-resolution FE volume conductor modeling in source analysis routine.

A full subtraction approach for finite element method based source analysis using constrained Delaunay tetrahedralisation

A mathematical dipole is widely used as a model for the primary current source in electroencephalography (EEG) source analysis. In the governing Poisson-type differential equation, the dipole leads to a singularity on the right-hand side, which has to be treated specifically. In this paper, we will present a full subtraction approach where the total potential is divided into a singularity and a correction potential. The singularity potential is due to a dipole in an infinite region of homogeneous conductivity. The correction potential is computed using the finite element (FE) method. Special care is taken in order to evaluate the right-hand side integral appropriately with the objective of achieving highest possible convergence order for linear basis functions. Our new approach allows the construction of transfer matrices for fast computation of the inverse problem for anisotropic volume conductors. A constrained Delaunay tetrahedralisation (CDT) approach is used for the generation of high-quality FE meshes. We validate the new approach in a four-layer sphere model with a highly conductive cerebrospinal fluid (CSF) and an anisotropic skull compartment. For radial and tangential sources with eccentricities up to 1 mm below the CSF compartment, we achieve a maximal relative error of 0.71% in a CDT-FE model with 360 k nodes which is not locally refined around the source singularity and therefore useful for arbitrary dipole locations. The combination of the full subtraction approach with the high quality CDT meshes leads to accuracies that, to the best of the author's knowledge, have not yet been presented before.

EEG source analysis of epileptiform activity using a 1 mm anisotropic hexahedra finite element head model

The major goal of the evaluation in presurgical epilepsy diagnosis for medically intractable patients is the precise reconstruction of the epileptogenic foci, preferably with non-invasive methods. This paper evaluates whether surface electroencephalography (EEG) source analysis based on a 1 mm anisotropic finite element (FE) head model can provide additional guidance for presurgical epilepsy diagnosis and whether it is practically feasible in daily routine. A 1 mm hexahedra FE volume conductor model of the patient's head with special focus on accurately modeling the compartments skull, cerebrospinal fluid (CSF) and the anisotropic conducting brain tissues was constructed using non-linearly co-registered T1-, T2- and diffusion-tensor-magnetic resonance imaging data. The electrodes of intra-cranial EEG (iEEG) measurements were extracted from a co-registered computed tomography image. Goal function scan (GFS), minimum norm least squares (MNLS), standardized low resolution electromagnetic tomography (sLORETA) and spatio-temporal current dipole modeling inverse methods were then applied to the peak of the averaged ictal discharges EEG data. MNLS and sLORETA pointed to a single center of activity. Moving and rotating single dipole fits resulted in an explained variance of more than 97%. The non-invasive EEG source analysis methods localized at the border of the lesion and at the border of the iEEG electrodes which mainly received ictal discharges. Source orientation was towards the epileptogenic tissue. For the reconstructed superficial source, brain conductivity anisotropy and the lesion conductivity had only a minor influence, whereas a correct modeling of the highly conducting CSF compartment and the anisotropic skull was found to be important. The proposed FE forward modeling approach strongly simplifies meshing and reduces run-time (37 ms for one forward computation in the model with 3.1 million unknowns), corroborating the practical feasibility of the approach.

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.

Detection of focal epileptiform events in the EEG by spatio-temporal dipole clustering

OBJECTIVE

Methods for the detection of epileptiform events can be broadly divided into two main categories: temporal detection methods that exploit the EEG's temporal characteristics, and spatial detection methods that base detection on the results of an implicit or explicit source analysis. We describe how the framework of a spatial detection method was extended to improve its performance by including temporal information. This results in a method that provides (i) automated localization of an epileptogenic focus and (ii) detection of focal epileptiform events in an EEG recording. For the detection, only one threshold value needs to be set.

METHODS

The method comprises five consecutive steps: (1) dipole source analysis in a moving window, (2) automatic selection of focal brain activity, (3) dipole clustering to arrive at the identification of the epileptiform cluster, (4) derivation of a spatio-temporal template of the epileptiform activity, and (5) template matching. Routine EEG recordings from eight paediatric patients with focal epilepsy were labelled independently by two experts. The method was evaluated in terms of (i) ability to identify the epileptic focus, (ii) validity of the derived template, and (iii) detection performance. The clustering performance was evaluated using a leave-one-out cross validation. Detection performance was evaluated using Precision-Recall curves and compared to the performance of two temporal (mimetic and wavelet based) and one spatial (dipole analysis based) detection methods.

RESULTS

The method succeeded in identifying the epileptogenic focus in seven of the eight recordings. For these recordings, the mean distance between the epileptic focus estimated by the method and the region indicated by the labelling of the experts was 8mm. Except for two EEG recordings where the dipole clustering step failed, the derived template corresponded to the epileptiform activity marked by the experts. Over the eight EEGs, the method showed a mean sensitivity and selectivity of 92 and 77%, respectively.

CONCLUSIONS

The method allows automated localization of the epileptogenic focus and shows good agreement with the region indicated by the labelling of the experts. If the dipole clustering step is successful, the method allows a detection of the focal epileptiform events, and gave a detection performance comparable or better to that of the other methods.

SIGNIFICANCE

The identification and quantification of epileptiform events is of considerable importance in the diagnosis of epilepsy. Our method allows the automatic identification of the epileptic focus, which is of value in epilepsy surgery. The method can also be used as an offline exploration tool for focal EEG activity, displaying the dipole clusters and corresponding time series.