Influence of metallic artifact filtering on MEG signals for source localization during interictal epileptiform activity

Objective. Medical intractable epilepsy is a common condition that affects 40% of epileptic patients that generally have to undergo resective surgery. Magnetoencephalography (MEG) has been increasingly used to identify the epileptogenic foci through equivalent current dipole (ECD) modeling, one of the most accepted methods to obtain an accurate localization of interictal epileptiform discharges (IEDs). Modeling requires that MEG signals are adequately preprocessed to reduce interferences, a task that has been greatly improved by the use of blind source separation (BSS) methods. MEG recordings are highly sensitive to metallic interferences originated inside the head by implanted intracranial electrodes, dental prosthesis, etc and also coming from external sources such as pacemakers or vagal stimulators. To reduce these artifacts, a BSS-based fully automatic procedure was recently developed and validated, showing an effective reduction of metallic artifacts in simulated and real signals (Migliorelli et al 2015 J. Neural Eng. 12 046001). The main objective of this study was to evaluate its effects in the detection of IEDs and ECD modeling of patients with focal epilepsy and metallic interference. Approach. A comparison between the resulting positions of ECDs was performed: without removing metallic interference; rejecting only channels with large metallic artifacts; and after BSS-based reduction. Measures of dispersion and distance of ECDs were defined to analyze the results. Main results. The relationship between the artifact-to-signal ratio and ECD fitting showed that higher values of metallic interference produced highly scattered dipoles. Results revealed a significant reduction on dispersion using the BSS-based reduction procedure, yielding feasible locations of ECDs in contrast to the other two approaches. Significance. The automatic BSS-based method can be applied to MEG datasets affected by metallic artifacts as a processing step to improve the localization of epileptic foci.Postprint (published version

An Application of Gaussian Processes on Ocular Artifact Removal from EEG

International audienceConsequences of eye movements are one of the main inferences that distort the brain EEG recordings. In this paper, a multi-modal approach is used to estimate the ocular artifacts in the EEG: both vertical and horizontal eye movement signals recoded by an eye tracker are used as a reference to denoise the EEG. A Gaussian process, i.e. a second order statistics method, is assumed to model the link between the eye tracker signals and the EEG signals. The proposed method is thus a non-linear extension of the well-known adaptive filtering and can be applied with a single EEG signal contrary to independent component analysis (ICA) which is extensively used. The results show the applicability and the efficiency of this model on the ocular artifact removal

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

Validating and improving the correction of ocular artifacts in electro-encephalography

For modern applications of electro-encephalography, including brain computer interfaces and single-trial Event Related Potential detection, it is becoming increasingly important that artifacts are accurately removed from a recorded electro-encephalogram (EEG) without affecting the part of the EEG that reflects cerebral activity. Ocular artifacts are caused by movement of the eyes and the eyelids. They occur frequently in the raw EEG and are often the most prominent artifacts in EEG recordings. Their accurate removal is therefore an important procedure in nearly all electro-encephalographic research. As a result of this, a considerable number of ocular artifact correction methods have been introduced over the past decades. A selection of these methods, which contains some of the most frequently used correction methods, is given in Section 1.5. When two different correction methods are applied to the same raw EEG, this usually results in two different corrected EEGs. A measure for the accuracy of correction should indicate how well each of these corrected EEGs recovers the part of the raw EEG that truly reflects cerebral activity. The fact that this accuracy cannot be determined directly from a raw EEG is intrinsic to the need for artifact removal. If, based on a raw EEG, it would be possible to derive an exact reference on what the corrected EEG should be, then there would not be any need for adequate artifact correction methods. Estimating the accuracy of correction methods is mostly done either by using models to simulate EEGs and artifacts, or by manipulating the experimental data in such a way that the effects of artifacts to the raw EEG can be isolated. In this thesis, modeling of EEG and artifact is used to validate correction methods based on simulated data. A new correction method is introduced which, unlike all existing methods, uses a camera to monitor eye(lid) movements as a basis for ocular artifact correction. The simulated data is used to estimate the accuracy of this new correction method and to compare it against the estimated accuracy of existing correction methods. The results of this comparison suggest that the new method significantly increases correction accuracy compared to the other methods. Next, an experiment is performed, based on which the accuracy of correction can be estimated on raw EEGs. Results on this experimental data comply very well with the results on the simulated data. It is therefore concluded that using a camera during EEG recordings provides valuable extra information that can be used in the process of ocular artifact correction. In Chapter 2, a model is introduced that assists in estimating the accuracy of eye movement artifacts for simulated EEG recordings. This model simulates EEG and eye movement artifacts simultaneously. For this, the model uses a realistic representation of the head, multiple dipoles to model cerebral and ocular electrical activity, and the boundary element method to calculate changes in electrical potential at different positions on the scalp. With the model, it is possible to simulate different data sets as if they are recorded using different electrode configurations. Signal to noise ratios are used to assess the accuracy of these six correction methods for various electrode configurations before and after applying six different correction methods. Results show that out of the six methods, second order blind identification, SOBI, and multiple linear regression, MLR, correct most accurately overall as they achieve the highest rise in signal to noise ratio. The occurrence of ocular artifacts is linked to changes in eyeball orientation. In Chapter 2 an eye tracker is used to record pupil position, which is closely linked to eyeball orientation. The pupil position information is used in the model to simulate eye movements. Recognizing the potential benefit of using an eye tracker not only for simulations, but also for correction, Chapter 3 introduces an eye movement artifact correction method that exploits the pupil position information that is provided by an eye tracker. Other correction methods use the electrooculogram (EOG) and/or the EEG to estimate ocular artifacts. Because both the EEG and the EOG recordings are susceptive to cerebral activity as well as to ocular activity, these other methods are at risk of overcorrecting the raw EEG. Pupil position information provides a reference that is linked to the ocular artifact in the EEG but that cannot be affected by cerebral activity, and as a result the new correction method avoids having to solve traditionally problematic issues like forward/backward propagation and evaluating the accuracy of component extraction. By using both simulated and experimental data, it is determined how pupil position influences the raw EEG and it is found that this relation is linear or quadratic. A Kalman filter is used for tuning of the parameters that specify the relation. On simulated data, the new method performs very well, resulting in an SNR after correction of over 10 dB for various patterns of eye movements. When compared to the three methods that performed best in the evaluation of Chapter 2, only the SOBI method which performed best in that evaluation shows similar results for some of the eye movement patterns. However, a serious limitation of the correction method is its inability to correct blink artifacts. In order to increase the variety of applications for which the new method can be used, the new correction should be improved in a way that enables it to correct the raw EEG for blinking artifacts. Chapter 4 deals with implementing such improvements based on the idea that a more advanced eye-tracker should be able to detect both the pupil position and the eyelid position. The improved eye tracker-based ocular artifact correction method is named EYE. Driven by some practical limitations regarding the eye tracking device currently available to us, an alternative way to estimate eyelid position is suggested, based on an EOG recorded above one eye. The EYE method can be used with both the eye tracker information or with the EOG substitute. On simulated data, accuracy of the EYE method is estimated using the EOGbased eyelid reference. This accuracy is again compared against the six other correction methods. Two different SNR-based measures of accuracy are proposed. One of these quantifies the correction of the entire simulated data set and the other focuses on those segments containing simulated blinking artifacts. After applying EYE, an average SNR of at least 9 dB for both these measures is achieved. This implies that the power of the corrected signal is at least eight times the power of the remaining noise. The simulated data sets contain a wide range of eye movements and blink frequencies. For almost all of these data sets, 16 out of 20, the correction results for EYE are better than for any of the other evaluated correction method. On experimental data, the EYE method appears to adequately correct for ocular artifacts as well. As the detection of eyelid position from the EOG is in principle inferior to the detection of eyelid position with the use of an eye tracker, these results should also be considered as an indicator of even higher accuracies that could be obtained with a more advanced eye tracker. Considering the simplicity of the MLR method, this method also performs remarkably well, which may explain why EOG-based regression is still often used for correction. In Chapter 5, the simulation model of Chapter 2 is put aside and, alternatively, experimentally recorded data is manipulated in a way that correction inaccuracies can be highlighted. Correction accuracies of eight correction methods, including EYE, are estimated based on data that are recorded during stop-signal tasks. In the analysis of these tasks it is essential that ocular artifacts are adequately removed because the task-related ERPs, are located mostly at frontal electrode positions and are low-amplitude. These data are corrected and subsequently evaluated. For the eight methods, the overall ranking of estimated accuracy in Figure 5.3, corresponds very well with the correction accuracy of these methods on simulated data as was found in Chapter 4. In a single-trial correction comparison, results suggest that the EYE corrected EEG, is not susceptible to overcorrection, whereas the other corrected EEGs are

Acute sleep deprivation induces a local brain transfer information increase in the frontal cortex in a widespread decrease context

Sleep deprivation (SD) has adverse effects on mental and physical health, affecting the cognitive abilities and emotional states. Specifically, cognitive functions and alertness are known to decrease after SD. The aim of this work was to identify the directional information transfer after SD on scalp EEG signals using transfer entropy (TE). Using a robust methodology based on EEG recordings of 18 volunteers deprived from sleep for 36 h, TE and spectral analysis were performed to characterize EEG data acquired every 2 h. Correlation between connectivity measures and subjective somnolence was assessed. In general, TE showed medium-and long-range significant decreases originated at the occipital areas and directed towards different regions, which could be interpreted as the transfer of predictive information from parieto-occipital activity to the rest of the head. Simultaneously, short-range increases were obtained for the frontal areas, following a consistent and robust time course with significant maps after 20 h of sleep deprivation. Changes during sleep deprivation in brain network were measured effectively by TE, which showed increased local connectivity and diminished global integration. TE is an objective measure that could be used as a potential measure of sleep pressure and somnolence with the additional property of directed relationships.Postprint (published version

Acute sleep deprivation induces a local brain transfer information increase in the frontal cortex in a widespread decrease context

Sleep deprivation (SD) has adverse effects on mental and physical health, affecting the cognitive abilities and emotional states. Specifically, cognitive functions and alertness are known to decrease after SD. The aim of this work was to identify the directional information transfer after SD on scalp EEG signals using transfer entropy (TE). Using a robust methodology based on EEG recordings of 18 volunteers deprived from sleep for 36 h, TE and spectral analysis were performed to characterize EEG data acquired every 2 h. Correlation between connectivity measures and subjective somnolence was assessed. In general, TE showed medium-and long-range significant decreases originated at the occipital areas and directed towards different regions, which could be interpreted as the transfer of predictive information from parieto-occipital activity to the rest of the head. Simultaneously, short-range increases were obtained for the frontal areas, following a consistent and robust time course with significant maps after 20 h of sleep deprivation. Changes during sleep deprivation in brain network were measured effectively by TE, which showed increased local connectivity and diminished global integration. TE is an objective measure that could be used as a potential measure of sleep pressure and somnolence with the additional property of directed relationships.Postprint (published version

Identification of audio evoked response potentials in ambulatory EEG data

Electroencephalography (EEG) is commonly used for observing brain function over a period of time. It employs a set of invasive electrodes on the scalp to measure the electrical activity of the brain. EEG is mainly used by researchers and clinicians to study the brain’s responses to a specific stimulus - the event-related potentials (ERPs). Different types of undesirable signals, which are known as artefacts, contaminate the EEG signal. EEG and ERP signals are very small (in the order of microvolts); they are often obscured by artefacts with much larger amplitudes in the order of millivolts. This greatly increases the difficulty of interpreting EEG and ERP signals.Typically, ERPs are observed by averaging EEG measurements made with many repetitions of the stimulus. The average may require many tens of repetitions before the ERP signal can be observed with any confidence. This greatly limits the study and useof ERPs. This project explores more sophisticated methods of ERP estimation from measured EEGs. An Optimal Weighted Mean (OWM) method is developed that forms a weighted average to maximise the signal to noise ratio in the mean. This is developedfurther into a Bayesian Optimal Combining (BOC) method where the information in repetitions of ERP measures is combined to provide a sequence of ERP estimations with monotonically decreasing uncertainty. A Principal Component Analysis (PCA) isperformed to identify the basis of signals that explains the greatest amount of ERP variation. Projecting measured EEG signals onto this basis greatly reduces the noise in measured ERPs. The PCA filtering can be followed by OWM or BOC. Finally, crosschannel information can be used. The ERP signal is measured on many electrodes simultaneously and an improved estimate can be formed by combining electrode measurements. A MAP estimate, phrased in terms of Kalman Filtering, is developed using all electrode measurements.The methods developed in this project have been evaluated using both synthetic and measured EEG data. A synthetic, multi-channel ERP simulator has been developed specifically for this project.Numerical experiments on synthetic ERP data showed that Bayesian Optimal Combining of trial data filtered using a combination of PCA projection and Kalman Filtering, yielded the best estimates of the underlying ERP signal. This method has been applied to subsets of real Ambulatory Electroencephalography (AEEG) data, recorded while participants performed a range of activities in different environments. From this analysis, the number of trials that need to be collected to observe the P300 amplitude and delay has been calculated for a range of scenarios

Eye blink characterization from frontal EEG electrodes using source separation and pattern recognition algorithms

Due to its major safety applications, including safe driving, mental fatigue estimation is a rapidly growing research topic in the engineering field. Most current mental fatigue monitoring systems analyze brain activity through electroencephalography (EEG). Yet eye blink analysis can also be added to help characterize fatigue states. It usually requires the use of additional devices, such as EOG electrodes, uncomfortable to wear, or more expensive eye trackers. However, in this article, a method is proposed to evaluate eye blink parameters using frontal EEG electrodes only. EEG signals, which are generally corrupted by ocular artifacts, are decomposed into sources by means of a source separation algorithm. Sources are then automatically classified into ocular or non-ocular sources using temporal, spatial and frequency features. The selected ocular source is back propagated in the signal space and used to localize blinks by means of an adaptive threshold, and then to characterize detected blinks. The method, validated on 11 different subjects, does not require any prior tuning when applied to a new subject, which makes it subject-independent. The vertical EOG signal was recorded during an experiment lasting 90 min in which the participants’ mental fatigue increased. The blinks extracted from this signal were compared to those extracted using frontal EEG electrodes. Very good performances were obtained with a true detection rate of 89% and a false alarm rate of 3%. The correlation between the blink parameters extracted from both recording modalities was 0.81 in average

Aberrant activity in conceptual networks underlies N400 deficits and unusual thoughts in schizophrenia.

BackgroundThe N400 event-related potential (ERP) is triggered by meaningful stimuli that are incongruous, or unmatched, with their semantic context. Functional magnetic resonance imaging (fMRI) studies have identified brain regions activated by semantic incongruity, but their precise links to the N400 ERP are unclear. In schizophrenia (SZ), N400 amplitude reduction is thought to reflect overly broad associations in semantic networks, but the abnormalities in brain networks underlying deficient N400 remain unknown. We utilized joint independent component analysis (JICA) to link temporal patterns in ERPs to neuroanatomical patterns from fMRI and investigate relationships between N400 amplitude and neuroanatomical activation in SZ patients and healthy controls (HC).MethodsSZ patients (n = 24) and HC participants (n = 25) performed a picture-word matching task, in which words were either matched (APPLE→apple) by preceding pictures, or were unmatched by semantically related (in-category; IC, APPLE→lemon) or unrelated (out of category; OC, APPLE→cow) pictures, in separate ERP and fMRI sessions. A JICA "data fusion" analysis was conducted to identify the fMRI brain regions specifically associated with the ERP N400 component. SZ and HC loading weights were compared and correlations with clinical symptoms were assessed.ResultsJICA identified an ERP-fMRI "fused" component that captured the N400, with loading weights that were reduced in SZ. The JICA map for the IC condition showed peaks of activation in the cingulate, precuneus, bilateral temporal poles and cerebellum, whereas the JICA map from the OC condition was linked primarily to visual cortical activation and the left temporal pole. Among SZ patients, fMRI activity from the IC condition was inversely correlated with unusual thought content.ConclusionsThe neural networks associated with the N400 ERP response to semantic violations depends on conceptual relatedness. These findings are consistent with a distributed network underlying neural responses to semantic incongruity including unimodal visual areas as well as integrative, transmodal areas. Unusual thoughts in SZ may reflect impaired processing in transmodal hub regions such as the precuneus, leading to overly broad semantic associations

Limitations of ICA for artefact removal

This paper reports analysis of the limitations of using independent component analysis (ICA) for biosignal analysis especially artefact removal. The possible difficulty is that there are limited number of electrodes (recordings) making it an overcomplete problem (non-square ICA). The other difficulty is the distribution of biosignal being close to Gaussian. These two properties of the signals may make these outside the standard ICA application. This paper reports that ICA is able to successfully separate the biosignals if the number of recordings are not less than the number of sources. If that is not the case, ICA separates artefact component only when the corresponding artefact is predominant. The experiments demonstrate that the results are not reliable and hence the authors recommend that caution should be exercised before using ICA for such application