Skip to main content

ORIGINAL RESEARCH article

Front. Hum. Neurosci., 21 September 2022
Sec. Brain-Computer Interfaces
Volume 16 - 2022 | https://doi.org/10.3389/fnhum.2022.1010760

Single-channel EEG signal extraction based on DWT, CEEMDAN, and ICA method

  • 1School of Computer Science and Engineering, Guilin University of Aerospace Technology, Guilin, China
  • 2School of Computer Science and Information Security, Guilin University of Electronic Technology, Guilin, China
  • 3School of Electronic Information and Automation, Guilin University of Aerospace Technology, Guilin, China

In special application scenarios, such as portable anesthesia depth monitoring, portable emotional state recognition and portable sleep monitoring, electroencephalogram (EEG) signal acquisition equipment is required to be convenient and easy to use. It is difficult to remove electrooculogram (EOG) artifacts when the number of EEG acquisition channels is small, especially when the number of observed signals is less than that of the source signals, and the overcomplete problem will arise. The independent component analysis (ICA) algorithm commonly used for artifact removal requires the number of basis vectors to be smaller than the dimension of the input data due to a set of standard orthonormal bases learned during the convergence process, so it cannot be used to solve the overcomplete problem. The empirical mode decomposition method decomposes the signal into several independent intrinsic mode functions so that the number of observed signals is more than that of the source signals, solving the overcomplete problem. However, when using this method to solve overcompleteness, the modal aliasing problem will arise, which is caused by abnormal events such as sharp signals, impulse interference, and noise. Aiming at the above problems, we propose a novel EEG artifact removal method based on discrete wavelet transform, complete empirical mode decomposition for adaptive noise (CEEMDAN) and ICA in this paper. First, the input signals are transformed by discrete wavelet (DWT), and then CEEMDAN is used to solve the overcomplete and mode aliasing problems, meeting the a priori conditions of the ICA algorithm. Finally, the components belonging to EOG artifacts are removed according to the sample entropy value of each independent component. Experiments show that this method can effectively remove EOG artifacts while solving the overcomplete and modal aliasing problems.

Introduction

EEG (Saeidi et al., 2021) is a non-linear and non-stationary electrophysiological signal of the central nervous system that contains rich information about brain activity. It is an important information source for human brain research, disease diagnosis and rehabilitation engineering. It also has rich rhythmic activity and can be used to characterize the dynamic changes in brain function. EEG signal acquisition equipment is required to be more portable and easier to use with the application of brain computer interfaces (BCIs). The fewer the number of channels is, the better in some special application scenarios, such as portable anesthesia depth monitoring, portable emotional state recognition (Wan et al., 2022) and portable sleep monitoring (Kwon et al., 2021), in which only a single channel is needed. However, EEG signals are easily contaminated by physiological and non-physiological artifacts, such as EOG (Miao et al., 2021), electromyogram (EMG) (Meng et al., 2022) and electrocardiogram (ECG) artifacts (Mourad, 2022). As a result, the brain function information in EEG signals is concealed by artifacts, which results in inaccurate classification. EOG artifacts have the highest amplitude and strongest randomness compared with the other two kinds of artifacts. The presence of artifacts makes EEG signals prone to obvious distortion (Jiang et al., 2019; Gu et al., 2021), interferes with the inherent information expression of neuronal electrical activity, weakens the signal-to-noise ratio and increases the difficulty of preprocessing in the recording process (Dora and Biswal, 2020; Liu et al., 2021; Sun et al., 2021).

In recent years, many classical methods, such as average artifact regression analysis (Semlitsch et al., 1986), principal component analysis (PCA) (Vigon et al., 2000), and ICA (Makeig et al., 1996; Jiang et al., 2020; Yuan et al., 2020), have been widely used to remove EOG artifacts from multichannel EEG signals. In the first method, it is assumed that the conductivity between the electrode collecting EOG signal and other electrodes remains unchanged. Then, the correlation between the EOG channel and other channels is estimated, and the EOG signals are removed from each channel according to the conductivity. In the PCA method, the EEG and EOG signals need to be recorded at the same time during the experiment, and the artifacts are removed by analyzing the main components. Different from the first method, the ICA method obtains each independent source signal and separates the features only according to its statistical characteristics and the observed signal when the source signal and transmission channel parameters are unknown (Li et al., 2013).

The above methods are all difficult to apply to a single-channel portable BCI system, as they require more EEG channels to achieve a better separation effect. Kumar proposed a method that removes single-channel EOG artifacts by using a wavelet transform soft threshold (Kumar et al., 2008). However, this method requires much prior knowledge and is highly subjective, and directly removing the wavelet coefficients will negatively impact the components of the EEG source signal. Mammone proposed a WT-ICA algorithm to remove EOG artifacts based on both wavelet transform (WT) and ICA (Mammone et al., 2012). This algorithm meets the ICA a priori condition by wavelet decomposition of the single-channel EEG signal. However, the wavelet transform not only increases the observation signals but also decomposes the source signal into several subsource signals, which leads to the overcomplete problem. The empirical mode decomposition (EMD) (Wu and Huang, 2011) proposed by Huang is an adaptive time-frequency decomposition method that can better deal with non-linear and non-stationary signals (Park et al., 2011). However, there will be intermittent phenomena due to abnormal events (such as tip signals, pulse interference and noise) in the EMD process, which results in mode aliasing and causes the IMF component to lose its specific physical significance.

Aiming at the problem that the ICA method cannot be used to solve the overcomplete and modal aliasing problems caused by using empirical mode decomposition to make the number of observed signals greater than that of the source signals in the process of single-channel EOG artifact removal, we propose a single-channel EOG artifact removal algorithm (DWT-CEEMDAN-ICA) based on DWT, CEEMDAN and ICA in this paper. First, the source signals are transformed by a discrete wavelet and decomposed by the CEEMDAN method so that they meet the a priori condition of the ICA algorithm, which solves the overcomplete and mode aliasing problems. Experiments show the effectiveness and stability of this algorithm compared with the other EOG artifact removal algorithms.

This paper is organized as follows. Related work is introduced in Section 2. In Section DWT-CEEMDAN-ICA algorithm, the proposed method is described in detail. Experiments are presented in Section 4, and some concluding remarks are given in the last section.

Related works

Discrete wavelet transform

The idea of wavelet transform is to gradually refine the signal at multiple scales through expansion and translation operations so that it can be subdivided according to time in the high-frequency domain and subdivided according to frequency in the low-frequency domain. The EEG signal after wavelet transform has better frequency resolution corresponding to the two different domains, which automatically meets the requirements of time-frequency signal analysis and focuses on any detail. The window function is given by:

ψa,b(t)=1aψ(t-ba)    (1)

Where a and b are the scale displacement and time displacement, respectively.

In the DWT method, the parameter α, τ of wavelet basis function ψ(α, τ) needs to be limited to discrete points. The basis function is:

ψj,k(t)=2-j2ψ(2-jt-kτ0)    (2)

Where j and k are the frequency resolution and time translation, respectively; then, the DWT at this time is:

Wψf(j,k)= f(t)ψj,k*(t)dt    (3)

Complete empirical mode decomposition

The idea of Empirical mode decomposition (EMD) is to decompose the signal step by step according to the fluctuation or trend to produce a series of data sequences with different characteristic scales. Each sequence is called an intrinsic mode function (IMF) (Boda et al., 2021) and meets the following conditions: The number of extreme points and zero crossings in the whole data segment is equal or has no more than one difference. The average value of the upper and lower packet routes formed by the local maximum and minimum points is zero and locally symmetrical about the time axis at any time. However, due to the large amount of noise, jumping change of the time scale and boundary effect in the actual signal, the phenomenon of mode aliasing will be caused in the process of EMD. The EMD formula is as follows:

x(t)=i=1Nan(t)+rn(t)    (4)

Where an(t) is the nth-order IMF, rn(t) is the remainder, and N is the number of IMFs.

To solve the problem of mode aliasing in the EMD method, Torres et al. proposed a complete empirical mode decomposition (CEEMDAN) algorithm that can adapt to noise (Xu et al., 2018). White noise is added to the residual value, and the mean value is calculated for each IMF component and then iterated step by step.

The method is described as follows:

EEG signals and EOG artifacts partially coincide in the time domain and frequency domain of the low-frequency band. If the artifacts are directly removed from the result of the discrete wavelet transform, part of the EEG signal will be lost, resulting in distortion. In addition, the a priori condition of the ICA method will no longer be met. Therefore, we use the Algorithm 1 (CEEMDAN method) to decompose the wavelet coefficients after wavelet transform into several IMFs.

ALGORITHM 1
www.frontiersin.org

Algorithm 1. CEEMDAN

Independent component analysis and sample entropy

The original signal collected through the electrode is the linear instantaneous mix of EEG signals and EOG artifacts, which are independent of each other. Therefore, the ICA method can be used to decompose the original signal into multiple independent component spaces to realize the separation of EEG signals and EOG artifacts. Let S = [S1, S2, …, SM] be the M mutually statistically independent source signal. X = [x1, x2, …, xN] is the n-dimensional observation signal generated by the linear mixing of S through an unknown matrix A, i.e., X = A × S. Under the condition that both A and S are unknown, the ICA method uses the assumption that X and S are statistically independent to find a linear transformation separation matrix W to make the output signal approach S as much as possible. The FastICA (Chen et al., 2018) algorithm takes the maximum negative entropy as the search direction and can sequentially extract independent sources. In addition, it adopts fixed-point iteration to make the convergence faster and more robust. Therefore, we use this method to process the IMF and obtain the source signal S.

The EEG signal comes from brain bioelectric activity and contains much physiological and pathological information, while the EOG signal only represents eye movement and blinking. Compared with EOG signals, EEG signals have more complex characteristics. The higher the complexity is, the higher the corresponding entropy. Therefore, the components with high entropy can be extracted as EEG signals. Compared with the approximate entropy (Pincus and Goldberger, 1994; Li et al., 2010), the sample entropy has a better estimation effect in the time-domain statistics aspect (e.g., mean and variance) and can be used to calculate the mixed signal composed of the determined signal and random signal. Therefore, we use Algorithm 2 (Calculate sample-entropy) to remove the EOG artifacts in this paper. The Algorithm 2 is as follows:

ALGORITHM 2
www.frontiersin.org

Algorithm 2. Calculate sample-entropy sampEn(m, r, N)

DWT-CEEMDAN-ICA algorithm

The innovation of this algorithm is that it can effectively solve the overcomplete and modal aliasing problem existing in the current single-channel EEG removal algorithm, and improve the retention of effective EEG information. The algorithm proposed in this paper is described in detail as follows:

1: The collected EEG signal is transformed by db4 wavelet to obtain a low-frequency approximate component A7and seven high-frequency detail components Di(i = 1, …, 7) corresponding to different layers decomposed by 7 layers.

2: These components are single-branch reconstructed and CEEMDAN decomposed to obtain several IMFs, and FastICA decomposition is carried out to calculate several independent components and their sample entropy.

3: The component corresponding to the sample entropy satisfying the threshold discriminant proposed by Gomez-Herrero et al. (2006) is regarded as an EOG artifact and set to zero. Then, the inverse ICA transform is carried out to obtain the EEG signal after the artifact is removed. The threshold discriminant proposed by Gomez Herrero is as follows:

φ(k+1)-φ(k)<φ(k)-φ(k-1)S.t. 1<k  [N/2]    (17)

where φ(k) represents the entropy value of the k-th independent component sorted in ascending order. The ICA components corresponding to the top k entropy values are regarded as EOG artifacts where k takes the smallest integer which satisfies the above conditions.

4: Repeat steps 2–3 for the remaining wavelet coefficients to obtain all the signals after removing artifacts and then carry out wavelet reconstruction to obtain a complete EEG signal without artifacts.

Figure 1 shows the flowchart of the proposed algorithm.

FIGURE 1
www.frontiersin.org

Figure 1. Flowchart of DWT-CEEMDAN-ICA.

Experiments

In this section, a series of experiments are conducted to evaluate the efficiency of the proposed algorithm compared with a number of state-of-the-art algorithms used to remove EOG artifacts. The algorithms in our comparisons include the following:

1) Wave transform (WT) algorithm. First, the original signals are decomposed by wavelet transform, then the wavelet coefficients are denoised by the soft threshold method, and finally the signal is reconstructed by wavelet reconstruction.

2) WT-ICA algorithm. First, the original signals are decomposed by wavelet transform, then all wavelet coefficients are decomposed into independent components by the ICA algorithm, and finally the signals are reconstructed by ICA inverse transform of these components.

3) EMD-ICA algorithm. The original signals are first decomposed by empirical mode decomposition to acquire a series of IMFs, and then they are decomposed into independent components. Finally, the signals are reconstructed by ICA inverse transformation of these components.

The original signals used in the experiments are collected from the FP1 channel by a single electrode EEG machine (Mindwave) of NeuroSky Technology Company (NeuroSky). The device runs TGAM EEG module and adopts Bluetooth plus BLE dual-mode transmission. The sampling rate is 512 Hz, and the sampling time is 2 min. First, the ears and forehead of the subject were treated to remove grease and cutin to reduce the components of EMG artifacts. Then, the subject was allowed to close his eyes and rest for 2 min. During the collection, the subject was instructed to remain calm and blink several times in the natural manner. As a result, the collected signals contain fewer other artifacts, which can be approximately considered to contain only EEG signals and EOG artifacts.

The root mean square error RMSE and correlation coefficient R are introduced to evaluate the performance of the algorithms. The smaller the RMSE is, the closer the signal is to the original signals after artifact removal. The better the removal effect is, the larger the correlation coefficient R and the more complete the effective information of the signals that is retained.

RMSE=1ni=1n[xi-yi]2    (18)
R=i=1n(xi-x¯)(yi-y¯)i=1n(xi-x¯)2i=1n(yi-y¯)2    (19)

Preprocessing

Signal jump and mechanical noise are generated in the collection process, resulting in a certain amount of signals that are logically unreasonable and exceed the normal range of the original EEG signals. These signals affect the analysis quality of EEG signal data. Therefore, it is necessary to check the consistency and remove the unreasonable information from the original signals and estimate and supplement the defective signal by using the mean difference method of surrounding signals to restore continuity and time-frequency characteristics. The original signals and the signals after the consistency check and interpolation are shown in Figures 2, 3, respectively.

FIGURE 2
www.frontiersin.org

Figure 2. Original signal.

FIGURE 3
www.frontiersin.org

Figure 3. Signal after consistency check and interpolation.

The frequency used to analyze the characteristics of EEG signals is mainly below 64 Hz, so it is necessary to filter and eliminate frequency interference by using a 0.05–64 Hz bandpass filter and a 50 Hz band-stop filter, respectively. Figures 4, 5 are the spectrum diagrams after bandpass and band-stop filtering, respectively.

FIGURE 4
www.frontiersin.org

Figure 4. Spectrogram after bandpass filtering.

FIGURE 5
www.frontiersin.org

Figure 5. Spectrogram after bandstop filtering.

Experimental results and analysis

It can be seen from the preprocessed signals that the artifacts are mainly concentrated on the sampling points in the range of 30,000–40,000. Therefore, the points in this range were used to compare the effectiveness of removing EOG artifacts. The retention degree for effective information is obtained by comparing the correlation coefficients of the sampling points in the range of 50,000–55,000 before and after being processed by the DWT-CEEMDAN-ICA algorithm.

Figure 6 shows the results of the preprocessed signals after wavelet decomposition. Then, CEEMDAN decomposition is carried out to obtain IMFi (i = 1.0.0.15) components for each wavelet coefficient. Take the high-frequency detail coefficient D5 as an example. Figure 7 shows the decomposition results. All the IMF components are decomposed by the ICA algorithm to obtain each independent component, as shown in Figure 8. Finally, the sample entropy of each ICA component is calculated, as shown in Table 1. From the table, we can see that the more complex the independent component is, the higher the sample entropy. Therefore, the component corresponding to the sample entropy satisfying the threshold discriminant is regarded as an EOG artifact and set to zero. Figure 9 shows the comparison result of D5 between its original value and the value after filtering the EOG artifacts by the DWT-CEEMDAN-ICA algorithm. We can see from the figure that our proposed algorithm has a very good effect in removing EOG artifacts.

FIGURE 6
www.frontiersin.org

Figure 6. Wavelet decomposition results.

FIGURE 7
www.frontiersin.org

Figure 7. CEEMDAN decomposition results on high-frequency detail coefficient D5.

FIGURE 8
www.frontiersin.org

Figure 8. ICA decomposition results.

TABLE 1
www.frontiersin.org

Table 1. Sample entroy.

FIGURE 9
www.frontiersin.org

Figure 9. D5 signal and clean D5 signal.

However, there is still a gap between the original signals and the signals reconstructed only from the D5 coefficient. To prevent distortion and retain more effective EEG signals, other wavelet coefficients also need to be processed by the DWT-CEEMDAN-ICA algorithm and then reconstructed by wavelet reconstruction to obtain the final clean EEG signals without EOG artifacts. The result is shown in Figure 10. From this figure, it is shown that the algorithm removed EOG artifacts well and had a high degree of fitting for the EEG signals.

FIGURE 10
www.frontiersin.org

Figure 10. Original signal and clean signal.

Figures 1113 show the results of the other three algorithms on EOG artifact removal. We can see that the WT algorithm can not only effectively remove EOG artifacts but can also remove many valuable EEG signals, resulting in serious signal distortion. Moreover, this algorithm is highly subjective and the threshold and basis function need to be manually selected. Compared with the WT algorithm, the WT-ICA algorithm can retain more original EEG signals. Because of the problem of overcompleteness and the subjectivity of judgment, the effect is different every time, and the fitting of the original EEG signals is not good. The EMD-ICA algorithm cannot effectively remove EOG artifacts because of mode aliasing and noise.

FIGURE 11
www.frontiersin.org

Figure 11. Raw EEG signal and result of WT.

FIGURE 12
www.frontiersin.org

Figure 12. Raw EEG signal and result of WT-ICA.

FIGURE 13
www.frontiersin.org

Figure 13. Original signal and EMD-ICA algorithm signal.

Figure 14 shows the correlation coefficient R and root mean square error RMSE calculated by using the four algorithms for sampling points in the range of 50,000–55,000. The algorithm proposed in this paper obtains the largest correlation coefficient and the smallest root mean square error, which proves that it not only solves the overcomplete and modal aliasing problems but also effectively removes the EOG artifacts and retains more valuable original information.

FIGURE 14
www.frontiersin.org

Figure 14. Comparison results between the errors and correlation coefficients of each algorithm.

Conclusion

As a result of the problems of overcompleteness and mode aliasing in the single-channel EOG signal removal algorithm, the single-channel EEG equipment is restricted by few acquisition channels and lack of a reference electrode, so EOG artifacts cannot be effectively removed. We propose a novel method in this paper that integrates the discrete wavelet transform, complete empirical mode decomposition for adaptive noise, independent component analysis and the sample entropy algorithm. We carry out a series of experiments to demonstrate its effectiveness. Compared with some existing methods, our method can effectively identify and remove EOG artifacts from original signals while solving the above problems.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by School of Computer Science and Engineering, Guilin University of Aerospace Technology. The patients/participants provided their written informed consent to participate in this study.

Author contributions

QH, ML, and YL: contributed to the design of this work. QH: methodology, validation, formal analysis, and writing—original draft. ML: software, data curation, and data analysis. YL: conceptualization, resources, writing—review and editing, and funding acquisition. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (81960324), General Project of Guangxi Natural Science Foundation (Guangdong-Guangxi Joint Fund, 2021GXNSFAA075037), and Guangxi University Middle-aged and Young Teachers' Basic Scientific Research Ability Improvement Project (2019KY0798, 2022KY0791, 2020KY21023, and 2019KY0822).

Acknowledgments

The authors would like to thank all the subjects for volunteering to participate in our experiments.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Boda, S., Mahadevappa, M., and Dutta, P. K. (2021). A hybrid method for removal of power line interference and baseline wander in ECG signals using EMD and EWT. Biomed. Signal Process. Control. 67, 102466. doi: 10.1016/j.bspc.2021.102466

CrossRef Full Text | Google Scholar

Chen, M. Q., Zhang, X., Chen, X., and Zhou, P. (2018). Automatic implementation of progressive FastICA peel-off for high density surface EMG decomposition. IEEE Trans. Neural Syst. Rehabil. Eng. 26, 144–152. doi: 10.1109/TNSRE.2017.2759664

PubMed Abstract | CrossRef Full Text | Google Scholar

Dora, C., and Biswal, P. K. (2020). An improved algorithm for efficient ocular artifact suppression from frontal eeg electrodes using vmd. Biocybern. Biomed. Eng. 40, 148–161. doi: 10.1016/j.bbe.2019.03.002

CrossRef Full Text | Google Scholar

Gomez-Herrero, G., De Clercq, W, Anwar, H., Kara, O., Egiazarian, K., Van, H. S., et al. (2006). “Automatic removal of ocular artifacts in the EEG without an EOG reference channel,” in Proceedings of the 7thNordic Signal Processing Symposium (NORSIG 2006) (Reykjavik: IEEE), 130–133. doi: 10.1109/NORSIG.2006.275210

CrossRef Full Text | Google Scholar

Gu, Y., Li, X., and Chen, S. (2021). Aoar: an automatic ocular artifact removal approach for multi-channel electroencephalogram data based on non-negative matrix factorization and empirical mode decomposition. J. Neural. Eng. 18, 056012. doi: 10.1088/1741-2552/abede0

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, L., Liu, J., Zhou, D., and Zhou, Q. (2020). Predicting the evolution of hot topics: a solution based on the online opinion dynamics model in social network. IEEE Trans. Syst. Man Cybern. Syst. 50, 3828–3840. doi: 10.1109/TSMC.2018.2876235

CrossRef Full Text | Google Scholar

Jiang, X., Bian, G. B., and Tian, Z. (2019). Removal of artifacts from eeg signals: a review. Sensors. 19, 987. doi: 10.3390/s19050987

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, P. S., Arumuganathan, R., Sivakumar, K., and Vimal, C. (2008). Removal of ocular artifacts in the eeg through wavelet transform without using an eog reference channel. Int. J. Open Probl. Comput. Math. 1, 188–200. doi: 10.1109/JBHI.2021.3131186

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwon, S., Kim, H., and Yeo, W. H. (2021). Recent advances in wearable sensors and portable electronics for sleep monitoring. iScience. 24, 102461. doi: 10.1016/j.isci.2021.102461

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Chen, W. T., Shao, X. M., and Wang, Z. Z. (2010). Analysis of amplitude-integrated EEG in the newborn based on approximate entropy. IEEE Transac. Biomed. Eng. 57, 2459–2466. doi: 10.1109/TBME.2010.2055863

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M. G., Cui, Y., and Yang, J. F. (2013). Automatic removal of ocular artifact from EEG with DWT and ICA method. Appl. Math. Inform. Sci. 2, 809–816. doi: 10.12785/amis/070252

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y. Z., Habibnezhad, M., Shayesteh, S., Jebelli, H., and Lee, S. H. (2021). Paving the way for future EEG studies in construction:dependent component analysis for automatic ocular artifact removal from brainwave signals. J. Construct. Eng. Manag. 147, 04021087. doi: 10.1061/(ASCE)CO.1943-7862.0002097

CrossRef Full Text | Google Scholar

Makeig, S., Bell, A. J., Jung, T. P., and Sejnowski, T. J. (1996). Independent component analysis of electroencephalographic data. Advances in Neural Information Processing Syatems. 8, 145–151

PubMed Abstract | Google Scholar

Mammone, N., La, F. F., and Morabito, F. C. (2012). Auto matic artifact rejection from multichannel scalp EEG by wavelet ICA. IEEE Sens. J. 12, 533–542. doi: 10.1109/JSEN.2011.2115236

CrossRef Full Text | Google Scholar

Meng, L., Chen, Q., Jiang, X, Y., Liu, X, Y., Fan, J, H., Dai, C. Y., et al. (2022). Evaluation of decomposition parameters for high-density surface electromyogram using fast independent component analysis algorithm. Biomed. Signal Process. Control. 75, 103615. doi: 10.1016/j.bspc.2022.103615

CrossRef Full Text | Google Scholar

Miao, M., Hu, W., Xu, B., and Zhang, J. (2021). Automated CCA-MWF algorithm for unsupervised identification and removal of EOG artifacts from EEG. IEEE J. Biomed. Health Inform. 26, 3607–3617. Available online at: https://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.502.6932&rep=rep1&type=pdf

PubMed Abstract | Google Scholar

Mourad, N. (2022). Ecg denoising based on successive local filtering. Biomed. Sig. Process. Control. 73, 103431. doi: 10.1016/j.bspc.2021.103431

CrossRef Full Text | Google Scholar

Park, C., Looney, D., Kidmose, P., Ungstrup, M., and Mandic, D. P. (2011). Time-frequency analysis of eeg asymmetry using bivariate empirical mode decomposition. IEEE Transac. Neural Syst. Rehabil. Eng. 19, 366–373. doi: 10.1109/TNSRE.2011.2116805

PubMed Abstract | CrossRef Full Text | Google Scholar

Pincus, S. M., and Goldberger, A. L. (1994). Physiological time-series analysis: what does regularity quantify? Am. J. Physiol. 266, 1643–1656. doi: 10.1152/ajpheart.1994.266.4.H1643

PubMed Abstract | CrossRef Full Text | Google Scholar

Saeidi, M., Karwowski, W., Farahani, F. V., Fiok, K., Taiar, R., Hancock, P. A., et al. (2021). Neural decoding of EEG signals with machine learning: a systematic review. Brain Sci. 11, 1525. doi: 10.3390/brainsci11111525

PubMed Abstract | CrossRef Full Text | Google Scholar

Semlitsch, H. V., Anderer, P., Schuster, P., and Presslich, O. (1986). A solution for reliable and valid reduction of ocular artifacts, applied to the P300 ERP. Psychophysiology. 23, 695–703. doi: 10.1111/j.1469-8986.1986.tb00696.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, R., Chan, C., Hsiao, J. H., and Tang, A. C. (2021). Validation of sobi-dans method for automatic identification of horizontal and vertical eye movement components from eeg. Psychophysiology. 58, e13731. doi: 10.1111/psyp.13731

PubMed Abstract | CrossRef Full Text | Google Scholar

Vigon, L., Saatchi, M. R., Mayhew, J., and Fernandes, R. (2000). Quantitative evaluation of techniques for ocular artefact filtering of EEG waveforms. IEE Proc. Sci. Measure. Technol. 147, 219–228. doi: 10.1049/ip-smt:20000475

CrossRef Full Text | Google Scholar

Wan, C. T., Chen, D. Y., Huang, Z. Q., and Luo, X. (2022). A wearable head mounted display bio-signals pad system for emotion recognition. Sensors. 22, 142. doi: 10.3390/s22010142

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Z., and Huang, N. E. (2011). Ensemble empirical mode decomposition: a noise-assisted data analysis method. Adv. Adapt. Data Anal. 1, 1–41. doi: 10.1142/S1793536909000047

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y., Luo, M. Z., Li, T., and Song, G. B. (2018). ECG signal de-noising and baseline wander correction based on CEEMDAN and wavelet threshold. Sensors. 17, 2754. doi: 10.3390/s17122754

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, P., Jiang, L., Liu, J., and Zhou, D. (2020). “Dual-level attention based on heterogeneous graph convolution network for aspect-based sentiment classification,” in Wireless Communications and Mobile Computing. doi: 10.1109/SmartCloud49737.2020.00022

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: electroencephalogram, discrete wavelet transform, empirical mode decomposition, independent component analysis, sample entropy

Citation: Hu Q, Li M and Li Y (2022) Single-channel EEG signal extraction based on DWT, CEEMDAN, and ICA method. Front. Hum. Neurosci. 16:1010760. doi: 10.3389/fnhum.2022.1010760

Received: 03 August 2022; Accepted: 29 August 2022;
Published: 21 September 2022.

Edited by:

Tianyou Yu, South China University of Technology, China

Reviewed by:

Jiang Lei, Hunan University of Science and Technology, China
Wang Zelin, Nantong University, China

Copyright © 2022 Hu, Li and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yunde Li, lyde2008@163.com

Download