IPM Template V5
IPM Template V5
Mehreen Naeem a,, Muhammad Jawad Khan a, , Kashif Javed a, Saddaf Rubab a
a
National University of Sciences and Technology (NUST), Islamabad, Pakistan
Abstract
This study presents a hybrid brain computer interface (BCI) system that achieves
better accuracy based on event related potential signals. Following system based on
the P300-SSVEP hybrid sequential BCI system to decode six reactive brain
commands using ensemble classifier. For P300 we have used dataset of six images
stimulus. The feature vector is extracted from the denoised waves after filtered
through least mean square (LMS) filters. Extracted feature samples are fed into
support vector machine ensemble classifier model for classification. To minimize the
false detection rate, output from ensemble classifier trigger respective SSVEP
frequency stimulus. EEG signal is recorded at occipital region. After classification of
SSVEP signal command is send to control quadcopter. For BCI application, a virtual
quadcopter environment is created and controlled by proposed hybrid BCI system.
Results of the proposed study using LMS filter demonstrate that the sequential P300
and SSVEP model yields the combine accuracy of 96.98%.
Corresponding authors
Email address: [email protected]
Email address: [email protected]
1. Introduction
2
rate [24]. Usually, stimulus frequencies are higher than 6Hz [25]. P300 is difficult to
detect so it still required several trials of stimulus to obtain a strong signal whereas
SSVEP is detected in single trial [14].
The rest of this study is organized as follows. Section 2 investigates the related
work and presents a review on the previous methods for P300/SSVEP acquisition.
Section 3 provides details of the dataset and describes the design details of the
preprocessing process and proposed methodology. Section 4 presents the
experimental results. Section 5 outlines the conclusion and discusses on the
limitations of the proposed method.
2. Related Work
In BCI, extraction of patterns containing different order information in the signal is
a challenging task [26]. There are several modalities for the extraction of features by
removing artifacts from the corrupted signal [27, 28]. In [29] feature vector for P300
is obtained by averaging a signal and spatial/ temporal filtering. In adaptive filter
output signal is subtracted from the raw EEG signal to eliminate the artifacts [30, 27].
Hybrid strategies for features extraction were used to exploit the advantage of each
method such as Empirical mode decomposition (EMD), Blind Source Separation
(BSS), Wavelet Transform, and Regression Methods [31]. Drawback of EMD is that
it is limited to numerical simulations and lack a formal mathematical framework [32].
In another related research, P300 signal is extracted using fusion of deep features with
the temporal features based on sparse auto encoder [33]. To obtain SSVEP features,
the power spectrum of the data was taken followed by the band-pass filtered using
FIR filter [34]. The identification of EEG signal relies on a classification algorithm
[12]. Selected feature samples are fed into a relative classifier that automatically
estimating the class of data. Mostly classifiers used in BCI systems are Linear
Discriminant Analysis (LDA), Support Vector Machine (SVM), Neural Networks
(NN) [35, 36], Multi-Layer Perceptron (MLP) [26] and deep learning methods [37]. A
data partitioning technique Ensemble classifier has been introduced to reduce the
classifier variability. The average classification accuracy of hemodynamic responses
using extreme machine learning is reported as 76.7% for two classes [38]. A
3
quadcopter in a 3D physical space was driven by means of a motor imagery EEG BCI
system [4] and accuracy acquired up to 90.5% of all valid targets [39].
Several Hybrid EEG techniques have been used in the BCI field to improve
accuracy by using simple machine learning algorithms [30, 40]. In most existing BCI,
the decoding of output from EEG relies on the single classification algorithm [41 –
43]. P300 signal is a weak signal, and it has required an average number of trials and
a strong classifier with more training for detection. The problem with a single
classifier is that it may cause over-fitting over extra training. So, we have used
ensemble classifiers to reduce overfitting.
Following research presents the performance of P300 classification in comparison
with the work of Turnip et al. [44]. In the present BCI research, we have developed a
sequential Hybrid P300/SSVEP system that will generate six commands accurately
and improve classification accuracy by using an ensemble machine learning
algorithm. Previously recorded P300 dataset of six images stimulus is used.
Preprocessing and filtration of signal is done using adaptive LMS filters. After which
features are extracted and then concatenating them to feed into ensemble SVM
(eSVM) classifier. The classified output is given to the SSVEP paradigm to flash the
stimulus. A person is asked to look at the stimulus for several seconds. Hence SSVEP
signal is recorded, and the features vector is extracted by taking power spectrum of a
signal after band-pass filtration. The classified output of the SSVEP signal is used to
control quadcopter in a virtual environment.
3. Methods
3.1 Datasets
The present study uses the dataset of Hoffman [45] for the six-choice signal
paradigm. The dataset is available on the website of the EPFL BCI group. The
experiment was performed with four disabled and four healthy. The disabled
subjects had limb muscle control abilities, and all were wheelchair-bound. Some
disabled subjects had varying communication, and some are suffered from
dysarthria. The healthy subjects were Ph.D. students and none of them had any
neurological deficits. In this experiment, users were asked to focus on the screen
4
containing six images stimulus. The healthy subjects were Ph.D. students and none
of them had any neurological deficits. In this experiment, users were asked to focus
on the screen containing six images stimulus. The stimulus was images of a
television, a lamp, a window, a telephone, a door, and a radio. The images were
flashed in a random order such that one image at a time. Stimulus flashed for
100ms and in the next 300ms, none of the images was flashed. The total inter-
stimulus interval (ISI) was 400ms. The data was recorded at a 4048 Hz sampling
frequency from 32 electrodes according to the 10-20 system. The dataset contained
four sessions for each subject. Each session was consisting of six runs such that one
run for each of the six stimuli. One run consisted of an average of 135 target and
non-target trials. The entire data for one subject was consist of 3240 trials. Fig.1
[45] presents the protocol of data recording. The dataset which we are using for
SSVEP dataset recorded by Vilic [46]. The data was taken on six flickering stimuli
and recorded by placing electrodes at Oz for signal, Fz as a reference, and Fpz as
ground using the 10-20 system. A six-choice flickering stimulus panel was used for
5
a collection of a dataset. Each stimulus flash with multiple frequencies of 6Hz,
6.5Hz, 7Hz, 7.5Hz, 8.2Hz, 9.3Hz. The experiment was performed with five
subjects. The data was recorded in two sessions for each of the subjects. Each
session had ten trials of sixteen seconds. In Fig. 2 recording scenario is explained.
Fig. 2. SSVEP dataset structure for five subjects. Each consists of two sessions and
further divided in ten trial each.
6
Preprocessing of SSVEP is quite simple as compared to P300. SSVEP is detected
in the frequency domain as the range of stimulus frequencies is known. According to
given stimulus frequencies, the EEG raw signal is filtered by bandpass filters between
4Hz to 11Hz. The filtered output that is in the time domain, is converted into the
frequency domain by taking Fast Fourier Transform (FFT) as shown in Fig 3. By
analyzing the power spectrum of the FFT signal, the dominant frequency having
maximum amplitude is considered as the frequency of the stimulus on which the
subject has focused on.
The feature extractions intend to extract the important information from brain
signals. We have used Adaptive Least Means Square (LMS) filters to extract the P300
features. By using preprocessed signal, LMS changes its filter coefficient based on the
desired or reference signal by reducing the least mean square error. In our case desired
reference signal is the known noise signal. This delayed version of the noisy signal is
extracted by finding common signal between target and non-target signal of P300.
The algorithm tries to find the common noise in both the reference signal and filter
output signal, by estimating the error between the two. In Fig. 4 x(k) is the one trial
signal of 1000ms, y(k) is the filtered output, d(k) represents the reference noise signal.
The whole process is repeated for M channels and N trails.
𝑦 = 𝑤 × 𝑥 (1)
𝑒 = 𝑑 − 𝑦 (2)
𝑤 = 𝑤 + (µ × 𝑒 × 𝑥) (3)
7
Fig. 4. Least Mean Square Filter (LMS) design for P300 features extraction followed
by pre-processing of signal. Bandpass filter, down-sampling, windsorizing and
channel selection methods are the steps involved in the pre-processing of the data.
In mentioned equations, w denotes the learning weights, e is the error between the
reference signal d and the filter output y, where the µ is the learning parameter When
error gets minimized, output y which is filtered signal is subtracted from the initial
input signal to get the desired signal as shown in Fig 5. The LMS filtered output is
then analysis by statistical methods. Following statistical measurement as a feature
vector are computed: positive peak, positive peak latency, negative peak, negative
peak latency, peak to peak, peak to peak latency, and root mean square (RMS) value.
SSVEP is a periodic response generated when a subject has focused on some
repetitive stimulus. It is confined to a specific set of stimulus frequencies, so instead
of the time domain, it is analyzed in the frequency domain. The power spectrum of a
signal is taken followed by FFT. The frequency having the highest power respective
to their amplitude is considered as a feature for the classification process as shown in
Fig 6. The red dot in a figure shows the highest power of the frequency in a signal.
8
Fig. 5. Filtered output from LMS filter for target and non-target P300
Fig. 6. Power spectrum of signal for SSVEP feature extraction. Red dot shows the
highest power of the frequency in a signal.
9
3.3 Classification
An average of trials is needed for the appropriate P300 signal detection. As P300 is
difficult to detect in a single trial, so strong classifier is needed for classification. We
are using ensemble SVM for the classification of P300 of a single trial. The
architecture of model is shown in Fig. 7. A good surface decision can be made by
using ensemble SVM if the input features are well organized. We have 3240 samples
for the single subject. After shuffling and standardization, data is divided into training
and test data in the ratio of 4:1. Training data is further split into three equal divisions
by the method of Bootstrap Aggregation commonly called bagging. In this method
samples from observation are selected randomly with replacement. After the splitting
of training data into three subsets, an SVM model has developed from each
Fig. 7. Ensemble Classification model for P300 detection. Training data split into
Boostrap samples
10
concatenated features dataset. SVM model with the polynomial kernel of order two is
used for training. Kernel function is expressed as:
𝐾(𝑥𝑖 , 𝑥𝑗 ) = (1 + 𝑥𝑖′ 𝑥𝑗 )
2 (4)
The output from all the standalone models of the same types is gathered. The final
prediction is based on the aggregation of prediction of each base model. For testing
and validation, the test data is passed through each three trained models of SVM and
predicted. By the voting method, final prediction can be made as:
Where M1, M2, M3 are the trained models of SVM, and x is the test sample.
SSVEP Classification is done by SVM with kernel. It is easy to detect, so with
little training of data can give remarkable results. In AVI SSVEP dataset, we have 100
samples, 80 random samples are separate out for training. Classification of SSVEP
followed by preprocessing is shown in Fig. 8.
11
classified command decoded from sequential hybrid BCI system is transmitted with a
fixed time interval to the quadcopter. The generated output is used to control the
virtual quadcopter in a 3D environment by MATLAB simulation. Fig. 9 provides a
workflow of the proposed scheme for quadcopter control. A Graphical User Interface
(GUI) is developed, that displays the control of quadcopter. The aim of our project is
to perform maximum operations by using lesser number of outputs. The total control
commands of quadcopter are eight out of which we have formulated six commands
for the flight of quadcopter: ascend/descend movement, forward/backward movement,
clockwise and counterclockwise rotation. Following mentioned commands cover all
the possible movements of flight. The remaining right and left movements are
covered by turning the quadcopter in the desired direction and then move forward.
The display panel on the GUI shows the respective commands for all the six outputs
generated from the hybrid model shown in Fig. 10.
4. Experiment
In this proposed scheme of sequential hybrid BCI model, P300 signal
preprocessing step includes the channel selection, windowing, standardization and
LMS filter. The filtered output for eight subjects with target and non-target P300
results are shown in Fig. 11. The positive deflection of P300 signal of target stimulus
is observed nearly at 300ms in graphs for subjects 1, 2, 3 and 8. For Subjects 4, 5, 6
Fig. 9. Structure of adopted methodology. Hybrid sequential P300 and SSVEP signal
architecture used to control flight of the quadcopter
12
Fig. 10. Graphical User Interface for quadcopter flight control. It consists of a panel
having list of commands and trajectory plot
13
Fig. 11. LMS filter results for eight subjects. Target signal was shown with blue and
non-target signal with red line
14
Fig. 12. ROC plots for three SVM models
15
the movement of a virtual quadcopter. The multiple trajectories were executed in
MATLAB stimulation. The trajectories followed by the quadcopter were shown in
Fig. 15.
Fig. 13. SSVEP Pre-processing at occipital region for frequencies 6Hz, 6.5Hz and 7 Hz.
Fig. 14. SSVEP Pre-processing at occipital region for frequencies 7.5Hz, 8.2H and, 9.3
Hz.
16
Fig. 15. Quadcopter trajectories using commands generated from hybrid model.
5. Discussion
In this research, we have decoded six active BCI commands by using sequential
hybrid P300/SSVEP. Available dataset of P300 and SSVEP was used. P300 dataset of
Hoffmann was composed of six images stimulus. If target P300 is truly classified by
17
Table 1. Accuracies reported in previous related research
18
with SSVEP [42], [47] shown in Table 1. Kunda et al. had implemented sparse
autoencoder for deep features extraction with an ensemble of SVMs on the BNCI
Horizon dataset and achieved an average accuracy of 95.5% [33]. In the present study,
an ensemble SVM was used with LMS filter for P300 classification with 96.63%
accuracy. The comparison of 8 subjects with eLDA model and eSVM model are
represented in a graph Fig. 16. For P300 single classifier may cause an overfitting
problem over more training. The ensemble classifier reduces the effect of overfitting
by making a prediction from the aggregation of multiple models [48].
BCI-based SSVEP has a high information transfer rate (ITR) with a less training
period. Ojha et al. had used empirical mode decomposition (EMD) decomposition
mode on AVI SSVEP dataset with a classification accuracy of 81.45% [49]. In
another research maximum of 87.5% accuracy on the same dataset was reported [50].
By using a frequency-based approach we have achieved 94.6% classification accuracy
by SVM learner. SSVEP required training for its detection. If we combine P300 and
SSVEP less training for its detection. If we combine P300 and SSVEP modalities
there is less chance of false detection.
Classification Accuracy
100 97.2% 96.7% 97.26% 97.43% 97.84% 97.2%
95.45%
94%
95
90
85.04% 84.61% 84.82%
85 83.7% 82.72% 83.33%
81.58% 82.24%
80
75
70
S1 S2 S3 S4 S5 S6 S7 S8
eLDA eSVM
19
6. Conclusion
We have employed a sequential hybrid P300/SSVEP system that generates six
commands accurately. Six images stimulus oddball paradigm is used to generate the
P300 signal. The LMS filter used with common noise as a reference provides a noise-
free signal. After LMS filtration followed by pre-processing, features are extracted
and fed into the ensemble SVM model. Classified output triggers the SSVEP
frequency signal respective to ensemble output. After signal acquisition and filtering,
extracted features are classified by the SVM algorithm. The generated output is used
to control the virtual quadcopter on MATLAB simulation. The collective accuracy for
Hybrid Model achieved using the ensemble classifier was 96.98%, relatively
improved than other conventional classifiers.
Acknowledgements
References
[1] Z.T. Al-qaysi, B.B. Zaidan, A.A. Zaidan, M.S. Suzani, Computer Methods
and Programs in Biomedicine A review of disability EEG based wheelchair
control system : Coherent taxonomy , open challenges and recommendations,
Comput. Methods Programs Biomed. 164 (2018) 221–237.
https://doi.org/10.1016/j.cmpb.2018.06.012.
20
[4] M. Aljalal, R. Djemal, S. Ibrahim, Robot Navigation Using a Brain Computer
Interface Based on Motor Imagery, J. Med. Biol. Eng. (2018).
https://doi.org/10.1007/s40846-018-0431-9.
[5] X. Duan, S. Xie, X. Xie, Y. Meng, Z. Xu, Quadcopter Flight Control Using a
Non-invasive Multi-Modal Brain Computer Interface, 13 (2019) 1–11.
https://doi.org/10.3389/fnbot.2019.00023.
[8] R. Sharma, EEG Signal Based Movement Control for Mobile Robots, (2019).
[12] J. Zhao, W. Li, M. Li, Comparative Study of SSVEP- and P300- Based
Models for the Telepresence Control of Humanoid Robots, (2015) 1–18.
https://doi.org/10.1371/journal.pone.0142168.
21
https://doi.org/10.5772/56135.This.
[14] M.R. Meshriky, S. Eldawlatly, G.M. Aly, An Intermixed Color Paradigm for
P300 Spellers : A Comparison with Gray-scale Spellers, (2017).
https://doi.org/10.1109/CBMS.2017.123.
[15] Y. Zhang, Q. Zhao, J. Jing, X. Wang, human faces A novel BCI based on
ERP components sensitive to configural processing of human faces, (2012).
https://doi.org/10.1088/1741-2560/9/2/026018.
[16] B.H. Kim, M. Kim, S. Joevaluated, Quadcopter fl ight control using a low-
cost hybrid interface with EEG-based classi fi cation and eye tracking,
Comput. Biol. Med. (2014) 1–11.
https://doi.org/10.1016/j.compbiomed.2014.04.020.
[20] J. Li, J. Pu, H. Cui, X. Xie, S. Xu, T. Li, Y. Hu, An Online P300 Brain –
Computer Interface Based on Tactile Selective Attention of Somatosensory
Electrical Stimulation, J. Med. Biol. Eng. (2018).
https://doi.org/10.1007/s40846-018-0459-x.
22
Imagery and P300 potential, 4320 (2017) 1–10.
https://doi.org/10.1109/TNSRE.2017.2766365.
[22] A.M. Norcia, L.G. Appelbaum, J.M. Ales, B.R. Cottereau, B. Rossion, The
steady-state visual evoked potential in vision research : A review, 15 (2015)
1–46. https://doi.org/10.1167/15.6.4.doi.
[23] Z. Wang, Y. Yu, M. Xu, Y. Liu, E. Yin, Z. Zhou, Towards a Hybrid BCI
Gaming Paradigm Based on Motor Imagery and SSVEP Towards a Hybrid
BCI Gaming Paradigm Based on Motor Imagery and SSVEP, Int. J. Human–
Computer Interact. 00 (2018) 1–9.
https://doi.org/10.1080/10447318.2018.1445068.
[24] M. Xu, H. Qi, B. Wan, T. Yin, Z. Liu, A hybrid BCI speller paradigm
combining P300 potential and the SSVEP, (2013).
https://doi.org/10.1088/1741-2560/10/2/026001.
23
https://doi.org/10.1088/1741-2552/aaf046.
[29] J. Long, J. Wang, T. Yu, An Efficient Framework for EEG Analysis with
Application to Hybrid Brain Computer Interfaces Based on Motor Imagery
and P300, 2017 (2017).
[31] X. Jiang, G. Bian, Removal of Artifacts from EEG Signals : A Review, (2019)
1–18. https://doi.org/10.3390/s19050987.
[34] Y. Chen, S. Chen, I.A.E. Zaeni, C. Wu, applied sciences Fuzzy Tracking and
Control Algorithm for an SSVEP-Based BCI System, (2016) 1–11.
https://doi.org/10.3390/app6100270.
[35] S. Naseem, K. Javed, M.J. Khan, S. Rubab, M.A. Khan, Y. Nam, Integrated
CWT-CNN for epilepsy detection using multiclass EEG dataset, Comput.
Mater. Contin. 69 (2021) 471–486.
https://doi.org/10.32604/cmc.2021.018239.
[36] S. Kundu, S. Ari, Biomedical Signal Processing and Control P300 based
character recognition using convolutional neural network and support vector
machine, Biomed. Signal Process. Control. 55 (2020) 101645.
24
https://doi.org/10.1016/j.bspc.2019.101645.
[37] A. Sarwar, K. Javed, M.J. Khan, S. Rubab, O.Y. Song, U. Tariq, Enhanced
accuracy formotor imagery detection using deep learning for BCI, Comput.
Mater. Contin. 68 (2021) 3825–3840.
https://doi.org/10.32604/cmc.2021.016893.
25
[45] A. Vilic, L. Revision, The AVI SSVEP Dataset, (2013) 1–6.
[49] M.K. Ojha, M.K. Mukul, Detection of Target Frequency from SSVEP Signal
Using Empirical Mode Decomposition for SSVEP Based BCI, Wirel. Pers.
Commun. 116 (2021) 777–789. https://doi.org/10.1007/s11277-020-07738-9.
[50] A.F. Demir, H. Arslan, I. Uysal, Bio-Inspired Filter Banks for Frequency
Recognition of SSVEP-Based Brain-Computer Interfaces, IEEE Access. 7
(2019) 160295–160303. https://doi.org/10.1109/ACCESS.2019.2951327.
26