Comparison of signal decomposition techniques for analysis of human cortical signals

Abstract

Abstract Objective. Conventional neural signal analysis methods assume that features of interest are linear, time-invariant signals confined to well-delineated spectral bands. However, new evidence suggests that neural signals exhibit important non-stationary characteristics with ill-defined spectral distributions. These features pose a need for signal processing algorithms that can characterize temporal and spectral features of non-linear time series. This study compares the effectiveness of four algorithms in extracting neural information for use in decoding cortical signals: Fast Fourier Transform bandpass filtering (FFT), principal spectral component analysis (PSCA), wavelet analysis (WA), and empirical mode decomposition (EMD). Approach. Electrocorticographic signals were recorded from the motor and sensory cortex of two epileptic patients performing finger movements. Each signal processing algorithm was used to extract beta (10–30 Hz) and gamma (66–114 Hz) band power to detect thumb movement and decode finger flexions, respectively. Naïve-Bayes (NB), support vector machine (SVM), and linear discriminant analysis (LDA) classifiers using each signal were validated using leave-one-out cross-validation. Main results. Decoders using all four signal decompositions achieved above 90% average accuracy in finger movement detection using beta power. When decoding individual finger flexion using gamma, the PSCA NB classifiers achieved 78 ± 4% accuracy while FFT, WA, and EMD analysis achieved accuracies of 73 ± 8%, 68 ± 7%, and 62 ± 3% respectively, with similar results using SVM and LDA. Significance. These results illustrate the relative levels of useful information contributed by each decomposition method in the case of finger movement decoding, which can inform the development of effective neural decoding pipelines. Further analyses could compare performance using more specific non-sinusoidal features, such as transients and phase-amplitude coupling.