Next Article in Journal
Social Milieu Oriented Routing: A New Dimension to Enhance Network Security in WSNs
Previous Article in Journal
A Secure Scheme for Distributed Consensus Estimation against Data Falsification in Heterogeneous Wireless Sensor Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid EEG—Eye Tracker: Automatic Identification and Removal of Eye Movement and Blink Artifacts from Electroencephalographic Signal

Department of Cogno-Mechatronics Engineering, Pusan National University, 2 Busandaehak-ro 63beon-gil Geumjeong-gu, Busan 609-735, Korea
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(2), 241; https://doi.org/10.3390/s16020241
Submission received: 28 November 2015 / Revised: 5 February 2016 / Accepted: 14 February 2016 / Published: 19 February 2016
(This article belongs to the Section Physical Sensors)

Abstract

:
Contamination of eye movement and blink artifacts in Electroencephalogram (EEG) recording makes the analysis of EEG data more difficult and could result in mislead findings. Efficient removal of these artifacts from EEG data is an essential step in improving classification accuracy to develop the brain-computer interface (BCI). In this paper, we proposed an automatic framework based on independent component analysis (ICA) and system identification to identify and remove ocular artifacts from EEG data by using hybrid EEG and eye tracker system. The performance of the proposed algorithm is illustrated using experimental and standard EEG datasets. The proposed algorithm not only removes the ocular artifacts from artifactual zone but also preserves the neuronal activity related EEG signals in non-artifactual zone. The comparison with the two state-of-the-art techniques namely ADJUST based ICA and REGICA reveals the significant improved performance of the proposed algorithm for removing eye movement and blink artifacts from EEG data. Additionally, results demonstrate that the proposed algorithm can achieve lower relative error and higher mutual information values between corrected EEG and artifact-free EEG data.

1. Introduction

In recent years, non-invasive neuro-imaging has become a valuable research tool to understand the underlying functionality of the brain [1,2,3,4,5]. Electroencephalogram (EEG) with the advantages of portability and high temporal resolution is a non-invasive brain-imaging technique used to measure different physiological states of the brain with amplitude typically in order of a few microvolts [6,7]. Unfortunately, measurements from EEG are highly contaminated with eye movement and blink artifacts which are several times higher in magnitude as compared to neuronal activity [8,9,10,11,12,13]. This issue has become a recurrent problem, for example in brain-computer interface (BCI) where it has been proved to decrease the classification accuracy [14].
Various methods to tackle this challenging task have been proposed in the literature. A straight forward way to reduce ocular artifacts is to discard the artifactual epochs from EEG data. However, this may cause a considerable data loss related to neuronal activity and it also requires a huge amount of time. In contrast, several automated methods have been proposed to detect and remove/reduce ocular artifacts in past. These methods can be divided into two main categories, i.e., regression based methods [9,15,16,17,18,19,20] and blind source separation techniques [21,22,23,24,25,26,27].
In regression algorithms, propagation coefficients of ocular artifacts are calculated to estimate the amount of electrooculography (EOG) signal present in EEG signal and subsequently the estimated EOG is subtracted from EEG to obtain artifact-free EEG [28,29,30,31,32]. Since EOG also captures the neuronal activity from prefrontal cortex, therefore the common neuronal activity in EEG and EOG might be lost in regression methods and due to this bidirectional contamination problem these algorithms were proved to be less effective [29]. Figure 1A shows the schematic diagram of regression algorithms.
On the other hand, the second class of methods are based on the assumption that the neuronal activity and artifactual activity are independent from each other. Most commonly used method in blind source separation is independent component analysis (ICA) [21,31,32]. In ICA based algorithms, EEG signal is decomposed into several independent components (ICs), which are then classified into neuronal-activity-related and artifactual ICs. The selected artifactual ICs are then removed to obtain artifact-free EEG. The main issue of this methodology is the selection of artifactual components, usually this can be done by the visual inspection of an expert but this approach might lead to misclassification of ICs and divergent results [33]. Recently, several automated criteria have been developed to tackle this issue [26,32,33,34] and they showed significant improvement in terms of artifact removal from EEG data but the distortion produced due the removal of those ICs in EEG signal was left unaddressed, since an artifactual IC may contain neuronal activity along with artifacts [35]. The schematic diagram illustrating the ICA based removal of ocular artifacts is shown in Figure 1B.
Recently, researchers proposed to combine these two methodologies to remove ocular artifacts from EEG by utilizing their advantageous features [29,36]. Although these methods proved to be effective in removing ocular artifacts but they always require simultaneous EOG recordings, which is not plausible in applications like BCI. To this extent, Kierkels and coauthors [37] proposed to use an eye tracker for eye movement artifact removal from EEG. They used eye tracker to generate eye positions, which were used as inputs to Kalman filter to remove eye movement artifacts. Although, their method achieved improved results over other artifact removal techniques but it is not able to deal with blink artifacts. Later, Noureddin and colleagues [38] used a high speed eye tracker to propose a regression based technique with recursive least square and H filters. Plöchl and co-authors [33] proposed a simple criteria for identifying artifactual ICs based on the events of eye tracker and they replaced the selected ICs with zero to get clean EEG.
In this paper, hybrid EEG and eye tracker system is used to develop a novel adaptive framework which automatically detects and removes eye movement and blink artifacts from EEG data. In contrast to the EOG based algorithms, the proposed algorithm is developed such that eye tracker and frontal EEG electrodes are used to detect and remove ocular artifacts. The proposed methodology combines the advantageous features of ICA and system identification to remove eye movement and blink artifacts. The first step of the proposed algorithm is to obtain ICs by ICA decomposition of EEG data. In the next step, we proposed to use composite multi-scale entropy and eye tracker events to automatically identify blink and eye movement artifacts related ICs, respectively. These artifactual ICs are filtered with median absolute deviation to remove high magnitude ocular activities and then processed to an auto-regressive exogenous model to completely remove eye movement and blink artifacts. The parameters of the auto-regressive exogenous model are estimated using affine projection algorithm. In the final step, inverse ICA is used to reconstruct EEG signal by back projecting all ICs. The performance of the proposed algorithm is demonstrated through results on experimental and standard EEG datasets. Furthermore, the proposed algorithm is compared with two state-of-the-art techniques namely ADJUST based ICA [34] and REGICA [30]. Relative error and mutual information are used as an evaluation indexes to measure the ability of removing ocular artifacts by the proposed algorithm. Paired t-test is used to validate the significant improvement in removing ocular artifacts by the proposed algorithm over previous methods. Results show that the proposed method is efficient for automatic detection and removal of eye movement and blink artifacts from EEG signals. The schematic diagram and summary of the proposed algorithm are shown in Figure 2 and Table 1, respectively.
The remainder of this paper is organized as follows. In Section 2, a detailed description of the datasets and proposed algorithm is presented. Evaluation indexes used in this study are formulated in Section 3 followed by the results and discussion sections. Finally, the main idea and conclusions are summarized in last section.

2. Materials and Methods

2.1. Materials

In this paper, we used experimental and standard EEG datasets to demonstrate the performance of the proposed algorithm.

2.1.1. Participants

EEG and eye tracker measurements were acquired from 5 healthy subjects, all male, mean age 28. Experiment was conducted under the Declaration of Helsinki. The experiment was approved by the Institutional Review Board of Pusan National University. The experiment was conducted in a confined room with dim light to avoid environmental disturbance.

2.1.2. Experimental Procedure

All the participants were seated in an armchair at a distance of about 1m from a monitor screen (Samsung, SyncMaster B1940, 19”). The participants were asked to perform a task involving different eye movements and blinks. The experiment starts with the blank screen for 5 s. During this task, a red dot appears on the screen for the duration of 2 s at nine different positions in a square region of 960 × 960 pixels. All the participants were instructed to follow the moving dot. Before every position change of the dot, a fixation cross appeared on the screen for 1 s, thus providing the cue for the subject to follow the dot. Each participant was instructed to blink when “blink” appears on the screen. A blank screen appears for 5 s after every blink. All participants performed three experimental blocks. Each experimental block consisted of 18 saccade trials and 18 fixation trials.

2.1.3. EEG Recordings

EEG data were recorded using BrainAmp DC amplifier with an ActiCap 32-channel active electrode system developed by Brain Products GmbH, Gilching, Germany. All the data were sampled at a rate of 500 Hz. All the electrodes were placed according to international 10–20 system as shown in Figure 3A. The impedance of all the electrodes were reduced below 5 kΩ.

2.1.4. Eye Tracker Recordings

The eye movements were recorded with a video eye tracking system Eyelink 1000 developed by SR research Ltd., Ottawa, ON, Canada. The sampling rate was 1000 Hz. The velocity threshold of 30 °/s was used to define saccades, an acceleration threshold of 8000 °/s2 and a minimum deflection threshold of 0.1°. Figure 3B shows the distribution of saccades amplitude.

2.1.5. Preprocessing

In order to synchronize EEG and eye tracker, eye tracker data was down-sampled at the sampling rate of 500 Hz. EEG and eye tracker data were then aligned by cutting them into trials according to the triggers that were simultaneously sent to both, the EEG and the eye tracking system [32]. The EEG data was band pass filtered between 0.5–40 Hz. All the processing and analysis were done in Matlab (Mathworks) and EEGLAB using Intel core i3, 2.4 GHz with 4 GB RAM laptop.

2.1.6. Standard Dataset

The principle measure to evaluate the performance of the proposed algorithm is to check its ability of removing artifactual activities from standard dataset. Preparatory to such an evaluation, publicly available contaminated EEG dataset with eye tracker signals was utilized [39]. In this experiment, participant read lists of five words from left to right. Their task was to report whether the list contained the name of an animal. Eye movements were recorded binocularly with an Eyelink 1000 tracker at 1000 Hz.
EEG data was recorded from 72 channels with Biosemi Active amplifiers at a rate of 512 Hz. In the preprocessing part, EEG data was filtered between 0.5 and 40 Hz and the baseline was removed from all of the data.

2.2. Methods

2.2.1. Independent Component Analysis

ICA is a statistical technique used for decomposing multichannel data into several ICs under the following assumptions [40]:
  • The number of ICs are less than or equal to the number of observed signals.
  • The artifactual and cerebral sources are linearly mixed and statistically independent.
  • Propagation delays through the missing medium (brain) are negligible.
The basic purpose of ICA is to consider the non-Gaussianity of the measured signal and to find their projections. Mathematical mode of ICA, for the observed EEG data can be represented as
x ( k ) = A s ( k ) , k = 1 , 2 , 3 , ... , N
where x ( k ) M × 1 is the measured EEG signal, s ( k ) M × 1 is the corresponding IC, A M × M is the full rank unknown mixing matrix, k is discrete time, N is the number of samples and M is the number of ICs. Since the total number of ICs contributing to EEG data are unknown, therefore in this study they are supposed to be equal to the number of electrodes used in EEG data acquisition. Given x ( k ) , the issue is how to estimate both A and s ( k ) . The ICs s ^ i ( k ) , i = 1 , 2 , 3 , ... , M can be represented as
s ^ i ( k ) = w i T x ( k ) ,     k = 1 , 2 , 3 , ... , N
where w i is a column vector. After estimation of each w i , the ICs can be obtained by using the following expression
s ^ ( k ) = W x ( k ) ,     W A 1
The algorithm implemented in this paper is infomax ICA with default parameters using runica function of the EEGLAB tool box (MATLAB, Torrance, CA, USA) [41]. These parameters involved pre-sphering of the data and to avoid training if weight change was less then 10−6.

2.2.2. Features Computation

Horizontal Eye Movements

Horizontal eye movements generate large amplitude fluctuations in frontal channels that are typically slower than those of blinks, therefore not efficiently identifiable by the composite multi-scale entropy. To identify horizontal eye movement related artifacts, all the ICs were portioned into saccade and fixation epochs [33]. Saccade epochs were defined as the time between horizontal eye movements start and end as given by the eye tracker. An additional interval of 5 ms before and 10 ms after was added to saccade epochs. In contrast, the fixation epochs were defined as the time between saccade epochs. Finally, the ratio of the mean variance for saccade and fixation epochs was calculated for all the ICs:
R a t i o = m e a n ( var  i a n c e s a c c a d e ) m e a n ( var  i a n c e f i x a t i o n )
If for given IC the ratio of mean variance defined in Equation (7) was greater than 1.1, the corresponding IC was selected as artifact related IC and subsequently processed for correction.

Vertical Eye Movements

Since the time course of artifacts caused by vertical eye movements is similar to the one generated by horizontal eye movements, the feature described in Equation (7) can be used to identify vertical eye movements related ICs.

2.2.3. Median Absolute Deviation

Once the artifactual ICs are identified, they are then processed for correction by a two-step methodology. In the first step, the ocular activities (outliers) that are set to zero are only those notable ones that are of a considerably high magnitude. In this way, the amount of neuronal activity in the ocular artifacts related components can be retained. In the present study, median absolute deviation was used to detect and remove high-magnitude ocular activities from the components [45]. The step-wise procedure for such removal is as follows:
(1)
Evaluate the median absolute deviation of the identified ocular activity among the identified artifactual ICs (median absolute deviation is defined as the median of the absolute deviation from the median)
M A D = b M ( | u i ( k ) M ( u i ) | )
where M A D is the median absolute deviation, M is the median, M ( u i ) is the median of the ith artifactual IC, b is a constant;
(2)
If u i ( j ) exceeds the criteria calculated using Equation (9), it is thresholded to zero:
M ( u i ) 3 * M A D < u i ( k ) < M ( u i ) + 3 * M A D
u i ( k ) M ( u i ) M A D > | ± 3 |

2.2.4. Auto-Regressive Exogenous Model

The procedure described above will only remove those artifactual components which can be clearly seen and detectable in ocular artifacts related ICs. Auto-regressive exogenous model is used to completely remove ocular artifacts from identified ICs. It is expected that the amount of neuronal activity included in identified components is much lower than that of present in contaminated EEG. Therefore, in the present study we applied auto-regressive exogenous model to ocular artifacts related ICs instead of EEG data. A linear auto-regressive exogenous model can be used to remove ocular artifacts from artifactual ICs by the following equation
y ( k ) + i = 1 p a i y ( k i ) = j = 1 q b j E F p 1 ( k j ) + l = 1 r c l E F p 2 ( k l ) + e ( k )
where y is the output, E F p 1 and E F p 2 are the inputs of the auto-regressive exogenous model representing Fp1 and Fp2 electrode of EEG at discrete time k, a i , b j and c l are the parameters to be estimated, p, q and r represent the order of the model and e ( k ) is the error assumed to be white-noise disturbance. Mathematically, Equation (11) can be rewritten in linear regression as below
y ( k ) = X T ( k ) β + e ( k )
where X and β have the following form
X ( k ) = [ y ( k 1 ) , , y ( k p ) , E F p 1 ( k 1 ) , , E F p 1 ( k q ) , E F p 2 ( k r ) , , E F p 2 ( k r ) ] T β = [ a 1 , , a p , b 1 , , b q , c 1 , , c r ] T }

2.2.5. Affine Projection Algorithm

The error signal can be obtained by subtracting the estimated output from the desired signal. Mathematically, it can be written as [46]
e ( k ) = y ( k ) X T ( k ) β ^ ( k 1 )
Then, the objective of the estimation problem is to minimize the squared Euclidean norm of
β ^ ( k + 1 ) β ^ ( k ) 2 2
Subject to constraints
y ( k ) X T ( k ) β ^ ( k ) = 0
Thus, an update equation is required such that the difference between two consecutive estimations of the unknown parameters is minimized. This can be achieved by using the method of Lagrange multipliers which converts the constrained minimization into an unconstrained one. Thus, the cost function can be defined as
J ( k ) = β ^ ( k ) β ^ ( k 1 ) 2 2 + [ y ( k ) X T ( k ) β ^ ( k ) ] θ ( k )
where θ is the Lagrange multipliers vector. Taking the gradient of J(k) with respect to β ^ ( k ) and equating the result to zero, we can find
β ^ ( k ) = β ^ ( k 1 ) + 1 2 X ( k ) θ ( k )
Using Equations (12) and (14) and the Lagrange method, we can write
θ ( k ) = 2 [ X T ( k ) X ( k ) ] 1 e ( k )
Substituting Equation (19) into Equation (18), the optimal change in the parameter vector can be written as
β ^ ( k ) = β ^ ( k 1 ) + X ( k ) [ X T ( k ) X ( k ) ] 1 e ( k )
A step size parameter μ and regularization constant δ can be used to modify the above equation for efficient updating of the parameter vector as
β ^ ( k ) = β ^ ( k 1 ) + μ X ( k ) [ δ I + X T ( k ) X ( k ) ] 1 e ( k )

3. Evaluation Index

Since it is not possible to exactly evaluate the performance of proposed methodology that how much artifacts from EEG data has been removed due to the unknown contributions of the neuronal activity and ocular activity. That is, it is not possible to measure signal to artifact ratio. But the performance of the proposed algorithm can be evaluated using data in the intervals before artifact contamination [47]. To test how well the proposed algorithm performed in comparison to the conventional methods, we asked an independent EEG expert to tag all intervals in our data that he considered as ocular artifacts related. The selection was done by visual inspection of the EEG time series. All the remaining data was considered as artifact-free EEG data and used to calculate the performance metrics. In this study, we used two performance measures to quantitatively verify the ability of the proposed algorithm in removing ocular artifacts and how much the EEG signals are distorted after the artifact rejection procedure.

3.1. Relative Error

In the present study, the relative error criteria was used to evaluate the proposed algorithm’s utility in removing ocular artifacts from EEG signals and comparing it with the results of the conventional methodologies. Relative error is defined as [47]
R E = | E E G f r e e E E G o u t | | E E G f r e e |
where R E represents the relative error, E E G f r e e is the artifact-free EEG signal as selected by EEG expert, E E G o u t is the artifact corrected EEG signal from the proposed algorithm and | | denotes the norm calculation for vector.

3.2. Mutual Information

The mutual information index was calculated to measure the mutual dependence of the artifact-free EEG signal and output EEG from the proposed method. Mathematically, it is found by using Kullback-Leibler divergence between the probability distribution function as [24]
M I = f ( a , b ) log f ( a , b ) f ( a ) f ( b ) d a d b
where M I is the mutual information, f ( a , b ) is the joint pdf and f ( a ) and f ( b ) are the marginal pdfs. Mutual information is calculated using an open source MATLAB function minfo.m developed by Dr. Jason Palmer [48]. If the mutual information between artifact-free EEG and output EEG from proposed method is large, it means they are closely related.

4. Results

This paper presents an automatic framework based on ICA and auto-regressive exogenous model to identify and remove ocular artifacts from EEG signals by combining EEG and eye tracker. The effectiveness of the proposed algorithm was demonstrated using experimental and standard EEG datasets. The performance of the proposed algorithm is compared with two conventional methods, i.e., ICA and REGICA to verify the significant improvement of results. In this study, ADJUST as implemented in EEGLAB toolbox is used to represent ICA based algorithms and REGICA is used to represent methods based on the combination of ICA and regression.
Five experimental EEG datasets were used to verify the performance of the proposed algorithm. Figure 4 plots the results of artifact removal by the proposed algorithm for one subject. Figure 4A shows the experimental EEG data for one subject, Figure 4B, the corresponding ICs obtained from ICA decomposition of the EEG data, and Figure 4C, comparison of the artifact-free EEG data obtained after implementation of the proposed algorithm and conventional algorithms. In Figure 4B, ICs 2 and 13 are blink related components and ICs 9 and 23 are eye movement related components as identified by the proposed algorithm. It is evident in Figure 4C that the ocular artifacts were significantly removed by the proposed algorithm, in contrast to Figure 4A. Furthermore, the comparison with the conventional algorithms show the improved performance of the proposed algorithm.
Figure 5 compares the proposed method and ADJUST for artifact removal for one subject at Fp1 and Fp2. Figure 5A shows the contaminated EEG data; Figure 5B,C are the corresponding artifact-free EEG after implementation of the proposed algorithm and ADJUST, respectively; Figure 5D compares the proposed algorithm and ADJUST with the contaminated EEG data. The two black boxes on the left and two on the right in Figure 5D highlight the uncontaminated and contaminated regions of the EEG signals at Fp1 and Fp2, respectively, which are partially enlarged in Figure 5E. Figure 5E shows that ADJUST causes distortion and loss of neuronal activity from the EEG data, whereas the proposed algorithm successfully preserves the neuronal-activity-related EEG signal intact. Indeed, in Figure 5E, it can be seen that the proposed algorithm performs better in removing ocular artifacts and reconstructing the EEG signal. Figure 6 illustrates the comparison of corrected EEG by proposed method with REGICA method. Figure 6D, meanwhile, compares the proposed algorithm with the REGICA at Fp1 and Fp2, respectively, and indicates that the proposed algorithm offers significantly better performance (Figure 6E).
Figure 7 and Figure 8 show the comparison results of the proposed algorithm with ADJUST and REGICA for all subjects at Fp1 and Fp2, respectively. It can be seen that the proposed algorithm outperforms conventional methods in removing ocular artifacts from EEG data as well as in preserving the neuronal activity related EEG signal (enlarged panels). In order to investigate the effect of the different ocular artifacts reduction methods on artifact-free data in frequency domain, the power spectral density (PSD) is computed and compared.
For this purpose, 10 s of artifact-free EEG data before artifact contamination is selected by an expert. Then, the PSD was computed using pwlech function in MATLAB. The resulting PSD of one subject on a frontal electrode Fp1 and a most occipital electrode Oz is shown in Figure 9. It can be seen in Figure 9A,B (upper panel) that the proposed algorithm outperforms the conventional methods (enlarged panels), whereas the PSD at Oz (bottom panel) for all the algorithms show good agreement with the PSD of the artifact-free EEG.
Furthermore, Table 2 lists the relative error indices calculated for all three algorithms with respect to the five EEG datasets. The relative error values show the superior performance of the proposed algorithm over ADJUST and REGICA for all of the datasets. A paired t-test was run to determine if the relative error values differed statistically among the proposed algorithm, ADJUST and REGICA. It can be seen from Table 2 that the proposed algorithm is highly statistically significant when compared with ADJUST and REGICA, except for one subject against REGICA (p < 0.17). For the purpose of further validation, mutual information index is used to evaluate the performance of the proposed algorithm against ADJUST and REGICA. Average mutual information values of all subjects for all electrodes are listed in Table 3. The results of mutual information index show that the proposed algorithm preserved more mutual information between the artifact-free EEG signal and the reconstructed EEG signal as compared to the conventional methods.
Finally, the proposed algorithm was tested on a standard EEG dataset to determine its utility for removal of ocular artifacts from EEG. Since the standard dataset does not contain EOG signals, so we only used this dataset to compare the proposed algorithm with ADJUST. Figure 10 compares the proposed method with ADJUST for artifact removal from standard EEG data at Fp1 and Fp2. Figure 10A shows the standard EEG data with the ocular contamination, Figure 10B,C, the artifact-free EEG after implementation of the proposed algorithm and ADJUST, respectively, and Figure 10D, a comparison of the proposed algorithm and ADJUST with contaminated standard EEG data. The black boxes in Figure 10D highlight the uncontaminated and contaminated EEG-signal regions, which are partially enlarged in Figure 10E. It can be visualized that the performance of the proposed algorithm in removing ocular artifacts and maintaining the neuronal-activity-related EEG signal intact is significantly better than ADJUST (Figure 10E).

5. Discussion

The analysis of EEG signal always requires the identification and removal of artifacts due to eye movements and blinking. In this study, a novel algorithm, based on the combination of ICA and auto-regressive exogenous model, is proposed for automatic identification and removal of ocular activities from EEG. Since using EOG is not plausible for applications like BCI, therefore in the proposed algorithm eye tracker and frontal EEG electrodes are used to remove ocular artifacts from EEG signal. In the literature, regression based techniques are the most commonly used methods to remove ocular artifacts. Although these techniques are simple and fast but due to bidirectional contamination they proved to be less effective [29]. To overcome this issue, the proposed algorithm reconstruct ICs instead of applying regression to EEG signal. Our assumption lies with the fact that the amount of the cerebral activity included in the contaminated ICs is much lower compared to that existing in the contaminated EEG signals. So, as long as the cerebral and ocular activities are derived from independent sources, the cerebral activity included in the artifactual ICs tends to be minimal. It can, therefore, be assumed that the artifactual components contain less cerebral activity common to EOG/frontal EEG. In this way, filtering artifactual ICs with auto-regressive exogenous model will cause less removal of cerebral activity. Thus, the corrected EEG contains more cerebral information in contrast to the conventional regression analysis.
Recent efforts on artifact removal in EEG signal has shown a great utility of ICA. Although the success of ICA is encouraging, it should be treated with care [49]. Existing studies have focused almost extensively on the important reduction of the typical artifacts in ICA corrected EEG signal [50,51,52], while distortion of the cerebral part of EEG signal introduced by the method as a side effect have been left unattended [35]. First, the effectiveness of ICA strongly depends on the quality of the signal decomposition. Not all signal sources may be isolated into separate components and there are no definite means to evaluate whether or not contributions of other sources confound a particular component [33]. Additionally, selection of ocular artifacts related components is another issue in ICA based methodologies. Usually this can be done by inspecting the time series and topographic maps of the ICs [31,53,54], thus relying on the subjective judgment of the experimenter. Usually, this approach leads to misclassification of ICs and divergent results. To overcome these issues, we proposed a procedure to identify eye artifact-related ICs, by composite multi-scale entropy and comparing their activations during saccade and during fixation intervals, as defined by high temporal resolution eye tracking. The identified ICs are then filtered with median absolute deviation and auto-regressive exogenous model, so that the underlying neuronal activity will be preserved and distortion in the cerebral part of EEG will be minimized. Furthermore, several BCIs have been proposed in literature by using hybrid EEG- eye tracker system. Thus, the proposed algorithm can also be used with those algorithm to improve EEG signal quality which can further be used to improve classification accuracy of the BCI. Our results suggest that the proposed algorithm outperforms ICA (Figure 5, Figure 7 and Figure 9, Table 2, Table 3, Table 4 and Table 5). It can be seen in Figure 7 and Figure 8 that artifact-free EEG by the proposed algorithm for Fp1 and Fp2 are slightly different even the input EEG is same. This might be due to the ICA decomposition of the EEG signal, since the weight matrix contain different weights for different electrodes.
Urigüen and Garcia-Zapirain [55] reviewed all previous methods for removing/reducing artifacts from EEG signal and they concluded that an optimal method should consist of combining various algorithms in cascade to enhance the quality of the signal by using multiple processing stages. The idea of recovering neuronal signal from artifactual ICs was firstly proposed by Castellanos and Makarov [35]. They proposed that artifactual ICs should not be simply replaced with zero as they might have leaked neural signal in it. Recently another methodology based on the same assumption was developed by Klados and colleagues [29]. In their method, they used regression based removal of ocular activity from ICs. Although these methods proved to be effective in removing ocular artifacts from EEG signal but they require the processing of all ICs and more computational cost which is not plausible for applications like BCI. Furthermore, processing of all ICs may cause to produce distortion in those ICs which are not related to ocular artifacts and results in distortion to EEG signal. In contrast to all these methodologies, the proposed algorithm filters only the automatically identified artifactual components using median absolute deviation and auto-regressive exogenous model. In this way more neuronal activity related information can be preserved and it requires less computational cost. Results enhance our hypothesis that instead of processing all ICs only artifactual ICs should be processed (Figure 6, Figure 8 and Figure 9, Table 2, Table 3, Table 4 and Table 5).
The performance of the proposed algorithm is compared with two state-of-the-art techniques, (1) ADJUST based ICA [34] and (2) REGICA [29]. Both conventional algorithms have been implemented through EEGLAB toolbox. In the current study, composite multi-scale entropy and eye tracker based criteria are proposed for use in automatic identification of blink and eye movement related artifactual components. To verify that the proposed algorithm can differentiate between artifactual ICs and neuronal-activity related ICs, a comparison with manual detection by two experienced experts and ADJUST was carried out. The criteria used to recognize artifactual components by EEG experts was in view of time course, topographic maps and the power spectrum plots of the ICs in EEGLAB. Notably, both experts’ selection of eye movement and blink related ICs was identical. The results of this comparison was listed in Table 4. The performance of the proposed algorithm is also statistically analyzed by calculating True positive (IC marked as artifactual both by the algorithm and visual inspection), False Positive (IC marked by the algorithm but not with the visual inspection), True Negative (IC marked neither by the algorithm nor with the visual inspection), and False Negative (IC not marked with the algorithm but with the visual inspection). The total count for each parameter for all subjects is tabulated in Table 5. Average sensitivity and average specificity for all subjects is calculated as follows [27]:
S e n s i t i v i t y = T P T P + F N × 100 %
S p e c i f i c i t y = T N T N + F P × 100 %
The agreement rate between visual inspection and the proposed algorithm and ADJUST was calculated using [27]:
A g r e e m e n t R a t e = T P + T N T P + T N + F P + F N
The agreement rate between the proposed algorithm and visual inspection is found to be 98.75%, whereas for ADJUST the agreement rate is found to be 96.25%. The result of this analysis suggest that the proposed algorithm can be used as a valuable tool for automatic identification of artrifactual ICs. Relative error is a measure which depicts the performance of each algorithm in both removing ocular artifacts as well as quantifying the amount of distortion introduced in the time-domain. Our results suggest that the proposed algorithm (0.0147 ± 0.0220) has a better performance in removing ocular artifacts from EEG signal, since it successfully removes ocular artifacts, while at the same time keeping the cerebral signal intact in the time domain when it was compared to ADJUST (0.1606 ± 0.1498) and REGICA (0.0512 ± 0.0637) algorithms. Moreover, paired t-test enhances the dominance of the proposed algorithm as the difference with ADJUST method is highly statistically significant for all subjects (p < 0.001). In case of REGICA the difference is highly statistically significant with p < 0.001 except for subject 2 (p < 0.011) and subject 5 (p < 0.17). Furthermore, mutual information index was adopted to investigate that how much information artifact-free EEG signal shares with the reconstructed EEG signal after the implementation of different methods. The average mutual information value for all datasets using the proposed algorithm (2.6461) against ADJUST (1.7488) and REGICA (2.2578) demonstrate improved performance of the proposed method. This analysis enhances our hypothesis that ICA and regression method are less effective in removing ocular artifacts from EEG signal. Finally, the contribution of the different component of the proposed algorithm is analyzed by eliminating different component one by one. After eliminating different components, the performance of the proposed algorithm in removing artifacts and keeping the neuronal activity related EEG signal intact is decreased and the relative error is increased. In case of eliminating composite multi-scale entropy and eye tracker based criteria for detection of artifactual ICs cause an increase in the relative error (0.4219 ± 0.0889). This enhances our hypothesis that only artifactual ICs should be identified and processed for artifact correction. Hence the utility and value of composite multi-scale entropy and eye tracker based criteria as a useful tool for identification of eye movement and blink related ICs. Furthermore, the ability of removing ocular artifacts by the proposed algorithm is decreased by the elimination of median absolute deviation. However, the elimination of median absolute deviation has very less effect on the neuronal activity related EEG signals in nonartifactual zone because median absolute deviation is only used to remove high magnitude ocular activities. Lastly, in the proposed algorithm auto-regressive exogenous and affine projection algorithm is used to remove remaining artifacts and to compensate for the possible neuronal activity loss due to median absolute deviation filtration. The elimination of auto-regressive exogenous and affine projection algorithm cause an increase in relative error (0.0253 ± 0.0156). This analysis show that the performance of the proposed algorithm is decreased by the elimination of different components. One striking feature of the proposed algorithm is that it does not require any calibration or pre-training. Also, this is one of the few studies concerning an artifact removal technique, in which statistical analysis is used to evaluate the performance of the proposed methodology.

6. Conclusions

This paper presents a novel algorithm using hybrid EEG and eye tracker system to automatically identify and remove ocular artifacts from EEG data by combining ICA and auto-regressive exogenous model. The performance of the proposed algorithm is demonstrated using experimental and standard EEG datasets. The proposed methodology enables the removal of ocular artifacts in the artifactual zone while keeping the neuronal activity related EEG signal intact in the non-artifactual zone. Additionally, results show that the proposed algorithm outperformed the two state of the art techniques based on ADJUST and REGICA. Relative error and mutual information are used as evaluation indexes to quantify the amount of distortion produced in corrected EEG by each algorithm. The statistical significance of the proposed algorithm is verified using paired t-test.

Acknowledgments

This research was supported by the research programs of the National Research Foundation of Korea (NRF) (NRF-2015R1A5A1037668) and “Leaders in Industry-university Cooperation” Project, supported by the Ministry of Education (MOE).

Author Contributions

M.M.N.M. researched the literatures, designed the algorithm, carried out the simulations, analyzed the simulations results and wrote the manuscript. S.K. conducted the experiment and participated in data processing. M.Y.J. supervised the overall study and suggested the theoretical aspects of the study. M.A.K. contributed to the algorithm coding and participated in revising the manuscript. All authors have approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jöbsis, F.F. Noninvasive infrared monitoring of cerebral and myocardial oxygen sufficiency and circulatory parameters. Science 1977, 198, 1264–1267. [Google Scholar] [CrossRef] [PubMed]
  2. Friston, K.J.; Jezzard, P.; Turner, R. Analysis of functional MRI time-series. Hum. Brain Mapp. 1994, 1, 153–171. [Google Scholar] [CrossRef]
  3. Hogervorst, M.A.; Brouwer, A.-M.; van Erp, J.B.F. Combining and comparing EEG, peripheral physiology and eye-related measures for the assessment of mental work load. Front. Neurosci. 2014, 8. [Google Scholar] [CrossRef]
  4. Kamran, M.A.; Hong, K.-S.; Mannan, M.N.M. Identification of fNIRS based Brain Activity Using Adaptive Algorithm. NUST J. Eng. Sci. 2012, 5, 7–10. [Google Scholar]
  5. Kamran, M.A.; Jeong, M.Y.; Mannan, M.N.M. Optimal hemodynamic response model for functional near-infrared spectroscopy. Front. Behav. Neurosci. 2015, 9. [Google Scholar] [CrossRef] [PubMed]
  6. Zeng, H.; Song, A.; Yan, R.; Qin, H. EOG artifact correction from EEG recording using stationary subspace analysis and empirical mode decomposition. Sensors 2013, 13, 14839–14859. [Google Scholar] [CrossRef] [PubMed]
  7. Sameni, R.; Gouy-Pailler, C. An iterative subspace denoising algorithm for removing electroencephalogram ocular artifacts. J. Neurosci. Methods, 2014, 225, 97–105. [Google Scholar] [CrossRef] [PubMed]
  8. Berg, P.; Scherg, M. Dipole modelling of eye activity and its application to the removal of eye artefacts from the EEG and MEG. Clin. Phys. Physiol. Meas. 1991, 12, 49–54. [Google Scholar] [CrossRef] [PubMed]
  9. Wallstrom, G.L.; Kass, R.E.; Miller, A.; Cohn, J.F.; Fox, N.A. Automatic correction of ocular artifacts in the EEG: A comparison of regression-based and component-based methods. Int. J. Psychophysiol. 2004, 53, 105–119. [Google Scholar] [CrossRef] [PubMed]
  10. Dimigen, O.; Sommer, W.; Hohlfeld, A.; Jacobs, A.; Kliegl, R. Co-registration of eye movements and EEG in natural reading: Analyses and review. J. Exp. Psychol. 2011, 140, 552–572. [Google Scholar] [CrossRef] [PubMed]
  11. Carl, C.; Hipp, J.F.; König, P.; Engel, A.K. Spectral signatures of saccade target selection. Brain Topogr. 2016, 29, 130–148. [Google Scholar] [CrossRef] [PubMed]
  12. Carl, C.; Açık, A.; König, P.; Engel, A.K.; Hipp, J.F. The saccade spike artifact in MEG. Neuroimage 2012, 59, 1657–1667. [Google Scholar] [CrossRef] [PubMed]
  13. Engbert, R.; Mergenthaler, K. Microsaccades are triggered by low retinal image slip. Proc. Natl. Acad. Sci. USA 2006, 103, 7192–7197. [Google Scholar] [CrossRef] [PubMed]
  14. Fatourechi, M.; Bashashati, A.; Ward, R.K.; Birch, G.E. EMG and EOG artifacts in brain computer interface systems: A survey. Clin. Neurophysiol. 2007, 118, 480–494. [Google Scholar] [CrossRef] [PubMed]
  15. Croft, R.J.; Barry, R.J. EOG correction: A new perspective. Electroencephalogr. Clin. Neurophysiol. 1998, 107, 387–394. [Google Scholar] [CrossRef]
  16. Croft, R.J.; Barry, R.J. EOG correction: A new aligned-artifact average solution. Electroencephalogr. Clin. Neurophysiol. 1998, 107, 395–401. [Google Scholar] [CrossRef]
  17. Croft, R.J.; Barry, R.J. Removal of ocular artifact from the EEG: A review. Clin. Neurophysiol. 2000, 30, 5–19. [Google Scholar] [CrossRef]
  18. Croft, R.J.; Barry, R.J. EOG correction: Which regression should we use? Psychophysiology 2000, 37, 123–125. [Google Scholar] [CrossRef] [PubMed]
  19. Croft, R.J.; Chandler, J.S.; Barry, R.J.; Cooper, N.R.; Clarke, A.R. EOG correction: A comparison of four methods. Psychophysiology 2005, 42, 16–24. [Google Scholar] [CrossRef] [PubMed]
  20. Romero, S.; Mañanas, M.A.; Barbanoj, M.J. Ocular reduction in EEG signals based on adaptive filtering, regression and blind source separation. Ann. Biomed. Eng. 2009, 37, 176–191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Barbati, G.; Porcaro, C.; Zappasodi, F.; Rossini, P.M.; Tecchio, F. Optimization of independent component analysis approach for artifact identification and removal in MEG signals. Clin. Neurophysiol. 2004, 115, 1220–1232. [Google Scholar] [CrossRef] [PubMed]
  22. Joyce, C.A.; Gorodnitsky, I.F.; Kutas, M. Automatic removal of eye movement and blink artifacts from EEG data using blind component separation. Psychophysiology 2004, 41, 313–325. [Google Scholar] [CrossRef] [PubMed]
  23. Hoffmann, S.; Falkenstein, M. The correction of eye blink artefacts in the EEG: A comparison of two prominent methods. PLoS ONE 2008, 3. [Google Scholar] [CrossRef] [PubMed]
  24. Ghandeharion, H.; Erfanian, A. A fully automatic ocular artifact suppression from EEG data using higher order statistics: Improved performance by wavelet analysis. Med. Eng. Phys. 2010, 32, 720–729. [Google Scholar] [CrossRef] [PubMed]
  25. Javidi, S.; Mandic, D.P.; Took, C.C.; Cichocki, A. Kurtosis-based blind source extraction of complex non-circular signals with application in EEG artifact removal in real-time. Front. Neurosci. 2011, 5. [Google Scholar] [CrossRef] [PubMed]
  26. Winkler, I.; Haufe, S.; Tangermann, M. Automatic classification of artifactual ICA-components for artifact removal in EEG signals. Behav. Brain Funct. 2011, 7. [Google Scholar] [CrossRef] [PubMed]
  27. Mahajan, R.; Morshed, B.I. Unsupervised Eye Blink Artifact Denoising of EEG Data with Modified Multiscale Sample Entropy, Kurtosis, and Wavelet-ICA. IEEE J. Biomed. Health Inform. 2015, 19, 158–165. [Google Scholar] [CrossRef] [PubMed]
  28. He, P.; Wilson, G.; Russel, C. Removal of ocular artifacts from electroencephalogram by adaptive filtering. Med. Biol. Eng. Comput. 2004, 42, 407–412. [Google Scholar] [CrossRef] [PubMed]
  29. Klados, M.A.; Papadelis, C.; Braun, C.; Bamidis, P.D. REG-ICA: A hybrid methodology combining blind source separation and regression techniques for the rejection of ocular artifacts. Biomed. Signal Process. Control 2011, 6, 291–300. [Google Scholar] [CrossRef]
  30. Peng, H.; Hu, B.; Shi, Q.; Ratcliffe, M.; Zhao, Q.; Qi, Y.; Gao, G. Removal of ocular artifacts in EEG—An improved approach combining DWT and ANC for portable applications. IEEE J. Biomed. Health Inform. 2013, 17, 600–607. [Google Scholar] [CrossRef] [PubMed]
  31. Iriarte, J.; Urrestarazu, E.; Valencia, M.; Alegre, M.; Malanda, A.; Viteri, C.; Artieda, J. Independent component analysis as a tool to eliminate artifacts in EEG: A quantitative study. J. Clin. Neurophysiol. 2003, 20, 249–257. [Google Scholar] [CrossRef] [PubMed]
  32. Kong, W.Z.; Zhou, Z.P.; Hu, S.Q.; Zhang, J.H.; Babiloni, F.; Dai, G.J. Automatic and direct identification of blink components from scalp EEG. Sensors 2013, 13, 10783–10801. [Google Scholar] [CrossRef] [PubMed]
  33. Plöchl, M.; Ossandón, J.P.; König, P. Combining EEG and eye tracking: Identification, characterization, and correction of eye movement artifacts in electroencephalographic data. Front. Hum. Neurosci. 2012, 6. [Google Scholar] [CrossRef] [PubMed]
  34. Mognon, A.; Jovicich, J.; Bruzzone, L.; Buiatti, M. ADJUST: An automatic EEG artifact detector based on the joint use of spatial and temporal features. Psychophysiology 2010, 48, 229–240. [Google Scholar] [CrossRef] [PubMed]
  35. Castellanos, N.P.; Makarov, V.A. Recovering EEG brain signals: Artifact suppression with wavelet enhanced independent component analysis. J. Neurosci. Methods 2006, 158, 300–312. [Google Scholar] [CrossRef] [PubMed]
  36. Ehinger, B.V.; König, P.; Ossandón, J.P. Predictions of visual content across eye movements and their modulation by inferred information. J. Neurosci. 2015, 35, 7403–7413. [Google Scholar] [CrossRef] [PubMed]
  37. Kierkels, J.; Riani, J.; Bergmans, J. Using an eye tracker for accurate eye movement artifact correction. IEEE Trans. Biomed. Eng. 2007, 54, 1257–1267. [Google Scholar] [CrossRef] [PubMed]
  38. Noureddin, B.; Lawrence, P.D.; Birch, G.E. Online removal of eye movement and blink EEG artifacts using a high-speed eye tracker. IEEE Trans. Biomed. Eng. 2012, 59, 2103–2110. [Google Scholar] [CrossRef] [PubMed]
  39. EYE-EEG: Eye tracking & EEG. Available online: http://www2.hu-berlin.de/eyetracking-eeg/testdata.html (accessed on 20 October 2015).
  40. Delorme, A.; Sejnowski, T.; Makeig, S. Enhanced detection of artifacts in EEG data using higher-order statistics and independent component analysis. Neuroimage 2007, 34, 1443–1449. [Google Scholar] [CrossRef] [PubMed]
  41. Delorme, A.; Makeig, S. EEGLAB: An open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 2004, 134, 9–21. [Google Scholar] [CrossRef] [PubMed]
  42. Costa, M.; Goldberger, A.L.; Peng, C.K. Multiscale entropy analysis of biological signals. Phys. Rev. E 2005, 71. [Google Scholar] [CrossRef] [PubMed]
  43. Bosl, W.; Tierney, A.; Tager-Flusberg, H.; Nelson, H. EEG complexity as a biomarker for autism spectrum disorder risk. BMC Med. 2011, 9. [Google Scholar] [CrossRef] [PubMed]
  44. Wu, S.-D.; Wu, C.-W.; Lin, S.-G.; Wang, C.-C.; Lee, K.-Y. Time Series Analysis Using Composite Multiscale Entropy. Entropy 2013, 15, 1069–1084. [Google Scholar] [CrossRef]
  45. Leys, C.; Ley, C.; Klein, O.; Bernard, P.; Licata, L. Detecting outliers: Do not use standard deviation around the mean, use absolute deviation around the median. J. Exp. Soc. Psychol. 2013, 49, 764–766. [Google Scholar] [CrossRef]
  46. Kamran, M.A.; Hong, K.-S. Linear parameter-varying model and adaptive filtering technique for detecting neuronal activities: An fNIRS study. J. Neural Eng. 2013, 10. [Google Scholar] [CrossRef] [PubMed]
  47. Wang, Z.; Xu, P.; Liu, T.; Tian, Y.; Lei, X.; Yao, D. Robust removal of ocular artifacts by combining independent component analysis and system identification. Biomed. Signal Process. Control 2014, 10, 250–259. [Google Scholar] [CrossRef]
  48. Swartz Center for Computational Neuroscience. Available online: http://sccn.ucsd.edu/~jason/minfo.m (accessed on 1 November 2015).
  49. Stone, J.V. Independent component analysis: An introduction. Trends Cogn. Sci. 2002, 6, 59–64. [Google Scholar] [CrossRef]
  50. Tran, Y.; Craig, A.; Boord, P.; Craig, D. Using independent component analysis to remove artifact from electroencephalographic measured during stuttered speech. Med. Bio. Eng. Comput. 2004, 42, 627–633. [Google Scholar] [CrossRef]
  51. Bian, N.Y.; Wang, B.; Cao, Y. Automatic removal of artifacts from EEG data using ICA and nonlinear exponential analysis. Acta Biophys. Sin. 2006, 22, 149–156. [Google Scholar]
  52. Mammone, N.; Morabito, F.C. Enhanced automatic artifact detection based on independent component analysis and Renyi’s entropy. Neural Netw. 2008, 21, 1029–1040. [Google Scholar] [CrossRef] [PubMed]
  53. Jung, T.P.; Makeig, S.; Humphries, C.; Lee, T.W.; McKeown, M.J.; Iragui, V.; Sejnowski, T.J. Removing electroencephalographic artifacts by blind source separation. Psychophysiology 2000, 37, 163–178. [Google Scholar] [CrossRef] [PubMed]
  54. Makeig, S.; Onton, J. ERP Features and EEG Dynamics: An ICA Perspective. In Oxford Handbook of Event-Related Potential Components; Luck, S.J., Kappenman, E.S., Eds.; Oxford University Press: New York, NY, USA, 2009; pp. 51–87. [Google Scholar]
  55. Urigüen, J.A.; Garcia-Zapirain, B. EEG artifact removal—State-of-the-art and guidelines. J. Neural Eng. 2015, 12. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic diagrams. (A) Regression method; (B) Independent component analysis.
Figure 1. Schematic diagrams. (A) Regression method; (B) Independent component analysis.
Sensors 16 00241 g001
Figure 2. Schematic diagram of the proposed algorithm.
Figure 2. Schematic diagram of the proposed algorithm.
Sensors 16 00241 g002
Figure 3. (A) Electroencephalogram (EEG) electrode configuration; (B) Distribution of saccade amplitude.
Figure 3. (A) Electroencephalogram (EEG) electrode configuration; (B) Distribution of saccade amplitude.
Sensors 16 00241 g003
Figure 4. Results on experimental dataset. (A) Experimental EEG data for one subject; (B) Independent components (ICs) obtained from independent component analysis (ICA) decomposition of EEG data; (C) Comparison of the corrected EEG by the proposed algorithm and conventional algorithms.
Figure 4. Results on experimental dataset. (A) Experimental EEG data for one subject; (B) Independent components (ICs) obtained from independent component analysis (ICA) decomposition of EEG data; (C) Comparison of the corrected EEG by the proposed algorithm and conventional algorithms.
Sensors 16 00241 g004
Figure 5. Comparison of the proposed algorithm with ADJUST using experimental data. (A) Contaminated experimental EEG data at Fp1 and Fp2; (B) Corrected EEG by the proposed algorithm; (C) Corrected EEG by ADJUST; (D) Comparison of corrected EEG with contaminated experimental EEG. (E) Partial enlargement of highlighted regions.
Figure 5. Comparison of the proposed algorithm with ADJUST using experimental data. (A) Contaminated experimental EEG data at Fp1 and Fp2; (B) Corrected EEG by the proposed algorithm; (C) Corrected EEG by ADJUST; (D) Comparison of corrected EEG with contaminated experimental EEG. (E) Partial enlargement of highlighted regions.
Sensors 16 00241 g005
Figure 6. Comparison of the proposed algorithm with REGICA using experimental data. (A) Contaminated experimental EEG data at Fp1 and Fp2; (B) Corrected EEG by the proposed algorithm; (C) Corrected EEG by REGICA; (D) Comparison of corrected EEG with contaminated experimental EEG; (E) Partial enlargement of highlighted regions.
Figure 6. Comparison of the proposed algorithm with REGICA using experimental data. (A) Contaminated experimental EEG data at Fp1 and Fp2; (B) Corrected EEG by the proposed algorithm; (C) Corrected EEG by REGICA; (D) Comparison of corrected EEG with contaminated experimental EEG; (E) Partial enlargement of highlighted regions.
Sensors 16 00241 g006
Figure 7. Comparison results of the proposed algorithm and ADJUST for all subjects at Fp1 and Fp2.
Figure 7. Comparison results of the proposed algorithm and ADJUST for all subjects at Fp1 and Fp2.
Sensors 16 00241 g007
Figure 8. Comparison results of the proposed algorithm and REGICA for all subjects at Fp1 and Fp2.
Figure 8. Comparison results of the proposed algorithm and REGICA for all subjects at Fp1 and Fp2.
Sensors 16 00241 g008
Figure 9. Comparison of the proposed algorithm with ADJUST and REGICA in frequency domain at Fp1 and Oz. (A) EEG spectra after applying filter 0.5–40 Hz; (B) EEG spectra after applying filter 0.5–20 Hz.
Figure 9. Comparison of the proposed algorithm with ADJUST and REGICA in frequency domain at Fp1 and Oz. (A) EEG spectra after applying filter 0.5–40 Hz; (B) EEG spectra after applying filter 0.5–20 Hz.
Sensors 16 00241 g009aSensors 16 00241 g009b
Figure 10. Comparison of the proposed algorithm with ADJUST using standard data. (A) Contaminated experimental EEG data at Fp1 and Fp2. (B) Corrected EEG by the proposed algorithm. (C) Corrected EEG by ADJUST. (D) Comparison of Corrected EEG with contaminated experimental EEG. (E) Partial enlargement of highlighted regions.
Figure 10. Comparison of the proposed algorithm with ADJUST using standard data. (A) Contaminated experimental EEG data at Fp1 and Fp2. (B) Corrected EEG by the proposed algorithm. (C) Corrected EEG by ADJUST. (D) Comparison of Corrected EEG with contaminated experimental EEG. (E) Partial enlargement of highlighted regions.
Sensors 16 00241 g010
Table 1. Summary of the proposed algorithm.
Table 1. Summary of the proposed algorithm.
Input: Contaminated EEG data, Eye tracker data
Output: Artifact-free EEG data
Synchronization of EEG and eye tracker
Decompose contaminated EEG data using ICA to get ICs
Portioning of ICs into saccade and fixation epochs
Calculate composite multi-scale entropy and ratio of mean variance to identify ocular artifacts related ICs
Apply median absolute deviation to remove high magnitude ocular activities from identified ICs
Filter ICs with auto-regressive exogenous model and affine projection algorithm
Artifact-free EEG data by back projecting all ICs using inverse ICA
Table 2. Relative error values results for the proposed algorithm against ADJUST and REGICA using artifact-free EEG data.
Table 2. Relative error values results for the proposed algorithm against ADJUST and REGICA using artifact-free EEG data.
SubjectProposedADJUST ICAp-valREGICAp-val
10.0001 ± 0.00010.0645 ± 0.0944<0.0010.0452 ± 0.0613<0.001
20.0185 ± 0.03050.1741 ± 0.1417<0.0010.0437 ± 0.0393<0.011
30.0146 ± 0.02300.1025 ± 0.1191<0.0010.0527 ± 0.06–<0.005
40.0134 ± 0.02440.1893 ± 0.1582<0.0010.0779 ± 0.1140<0.005
50.0273 ± 0.03200.2727 ± 0.2358<0.0010.0366 ± 0.0389<0.17
Average0.0147 ± 0.02200.1606 ± 0.1498 0.0512 ± 0.0637
Table 3. Average mutual information values results of all subjects for the proposed algorithm against ADJUST and REGICA for all electrodes.
Table 3. Average mutual information values results of all subjects for the proposed algorithm against ADJUST and REGICA for all electrodes.
Electrode LocationProposedADJUSTREGICA
Fp12.71480.66801.4023
Fp22.62860.75781.4563
F72.71691.17961.9115
F32.71651.27892.0176
Fz2.69731.39321.9717
F22.64341.37892.0190
F82.55381.42371.8863
FC52.69461.29192.3216
FC32.52621.56922.3181
FC22.57201.87422.3239
FC62.60341.80462.1482
T72.71581.45712.3656
C32.66851.47542.3942
Cz2.74352.06692.4309
C42.60751.98132.3062
T82.66111.85242.3069
TP92.65151.58512.3234
CP52.68881.55662.4402
CP12.66692.09102.4864
CP22.71862.15922.4384
CP62.58392.02362.3601
TP102.63071.95412.3272
P72.56001.85882.4463
P32.62752.17762.4715
Pz2.60762.18112.3613
P22.62302.17392.3948
P82.62902.07962.4958
PO92.65402.21942.4626
O12.65972.13802.4375
Oz2.67582.25842.4655
O22.58361.91912.3257
PO102.65142.13372.4338
Average2.64611.74882.2578
Table 4. Comparison results for identification of artifactual ICs.
Table 4. Comparison results for identification of artifactual ICs.
SubjectEEG ExpertsPropsoed AlgorithmADJUST
Vertical and HorizontalBlinkVertical and HorizontalBlinkVertical and HorizontalBlink
1111101
2111202
3122222
4212121
5111131
Table 5. Performance evaluation of the proposed algorithm and ADJUST for all subjects.
Table 5. Performance evaluation of the proposed algorithm and ADJUST for all subjects.
MethodTrue Positive (TP)False Positive (FP)True Negative (TN)False Negative (FN)Average SensitivityAverage Specificity
Proposed1221460100%98.64%
ADJUST104144283.33%97.29%

Share and Cite

MDPI and ACS Style

Mannan, M.M.N.; Kim, S.; Jeong, M.Y.; Kamran, M.A. Hybrid EEG—Eye Tracker: Automatic Identification and Removal of Eye Movement and Blink Artifacts from Electroencephalographic Signal. Sensors 2016, 16, 241. https://doi.org/10.3390/s16020241

AMA Style

Mannan MMN, Kim S, Jeong MY, Kamran MA. Hybrid EEG—Eye Tracker: Automatic Identification and Removal of Eye Movement and Blink Artifacts from Electroencephalographic Signal. Sensors. 2016; 16(2):241. https://doi.org/10.3390/s16020241

Chicago/Turabian Style

Mannan, Malik M. Naeem, Shinjung Kim, Myung Yung Jeong, and M. Ahmad Kamran. 2016. "Hybrid EEG—Eye Tracker: Automatic Identification and Removal of Eye Movement and Blink Artifacts from Electroencephalographic Signal" Sensors 16, no. 2: 241. https://doi.org/10.3390/s16020241

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop