Metaheuristic Optimization Algorithm for Signals Classification of Electroencephalography Channels

Marwa M.Eid,Fawaz Alassery,Abdelhameed Ibrahim and Mohamed Saber

1Department of Communications and Electronics,Delta Higher Institute of Engineering and Technology,Mansoura,35111,Egypt

2Department of Computer Engineering,College of Computers and Information Technology,Taif University,Taif,21944,Saudi Arabia

3Computer Engineering and Control Systems Department,Faculty of Engineering,Mansoura University,Mansoura,35516,Egypt

4Department of Electronics and Communications Engineering,Faculty of Engineering,Delta University for Science and Technology,Mansoura,11152,Egypt

Abstract: Digital signal processing of electroencephalography (EEG) data is now widely utilized in various applications, including motor imagery classification,seizure detection and prediction,emotion classification,mental task classification, drug impact identification and sleep state classification.With the increasing number of recorded EEG channels, it has become clear that effective channel selection algorithms are required for various applications.Guided Whale Optimization Method (Guided WOA), a suggested feature selection algorithm based on Stochastic Fractal Search(SFS)technique,evaluates the chosen subset of channels.This may be used to select the optimum EEG channels for use in Brain-Computer Interfaces (BCIs), the method for identifying essential and irrelevant characteristics in a dataset, and the complexity to be eliminated.This enables (SFS-Guided WOA) algorithm to choose the most appropriate EEG channels while assisting machine learning classification in its tasks and training the classifier with the dataset.The(SFSGuided WOA) algorithm is superior in performance metrics, and statistical tests such as ANOVA and Wilcoxon rank-sum are used to demonstrate this.

Keywords: Signals; metaheuristics optimization; feature selection; multilayer perceptron;support vector machines

1 Introduction

Digital signal processing is critical for several applications,including seizure detection/prediction,sleep state classification,and categorization of motor imagery.As shown in Fig.1,digital EEG signal processing consists of three components: a signal collection unit, a feature extraction unit, and a decision algorithm.The EEG signal collected from the scalp,brain surface,or brain interior is used as the system’s input.Electrodes, whether invasive or non-invasive, are used to represent the signal acquisition unit.The feature extraction unit is a signal processing device responsible for extracting distinguishing characteristics from a channel (s).For example, in a brain-computer interface (BCI),the decision unit is a hybrid unit that performs categorization,decision-making,and decision-passing to external devices that output the subject’s intention[1].

Figure 1:Processing of EEG signals

As previously stated,the interface between the brain and the computer(or another device)may be intrusive or non-invasive.While invasive technologies have recently demonstrated some promise in a variety of applications due to their high accuracy and low noise [2], noninvasive technologies are still widely used for safety applications with some additional signal processing tasks to compensate for noise and resolution limitations.Scalp EEG acquisition devices are usually chosen because they are inexpensive, simple to use, portable, and provide excellent temporal resolution.The scalp EEG waves may be recorded in a variety of modalities,including unipolar and bipolar.The former mode records the voltage differential between all electrodes and a reference electrode,where each electrodereference pair forms a channel.In contrast,the bipolar mode records the voltage differences between two designated electrodes,each pair forming a channel.The International Federation of Societies for Electroencephalography and Clinical Neurophysiology (IFSECN) proposed an electrode placement method on the scalp dubbed the International 10–20 system [3].Fig.2 illustrates the 10–20 EEG electrode locations for electrode insertion on the left and top of the head.These electrodes(channels)depict the activity of several brain regions.

Figure 2:EEG 10–20 electrode placement[3]

The brain regions are shown in Fig.3.Most of the relevant information about the functioning condition of the human brain is contained in five main brain waves,each with its own distinct frequency band.Delta band within(0–4 Hz),theta band within(3.5–7.5 Hz),alpha band within(7.5–13 Hz),beta band within (13–26 Hz), and gamma band within (26–70 Hz) are these frequency bands [4].Delta waves are associated with profound slumber.Theta waves are associated with the most meditative state(body asleep/mind awake).Alpha waves are associated with dreams and relaxation.Beta waves are the most prevalent during the waking state of intense concentration.Gamma waves are intimately connected with the brain’s decision-making process.When dealing with mental disease situations,unanticipated changes in brain waves occur, necessitating a significant amount of signal processing to diagnose aberrant conditions[4].The frequency range,speed,mental state,and waveforms of the EEG are shown in Tab.1.

Figure 3:Human brain and various lobes[4]

Table 1: EEG frequency band,speed,state,and waveforms[4]

The EEG signals collected are often multi-channel in nature.For example,we have two options while classifying these signals:work on a subset of channels chosen based on specific criteria or work on all channels [5].The method of EEG data categorization based on channel selection is shown in Fig.4.Reduce the number of channels in this signal processing environment since the setup procedure with many channels is time-consuming and inconvenient for the subject.Additionally,it increases the computational complexity of the system,which some applications need to be minimal.

Figure 4:General process of EEG signal classification

Seizure prediction and detection is another area where channel reduction may be helpful.The scientific and industry community are particularly interested in medical support systems’portable development.This can detect the onset of epileptic seizures early or even hours in advance by incorporating algorithms, thereby avoiding injury [6,7].Developing such portable systems based on computationally efficient prediction algorithms that use the fewest possible channels to reduce system power consumption is critical to system longevity.In the processing of EEG data,numerous channel selection techniques have been investigated.

The current work contribution can be summarized as follow.

1.A continuous version of the Guided Whale Optimization based on Stochastic Fractal Search algorithm(Continuous SFS-Guided WOA)is presented.

2.A binary version of the SFS-Guided WOA algorithm (Binary SFS-Guided WOA) is also presented.

3.Two publicly accessible datasets for electroencephalogram (EEG) signal processing, named BCI Competition IV-dataset 2a and BCI Competition IV-data set III,are utilized to test the suggested method.

4.The SFS-Guided WOA algorithm is employed to evaluate the chosen subset of EEG channels of the two datasets.

5.This is used to select the optimum EEG channels for use in Brain-Computer Interfaces(BCIs).

6.Statistical tests such as ANOVA and Wilcoxon rank-sum are used to demonstrate the presented method’s performance.

2 Literature Review

Feature selection methods may be categorized as filter-based, wrapper-based, or hybrid-based[8–10].The advantage of filter-based completely characteristic selection methods over traditional characteristic selection strategies is their speed and capacity to expand to large datasets.

2.1 Related Work

The optimization technique is widely used in various fields of study,including computer science,engineering[11],health[12],agriculture,and feature selection[13].The primary aim of optimization is to choose the optimal solution to a given problem among the available solutions that match the problem description.Additionally, optimization algorithms have a goal that must be reduced or maximized under the addressed problem[14–16].

Recently,numerous studies have used optimization to resolve given problems,such as the Whale Optimization Algorithm(WOA).WOA was used to locate the optimal weights for training the neural community and developed a multi-objective model of WOA,which was then applied to the problem of forecasting wind speed.Additionally,WOA was widely employed to determine the final location and length of capacitors used inside the radial system[17].Additionally,they used WOA to circumvent the difficulty of determining the final length used by a distributed generator[18],and they benefited from the use of WOA for image segmentation[19,20].

2.2 EEG Signal

The nature of EEG alerts may be very complicated since they are no longer linked, however random.The EEG dimension is determined by various factors, most notably the individual’s age,gender,psychological state,and intellectual state of the issue[21].Thus,comprehending the behavior and movement of brain cells involves various linear and nonlinear signal-processing methods that result in the physiological state and circumstances of the issue.Numerous ways are advocated for capturing the dynamic capabilities and sudden changes that may occur.The first step is preprocessing,which includes recording warnings, removing artifacts, signal averaging, output thresholding, and signal enhancement.The second stage is the function extraction technique within the procedure,which determines a feature vector from an ordinary vector[22].

2.3 Traditional Optimizer

Optimizer of Genetic Algorithm(GA)is inspired by biology(survival of the fittest).Initialization is a critical GA process.Alternatively,other genetic operators,such as elitism,may be used[23].The advantages of this optimizer are its simplicity and ability to deal with noisy fitness functions.Due to delayed convergence,premature convergence,and parameter change,complexity is not scalable.This technique is utilized in the construction of image processing filters as well as antennas.

Particle Swarm Optimization(PSO)is another optimization technique that simulates the motions and interactions of individuals in a flock of birds or a school of fish[24].Every particle is guided by its best-known location and the swarm’s optimal position.It is stable,simple to implement,and a suitable model of collaboration, but starting settings are elusive.It has a long convergence time and a high computational cost.Gene clustering,antenna design,vehicle routing,control design,and dimension reduction are only a few examples of uses.

The Grey Wolf Optimizer (GWO) is an algorithm that mimics grey wolf leadership, social structure, and hunting behavior.Encircling and assaulting the victim are the first two phases.This optimizer possesses Exploration and exploitation must be conducted in a balanced manner.While high search accuracy is simple to implement,it results in premature convergence due to the fluctuating positions of the three leaders.The greater the number of variables,the lower performance is achieved.It is utilized in feature selection, parameter adjustment of PID controllers, clustering, robotics, and route finding[25].

The foraging behaviors of humpback whales inspired the Whale Optimization Algorithm(WOA).They catch fish with bubbles as they swirl around a school of fish.It is a simple method for exploring a vast search space that is sluggish to convergence, prone to local optima stagnation, and computationally costly.WOA is applied in route planning, voltage offset reduction, and precision control of laser sensor systems[26].

3 Suggested(SFS-Guided WOA)Algorithm

The Guided WOA is a variant of the standard WOA.In the Guided WOA technique,to address the main disadvantage of this method,the search strategy for a single random whale may be substituted with an advanced design capable of quickly moving the whales toward the optimal solution or prey.The original WOA compels whales to travel randomly around one another,comparable to the global search.A whale may follow three random whales rather than one to improve exploration performance in the modified WOA(Guided WOA)[27].This may encourage whales to do more exploration while remaining unaffected by the leading position.

The Stochastic Fractal Search (SFS) technique’s diffusion process may generate a sequence of random walks around the optimum solution.This enhances the Guided WOA’s exploration capacity by using this diffusion process to find the optimal solution.Gaussian random walks are used as a component of the diffusion process that occurs around the updated optimum position.Algorithm 1 shows the continuous version of the SFS-Guided WOA algorithm.The binary conversion of the algorithm is shown in Algorithm 2, which explains step by step how to convert the continuous algorithm to a binary one to be applied for the tested EEG problem.

Algorithm 1:Continuous SFS-Guided WOA Algorithm 1:Initialize Guided WOA population Gi(i=1,2,...,n)with size n,maximum iterations Maxiter,fitness function Fn.2:Initialize Guided WOA parameters a,A,C,l,r1,r2,r3 3:Initialize Guided WOA parameters w1,w2,w3 4:Set t=1 5:Calculate fitness function Fn for each Gi 6:Find best individual G*7:while t ≤Maxiter(Termination condition)do 8: for(i=1:i ＜n+1)do 9: if(r3 ＜0.5)then 10: if(|A|＜1)then 11: Update position of current search agent as G(t+1)=G*(t)-A.D 12: else 13: Select three random search agents Grand1,Grand2,and Grand3 14: Update(z)by the exponential form of z=1-images/BZ_451_632_1924_663_1970.png t Maxiterimages/BZ_451_804_1924_835_1970.png2 15: Update position of current search agent as G(t+1)=w1*Grand1+z*w2*(Grand2-Grand3)+(1-z)*w3*(G-Grand1)16: end if 17: else 18: Update position of current search agent as G(t+1)=D′.ebl.cos(2πl(wèi))+G*(t)19: end if 20: end for 21: for(i=1:i ＜n+1)do 22: Calculate G′*i =Gaussion(μG*,σ)+(η×G*-η′×Pi)23: end for 24: Update a,A,C,1,r3 25: Calculate fitness function Fn for each Gi 26: Find best individual G*(Continued)

Algorithm 1:Continued 27: Set t=t+1 28: end while 29: return G*Algorithm 2:Binary SFS-Guided WOA Algorithm 1:Initialize SFS-Guided WOA algorithm configuration,including population and parameters 2:Change current solutions to binary solution(0 or 1)3:Evaluate fitness function and determine the best solution 4:Train k-NN based model and then calculate error 5:while t ≤itersmax do 6: Apply SF5-Guided WOA algorithm 7: Change updated solution to binary solution(0 or 1)based on the following equation G(t+1)d =■■ ■1 if 1 1+e-10(G*-0.5) ≥0.5 0 otherwise 8: Evaluate fitness function for each agent 9: Update parameters 10: Update best solution 11:end while 12:Return optimal solution

4 Results and Discussion

This section discusses the experimental results.The data preprocessing process is explained,including the description and the correlation matrix of tested EEG datasets.Configuration of the suggested algorithm is also discussed.Performance metrics and results discussion are described in detail in this part.

4.1 Data Preprocessing

Two publicly accessible datasets for electroencephalogram(EEG)signal processing are utilized in this work to test our suggested method.Tab.2 shows the description of the dataset.The Statistics of the EEG Dataset is discussed in Tab.3 and Fig.5 shows the correlation matrix of EEG dataset.The BCI Competition IV dataset is the first.The fourth BCI competition was held in 2008 at Austria’s Graz University of Technology.For the sake of this research,we will analyze dataset 2a from the competition mentioned above.This dataset is freely accessible through [28].The dataset contains the EEG data of nine healthy individuals.The subjects were healthy and ordinary people.They were instructed to complete the motor imagery activities while seated in a comfy armchair in front of an LCD display.To finish all charges,subjects used four distinct kinds of motor imagery.These activities required the use of one’s imagination to move the left or right hand,foot,or tongue.To initiate the experimental paradigm, a brief auditory beep was played.Then, after two seconds, a fixation cross appeared on the LCD and was replaced by an arrow pointing up,down,right,or left.The participants completed one of the imaging tasks involving the mouth,feet,and left or right-hand motions,depending on the orientation of the needle.The performance subject retained the chosen item’s imagination for about three seconds until the fixation cross vanished and the LCD became completely dark.Then, after a brief pause of about two seconds,the next job was resumed.This procedure was repeated 72 times for each of the four activities,totaling 288 instances of motor imagery per participant.

Table 2: Datasets description

Table 3: Description statistics of EEG dataset

Figure 5:Correlation matrix of EEG dataset

4.2 Configuration

Each dataset is subdivided into three equal-sized segments at random: training, validation, and test.During the learning phase, training is utilized to fine-tune the KNN classifier.Validation is a technique for testing.When determining the fitness function of a particular solution.Normalize data to ensure that all features are contained within the same limits and are handled equally by the machine learning model.One of the simplest methods for scaling data is to use the min-max scaler,which scales and bounds data features between 0 and 1.

4.3 Evaluation Metrics

The evaluation metrics of the suggested method and compared algorithms are shown in Tab.4.The used variables in Tab.4 are the number of optimizer’s runs,M,the best solution at the run numberj,In addition,vector size,.A number of tested points,N.Classifier label of the output for a pointi,Ci,and label of the class for a pointi,Li.Finally,the total number of features,D,andMatch,is used to calculate the matching between two inputs.

Table 4: Evaluation metrics

4.4 Results Evaluation

Results of the experimental for the two tested datasets,D1 and D2,based on the suggested and compared methods are shown in Tab.5.The results are compared to GWO, GA, WOA, and PSO algorithms.The average error of (0.161956522) for D1 and of (0.027467811) is much better based on the suggested method.The average select error of (0.385714286) and (0.61875) for D1 and D2,respectively,show the performance of the SFS-Guided WOA algorithm.Average,best,worst fitness and standard deviation show the quality of the suggested method compared to other optimization techniques.

Table 5: Experimental results of the suggested and compared methods

ANOVA and Wilcoxon Signed Rank tests are performed to confirm the suggested method compared to other algorithms.Tabs.6 and 7 show the ANOVA test results of the tested algorithms based on the first dataset(D1)and the second dataset(D2),respectively.The results indicated that the p-value is less than 0.05.Wilcoxon Signed-Rank test results based on ten runs for the first dataset(D1)and the second dataset(D2)using the suggested and compared algorithms are shown in Tabs.8 and 9,respectively.The statistical tests results confirm the performance of the SFS-Guided WOA algorithm for the EEG datasets.

Table 6: ANOVA test of the first dataset(D1)

Table 7: ANOVA test of the second dataset(D2)

Table 8: Wilcoxon signed rank test of the first datas et(D1)

Table 9: Wilcoxon signed rank test of the first dataset(D1)

The average error of the suggested (bSFS-Guided WOA) and compared algorithms (bPSO,bWOA, bGA and bGWO) over the two tested datasets (D1 and D2) is shown in Fig.6.The figure indicates the performance of the suggested method over the tested datasets.Residual,Homoscedasticity,QQ plots and heat map of the suggested and compared algorithms over the first tested dataset(D1)and the second dataset(D2)are shown in Figs.7 and 8,respectively.

Figure 7: Continued

Figure 7: Residual, Homoscedasticity, QQ plots and heat map of the suggested and compared algorithms over the first tested dataset(D1)

Figure 8: Residual, Homoscedasticity, QQ plots and heat map of the suggested and compared algorithms over the first tested dataset(D2)

5 Conclusion

In this work, the Guided Whale Optimization Method (Guided WOA) algorithm based on Stochastic Fractal Search (SFS) technique is used to evaluate the chosen subset of channels for EEG datasets.This method is used to select the optimum EEG channels for use in Brain-Computer Interfaces (BCIs).The (SFS-Guided WOA) algorithm is superior in terms of performance metrics,and statistical tests such as ANOVA and Wilcoxon rank-sum are used to demonstrate this.The results for the two tested datasets based on the suggested and compared methods(GWO,GA,WOA,and PSO algorithms) show the quality of the recommended method.The average error and average select error confirm the performance of the SFS-Guided WOA algorithm.Other metrics, such as average, best, worst fitness and standard deviation, also show the quality of the suggested method compared to other optimization techniques.The average error of the presented(bSFS-Guided WOA)algorithm and compared algorithms (bPSO, bWOA, bGA and bGWO) indicates the performance of the recommended method over the tested datasets.Residual,Homoscedasticity,QQ plots and heat map of the suggested and compared algorithms are also tested over the two datasets.The recommended method in this work will be tested for other datasets in the future.

Acknowledgement:The authors thank Taif University Accessibility Center for the study participants.We deeply acknowledge Taif University for supporting this study through Taif University Researchers Supporting Project Number(TURSP-2020/150),Taif University,Taif,Saudi Arabia.

Funding Statement:Funding for this study is received from Taif University Researchers Supporting Project No.(Project No.TURSP-2020/150),Taif University,Taif,Saudi Arabia.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.