Adaptive CSP with Subspace Alignment for Subject-To-Subject Transfer in Motor Imagery Brain-Computer Interfaces

Total Page:16

File Type:pdf, Size:1020Kb

Adaptive CSP with Subspace Alignment for Subject-To-Subject Transfer in Motor Imagery Brain-Computer Interfaces IEEE Proceedings of the 6th International Winter Conference on Brain-Computer Interface Adaptive CSP with subspace alignment for subject-to-subject transfer in motor imagery brain-computer interfaces Yiming Jin Mahta Mousavi Virginia R. de Sa Computer Science and Engineering Electrical and Computer Engineering Department of Cognitive Science Shanghai Jiao Tong University UC San Diego UC San Diego Shanghai, China La Jolla, CA, USA La Jolla, CA, USA [email protected] [email protected] [email protected] Abstract—In brain-computer interfaces, adapting a classifier signals [5], [6]. Previous adaptive CSP methods rarely check from one user to another is challenging but essential to reduce the distribution of extracted features after updating the covari- training time for new users. Common Spatial Patterns (CSP) ance matrix. A subtle change in the covariance matrix may is a widely used method for learning spatial filters for user specific feature extraction but the performance is degraded when lead to significant difference between source and target data applied to a different user. This paper proposes a novel Adaptive filtered by CSP. To reduce domain variance of output data, we Selective Common Spatial Pattern (ASCSP) method to update the propose a novel Adaptive Selective Common Spatial Pattern covariance matrix using selected candidates. Subspace alignment (ASCSP) method that selects the most appropriate candidates is then applied to the extracted features before classification. The to update the covariance matrix without resulting in a great proposed method outperforms the standard CSP and adaptive CSP algorithms previously proposed. Visualization of extracted change to the distribution of extracted features. Subspace features is provided to demonstrate how subspace alignment alignment [7] is then applied to the extracted features before contributes to reduce the domain variance between source and final classification to further reduce the domain variance. target domains. Keywords—BCI, motor imagery, transfer learning, CSP, sub- II. DATASET AND METHOD space alignment A. Dataset I. INTRODUCTION To test our methodology, we used data from 6 participants from the study previously published in [8]. In this study, par- A brain-computer interface (BCI) is a platform that can ticipants were instructed to perform kinesthetic motor imagery interpret a user’s brain signals, such as those reflected in of their left or right hand to control a cursor to hit a target on electroencephalogram (EEG) signals, and send commands to a screen in front of them. In each trial, a cursor was presented the external world. Performing real or imagined movements at the center of the monitor and the target at either end - the will elicit different neural patterns which can be detected center was three cursor steps away from each end. After 1.5 by BCI. BCI-based user specific motor imagery has been seconds the target would vanish to reduce distraction and then widely investigated and common spatial patterns (CSP) has the cursor began to move at the speed of one step per second. been successfully applied to this problem [1]. Participants were led to believe that they controlled the cursor; However, due to the distortion of temporal and local infor- however, the cursor moved based on a pre-determined order mation of brain signals and variance between subjects, it is for the purposes of the original study [8]. difficult to classify EEG patterns in one user using a classifier Data were collected using a 64-channel BrainAmp system trained in another user. To overcome these drawbacks, various (Brain Products GmbH). The electrodes were distributed on modified CSP methods have been proposed to obtain more ro- the cap based on the International 10-20 system. Data were bust filters (see e.g. [2]). There are two main ways to approach collected at 5000 Hz sampling rate but were downsampled the problem: One is to regularize the objective function by to 500 Hz for further analysis. Data were analyzed offline adding a penalty term. Such regularization uses priors to guide in MATLAB and EEGLAB [9]. For further details on the the optimization process toward robust spatial filters. The other pre-processing of the data please refer to [8]. In this study, is to incorporate a weighted average of covariance matrices we use data 150 ms to 950 ms after the cursor movements from other subjects to regularize the current covariance matrix towards the target and only consider the frequency band of [3], [4]. 7-12 Hz as this covers the mu band for motor imagery. Since Inspired by adaptive learning, spatial filters learned by CSP in the pre-processing phase for each participant up to five can be adapted recursively by new incoming unlabeled EEG channels contaminated with muscle and other artifacts were Funding: NSF: SMA1041755, IIS1528214, NIH:NR013500, IBM, Adobe. removed, the common channels shared by all participants were selected, as the same channels are required for subject-to- Based on this equation, the optimal M ∗ is obtained as ∗ 0 subject transfer. M = X XT . This implies that the new coordinate system is S 0 B. Common Spatial Patterns equivalent to Xa = XSXSXT . Subspace alignment domain adaptation algorithm [7] is presented in Algorithm 1. The common spatial pattern (CSP) algorithm was suggested for motor imagery classification by Ramoser et al. [1]. CSP Algorithm 1: Subspace Alignment [7] learns optimal spatial filters that maximize filtered variance for one class and simultaneously minimize filtered variance for Input: Source data S, target data T , reduced dimension d the other class. The ith trial of the pre-processed EEG data Output: Subspace aligned data Sa;TT 1 XS PCA(S; d) for each class is denoted as an M × N matrix Xi(y) with 2 XT PCA(T; d) class label y 2 f1; 2g, where M is the number of channels 0 and N is the number of samples. The averaged normalized 3 Xa XSXSXT 4 S = SX spatial covariance matrix Cy is computed for each class as: a a 5 TT = TXT ny T 1 X Xi(y)X (y) 6 return Sa;TT C = i ; (1) y n trace(X (y)XT (y)) y i=1 i i where ny is the number of EEG trials in class y and T D. Adaptive Selective CSP represents the matrix transpose operator. The optimal set of CSP filters can be found to maximize the following Rayleigh With the help of CSP and SA, a novel adaptive selective quotient: CSP (ASCSP) is proposed. Current adaptive CSP methods T WCyW ignore the change of the distribution of the CSP filtered fea- max T : (2) tures extracted after updating the covariance matrices. ASCSP W (C1 + C2)W checks the variance of CSP filtered data before updating the This problem can be solved by the generalized eigenvalue covariance matrix and selects the probable candidate trials problem in the form: WC = ΛWC , where Λ is a 1 2 from target data that maintain the similarity of source and diagonal matrix containing eigenvalues of C−1C sorted in 2 1 target distribution. descending order and the matrix W consists of corresponding Given EEG trials Xs = fXs1 [ Xs2 g from source subject eigenvectors. With the transformation matrix W , X(y) is labeled with left (class 1) and right (class 2) respectively spatially filtered to obtain: Z = WX(y). First and last m and unlabeled trials Xt from target subject, C and C are rows of Z are selected to discriminate the two classes. Feature 1 2 initialized as averaged normalized covariance matrices from vector of rth spatial filter is constructed as: Xs1 and Xs2 using Eq. (1). var(z ) f = log r ; (3) The first step of ASCSP is to select candidate trial xi from r P2m t j=1 var(zj) X for class 1. CSP filters are trained from source labeled data and applied to both source (Xs) and target (Xt) data to extract where var() is the variance calculator and z (r=1,..,2m) is r logarithmic features. Subspace alignment is then applied to the rth row of Z. The logarithmic transformation makes the the logarithmic features to reduce domain discrepancy. Linear distribution of f more similar to Gaussian. r Discriminant Analysis (LDA) is trained upon the extracted C. Subspace Alignment features in source aligned subspace and then predicts the Fernando et al. [7] proposed a domain adaptation algorithm label in the target aligned subspace. If the prediction is class based on unsupervised subspace alignment (SA). This algo- 1 then this trial is selected as the candidate to update the rithm is part of our proposed method to reduce the variance covariance matrix later, otherwise another trial is picked and between source and target domains. Given source data yS and tested sequentially until a label of class 1 is obtained. target data yT , source subspace XS and target subspace XT After selecting the candidate xi for class 1, candidate xj for are generated from the eigenvectors of these input data by class 2 is chosen from remaining Xt. New covariance matrices ~ ~ PCA. Then the given data are projected to the subspaces by C1 and C2 are obtained from: the operations ySXS and yT XT . A linear transformation M is learned to map the source subspace to the target subspace new C1n1 cov(xi) by minimizing a Bregman matrix divergence: C1 + ; n1 + 1 trace(cov(xi))(n1 + 1) ∗ 2 M = argminM jjXSM − XT jjF ; (4) new C2n2 cov(xj) C2 + ; (6) 2 n2 + 1 trace(cov(xj))(n2 + 1) where jj · jjF is the Frobenius norm. Because the Frobenius norm is invariant to orthonormal operations, the objective where cov() calculates the covariance matrix. CSP is trained s function can be rewritten as: with C~1 and C~2 and applied to both source X and target 0 0 t ~ ∗ 2 X data.
Recommended publications
  • Invariant Common Spatial Patterns: Alleviating Nonstationarities in Brain-Computer Interfacing
    Invariant Common Spatial Patterns: Alleviating Nonstationarities in Brain-Computer Interfacing Benjamin Blankertz1,2 Motoaki Kawanabe2 Ryota Tomioka3 Friederike U. Hohlefeld4 Vadim Nikulin5 Klaus-Robert Müller1,2 1TU Berlin, Dept. of Computer Science, Machine Learning Laboratory, Berlin, Germany 2Fraunhofer FIRST (IDA), Berlin, Germany 3Dept. Mathematical Informatics, IST, The University of Tokyo, Japan 4Berlin School of Mind and Brain, Berlin, Germany 5Dept. of Neurology, Campus Benjamin Franklin, Charité University Medicine Berlin, Germany {blanker,krm}@cs.tu-berlin.de Abstract Brain-Computer Interfaces can suffer from a large variance of the subject condi- tions within and across sessions. For example vigilance fluctuations in the indi- vidual, variable task involvement, workload etc. alter the characteristics of EEG signals and thus challenge a stable BCI operation. In the present work we aim to define features based on a variant of the common spatial patterns (CSP) algorithm that are constructed invariant with respect to such nonstationarities. We enforce invariance properties by adding terms to the denominator of a Rayleigh coefficient representation of CSP such as disturbance covariance matrices from fluctuations in visual processing. In this manner physiological prior knowledge can be used to shape the classification engine for BCI. As a proof of concept we present a BCI classifier that is robust to changes in the level of parietal α-activity. In other words, the EEG decoding still works when there are lapses in vigilance. 1 Introduction Brain-Computer Interfaces (BCIs) translate the intent of a subject measured from brain signals di- rectly into control commands, e.g. for a computer application or a neuroprosthesis ([1, 2, 3, 4, 5, 6]).
    [Show full text]
  • Canonical Correlation Approach to Common Spatial Patterns
    6th Annual International IEEE EMBS Conference on Neural Engineering San Diego, California, 6 - 8 November, 2013 Canonical Correlation Approach to Common Spatial Patterns Eunho Noh1 and Virginia R. de Sa2 Abstract— Common spatial patterns (CSPs) are a way of spa- In this paper, we propose a method called canonical tially filtering EEG signals to increase the discriminability be- correlation approach to common spatial patterns (CCACSP) tween the filtered variance/power between the two classes. The which incorporates the temporal structure of the data to proposed canonical correlation approach to CSP (CCACSP) utilizes temporal information in the time series, in addition to extract discriminative and uncorrelated sources. The method exploiting the covariance structure of the different classes, to was applied to a simulated dataset as well as two BCI find filters which maximize the bandpower difference between competition datasets to compare our approach to the standard the classes. We show with simulated data, that the unsupervised CSP algorithm. canonical correlation analysis (CCA) algorithm is better able to extract the original class-discriminative sources than the CSP II. METHODS algorithm in the presence of large amounts of additive Gaussian noise (while the CSP algorithm is better at very low noise A. CSP algorithm levels) and that our CCACSP algorithm is a hybrid, yielding CSP finds spatial filters that maximize the variance/power good performance at all noise levels. Finally, experiments on data from the BCI competitions confirm the effectiveness of of spatially filtered signals under one condition while min- the CCACSP algorithm and a merged CSP/CCACSP algorithm imizing it for the other condition.
    [Show full text]
  • Major Depression Detection from EEG Signals Using Kernel Eigen-Filter-Bank Common Spatial Patterns
    sensors Article Major Depression Detection from EEG Signals Using Kernel Eigen-Filter-Bank Common Spatial Patterns Shih-Cheng Liao 1,†, Chien-Te Wu 1,2,†, Hao-Chuan Huang 3, Wei-Teng Cheng 4 and Yi-Hung Liu 3,5,* 1 Department of Psychiatry, National Taiwan University Hospital, Taipei 10051, Taiwan; [email protected] (S.-C.L.); [email protected] (C.-T.W.) 2 School of Occupational Therapy, College of Medicine, National Taiwan University, Taipei 10051, Taiwan 3 Graduate Institute of Mechatronics Engineering, National Taipei University of Technology, Taipei 10608, Taiwan; [email protected] 4 Department of Mechanical Engineering, Chung Yuan Christian University, Chungli 32023, Taiwan; [email protected] 5 Department of Mechanical Engineering, National Taipei University of Technology, Taipei 10608, Taiwan * Correspondence: [email protected]; Tel.: +886-2-2771-2171 (ext. 2066) † These authors contributed equally to this work. Academic Editor: Ioannis Kompatsiaris Received: 14 April 2017; Accepted: 10 June 2017; Published: 14 June 2017 Abstract: Major depressive disorder (MDD) has become a leading contributor to the global burden of disease; however, there are currently no reliable biological markers or physiological measurements for efficiently and effectively dissecting the heterogeneity of MDD. Here we propose a novel method based on scalp electroencephalography (EEG) signals and a robust spectral-spatial EEG feature extractor called kernel eigen-filter-bank common spatial pattern (KEFB-CSP). The KEFB-CSP first filters the multi-channel raw EEG signals into a set of frequency sub-bands covering the range from theta to gamma bands, then spatially transforms the EEG signals of each sub-band from the original sensor space to a new space where the new signals (i.e., CSPs) are optimal for the classification between MDD and healthy controls, and finally applies the kernel principal component analysis (kernel PCA) to transform the vector containing the CSPs from all frequency sub-bands to a lower-dimensional feature vector called KEFB-CSP.
    [Show full text]
  • MATHEMATICAL ENGINEERING TECHNICAL REPORTS Spectrally
    MATHEMATICAL ENGINEERING TECHNICAL REPORTS Spectrally weighted Common Spatial Pattern algorithm for single trial EEG classification Ryota TOMIOKA, Guido DORNHEGE, Guido NOLTE, Benjamin BLANKERTZ, Kazuyuki AIHARA, and Klaus-Robert MULLER¨ METR 2006–40 July 2006 DEPARTMENT OF MATHEMATICAL INFORMATICS GRADUATE SCHOOL OF INFORMATION SCIENCE AND TECHNOLOGY THE UNIVERSITY OF TOKYO BUNKYO-KU, TOKYO 113-8656, JAPAN WWW page: http://www.i.u-tokyo.ac.jp/mi/mi-e.htm The METR technical reports are published as a means to ensure timely dissemination of scholarly and technical work on a non-commercial basis. Copyright and all rights therein are maintained by the authors or by other copyright holders, notwithstanding that they have offered their works here electronically. It is understood that all persons copying this information will adhere to the terms and constraints invoked by each author’s copyright. These works may not be reposted without the explicit permission of the copyright holder. Spectrally weighted Common Spatial Pattern algorithm for single trial EEG classification Ryota TOMIOKA1;2, Guido DORNHEGE2, Guido NOLTE2, Benjamin BLANKERTZ2, Kazuyuki AIHARA1, and Klaus-Robert MULLER¨ 2;3 1 Dept. Mathematical Informatics, IST, The University of Tokyo, Japan 2 Fraunhofer FIRST.IDA, Berlin, Germany 3 Technical University Berlin, Germany [email protected] July 5, 2006 Abstract We propose a simultaneous spatio-temporal filter optimization algo- rithm for the single trial ElectroEncephaloGraphy (EEG) classification problem. The algorithm is a generalization of the Common Spatial Pat- tern (CSP) algorithm, which incorporates non-homogeneous weighting of the cross-spectrum matrices. We alternately update the spectral weighting coefficients and the spatial projection.
    [Show full text]
  • Regularizing Common Spatial Patterns to Improve BCI Designs: Theory and Algorithms Fabien Lotte, Cuntai Guan
    Regularizing Common Spatial Patterns to Improve BCI Designs: Theory and Algorithms Fabien Lotte, Cuntai Guan To cite this version: Fabien Lotte, Cuntai Guan. Regularizing Common Spatial Patterns to Improve BCI Designs: Theory and Algorithms. 2010. inria-00476820v1 HAL Id: inria-00476820 https://hal.inria.fr/inria-00476820v1 Preprint submitted on 4 May 2010 (v1), last revised 24 Sep 2010 (v4) HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. X, NO. Y, MONTH 2010 1 Regularizing Common Spatial Patterns to Improve BCI Designs: Theory and Algorithms Fabien LOTTE∗, Member, IEEE, Cuntai GUAN, Senior Member, IEEE Abstract—One of the most popular feature extraction algo- expressed with a different formulation and therefore lack a rithms for Brain-Computer Interfaces (BCI) is the Common unifying regularization framework. Moreover, they were only Spatial Patterns (CSP) algorithm. Despite its known efficiency compared to standard CSP, and typically with only 4 or 5 and widespread use, CSP is also known to be very sensitive to noise and prone to overfitting. To address this issue, some groups subjects [12][11][13], which prevents us from assessing their have recently proposed to regularize CSP.
    [Show full text]
  • Masking Covariance for Common Spatial Pattern As Feature Extraction
    Masking Covariance for Common Spatial Pattern as Feature Extraction H. Asyraf1, M. I. Shapiai1, 2, N. A. Setiawan3, W. S. N. S. Wan Musa1 1Electronic System Engineering, MJIIT. 2Centre of Artificial Intelligence and Robotics (CAIRO), UTM Kuala Lumpur, Malaysia. 3Faculty of Engineering, Universitas Gadjah Mada Yogyakarta, Indonesia. [email protected] Abstract—Brain Computer Interface (BCI) requires accurate problem lead to the study of numerous feature extractions such and reliable discrimination of EEG’s features. One of the most as Common Spatial Pattern (CSP) [2]. Laplacian filter [3] and common method used in feature extraction for BCI is a common common average reference [4] to obtain only usable feature to spatial pattern (CSP). CSP is a technique to distinguish two distinguish the EEG signal into its category. This paper will be opposite features by computing the spatial pattern of the focusing on the implementation of CSP method to extract the measured EEG channels. The existing CSP is known to be prone to the over-fitting problem. Covariance estimation is an relevant feature from EEG Signal. important process in obtaining spatial pattern using CSP. The CSP method is an algorithm commonly used to extract the empirical covariance estimation ignores the spurious information most significant discriminative information from EEG signals. between channels of EEG device. This may cause inaccurate It was first suggested for binary classification of EEG trials [5]. covariance estimation, which results to lower accuracy CSP algorithm performs covariance estimation process where performance. In this study, a masking covariance matrix is it computes the projection of most differing power or variance introduced based on the functionality of brain region.
    [Show full text]
  • Robust Common Spatial Filters with a Maxmin Approach
    1 Robust Common Spatial filters with a Maxmin Approach Motoaki Kawanabe1†∗, Wojciech Samek2y, Klaus-Robert Muller¨ 2;3∗, and Carmen Vidaurre2 1ATR Brain Information Communication Research Laboratory Group, 2-2-2 Hikaridai Seika-cho, Soraku-gun, Kyoto 619-0288, Japan. 2Department of Machine Learning, Berlin Institute of Technology (TU Berlin), Marchstr. 23, 10587 Berlin, Germany. 3Department of Brain and Cognitive Engineering, Korea University, Anam-dong, Seongbuk-gu, Seoul 136-713, Korea. y M. Kawanabe and W. Samek contributed equally to this work. Keywords: Common Spatial Patterns, Brain-Computer Interfacing, Robustness. Abstract Electroencephalographic signals are known to be non-stationary and easily affected by artifacts, therefore their analysis requires methods that can deal with noise. In this work we present a way to robustify the popular Common Spatial Patterns (CSP) algorithm under a maxmin approach. In contrast to standard CSP that maximizes the variance ratio between two conditions based on a single estimate of the class covariance matrices, we propose to robustly compute spatial filters by maximizing the minimum vari- ance ratio within a prefixed set of covariance matrices called the tolerance set. We show that this kind of maxmin optimization makes CSP robust to outliers and reduces its tendency to overfit. We also present a data driven approach to construct a tolerance set that captures the actual variability of the covariance matrices over time and show its ability to reduce non-stationarity of the extracted features and to sig- nificantly improve the classification accuracy. We test the spatial filters derived with this approach and compare them to standard CSP and to a state-of-the-art method on a real world brain-computer interface (BCI) data set in which we expect substantial fluctuations caused by environmental differences.
    [Show full text]
  • Augmented Complex Common Spatial Patterns for Classification of Noncircular Eeg from Motor Imagery Tasks 3
    IEEE TRANSACTIONS ON NEURAL SYSTEMS AND REHABILITATION ENGINEERING, VOL. 22, NO. 1, JANUARY 2014 1 Augmented Complex Common Spatial Patterns for Classification of Noncircular EEG From Motor Imagery Tasks Cheolsoo Park, Clive Cheong Took, Senior Member, IEEE, and Danilo P. Mandic, Fellow, IEEE Abstract—A novel augmented complex-valued common spatial beta/gamma (above 20 Hz) in EEG, that have been observed pattern (CSP) algorithm is introduced in order to cater for general when subjects plan and execute their hand or finger movement complex signals with noncircular probability distributions. This is [6]–[9]. Nikouline et al. [10] demonstrated that somatosensory a typical case in multichannel electroencephalogram (EEG), due to the power difference or correlation between the data channels, yet stimuli suppressed rhythms at somatosensory cortex (SI), current methods only cater for a very restrictive class of circular and Yuan et al. [9] suggested that the and rhythms are data. The proposed complex-valued CSP algorithms account for changed by the reflection of phase coherence in thalamocortical the generality of complex noncircular data, by virtue of the use of circuits. In addition, positron emission tomography (PET) or augmented complex statistics and the strong-uncorrelating trans- functional magnetic resonance imaging (fMRI) imaging studies form (SUT). Depending on the degree of power difference of com- plex signals, the analysis and simulations show that the SUT based have established that several active brain regions are observed algorithm maximizes the inter-class difference between two motor during motor imagery tasks, including supplementary motor imagery tasks. Simulations on both synthetic noncircular sources area, superior and inferior parietal lobule, dorsal and ventral lat- and motor imagery experiments using real-world EEG support the eral pre-motor cortices, pre-frontal areas, inferior frontal gyrus, approach.
    [Show full text]
  • Arxiv:1802.09046V2 [Cs.NE] 7 Mar 2021 Boosting Learning Abilities of Children Who Have Difficulty with This
    Multiclass Common Spatial Pattern for EEG based Brain Computer Interface with Adaptive Learning Classifier Hardik Meisheri1, Nagaraj Ramrao1, Suman Mitra1 1DA-IICT, Gandhinagar, India Abstract In Brain Computer Interface (BCI), data generated from Electroencephalogram (EEG) is non-stationary with low signal to noise ratio and contaminated with artifacts. Common Spatial Pattern (CSP) algorithm has been proved to be effective in BCI for extracting features in motor imagery tasks, but it is prone to overfitting. Many algorithms have been devised to regularize CSP for two class problem, however they have not been effective when applied to multiclass CSP. Outliers present in data affect extracted CSP features and reduces performance of the system. In addition to this non-stationarity present in the features extracted from the CSP present a challenge in classification. We propose a method to identify and remove artifact present in the data during pre-processing stage, this helps in calculating eigenvectors which in turn generates better CSP features. To handle the non-stationarity, Self-Regulated Interval Type-2 Neuro-Fuzzy Inference System (SRIT2NFIS) was proposed in the literature for two class EEG classification problem. This paper extends the SRIT2NFIS to multiclass using Joint Approximate Diag- onalization (JAD). The results on standard data set from BCI competition IV shows significant increase in the accuracies from the current state of the art methods for multiclass classification. 1 Introduction There is on-going research to connect brain signals with computers/machines in order to create a comfortable communication pathway for disabled people with external environment. Brain Computer Interface (BCI) aims to provide non-muscular communication between brain and computer/machine by interpreting the thoughts in the brain.
    [Show full text]
  • Common Spatial Pattern for the Classification of Imagined
    Common Spatial Pattern for the Classification of Imagined Geometric Objects Fabio R. Llorella Costa1, Gustavo Patow1 and Jose´ M. Azor´ın2 1VIRViG, Universitat de Girona, Girona, Spain 2BMI-LAB, Universidad Miguel Hernandez,´ Elche, Spain [email protected], [email protected], [email protected] Keywords: Brain-Computer Interface, Common Spatial Pattern, Support Vector Machine, Visual Imagery. Abstract: Electroencephalographic (EEG) signals contain cognitive information, which can be used by Brain-Computer Interface (BCI) systems to control devices through thought. In this work we study the possibility of detecting the visual imagination of seven different geometric objects (triangle, circle, square, pentagon, line, hexagon and parallelogram). The power spectral density in the a band were compared offline with using common spatial pattern (CSP) and the variance of each channel, obtaining as a best result the calculation of the CSP plus variance in the a band and classifying the vector of features with a support vector machine (SVM), obtaining an average result of 52% accuracy and a kappa value of 0.43 in the classification of the seven geometrical shapes, reaching up to 83% and a kappa value of 0.78 for a single user. 1 INTRODUCTION use techniques widely used in motor imagery classi- fication, such as the common spatial pattern, which Classification of electroencephalographic (EEG) sig- has offered good results (Wang et al., 2005; DaSalla nals is not a simple task as they show a great variabil- et al., 2009) and has been used to classify support vec- ity in time, that is, they are non-stationary signals (Lo tor machines, also widely used in the classification of et al., 2009).
    [Show full text]
  • Motor-Imagery Classification Using Riemannian Geometry with Median
    electronics Article Motor-Imagery Classification Using Riemannian Geometry with Median Absolute Deviation Abu Saleh Musa Miah 1 , Md Abdur Rahim 2 and Jungpil Shin 2,* 1 Bangladesh Army University of Science and Technology, Saidpur 5311, Bangladesh; [email protected] 2 School of Computer Science and Engineering, The University of Aizu, Aizuwakamatsu, Fukushima 965-8580, Japan; [email protected] * Correspondence: [email protected] Received: 1 September 2020; Accepted: 24 September 2020; Published: 27 September 2020 Abstract: Motor imagery (MI) from human brain signals can diagnose or aid specific physical activities for rehabilitation, recreation, device control, and technology assistance. It is a dynamic state in learning and practicing movement tracking when a person mentally imitates physical activity. Recently, it has been determined that a brain–computer interface (BCI) can support this kind of neurological rehabilitation or mental practice of action. In this context, MI data have been captured via non-invasive electroencephalogram (EEGs), and EEG-based BCIs are expected to become clinically and recreationally ground-breaking technology. However, determining a set of efficient and relevant features for the classification step was a challenge. In this paper, we specifically focus on feature extraction, feature selection, and classification strategies based on MI-EEG data. In an MI-based BCI domain, covariance metrics can play important roles in extracting discriminatory features from EEG datasets. To explore efficient and discriminatory features for the enhancement of MI classification, we introduced a median absolute deviation (MAD) strategy that calculates the average sample covariance matrices (SCMs) to select optimal accurate reference metrics in a tangent space mapping (TSM)-based MI-EEG.
    [Show full text]
  • Spatial Filtering for Brain Computer Interfaces
    Spatial Filtering for Brain Computer Interfaces: A Comparison between the Common Spatial Pattern and Its Variant He He and Dongrui Wu, Senior Member, IEEE Key Laboratory of the Ministry of Education for Image Processing and Intelligent Control, School of Automation, Huazhong University of Science and Technology, Wuhan, China E-mail: [email protected], [email protected] Abstract—The electroencephalogram (EEG) is the most popu- CSP finds a set of spatial filters that can achieve good lar form of input for brain computer interfaces (BCIs). However, discrimination among different classes. At the same time, it it can be easily contaminated by various artifacts and noise, can reduce the dimensionality of the EEG signals. Different e.g., eye blink, muscle activities, powerline noise, etc. Therefore, the EEG signals are often filtered both spatially and temporally approaches have been proposed in the literature to compute to increase the signal-to-noise ratio before they are fed into a them [2], [8], [12], [17]. In this paper we review some typical machine learning algorithm for recognition. This paper considers approaches, and propose a new approach by solving a slightly spatial filtering, particularly, the common spatial pattern (CSP) different objective function in the optimization. We then use filters for EEG classification. In binary classification, CSP seeks two motor imagery datasets in BCI Competition IV to compare a set of filters to maximize the variance for one class while minimizing it for the other. We first introduce the traditional these approaches. solution, and then a new solution based on a slightly different The rest of this paper is organized as follows: Section II objective function.
    [Show full text]