Open access peer-reviewed chapter

Empirical Mode Decomposition of EEG Signals for the Effectual Classification of Seizures

Written By

Fasil OK and Reghunadhan Rajesh

Reviewed: August 5th, 2019 Published: September 9th, 2020

DOI: 10.5772/intechopen.89017

Chapter metrics overview

711 Chapter Downloads

View Full Metrics


Empirical mode decomposition (EMD) is a remarkable method for the analysis of nonlinear and non-stationary data. EMD will breakdown the given signal into intrinsic mode functions (IMFs), which can represent natural signals effectively. In this work, the competence of EMD with traditional features to classify the seizure and non-seizure EEG signals is studied. Due to the complex nature of human brain, the EEG signals which are recorded from different regions of brain are non-stationary in nature. Different features such as entropy features (approximate entropy (ApEn), sample entropy (SmEn), Shannon entropy (ShEn), Rényi entropy (RnEn)), fractal dimension features (Petrosian fractal dimension, Higuchi fractal dimension, Katz fractal dimension), statistical features (mean, standard deviation and energy) and exponential energy features are extracted from IMFs and fed to a SVM classifier. The performances of extracted features are studied independently. The result shows that, the EMD method is well suited for complex seizure EEG signal classification.


  • seizures
  • EEG
  • empirical mode decomposition
  • intrinsic mode functions

1. Introduction

Seizures are characterized as unexpected, unprovoked and uncontrolled explosion of electrical impulses in brain [1]. During the seizure, the patient may experiences changes in behavior, loss of consciousness, unusual movements and unusual feelings [2, 3]. The recurrent and unprovoked seizure leads to epilepsy disorder which is a prevalent neurological disorder. Epilepsy disorder will tamper the patients way of life with social stigma, work productivity lose and premature death [4].

Electroencephalogram (EEG) is one of the traditional and easiest tool for the identification and diagnosis of seizures [5]. The availability of EEG for common people within their budgetary limits made it a typical method. Due to the sophisticated nature of brain system, the EEG signals acquired from the brain are also complicated. Automated analysis of EEG signals using modern signal processing techniques might be effortless and precise for the diagnosis of seizures rather than manual approach [6].

Out of modern signal processing techniques, empirical mode decomposition (EMD) is one of the widely used techniques for the efficient interpretation of signals and images. After the introduction of EMD by Huang [7] in 1998, several studies utilized the EMD for various applications. In Nunes et al. [8] used the EMD for texture analysis and image filtering. They have used bi-dimensional EMD in their method. In another work Zeng et al. [9] applied EMD for the effective classification of gait patterns between patients with Parkinson disease and healthy subjects. In another work Hasan et al. [10] combined the deep learning methods with EMD to classify cardiovascular disease. Xiwei et al. [11] utilized the advantages of EMD in a wind speed prediction model, in which, authors used EMD for the extraction of fluctuation features of wind speed data. Another important study by Thilagaraj et al. [12] also used EMD for the identification of alcoholism.

The usefulness of the empirical mode composition for the effective understanding of the EEG signal is proven in many works in the literature. In [13], authors classified the level of autism severity from EEG with the help of EMD. They have used artificial neural network for the classification of extracted feature from intrinsic mode functions (IMFs). Two-class motor imagery EEG signals are classified in another important study based on EMD [14]. Similarly Gaur et al. [15] used multivariate empirical mode decomposition for the effective classification of multi-class BCI by analyzing EEG signals.

In this work, we have studied the effectiveness of empirical mode decomposition for the classification of seizures by analyzing EEG signals. The filtered EEG signals are segmented into 10 non-overlapping segments and decomposed into IMFs using EMD. First four IMFs are used for the feature extraction. Various features such as approximate entropy, sample entropy, Shannon entropy, Rényi entropy, exponential energy, fractal dimensional features and statistical features (mean, standard deviation and energy) are extracted from the IMFs. Support vector machine (SVM) with RBF kernel is used for classifying the seizure.

The remaining sections of the paper are as follows. A short description of EMD and algorithm is explained in Section 2. Section 3 explains the details of the dataset used in this study and in Section 4; various feature extraction methods are mentioned. In Section 5 experimental setup and results are explained. A detailed discussion of achieved results is given in Section 6 and Section 7 concludes the paper.


2. Empirical mode decomposition (EMD)

Empirical mode decomposition is a data-driven decomposition method proposed by Huang et al. for the analysis of nonlinear and non-stationary data [7], which will decomposes the signal into finite and smaller number of intrinsic mode functions (IMFs). A non-stationary signal can be represented as sum of IMFs and each IMFs should follow two conditions: (1) the number of extrema and number of zero crossing of the IMFs should be equal or differ at most by one and (2) the mean value of two envelopes defined by local maxima and local minima should be zero [16].

IMFs can be extracted from a signal through a iterative method known as shifting process as follows:

  1. Use cubic spline interpolation method to construct upper (emax) and lower (emin) envelops by connecting detected maxima and minima individually from the signal x(t).

  2. Calculate the mean mt=emax+emin2.

  3. Extract the difference d(t) between signal x(t) and calculated m1(t), dt=xtmt.

  4. Check whether d(t) is satisfying the IMFs basic conditions. Repeat step 1 to 3 until d(t) satisfying the IMFs conditions.

  5. if d(t) satisfies IMFs condition once, define the first IMF as IMF1 = d(t).

  6. The next IMFs can be obtained by generating residue r(t) as rt=xtIMF1 and use these residue as the original data for the next iteration.

  7. Iteration will stop when final residue is a function which cannot produce any more IMFs or final residue is constant/monotonic function.

    The original signal can be represented as the sum of all IMFs and final residual.


where IMFi is the ith IMF, K is the number of IMF and rK(t) is the final residual.


3. Dataset

A benchmark data set named as Bern-Barcelona EEG dataset is used in this study. The dataset includes two class EEG signals such as focal and non-focal. Each class contains 3750 pairs of signals. EEG signals in the focal class are collected from the epileptic area of the brain and non-focal signals are collected from non-epileptic area of the brain. The signals are 20 s duration with 10,240 samples in each. The signals are sampled at 512 Hz sampling rate. In our study we have used 50 signals from each class as did in many other studies [17, 18, 19].


4. Feature extraction

Feature extraction is one of the important tasks in any machine learning application. An effective and unbiased feature will provide the best results. There are several features, which are traditionally used for various EEG related studies.

Entropy features are widely used for the analysis of various non-stationary bio-signals [20, 21, 22]. Different verities of entropy are introduced in past years. In this work we have used four verities of entropy features, namely approximate entropy (ApEn), sample entropy (SmEn), Shannon entropy (ShEn) and Rényi entropy (RnEn). Among considered entropy features, approximate entropy introduced by Pincus [23] is a good measure of complexity for non-stationary signals. One of the study proposed by Hozinger et al. [24], extracted approximate entropy from ECG time-series for better understanding of electrocardiogram (ECG) signals. Another study by Ahmed et al. [25] utilized approximate entropy for surface electromyogram (EMG) signal classification. Similar to [24], they also extracted approximate entropy from direct signals with no transformation. Also, other entropy measures such as sample entropy [27, 28, 29], Shannon entropy [30, 31] and Rényi entropy [32, 33] are used in many studies.

Fractal dimension based feature are also got wide attention of researchers in recent years. The fractal dimensions are better measures of complexity of a non-linear or non-stationary data [35]. In this work we extracted three different fractal dimension features such as Petrosian fractal dimension, Higuchi fractal dimension and Katz fractal dimension. These measures are used in various EEG related studies in the literature. In a study of drowsiness detection [36], authors extracted Petrosian and Higuchi fractal dimensions from EEG time domain signals. Similarly in another work, Acharya et al. [37] extracted Katz fractal dimension with other features for the classification of various sleep stages. We have also extracted one of the newly introduced feature, namely exponential energy by Fasil and Rajesh [26]. Some of the statistical features (mean, standard deviation and energy) are also tested in this work.


5. Experiments and results

In this work, seizure EEG signals and non-seizure EEG signals are classified by decomposing the EEG signal into IMFs using empirical mode decomposition. The frequencies beyond 60 Hz are irrelevant in the EEG analysis due to the non-availability of proper information in higher frequencies [34]. A sixth order butter-worth filter is used to remove frequencies beyond 60 Hz. The signals are further segmented into 10 non-overlapping segments. Empirical mode decomposition is applied on the segmented EEG signals and first four IMFs are obtained. Feature are extracted from four IMFs and averaged across the segments. Support vector machine with RBF kernel is used for the classification task. An overall diagram of the work is given in Figure 1.

Figure 1.

Block diagram of the proposed seizure classification method.

The empirical mode decomposition produces six IMFs in total, though we have considered only first four IMFs. The reason behind this selection procedure is the non-availability of useful information in last IMFs. In this work we have extracted various features such as approximate entropy (ApEn), sample entropy (SmEn), Shannon entropy (ShEn), Rényi entropy (RnEn), Petrosian fractal dimension, Higuchi fractal dimension, Katz fractal dimension, exponential energy from four IMFs and statistical feature (mean, standard deviation and energy).

Sample focal EEG signal and IMFs obtained from focal EEG signals are shown in Figure 2. Similarly sample non-focal EEG signal and IMFs obtained from non-focal EEG signals are shown in Figure 3.

Figure 2.

Focal EEG signals and six IMFs obtained from focal EEG signal.

Figure 3.

Non-focal EEG signals and six IMFs obtained from non-focal EEG signal.

Each record in the dataset contains a pair of signals denoted as ‘x’ and ‘y’. EMD applied on both signal separately and total 8 IMFs are obtained (4 from ‘x’ and 4 from ‘y’). To investigate the ability of each feature to classify seizures, experiments are conducted on all features individually. In classification task, the capacity of RBF kernel in support vector machine is already proved in various seizure studies [26, 34]. In this work we have used RBF kernel in support vector machine for the classification task.

We have used k-fold cross-validation with k=5 for testing the ability of extracted features from IMFs. The results of the experiments are calculated with three benchmark measures, such as accuracy, sensitivity and specificity. The measured results are tabulated in Table 1.

FeatureAccuracy (%)Sensitivity (%)Specificity (%)
Shannon entropy582690
Statistical features685482
Higuchi fractal dimension717270
Sample entropy738360
Approximate entropy758268
Petrosian fractal dimension756882
Rényi entropy798276
Katz fractal dimension808872
Exponential energy848484

Table 1.

Results of various features extracted from the IMFs.

The results in Table 1 indicates that the features extracted from the IMFs of empirical mode decomposition gives promising results. Among the tested features, exponential energy feature provided better accuracy with 84%. A box-plot of extracted exponential energy is shown in Figure 4. Katz Fractal Dimension also provides better accuracy of 80% followed by Rényi entropy with 79%. Statistical features and Shannon entropy gives less accuracy out of all. It is noted that Shannon entropy giving very low sensitivity value and very high specificity, which indicates that more number of seizure signals are miss-classified as non-seizure signals.

Figure 4.

Boxplot of the extracted exponential energy feature of four IMFs of focal (green color box, labeled as F_IMF) and non-focal (blue color box, labeled as NF_IMF) EEG signals.


6. Discussion

The study of EEG signals using empirical mode decomposition (EMD) gives an insight into the effectiveness of EMD method to analyze EEG signal for seizure classification. The features (includes four types of entropy features, three types of fractal dimensions, statistical features and exponential energy) considered in this work, produces better classification accuracy when it is extracted from decomposed IMFs.

Empirical mode decomposition method decomposes the signals into various intrinsic mode functions (IMFs). Since, IMFs carries more detailed information of a signal, the features extracted from these IMFs leads to better classification.

Similar to EMD, discrete wavelet transformation (DWT) is a method, which decomposes the signal into various sub-bands [38, 39, 40, 42]. Many EEG related studies used DWT method for various analysis. Li et al. [41] combined DWT method with envelope analysis for the effective feature extraction to classify epileptic signal. In another work, Kumar et al. [42] extracted fuzzy entropy from the sub-bands of DWT for seizure detection. Similarly Liu et al. [43], Mohammadi et al. [44] and Silveira et al. [45] also used DWT method to analyze EEG signals for various purposes. Though, EMD is more better than the DWT method.

A comparison of EMD method with DWT is also carried out in this work. The same features which are extracted from the IMFs are also extracted from the DWT sub-bands and classified with same classifier. A bar chart of the comparison of classification accuracy is given in Figure 5. The comparison results show that the EMD based feature produces better classification results than DWT based features. EMD based method produced an average accuracy of 73.66%. In case of DWT the average accuracy is 68%. Although, DWT methods shows a slight improvements in results for approximate entropy and Shannon entropy features.

Figure 5.

A comparison of classification accuracy between empirical mode decomposition (EMD) and discrete wavelet transform (DWT). Red dashed vertical line indicates the average accuracy of all DWT features and green dashed vertical line indicates the average accuracy of all EMD features.

Among various entropy features, EMD-Rényi combination (79% accuracy) provides higher classification accuracy. Approximate entropy extracted from IMFs produced an accuracy of 75%. Shannon entropy with EMD is not a good choice of feature for epileptic seizure detection. The classification accuracy produced by Shannon entropy is only 58%. Complexity of EEG data is the reason for less percentage of accuracy.

Three fractal dimensions (Petrosian fractal dimension, Higuchi fractal dimension, Katz fractal dimension) used in this work also produce promising results when they are extracted from IMFs. In this study, EMD based Katz fractal dimension produces higher (80%) classification accuracy than Petrosian (75%) fractal dimension and Higuchi (71%) fractal dimension.

EMD based statistical features did not produce promising results for classification of epileptic EEG signals. But the results are comparatively better than the features from time domain and DWT domain. The highest classification accuracy (84%) reported in this study is with newly introduced exponential energy feature by Fasil and Rajesh [26]. Exponential energy feature utilizes the detailed information available in IMFs to classify epileptic EEG signals effectively. The achieved results show the effectiveness of empirical mode decomposition (EMD) as major step in epilepsy classification.


7. Conclusions

The scope of the empirical mode decomposition of EEG signals in effectual classification of seizure is studied in this work. Four intrinsic mode functions (IMFs) are obtained by applying EMD on filtered EEG signals. Widely used features such as entropy features, fractal dimension features, statistical features and exponential energy features are extracted and its discriminating power is studied. SVM with RBF kernel is used for the classification task. Exponential energy feature provided better results for the seizure classification.

Seizure identification is a challenging and risk bearing activity, which require better accuracy. In future, authors will concentrate on improving the results by incorporating other signal transformation methods with EMD.


Conflict of interest

Authors declare no conflict of interest.


  1. 1. Falco-Walter JJ, Scheffer IE, Fisher RS. The new definition and classification of seizures and epilepsy. Epilepsy Research. 2018;139:73-79
  2. 2. Saab ME, Gotman J. A system to detect the onset of epileptic seizures in scalp EEG. Clinical Neurophysiology. 2005;116(2):427-442
  3. 3. Tzimourta KD, Tzallas AT, Giannakeas N, Astrakas LG, Tsalikakis DG, Tsipouras MG. Epileptic seizures classification based on long-term EEG signal wavelet analysis. In: Precision Medicine Powered by Health and Connected Health. Singapore: Springer; 2018. pp. 165-169
  4. 4. Annegers JF, Coan SP. The risks of epilepsy after traumatic brain injury. Seizure. 2000;9(7):453-457
  5. 5. Gupta A, Singh P, Karlekar M. A novel signal modeling approach for classification of seizure and seizure-free EEG signals. IEEE Transactions on Neural Systems and Rehabilitation Engineering. 2018;26(5):925-935
  6. 6. Kumar TS, Kanhangad V, Pachori RB. Classification of seizure and seizure-free EEG signals using local binary patterns. Biomedical Signal Processing and Control. 2015;15:33-40
  7. 7. Huang NE, Shen Z, Long SR, Wu MC, Shih HH, Zheng Q, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proceedings of the Royal Society of London, Series A: Mathematical, Physical and Engineering Sciences. 1998;454(1971):903-995
  8. 8. Nunes JC, Bouaoune Y, Delechelle E, Niang O, Bunel P. Image analysis by bidimensional empirical mode decomposition. Image and Vision Computing. 2003;21(12):1019-1026
  9. 9. Zeng W, Yuan C, Wang Q, Liu F, Wang Y. Classification of gait patterns between patients with Parkinsons disease and healthy controls using phase space reconstruction (PSR), empirical mode decomposition (EMD) and neural networks. Neural Networks. 2019;111:64-76
  10. 10. Hasan NI, Bhattacharjee A. Deep learning approach to cardiovascular disease classification employing modified ECG signal from empirical mode decomposition. Biomedical Signal Processing and Control. 2019;52:128-140
  11. 11. Mi X, Liu H, Li Y. Wind speed prediction model using singular spectrum analysis, empirical mode decomposition and convolutional support vector machine. Energy Conversion and Management. 2019;180:196-205
  12. 12. Thilagaraj M, Rajasekaran MP. An empirical mode decomposition (EMD)-based scheme for alcoholism identification. Pattern Recognition Letters. 2019;125:133-139
  13. 13. Hadoush H, Alafeef M, Abdulhay E. Automated identification for autism severity level: EEG analysis using empirical mode decomposition and second order difference plot. Behavioural Brain Research. 2019;362:240-248
  14. 14. Ghritlahare R, Sahu M, Kumar R. Classification of two-class motor imagery EEG signals using empirical mode decomposition and Hilbert-Huang transformation. In: Computing and Network Sustainability. Singapore: Springer; 2019. pp. 375-386
  15. 15. Gaur P, Pachori RB, Wang H, Prasad G. A multi-class EEG-based BCI classification using multivariate empirical mode decomposition based filtering and Riemannian geometry. Expert Systems with Applications. 2018;95:201-211
  16. 16. Martis RJ, Acharya UR, Tan JH, Petznick A, Yanti R, Chua CK, et al. Application of empirical mode decomposition (EMD) for automated detection of epilepsy using EEG signals. International Journal of Neural Systems. 2012;22(06):1250027
  17. 17. Zhu G, Li Y, Wen PP, Wang S, Xi M. Epileptogenic focus detection in intracranial EEG based on delay permutation entropy. In: AIP Conference Proceedings; AIP; Vol. 1559, No. 1; 2013. pp. 31–36
  18. 18. Sharma R, Pachori RB, Gautam S. Empirical mode decomposition based classification of focal and non-focal EEG signals. In: 2014 International Conference on Medical Biometrics; IEEE; 2014. pp. 135–140
  19. 19. Sharma R, Pachori R, Acharya U. An integrated index for the identification of focal electroencephalogram signals using discrete wavelet transform and entropy measures. Entropy. 2015;17(8):5218-5240
  20. 20. Chen W, Wang Z, Xie H, Yu W. Characterization of surface EMG signal based on fuzzy entropy. IEEE Transactions on Neural Systems and Rehabilitation Engineering. 2007;15(2):266-272
  21. 21. Lv Z, Wu XP, Li M, Zhang D. A novel eye movement detection algorithm for EOG driven human computer interface. Pattern Recognition Letters. 2010;31(9):1041-1047
  22. 22. Li T, Zhou M. ECG classification using wavelet packet entropy and random forests. Entropy. 2016;18(8):285
  23. 23. Pincus SM. Approximate entropy as a measure of system complexity. Proceedings of the National Academy of Sciences. 1991;88(6):2297-2301
  24. 24. Holzinger A, Stocker C, Bruschi M, Auinger A, Silva H, Gamboa H, Fred A. On applying approximate entropy to ECG signals for knowledge discovery on the example of big sensor data. In: International Conference on Active Media Technology; Springer: Berlin, Heidelberg; 2012. pp. 646–657
  25. 25. Ahmad SA, Chappell PH. Surface EMG classification using moving approximate entropy. In: 2007 International Conference on Intelligent and Advanced Systems; IEEE. 2007. pp. 1163–1167
  26. 26. Fasil OK, Rajesh R. Time-domain exponential energy for epileptic EEG signal classification. Neuroscience Letters. 2019;694:1-8
  27. 27. Richman JS, Moorman JR. Physiological time-series analysis using approximate entropy and sample entropy. American Journal of Physiology. Heart and Circulatory Physiology. 2000;278(6):H2039-H2049
  28. 28. Jie X, Cao R, Li L. Emotion recognition based on the sample entropy of EEG. Bio-medical Materials and Engineering. 2014;24(1):1185-1192
  29. 29. Absolo D, Hornero R, Espino P, Alvarez D, Poza J. Entropy analysis of the EEG background activity in Alzheimer’s disease patients. Physiological Measurement. 2006;27(3):241
  30. 30. Cao C, Slobounov S. Application of a novel measure of EEG non-stationarity as Shannon-entropy of the peak frequency shifting for detecting residual abnormalities in concussed individuals. Clinical Neurophysiology. 2011;122(7):1314-1321
  31. 31. Yulmetyev RM, Emelyanova NA, Gafarov FM. Dynamical Shannon entropy and information Tsallis entropy in complex systems. Physica A: Statistical Mechanics and its Applications. 2004;341:649-676
  32. 32. Liang Z, Wang Y, Sun X, Li D, Voss LJ, Sleigh JW, et al. EEG entropy measures in anesthesia. Frontiers in Computational Neuroscience. 2015;9:16
  33. 33. Inuso G, La Foresta F, Mammone N, Morabito FC. Brain activity investigation by EEG processing: Wavelet analysis, kurtosis and Renyi’s entropy for artifact detection. In: 2007 International Conference on Information Acquisition; IEEE; 2007. pp. 195–200
  34. 34. Fasil OK, Rajesh R, Thasleema TM. Influence of differential features in focal and non-focal EEG signal classification. In: 2017 IEEE Region 10 Humanitarian Technology Conference (R10-HTC); IEEE; 2017. pp. 646–649
  35. 35. Accardo A, Affinito M, Carrozzi M, Bouquet F. Use of the fractal dimension for the analysis of electroencephalographic time series. Biological Cybernetics. 1997;77(5):339-350
  36. 36. Mardi Z, Ashtiani SN, Mikaili M. EEG-based drowsiness detection for safe driving using chaotic features and statistical tests. Journal of Medical Signals and Sensors. 2011;1(2):130
  37. 37. Acharya R, Faust O, Kannathal N, Chua T, Laxminarayan S. Non-linear analysis of EEG signals at various sleep stages. Computer Methods and Programs in Biomedicine. 2005;80(1):37-45
  38. 38. Ocak H. Automatic detection of epileptic seizures in EEG using discrete wavelet transform and approximate entropy. Expert Systems with Applications. 2009;36(2):2027-2036
  39. 39. Kumar Y, Dewal ML, Anand RS. Epileptic seizures detection in EEG using DWT-based ApEn and artificial neural network. Signal, Image and Video Processing. 2014;8(7):1323-1334
  40. 40. Guo L, Rivero D, Seoane JA, Pazos A. Classification of EEG signals using relative wavelet energy and artificial neural networks. In: Proceedings of the First ACM/SIGEVO Summit on Genetic and Evolutionary Computation; ACM; 2009. pp. 177–184
  41. 41. Li M, Chen W, Zhang T. Classification of epilepsy EEG signals using DWT-based envelope analysis and neural network ensemble. Biomedical Signal Processing and Control. 2017;31:357-365
  42. 42. Kumar Y, Dewal ML, Anand RS. Epileptic seizure detection using DWT based fuzzy approximate entropy and support vector machine. Neurocomputing. 2014;133:271-279
  43. 43. Liu Y, Zhou W, Yuan Q, Chen S. Automatic seizure detection using wavelet transform and SVM in long-term intracranial EEG. IEEE Transactions on Neural Systems and Rehabilitation Engineering. 2012;20(6):749-755
  44. 44. Mohammadi Z, Frounchi J, Amiri M. Wavelet-based emotion recognition system using EEG signal. Neural Computing and Applications. 2017;28(8):1985-1990
  45. 45. Silveira TD, Kozakevicius AD, Rodrigues CR. Drowsiness detection for single channel EEG by DWT best m-term approximation. Research on Biomedical Engineering. 2015;31(2):107-115

Written By

Fasil OK and Reghunadhan Rajesh

Reviewed: August 5th, 2019 Published: September 9th, 2020