 Research
 Open Access
 Published:
Epileptic seizure classifications using empirical mode decomposition and its derivative
BioMedical Engineering OnLine volume 19, Article number: 10 (2020)
Abstract
Background
Epilepsy is one of the most common neurological disorders associated with disruption of brain activity. In the classification and detection of epileptic seizures, electroencephalography (EEG) measurements, which record the electrical activities of the brain, are frequently used. Empirical mode decomposition (EMD) and its derivative, ensemble EMD (EEMD) are recently developed methods used to decompose nonstationary and nonlinear signals such as EEG into a finite number of oscillations called intrinsic mode functions (IMFs). Our main objective in this study is to present a hybrid IMF selection method combining four different approaches (energy, correlation, power spectral distance, and statistical significance measures), and investigate the effect of selected IMFs extracted by EMD and EEMD on the classification. We have applied the proposed IMF selection approach on the classification of EEG signals recorded from epilepsy patients who are under treatment at our collaborator hospital. Multichannel EEG signals collected from epilepsy patients are decomposed into IMFs, and then IMF selection was performed. Finally, time and spectraldomain, and nonlinear features are extracted and feature sets are created for the classification.
Results
The maximum classification accuracies obtained using various combinations of IMFs were 94.56%, 95.63%, 96.8%, and 96.25% for SVM, KNN, naive Bayes, and logistic regression classifiers, respectively, by using EMD analysis; whereas, the EEMD approach has provided maximum classification accuracies of 96.06%, 97%, 97%, and 96.25% for SVM, KNN, naive Bayes, and logistic regression, respectively. Classification performance with the same features obtained using direct EEG signals instead of the decomposed IMFs was worse than the aforementioned 2 approaches for every combination.
Conclusion
Simulation results demonstrate that the proposed IMF selection approach affects the classification results. Also, EEMD provides a robust method for feature extraction from EEG signals in order to classify preseizure and seizure segments.
Background
Epilepsy is one of the neurological disorders associated with disruption of brain activity that affects approximately 50 million people of the world’s population [1, 2]. Detection of epileptic seizures is performed by neurologists by a visual examination of longterm electroencephalogram (EEG) signals. However, this method is very timeconsuming and generally yields incorrect results. On the other hand, epileptic seizures are initiated in different brain lobes of different individuals, so it is not possible to determine a standard focus center for the studies. Therefore longterm EEG recordings are needed to detect epileptic seizures and determine focus center [2,3,4,5]. Since visual examination of longterm EEG data makes it difficult to diagnose the disease, automatic seizure detection has become a very popular research area and various signal processing methods have been applied to solve this problem [2, 5, 6].
Many types of seizure detection and classification algorithms have been proposed in the literature [5]. These studies will be briefly discussed in "Related studies" section. In this present study, empirical mode decomposition (EMD) and its derivative, ensemble EMD (EEMD) based classification model for epileptic EEG data is introduced. Our aim is to distinguish preseizure and seizure epileptic EEG signals by classifying the features extracted from selected IMFs of EMD, or EEMD. Simulations are performed to evaluate the effectiveness of selecting the IMFs based on some metrics as opposed to using first several IMFs for the classification.
The rest of the paper is organized as follows. The review of some of the previous related work is given in "Related studies" section. Experimental results of the proposed method are shown in "Results" section. Discussion of the results is reported in "Discussion" section. The description of the data set, EMD algorithm, EEMD algorithm and the details of the proposed methodology are discussed in "Materials and methods" section.
Related studies
Epileptic seizure detection and classification studies have been reported frequently in the literature using various signal processing and classification methods. A variety of features such as temporal, spectral, statistical and nonlinear features are exploited to improve the detection and classification performance.
Several methods have been presented for the detection and classification of seizure and seizurefree EEG segments by using time and frequency domain features such as energy [7], exponential energy [8], matrix determinant [2], spectral power of Hjorth’s mobility components [9], crosscorrelation, power spectral density [10], subband spectral powers [11], average value, maximum value, and minimum value [5]. Furthermore, several studies may be found in the literature using the wavelet transform and its derivative approaches [6, 12].
Weighted multiscale Renyi permutation entropy (WMRPE), weighted permutation entropy (WPE), fuzzy entropy (FuzzyEn), a sigmoid entropy, approximate entropy (ApEn) based methods have also been frequently applied to this problem [13,14,15]. Additionally, nonlinear parameters such as fractal dimension, scaling exponent obtained with detrended fluctuation analysis (DFA), Hurst’s exponent have been utilized in many studies and successful results have been obtained for the detection and classification of seizure and seizurefree epileptic EEG signals [16, 17].
Time–frequency analysis methods such as EMD, EEMD, multivariate empirical mode decomposition (MEMD), complete ensemble empirical mode decomposition (CEEMD) which are developed for the analysis of nonlinear and nonstationary signals, have been successfully applied into detection or classification of seizure and seizurefree epileptic EEG signals in many studies [1, 18,19,20,21,22,23,24,25,26,27,28]. These methods decompose a given signal into a finite number of zero–mean oscillations called intrinsic mode functions (IMFs). One of the major problems while using EMD and other similar decomposition methods is how to choose which IMFs to be used in the classification algorithms. In most studies, the first several IMFs, known to have highfrequency oscillations, are automatically selected for feature extraction [19,20,21,22]. It may be discussed that there is a lack of methods in the literature for the selection of best IMFs while using EMD and other similar decompositions.
After a brief investigation, it may be observed that successful classification results have been obtained by using classification algorithms such as support vector machine (SVM) [5, 13,14,15, 17, 19, 21, 27], Artificial Neural Networks (ANN) [12, 20, 24], Knearest neighbor (KNN) algorithm [26], Extreme Learning Machine (ELM) [6, 16], Multilayer Perceptron Neural Network (MLPNN) [25], etc.
Results
EEG signals including preseizure and seizure segments obtained from 10channel EEG recordings of 16 epilepsy patients who are under treatment at Izmir Katip Celebi University School of Medicine, Department of Neurology, were analyzed using EMD, and EEMD approaches and various classifiers. The hybrid IMF selection process including energy, correlation, power spectral distance, and statistical significance measures was carried out for EMD and EEMD approaches in order to identify the IMFs that best represent the original signal as described in "Selection of intrinsic mode functions (IMFs)" section. After the IMF selection process, timedomain (energy, mean value, skewness, and kurtosis) and spectraldomain (total power, spectral entropy, 1st, 2nd, and 3rd moments), and nonlinear (Hurst exponent and Higuchi fractal dimension) featuresets were created using the selected three IMFs (IMF1, IMF3, IMF2) obtained by EMD, and EEMD approaches, and the EEG signal itself. In addition, we also performed simulations to compare the performance of our proposed approach with that of Discrete Wavelet Transforms (DWT). Since three selected IMFs of EMD and EEMD approaches are used for feature extraction and classification, threelevel decomposition is used for DWT utilizing Daubechies4 (db4) mother wavelet function [23]. Finally, SVM, KNN, naive Bayes, and logistic regression classifiers are used for the classification, and the results are evaluated.
Performance evaluation results of our proposed approach are given in Tables 1, 2, 3, 4. In these tables IMF1, IMF2, or IMF3; show that the features for classifications are calculated by using the corresponding IMF; IMF 1–3 denotes that the features are extracted using all three IMFs. On the other hand IMF1– IMF2 shows that the features are extracted from IMF1 and IMF2. Additionally, AC+DC13 show that the features are extracted from approximation coefficient (AC) and 3 detail coefficients (DC) of DWT. Furthermore, the italicized numbers in table cells indicate the best performance in accuracy for each approach (Tables 1, 2, 3) and classifier (Table 4).
Table 1 summarizes the performance evaluation of timedomain features used for classification. Using the timedomain features calculated from the IMF1–IMF3 (the most favorable two IMFs) of EMD, we obtain 97.18% classification accuracy and 97.14% Fscore using the logistic regression classifier. While the logistic regression algorithm yields the highest accuracy (98.13%) and Fscore (98.13%) values by using the timedomain features calculated from IMF1–IMF3 of EEMD, the SVM algorithm performs the worst (ACC: 62.44%, F1score: 60.80%) for the same features calculated from IMF2. When the same features calculated from the subbands obtained using DWT, we achieved 94.25% accuracy and 94.31% Fscore for the KNN classifier. To reveal the effect of decomposition, we analyzed the EEG signal itself and repeated the above feature extractions and classification. Using the timedomain features and KNN classifier, we obtain 89.75% accuracy and 89.96% Fscore, where the SVM performed very poorly (ACC: 53.94% and Fscore: 45.26%). Results of all classification using timedomain features are provided in Table 1.
We give the performance metrics for spectral features used in classification for different IMF combinations in Table 2. We observe that naive Bayes provides 96.88% accuracy and 96.77% Fscore using spectral features calculated from IMF1–IMF3 of EMD. However, higher classification performance is obtained by the same features calculated from IMF2–IMF3 of EEMD with logistic regression. While 95% accuracy and 94.87% Fscore were obtained from the spectral feature of DWT using naive Bayes classifier; 93.31% accuracy and 93.37% Fscore were achieved using the same feature obtained from EEG signals itself.
Classification results using nonlinear features are given in Table 3. The results suggest that the nonlinear features extracted from IMF1–4 of EMD provided classification performance with 95% accuracy and 95.01% Fscore using KNN and SVM. However, EEMD approach provided 92.94% accuracy and 92.90% Fscore using the same features with SVM. Using the features obtained from the EEG signal itself, accuracy and Fscore are obtained 69.38% and 68.95%, respectively, with KNN. On the other hand, 87.50% accuracy and 87.42% Fscore were obtained using the nonlinear feature of the DWT approach by the logistic regression classifier.
In order to determine the effect of IMF selection on the classification performance and to compare the approaches, the classification is performed with the combination of time, spectral, and nonlinear features. The classification results are shown in Table 4. In EMD approach, the SVM provided the maximum classification accuracy (94.56%) using combined features of IMF1–IMF2. However, KNN (95.63%), naive Bayes (96.88%), and logistic regression (96.25%) classifiers resulted in the highest accuracies using combined features of IMF1–IMF3.
On the other hand, in the EEMD approach SVM (96.06%) and logistic regression (96.25%) classifiers provided the highest classification accuracy for the combined features of IMF1. While KNN (97.06 %) achieves the best performance using combined features of IMF1–3, naive Bayes (97%) yielded maximum classification accuracy using the combined feature of IMF1–IMF3.
DWT approach provided maximum classification accuracy of 94.56% with naive Bayes classifier for the combined features of subbands. Notice that by using the same features extracted from the EEG signal (the last row), KNN (93.25%) provides the best classification performance. We also observed that the classification performance of the combined featureset created by using the EEG signal is worse than the EMD and EEMD approaches. Furthermore, the highest classification performance for all classifiers is achieved using features extracted by EEMD approach. Apart from the selected first 3 IMF, the success of the classification was not improved when the features obtained using the 4th IMF were included in the classification process.
In order to investigate the channelbased performance of our approaches, the classification is performed for 10 channels separately using total features of IMF1–3. The average mean classification accuracies for the channels in the left (Fp1–F7, F7–T1, T1–T3, T3–T5, Fp1–F3 channels) and right (Fp2–F8, F8–T2, T2–T4, T4–T6, Fp2–F4 channels) hemispheres are calculated. The classification accuracy of EEMD and EEGsignal based approaches are higher in the left hemisphere for all four classifiers (shown in Fig. 1b, c). These results are supported by the clinical information about epileptic focus areas of patients in our study, shown in Table 5. However, in the EMDbased approach, the classification accuracy is higher for the left hemisphere only for KNN and naive Bayes classifier (shown in Fig. 1a).
Discussion
In our proposed study, the main objective is to present a hybrid IMF selection method and explore the effect of selected IMFs extracted by EMD and EEMD, on the classification performance. Our approach investigates the advantage of using EEMD, where noiseadded versions of the signal are decomposed to eliminate the wellknown, modemixing problem of EMD. The problem of mode mixing can be described as the occurrence of very different oscillations in one mode, or very similar oscillations in different modes. EEMD method has been developed to overcome this shortcoming of EMD. As such, in our experiments we included EEMD as well as EMD to compare their classification performance.
We have applied the proposed IMF selection approach on the classification of EEG signals recorded from epilepsy patients who are under treatment at our collaborator hospital. We have used 10channel EEG signals recorded from 16 patients, providing a total of 160 preseizure, and 160 seizure (320 total) EEG segments. In addition, 4 timedomain, 5 frequency domain, and 2 nonlinear features are extracted from each selected IMF of those EEG segments. The timedomain, spectraldomain, and nonlinear features obtained from the selected three IMFs (IMF1, IMF3 and IMF2; in this order) were classified using support vector machine (SVM), Knearest neighbor (KNN), naive Bayes, and logistic regression classifiers, and the performances of EMD and EEMD approaches were compared. Then by using this selection approach, we explore the advantages of IMF selection in either EMD or EEMD approaches as opposed to using first several IMFs (IMF1–4). In order to reveal the advantages of using EMD or EEMD approaches, the same features were extracted from the EEG signal itself, and the subbands obtained by the DWT approach, and classification processes is repeated.
Performance of SVM classifier with time featureset was found to be poor for both approaches. When nonlinear featureset was used, the success of four classifiers was found to be low in both approaches. Using the spectral featureset, we obtain higher accuracies for all classifiers except logistic regression. This suggests that epileptic seizures cause distinctive changes in the frequency domain. In addition, when IMFbased classification results were evaluated, we notice that the success of classification performed only by the features obtained from the combination of selected IMFs was higher or similar to randomly selected first 4 IMFs (except nonlinear feature set). This shows that the IMF selection process helps improve the classification performance as selected IMFs carry the most useful information for the discrimination between the seizure and preseizure segments of EEG signals. The classification accuracy obtained using EMD or EEMD approaches using each featureset is higher than that of the features obtained directly from EEG signals, and subbands of DWT, for all four classifiers. The computational complexity of EMD and its derivative, over classical approaches such as DWT, and fast Fourier transform (FFT) is generally considered as a disadvantage. Contrary to common knowledge, if the number of sifting steps in the EMD algorithm is equal to 10, the computational complexity is given as \(\mathcal {O}(N\text{log}{N}),\) which is same as the computational complexity of FFT, where \(\mathcal {O}\) denotes the order of computation, and N shows the signal sample size. In addition to EMD, the number of ensembles is added to the computational complexity in the EEMD approach [29]. Therefore, in signal processing applications, EMDbased approaches may be preferred considering the tradeoff between the performance and computational cost.
Evaluating the channelbased classification performances, the classification success of the features obtained by EEMD approach was found to be higher than other approaches for all 4 classifiers (shown in Fig. 1).
The innovative contributions of our study can be highlighted as follows:
We propose a hybrid IMF selection method considering different approaches such us energy, correlation, power spectral distance, and statistical significance test.
We demonstrate the advantages of using selected IMFs by the proposed approach of either EMD or EEMD approaches as opposed to randomly selecting first several IMFs.
We investigate the performance improvement by using ensemble EMD in the classification of epileptic seizures as compared to traditional EMD, the EEG signal itself, and DWTbased approaches.
Conclusion
There are many studies in the literature for the detection and classification of epileptic seizures. Many studies have been performed in this field by using EMD and derivative approaches used in our study [1, 18,19,20,21,22,23,24,25,26,27,28]. EMD and its extensions (ensemble, multivariate and other) are suitable for the analysis of nonlinear and nonstationary signals such as EEG. In these methods, EEG signals are decomposed into IMFs which are zero–mean oscillations. Determining which of these IMFs contain useful information is vital for the success of the analysis. In most of the previous studies, the first 5 IMFs [19, 22] or first 4 IMFs [1, 17, 20, 25] have been selected, because they contain highfrequency information. In other words, no IMF selection process was performed in the initial stage of these studies. On the other hand, there are several IMF selection procedures presented in the literature based on energy, correlation coefficient, power spectrum, and statistical significance [24, 30,31,32,33]. If the signal to be analyzed contains noise, the energy and correlation coefficient of the IMFs where the noise component is dominant, will be high and misleading [30]. Therefore, the use of these IMF selection methods alone is not sufficient to determine the appropriate IMFs.
In our study, we propose a hybrid IMF selection approach considering energy, correlation, power spectral distance, and statistical significance measures. We explore the advantages of the proposed IMF selection in either EMD or EEMD approaches as opposed to using randomly selected IMFs. In our epileptic EEG classification experiments, the proposed EMD and EEMDbased approaches outperformed the EEGbased and DWTbased approaches for all classifiers and feature sets we used. The selection algorithm for both EMD and EEMD suggests IMF1, IMF3 and IMF2 in this order. We use these IMFs separately and their combinations for feature extraction and evaluate the classification performance. The classification performance of selected IMFs and their combinations was generally higher than the classification success of randomly selected IMF1–4. It is obvious that in another signal processing problem, the selection algorithm may yield a completely different set of IMFs. Hence the use of first k IMFs in the classification process, as generally done in previous studies, is not the best approach. In our simulations, highest classification accuracies were obtained by using the EEMD approach where the discriminative information about epileptic seizures in the channels may be revealed more clearly (shown in Fig. 1). Note that, working with 3 or more IMFs increases both the computational load and processing time. It may be concluded that performing an IMF selection procedure before obtaining the features directly affects the success and computational load of the study.
Materials and methods
Proposed approach
In this study, we present a method for preseizure and seizure classification algorithm using EMD and EEMDbased feature extraction methods and various classifiers as depicted in Fig. 2. EEG data recorded from diagnosed epilepsy patients are labeled by physicians, and divided into preseizure and seizure sections. These EEG segments are decomposed into intrinsic mode functions (IMFs) using both EMD and EEMD methods for each EEG channel separately. Subsequently, optimum IMFs that best represent the signal are selected by combining several selection approaches. Following the IMF selection process, temporal, spectral, statistical and nonlinear features were calculated from the selected IMFs. Finally, the extracted features were classified by using naive Bayes, Knearest neighbor (KNN), support vector machine (SVM), and logistic regression methods.
Data set
Epileptic EEG data of 16 epilepsy patients recorded using surface electrodes in Izmir Katip Celebi University, School of Medicine, Neurology Department were used in this study. EEG data were recorded using the Neurofax EEG device, from 18 different channels and a sampling frequency of 100 Hz. Surface EEG data were recorded from, Fp1–F7, F7–T1, T1–T3, T3–T5, T5–O1, Fp1–F3, F3–C3, C3–P3, P3–O1, Fp2–F8, F8–T2, T2–T4, T4–T6, T6–O2, Fp2–F4, F4–C4, C4–P4, P4–O2, electrode positions, according to the International 10–20 electrode placement system. In order to use this EEG data within the scope of our study, Izmir Katip Celebi University NonInvasive Clinical Research Ethics Committee was applied and Ethical Approval dated 08.08.2019 and numbered 296 was obtained. As discussed in [34], EEG signals recorded from the temporal and frontal lobeweighted 10 channels (Fp1–F7, F7–T1, T1–T3, T3–T5, Fp1–F3, Fp2–F8, F8–T2, T2–T4, T4–T6, Fp2–F4) are used in the study.
Oneminute preseizure and seizure epochs were marked by neurologist in the epileptic EEG signals recorded from selected channels. A total of 2 EEG epochs, one preseizure, and one seizure EEG epoch were used for each patient for our study. Thus, a total of 32 EEG epochs (containing 10 channels, for 1 min) were analyzed. Summary of the EEG data set used in the proposed study is presented in Table 5.
Analysis of EEG signals using EMD and EEMD methods
We applied empirical mode decomposition (EMD) and ensemble EMD methods for the analysis of EEG signals in our study. In the following, we present a brief introduction to these decomposition methods.
Empirical mode decomposition (EMD)
Empirical mode decomposition which produces a collection of intrinsic mode functions (IMF) with zero–mean oscillations, is used as an adaptive time–frequency signal analysis method. In nonlinear and nonstationary processes, it is applied as a feature extraction and noise reduction method in signal processing applications. It is the most important rule of the EMD method that the sum of these obtained IMFs give the original signal. It is essential for the IMF to satisfy two conditions: (1) the number of zero crossing and extrema should be equal or it varies with one, (2) the mean value of the upper and lower envelopes should be zero. The process of the EMD algorithm is to extract IMF, also called Sifting, can be performed as shown in Algorithm 1 [19, 24].
Ensemble empirical mode decomposition (EEMD)
Although the standard EMD algorithm provides successful results in signal processing applications as a time–frequency analysis method, it suffers from a problem called “mode mixing”. The problem of mode mixing can be described as the occurrence of very different oscillations in one mode, or very similar oscillations in different modes. The ensemble empirical mode decomposition (EEMD) method has been developed to overcome this problem. In the EEMD method, Gaussian white noise is added to the signal to be analyzed and the signal is decomposed into the intrinsic mode functions (IMF) using the EMD method. Due to the statistical properties of Gaussian white noise, the continuity of the signal is obtained in different frequency regions, so that the problem of mode mixing is reduced. The process of the EEMD algorithm is demonstrated in Algorithm 2 [28].
In the proposed method, we had a 10channel and twoepoch EEG signal for each patient (total number of patients is 16). Hence the size of the preseizure and seizure EEG data set was 16 × 10. Maximum numbers of obtained IMFs after applying the EMD and EEMD were 16 and 15, respectively. Therefore, since it would be timeconsuming and meaningless to obtain features from all IMFs, IMF selection process was performed before the feature extraction.
Discrete wavelet transforms (DWT) has widely been used for the analysis of nonstationary signals [23]. In our study, we use the DWTbased approach for feature extraction and classification of epileptic EEG segments to investigate the advantages of proposed EMD and EEMDbased approaches. DWT decomposes a given signal x[n] into detail and approximation coefficients by using a set of mother wavelet function [23, 35]. In our study, Daubechies4 (db4) mother wavelet and 3level subband decomposition are used.
Selection of intrinsic mode functions (IMFs)
In this study, we propose a hybrid IMF selection method by using energybased, correlationbased, PSD distancebased, and ttestbased approaches. Preseizures and seizures epileptic EEG data of 16 patients recorded from 10 channels were decomposed into the IMFs using both EMD and EEMD approaches (example signals are shown in Fig. 3), then the proposed IMF selection procedure in the following described is executed.
Energybased selection method
The energies of each IMFs are calculated as shown in Eq. (1). Since the higherenergy IMF is considered to be the best representative of the original signal, the IMFs were ranked from the high to the low energy IMF [30].
Here, \(\text{IMF}_{i}\) is the \(i{\text{th}}\;{\text{IMF}}\) and \(E_{\text{IMF}_{i}}\) is the energy of this IMF.
The correlationbased selection method
The correlation coefficient of each IMFs are calculated as shown in Eq. (2). Since the IMF with high correlation coefficient is considered to be a good representative IMF of the original signal, the IMFs are ranked from the high to low correlation coefficient IMF [31].
Here, \(C_{x,\text{IMF}_{i}}\) is the crosscovariance of the original signal and \(i{\text{th}}\;{\text{IMF}}\), \(\sigma _{x}\), and \(\sigma _{\text{IMF}_{i}}\) are the standard deviations of the original signal and \(\text{IMF}_i,\) respectively, and\(\rho\) denotes the correlation coefficient.
The PSD distancebased selection method
Another IMF selection method, based on power spectral densities (PSD) was also utilized by using the power spectral densities of the original signal and IMFs. The distances between the estimated PSDs are calculated using the Kullback Liebler distance (KLD) method as shown in Eq. (3). If the distance between the PSDs of original signal and an IMF is minimum, that IMF is considered to be the best representative IMF of the original signal. Hence, the IMFs are ranked from the low to the high PSD distance IMF [32, 33].
where \(S_x(.)\) is the power spectrum of the original signal, \(S_{\text{IMF}_{i}}(.)\) is the power spectrum of the \(i{\text{th}}\;\text{IMF}\), the \({\text{dis}}_{\text{KLD}}(x,\text{IMF}_{i})\) shows the KLD between the power spectra of the \(i{\text{th}}\;{\text{IMF}}\) and that of the original signal.
Statistical significancebased selection method
We also use the ttest statistical significance measure for the selection of best IMFs. The ttest is based on the principle of generating a null hypothesis that a single sample data set comes from a normal distribution. In this statistical significance test, test statistic values; hvalue and pvalue are calculated. Here, the hvalue indicates whether the distribution of data is normal, and the pvalue indicates the statistical significance of the data. If a pvalue greater than the specified threshold of \(\alpha\) (often chosen as 0.05 or 5% in the literature), the distribution of data can be interpreted as normal (null hypothesis is satisfied, hvalue = 0). Otherwise, if this pvalue is less than that threshold, the distribution of data may not be interpreted as normal (null hypothesis is not satisfied, hvalue = 1). The pvalues of the data whose distribution is known to be normal (hvalue = 0) can be used as a statistical significance measure. It has previously been recommended to select the IMFs with high pvalues in order to create a feature set with improved classification performance [24]. As such, we calculate the pvalue for every IMFs by applying the ttest. Since the pvalue obtained here shows the statistical significance of IMFs, the IMFs are ranked from the high to low pvalue IMF.
Table 6 shows the results of the above four selection approaches for one of the patients and one EEG channel.
These procedures were applied to the preseizure and seizure EEG data of 10 different channels of each patient separately. Finally, 40 metrics for 10 channels are calculated for each patient. All ranking matrices were combined and a 1280 × 16 dimensional ranking matrix for all preseizure and seizure EEG data was obtained. To determine the first priority selected IMFs for all signals, the histogram of the 1st column of the ranking matrix was calculated. The resulting histogram is shown in Fig. 4.
Examining the histogram shown in Fig. 4, we observe that the IMF1 is the first priority selected IMF, IMF2 is the third, and IMF3 is the second priority selected IMF. In our simulation, we choose these three IMFs (IMF1, IMF3, IMF2) for feature extraction.
The histogram shown in Fig. 4 suggests IMF1, IMF3 and IMF2 in this order.
Classification of preseizure and seizure EEG segments
In this section, we present a method to classify the preseizure and seizure segment of EEG signals collected from epilepsy patients. These EEG signals are detailed introduced in "Data set" section. We use the selected best IMFs represented the EEG signals, we extract a set of feature.
Feature extraction
Timedomain, spectral, and nonlinear features were obtained using the selected IMFs and original EEG signals to obtain feature sets.
Timedomain feature set: after the IMF selection process was carried out, the timedomain feature data set was created, using directly the EEG signals, using the first three of the IMFs obtained by EMD and EEMD methods, and using the subbands of DWT. Energy, mean value, skewness, and kurtosis values were calculated for 3 IMFs, DWT subbands, and EEG signals in the timedomain [8, 23].
$$E= \sum _{n=0}^{N1}X[n]^{2}$$(4)$$\mu=\frac{1}{N}{\sum _{n=0}^{N1}{X[n]}}$$(5)$$S = \frac{\frac{1}{N}{\sum _{n=0}^{N1}\left(X[n]\mu \right)^{3}}}{\left(\sqrt{\frac{1}{N}{\sum _{n=0}^{N1}\left([n]\mu \right)^{2}}}\right)^3}$$(6)$$K= \frac{\frac{1}{N}{\sum _{n=0}^{N1}(X[n]\mu )^{4}}}{\left(\frac{1}{N} {\sum _{n=0}^{N1}(X[n]\mu )^{2}}\right)^2}$$(7)
In the above equations, X[n] indicates the EEG signal or IMFs, N is the size of the signal or IMFs. E denotes the energy, \(\mu\) is the mean value; S indicates the skewness, K is the kurtosis value.
In the EMD and EEMDbased approaches a total of 320 × 12 size, and DWT based approach a total of 320 × 16 size feature sets were obtained. Applying the same procedure to the EEG signal itself, a total of 320 × 14 size feature set for preseizure and seizure EEG data was obtained.
Spectraldomain feature set: to generate this feature data set, the spectrum of the signal or IMF calculated by the periodogram method was used. Total power, spectral entropy, 1st, 2nd, and 3rd moments were calculated using the spectrum of signals [10, 26].
$$S(\omega _k) = \frac{1}{N}X(\omega _k)^{2}$$(8)$$S_{\text{T}}= \sum _{k=0}^{N1}{S(\omega _k)}$$(9)$$M_j= \sum _{k=0}^{N1}(\omega _k)^jS(\omega _k), \quad j= 1,2,3$$(10)$$H= \sum _{k=0}^{N1}{P(\omega _k)\log _2{P(\omega _k)}}$$(11)
Here, in Eqs. (8) and (9), \(S(\omega _k)\) denotes the power spectral density of the signal estimated by periodogram method, \(X(\omega _k)\) is the discrete Fourier transform of the signal x[n] [10], and \(S_{\text{T}}\) is the total power. In addition, N indicates the size of the corresponding signal and \(\omega _k= \frac{2\uppi }{N}k\); \(M_j\) given in Eq. (10), indicate the higher order spectral moments of the corresponding signal. H shown in Eq. (11) denotes the spectral entropy of the signal, and \(P(\omega _k)=\frac{S(\omega _k)}{S_{\text{T}}}\) indicates the normalized power spectral distribution [26].
In the EMD and EEMDbased approaches a total of 320 × 15 size, and DWT based approach a total of 320 × 20 size feature sets were obtained. Applying the same procedure to the EEG signal itself, a total of 320 × 5 size feature set for preseizure and seizure EEG data was obtained.

Nonlinear feature set: nonlinear features such as the Hurst exponent and Higuchi fractal dimension were computed to obtain this feature data set. These nonlinear features are used to analyze the complexity and selfsimilarity of brain recordings and other biological signals. Calculation of Hurst exponent and Higuchi fractal dimension were given in Eqs. (12), (13), (14), and (15); (16), (17), and (18), respectively.
$$X[n]= \{X[1], X[2],\ldots , X[N]\}$$(12)$$X_\text{A}[n]= \sum _{i=1}^n{X[i]\mu }, \quad n=1,\ldots ,N$$(13)$$\begin{aligned} R[m] & = \text{max}(\{X_A[1], X_{\text{A}}[2],\ldots , X_{\text{A}}[m]\})\text{min}(\{X_{\text{A}}[1], X_{\text{A}}[2],\ldots ,X_{\text{A}}[m]\}), \\ S[m] & = \sqrt{\frac{1}{m}{\sum _{k=1}^m(X[k]\bar{X}_m})^{2}},\quad m = 1,\ldots , N \end{aligned}$$(14)$${\text{LN}}= \ln {\frac{R(k)}{S(k)}},\quad k = 1,\ldots , N,$$(15)where X[n] given in Eq. (12) shows the EEG signal or the IMFs to be analyzed and \(\mu\) indicates the mean value of this signal. The \(X_{\text{A}}[n]\) shown in Eq. (13) indicates the accumulated deviation value of X[n]. Equations (14) shows the range series R[m] and the standard deviation S[m] of the timeseries X[n], and \(\bar{X}_m\) is the mean value from X[1] to X[m]. In Eq. (15), LN shows the logarithmic value. The Hurst exponent is calculated as the slope of the line where LN is plotted with respect to ln(k).
The value of Hurst Exponent (HE) ranges from 0 to 1. If there is no correlation in the timeseries, \(\text{HE} =0.5\); if timeseries has longrange anticorrelations, \(0< \text{HE} < 0.5\) and if there is longrange correlations in the timeseries, \(0.5< \text{HE} < 1\) [16].
Higuchi fractal dimension (HFD) used to calculate the fractal dimension (FD) directly from timeseries signals. The most important parameter that must be determined for the calculation of Higuchi fractal dimension is \(k{(\text{max})}\). The HFD values calculated in a given \(k{(\text{max})}\) range are plotted against this range in order to determine the optimal value for the \(k{(\text{max})}\) parameter. The k value that the obtained curve reaches the saturation point is determined as \(k{(\text{max})}\) [17, 36].
In Eq. (16) , X indicates the onedimensional timeseries EEG signal or the IMFs and \({X_{k}}^m\) indicates the new timeseries. Here, k and m are integers and the int(.) operation indicates the integer part of the \((Nm)/k\) value, N is the length of the original signal. The L[m, k] calculated in Eq. (17) indicates the the size of the new timeseries signals. The L[k] calculated by using the average of the L(m, k) values in Eq. (18) indicates the length of the curve for the k new time interval. HFD is calculated as the slope of the line where L[k] is plotted with respect to \(\text{ln}[1/k],k=1,2,\ldots ,k_{\text{max}}\).
In our study, HFD values calculated against different k (max) values were plotted and a graph was obtained. It was observed that this graph reached saturation point when \(k_{\text{max}} = 30.\)
In the EMD and EEMDbased approaches a total of 320 × 6 size, and DWT based approach a total of 320 × 8 size feature sets were obtained. Applying the same procedure to the EEG signal itself, a total of 320 × 2 size feature set for preseizure and seizure EEG data was obtained.
Classification
Features extracted from the selected IMFs of the EEG signals are used to discriminate the preseizure and seizure segments of the EEG by using the support vector machines (SVM), Knearest neighbor (KNN), and naive Bayes classifiers. In the following, we present the fundamentals of these classification methods.

Support vector machine (SVM): support vector machine (SVM), a supervised machine learning algorithm, is a successful algorithm that is frequently used in both classification and regression studies. In this algorithm, the elements of the data set containing n features are placed as elements of the coordinate system in an ndimensional space. Then, the classification is performed by finding the hyperplane that separates the classes best. There are many possible hyperplanes that can separate the two classes. What is important here is to choose the hyperplane from which the highest classification performance may be achieved.
Let \((x_k,y_k)\) be given as a separable sample example. Here, k indicates the size of the feature set and \(y\in \{1,1\}\) indicates the class label. Thence, separating hyperplane can be formulated with \(f(x)=\vec{w}x+c.\) Here, \(\vec{w}\) indicates the hyperplane parameters and c indicates the offset. The hyperplanes that can separate the two classes from each other with minimum error provide \(y_k[(\vec{w}x_k)+c]1\ge 0, k=1, 2,\ldots , n\) condition. The main purpose here is to achieve the maximum margin. Here, the margin is the distance between the support vectors belonging to two different class. Finally, the data falling on different sides of the hyperplane is assigned as an element of a different class [13, 14, 18, 19, 26].

Knearest neighbor (KNN): it is one of the learningbased pattern recognition methods. The data set is divided into two parts as training and tests then the learning process is performed according to the data in the training set. First, the distance between the sample to be classified and all the data in the training set is calculated. Then, the Knearest neighbors that have minimum distance is determined. Finally, the most common class among these Knearest neighbors is selected as the class of the new sample. Various distance measurement methods such as Euclidean, Manhattan, Minkowski, and Hamming can be used for distance calculation [26, 35, 37]. In our study, the most commonly used Euclidean distance calculation method is used [shown in Eq. (19)] and k value is chosen as 5.
$$\text{ED}=\sqrt{\sum _{m=1}^n{(x_my_m)^2}}$$(19) 
Naive Bayes: it is one of the probabilistic classifier based on Bayes theorem in which classification is performed according to probability basics. The classification process is performed by calculating the membership probability of a sample to all classes in the data set.
Let \(X=\{x_1,x_2,\ldots ,x_n\}\) be given. Here, n is the number of features, X indicates the sample in the featureset. In addition, \(\{M_1, M_2,\ldots ,M_m\}\) represents classes, here m is the number of classes. The probability that each X data in the data set is a member of the \(M_i\) class is calculated as given in Eq. (20):
Then the X data is assigned to the class in which class membership is highest. Here, X data is assigned to the \(M_i\) class, where \(P(M_i)\) indicates the class prior probabilities, P(X) indicates the prior probability of sample X, \(P(X/M_i)\) indicates the probability of X conditioned on \(M_i\) and \(P(M_i/X)\) indicates the probability of \(M_i\) conditioned on X [35, 37].

Logistic regression: logistic regression (LR) is a frequently used statistical classification technique in which the probability (P1), of dichotomous outcome event limited to two values such as yes/no, on/off, or 1/0, is related to a set of independent variables, and given in Eq. (21):
$$\text{logit}(P_1)= \ln \left( \frac{P_1}{1P_1}\right) =\beta _0+ \beta _1X_1+ \ldots + \beta _nX_n.$$(21)Here, \(\beta _0\) is the intercept and \(\{\beta _1X_1+ \cdots + \beta _nX_n\}\) are the coefficients associated with the independent variable \(\{X_1, X_2,\ldots , X_n\}\). Generally, in the logistic regression method, the maximum likelihood estimation (MLE) method is used to calculate the coefficients \(\{\beta _1X_1+ \cdots + \beta _nX_n\}\).
The probability of an event existing as a function of the independent variables is nonlinear as extracted from Eq. (22) [38]:
Here, \(P_1\in \{0,1\}\) indicates the probability value.
If the result of our Eq. (22) is \(\infty\), the probability is 0 (\(P_1=0\)), and if the result of this equation is \(\infty\), our probability is 1.
Performance evaluation
In this study, accuracy (ACC), sensitivity (SEN), selectivity (SPE), and precision (PRE) expressed as the performance criteria and Fscore values that is the combination of previous parameters were used for performance evaluation. Fivefold crossvalidation (CV) method has been used to establish the performances of the classifiers.
The feature set used in the kfold CV method is randomly separated into k different folds with the same size. Of these k folds, (k − 1) folds are used for training and the other one (1) fold is used for testing. No fold is used for validation processes. This process is repeated k times and the accuracy value is calculated separately for each iteration. After k iterations, the average accuracy value is obtained. This average accuracy obtained is accepted as CV accuracy [21, 23].
In Eq. 23 , falsepositive (FP) indicates the number of samples for class 0, but is mistaken for class 1 by the algorithm. Falsenegative (FN) denotes the number of samples for class 1, but is mistaken for class 0 by the algorithm. Truepositive (TP) (the number of samples for class 1) and truenegative (TN) (the number of samples for class 0) indicate the numbers of samples that are exactly classified by the algorithm [13, 14].
Availability of data and materials
The datasets generated during and/or analyzed during the current study are available from the corresponding author at reasonable request.
Abbreviations
 EEG:

Electroencephalogram
 EMD:

Empirical mode decomposition
 EEMD:

Ensemble empirical mode decomposition
 IMF:

Intrinsic mode function
 SVM:

Support vector machine
 KNN:

Knearest neighbors
References
 1.
Oweis RJ, Abdulhay EW. Seizure classification in EEG signals utilizing Hilbert–Huang transform. Biomed Eng Online. 2011;10(1):38.
 2.
Raghu S, Sriraam N, Hegde AS, Kubben PL. A novel approach for classification of epileptic seizures using matrix determinant. Expert Syst Appl. 2019;127:323–41.
 3.
Faust O, Acharya UR, Adeli H, Adeli A. Waveletbased EEG processing for computeraided seizure detection and epilepsy diagnosis. Seizure. 2015;26:56–64.
 4.
Fürbass F, Kampusch S, Kaniusas E, Koren J, Pirker S, Hopfengärtner R, Stefan H, Kluge T, Baumgartner C. Automatic multimodal detection for longterm seizure documentation in epilepsy. Clin Neurophysiol. 2017;128(8):1466–72.
 5.
Mahmoodian N, Boese A, Friebe M, Haddadnia J. Epileptic seizure detection using crossbispectrum of electroencephalogram signal. Seizure. 2019;66:4–11.
 6.
Yuan Q, Zhou W, Zhang L, Zhang F, Xu F, Leng Y, Wei D, Chen M. Epileptic seizure detection based on imbalanced classification and wavelet packet transform. Seizure. 2017;50:99–108.
 7.
Khan NA, Ali S. A new feature for the classification of nonstationary signals based on the direction of signal energy in the time–frequency domain. Comput Biol Med. 2018;100:10–6.
 8.
Fasil O, Rajesh R. Timedomain exponential energy for epileptic EEG signal classification. Neurosci Lett. 2019;694:1–8.
 9.
Kang JH, Chung YG, Kim SP. An efficient detection of epileptic seizure by differentiation and spectral analysis of electroencephalograms. Comput Biol Med. 2015;66:352–6.
 10.
Iscan Z, Dokur Z, Demiralp T. Classification of electroencephalogram signals with combined time and frequency features. Expert Syst Appl. 2011;38(8):10 499–505.
 11.
Bandarabadi M, Teixeira CA, Rasekhi J, Dourado A. Epileptic seizure prediction using relative spectral power features. Clin Neurophysiol. 2015;126(2):237–48.
 12.
Bandil MK, Wadhwani A. Multiresolution EEG and EEG subband features optimization for epileptic classification using hybrid evolutionary computing technique. Procedia Comput Sci. 2019;152:243–51.
 13.
Tawfik NS, Youssef SM, Kholief M. A hybrid automated detection of epileptic seizures in EEG records. Comput Electr Eng. 2016;53:177–90.
 14.
Xiang J, Li C, Li H, Cao R, Wang B, Han X, Chen J. The detection of epileptic seizure signals based on fuzzy entropy. J Neurosci Methods. 2015;243:18–25.
 15.
Raghu S, Sriraam N, Temel Y, Rao SV, Hegde AS, Kubben PL. Performance evaluation of DWT based sigmoid entropy in time and frequency domains for automated detection of epileptic seizures using SVM classifier. Comput Biol Med. 2019;110:127–43.
 16.
Yuan Q, Zhou W, Li S, Cai D. Epileptic EEG classification based on extreme learning machine and nonlinear features. Epilepsy Res. 2011;96(1–2):29–38.
 17.
Sharma M, Pachori RB, Acharya UR. A new approach to characterize epileptic seizures using analytic time–frequency flexible wavelet transform and fractal dimension. Pattern Recognit Lett. 2017;94:172–9.
 18.
DuqueMuñoz L, EspinosaOviedo JJ, CastellanosDominguez CG. Identification and monitoring of brain activity based on stochastic relevance analysis of shorttime EEG rhythms. Biomed Eng Online. 2014;13(1):123.
 19.
Li S, Zhou W, Yuan Q, Geng S, Cai D. Feature extraction and recognition of ictal EEG using EMD and SVM. Comput Biol Med. 2013;43(7):807–16.
 20.
Pachori RB, Patidar S. Epileptic seizure classification in EEG signals using secondorder difference plot of intrinsic mode functions. Comput Methods Progr Biomed. 2014;113(2):494–502.
 21.
Sharma R, Pachori RB. Classification of epileptic seizures in EEG signals based on phase space representation of intrinsic mode functions. Expert Syst Appl. 2015;42(3):1106–17.
 22.
Zheng Y, Wang G, Li K, Bao G, Wang J. Epileptic seizure prediction using phase synchronization based on bivariate empirical mode decomposition. Clin Neurophysiol. 2014;125(6):1104–11.
 23.
Alickovic E, Kevric J, Subasi A. Performance evaluation of empirical mode decomposition, discrete wavelet transform, and wavelet packed decomposition for automated epileptic seizure detection and prediction. Biomed Signal Process Control. 2018;39:94–102.
 24.
Zahra A, Kanwal N, Ur Rehman N, Ehsan S, McDonaldMaier KD. Seizure detection from EEG signals using multivariate empirical mode decomposition. Comput Biol Med. 2017;88:132–41.
 25.
Djemili R, Bourouba H, Korba MA. Application of empirical mode decomposition and artificial neural network for the classification of normal and epileptic EEG signals. Biocybern Biomed Eng. 2016;36(1):285–91.
 26.
Das AB, Bhuiyan MIH. Discrimination and classification of focal and nonfocal EEG signals using entropybased features in the EMDDWT domain. Biomed Signal Process Control. 2016;29:11–21.
 27.
Sikdar D, Roy R, Mahadevappa M. Epilepsy and seizure characterisation by multifractal analysis of EEG subbands. Biomed Signal Process Control. 2018;41:264–70.
 28.
Bizopoulos PA, Tsalikakis DG, Tzallas AT, Koutsouris DD, Fotiadis DI. EEG epileptic seizure detection using kmeans clustering and marginal spectrum based on ensemble empirical mode decomposition. In: 13th IEEE international conference on bioinformatics and bioengineering. New York: IEEE; 2013. p. 1–4.
 29.
Wang YH, Yeh CH, Young HWV, Hu K, Lo MT. On the computational complexity of the empirical mode decomposition algorithm. Phys A Stat Mech Appl. 2014;400:159–67.
 30.
Junsheng C, Dejie Y, Yu Y. Research on the intrinsic mode function (IMF) criterion in EMD method. Mech Syst Signal Process. 2006;20(4):817–24.
 31.
Peng Z, Peter WT, Chu F. A comparison study of improved Hilbert–Huang transform and wavelet transform: application to fault diagnosis for rolling bearing. Mech Syst Signal Process. 2005;19(5):974–88.
 32.
Lozano M, Fiz JA, Jané R. Performance evaluation of the Hilbert–Huang transform for respiratory sound analysis and its application to continuous adventitious sound characterization. Signal Process. 2016;120:99–116.
 33.
Komaty A, Boudraa AO, Augier B, DaréEmzivat D. EMDbased filtering using similarity measure between probability density functions of IMFs. IEEE Trans Instrum Meas. 2013;63(1):27–34.
 34.
Grewal S, Gotman J. An automatic warning system for epileptic seizures recorded on intracerebral EEGs. Clin Neurophysiol. 2005;116(10):2460–72.
 35.
Sharmila A, Geethanjali P. DWT based detection of epileptic seizure from EEG signals using naive bayes and kNN classifiers. IEEE Access. 2016;4:7716–27.
 36.
Gómez C, Mediavilla Á, Hornero R, Abásolo D, Fernández A. Use of the Higuchi’s fractal dimension for the analysis of MEG recordings from Alzheimer’s disease patients. Med Eng Phys. 2009;31(3):306–13.
 37.
Acharya UR, Molinari F, Sree SV, Chattopadhyay S, Ng KH, Suri JS. Automated diagnosis of epileptic EEG using entropies. Biomed Signal Process Control. 2012;7(4):401–8.
 38.
Alkan A, Koklukaya E, Subasi A. Automatic seizure detection in EEG using logistic regression and artificial neural network. J Neurosci Methods. 2005;148(2):167–76.
Acknowledgements
This study was supported by Izmir Katip Celebi University Scientific Research Projects Coordination Unit: Project numbers: 2019GAPMÜMF0003 and 2019TDRFEBE0005.
Author information
Affiliations
Contributions
EEG recordings collected from patients who were followed up by HST were used. Oneminute preseizure and seizure epochs were marked by SKA for each EEG data. AA has determined the signal processing methods to be used. OKC performed data analysis using the specified signal processing methods. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Dated 08.08.2019 and numbered 296, NonInvasive Clinical Research Ethical Approval was obtained from the Izmir Katip Celebi University Hospital.
Consent for publication
When EEG recordings were obtained from patients, prior approval was obtained for future use of these records.
Competing interests
The authors declare that they have no competing interests.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Karabiber Cura, O., Kocaaslan Atli, S., Türe, H.S. et al. Epileptic seizure classifications using empirical mode decomposition and its derivative. BioMed Eng OnLine 19, 10 (2020). https://doi.org/10.1186/s129380200754y
Received:
Accepted:
Published:
Keywords
 Electroencephalogram (EEG)
 Epilepsy
 Epileptic seizure classification
 Empirical mode decomposition
 Ensemble empirical mode decomposition
 Intrinsic mode function selection