Next Article in Journal
Modeling and Optimization of Laser Cladding Fixation Process for Optical Fiber Sensors in Harsh Environments
Previous Article in Journal
Long-Term Polygraphic Monitoring through MEMS and Charge Transfer for Low-Power Wearable Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

DF-SSmVEP: Dual Frequency Aggregated Steady-State Motion Visual Evoked Potential Design with Bifold Canonical Correlation Analysis

1
Department of Electrical and Computer Engineering, Concordia University, 1455 De Maisonneuve Blvd. W. EV-009.187, Montreal, QC H3G 1M8, Canada
2
Concordia Institute for Information System Engineering, Concordia University, 1455 De Maisonneuve Blvd. W. EV-009.187, Montreal, QC H3G 1M8, Canada
3
Department of Electrical Engineering and Computer Science, York University, Toronto, ON M3J 1P3, Canada
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(7), 2568; https://doi.org/10.3390/s22072568
Submission received: 20 January 2022 / Revised: 19 March 2022 / Accepted: 23 March 2022 / Published: 27 March 2022
(This article belongs to the Section Biomedical Sensors)

Abstract

:
Recent advancements in Electroencephalographic (EEG) sensor technologies and signal processing algorithms have paved the way for further evolution of Brain Computer Interfaces (BCI) in several practical applications, ranging from rehabilitation systems to smart consumer technologies. When it comes to Signal Processing (SP) for BCI, there has been a surge of interest on Steady-State motion Visual Evoked Potentials (SSmVEP), where motion stimulation is used to address key issues associated with conventional light flashing/flickering. Such benefits, however, come with the price of being less accurate and having a lower Information Transfer Rate (ITR). From this perspective, this paper focuses on the design of a novel SSmVEP paradigm without using resources such as trial time, phase, and/or number of targets to enhance the ITR. The proposed design is based on the intuitively pleasing idea of integrating more than one motion within a single SSmVEP target stimuli, simultaneously. To elicit SSmVEP, we designed a novel and innovative dual frequency aggregated modulation paradigm, called the Dual Frequency Aggregated Steady-State motion Visual Evoked Potential (DF-SSmVEP), by concurrently integrating “Radial Zoom” and “Rotation” motions in a single target without increasing the trial length. Compared to conventional SSmVEPs, the proposed DF-SSmVEP framework consists of two motion modes integrated and shown simultaneously each modulated by a specific target frequency. The paper also develops a specific unsupervised classification model, referred to as the Bifold Canonical Correlation Analysis (BCCA), based on two motion frequencies per target. The corresponding covariance coefficients are used as extra features improving the classification accuracy. The proposed DF-SSmVEP is evaluated based on a real EEG dataset and the results corroborate its superiority. The proposed DF-SSmVEP outperforms its counterparts and achieved an average ITR of 30.7 ± 1.97 and an average accuracy of 92.5 ± 2.04, while the Radial Zoom and Rotation result in average ITRs of 18.35 ± 1 and 20.52 ± 2.5, and average accuracies of 68.12 ± 3.5 and 77.5 ± 3.5, respectively.

1. Introduction

Recent advancements in Electroencephalographic (EEG) sensor technologies and Signal Processing (SP) algorithms have paved the way for further evolution of Brain Computer Interface (BCI) systems [1]. The ultimate goal of a BCI system is to establish a robust communication channel with high throughput and accuracy between brain and the outer world. BCI systems have found several practical applications ranging from rehabilitation/assistive systems to diagnosis/prognoses of neurological disorders [2,3,4]. Recent technology trends show that leading technology companies are racing to develop advanced BCI systems coupled with Augmented Reality (AR) visors. It is widely expected that AR coupled with BCI would be the next era of computing. In this regard and for integration within an AR environment.
Literature Review: Generally speaking, there are two main visual BCI Paradigms, (1) Steady-State Visually Evoked Potential (SSVEP) [5,6,7,8,9,10,11,12,13,14], where light flashing (flickering) visual stimulus is used to induce evoked potentials in the EEG signals, and; (2) Steady-State motion Visual Evoked Potentials (SSmVEP) [15,16,17,18,19], where instead of using flickering, some form of graphical motion is used to evoke potentials. The former category (SSVEP) has been the main research theme due to its high achievable Information Transfer Rate (ITR), minimal requirement for user training, and excellent interactive potentials, such as high tolerance to artifacts and robust performance across users. However, flickering light causes extensive mental stress. Continuous use of SSVEPs (looking at flickering patterns for a long period of time), therefore, may cause seizure or eye fatigue.The second category (SSmVEP) is introduced to address these issues while keeping all the aforementioned benefits of the SSVEPs.
One issue that is shared by both categories is the challenge of coding more targets with available resources. Within the context of SSVEPs, the following research works were conducted to address this issue: Reference [12] introduced simultaneous phase and frequency modulations. Reference [13] used modulation in time, i.e., lengthening of the trial duration. Two different flickering-frequencies are shown consecutively for each target. Following a similar path, Reference [10] focused on using more than one frequency in a single target (but not at the same time) via Frequency-Shift Keying (FSK) modulation, also referred to as code modulation, i.e., trial time is again used for modulation purposes. Similarly, Reference [11] used code modulation with a single frequency for each target but with different phase shifts over one trial to enhance the system. Such code VEPs [10,11,14] and phase modulation [11,12] techniques are, however, very sensitive to synchronization, and as a trial time is used for modulation purposes, ITR will be compromised. When it comes to SSmVEPs, the issue of coding more targets with enhanced ITR has not yet been considered; this paper addresses this gap.
Contributions: To address the aforementioned problem, we focus on designing a novel SSmVEP paradigm without using additional resources such as trial time, phase, and/or number of targets to enhance the ITR. The proposed design is based on the intuitively pleasing idea of using more than one simultaneous motion within a single SSmVEP target stimuli. More specifically, as shown in Figure 1, to elicit SSmVEP we designed a novel and innovative dual frequency aggregated modulation paradigm, referred to as the DF-SSmVEP, by concurrently integrating “Radial Zoom” and “Rotation” motions in a single target without increasing the trial length. Figure 2i visually compares four different paradigms: Conventional SSVEP frequency modulation is shown in Figure 2(ia), where two target frequencies, “F1” and “F2”, are evoked in two different trials via flickering. Figure 2(ib) is similar to Figure 2(ia) where now two target frequencies are used together, one after another by increasing the trial time. Figure 2(ic) illustrates two SSmVEP modulations similar to Figure 2(ia), but target frequencies are evoked now via motion of the circle. Figure 2(id) shows the proposed DF-SSmVEP design where now two target frequencies are used together simultaneously eliminating the need to increase (sacrifice) the trial length for achieving higher accuracy as is the case in code/frequency modulated SSVEPs [10,11].
The paper also develops a specific unsupervised classification model adopted to the proposed innovative DF-SSmVEP paradigm. More specifically, in contrast to existing works, we propose an unsupervised SSmVEP detection technique, called the Bifold Canonical Correlation Analysis (BCCA) using unique characteristics of the proposed dual aggregated frequency design. The BCCA exploits the availability of two motion frequencies for each target and separately considers each single frequency of the targets as a reference. The corresponding covariance coefficients are then used as extra features advancing the classification accuracy. The proposed DF-SSmVEP is evaluated based on a real EEG dataset.

2. The Proposed DF-SSmVEP

In this section, we present the proposed DF-SSmVEP framework. The designed DF-SSmVEP stimulation paradigm includes a green and black circle with two motion modes. The first motion is the “Radial Zoom” in which the size of the circle changes periodically. The second mode is the “Reciprocal Rotation” of the circle between −45 and 75 . Radial zoom motion and rotation motion are selected as candidates for integration following previous evaluations [19,20]. Generally speaking, the motion inversion frequency is defined as the frequency with which the direction of the current motion is changing within the reciprocal motion. Furthermore, the SSmVEP frequency is equal to the stimulation’s frequency, which is related to the aforementioned motion inversion frequency. When it comes to the design of the DF-SSmVEP paradigm, we had to decide on the motion choices. In this regard, we followed the intuition provided by Reference [20], i.e., to create a stimuli for evoking SSmVEPs, we can use any paradigm that has periodically changing motion. The two designs are integrated such that the focal point of one paradigm is overlaid with that of the second one. In the proposed DF-SSmVEP paradigm, the focal point will be the center of the black-green circle, i.e., the center of oscillation for the two segments of the design (resizing and oscillation of circle). To make the proposed design as efficient as possible and to reduce fatigue [21], the proposed DF-SSmVEP design does not include high contrast colors, improving the practical applicability of the SSmVEP stimuli. As shown in Figure 2(ib), the Luminance Contrast Ratio (LCR) associated with our green-black paradigm is lower than that of the conventional black-white design.
In the design stage and for selecting the underlying targeted frequencies, one needs to take into account the refresh rate of the monitor over which the SSmVEP stimuli is shown. Essentially, the monitor’s refresh rate is a restricting factor limiting the number of frames that can be shown within each motion direction change. According to Reference [21] and for achieving flexibility in the implementation of target frequencies, the underlying paradigms should have a different (variable) number of frames within each cycle. In the implemented setting, the refresh rate of the monitor was 60 Hz, based on which initially we designed the number of frames within each half of a cycle in an asymmetric fashion. We refer to this initial step as binary stimulus sequencing. More specifically, within each half of a cycle of the DF-SSmVEP, the number of frames is constructed as follows
S ( f ( T ) , i ) = square 2 π f ( T ) ( i R r ) ,
where the target frequency is denoted by f; Refresh rate of the monitor is denoted by R r ; i is used as the frame index. Finally, the stimulus sequences corresponding to f are represented by S ( f , i ) . It is worth noting that within each half of a cycle, an expansion or contraction is performed, which is equivalent to half of the reciprocal motion. Therefore, we can generate stimulus with oscillations up to f R r / k around the SSmVEP’s target frequency. More specifically, we need a minimum number of k frames over half of a cycle to have an understandable and comfortable SmVEP paradigm.
Coding Algorithm: Assume that the maximum number of targets (objects shown on the screen simultaneously) is denoted by N ( T ) . In other words, N ( T ) target frequencies are selected within the limited frequency spectrum of [ f min , f max ] available for constructing the stimuli. f i , for ( 1 i N ( T ) ), represents the target frequency for the i th target/object assumed to be sorted in an ascending order, i.e., ( f min f 1 < f 2 < < f N ( T ) f max ). These N ( T ) target frequencies need to be derived in an intelligent fashion such that the best performance among all the susceptible frequencies is achieved (i.e., achieve accuracy improvements without reducing the ITR). In the proposed dual aggregated design, each target includes two motions with two distinct frequencies. These frequencies are assigned to targets in which no two pairs of targets have more than one adjacent frequencies. More specifically, the objective of the coding algorithm is to find these two underlying target frequencies in such a way that each pair of objects at most have one adjacent target frequency. For each consecutive frequencies f i and f i + 1 , g i is defined as
g i = f i + f i + 1 2 , i [ 1 , N 1 ] f i + M , i = N
where M = min [ f i + 1 f i 2 ] i [ 1 , N 1 ] . Each g i is adjacent to f i and f i + 1 . More specifically, consider P = { ( a i , b i ) } , for ( 1 i N ), representing the set of N 5 target pairs where a i and b i are the new SSmVEP frequencies used for the ith object.s a i and b i in P are defined as follows:
a i = f 1 , i = 1 f i + 1 , i [ 2 , N 1 ] f 2 , i = N   b i = g N 1 , i = 1 g i 1 , i [ 2 , N 1 ] g N , i = N .
Pre-processing: The proposed SSmVEP paradigm is implemented via a BCI system for real EEG data collection. In this regard, the first step is pre-processing of EEG signals associated with the proposed SSmVEP paradigm, as collected EEG signals are, typically, contaminated by high/low noises and/or artifacts. In this context, one needs to apply time-domain and/or spatial filters on the EEG signals to remove the artificats/noises and allow efficient extraction of SSmVEP related features. To filter high-frequency noises and artifacts, we applied a zero-phase Chebyshev Type I bandpass filter, which passes frequencies within the rage of 2–40 Hz.

Proposed BCCA Paradigm

To learn potential linear relations that exist between two multi-dimensional EEG signals, typically, the Canonical Correlation Analysis (CCA) is used. The CCA is a statistical approach that aims at extracting two linear projection vectors in such a way that the resulting linear combination of the two original EEG signals has maximum correlation coefficient. More specifically, consider that the two original EEG signals are represented by matrices X and Y. Then, the objective is to compute projection vectors w x and w y that results in a maximized correlation coefficient, denoted by ρ between w x T X and w y T Y . This process can be performed as follows:
ρ = max E ( w x T X Y T w y ) E ( w x T X Y T w x ) E ( w y T X Y T w y ) .
Reference signals are constructed at the stimulation frequency f i as
y i = [ cos 2 π f i t , sin 2 π f i t , , cos 2 π N h f i t , sin 2 π N h f i t ] T ,
where t = 1 f s , , m f s with m denoting the sample point; The sampling rate is represented by f s , and; The number of harmonics is represented by N h . It is worth noting that the value of N h depends on the designed paradigm and varies from one scenario to another. Therefore, we used the Welch Power Spectrum of the underlying signals and computed N h experimentally. Consider X as the matrix of the EEG signals collected from K different channels. The CCA finds linear combination of coefficients with the largest correlation between X and Y . We would like to mention that use of CCA for comparing SSVEP and SSmVEP paradigms, which evoke specific frequency range and harmonics within EEG signals, is common in the literature [19,22]. These targeted frequencies are dominant in the signal harmonics and other components can be approximately modeled by linear noise. In this model, the collected signals from the involved regions, Occipital lobe, contain both evoked SSVEP harmonics and other artifacts. Moreover, the goal of CCA is increasing the Signal to Noise Ratio (SNR) and removing background EEG activities to detect evoked Signals. The standard CCA method performs canonical correlation analysis between multi-channel EEG signals and predefined sinusoidal reference signals at stimulation frequencies and then identifies the target frequency based on the canonical correlation values.
In the BCCA fusion, there is a feature vector for each sample concerning each target frequency. Contemplating the Power Spectral Density (PSD) of the EEG signal, collected during the aggregated paradigm shows that only one of the peaks is significant for some trials, i.e., one of the two modulated frequencies has more impact on the visual pathways of the brain. To capitalize on this unique property and enhance the DF-SSmVEP, the following three references are incorporated to create the feature vector
y 1 = [ cos ( 2 π f i , 1 t ) , sin ( 2 π f i , 1 t ) , , cos ( 2 π N h f i , 1 t ) , sin ( 2 π N h f i , 1 t ) ] T
y 2 = [ cos ( 2 π f i , 2 t ) , sin ( 2 π f i , 2 t ) , , cos ( 2 π N h f i , 2 t ) , sin ( 2 π N h f i , 2 t ) ] T
y c = C ( y 1 , y 2 , [ cos ( 2 π ( f i , 1 + f i , 2 ) t ) , sin ( 2 π ( f i , 1 + f i , 2 ) t ) ] T , 2 )
where f i , j represents the jth stimulation frequency of ith target, and the operator C ( a , b , c , 2 ) concatenates three matrixes a, b and c vertically. The projection of each vector is separately calculated leading to three different weight vectors between test signal X and: (i) ( w y 1 , w X 1 ) sine/cosine reference of first frequencies; (ii) ( w y 2 , w X 2 ) reference of second frequency of targets, and; (iii) ( w y c , w X c ) sine/cosine reference of both frequencies of targets. The feature vector v is
v = ρ 1 , ρ 2 , ρ c T , and ρ a = ρ 1 + ρ 2 + ρ c 3 ,
where ρ a is used as the final value to represent the correlation between the unknown sample and the frequencies of a target. It is worth mentioning that the proposed BCCA is an unsupervised technique as such there is no need for a separate training step. Consequently, all the available trials of sessions are used in the testing stage.

3. Experimental Results

In this section, we present different experimental evaluations based on a real EEG dataset collected for this study. We aim to evaluate the performance and effectiveness of the proposed DFSSmVEP paradigm and also compare it with its conventional SSmVEP counterparts, i.e., Radial Zoom and Rotation, when applied independently. We would like to mention that another possibility would be to consider comparisons with conventional SSVEP (flickering) paradigms. However, SSVEPs are inherently categorized as a different group of VEPs because of the induced level of mental load and fatigue. As such and following the common practice, we have not included results for conventional SSVEPs.

3.1. Experimental Setup

To evaluate the performance of the proposed DF-SSmVEP framework, we created a real EEG dataset from ten healthy individuals (five men and five women). Among the ten participants (each between 20 and 27 years old), none had any evidence of color or visual-recognition ailment issues, but five had prior experience with use of EEG signals and have participated previously in BCI experiments. We would like to mention that, for data collection, we followed the common, standard, and accepted approach in the BCI domain (e.g., [20]) where typically 9–11 subjects are used, and the trial length ranges from 2 s to 6 s. The data were collected with the policy certification number 3007997 of Ethical acceptability for research involving human subjects approved by Concordia University. We used g.Nautilus, which is a wireless and portable bio-signal acquisition headset from g.tech Medical Engineering to collect the required EEG signals from the ten participants. The used g.Nautilus headset has sampling rate of 500 Hz and consists of 32 active, bipolar, and wet electrodes. Following the conventional approach, we placed the reference electrode at the earlobe, and the ground electrode at the frontal position (Fpz). To collect the EEG signals, we focused on the parietal-occipital area of the brain and used signals from the following electrodes: P o 3 , P o 4 , P o 7 , P o 8 , P z , and O z . For illustrating the stimuli, we used a 21.5 -inch LED screen with resolution of 1920 × 1080 , and refresh rate of 60 Hz. The stimuli is shown over a white background while each participant was sitting at a distance of 70 cm from the screen.
Experimental Protocol: We used separate runs to test performance of the proposed DF-SSmVEP against the following two different paradigms: (i) Radial zoom, and; (ii) Rotation. Please refer to Reference [23] for videos of the test settings. Each video consists of five targets oscillating with different frequencies. For individual videos, target frequencies were 5, 6, 7, 8, and 9 Hz. For aggregated videos, target frequencies of the radial zoom pattern were 5, 6, 7, 8, and 9 Hz, and target frequencies of the rotation pattern were 5.5, 6.5, 7.5, 8.5, and 9.5 Hz. In each run, we included two consecutive sessions, in each of which the participants are asked to look directly at the target object identified with a pointer. During each session, we point to each of the underlying targets over four trials, where each of these four trials lasts for 3.5 s. A 2.5 s break is considered between each run of the consecutive trials. For evaluation purposes, ITR is used, which assesses the speed of a BCI systems as
ITR = 60 T log 2 K + σ log 2 σ + ( 1 σ ) log 2 ( 1 σ K 1 ) ,
where σ denotes the recognition accuracy; K denotes the number of stimuli, and; T represents a value computed by adding the time of the resting state and the trial time. Evaluation via the One-way Analysis Of Variance (ANOVA) [24] with Tukey post hoc analysis is also used to confirm that responses to the proposed DF-SSmVEP-stimuli is statistically meaningful (p < 0.05). Some of the functions of [25] are used in the pipeline code. We would like to mention that one-way ANOVA (which can be considered to be the parametric equivalent of the Kruskal–Wallis test) is the common approach used in the literature to measure the significance of the VEP stimuli as such and for consistency is used in our study. Alternative non-parametric tests could be used for further analysis but heteroscedasticity of data (i.e., different groups have different standard deviations) should be carefully taken into consideration.

3.2. Results

As shown in Figure 3, accuracy and ITR are measured for different time windows for each trial ranging from 0.5 to 4 s with an interval of 500 ms. The highest rate of transmission (ITR) belongs to aggregated motion (30.7 ± 1.97), which also achieved the best accuracy (92.5 ± 2.04). We used the ANOVA test following the existing literature [18,20] that used ANOVA test for significance comparison between accuracies and/or ITRs of different paradigms. It is worth noting that based on the Central Limit Theorem, we can safely assume that the samples have a normal distribution. The one-way ANOVA on accuracies of DF-SSmVEP paradigm reveals that there is no significant effect of frequencies (classes) on accuracies (F = 2.78, p = 0.065), so all target frequencies of DF-SSmVEP are feasible in BCI systems. The Tukey post hoc test on the accuracy and the ITR shows significant differences between the performance of DF-SSmVEP with BCCA and other paradigms with corresponding classifiers for which the highest accuracies are acquired, i.e., (Accuracy: p D F R = 0.042, p D F R Z = 0.002, p R R Z = 0.218; ITR: p D F R = 0.011, p D F R Z = 0.001, p R R Z = 0.264).
Table 1 shows the overall recognition accuracies and ITRs. It can be observed that aggregating two SSmVEP paradigms using the proposed DF-SSmVEP results in compelling performance improvement. The results also show superiority of the BCCA as the best unsupervised target detector among its counterparts. Figure 4 illustrates Welch PSD of two targets of aggregated motion after spatial filtering for Subject 1. Two significant peaks of SSmVEP frequency around the two frequencies of each target is observable. Generally speaking, the Welch PSD is used to estimate random signals by dividing the data with a length of N into M segments of length L. The formula to compute the used window averaged period is given by P ( w ) = 1 / M i = 1 M [ 1 L P 0 | n = 1 L w ( n ) x i ( n ) e j w n | 2 ) ] , where p 0 refers to the power of window w ( n ) , which is given by P 0 = 1 / L n = 1 L | w ( n ) | 2 . According to [26,27], using welch PSD for detecting SSVEP and SSMVEP harmonics within EEG signal is more effective compared to Short-Time Fourier Transform (STFT) and filter banks. The Welch method reduces noise in the estimated power spectra in exchange for reducing the frequency resolution. Due to the presence of noise caused by EEG artifacts, the noise reduction from Welch’s method is indeed desired. Although the length of a trial is short, the target frequency and visual frequency band is limited to 20 Hz, which means a complete cycle of oscillation takes at most 50 ms that is much less than the trial length. Hence, the band-power of signal contains signals’ power, therefore, the superposition of sub powers will intensify the signal power and reduce noise power.
Mean and standard deviation based on four different Performance Indices (PI), i.e., precision, sensitivity, specificity, and accuracy across all runs for each class, are also shown in Table 2. In these calculations, the time window is set to 3.5 s. Each row of Table 2 corresponds to the dual frequency of the proposed DF-SSmVEP paradigm in each class. To investigate the robustness of the proposed methodology, we follow the common approach, which is to show paradigms repeatedly to different subjects (in contrast to using EEG signals of one trial several times) and average the performance across Runs and Subjects. As mentioned in the experimental protocol, in each run, we included two consecutive sessions, in each of which the participants are asked to look directly at the target object identified with a pointer. During each session, we point to each of the underlying targets over four trials. Hence, for each subject, the dataset includes 40 multidimensional signals associated with five classes. Each class represents target frequencies used in the corresponding motion paradigm where the cursor points to during data session. For example, if a subject stares at the rotation paradigm oscillating with 5 Hz and the cursor points to it, the EEG signal collected in this period of time will be labeled by 5 Hz. These labels are considered to be the ground truth for the dataset. BCCA and CCA also predict target frequencies used within a modulation for a given EEG signals. In this prediction task, the ratio of true positives (EEG signals predicted by correct target frequencies) over total records represents the accuracy. More specifically, in Table 2, specificity, sensitivity, accuracy, and precision of prediction for each class are measured as follows:
Sensitivity = True Positive True Positive + False Negative
Specificity = True Negative True Negative + False Positive
Precision = True Positive True Positive + False Positive
Precision = True Positive + True Negative Total number of predictions
where “True Positive” for a given class is the number of EEG signals labeled with name of that class correctly; “True Negative” for a given class is the number of EEG signals labeled with the name of other classes correctly; “False Positive” for a given class is the number of EEG signals labeled with the name of that class incorrectly, and; “False Negative” for a given class is the number of EEG signals labeled with the name of the other classes incorrectly.

4. Conclusions

To address lower accuracy and ITR of SSmVEP designs, this paper proposed an intuitively pleasing, novel, and innovative dual frequency aggregated modulation paradigm. Referred to as the DF-SSmVEP, the novel design is constructed by concurrently integrating “Radial Zoom” and “Rotation” motions in a single target without increasing the trial length. The paper also develops a specific unsupervised classification model, referred to as the BCCA, which uses availability of two motion frequencies per each target. While conventional Radial Zoom [19] and Rotation paradigms [19] have been used previously to design SSmVEPs, the proposed DFSSmVEP is developed by simultaneous combination (use) of these two conventional paradigms jointly. In other words, the DFSSmVEP is the first SSmVEP paradigm of its nature to concurrently use two different motions in one paradigm. Additionally, the frequency modulation and allocation of harmonics within accessible frequency bandwidth (visual bandwidth) is another novel characteristic of the DFSSmVEP. This is a unique way of associating different harmonics to different motion modes in one motion paradigm, which can be considered to be a step forward toward combining more SSmVEP patterns to develop new SSmVEP-based stimuli for evoking harmonic signals more strongly. The proposed DF-SSmVEP is evaluated via a real EEG dataset achieving average ITR of 30.07 ± 1.97 and average accuracy of 92.5 ± 2.04. As a final note, we would like to mention that it is well-known that EEG signals have a highly non-linear nature, the nature of non-linearity may even vary from subject to subject or even session to session. Although CCA has been widely used for comparing SSmVEP paradigms, its non-linear counterparts could be considered alternatively to relax the assumption of linear noise for removing artifacts to potentially improve the overall performance of the paradigm. In particular, advanced Deep Neural Networks (DNN) (e.g., similar in nature to the ones developed in our prior work [28]) can be designed to learn the type of non-linearity for each subject and remove it accordingly. However, such an approach, typically, is supervised needing a training phase making unsupervised CCA approach, although less accurate, more appealing. Development of semi-supervised or self-supervised DNN-based methods seem to be an interesting venue to pursue as a future direction.

Author Contributions

Conceptualization, R.K. and A.M.; methodology, R.K. and A.M.; software, R.K.; validation, R..K., A.M., H.B. and A.A.; writing—original draft preparation, R..K. and A.M.; writing—review and editing, R.K., A.M. and H.B.; supervision, A.M. and A.A.; funding acquisition, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the Natural Sciences and Engineering Research Council (NSERC) of Canada through the NSERC Discovery Grant RGPIN-2016-04988.

Institutional Review Board Statement

The data were collected with the policy certification number 3007997 of Ethical acceptability for research involving human subjects approved by Concordia University.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The anonymized dataset created for this research work is available through the following page: https://github.com/raykakarimi/DFSSMVEP (accessed on 19 January 2022).

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Chen, X.; Xu, X.; Liu, A.; Lee, S.; Chen, X.; Zhang, X.; McKeown, M.J.; Wang, Z.J. Removal of Muscle Artifacts From the EEG: A Review and Recommendations. IEEE Sens. J. 2019, 19, 5353–5368. [Google Scholar]
  2. Shahtalebi, S.; Mohammadi, A. Bayesian Optimized Spectral Filters Coupled with Ternary ECOC for Single-Trial EEG Classification. IEEE Trans. Neural Syst. Rehabil. Eng. 2018, 26, 2249–2259. [Google Scholar]
  3. Samanta, K.; Chatterjee, S.; Bose, R. Cross-Subject Motor Imagery Tasks EEG Signal Classification Employing Multiplex Weighted Visibility Graph and Deep Feature Extraction. IEEE Sens. Lett. 2020, 4, 1–4. [Google Scholar]
  4. Dagois, E.; Khalaf, A.; Sejdic, E.; Akcakaya, M. Transfer Learning for a Multimodal Hybrid EEG-fTCD Brain–Computer Interface. IEEE Sens. Lett. 2019, 3, 1–4. [Google Scholar]
  5. Zhang, Y.; Yin, E.; Li, F.; Zhang, Y.; Guo, D.; Yao, D.; Xu, P. Hierarchical Feature Fusion Framework for Frequency Recognition in SSVEP-based BCIs. Neural Netw. 2019, 119, 1–9. [Google Scholar]
  6. Kubacki, A. Use of Force Feedback Device in a Hybrid Brain-Computer Interface Based on SSVEP, EOG and Eye Tracking for Sorting Items. Sensors 2021, 21, 7244. [Google Scholar]
  7. Ikeda, A.; Washizawa, Y. Steady-State Visual Evoked Potential Classification Using Complex Valued Convolutional Neural Networks. Sensors 2021, 21, 5309. [Google Scholar]
  8. De la Cruz-Guevara, D.R.; Alfonso-Morales, W. Caicedo-Bravo, E. Solving the SSVEP Paradigm Using the Nonlinear Canonical Correlation Analysis Approach. Sensors 2021, 21, 5308. [Google Scholar]
  9. Chen, Y.-J.; Chen, P.-C.; Chen, S.-C.; Wu, C.-M. Denoising Autoencoder-Based Feature Extraction to Robust SSVEP-Based BCIs. Sensors 2021, 21, 5019. [Google Scholar]
  10. Zhao, X.; Zhao, D.; Wang, X.; Hou, X. A SSVEP Stimuli Encoding Method Using Trinary Frequency-shift Keying Encoded SSVEP (TFSK-SSVEP). Front. Hum. Neurosci. 2017, 11, 278. [Google Scholar]
  11. Wei, Q.; Feng, S.; Lu, Z. Stimulus Specificity of Brain-Computer Interfaces based on Code Modulation Visual Evoked Potentials. PLoS ONE 2016, 11, e0156416. [Google Scholar]
  12. Nakanishi, M.; Wang, Y.; Wang, Y.T.; Mitsukura, Y.; Jung, T.-P. A High-Speed Brain Speller using Steady-State Visual Evoked Potentials. Int. J. Neural Syst. 2014, 24, 1450019. [Google Scholar]
  13. Zhang, Y.; Xu, P.; Liu, T.; Hu, J.; Zhang, R.; Yao, D. Multiple Frequencies Sequential Coding for SSVEP-based Brain-Computer Interface. PLoS ONE 2012, 7, e29519. [Google Scholar]
  14. Kadioglu, B.; Yildiz, I.; Closas, P.; Fried-Oken, M.B.; Erdogmus, D. Robust Fusion of c-VEP and Gaze. IEEE Sens. Lett. 2019, 3, 1–4. [Google Scholar]
  15. Zhang, X.; Hou, W.; Wu, W.; Chen, L.; Jiang, N. Enhancing Detection of SSMVEP Induced by Action Observation Stimuli Based on Task-Related Component Analysis. Sensors 2021, 21, 5269. [Google Scholar]
  16. Beveridge, R.; Wilson, S.; Callaghan, M.; Coyle, D. Neurogaming with Motion-onset Visual Evoked Potentials (MVEPs): Adults versus Teenagers. IEEE Trans. Neural Syst. Rehabil. Eng. 2019, 27, 572–581. [Google Scholar]
  17. Yan, W.; Xu, G.; Chen, L.; Zheng, X. Steady-state Motion Visual Evoked Potential (SSmVEP) Enhancement Method based on Time-Frequency Image Fusion. Comput. Intell. Neurosci. 2019, 2019, 9439407. [Google Scholar]
  18. Han, C.; Xu, G.; Xie, J.; Chen, C.; Zhang, S. Highly Interactive Brain-Computer Interface based on Flicker-Free Steady-State Motion Visual Evoked Potential. Sci. Rep. 2018, 8, 5835. [Google Scholar]
  19. Chai, X.; Zhang, Z.; Guan, K.; Liu, G.; Niu, H. A Radial Zoom Motion-based Paradigm for Steady State Motion Visual Evoked Potentials. Front. Hum. Neurosci. 2019, 13, 127. [Google Scholar]
  20. Yan, W.; Xu, G.; Xie, J.; Li, M.; Dan, Z. Four Novel Motion Paradigms based on Steady-state Motion Visual Evoked Potential. IEEE Trans. Biomed. Eng. 2018, 65, 1696–1704. [Google Scholar]
  21. Nakanishi, M.; Wang, Y.; Mitsukura, Y.; Jung, T.P. An Approximation Approach for Rendering Visual Flickers in SSVEP-based BCI using Monitor Refresh Rate. In Proceedings of the 2013 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Osaka, Japan, 3–7 July 2013; pp. 2176–2179. [Google Scholar]
  22. Wang, Y.; Jung, T. Visual stimulus design for high-rate SSVEP BCI. Electron. Lett. 2010, 46, 1057–1058. [Google Scholar]
  23. Available online: https://github.com/raykakarimi/DF-SSmVEP-videos (accessed on 20 December 2021).
  24. Xie, J.; Xu, G.; Wang, J.; Zhang, F.; Zhang, Y. Steady-state Motion Visual Evoked Potentials Produced by Oscillating Newton’s Rings: Implications for Brain-computer Interfaces. PLoS ONE 2012, 7, e39707. [Google Scholar]
  25. Zhang, Y.; Zhou, G.; Jin, J.; Wang, X.; Cichocki, A. Frequency recognition in SSVEP-based BCI using multiset canonical correlation analysis. Int. J. Neural Syst. 2014, 24, 1450013. [Google Scholar]
  26. Carvalho, S.N.; Costa, T.B.S.; Uribe, L.F.S.; Soriano, D.C.; Yared, G.F.G.; Coradine, L.C.; Attux, R. Comparative analysis of strategies for feature extraction and classification in SSVEP BCIs. Biomed. Signal Process. Control 2015, 21, 34–42. [Google Scholar]
  27. Yan, W.; Xu, G.; Li, M.; Xie, J.; Han, C.; Zhang, S.; Luo, A.; Chen, C. Steady-state motion visual evoked potential (SSMVEP) based on equal luminance colored enhancement. PLoS ONE 2017, 12, e0169642. [Google Scholar]
  28. Karimi, R.; Mohammadi, A.; Rosero, L.; Asif, A. Deep Video Canonical Correlation Analysis for Steady State motion Visual Evoked Potential Feature Extraction. In Proceedings of the 2020 28th European Signal Processing Conference (EUSIPCO), Amsterdam, The Netherlands, 18–21 January 2021; pp. 1130–1134. [Google Scholar]
Figure 1. Proposed DF-SSmVEP paradigm developed by concurrent inclusion of two types of the motion (rotation and resizing).
Figure 1. Proposed DF-SSmVEP paradigm developed by concurrent inclusion of two types of the motion (rotation and resizing).
Sensors 22 02568 g001
Figure 2. (i) Comparison between existing SSVEP frequency modulation schemes (a,b) with SSmVEP (c), and the proposed DF-SSmVEP (d). (ii) Luminance contrast ratio of the colors used in the designed DF-SSmVEP.
Figure 2. (i) Comparison between existing SSVEP frequency modulation schemes (a,b) with SSmVEP (c), and the proposed DF-SSmVEP (d). (ii) Luminance contrast ratio of the colors used in the designed DF-SSmVEP.
Sensors 22 02568 g002
Figure 3. Mean and standard deviation across all subjects for each time window: (a) Accuracy comparisons. (b) ITR comparisons.
Figure 3. Mean and standard deviation across all subjects for each time window: (a) Accuracy comparisons. (b) ITR comparisons.
Sensors 22 02568 g003
Figure 4. (a) The PSD plot based on aggregated SSmVEP with 6–9.5 Hz target frequencies. (b) Similar to (a) but with 8–6.5 Hz target frequencies.
Figure 4. (a) The PSD plot based on aggregated SSmVEP with 6–9.5 Hz target frequencies. (b) Similar to (a) but with 8–6.5 Hz target frequencies.
Sensors 22 02568 g004
Table 1. Mean accuracy (%) and mean ITR (bits/min) comparison between four methods of spatial filtering across the three motions. The best ITR among different time windows are reported for each filter. Maximum Contrast Fusion (MCF) [20] and T-F Image Fusion [17] are two types of spatial filtering.
Table 1. Mean accuracy (%) and mean ITR (bits/min) comparison between four methods of spatial filtering across the three motions. The best ITR among different time windows are reported for each filter. Maximum Contrast Fusion (MCF) [20] and T-F Image Fusion [17] are two types of spatial filtering.
ParadigmsRadial ZoomRotationDF-SSmVEP
Filters
ACCITRACCITRACCITR
MCF + CCA68.1218.3577.520.5281.8821.89
T-F Image Fusion + CCA59.313.3968.7513.7363.2513.93
CCA Fusion63.517.2476.1718.0584.3823.73
 BCCA Fusion----92.530.7
Table 2. Comparison between mean and standard deviation values associated with performance indices (precision, sensitivity, specificity, and accuracy across) all run for three paradigms, i.e., Rotation (R), Radial Zoom (RZ), and the proposed DF-SSmVEP (denoted by DF).
Table 2. Comparison between mean and standard deviation values associated with performance indices (precision, sensitivity, specificity, and accuracy across) all run for three paradigms, i.e., Rotation (R), Radial Zoom (RZ), and the proposed DF-SSmVEP (denoted by DF).
PISpecificitySensitivityPrecisionAccuracy
Classes
RRZDFRRZDFRRZDFRRZDF
9 Hz or0.9120.9340.9930.9000.6000.8750.7350.7100.9770.9100.8670.970
(9 Hz, 7.5 Hz)±0.052±0.037±0.013±0.098±0.226±0.083±0.137±0.061±0.047±0.054±0.020±0.01
6 Hz or0.9400.8930.9840.7250.8750.9620.7570.7020.9470.8970.8900.980
(6 Hz, 9.5 Hz)±0.023±0.062±0.027±0.098±0±0.060±0.091±0.167±0.870±0.027±0.050±0.023
5 Hz or0.9460.8750.9400.6250.9370.9750.7370.6630.8290.8820.88750.952
(5 Hz, 8.5 Hz)±0.025±0.044±0.033±0.220±0.106±0.053±0.118±0.068±0.100±0.054±0.017±0.029
7 Hz or0.9810.9930.9930.7870.31210.9370.9330.9770.9420.8570.995
(7 Hz, 5.5 Hz)±0.030±0.013±0.013±0.177±0.135±0±0.100±0.140±0.047±0.026±0.031±0.011
8 Hz or0.9400.9060.9870.8500.6870.8120.8090.6490.9390.9220.8620.952
(8 Hz, 6.5 Hz)±0.047±0.062±0.016±0.098±0.244±0.135±0.131±0.221±0.079±0.027±0.095±0.036
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Karimi, R.; Mohammadi, A.; Asif, A.; Benali, H. DF-SSmVEP: Dual Frequency Aggregated Steady-State Motion Visual Evoked Potential Design with Bifold Canonical Correlation Analysis. Sensors 2022, 22, 2568. https://doi.org/10.3390/s22072568

AMA Style

Karimi R, Mohammadi A, Asif A, Benali H. DF-SSmVEP: Dual Frequency Aggregated Steady-State Motion Visual Evoked Potential Design with Bifold Canonical Correlation Analysis. Sensors. 2022; 22(7):2568. https://doi.org/10.3390/s22072568

Chicago/Turabian Style

Karimi, Raika, Arash Mohammadi, Amir Asif, and Habib Benali. 2022. "DF-SSmVEP: Dual Frequency Aggregated Steady-State Motion Visual Evoked Potential Design with Bifold Canonical Correlation Analysis" Sensors 22, no. 7: 2568. https://doi.org/10.3390/s22072568

APA Style

Karimi, R., Mohammadi, A., Asif, A., & Benali, H. (2022). DF-SSmVEP: Dual Frequency Aggregated Steady-State Motion Visual Evoked Potential Design with Bifold Canonical Correlation Analysis. Sensors, 22(7), 2568. https://doi.org/10.3390/s22072568

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