Neural decoding of electrocorticographic signals using dynamic mode decomposition

Objective. Brain-computer interfaces (BCIs) using electrocorticographic (ECoG) signals have been developed to restore the communication function of severely paralyzed patients. However, the limited amount of information derived from ECoG signals hinders their clinical applications. We aimed to develop a method to decode ECoG signals using spatiotemporal patterns characterizing movement types to increase the amount of information gained from these signals. Approach. Previous studies have demonstrated that motor information could be decoded using powers of specific frequency bands of the ECoG signals estimated by fast Fourier transform (FFT) or wavelet analysis. However, because FFT is evaluated for each channel, the temporal and spatial patterns among channels are difficult to evaluate. Here, we used dynamic mode decomposition (DMD) to evaluate the spatiotemporal pattern of ECoG signals and evaluated the accuracy of motor decoding with the DMD modes. We used ECoG signals during three types of hand movements, which were recorded from 11 patients implanted with subdural electrodes. From the signals at the time of the movements, the modes and powers were evaluated by DMD and FFT and were decoded using support vector machine. We used the Grassmann kernel to evaluate the distance between modes estimated by DMD (DMD mode). In addition, we decoded the DMD modes, in which the phase components were shuffled, to compare the classification accuracy. Main results. The decoding accuracy using DMD modes was significantly better than that using FFT powers. The accuracy significantly decreased when the phases of the DMD mode were shuffled. Among the frequency bands, the DMD mode at approximately 100 Hz demonstrated the highest classification accuracy. Significance. DMD successfully captured the spatiotemporal patterns characterizing the movement types and contributed to improving the decoding accuracy. This method can be applied to improve BCIs to help severely paralyzed patients communicate.


Introduction
A brain-computer interface (BCI), which enables communication via brain activity, is the only way to express thoughts for severely paralyzed patients who are completely unable to move their body [1]. BCIs are in high demand for these patients [2], who choose to undergo invasive procedures involving the placement of electrodes on the surface [3] or inside the cortex [4,5]. An electrocorticographic (ECoG) signal, which is measured by a surface electrode, is a promising signal for use in BCIs [6] due to its stability [7] and less invasive nature compared to penetrating electrodes [8,9]. An ECoG-BCI was successfully applied to a severely paralyzed patient with amyotrophic lateral sclerosis (ALS) and enabled the patient to stably communicate for months [10]. However, much less information is derived from an ECoG signal than from penetrating electrodes [11]. The patient with ALS was able to type faster using unit activity recorded by a penetrating electrode, at more than 30 characters per minute [11], than using ECoG signals [10]. Thus, for the clinical application of ECoG-BCIs, some novel technologies to extract information from ECoG signals are necessary.
To decode the ECoG signals, the oscillation component is considered to have important information, and various power spectral density (PSD) estimations [12][13][14] have been used as features for machine learning [15][16][17][18][19]. Among the methods to evaluate the PSD [18], fast Fourier transform (FFT) is often used for online decoding because of the low calculation cost [20]. FFT evaluates the powers of some frequency bands for each channel in a time series and for individual time points in space. However, because these features are independently evaluated, FFT cannot evaluate the patterns of ECoG signals in both space and time simultaneously. Therefore, a method to characterize the nonlinear dynamics as the linear summation of coherent patterns across both space and time (e.g. frequency and phase) is necessary.
Dynamic mode decomposition (DMD) has recently attracted attention as a way of obtaining modal representations of nonlinear dynamics from general multivariate time series data without explicit prior knowledge about the dynamics [21,22]. DMD is closely related to spectral analysis of the Koopman operator (see supplementary information (available online at stacks.doi.org/JNE/17/036009/mmedia)) and has been successfully applied for the extraction of spatiotemporal patterns, broadly in fluid mechanics [21][22][23][24], but rarely in neuroscience, e.g. ECoG and functional magnetic resonance imaging (fMRI) [25][26][27][28]. Although a previous study using human ECoG signals demonstrated that the DMD was able to capture some characteristic spatiotemporal patterns related to spindles [27], the spatiotemporal patterns were not used for neural decoding. It has not been revealed whether the spatial and temporal coherent patterns identified by DMD improve the accuracy of the neural decoding used in BCI technology.
A problem in the application of DMD for neural decoding is that DMD transforms the spatial and temporal signals into a set of vectors representing the characteristic spatial mode, which has a characteristic frequency in the time domain. Notably, because the characteristic frequency of the DMD mode is estimated from single trial data, the frequency is not fixed across trials and varies depending on the original signals. Therefore, to classify the signals by the vectors estimated by DMD, proper metrics are needed to evaluate how the sets of vectors are similar among trials. Here, we used the principal angle as a metric to evaluate the relationship among the sets of vectors for each trial [29]. A Grassmann kernel was used with support vector machine (SVM) to classify the set of vectors. In this study, we classified upper limb movements using the DMD components with SVM. The classification accuracy using the DMD component was compared with that using the FFT powers and that using the DMD components whose phases were shuffled.

Subjects
11 patients (four females, seven males; age range, 13-66 years) with subdural electrodes participated in this study. Eight patients were the same patients from our previous report [17]. Five patients had different degrees of motor dysfunction and sensory disturbances in their upper limbs due to strokes without damage to the sensorimotor cortex, previous surgery on the motor cortex, or brain tumor on the sensorimotor cortex (table 1). Six patients with epilepsy had no motor dysfunctions or sensory disturbances. All participants or their guardians gave written informed consent to participate in the study, which was approved by the ethics committee of Osaka University Hospital (No. 08061, No. 19257).

Movement tasks
Experiments were performed approximately 1 week after electrode placement. The patients performed three types of upper limb movements, selecting among grasping, pinching, hand opening, thumb flexion, and elbow flexion [17]. We presented the patients with three sound cues at 1 s intervals. The patients performed one of three prespecified movements after the third sound cue. The patient selected the movement type for each trial, so the order and number of exercises were random. The numbers of the selected movement types are shown in table 2. The patients were instructed to perform the selected movement once, immediately after the third cue, and to return their hands to the resting position, which was an immediate hand posture among the performed movements. For the resting position, patients were instructed to relax their hands or elbows with slightly flexed joints. Each type of movement was performed approximately 30 to 100 times. The movement instructions were delivered using a PC monitor controlled by ViSaGe (Cambridge Research System, Rochester, UK) placed in front of the patients.

ECoG recording and preprocessing
ECoG signals were recorded at 1 kHz by an EEG-1200 system (Nihon Koden, Tokyo, Japan). The subdural  Grasp, thumb, elbow 71 69 39 32 electrodes were placed on the sensorimotor cortex based on clinical necessity. For each patient, 15-60 electrodes were implanted. All signals were subtracted by the average of the signals of all electrodes at each time point, a well-known method for the calculation of a common average reference [30,31]. An epileptologist identified the channels with severe noise, and these were excluded from the analysis. Three channels for patient 7 and one channel for patient 10 were removed (supplementary figures 3, 5). For the following analysis, we used the ECoG signals at 0.5 s after the third sound cue, S execute , and at 1.0 s before the first sound cue, S norm (figure 1). Here, S execute and S norm were composed of n channels for 500 and 1000 time points, respectively. Here, we used 1.0 s of ECoG signals for normalization (S norm ) to decrease the variability in S execute . Notably, we also used 0.5 s of ECoG signals before the first sound cue for the S norm , but the classification accuracy using the FFT feature did not significantly change with the normalization of 0.5 s (supplementary information, table 1). All processes were implemented by MATLAB (R2015b).

FFT features for neural decoding
As a conventional approach for the neural decoding of ECoG signals, FFT was used to create the FFT features of the signals [32]. FFT was applied as 512 points and 1024 points for S execute (500 points) and S norm (1000 points), respectively, by zero-padding with a hamming window. S norm was only used to normalize S execute to reduce variation among trials. The power features of four frequency bands were obtained by averaging the estimated power of the four frequency bands: theta (4-8 Hz), alpha (8-13 Hz), beta (13)(14)(15)(16)(17)(18)(19)(20)(21)(22)(23)(24)(25)(26)(27)(28)(29)(30), and high gamma (80-150 Hz). In addition to these powers, we also calculated the average 0 Hz power for S execute and S norm . Finally, these five features of S execute were normalized by the mean power p norm and standard deviations σ norm of the corresponding features derived from S norm to create the FFT features, which consisted of five features by n channels [17,33]. The features were obtained by p execute −p norm σnorm (p execute is the power feature vector of S execute ). Here, the FFT powers of S execute were normalized because normalized features were shown to be decoded with better accuracy than FFT powers of S execute without normalization [15].

Computation of the dynamic mode decomposition.
Assuming that the sensorimotor cortex during exercise is one dynamic system, the system can be described using its observed variable x (ECoG signals) as follows: Figure 1. ECoG signals during hand movements. After selecting one of the three movement types, the patient performed the selected movements once after the third sound cue (shown as black speaker). The ECoG signals during the 0.5 s before the first sound cue and those after the last sound cue were used for the analysis.
where t and µ are the time and system parameters. Because the actual recorded signals are values at discrete times, x k = x (k∆t) is delimited by ∆t, and the dynamic system is described as Although F is difficult to obtain analytically, we can obtain linearly approximated dynamics x k+1 = Ax k from the recorded signals. The DMD can be used to estimate A of unknown dynamics. First, we briefly explain how the DMD was evaluated for column vector x k , which is the observations across all electrodes at the k-th time sample [27]. Gathering measurements from m snapshots in time, we constructed two n × (m − 1) raw data matrices, X 1 and X 2 : Here, the DMD of the data matrix pair X 1 and X 2 was given by the eigendecomposition of the transition matrix A, which is an n × n matrix satisfying X 2 = AX 1 .
Using the singular value decomposition (SVD) of the data matrix X 1 However, the eigenvalue decomposition of A requires a tremendous amount of calculation. Therefore, we used the orthogonal matrix U to was set to obtain the eigenvector W by applying the eigenvalue decomposition AW = W Λ. Then, by putting W ≜ −1 W, the DMD mode was obtained as The absolute values and phases of Λ are the decay rates and frequencies, respectively.
Finally, we were able to write an approximation of the observed data as a simple dynamic model X (t), where Ω = log (Λ) /∆t, t is time, and z is a set of weights to match the first time point measured such that X (1) = Φz. It should be noted that Φ and Ω are sets of vectors with complex values (figure 2).
DMD was originally used to evaluate the X 1 of n ≫ m. However, in our data, the ECoG signal has approximately 100 channels for 500 ms. Therefore, to calculate the DMD, we constructed a new augmented data matrix X 1 aug as follows: Figure 2. DMD of ECoG signals during movements. The ECoG signals of n channels (Sexecute) were converted into the n × 300 matrix of magnitude and phases that represents the complex values of the DMD modes Φj. Λ real and Λ imagin represent the decay rates and frequencies.
Here, h is a minimum integer that satisfies h < m+1 n+1 . X 1 aug is an array generated from X 1 of m length by stacking them in the channel dimension by sliding the step one by one. We constructed X 2 aug using the same procedure. The DMD features were obtained from X 1 aug and X 2 aug by the method previously mentioned. Notably, we generated X 1 aug and X 2 aug using only S execute without any normalization by S norm to avoid damaging the spatial relations among the channels of S execute by subtracting the averages and dividing the standard deviations of S norm for each channel. In fact, the DMD mode Φ will be changed by the normalization of X 1 aug and X 2 aug by the same S norm .
Here, the number of channels, n, and the size of the X 1 aug were different among the subjects. Therefore, to equalize the number of DMD modes, we selected the 300 modes with the largest singular values. Because the DMD mode φ i aug yielded from X 1 aug and X 2 aug was a vector of hn × 1, we used the first n vector of φ i aug corresponding to φ i , which is the mode for the original data X 1 and X 2 for the decoding task. In addition, because the size of φ i generally decreases with the higher frequency range (imagery part of the eigenvalue λ i ), we normalized each mode by the L2 norm of the mode φ i . Notably, the principal angle requires the orthonormal matrix to be equalized among the permutation of the composed vectors. To simplify the calculation, we aligned the mode magnitude by the L2 norm instead of the orthogonalization.

Neural decoding by the Grassmann kernel.
The FFT feature is represented as a set of vectors of n × 5. The DMD features are composed of eigenvectors Φ j corresponding to a characteristic frequency ω j (j = 1 . . . 300).
Therefore, the DMD feature is represented as a set of vectors n × 300. Notably, the DMD feature for classification is a subspace of the representation. To compare the subspaces for the classification task, we evaluated a principal angle, which measures the distance between the subspaces based on the set of vectors. Here, based on the principal angle, we used the Grassmann kernel to classify the set of vectors by SVM [29]. Y 1 and Y 2 were set as two orthonormal matrices of size D by m. The principal angles 0 ≦ θ 1 ≦…≦ θ m ≦ π/2 between two subspaces, span (Y 1 ) and span (Y 2 ), were defined recursively by the following: The principal angles were computed from the SVD of Y 1 ′ Y 2 as follows: , and cos Θ is the diagonal matrix cos Θ = diag(cos θ 1 . . . cos θ m ).
Based on the principal angles, the projection kernel of subspaces Y 1 and Y 2 was determined as follows: Here, we applied the projection kernel to the DMD features Φ i and Φ j to generate the Gram matrix k p Φ i , Φ j . To apply the Gram matrix to SVM, each component of the Gram matrix was divided by the average of all components of the Gram matrix. For comparison, the FFT features were also classified using SVM with the projection kernel. We also calculated the inner product matrix correlation [34] defined as r in (A, B) [tr(A ′ A)tr(B ′ B)] 1/2 of two i × j matrices A and B to examine the independent information contained in the DMD and FFT features.

T-distributed stochastic neighbor embedding analysis of the Gram matrix.
Based on the Grassmann kernel, the distance between each DMD mode is represented in the Grassmann space, in which the distance between each mode of the Gram matrix G i, j corresponding to the inner product of Φ i and Φ j was defined as To visualize the distribution of each DMD mode in the space, we applied tdistributed stochastic neighbor embedding (t-SNE) to the space [35].

Decoding analysis
The classifier was independently trained for each patient by SVM. The classification accuracy was evaluated by 10-fold nested cross validation. In the nested loop, the cost parameter was optimized in the range of 10 −1 > cost > 10 8 . The classification accuracies were evaluated as balanced accuracy because the number of each movement type was uneven [36]. The balanced accuracy is defined as the average of recall obtained in each category.
We evaluated how the classification accuracies depend on the characteristic frequency. We selected 10 modes Φ j (j = 1 . . . 10), which have the characteristic frequency ω j closest to 29 selected frequencies: 10, 20, …, 290 Hz. Then, the classification was performed with 10 selected modes for each frequency.
Φ j is the vector of complex values. Each complex value is expressed as a radius and a phase in a complex plane. To evaluate how much the phase information contributes to the decoding, we also evaluated the classification accuracy using the DMD features in which the phases of each Φ j were randomized.

Statistical tests
We compared the classification accuracies evaluated by the two types of features, DMD features and FFT features, by Student's paired t-test, with p < 0.05 indicating significance.
We also compared the correctness of prediction by SVM trained with DMD and FFT features by Welch's t-test, with p < 0.05 indicating significance.
The classification accuracy of each frequency was compared to the classification accuracy using the same features with shuffled labels by using Bonferroni-corrected Student's paired t-test, with p < 0.05 indicating significance.
The phases of the DMD mode of each electrode were evaluated by the Rayleigh test to show some biases in the phase distribution of the electrode among each movement task.

Gram matrix based on projection kernel with DMD mode and clustering of dynamics
For illustrative patient 1, we calculated the projection kernel values G i, j , which correspond to the inner products between each DMD mode Φ i and Φ j . The projection kernel values demonstrated some clusters according to the movement types ( figure 3(A)). Based on the projection kernel, we further calculated the distance matrix between each DMD mode, which represents the Grassmann space for the SVM classification. To visualize the Grassmann space, t-SNE was applied for the distance matrix to embed each trial into a two-dimensional space ( figure 3(B)). Each movement type was classified into different clusters. The Grassmann space using the projection kernel seemed to be successful at classifying the movement types.

Neural decoding using DMD features
Using the Grassmann kernel, the classification of the movement types was evaluated. The classification accuracy significantly varied among the different features: DMD feature, FFT feature, DMD feature with shuffled phase, and DMD features with shuffled labels (p = 9.12 × 10 -11 , one-way ANOVA; figure 4). The three types of hand movements were successfully classified using the DMD features (79.8 ± 6.28%; mean ± 95% confidence interval), the FFT features (67.3 ± 8.56%), and the DMD features with shuffled phase (66.7 ± 9.96%), with accuracies exceeding the accuracy using the DMD features with shuffled labels (33.1 ± 3.90%, i.e. chance level). Among these classification schemes, the classification accuracy using the DMD features was significantly superior to that using the FFT features and the DMD features with shuffled phase (p < 0.05, paired Student's t-test, uncorrected).
We also examined how the prediction of the DMD and FFT features were correlated. First, we evaluated the inner product matrix correlation of the Gram matrices of the DMD features and the FFT features. For 10 out of 11 patients, the correlation was more than 0.8 (supplementary information, table 2), suggesting that the Gram matrices were similar to each other. However, for patient 5, whose correlation was 0.359, the classification accuracy was higher with the DMD features for this patient. Second, we compared the predictions using two features. We divided the classification results with the two features into three  groups for each patient: Group 1, the predictions of both features were the same; Group 2, the prediction of DMD was correct, but that of FFT was incorrect; and Group 3, the prediction of DMD was incorrect, but that of FFT was correct. Among all patients, the ratio of the groups was significantly high for Group 1 (66.5 ± 8.02%) (figure 5). Moreover, for the conflicting predictions, Group 2 was significantly larger than Group 3. Therefore, these two features mostly result in the same predictions, but in the case of conflicting predictions, the DMD features predicted significantly better than the FFT features.
We compared the classification accuracy using subsets of the DMD features to evaluate which frequency band modes contributed to the classification. The classification accuracies using DMD features of the characteristic frequency close to the selected frequency significantly exceeded the accuracy by chance at 60, 80, 100, 110, and 150 Hz (p < 0.01, Bonferronicorrected paired Student's t-test; figure 6). The peak mean accuracy was 50.4 ± 2.42% (mean ± 95% confidence interval) at 100 Hz, which corresponds to the high gamma band.

DMD modes during movement
The DMD mode captured spatiotemporal patterns characterizing movement type. For the high gamma band frequency (80-150 Hz) in patient 1, the averaged amplitudes and phases of the DMD mode were evaluated and color-coded according to electrode position ( figure 7(A)). The magnitude of the DMD mode was high at the electrodes on the primary motor cortex and demonstrated different patterns in accordance with movement type for patient 1. In addition, the relative phases of each electrode were significantly biased around the motor cortex, and the averaged relative phases were different among movement types. Similarly, the magnitudes and phases of the DMD modes changed based on the movement types for all other patients (supplementary information). Figure 5. Comparison of correctness of prediction by SVM trained with DMD and FFT features. The mean and 95% confidence interval of prediction ratio where the two sets of features had the same prediction label, where DMD had correct prediction labels and FFT had incorrect prediction labels, and where DMD had incorrect prediction labels and FFT had correct prediction labels ( * * p < 0.01, Bonferroni-corrected Welch's t-test). Figure 6. Classification accuracy for each frequency mode. The classification accuracies using the DMD features whose frequency components were close to each frequency point (blue line) were compared to those using the same features with shuffled labels (chance level, orange line). The accuracies at approximately 100 Hz were significantly higher than those at chance level ( * p < 0.05, * * p < 0.01, paired Student's t-test, Bonferroni corrected).
On the other hand, although the FFT powers of each electrode averaged among the same frequency band demonstrated some spatial pattern characterizing the movement type, the distribution of the high powers was different than that of the magnitude of the DMD mode and not restricted around the primary motor cortex ( figure 7(B)). The DMD mode captured some characteristic spatiotemporal patterns of the phases and magnitudes.

Discussion
Our analyses indicated that the DMD mode with projection kernel provided better accuracy of the neural decoding than the FFT power, thereby demonstrating that DMD succeeded in capturing neural information better than FFT. Moreover, the projection kernel was suggested to represent the characteristic differences in ECoG signals among different movement types. Our proposed method using DMD with a projection kernel is a promising technique for neural decoding. DMD succeeded in extracting some spatial and temporal patterns characterizing three types of hand movements. Although a previous study demonstrated that DMD extracts some characteristic spatial patterns in ECoG signals [27], it remained unclear whether these patterns included more information about the movement types than the patterns evaluated by FFT. In our study, the decoding accuracy using DMD features was higher than that using FFT and DMD features with shuffled phases. Notably, both DMD and FFT were able to capture spatial and temporal patterns of each frequency range, as illustrated in figure 4. However, FFT was not able to maintain the information of phases between the dynamics of different channels. Indeed, the classification accuracy using the DMD mode significantly decreased when the phases of the DMD modes were shuffled and were similar to the accuracy using the FFT. Moreover, comparing the predictions of the DMD features and the FFT features, these predictions mostly overlapped, but the DMD features had more useful information than the FFT features. Therefore, our data suggest that the DMD mode improved the classification accuracy compared to the FFT by capturing some spatially and temporally coherent patterns in the ECoG signals to obtain motor information.
Among the frequency features examined, the DMD modes at approximately 100 Hz demonstrated higher classification accuracy, implying that they contained a larger amount of motor information. Some studies using FFT powers have also shown that the frequency powers at approximately 100 Hz provide more motor information [15,37]. Here, DMD was suggested to extract the spatial and temporal pattern relating to the oscillation at 100 Hz. Moreover, our results demonstrated that the DMD mode close to 60 Hz was also decoded with accuracy exceeding the chance level. Although a previous study using FFT signals demonstrated that the powers at approximately 60 Hz were difficult to decode due to some electrical power noise [37], our results suggested that the DMD mode is useful for extracting neural information from the ECoG signals, even with noise.
There are some problems in applying the proposed decoding method to a BCI. First, the DMD mode calculation must be made within at least several hundred milliseconds, which would require sufficient machine power and an appropriate approximation of the calculation. An incremental SVD might be useful to obtain such an appropriate approximation [38][39][40]. Second, we might need to select an appropriate DMD mode for online use during which the activity patterns of the brain continuously change. Because DMD assumes a stationary nonlinear property over a certain period of time, even for nonstationary dynamics, adapting the method to nonstationary data might be difficult. Therefore, some other machine learning techniques, such as a hidden Markov model, may be needed to address nonstationary states [41].
Although some technical improvements (e.g. calculation speed, adapting to nonstationary data) are necessary, improvements in the decoding accuracy of the performed movement types will improve the BCI as a communication device for the patient with limited ability to move their bodies. Our previous study demonstrated that the ECoG signals between the performed movement and attempted movements of paralyzed patients were similar and were decoded similarly. Our proposed method is also expected to improve the performance of attempted movements. Therefore, our proposed method with DMD will improve the decoding accuracy for paralyzed patients. Using precise decoding of attempted movements, paralyzed patients will be able to control some external devices to communicate with others [10]. Our method will improve BCI performance to realize better communication using BCI for paralyzed patients.
Practically, our method can be applied for to the functional mapping of movements or other various cortical functions, such as language, which has been developed with power spectral analysis of ECoG signals [42]. DMD will reveal the spatially and temporally coherent pattern of the ECoG signals in the characteristic frequency range with the highest information regarding particular cortical functions and might improve the accuracy of the functional mapping compared to that using power spectral analysis without phase coherence among electrodes. Similarly, our method can be applied to identify particular characteristic ECoG signals related to epileptic activity. Our neural decoding with DMD will be useful to identify spatially and temporally coherent patterns that represent particular types of information.
This study had some limitations. To simplify and shorten the calculation, we aligned the mode magnitude by the L2 norm instead of the orthogonalization. Therefore, the projection kernels between two matrices were not guaranteed to be the same when the rows of the matrix were exchanged with each other [29]. However, our results demonstrated that the simplified method succeeded in decoding three movement types with accuracy that exceeded that using FFT. The proposed method with a projection kernel was shown to be useful for motor decoding. Our future work should evaluate how orthogonalization affects the decoding accuracy. In addition, our current method ignored the frequency information. Although we used the DMD modes at specific frequencies, we did not use the frequency components for the decoding. The decoding accuracy is expected to improve with the addition of the frequency components into the DMD features. Notably, we compared the classification accuracy of three types of movements that were different among patients. The differences in the movement types might have affected the accuracy. However, by comparing the accuracy using FFT features and the DMD features with Grassmann kernels for each patient, we eliminated the effects of the different sets of movement types. Moreover, the ECoG signals in this study came from two types of patients with and without paralysis in their hands. These patients had different diseases and were implanted with different numbers of electrodes but implantations on their sensorimotor cortices was common among the patients. Notably, those differences induced no apparent differences in the decoding accuracies (supplementary information). However, the anatomical pattern of the DMD features might be affected by ischemic stroke and brain tumors in the motor cortex. We have demonstrated that the classification accuracy was significantly improved by DMD features compared to FFT features using the same classification method for the ECoG signals from the sensorimotor cortex in various conditions.