Next Article in Journal
Transformer-Based Network with Optimization for Cross-Subject Motor Imagery Identification
Previous Article in Journal
Optimized Solutions of Electrocardiogram Lead and Segment Selection for Cardiovascular Disease Diagnostics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

OptEF-BCI: An Optimization-Based Hybrid EEG and fNIRS–Brain Computer Interface

1
Department of Intelligent Mechatronics Engineering, Sejong University, Seoul 05006, Republic of Korea
2
Department of Scientific Computing, Pukyong National University, Busan 48513, Republic of Korea
3
Interdisciplinary Biology Laboratory (iBLab), Division of Biological Science, Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
4
Department of Robotics & Artificial Intelligence (R&AI), School of Mechanical and Manufacturing Engineering (SMME), National University of Sciences and Technology (NUST) H-12, Islamabad 44000, Pakistan
5
Department of Precision Medicine, School of Medicine, Sungkyunkwan University, Suwon 16419, Republic of Korea
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Bioengineering 2023, 10(5), 608; https://doi.org/10.3390/bioengineering10050608
Submission received: 21 April 2023 / Revised: 10 May 2023 / Accepted: 17 May 2023 / Published: 18 May 2023
(This article belongs to the Section Biosignal Processing)

Abstract

:
Multimodal data fusion (electroencephalography (EEG) and functional near-infrared spectroscopy (fNIRS)) has been developed as an important neuroimaging research field in order to circumvent the inherent limitations of individual modalities by combining complementary information from other modalities. This study employed an optimization-based feature selection algorithm to systematically investigate the complementary nature of multimodal fused features. After preprocessing the acquired data of both modalities (i.e., EEG and fNIRS), the temporal statistical features were computed separately with a 10 s interval for each modality. The computed features were fused to create a training vector. A wrapper-based binary enhanced whale optimization algorithm (E-WOA) was used to select the optimal/efficient fused feature subset using the support-vector-machine-based cost function. An online dataset of 29 healthy individuals was used to evaluate the performance of the proposed methodology. The findings suggest that the proposed approach enhances the classification performance by evaluating the degree of complementarity between characteristics and selecting the most efficient fused subset. The binary E-WOA feature selection approach showed a high classification rate (94.22 ± 5.39%). The classification performance exhibited a 3.85% increase compared with the conventional whale optimization algorithm. The proposed hybrid classification framework outperformed both the individual modalities and traditional feature selection classification (p < 0.01). These findings indicate the potential efficacy of the proposed framework for several neuroclinical applications.

1. Introduction

Brain–computer interfaces (BCIs) enable interaction with the environment using brain signals, without the involvement of muscles or the peripheral nervous system [1,2,3]. BCI have flourished and developed over the past few years. One of the most widely used paradigms in BCI research is motor imagery (MI) [4], which involves having the subject imagine completing a task such as gripping [5] or lifting [6]. MI-BCIs are often used to assist patients with motor function issues [7,8], either for assistance with daily living or for rehabilitation training.
According to Jeon et al. [9], MI tasks cause event-related desynchronization and synchronization in electroencephalography (EEG) [10,11]. Several feature extraction techniques have been developed to identify these activities [12,13,14]. However, because EEG is nonstationary, it is regarded as a biological signal with a low spatial resolution, high signal-to-noise ratio, and distortion [15]. Owing to these disadvantages, several studies have incorporated functional near-infrared spectroscopy (fNIRS) data to enhance the effectiveness of BCIs [16,17]. fNIRS detects the oxygen content in the blood to determine the activity level of the cerebral cortex [18]. It provides a higher spatial resolution and tolerates more motion artifacts than EEG. A 2004 study by Coyle et al. [19] used fNIRS signals to improve BCIs for the first time. The last several decades have seen a tremendous increase in the use of fNIRS in cognitive neuroscience owing to advances in brain research and cognitive neuroscience technology.
Inadequate temporal resolution and a lack of feature extraction techniques frequently result in the subpar performance of fNIRS-based BCI systems, but EEG features can be extracted efficiently using common space patterns (CSP) [20,21]. However, the fNIRS signal is distinct from the EEG signal. EEG reflects a range of high-frequency variations and the electrophysiology of the brain. In contrast, fNIRS uses a signal that varies gradually to evaluate the effect on brain hemodynamics. Integration of the complementary qualities of EEG and fNIRS may result in a high-performance BCI because the inner edge light signal and the electrophysiological signal do not interact [22,23]. Therefore, the use of the complementary information provided by both modalities, together with the appropriate feature fusion algorithms, can lead to a higher classification accuracy compared to single-modality BCI.
According to the literature [24], two categories of EEG–fNIRS fusion frameworks exist. In the first category, known as the decision level, features of both modalities are independently entered into a classifier. For instance, Fazli et al. [17] used the band power of the EEG signal, oxyhemoglobin (ΔHbO), and deoxyhemoglobin (∆HbR) as features that were individually categorized for MI tasks. A meta-classifier refined with a feedback loop was used to integrate the results of all the classifiers. The global peak cross-validation accuracy of each classifier was used in the feedback loop to identify the optimal strategy for combining their outputs. Similarly, the spectral and temporal features of hybrid EEG–fNIRS were used to classify the auditory and visual activities [25]. Both studies reported an average accuracy improvement of 5% compared with the single modality. Shargie et al. [26] trained two support vector machine (SVM) classifiers for the two modalities, respectively, and performed a decision-level fusion. They also reported a significant improvement in the classification accuracy.
The second category, known as feature-level fusion, entails selecting, altering, or concatenating information in the best possible manner prior to model training. Studies on the straightforward fusion of features have shown considerable improvement as compared to the results of a single modality. Shin et al. [27] fused EEG features (the log variance of the first and last three CSP components) and fNIRS features (average value of the time course and average slope) to classify MI (left-hand motor imagery (LHMI) vs. right-hand motor imagery (RHMI)) and mental arithmetic (MA) tasks (MA vs. baseline). This research group also designed a hybrid BCI model to categorize three activities (MA, RHMI, and resting-state activities) by combining features extracted from two EEG channels and two fNIRS optodes [28]. Their model correctly identified the three tasks with a 77.66 ± 12.1% accuracy. In a previous study [29], the authors investigated the performance of hybrid EEG and fNIRS signals in classifying a set of motor tasks. The study showed that hybrid EEG–fNIRS yielded a higher accuracy compared to that obtained using a single modality.
Sun et al. [30] developed a framework for feature fusion by employing linear, tensor, and p-th-order polynomials. For the MI and MA tasks, their proposed model attained classification accuracies of 77.53% and 90.19%, respectively. Their findings indicated that feature-level fusion improves the performance of hybrid BCIs. Choosing the best features via visual examination may be difficult, in particular when data from all the channels are utilized for feature extraction. According to the literature [31,32], selecting the right activity channels can help to improve the categorization ability. In addition, several studies have demonstrated the effectiveness of feature selection for BCIs [33,34,35]. Feature selection helps to reduce the dimensionality of the dataset, increase processing efficiency, and improve the classification accuracy. In a recent study, Zafar et al. [36] applied a graph convolutional network to determine the correlated activity channels in fNIRS data. A filter-based approach was employed to determine the optimal feature set for SVM training. The authors showed that their proposed method significantly improved the results compared to the full channel features. However, this methodology consists of two steps, channel selection and feature selection, making it less time-efficient.
The use of EEG–fNIRS multimodal techniques to enhance the classification ability has yielded encouraging results. However, powerful computational strategies for the successful incorporation of these modalities are still lacking. The techniques must consider the various feature properties that must be developed in order to maximize the integration of multimodalities. Such initiatives will likely result in significant breakthroughs in packages aimed at measuring brain activity, with an increased performance compared with a single modality. The following are the main highlights of the framework presented here:
  • First, the data acquired from both modalities were preprocessed to filter them and remove artifacts.
  • Second, the statistical temporal features of both modalities were extracted with a 10 s interval.
  • The features were fused, and the binary enhanced whale optimization algorithm (E-WOA) was used to select the optimal/efficient fused feature subset and to improve the efficiency of the multimodal characteristics by increasing their complementarity, redundancy, and significance in relation to classification labels.
  • The support-vector-machine-based cost function was used to classify the data.
  • An online MI dataset of 29 healthy individuals was used for validation.
  • Finally, the performance results of the proposed approach were compared with those of conventional WOA, other optimization algorithms, and the published literature using the same dataset.

2. Proposed Framework

Figure 1 depicts the architecture for data acquisition, preprocessing, the extraction of temporal features, multimodal feature fusion, optimal feature selection, and model training for hybrid BCI applications. The following sections explain each step in detail.

2.1. Data Acquisition

This study used EEG and fNIRS datasets available online to validate the proposed framework [27]. The dataset contained EEG and fNIRS data from 29 healthy participants. The average age of the participants was 28.5 ± 3.7 years (14 men and 15 women). None of the participants reported having any mental, neurological, or other brain disorders. The dataset included LHMI, RHMI, MA, and baseline mental tasks. The data were collected using EEG signals from the prefrontal, motor, parietal, and occipital brain areas, utilizing Fz as a ground electrode at a sampling rate of 1000 Hz. The 10-5 international system was used for the placement of the 30 EEG electrodes. Similarly, 36 fNIRS channels were set up around Fp1, Fpz, Fp2, C3, Oz, and C4, employing 14 sources and 16 detectors at a 3 cm distance using the 10-5 international system. fNIRS data were acquired at a frequency of 2.5 Hz. The dataset consisted of triggered, fNIRS, and EEG data from 6 sessions of each of 10 trials (i.e., 30 trials per task). At the start of each session, there was a 60 s rest time, followed by 20 trials (10 per task) and another 60 s rest interval. The experiment consisted of 2 s visual instructions, a 10 s activity phase, and a randomly allotted pre-rest period of 15–17 s. Figure 2 shows the experimental paradigm and the positions of the fNIRS and EEG optodes and electrodes, respectively [27].
This study used EEG and fNIRS MI datasets to validate the proposed approach. For the MI-related activities, subjects were instructed to envision their hands opening and closing while grasping a ball. The EEG and fNIRS data were down-sampled to 200 Hz and 10 Hz, respectively, for further processing.

2.2. Preprocessing

In the preprocessing phase, the common average reference was used to re-reference the EEG data [37]. The EOG and EEG data were filtered between 0.5 and 5 Hz and between 0.5 and 50 Hz after re-referencing, respectively. The EOG artifacts were eliminated using a hybrid independent component analysis–regression [38]. A third-order Butterworth band-pass filter with an 8–30 Hz cutoff frequency was applied to filter the EEG data [39]. The fNIRS data were preprocessed to eliminate physiological noise such as breathing, heartbeats, and low-frequency drift. A third-order Butterworth low-pass filter with a cutoff frequency of 0.1 Hz and a Butterworth high-pass filter with a cutoff frequency of 0.01 Hz were used to preprocess the acquired fNIRS data [40,41].

2.3. Feature Extraction

After EEG and fNIRS data preprocessing, the features for a 10 s interval were extracted. This study focused only on the most commonly used temporal statistical features, such as the slope, peak, mean, kurtosis, and skewness [42,43,44]. The peak was defined as the highest value, and curve fitting was used to determine the slope. The following equations were used to obtain the mean ( μ , Equation (1)), skewness ( S x , Equation (2)), and kurtosis ( K x , Equation (3)):
μ = 1 N k = k 1 k 2 Y ( k )
S x = E x ( Y x μ x ) 3 σ 3
K x = E x ( Y x μ x ) 4 σ 4
where Y corresponds to the input signal (EEG and fNIRS), σ to the standard deviation, and E x to the statistical expectation.
The features were extracted using all EEG and fNIRS channels. In total, 180 features (36 channels × 5 features) were extracted from the fNIRS data, whereas the EEG data contained 150 features (30 channels × 5 features). Subsequently, all features were normalized separately for both modalities. Finally, all the features from both modalities were concatenated (ffNIRS:fEEG). The final fused feature vector comprised 330 features (180 (ffNIRS) + 150 (fEEG)). An optimal feature selection approach is required to obtain the most reliable information for a robust and highly accurate EEG-fNRIS-based BCI.

2.4. Optimal Feature Selection Approach

Feature selection is a crucial step in EEG–fNIRS-based BCI systems. External devices can be controlled and environmental communication can be performed through EEG and fNIRS by measuring the electrical or hemodynamic activity of the brain. Owing to the complexity, high dimensionality, noise, and abundance of redundant and irrelevant information present in EEG and fNIRS data, feature selection is crucial. The accuracy and speed of a BCI system can be improved using feature selection, which also helps to reduce computational complexity and improve the interpretability of the results.
Filter and wrapper methods are the two primary divisions of feature selection approaches. Filter methods rank features according to their relevance and redundancy using statistical or information-theoretic metrics and then select the top-ranked features for further analysis. Wrapper approaches analyze the performance of a BCI system using several feature subsets and select the feature subset that provides the best performance using machine learning algorithms. In this study, a binary E-WOA wrapper-based feature selection method was applied to enhance the efficiency of the hybrid BCI.

2.4.1. Whale Optimization Algorithm (WOA)

The WOA is a metaheuristic optimization technique that was introduced in 2016 by Mirjalili [45]. It is based on the hunting behavior of humpback whales and is intended to address various optimization challenges. The bubble net feeding technique employed by humpback whales, whereby they ascend near to the surface and swim in a spiral pattern to catch their prey in a net of bubbles, serves as the model for the WOA [45]. The WOA keeps track of a population of potential solutions, or “whales”. The algorithm updates the whales’ location using three forms of movement: encircling prey, spiral bubble net attack, and searching for prey. The encircling of prey refers to the act of traveling in a circle around the optimal solution so as to converge on it. Spiral bubble net attack entails traveling toward the best solution available in the given moment, and searching for prey involves random movement.
Assuming that X ( t ) = ( x 1 ( t ) , x 1 ( t ) , x m ( t ) ) represents the population of whales in an N-dimensional search space, the position of each whale is represented by X i ( t ) = ( x i , 1 ( t ) , x i , 1 ( t ) , x i , D ( t ) ) . The matrix X ( 1 ) is randomly initialized for the first and subsequent iterations in the problem space. The vector X ( t ) values update themselves using the three previously defined WOA modeling methods. During optimization, WOA considers the rate of probability ( ρ ) that each X i ( t ) will alternate between encircling/searching for prey and spiral bubble net attacking strategies using the coefficient vector ( A i ( t ) ). The updated position of the whale can be computed using Equation (4):
X i ( t ) = { Encircling   prey   ( ρ i ( t ) < 0.5 )   and   ( | A i ( t ) | < 1 ) Search   for   prey   ( ρ i ( t ) < 0.5 )   and   ( | A i ( t ) | 1 ) Spiral   bubble - net   attacking   ( ρ i ( t ) 0.5 ) , 0 < ρ i ( t ) < 1
A i ( t ) = 2 × a i ( t ) × r a n d a i ( t )
where   a i ( t ) represents the linearly decreasing variable and can be computed using Equation (6):
a i ( t ) = 2 t × ( 2 M a x I t )
The mathematical model of the method of encircling prey is presented in the following equations:
                                  X i ( t + 1 ) = X b e s t ( t ) A i ( t ) × D ( t )                                                           D ( t ) = | C i ( t ) × X b e s t ( t ) X i ( t ) | C i ( t ) = 2 × r a n d }
where D ( t ) represents the distance between the current and optimal positions of the whale, and C i ( t ) represents the coefficient vector at the t iteration. The search for prey can be modeled as follows:
X i ( t + 1 ) = X r n d ( t ) A i ( t ) × D ( t )                         D ( t ) = | C i ( t ) × X r n d ( t ) X i ( t ) | }
The spiral bubble net attacking method can be modeled in turn using Equation (9), where C represents the logarithmic spiral shape:
                                              X i ( t + 1 ) = D ( t ) × exp c l × cos ( 2 π l ) + X b e s t ( t ) D ( t ) = | X b e s t ( t ) X i ( t ) | } , 1 l 1
Although the WOA is a well-known optimization method, it is still affected by early convergence, limited population diversity, and the misalignment of search strategies [46]. Therefore, Shahraki et al. [47] proposed enhanced WOA (E-WOA) and binary E-WOA variants for feature selection.

2.4.2. Enhanced WOA (E-WOA)

In E-WOA, a pooling technique and three effective search tactics—migration, preferential selection, and enriched surrounding prey—enhance the effectiveness of conventional WOAs [47]. The pooling mechanism maintains population diversity by fusing unsuccessful solutions from each iteration with a successful solution. Advanced search techniques are also included to strengthen the conventional WOA search techniques.
The E-WOA pool matrix ( P ( 1 ) , P ( 2 ) , , P ( k ) ) with the members P i = P i ( 1 ) , P i ( 2 ) , , P i ( k ) is computed using Equation (10) in the pooling mechanism at each iteration:
P i ( t ) = B i ( t ) × X b r n d ( t ) × B ¯ i ( t ) + X w o r s t ( t )
where X b r n d ( t ) represents a random position near the best whale X b e s t ( t ) , and X w o r s t ( t ) represents the worst whale at a given iteration. B i ( t ) and B ¯ i ( t ) are the random and reverse vectors, meaning that the corresponding values of zero-valued elements in B ¯ i ( t ) have a value of one and the corresponding values of non-zero-valued elements in B i ( t ) have a value of zero in B ¯ i ( t ) . The pooling technique uses a crossover operator to improve the diversity and combines the least viable solution with a promising solution. A new solution is updated using the existing pool member.
The migrating search technique randomly separates a proportion of the whales from the group using Equation (11) to cover previously unexplored areas and enhance exploration. Moreover, it is anticipated that the separated whales will boost population variety, which will help to reduce local optima trapping:
  X i ( t + 1 ) = X r n d ( t ) X b r n d ( t )         X r n d ( t ) = r a n d × ( δ max δ min ) + δ min       X b r n d ( t ) = r a n d ( δ b e s t _ max δ b e s t _ min ) + δ b e s t _ min }
where δ b e s t _ max and δ b e s t _ min represent the upper and lower boundaries of X b e s t ( t ) .
The preferential selection strategy enhances the search for prey approach. The preferential selection approach is modeled using Equation (12):
X i ( t + 1 ) = X i ( t ) + A i ( t ) × ( C i ( t ) × P r n d 1 ( t ) P r n d 2 ( t ) )
where P r n d 1 ( t ) and P r n d 2 ( t ) are randomly selected using ( P ( 1 ) , P ( 2 ) , , P ( k ) ) .
The encircling prey method from the conventional WOA is updated using the following equation:
X i ( t + 1 ) = X b e s t ( t ) A i ( t ) × D ( t )                             D ( t ) = | C i ( t ) × X b e s t ( t ) P r n d 3 ( t ) | }
Using a matrix pool, P r n d 3 ( t ) can be randomly selected. The flowchart corresponding to E-WOA is shown in Figure 3.
Shahraki et al. [47] also proposed the binary E-WOA for feature selection. The most significant/optimal features or variables associated with certain medical conditions can be successfully identified using binary optimization feature selection algorithms, making them useful for medical applications. Binary optimization feature selection algorithms assist in increasing the accuracy and efficacy of medical diagnoses and treatments by limiting the number of features or variables to those that are most important. The pseudo-code for the binary E-WOA is shown in Figure 4.

2.5. Support-Vector-Machine-Based Classification

The supervised classifier SVM has been used in several studies [42,48,49]. The main goal of SVM is to identify the optimal discriminative hyperplane that maximizes the margin between two classes. The so-called support vectors, or vectors on the margin, can be utilized to describe the ideal hyperplane. The hyperplane for a 2D feature space can be defined as follows:
f ( x ) = w x + b
where x denotes the input test data, w denotes the weight vector, and b denotes the bias factor. Projecting the test data from the input onto w, one can forecast the corresponding class output (i.e., response variable).

3. Results

As previously mentioned, this study used an open dataset of MI brain tasks to validate the proposed approach. Each trial was 10 s long and consisted of 30 trials for each task and each participant (30 in total). The performance of the temporal features of the unimodal data (EEG and ΔHbO-fNIRS) for all the channels using a 10 s interval is presented in Table 1. In total, 80% of the data were used for training, whereas the remaining 20% were reserved for the testing of the developed model. The number of features and the classification accuracy were used as comparison metrics. Equation (15) was used to compute the classification accuracy:
Classification   Accuracy   ( % ) = Correctly   classified   trail Total   no .   of   trails × 100
After a critical evaluation of the findings, it was revealed that the EEG and fNIRS data showed the highest classification accuracies (68.33% and 82.5%) for Subjects 26 and 28, respectively. The average accuracies of EEG and fNIRS were 53.59 ± 7.88 and 58.33 ± 10.13, respectively. The average accuracies of the two modalities were very low. Therefore, we applied the conventional WOA and binary E-WOA approaches to check the performance of the fused dataset. The fused feature vector contained 330 features (150 (EEG) + 180 (fNIRS)). The population selected for both optimizations was 10, with a maximum of 50 iterations. The results of the 10 runs are reported in Table 2.
As shown in Table 2, both optimization algorithms significantly enhanced the classification performance of the hybrid features from both modalities. Both wrapper-based optimization algorithms tested the effectiveness of the features by training a model with a subset of features. It is evident from Table 2 that the conventional WOA and binary E-WOA approaches enhanced the classification accuracy of Subject 1 using only 97 and 55 features, respectively. The binary E-WOA increased the classification accuracy by almost 9% compared to the conventional WOA. The classification accuracy of Subject 28 was 99.17 ± 2.64, nearly 16.5% higher than that obtained with fNIRS. Subject 28 used only 26.6 ± 13.06 optimal features from both modalities out of the total number of features (330). The box plot shown in Figure 5 summarizes the results in a simple manner.
As shown in Figure 5, the hybrid modalities outperformed the solo modalities using the conventional WOA and binary E-WOA. The conventional WOA and binary E-WOA demonstrated a high classification accuracy (90.37 ± 7.66% and 94.22 ± 5.39%, respectively). The proposed approach (binary E-WOA) only required 29.62 ± 20.78 features to train the model. In contrast, the conventional WOA required 70.81 ± 11.63 features, as depicted in Figure 5. Furthermore, the proposed approach (EEG–fNIRS-based binary E-WOA) was compared to other wrapper-based techniques [50,51,52,53,54,55]. These findings are shown in Figure 6.
The comparison of the results revealed that the EEG–fNIRS-based binary E-WOA outperformed all the other optimization algorithms in terms of the classification rate. A two-sample t-test also proved the reliability and adaptability of the results, as shown in Figure 6. Compared with the outcomes of the other optimization techniques, the binary E-WOA findings were highly accurate (p < 0.01). The proposed approach was also compared with those previously described in the literature, and the results are listed in Table 3.

4. Discussion

Multimodal hybrid techniques outperform unimodal data in terms of classification performance using conventional machine learning or deep learning techniques [59]. Although there is variability between EEG and fNIRS data, it is not as significant as initially believed based on the signal sources. Moreover, the adoption of certain effective techniques for multimodal fusion improves BCI performance. Owing to the different acquisition methods, noise sensitivities, and temporal and spatial resolutions, information from several modalities is not directly compatible. One of the main challenges in the development of hybrid BCIs is finding methods to integrate complementary information, maximize its advantages, and overcome the limitations of a single model.
Shin et al. [27] computed common spatial patterns for EEG and two statistical features for fNIRS and utilized linear discriminant analysis for MI task classification. The reported accuracy for the hybrid BCI was only 67.5%. Similarly, in another study [56], the authors developed an independent decision path fusion approach. He et al. [58] developed a neural network for a hybrid BCI model. However, the computational complexity of the model was high.
In this study, the statistical features of the two modalities were merged. Wrapper-based and binary E-WOAs were applied to select the optimal features. The results show that the binary E-WOA outperformed the conventional WOA in terms of accuracy (Table 2 and Figure 5). The ability of the binary E-WOA to detect and select more appropriate features from the hybrid EEG and fNIRS data for the classification task can be attributed to its effective search strategy. Feature selection in the conventional WOA is entirely determined by the objective function, which cannot accurately reflect the complexity and variety of brain MI datasets. This may lead to the selection of irrelevant or less-than-ideal features, which may, in turn, lead to subpar classification results, as shown in Figure 5. On the other hand, the binary E-WOA utilizes a more complex feature selection technique that incorporates the objective function and a collection of three advanced search strategies to assist in the search for optimal features, as discussed in Section 2.4.2. The three effective search strategies allow the binary E-WOA to explore the feature set more effectively and discover a more diverse collection of optimal features. This results in an improved classification accuracy and a more robust classification model. Moreover, when the results of the proposed approach were compared with those of other feature selection approaches, the binary E-WOA outperformed them and provided more robust and accurate classification results, further validating the reliability of the proposed framework. Therefore, it can be concluded that the proposed binary E-WOA hybrid EEG–fNIRS framework can be used for BCI applications.

5. Conclusions

In this study, an enhanced binary variant of the WOA was applied to improve the performance of an EEG–fNIRS-based BCI. The binary E-WOA was used to select the optimal statistical features for both modalities. The three novel search strategies enhanced the classification performance of the conventional WOA. After extensive training and testing, the proposed framework classified the MI task with a 94.22 ± 5.39% accuracy using features from both modalities. The results also suggested that the improved binary variant increased the classification accuracy for MI tasks by more than 3.85% compared to the conventional WOA. Furthermore, the results were compared with those obtained using other metaheuristic algorithms, validating the high performance of the proposed framework for classification. To demonstrate the statistical significance of the results, a two-sample t-test was performed, which revealed that the differences in accuracy, as compared with all the other approaches, were significant (p < 0.01). Therefore, the framework presented here may be useful for hybrid EEG–fNIRS BCI applications. In summary, the present study indicates that the binary E-WOA hybrid EEG–fNIRS could be used as a potential tool for various applications, including brain-controlled devices, communication, and neurorehabilitation.

Author Contributions

Conceptualization, M.U.A., K.S.K. and A.Z.; formal analysis, K.D.K. and A.Z.; funding acquisition, K.S.K. and S.W.L.; methodology, M.U.A. and K.S.K.; resources, S.W.L.; software, K.D.K.; supervision, S.W.L.; validation, A.Z.; writing—original draft, M.U.A. and K.S.K.; writing—review and editing, K.D.K., A.Z. and S.W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (2022R1C1C2003637) (to K.S.K.). This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MOE) (NRF2021R1I1A2059735) (S.W.L).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are included in the article.

Acknowledgments

This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (2022R1C1C2003637) (to K.S.K.). This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MOE) (NRF2021R1I1A2059735) (S.W.L).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nicolas-Alonso, L.F.; Gomez-Gil, J. Brain computer interfaces, a review. Sensors 2012, 12, 1211–1279. [Google Scholar] [CrossRef] [PubMed]
  2. Khosrowabadi, R.; Quek, C.; Ang, K.K.; Tung, S.W.; Heijnen, M. A Brain-Computer Interface for classifying EEG correlates of chronic mental stress. In Proceedings of the 2011 International Joint Conference on Neural Networks, San Jose, CA, USA, 31 July–5 August 2011; pp. 757–762. [Google Scholar]
  3. Hekmatmanesh, A.; Nardelli, P.H.; Handroos, H. Review of the state-of-the-art of brain-controlled vehicles. IEEE Access 2021, 9, 110173–110193. [Google Scholar] [CrossRef]
  4. Kaiser, V.; Kreilinger, A.; Müller-Putz, G.R.; Neuper, C. First steps toward a motor imagery based stroke BCI: New strategy to set up a classifier. Front. Neurosci. 2011, 5, 86. [Google Scholar] [CrossRef] [PubMed]
  5. Herath, H.M.K.K.M.B.; de Mel, W.R. Controlling an Anatomical Robot Hand Using the Brain-Computer Interface Based on Motor Imagery. Adv. Hum.-Comput. Interact. 2021, 2021, 5515759. [Google Scholar] [CrossRef]
  6. Kasemsumran, P.; Boonchieng, E. EEG-based motor imagery classification using novel adaptive threshold feature extraction and string grammar fuzzy K-nearest neighbor classification. J. Comput. 2019, 30, 27–40. [Google Scholar]
  7. Ang, K.K.; Guan, C.T.; Chua, K.S.G.; Ang, B.T.; Kuah, C.; Wang, C.C.; Phua, K.S.; Chin, Z.Y.; Zhang, H.H.; IEEE. Clinical study of neurorehabilitation in stroke using EEG-based motor imagery brain-computer interface with robotic feedback. In Proceedings of the 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Buenos Aires, Argentina, 31 August–4 September 2010; pp. 5549–5552. [Google Scholar]
  8. Lule, D.; Diekmann, V.; Kassubek, J.; Kurt, A.; Birbaumer, N.; Ludolph, A.C.; Kraft, E. Cortical plasticity in amyotrophic lateral sclerosis: Motor imagery and function. Neurorehabilit. Neural Repair 2007, 21, 518–526. [Google Scholar] [CrossRef]
  9. Jeon, Y.; Nam, C.S.; Kim, Y.J.; Whang, M.C. Event-related (De)synchronization (ERD/ERS) during motor imagery tasks: Implications for brain-computer interfaces. Int. J. Ind. Ergon. 2011, 41, 428–436. [Google Scholar] [CrossRef]
  10. Hekmatmanesh, A.; Wu, H.; Motie-Nasrabadi, A.; Li, M.; Handroos, H. Combination of discrete wavelet packet transform with detrended fluctuation analysis using customized mother wavelet with the aim of an imagery-motor control interface for an exoskeleton. Multimed. Tools Appl. 2019, 78, 30503–30522. [Google Scholar] [CrossRef]
  11. Hekmatmanesh, A.; Asl, R.M.; Wu, H.; Handroos, H. EEG control of a bionic hand with imagination based on chaotic approximation of largest Lyapunov exponent: A single trial BCI application study. IEEE Access 2019, 7, 105041–105053. [Google Scholar] [CrossRef]
  12. Dagdevir, E.; Tokmakci, M. Optimization of preprocessing stage in EEG based BCI systems in terms of accuracy and timing cost. Biomed. Signal Process. Control. 2021, 67, 102548. [Google Scholar] [CrossRef]
  13. Sadiq, M.T.; Yu, X.J.; Yuan, Z.H.; Fan, Z.M.; Rehman, A.U.; Ullah, I.; Li, G.Q.; Xiao, G.X. Motor Imagery EEG Signals Decoding by Multivariate Empirical Wavelet Transform-Based Framework for Robust Brain-Computer Interfaces. IEEE Access 2019, 7, 171431–171451. [Google Scholar] [CrossRef]
  14. Hekmatmanesh, A.; Wu, H.; Jamaloo, F.; Li, M.; Handroos, H. A combination of CSP-based method with soft margin SVM classifier and generalized RBF kernel for imagery-based brain computer interface applications. Multimed. Tools Appl. 2020, 79, 17521–17549. [Google Scholar] [CrossRef]
  15. Hallez, H.; Vanrumste, B.; Grech, R.; Muscat, J.; De Clercq, W.; Vergult, A.; D’Asseler, Y.; Camilleri, K.P.; Fabri, S.G.; Van Huffel, S.; et al. Review on solving the forward problem in EEG source analysis. J. Neuroeng. Rehabil. 2007, 4, 46. [Google Scholar] [CrossRef]
  16. Quaresima, V.; Ferrari, M. A Mini-Review on Functional Near-Infrared Spectroscopy (fNIRS): Where Do We Stand, and Where Should We Go? Photonics 2019, 6, 87. [Google Scholar] [CrossRef]
  17. Fazli, S.; Mehnert, J.; Steinbrink, J.; Curio, G.; Villringer, A.; Muller, K.R.; Blankertz, B. Enhanced performance by a hybrid NIRS-EEG brain computer interface. Neuroimage 2012, 59, 519–529. [Google Scholar] [CrossRef] [PubMed]
  18. Ferrari, M.; Quaresima, V. A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application. Neuroimage 2012, 63, 921–935. [Google Scholar] [CrossRef]
  19. Coyle, S.; Ward, T.; Markham, C.; McDarby, G. On the suitability of near-infrared (NIR) systems for next-generation brain–computer interfaces. Physiol. Meas. 2004, 25, 815. [Google Scholar] [CrossRef] [PubMed]
  20. Spüler, M.; Walter, A.; Rosenstiel, W.; Bogdan, M. Spatial filtering based on canonical correlation analysis for classification of evoked or event-related potentials in EEG data. IEEE Trans. Neural Syst. Rehabil. Eng. 2013, 22, 1097–1103. [Google Scholar] [CrossRef]
  21. Wu, D.; King, J.-T.; Chuang, C.-H.; Lin, C.-T.; Jung, T.-P. Spatial filtering for EEG-based regression problems in brain–computer interface (BCI). IEEE Trans. Fuzzy Syst. 2017, 26, 771–781. [Google Scholar] [CrossRef]
  22. Pfurtscheller, G.; Allison, B.Z.; Bauernfeind, G.; Brunner, C.; Solis Escalante, T.; Scherer, R.; Zander, T.O.; Mueller-Putz, G.; Neuper, C.; Birbaumer, N. The hybrid BCI. Front. Neurosci. 2010, 4, 3. [Google Scholar] [CrossRef] [PubMed]
  23. Dähne, S.; Biessmann, F.; Samek, W.; Haufe, S.; Goltz, D.; Gundlach, C.; Villringer, A.; Fazli, S.; Müller, K.-R. Multivariate machine learning methods for fusing multimodal functional neuroimaging data. Proc. IEEE 2015, 103, 1507–1530. [Google Scholar] [CrossRef]
  24. Deligani, R.J.; Borgheai, S.B.; McLinden, J.; Shahriari, Y. Multimodal fusion of EEG-fNIRS: A mutual information-based hybrid classification framework. Biomed. Opt. Express 2021, 12, 1635–1650. [Google Scholar] [CrossRef]
  25. Putze, F.; Hesslinger, S.; Tse, C.-Y.; Huang, Y.; Herff, C.; Guan, C.; Schultz, T. Hybrid fNIRS-EEG based classification of auditory and visual perception processes. Front. Neurosci. 2014, 8, 373. [Google Scholar] [CrossRef] [PubMed]
  26. Al-Shargie, F.; Tang, T.B.; Kiguchi, M. Stress Assessment Based on Decision Fusion of EEG and fNIRS Signals. IEEE Access 2017, 5, 19889–19896. [Google Scholar] [CrossRef]
  27. Shin, J.; von Lühmann, A.; Blankertz, B.; Kim, D.-W.; Jeong, J.; Hwang, H.-J.; Müller, K.-R. Open access dataset for EEG+ NIRS single-trial classification. IEEE Trans. Neural Syst. Rehabil. Eng. 2016, 25, 1735–1745. [Google Scholar] [CrossRef]
  28. Kwon, J.; Shin, J.; Im, C.-H. Toward a compact hybrid brain-computer interface (BCI): Performance evaluation of multi-class hybrid EEG-fNIRS BCIs with limited number of channels. PLoS ONE 2020, 15, e0230491. [Google Scholar] [CrossRef] [PubMed]
  29. Buccino, A.P.; Keles, H.O.; Omurtag, A. Hybrid EEG-fNIRS asynchronous brain-computer interface for multiple motor tasks. PLoS ONE 2016, 11, e0146610. [Google Scholar] [CrossRef]
  30. Sun, Z.; Huang, Z.; Duan, F.; Liu, Y. A Novel Multimodal Approach for Hybrid Brain–Computer Interface. IEEE Access 2020, 8, 89909–89918. [Google Scholar] [CrossRef]
  31. Zafar, A.; Ghafoor, U.; Yaqub, M.A.; Hong, K.-S. Initial-dip-based classification for fNIRS-BCI. In Proceedings of the Neural Imaging and Sensing 2019; SPIE: Bellingham, WA, USA, 2019; pp. 116–124. [Google Scholar]
  32. Asam, M.; Khan, S.H.; Akbar, A.; Bibi, S.; Jamal, T.; Khan, A.; Ghafoor, U.; Bhutta, M.R. IoT malware detection architecture using a novel channel boosted and squeezed CNN. Sci. Rep. 2022, 12, 15498. [Google Scholar] [CrossRef]
  33. Qi, Y.; Ding, F.; Xu, F.; Yang, J. Channel and Feature Selection for a Motor Imagery-Based BCI System Using Multilevel Particle Swarm Optimization. Comput. Intell. Neurosci. 2020, 2020, 8890477. [Google Scholar] [CrossRef]
  34. Naseer, N.; Hong, K.-S. fNIRS-based brain-computer interfaces: A review. Front. Hum. Neurosci. 2015, 9, 3. [Google Scholar] [CrossRef]
  35. Nazeer, H.; Naseer, N.; Khan, R.A.; Noori, F.M.; Qureshi, N.K.; Khan, U.S.; Khan, M.J. Enhancing classification accuracy of fNIRS-BCI using features acquired from vector-based phase analysis. J. Neural Eng. 2020, 17, 056025. [Google Scholar] [CrossRef]
  36. Zafar, A.; Dad Kallu, K.; Atif Yaqub, M.; Ali, M.U.; Hyuk Byun, J.; Yoon, M.; Su Kim, K. A Hybrid GCN and Filter-Based Framework for Channel and Feature Selection: An fNIRS-BCI Study. Int. J. Intell. Syst. 2023, 2023, 8812844. [Google Scholar] [CrossRef]
  37. Dien, J. Issues in the application of the average reference: Review, critiques, and recommendations. Behav. Res. Methods Instrum. Comput. 1998, 30, 34–43. [Google Scholar] [CrossRef]
  38. Mannan, M.M.N.; Jeong, M.Y.; Kamran, M.A. Hybrid ICA—Regression: Automatic identification and removal of ocular artifacts from electroencephalographic signals. Front. Hum. Neurosci. 2016, 10, 193. [Google Scholar] [CrossRef]
  39. Suk, H.I.; Lee, S.W. Subject and class specific frequency bands selection for multiclass motor imagery classification. Int. J. Imaging Syst. Technol. 2011, 21, 123–130. [Google Scholar] [CrossRef]
  40. Zafar, A.; Hong, K.-S. Neuronal activation detection using vector phase analysis with dual threshold circles: A functional near-infrared spectroscopy study. Int. J. Neural Syst. 2018, 28, 1850031. [Google Scholar] [CrossRef] [PubMed]
  41. Zafar, A.; Hong, K.-S. Reduction of onset delay in functional near-infrared spectroscopy: Prediction of HbO/HbR signals. Front. Neurorobot. 2020, 14, 10. [Google Scholar] [CrossRef]
  42. Hong, K.-S.; Khan, M.J.; Hong, M.J. Feature extraction and classification methods for hybrid fNIRS-EEG brain-computer interfaces. Front. Hum. Neurosci. 2018, 12, 246. [Google Scholar] [CrossRef] [PubMed]
  43. Hasan, M.A.H.; Khan, M.U.; Mishra, D. A Computationally Efficient Method for Hybrid EEG-fNIRS BCI Based on the Pearson Correlation. BioMed Res. Int. 2020, 2020, 1838140. [Google Scholar] [CrossRef]
  44. Zafar, A.; Hussain, S.J.; Ali, M.U.; Lee, S.W. Metaheuristic Optimization-Based Feature Selection for Imagery and Arithmetic Tasks: An fNIRS Study. Sensors 2023, 23, 3714. [Google Scholar] [CrossRef] [PubMed]
  45. Mirjalili, S.; Lewis, A. The Whale Optimization Algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar] [CrossRef]
  46. Gharehchopogh, F.S.; Gholizadeh, H. A comprehensive survey: Whale Optimization Algorithm and its applications. Swarm Evol. Comput. 2019, 48, 1–24. [Google Scholar] [CrossRef]
  47. Nadimi-Shahraki, M.H.; Zamani, H.; Mirjalili, S. Enhanced whale optimization algorithm for medical feature selection: A COVID-19 case study. Comput. Biol. Med. 2022, 148, 105858. [Google Scholar] [CrossRef]
  48. Aydin, E.A. Subject-Specific feature selection for near infrared spectroscopy based brain-computer interfaces. Comput. Methods Programs Biomed. 2020, 195, 105535. [Google Scholar] [CrossRef]
  49. Hong, K.-S.; Zafar, A. Existence of initial dip for BCI: An illusion or reality. Front. Neurorobot. 2018, 12, 69. [Google Scholar] [CrossRef]
  50. Bairathi, D.; Gopalani, D. A novel swarm intelligence based optimization method: Harris’ hawk optimization. In Proceedings of the Intelligent Systems Design and Applications: 18th International Conference on Intelligent Systems Design and Applications (ISDA 2018), Vellore, India, 6–8 December 2018; Volume 2, pp. 832–842. [Google Scholar]
  51. Kamel, S.; Hamour, H.; Ahmed, M.H.; Nasrat, L. Atom Search optimization Algorithm for Optimal Radial Distribution System Reconfiguration. In Proceedings of the 2019 International Conference on Computer, Control, Electrical, and Electronics Engineering (ICCCEEE), Khartoum, Sudan, 21–23 September 2019; pp. 1–5. [Google Scholar]
  52. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef]
  53. Fister, I.; Fister Jr, I.; Yang, X.-S.; Brest, J. A comprehensive review of firefly algorithms. Swarm Evol. Comput. 2013, 13, 34–46. [Google Scholar] [CrossRef]
  54. Yang, X.-S. (Ed.) Chapter 11—Bat Algorithms. In Nature-Inspired Optimization Algorithms, 2nd ed.; Academic Press: Cambridge, MA, USA, 2021; pp. 157–173. [Google Scholar]
  55. Ali, M.U.; Hussain, S.J.; Zafar, A.; Bhutta, M.R.; Lee, S.W. WBM-DLNets: Wrapper-Based Metaheuristic Deep Learning Networks Feature Optimization for Enhancing Brain Tumor Detection. Bioengineering 2023, 10, 475. [Google Scholar] [CrossRef]
  56. Jiang, X.; Gu, X.; Xu, K.; Ren, H.; Chen, W. Independent decision path fusion for bimodal asynchronous brain–computer interface to discriminate multiclass mental states. IEEE Access 2019, 7, 165303–165317. [Google Scholar] [CrossRef]
  57. Wang, Y.; Yang, Z.; Ji, H.; Li, J.; Liu, L.; Zhuang, J. Cross-Modal Transfer Learning from EEG to Functional Near-Infrared Spectroscopy for Classification Task in Brain-Computer Interface System. Front. Psychol. 2022, 13, 833007. [Google Scholar] [CrossRef] [PubMed]
  58. He, Q.; Feng, L.; Jiang, G.; Xie, P. Multimodal Multitask Neural Network for Motor Imagery Classification With EEG and fNIRS Signals. IEEE Sens. J. 2022, 22, 20695–20706. [Google Scholar] [CrossRef]
  59. Li, Y.; Zhang, X.; Ming, D. Early-stage fusion of EEG and fNIRS improves classification of motor imagery. Front. Neurosci. 2023, 16, 1062889. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Proposed multimodality-based framework for mental task classification.
Figure 1. Proposed multimodality-based framework for mental task classification.
Bioengineering 10 00608 g001
Figure 2. (a) Position of EEG and fNIRS electrodes and optodes, respectively. (b) Experimental paradigm.
Figure 2. (a) Position of EEG and fNIRS electrodes and optodes, respectively. (b) Experimental paradigm.
Bioengineering 10 00608 g002
Figure 3. Flowchart of E-WOA.
Figure 3. Flowchart of E-WOA.
Bioengineering 10 00608 g003
Figure 4. Pseudo-code of binary E-WOA for feature selection [47].
Figure 4. Pseudo-code of binary E-WOA for feature selection [47].
Bioengineering 10 00608 g004
Figure 5. Performance of a single modality and hybrid modalities with conventional WOA and binary E-WOA (a) in terms of the number of features and (b) in terms of accuracy.
Figure 5. Performance of a single modality and hybrid modalities with conventional WOA and binary E-WOA (a) in terms of the number of features and (b) in terms of accuracy.
Bioengineering 10 00608 g005
Figure 6. Comparison of various wrapper-based approaches. HHO: Harris hawk optimization [50], ASO: atom search optimization [51]; GWO: grey wolf optimization [52], FA: firefly algorithm [53], and BA: bat algorithm [54]. * p < 0.01.
Figure 6. Comparison of various wrapper-based approaches. HHO: Harris hawk optimization [50], ASO: atom search optimization [51]; GWO: grey wolf optimization [52], FA: firefly algorithm [53], and BA: bat algorithm [54]. * p < 0.01.
Bioengineering 10 00608 g006
Table 1. Performance of EEG and fNIRS for LHMI vs. RHMI.
Table 1. Performance of EEG and fNIRS for LHMI vs. RHMI.
SubjectEEGfNIRS
No. of FeaturesAccuracy (%)No. of FeaturesAccuracy (%)
115045.83 ± 11.9518050 ± 12.42
251.67 ± 12.9166.67 ± 14.16
352.5 ± 11.8271.67 ± 18.92
463.33 ± 12.5562.5 ± 9
557.5 ± 10.7246.67 ± 9.78
652.5 ± 6.8660 ± 12.91
750.83 ± 14.4151.67 ± 10.97
845 ± 11.9258.33 ± 18.43
961.67 ± 10.5453.33 ± 17.66
1042.5 ± 12.0865.83 ± 10.72
1155.83 ± 12.4555 ± 11.92
1247.5 ± 11.8258.33 ± 11.11
1351.67 ± 12.9165 ± 15.61
1455 ± 18.5153.33 ± 14.27
1545.83 ± 7.0851.67 ± 16.57
1658.33 ± 15.7154.17 ± 9
1742.5 ± 12.0844.17 ± 11.15
1851.67 ± 14.5964.17 ± 10.43
1957.5 ± 10.7237.5 ± 10.58
2061.67 ± 13.7265 ± 15.61
2151.67 ± 12.365.83 ± 12.7
2236.67 ± 5.8338.33 ± 9.78
2364.17 ± 14.7251.67 ± 10.24
2450.83 ± 13.2955.83 ± 4.03
2566.67 ± 14.1660.83 ± 11.15
2668.33 ± 6.5766.67 ± 14.7
2763.33 ± 5.8374.17 ± 17.32
2855.83 ± 9.6682.5 ± 7.3
2945.83 ± 15.3460.83 ± 13.64
Average 53.59 ± 7.88 58.33 ± 10.13
Table 2. Performance of the conventional WOA and binary E-WOA for the hybrid (EEG–fNIRS) dataset.
Table 2. Performance of the conventional WOA and binary E-WOA for the hybrid (EEG–fNIRS) dataset.
SubjectConventional WOABinary E-WOA
No. of FeaturesAccuracy (%)No. of FeaturesAccuracy (%)
196.8 ± 28.1182.5 ± 13.8654.4 ± 23.8991.67 ± 5.56
269.3 ± 10.3395 ± 4.322.5 ± 26.4497.5 ± 7.91
381 ± 31.0690.83 ± 10.7237.5 ± 19.4596.67 ± 4.3
477.9 ± 12.1287.5 ± 9.8235.6 ± 34.790 ± 6.57
574.2 ± 12.5582.5 ± 9.1733.6 ± 21.1194.17 ± 5.62
663.3 ± 9.6585.83 ± 11.8212.9 ± 9.2492.5 ± 4.73
767.1 ± 10.2488.33 ± 5.8323.1 ± 23.4893.33 ± 3.51
877.5 ± 14.9793.33 ± 8.6131.7 ± 15.3896.67 ± 4.3
974 ± 12.4494.17 ± 5.6258.2 ± 38.3195.83 ± 4.39
1073 ± 8.6290 ± 7.6642.6 ± 22.7892.5 ± 8.29
1173.3 ± 9.0789.17 ± 7.9118.8 ± 21.0992.5 ± 4.73
1261.4 ± 10.582.5 ± 7.313.2 ± 9.5390 ± 7.66
1374.7 ± 10.8695.83 ± 5.8929.8 ± 23.3895 ± 5.83
1459.7 ± 8.2692.5 ± 6.1522.1 ± 17.4995.83 ± 5.89
1568 ± 6.4187.5 ± 8.121.9 ± 14.990.83 ± 7.3
1672.3 ± 10.7989.17 ± 7.9130.7 ± 24.9191.67 ± 6.8
1775.5 ± 10.6295 ± 5.8329 ± 21.2994.17 ± 5.62
1865.1 ± 6.392.5 ± 4.7336.2 ± 15.595.83 ± 5.89
1969.1 ± 10.2993.33 ± 6.5730.4 ± 26.8896.67 ± 4.3
2069.3 ± 10.8894.17 ± 5.6235.8 ± 15.4595 ± 5.83
2166.9 ± 8.5490.83 ± 6.1528.9 ± 17.0794.17 ± 6.86
2266.1 ± 12.7276.67 ± 13.4922.2 ± 25.6687.5 ± 8.1
2368.7 ± 11.2195 ± 8.0534.2 ± 36.0195 ± 5.83
2473.4 ± 10.3888.33 ± 8.9632.2 ± 17.8689.17 ± 6.86
2568.2 ± 14.3192.5 ± 9.9830.9 ± 27.9495.83 ± 5.89
2669.5 ± 8.5196.67 ± 5.8327.3 ± 17.55100 ± 0
2759.9 ± 8.8197.5 ± 4.0322.5 ± 9.57100 ± 0
2870.5 ± 10.3298.33 ± 3.5126.6 ± 13.0699.17 ± 2.64
2968 ± 8.683.33 ± 8.7814.3 ± 12.7993.33 ± 5.27
Average 90.37 ± 7.66 94.22 ± 5.39
Table 3. Accuracy comparison with various studies.
Table 3. Accuracy comparison with various studies.
ReferenceYearApproachAccuracy (%)
Shin et al. [27]2016Common spatial pattern, mean, slope, shrinkage LDA67.5
Sun et al. [30]2020p-th-Order polynomial fusion77.53
Jiang et al. [56]2019Independent decision path fusion78.56
Wang et al. [57]2022R-CSP-E transfer and ensemble learning66.83
He et al. [58]2022End-to-end multimodal multitask neural network82.11
Present Study2023Temporal features and binary E-WOA94.22 ± 5.39
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ali, M.U.; Kim, K.S.; Kallu, K.D.; Zafar, A.; Lee, S.W. OptEF-BCI: An Optimization-Based Hybrid EEG and fNIRS–Brain Computer Interface. Bioengineering 2023, 10, 608. https://doi.org/10.3390/bioengineering10050608

AMA Style

Ali MU, Kim KS, Kallu KD, Zafar A, Lee SW. OptEF-BCI: An Optimization-Based Hybrid EEG and fNIRS–Brain Computer Interface. Bioengineering. 2023; 10(5):608. https://doi.org/10.3390/bioengineering10050608

Chicago/Turabian Style

Ali, Muhammad Umair, Kwang Su Kim, Karam Dad Kallu, Amad Zafar, and Seung Won Lee. 2023. "OptEF-BCI: An Optimization-Based Hybrid EEG and fNIRS–Brain Computer Interface" Bioengineering 10, no. 5: 608. https://doi.org/10.3390/bioengineering10050608

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