Computational and Mathematical Methods in Medicine

Volume 2013 (2013), Article ID 591216, 7 pages

http://dx.doi.org/10.1155/2013/591216

Research Article

Local Temporal Correlation Common Spatial Patterns for Single Trial EEG Classification during Motor Imagery

Rui Zhang, Peng Xu, Tiejun Liu, Yangsong Zhang, Lanjin Guo, Peiyang Li, and Dezhong Yao

Key Laboratory for Neuroinformation of Ministry of Education, School of Life Science and Technology, University of Electronic Science and Technology of China, Chengdu 610054, China

Received 1 August 2013; Revised 10 October 2013; Accepted 10 October 2013

Academic Editor: Shiliang Sun

Copyright © 2013 Rui Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Common spatial pattern (CSP) is one of the most popular and effective feature extraction methods for motor imagery-based brain-computer interface (BCI), but the inherent drawback of CSP is that the estimation of the covariance matrices is sensitive to noise. In this work, local temporal correlation (LTC) information was introduced to further improve the covariance matrices estimation (LTCCSP). Compared to the Euclidean distance used in a previous CSP variant named local temporal CSP (LTCSP), the correlation may be a more reasonable metric to measure the similarity of activated spatial patterns existing in motor imagery period. Numerical comparisons among CSP, LTCSP, and LTCCSP were quantitatively conducted on the simulated datasets by adding outliers to Dataset IVa of BCI Competition III and Dataset IIa of BCI Competition IV, respectively. Results showed that LTCCSP achieves the highest average classification accuracies in all the outliers occurrence frequencies. The application of the three methods to the EEG dataset recorded in our laboratory also demonstrated that LTCCSP achieves the highest average accuracy. The above results consistently indicate that LTCCSP would be a promising method for practical motor imagery BCI application.

1. Introduction

Brain-computer interfaces (BCI) use electroencephalographic signals or other electrophysiological measures of brain activity to provide a new nonmuscular channel for sending messages and commands to the external world. According to the different electrophysiological signals which they use, BCI can fall into 6 groups [1]: visual evoked potentials (VEP) based BCI [2]; slow cortical potentials (SCP) based BCI [3]; evoked potentials P300 based BCI [4]; mu and beta rhythms (ERD/ERS) based BCI [5], cortical neuronal action potentials based BCI [6] and hybrid BCI [7]. Among them, ERD/ERS based BCI has received a lot of attentions in recent years due to its potential application in motor rehabilitation and its assisting for the motor function impaired patients [8–10].

Feature extraction and classification algorithms play important roles for the performance of ERD/ERS based BCI, and there are various methods have been proposed to extract ERD/ERS related features [11, 12], such as the laplacian method [13], autoregressive spectral analysis [14], common spatial pattern (CSP) [15], discriminative spatial patterns [16], bispectrum analysis [17], and multivariate empirical mode decomposition [18]. Currently, CSP is one of the most popular feature extraction methods for ERD/ERS based BCI, its efficiencies have been proved by the BCI competitions [19, 20] and various ERD/ERS based BCI studies [21–23]. Besides BCI discipline, CSP has been also applied to normal versus abnormal EEGs classification [24] and EEG source localizations [25].

The aim of CSP is to learn the optimal spatial filters which maximize the variance of one class while minimizing the variance of the other class simultaneously [26, 27]. Mathematically, CSP relies on the simultaneous diagonalization of two covariance matrices. However, there are two inherent drawbacks for the estimation of the covariance matrices in using the conventional strategy; on one hand, it is prone to be influenced by outlier noise; even one outlier may make the spatial filters obtained meaningless [28]; on the other hand, the temporal information is neglected, while the time-dependent local variances may deliver more discriminant power than the global variances [29]. To deal with above shortcomings in conventional CSP, Local Temporal Common Spatial Patterns (LTCSP) propose to use the Euclidean distance between and (both and are -channel EEG recording vectors at two given time points and ) as a weight to emphasize the covariance matrices [29]. From the perspective of neurophysiology, the task related signal is usually hidden in strong spontaneous brain background activity, which is nonstationary  [30, 31]. Therefore, the scalp-recorded signal is very noisy even if the recordings are from high-performance EEG amplifier with suitable reference strategy [32, 33] and in a finely shielded room with a very cooperative subject. In this case, the Euclidean distance measure may give unsuitable weight coefficient due to the noise effect. Besides, the Euclidean distance is not convergent with distribution in a wide range, which may reduce the robustness of LTCSP. In this work, the correlation measure was newly proposed to replace the Euclidean distance in order to solve the above problems. In fact, when Müller-Gerking et al. firstly introduced CSP into BCI areas [26], the correlations of neighboring electrodes had been argued as an important information which implies the need to pay attention for them.

The framework of this paper is arranged as follows. Section 2 gives details about LTCSSP; In Section 3, CSP, LTCSP, and LTCCSP are evaluated on simulated outliers influenced BCI Competition datasets. Then the comparison results of the above three methods on the actual EEG dataset recorded in our lab are given in Section 4; Sections 5 and 6 include discussion and final conclusion, respectively. Besides the above sections, an appendix was provided to describe the mathematical details of CSP.

2. Materials and Methods

2.1. Principles of LTCCSP

LTCCSP is an extension of the conventional CSP (See Appendix for detail). The spatial filter matrix of conventional CSP is obtained by maximizing the variance of one class while minimizing the variance of the other class simultaneously. We denote matrix and as the EEG data under task is 1 and 2, with    being the number of channels and being the number of samples in each trial. Formally, the object function of CSP could be expressed as [34]

where and are the average normalized spatial covariance matrices of tasks 1 and 2 and is the spatial filter matrix. Generally, the first and last few columns of are served as optimized spatial filters.

Then taking the first columns of for an example, (1) could be transformed to

where and are, respectively, the numbers of trials under each task, denotes the th trial, and is the th column of the matrix . The last columns could be transformed as above, too.

Using the dimension reduction strategy in [35], the quadratic forms   could be transformed to

where and are the data vectors at time points and , respectively. Obviously, (3) mainly focuses on the global information of trials. However, in actual BCI situation, the local information may be helpful for task recognition. Therefore, the weight matrix can be added to emphasize the local information as follows [29]:

where is the weight matrix. The fundamental derivation of weight matrix is that if the two concerned patterns and are close, it will impose a large coefficient in the weight matrix. Hence, we will use the correlation coefficient to define the corresponding weight matrix as follows:

where denotes the correlation coefficient operator and is the local temporal range. Obviously, the weight matrix will impose a relatively larger coefficient on the similar patterns in the concerned local temporal range .

Moreover, (4) will be further expressed as

where the Laplacian matrix and is a diagonal matrix with the diagonal elements being the row sums of   ; that is, . Equation (6) is also held for EEG trials under task 2. Using (6), (2) can be converted as follows:

where and are used to replace    and , respectively. Let and be the average normalized local temporal correlation covariance matrices of two classes; (7) can be simplified to

Noting that (8) has the same structure as (1), so the solution of maximizing (8) subject to is similar to (1), and the spatial filter matrix could be given by

where is the matrix of eigenvectors of with being the diagonal matrix of associated eigenvalues and is the eigenvectors matrix of .

2.2. Online Implementation of LTCCSP

In the online application, when the spatial filter matrix is obtained from the training dataset, there also exists a little difference between LTCCSP and conventional CSP for feature extraction. Take as a generic symbol for and ; since is a semipositive definite matrix, it can be decomposed as . Then EEG data   from a trial  should be spatially filtered as . At last, the logarithm variances of the first and last   rows of are served as final features and will be sent to classifier for classification.

3. Evaluations on Simulated Outliers Influenced EEG Datasets

3.1. EEG Datasets Description

(1) Dataset IVa of BCI Competition III. This dataset contains EEG signals recorded from five subjects by using 118 electrodes [36]. In each trial, a visual cue was shown for 3.5 s, during which three kinds of motor imageries were performed, that is, left hand, right hand, and right foot. The motor imageries of right hand and foot were needed to be classified. The total number of EEG trials for each subject is 280. In particular, 168, 224, 84, 56, and 28 trials are used as training data corresponding to the five subjects: aa, al, av, aw, and ay, respectively. The data were band-pass filtered between 0.05 and 200 Hz and downsampled at 100 Hz for subsequent analysis.

(2) Dataset IIa of BCI Competition IV. The EEG data was recorded from nine subjects, who were asked to perform four different kinds of motor imagery tasks, that is, left hand, right hand, foot, and tongue [37]. Two sessions recorded on different days were made available for each individual, and each session consisted of 288 trials with 72 trials per class. The EEG signals, measured by 22 electrodes, were sampled at 250 Hz. In our experiment, one session containing the class labels for all trials is used as training set and the unlabelled session as test set. We focus on two-class classification scenario, that is, classifying EEG signals belonging to left and right hand motor imageries.

3.2. Preprocessing

Following the winner of BCI competition IV and [34], the EEG segments recorded from 0.5 s to 3.75 s after the visual cue were used for subsequent analysis on the first dataset, and the time interval from 0.5 s to 2.5 s was used on the second dataset. Then the EEG segments were band-pass filtered between 8 Hz and 30 Hz.

3.3. Introducing Outliers

We introduced outliers into the two datasets from BCI Competitions with the aim to simulate the actual recordings with outliers contamination and then quantitatively compare the classification accuracies of the conventional CSP and LTCSP methods and the proposed LTCCSP method. For each subject, the outliers were generated by the one-dimensional Gaussian distribution , where and are the mean and standard deviation of the training EEG segments, respectively. The generated outliers were added to the training EEG data to simulate polluted EEG signals. The number of outliers is varying from 0 to with step , where is the number of trials in training set. We randomly selected the time positions to add the outliers.

3.4. Results

We used CSP, LTCSP, and LTCCSP to extract the motor imagery-related features, respectively. Comparing with CSP, there are two parameters in LTCSP and one parameter in LTCCSP that need to be configured. We set in using LTCSP, where is the standard deviation of the squared norms of the training samples, as the highest average classification accuracy was achieved under this parameter [29]. For the local temporal range that exists in both LTCSP and LTCCSP, we used 10-fold cross-validation method to compare the average accuracies of within on the training set and then set as the one corresponding to the highest average accuracy.

Three pairs of spatial filters were used for feature extraction for all of the three CSP-based methods, as recommended in [15]. Then the log-variances of the spatially filtered EEG signals were used as input features for a classifier. Support Vector Machine (SVM), one of the most popular classifiers in BCI application, was used for classification [38], and the default parameters were set for SVM. The experiments were repeated ten times for each occurrence frequency of outliers; the average accuracies and the mean accuracies for all of the subjects were reported in Table 1.

Table 1: Average classification accuracies (%) of the Dataset IVa of BCI Competition III and Dataset IIa of BCI Competition IV for CSP, LTCSP, and LTCCSP with increasing occurrence frequencies of outliers.

Freq

Dataset IVa of BCI Competition III
Dataset IIa of BCI Competition IV

Subject
aa
al
av
aw
ay
1
2
3
4
5
6
7
8
9
Mean

0
CSP

LTCSP

LTCCSP

0.1
CSP

LTCSP

LTCCSP

0.2
CSP

LTCSP

LTCCSP

0.3
CSP

LTCSP

LTCCSP

0.4
CSP

LTCSP

LTCCSP

Paired t-test between two concerned methods, that is, LTCSP versus CSP and LTCCSP versus CSP. The bold values indicate the best performance among the three methods.

All three CSP-based methods achieved the highest classification accuracies on the clean EEG data compared to those cases with various occurrence frequencies of outliers added. When outliers were introduced, the performance of all the three methods was deteriorated with the increasing occurrence frequencies of outliers for most of the subjects except some subjects, that is, Sub ay, Sub 1, and Sub 2. For Sub 2, all the three methods achieved approximately 50% accuracies, which may indicate BCI illiteracy phenomenon [39], that is, the subject cannot control BCI effectively by motor imagery.

Table 1 also revealed that both LTCCSP and LTCSP achieve better performance than the conventional CSP under various outliers conditions, where LTCCSP showed the statistically significant improvement but not for LTCSP. Among the three methods, the best performance was achieved by LTCCSP on both clean EEG datasets and outliers influenced EEG datasets.

4. Evaluation on Actual EEG Dataset

4.1. EEG Datasets Description

Dataset Recorded in Our Laboratory. This dataset consists of EEG data from 13 subjects (11 males and 2 females, right handed, 19–25 years old). The experimental protocol was approved by the Institution Research Ethics Board at University of Electronic Science & Technology of China. All participants were asked to read and sign an informed consent form before participating in the study. After experiment, all the participants received a monetary compensation for their time and effort. Subjects sat in a comfortable armchair in front of a computer screen; they were asked to perform motor imagery with left hand or right hand according to the instructions appeared on the screen. Motor imagery lasted for 5 seconds, and followed a 5 seconds rest. Fifteen Ag/AgCl electrodes covers sensorimotor area were used for EEG recordings with Symtop Amplifier (Symtop Instrument, Beijing, China), and the signals were sampled with 1000 Hz and band-pass filtered between 0.5 Hz and 45 Hz. Four runs were collected; each run consisted of 50 trials, 25 trials for each class, and there was a 3-minute break between the consecutive two runs. The first two runs are served as training dataset, and the last two runs are served as test dataset.

4.2. Preprocessing

All the EEG segments during motor imagery period were selected, that is, from 0 s to 5 s. Then the selected EEG segments were band-pass filtered between 8 Hz and 30 Hz.

4.3. Results

We used CSP, LTCSP, and LTCCSP to extract the motor imagery related features, respectively, and the parameters of LTCCSP and LTCSP were the same as we mentioned in Section 3.4. SVM with the default parameters was used for classification. The classification accuracy for each subject and the mean accuracy were reported in Table 2.

Table 2: Classification accuracies (%) of the Dataset recorded in our laboratory for CSP, LTCSP, and LTCCSP.

Subject
TCY
LPY
GK
WCF
WZQ
CR
HYR
ZB
PKH
FNX
XJP
WXY
DT
Mean

CSP

LTCSP

LTCCSP

Paired t-test between two concerned methods, that is, LTCSP versus CSP and LTCCSP versus CSP. The bold values indicate the best performance among the three methods.

On the EEG dataset recorded in our lab, LTCCSP achieved the highest mean accuracy across the three methods. As revealed in Table 2, the highest accuracies were achieved by LTCCSP in 11 of 13 subjects. Furthermore, when we set the accuracies obtained by CSP as the baseline; the accuracies obtained by LTCCSP are higher or at least equal to the baseline, while the accuracies obtained by LTCSP are not consistently improved across the 13 subjects. Paired -test showed that there is statistical difference between CSP and LTCCSP, while no statistical difference existed between CSP and LTCSP.

5. Discussions

Both LTCCSP and LTCSP introduce local temporal information to the covariance matrices estimation procedures in CSP, and our experiments indicate that the local temporal information is effective in handling outliers. If we take the performance of CSP as baseline, Table 1 indicates that LTCCSP consistently improves the performance for most of subjects under outliers condition, but the improvement of LTCSP under the same situation is not stable. Table 1 also tells us that the performance of LTCSP is worse than that of CSP on the two clean BCI Competition datasets without outliers introduced, but Table 2 shows us that LTCSP achieves better performance than CSP on the EEG dataset recorded in our lab.

We think that the performance difference between LTCCSP and LTCSP is mainly caused by the definition of the weight matrix. LTCCSP uses correlation measure to construct the weight matrix, while LTCSP uses the Euclidean distance to construct the weight matrix. In our opinion, the Euclidean distance measure is not good enough because sometimes it will wrongly treat normal signal as noise signal. An extreme case is that, for example, and , where , are EEG data vectors selected from time points and , respectively. When we calculate the Euclidean distance of them, we will get the Euclidean norm of the constant vector. Obviously, if the constant vector is of a large norm, the Euclidean distance of them will be large too, then, in the framework of LTCSP the corresponding EEG information will be strongly suppressed by imposing a very small weight on the estimation of the covariance matrix. However, no matter how big is the norm of the constant vector, the correlation coefficient of them is fixed at 1. In practical application, similar situations may appear in natural EEG signals due to the transient baseline shift. In addition, since the crucial concern of CSP is based on the spatial pattern [26], the data vectors that have similar spatial patterns should be imposed with a relatively large weight. Therefore, correlation coefficient may be a more suitable measure to construct the weight matrix.

In summary, compared with the Euclidean distance-based LTCSP, correlation-based LTCCSP may provide the following merits.

(1)The correlation coefficient distributes in , while the Euclidean distance distributes in a wide range, which makes LTCCSP more stable than LTCSP.(2)Use of correlation may reduce the possibility of imposing small weight to natural EEG signal. (3)There is only one parameter in LTCCSP, while in LTCSP there are two parameters need to configure, and , which makes LTCCSP easier to implement in practical use.

It should be noted that several methods have also been proposed to improve CSP’s performance. For example, adaptive CSP is developed to deal with time-varying signals [40]; nonlinear CSP overcomes the linearity restriction [41]; L1 norm CSP enhances its robustness [42]. Besides, an excellent method named extreme energy ratio (EER) has been recently proposed [43, 44], which also relies on covariance matrices estimation and eigenvalue decomposition, but it aims at maximizing or minimizing the disparity of energy features between two classes of EEG signals. Afterwards, semisupervised EERs are developed to solve the small training set and time-varying problems in BCI application [45]. The algorithm proposed in this paper is mainly to establish robust covariance matrices estimation by introducing local temporal correlation information, and it could be complementary to the above existing methods to develop the more competitive method in the future.

6. Conclusions

In this paper, we proposed a practical feature extraction method named LTCCSP, which considers local temporal correlation information in the learning process of the conventional CSP, for optimizing spatial filters. The current results confirmed that LTCCSP has the ability to obtain meaningful spatial filters from natural EEG data and noise influenced EEG data. Furthermore, LTCCSP is simple for application as there is only one parameter that needs to be configured. However, it should be noted that though LTCCSP shows its effective ability for motor imagery-related feature extraction, it still cannot completely suppress the introduced noise, and much work is still needed to promote for a more robust feature extraction.

Appendix

The Method of Common Spatial Patterns

We describe here the mathematical part of the method of CSP as used in the present paper [26].

Let matrix denote the filtered data of a trial under task 1, with being the number of channels and being the number of samples in time. Thus, the recording at a given time point can be represented as a point in -dimensional Euclidean space and also can be seen as a spatial pattern. The normalized spatial covariance matrix is calculated as

Likewise, let denote the corresponding normalized spatial covariance matrix of a trial under task 2. Then, the normalized spatial covariance matrices are averaged over trials, and are obtained for each task. Next, whiten the composite spatial covariance matrix ; that is, determine a matrix such that

where is an identity matrix. The whitening matrix is formed as , where is the eigenvectors matrix of and is the diagonal matrix of associated eigenvalues. After that, let and , respectively; then and share the same eigenvectors matrix :

This decomposition can be accomplished due to . Make the eigenvalues mentioned are sorted in descending order; in consequence, the final optimal spatial filter is given by

Using this projection matrix , the EEG signals and are projected as

Since the sum of the corresponding eigenvalues is always one, the variances of first and last few rows of and are suitable features for classification.

Acknowledgments

This work was supported in part by Grants from 973 programs 2011CB707803, the National Nature Science Foundation of China (no. 61175117, no. 31070881, and no. 31100745), the program for New Century Excellent Talents in University (no. NCET-12-0089), the 863 project 2012AA011601, and the National Science & Technology Pillar Program 2012BAI16B02.

References

J. R. Wolpaw, N. Birbaumer, D. J. McFarland, G. Pfurtscheller, and T. M. Vaughan, “Brain-computer interfaces for communication and control,” Clinical Neurophysiology, vol. 113, no. 6, pp. 767–791, 2002.

Y. Zhang, P. Xu, T. Liu, J. Hu, R. Zhang, and D. Yao, “Multiple frequencies sequential coding for SSVEP-based brain-computer interface,” PLoS ONE, vol. 7, no. 3, Article ID e29519, 2012.

N. Birbaumer, N. Ghanayim, T. Hinterberger et al., “A spelling device for the paralysed,” Nature, vol. 398, no. 6725, pp. 297–298, 1999.

S. Halder, E. M. Hammer, S. C. Kleih, et al., “Prediction of auditory and visual p300 brain-computer interface aptitude,” PLoS ONE, vol. 8, Article ID e53513, 2013.

J. R. Wolpaw and D. J. McFarland, “Control of a two-dimensional movement signal by a noninvasive brain-computer interface in humans,” Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 51, pp. 17849–17854, 2004.

G. Schalk, K. J. Miller, N. R. Anderson et al., “Two-dimensional movement control using electrocorticographic signals in humans,” Journal of Neural Engineering, vol. 5, no. 1, pp. 75–84, 2008.

Y. Li, J. Long, T. Yu et al., “An EEG-based BCI system for 2-D cursor control by combining Mu/Beta rhythm and P300 potential,” IEEE Transactions on Biomedical Engineering, vol. 57, no. 10, pp. 2495–2505, 2010.

B. Varkuti, C. Guan, Y. Pan, et al., “Resting state changes in functional connectivity correlate with movement recovery for BCI and robot-assisted upper-extremity training after stroke,” Neurorehabilitation and Neural Repair, vol. 27, pp. 53–62, 2013.

M. Takahashi, K. Takeda, Y. Otaka, et al., “Event related desynchronization-modulated functional electrical stimulation system for stroke rehabilitation: a feasibility study,” Journal of NeuroEngineering and Rehabilitation, vol. 9, p. 56, 2012.

S. R. Soekadar, M. Witkowski, J. Mellinger, A. Ramos, N. Birbaumer, and L. G. Cohen, “ERD-based online brain-machine interfaces (BMI) in the context of neurorehabilitation: optimizing BMI learning and performance,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 19, no. 5, pp. 542–549, 2011.

F. Lotte, M. Congedo, A. Lécuyer, F. Lamarche, and B. Arnaldi, “A review of classification algorithms for EEG-based brain-computer interfaces,” Journal of Neural Engineering, vol. 4, no. 2, article, pp. R1–R13, 2007.

A. Bashashati, M. Fatourechi, R. K. Ward, and G. E. Birch, “A survey of signal processing algorithms in brain-computer interfaces based on electrical brain signals,” Journal of Neural Engineering, vol. 4, no. 2, article, pp. R32–R57, 2007.

D. J. McFarland, L. M. McCane, S. V. David, and J. R. Wolpaw, “Spatial filter selection for EEG-based communication,” Electroencephalography and Clinical Neurophysiology, vol. 103, no. 3, pp. 386–394, 1997.

D. J. McFarland and J. R. Wolpaw, “Sensorimotor rhythm-based brain-computer interface (BCI): model order selection for autoregressive spectral analysis,” Journal of Neural Engineering, vol. 5, no. 2, pp. 155–162, 2008.

B. Blankertz, R. Tomioka, S. Lemm, M. Kawanabe, and K.-R. Müller, “Optimizing spatial filters for robust EEG single-trial analysis,” IEEE Signal Processing Magazine, vol. 25, no. 1, pp. 41–56, 2008.

X. Liao, D. Yao, D. Wu, and C. Li, “Combining spatial filters for the classification of single-trial EEG in a finger movement task,” Transactions on Biomedical Engineering, vol. 54, pp. 821–831, 2007.

S. Shahid and G. Prasad, “Bispectrum-based feature extraction technique for devising a practical brain-computer interface,” Journal of Neural Engineering, vol. 8, no. 2, Article ID 025014, 2011.

C. Park, D. Looney, N. Rehman, A. Ahrabian, and D. P. Mandic, “Classification of motor imagery BCI using multivariate empirical mode decomposition,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 21, no. 1, pp. 10–22, 2012.

B. Blankertz, K.-R. Müller, G. Curio et al., “The BCI competition 2003: progress and perspectives in detection and discrimination of EEG single trials,” IEEE Transactions on Biomedical Engineering, vol. 51, no. 6, pp. 1044–1051, 2004.

B. Blankertz, K.-R. Müller, D. J. Krusienski et al., “The BCI competition III: validating alternative approaches to actual BCI problems,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 14, no. 2, pp. 153–159, 2006.

P. Xu, P. Yang, X. Lei, and D. Yao, “An enhanced probabilistic LDA for multi-class brain computer interface,” PLoS ONE, vol. 6, no. 1, Article ID e14634, 2011.

X. Lei, P. Yang, and D. Yao, “An empirical bayesian framework for brain-computer interfaces,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 17, no. 6, pp. 521–529, 2009.

J. Long, Y. Li, H. Wang, T. Yu, J. Pan, and F. Li, “A hybrid brain computer interface to control the direction and speed of a simulated or real wheelchair,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 20, no. 5, pp. 720–729, 2012.

Z. J. Koles, J. C. Lind, and P. Flor-Henry, “Spatial patterns in the background EEG underlying mental disease in man,” Electroencephalography and Clinical Neurophysiology, vol. 91, no. 5, pp. 319–328, 1994.

Z. J. Koles and A. C. K. Soong, “EEG source localization: implementing the spatio-temporal decomposition approach,” Electroencephalography and Clinical Neurophysiology, vol. 107, no. 5, pp. 343–352, 1998.

J. Müller-Gerking, G. Pfurtscheller, and H. Flyvbjerg, “Designing optimal spatial filters for single-trial EEG classification in a movement task,” Clinical Neurophysiology, vol. 110, no. 5, pp. 787–798, 1999.

C. Guger, H. Ramoser, and G. Pfurtscheller, “Real-time EEG analysis with subject-specific spatial patterns for a brain-computer interface (BCI),” IEEE Transactions on Rehabilitation Engineering, vol. 8, no. 4, pp. 447–456, 2000.

M. Grosse-Wentrup, C. Liefhold, K. Gramann, and M. Buss, “Beamforming in noninvasive brain-computer interfaces,” IEEE Transactions on Biomedical Engineering, vol. 56, no. 4, pp. 1209–1219, 2009.

H. Wang and W. Zheng, “Local temporal common spatial patterns for robust single-trial EEG classification,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 16, no. 2, pp. 131–139, 2008.

S. Lemm, B. Blankertz, G. Curio, and K.-R. Müller, “Spatio-spectral filters for improving the classification of single trial EEG,” IEEE Transactions on Biomedical Engineering, vol. 52, no. 9, pp. 1541–1548, 2005.

K.-R. Müller, C. W. Anderson, and G. E. Birch, “Linear and nonlinear methods for brain-computer interfaces,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 11, no. 2, pp. 165–169, 2003.

D. Yao, “Electric potential produced by a dipole in a homogeneous conducting sphere,” IEEE Transactions on Biomedical Engineering, vol. 47, no. 7, pp. 964–966, 2000.

D. Yao, “A method to standardize a reference of scalp EEG recordings to a point at infinity,” Physiological Measurement, vol. 22, no. 4, pp. 693–711, 2001.

F. Lotte and C. Guan, “Regularizing common spatial patterns to improve BCI designs: unified theory and new algorithms,” IEEE Transactions on Biomedical Engineering, vol. 58, no. 2, pp. 355–362, 2011.

Y. Koren and L. Carmel, “Robust linear dimensionality reduction,” IEEE Transactions on Visualization and Computer Graphics, vol. 10, no. 4, pp. 459–470, 2004.

G. Dornhege, B. Blankertz, G. Curio, and K.-R. Müller, “Boosting bit rates in noninvasive EEG single-trial classifications by feature combination and multiclass paradigms,” IEEE Transactions on Biomedical Engineering, vol. 51, no. 6, pp. 993–1002, 2004.

M. Naeem, C. Brunner, R. Leeb, B. Graimann, and G. Pfurtscheller, “Seperability of four-class motor imagery data using independent components analysis,” Journal of Neural Engineering, vol. 3, no. 3, article, pp. 208–216, 2006.

C.-C. Chang and C.-J. Lin, “LIBSVM: a Library for support vector machines,” ACM Transactions on Intelligent Systems and Technology, vol. 2, no. 3, article 27, 2011.

C. Vidaurre and B. Blankertz, “Towards a cure for BCI illiteracy,” Brain Topography, vol. 23, no. 2, pp. 194–198, 2010.

S. Sun and C. Zhang, “Adaptive feature extraction for EEG signal classification,” Medical and Biological Engineering and Computing, vol. 44, no. 10, pp. 931–935, 2006.

S. Sun and C. Zhang, “An optimal kernel feature extractor and its application to EEG signal classification,” Neurocomputing, vol. 69, no. 13-15, pp. 1743–1748, 2006.

P. Li, P. Xu, R. Zhang, L. Guo, and D. Yao, “L1 norm based common spatial patterns decomposition for scalp EEG BCI,” BioMedical Engineering Online, vol. 12, p. 77, 2013.

S. Sun, “The extreme energy ratio criterion for EEG feature extraction,” in Proceedings of the 18th International Conference on Artificial Neural Networks (ICANN '08), vol. 5164 of Lecture Notes in Computer Science, pp. 919–928, Springer, 2008.

S. Sun, “Extreme energy difference for feature extraction of EEG signals,” Expert Systems with Applications, vol. 37, no. 6, pp. 4350–4357, 2010.

W. Tu and S. Sun, “Semi-supervised feature extraction for EEG classification,” Pattern Analysis and Applications, vol. 16, pp. 213–222, 2013.