Predictability of epileptic seizures by fusion of scalp EEG and fMRI

The systems for prediction of epileptic seizure investigated in recent years mainly rely on the traditional nonlinear analysis of the brain signals from intracranial electroencephalograph (EEG) recordings. The overall objective of this work focuses on investigation of the predictability of seizure from the scalp signals by applying effective blind source separation (BSS) techniques to scalp EEGs, in which the epileptic seizures are considered as independent components of the scalp EEGs. The ultimate goal of the work is to pave the way for epileptic seizure prediction from the scalp EEG. The main contributions of this research are summarized as follows. Firstly, a novel constrained topographic independent component analysis (CTICA) algorithm is developed for the improved separation of the epileptic seizure signals. The related CTICA model is more suitable for brain signal separation due to the relaxation of the independence assumption, as the source signals geometrically close to each other are assumed to have some dependencies. By incorporating the spatial and frequency information of seizure signals as the constraint, CTICA achieves a better performance in separating the seizure signals in comparison with other conventional ICA methods. Secondly, the predictability of seizure is investigated. The traditional method for quantification of the nonlinear dynamics of time series is employed to quantify the level of chaos of the estimated sources. The simultaneously recorded intracranial and scalp EEGs are used for the comparison of the results. The experiment results demonstrate that the separated seizure sources have a similar transition trend as those achieved from the intracranial EEGs. Thirdly, simultaneously recorded EEG and functional Magnetic Resonance Imaging (fMRI) is studied in order to validate the activated area of the brain related to the seizure sources. An effective method to remove the fMRI scanner artifacts from the scalp EEG is established by applying the blind source extraction (BSE) algorithm. The results show that the effect of fMRI scanner artifacts has been reduced in scalp EEG recordings. Finally, a data driven model, spatial ICA (SICA) subject to EEG as the temporal constraint is proposed in order to detect the Blood Oxygen-Level Dependence (BOLD) from the seizure fMRI. In contrast to the popular model driven method General Linear Model (GLM), SICA does not rely on any predefined hemodynamic response function. It is based on the fact that brain areas executing different tasks are spatially independent. Therefore SICA works perfectly for non-event-related fMRI analysis such as seizure fMRI. By incorporating the temporal information existing within the EEG as the constraint, the superiority of the proposed constrained SICA is validated in terms of better algorithm convergence and a higher correlation between the time courses of the component and the seizure EEG signals as compared to SICA

Methods for noninvasive localization of focal epileptic activity with magnetoencephalography

Magnetoencephalography (MEG) is a noninvasive brain signal acquisition technique that provides excellent temporal resolution and a whole-head coverage allowing the spatial mapping of sources. These characteristics make MEG an appropriate technique to localize the epileptogenic zone (EZ) in the preoperative evaluation of refractory epilepsy. Presurgical evaluation with MEG can guide the placement of intracranial EEG (iEEG), the current gold standard in the clinical practice, and even supply sufficient information for a surgical intervention without invasive recordings, reducing invasiveness, discomfort, and cost of the presurgical epilepsy diagnosis. However, MEG signals have low signal-to-noise ratio compared with iEEG and can sometimes be affected by noise that masks or distorts the brain activity. This may prevent the detection of interictal epileptiform discharges (IEDs) and high-frequency oscillations (HFOs), two important biomarkers used in the preoperative evaluation of epilepsy. In this thesis, the reduction of two kinds of interference is aimed to improve the signal-to-noise ratio of MEG signals: metallic artifacts mask the activity of IEDs; and the high-frequency noise, that masks HFO activity. Considering the large number of MEG channels and the long duration of the recordings, reducing noise and marking events manually is a time-consuming task. The algorithms presented in this thesis provide automatic solutions aimed at the reduction of interferences and the detection of HFOs. Firstly, a novel automatic BSS-based algorithm to reduce metallic interference is presented and validated using simulated and real MEG signals. Three methods are tested: AMUSE, a second-order BSS technique; and INFOMAX and FastICA, based on high-order statistics. The automatic detection algorithm exploits the known characteristics of metallic-related interferences. Results indicate that AMUSE performes better when recovering brain activity and allows an effective removal of artifactual components.Secondly, the influence of metallic artifact filtering using the developed algorithm is evaluated in the source localization of IEDs in patients with refractory focal epilepsy. A comparison between the resulting positions of equivalent current dipoles (ECDs) produced by IEDs is performed: without removing metallic interference, rejecting only channels with large metallic artifacts, and after BSS-based reduction. The results show that a significant reduction on dispersion is achieved using the BSS-based reduction procedure, yielding feasible locations of ECDs in contrast to the other approaches. Finally, an algorithm for the automatic detection of epileptic ripples in MEG using beamformer-based virtual sensors is developed. The automatic detection of ripples is performed using a two-stage approach. In the first step, beamforming is applied to the whole head to determine a region of interest. In the second step, the automatic detection of ripples is performed using the time-frequency characteristics of these oscillations. The performance of the algorithm is evaluated using simultaneous intracranial EEG recordings as gold standard.The novel approaches developed in this thesis allow an improved noninvasive detection and localization of interictal epileptic biomarkers, which can help in the delimitation of the epileptogenic zone and guide the placement of intracranial electrodes, or even to determine these areas without additional invasive recordings. As a consequence of this improved detection, and given that interictal biomarkers are much more frequent and easy to record than ictal episodes, the presurgical evaluation process can be more comfortable for the patient and in a more economic way.La magnetoencefalografía (MEG) es una técnica no invasiva de adquisición de señales cerebrales que proporciona una excelente resolución temporal y una cobertura total de la cabeza, permitiendo el mapeo espacial de las fuentes cerebrales. Estas características hacen del MEG una técnica apropiada para localizar la zona epileptogénica (EZ) en la evaluación preoperatoria de la epilepsia refractaria. La evaluación prequirúrgica con MEG puede orientar la colocación del EEG intracraneal (iEEG), el actual modelo de referencia en la práctica clínica, e incluso suministrar información suficiente para una intervención quirúrgica sin registros invasivos; reduciendo la invasividad, la incomodidad y el costo del diagnóstico de la epilepsia prequirúrgica. Sin embargo, las señales MEG tienen baja relación señal ruido en comparación con el iEEG pudiendo imposibilitar la detección de descargas epileptiformes interictales (IEDs) y oscilaciones de alta frecuencia (HFOs), dos importantes biomarcadores utilizados en la evaluación preoperatoria de la epilepsia.En esta tesis, la reducción de dos tipos de interferencia está dirigida a mejorar la relación señal-ruido de la señal MEG: los artefactos metálicos que enmascaran la actividad de las IEDs; y el ruido de alta frecuencia, que enmascara la actividad de las HFOs. Debido al gran número de canales MEG y la larga duración de los registros, tanto reducir el ruido como seleccionar los biomarcadores manualmente es una tarea que consume mucho tiempo. Los algoritmos presentados en esta tesis aportan soluciones automáticas dirigidas a la reducción de interferencias y la detección de HFOs. En primer lugar, se presenta y valida un nuevo algoritmo automático basado en BSS para reducir interferencias metálicas mediante señales simuladas y reales. Se prueban tres métodos: AMUSE, una técnica BSS de segundo orden; y INFOMAX y FastICA, basados en estadísticos de orden superior. El algoritmo de detección automático utiliza las características conocidas de la señal producida por la interferencia metálica. Los resultados indican que AMUSE recupera mejor la actividad cerebral y permite una eliminación efectiva de componentes artefactuales.Posteriormente, se evalúa la influencia del filtrado de artefactos metálicos en la localización de IEDs en pacientes con epilepsia focal refractaria. Se realiza una comparación entre las posiciones resultantes de dipolos de corriente equivalentes (ECDs) producidos por IEDs: sin eliminar interferencias metálicas, rechazando solamente canales con elevados artefactos metálicos y, por último, después de una reducción utilizando el algoritmo BSS desarrollado. Los resultados muestran que se logra una reducción significativa en la dispersión utilizando el procedimiento de reducción basado en BSS, lo que produce ubicaciones factibles de los dipolos en contraste con los otros enfoques.En segundo lugar, se desarrolla un algoritmo para la detección automática ripples epilépticos en MEG utilizando sensores virtuales basados en la técnica de beamformer. La detección de ripples se realiza mediante un enfoque en dos etapas. Primero, se determina el área de interés usando beamformer. Posteriormente, se realiza la detección automática de ripples utilizando las características en tiempo-frecuencia. El rendimiento del algoritmo se evalúa utilizando registros iEEG simultáneos.Los nuevos enfoques desarrollados en esta tesis permiten una detección no invasiva mejor de los biomarcadores interictales, que pueden ayudar a delimitar la zona epileptogénica y guiar la colocación de electrodos intracraneales, o incluso determinar estas áreas sin este tipo de registros. Como consecuencia de esta mejora en la detección, y dado que los biomarcadores interictales son mucho más frecuentes y fáciles de registrar que los episodios ictales, la evaluación prequirúrgica puede ser más cómoda y menos costosa para el paciente.Postprint (published version

Detection and removal of eyeblink artifacts from EEG using wavelet analysis and independent component analysis

Electrical signals generated by brain activity that are measured by the electroencephalogram can be distorted by electrical activity originating from eyeblinks and eye movements. This thesis proposes a new technique to identify and remove eyeblink artifacts from EEG data. An algorithm using a combination of wavelet analysis and independent component analysis (ICA) is implemented to detect the temporal location of the eyeblink artifact and eliminate it without compromising the integrity of the primary EEG data. The discrete wavelet transform is performed on 10 second epochs of data to detect the occurrence of ocular artifact. ICA is used to separate out the independent components within the data and the temporal locations of the eyeblink are used to remove the artifact and reconstruct the EEG data without that source of distortion. The results obtained indicate that the technique implemented may be robust enough to effectively process EEG data and is capable of removing eyeblink artifacts successfully when they are prominent and the data does not contain a great deal of movement artifact. The results show an 88.68% detection rate, a false positive rate of 4.03%, and an 87.23% removal rate for all eyeblinks that were accurately detected. The statistics obtained compared favorably with work done by others in this field of investigation

A Computational Framework to Support the Automated Analysis of Routine Electroencephalographic Data

Epilepsy is a condition in which a patient has multiple unprovoked seizures which are not precipitated by another medical condition. It is a common neurological disorder that afflicts 1% of the population of the US, and is sometimes hard to diagnose if seizures are infrequent. Routine Electroencephalography (rEEG), where the electrical potentials of the brain are recorded on the scalp of a patient, is one of the main tools for diagnosing because rEEG can reveal indicators of epilepsy when patients are in a non-seizure state. Interpretation of rEEG is difficult and studies have shown that 20-30% of patients at specialized epilepsy centers are misdiagnosed. An improved ability to interpret rEEG could decrease the misdiagnosis rate of epilepsy. The difficulty in diagnosing epilepsy from rEEG stems from the large quantity, low signal to noise ratio (SNR), and variability of the data. A usual point of error for a clinician interpreting rEEG data is the misinterpretation of PEEs (paroxysmal EEG events) ( short bursts of electrical activity of high amplitude relative to the surrounding signals that have a duration of approximately .1 to 2 seconds). Clinical interpretation of PEEs could be improved with the development of an automated system to detect and classify PEE activity in an rEEG dataset. Systems that have attempted to automatically classify PEEs in the past have had varying degrees of success. These efforts have been hampered to a large extent by the absence of a \gold standard\u27 data set that EEG researchers could use. In this work we present a distributed, web-based collaborative system for collecting and creating a gold standard dataset for the purpose of evaluating spike detection software. We hope to advance spike detection research by creating a performance standard that facilitates comparisons between approaches of disparate research groups. Further, this work endeavors to create a new, high performance parallel implementation of ICA (independent component analysis), a potential preprocessing step for PEE classification. We also demonstrate tools for visualization and analysis to support the initial phases of spike detection research. These tools will first help to develop a standardized rEEG dataset of expert EEG interpreter opinion with which automated analysis can be trained and tested. Secondly, it will attempt to create a new framework for interdisciplinary research that will help improve our understanding of PEEs in rEEG. These improvements could ultimately advance the nuanced art of rEEG interpretation and decrease the misdiagnosis rate that leads to patients suering inappropriate treatment

Enhancing brain-computer interfacing through advanced independent component analysis techniques

A Brain-computer interface (BCI) is a direct communication system between a brain and an external device in which messages or commands sent by an individual do not pass through the brain’s normal output pathways but is detected through brain signals. Some severe motor impairments, such as Amyothrophic Lateral Sclerosis, head trauma, spinal injuries and other diseases may cause the patients to lose their muscle control and become unable to communicate with the outside environment. Currently no effective cure or treatment has yet been found for these diseases. Therefore using a BCI system to rebuild the communication pathway becomes a possible alternative solution. Among different types of BCIs, an electroencephalogram (EEG) based BCI is becoming a popular system due to EEG’s fine temporal resolution, ease of use, portability and low set-up cost. However EEG’s susceptibility to noise is a major issue to develop a robust BCI. Signal processing techniques such as coherent averaging, filtering, FFT and AR modelling, etc. are used to reduce the noise and extract components of interest. However these methods process the data on the observed mixture domain which mixes components of interest and noise. Such a limitation means that extracted EEG signals possibly still contain the noise residue or coarsely that the removed noise also contains part of EEG signals embedded. Independent Component Analysis (ICA), a Blind Source Separation (BSS) technique, is able to extract relevant information within noisy signals and separate the fundamental sources into the independent components (ICs). The most common assumption of ICA method is that the source signals are unknown and statistically independent. Through this assumption, ICA is able to recover the source signals. Since the ICA concepts appeared in the fields of neural networks and signal processing in the 1980s, many ICA applications in telecommunications, biomedical data analysis, feature extraction, speech separation, time-series analysis and data mining have been reported in the literature. In this thesis several ICA techniques are proposed to optimize two major issues for BCI applications: reducing the recording time needed in order to speed up the signal processing and reducing the number of recording channels whilst improving the final classification performance or at least with it remaining the same as the current performance. These will make BCI a more practical prospect for everyday use. This thesis first defines BCI and the diverse BCI models based on different control patterns. After the general idea of ICA is introduced along with some modifications to ICA, several new ICA approaches are proposed. The practical work in this thesis starts with the preliminary analyses on the Southampton BCI pilot datasets starting with basic and then advanced signal processing techniques. The proposed ICA techniques are then presented using a multi-channel event related potential (ERP) based BCI. Next, the ICA algorithm is applied to a multi-channel spontaneous activity based BCI. The final ICA approach aims to examine the possibility of using ICA based on just one or a few channel recordings on an ERP based BCI. The novel ICA approaches for BCI systems presented in this thesis show that ICA is able to accurately and repeatedly extract the relevant information buried within noisy signals and the signal quality is enhanced so that even a simple classifier can achieve good classification accuracy. In the ERP based BCI application, after multichannel ICA the data just applied to eight averages/epochs can achieve 83.9% classification accuracy whilst the data by coherent averaging can reach only 32.3% accuracy. In the spontaneous activity based BCI, the use of the multi-channel ICA algorithm can effectively extract discriminatory information from two types of singletrial EEG data. The classification accuracy is improved by about 25%, on average, compared to the performance on the unpreprocessed data. The single channel ICA technique on the ERP based BCI produces much better results than results using the lowpass filter. Whereas the appropriate number of averages improves the signal to noise rate of P300 activities which helps to achieve a better classification. These advantages will lead to a reliable and practical BCI for use outside of the clinical laboratory