Decoding of Motor Coordination Imagery Involving the Lower Limbs by the EEG-Based Brain Network

Compared with the efficacy of traditional physical therapy, a new therapy utilizing motor imagery can induce brain plasticity and allows partial recovery of motor ability in patients with hemiplegia after stroke. Here, we proposed an updated paradigm utilizing motor coordination imagery involving the lower limbs (normal gait imagery and hemiplegic gait imagery after stroke) and decoded such imagery via an electroencephalogram- (EEG-) based brain network. Thirty subjects were recruited to collect EEGs during motor coordination imagery involving the lower limbs. Time-domain analysis, power spectrum analysis, time-frequency analysis, brain network analysis, and statistical analysis were used to explore the neural mechanisms of motor coordination imagery involving the lower limbs. Then, EEG-based brain network features were extracted, and a support vector machine was used for decoding. The results showed that the two employed motor coordination imageries mainly activated sensorimotor areas; the frequency band power was mainly concentrated within theta and alpha bands, and brain functional connections mainly occurred in the right forehead. The combination of the network attributes of the EEG-based brain network and the spatial features of the adjacency matrix had good separability for the two kinds of gait imagery (p < 0.05), and the average classification accuracy of the combination feature was 92.96% ± 7.54%. Taken together, our findings suggest that brain network features can be used to identify normal gait imagery and hemiplegic gait imagery after stroke.

For the decoding of motor coordination imagery, conventional features have mainly adopted time-domain features [35], frequency-domain features [36], time-frequency features [37], spatial-domain features [38], and fusion of time-frequency-spatial multidomain features [39]. ese traditional features mainly reflect the activation level of the corresponding brain area. However, motor coordination imagery of the lower limbs involves not only the activation of related brain regions but also the coordinated control among related brain regions [36]. erefore, it is necessary to use brain network features to more comprehensively analyze the neural correlates of motor coordination imagery.
In the present study, we proposed a new paradigm (normal gait imagery and hemiplegic gait imagery after stroke) of motor coordination imagery involving the lower limbs. In our paradigm, normal gait imagery is abbreviated as motor coordination imagery (MCI), and poststroke hemiplegic gait imagery is abbreviated as imaging exceptional gait (IEG). In addition to traditional time-domain analysis, power spectrum analysis, and time-frequency analysis, a method for brain network analysis was also used to investigate the neural mechanisms and decoding of motor coordination imagery involving the lower limbs. e rest of the paper is organized as follows. Details of the method are illuminated in Section 2 including experiments configuration and data analysis. Section 3 displays the brain topographic maps, functional brain networks during MCI, IEG, and REST, and comparison results of the average classification accuracy of MCI and IEG. In Section 4, a brief discussion of the results is given. Finally, Section 5 concludes the paper.

Subjects.
Previous studies have shown that brain activity patterns evoked by motor imagery and actual exercise are similar [33]. Hence, our present study used the Kinesthetic and Visual Imagery Questionnaire (KVIQ) [38][39][40][41][42][43][44][45][46][47][48][49][50] to recruit 30 healthy subjects (numbered s1-s30) with strong motor imagery abilities (questionnaire score ≥ 70 points) to serve as participants for data collection. All recruited subjects were male, right-handed, without psychological disease and were 25 ± 1 years old. is study was approved by the Medical Ethics Committee of Kunming University of Science and Technology. Each subject signed an informed written consent before the start of the experiment.

Experimental Paradigm.
Coordinated movement is a type of exercise to restore and strengthen motor coordination. It usually involves the coordinated movement of multiple joints and limbs. Specifically, it includes motor coordination of the following: the upper and lower limbs, the limb trunk, and both symmetrical coordination and asymmetrical coordination of the left and right limbs. Among these components, coordinated movement of the upper limbs and hands can train the accuracies, reaction speeds, and rhythms of actions, while coordinated movement of the lower limbs can train correct gaits and motor coordination of the upper and lower limbs. In order to establish a method for recovering motor abilities of poststroke hemiplegic patients, the present study focused on motor coordination imagery of the lower limbs. We designed MCI and IEG. e timing of each trial is shown in Figure 1. e screen at the beginning of the experiment showed a fixed cross for 5 s, which prompted subjects to prepare for the experiment. When the fixed cross disappeared, the screen showed a prompt for normal gait imagery for 2 s, which prompted subjects to prepare for normal gait imagery. Once the prompt disappeared, the screen became blank and required the subjects to imagine a normal gait for 3 s. When the imagery period ended, the screen showed a star prompt for 5 s, indicating the subjects to rest. After the star prompt disappeared, the screen showed a prompt for poststroke hemiplegic gait imagery for 2 s, which prompted subjects to prepare for poststroke hemiplegic gait imagery. After the prompt disappeared, the screen became blank and required subjects to imagine a poststroke hemiplegic gait for 3 s. At the end of the imagery period, the screen showed a star prompt for 5 s, prompting the subjects to rest. is entire sequence comprised one experimental block, and each subject completed a total of 50 blocks.

Experimental Setup.
During the experiment, subjects were approximately 70 cm away from the prompt screen, and their hands were placed on the desktop. A TCL 24-inch LCD monitor was used for task prompts. A Lenovo inkPad computer running Windows 10 was used for data acquisition and processing via MATLAB software. e electroencephalogram (EEG) amplifier consisted of an NT9200 series (Beijing Zhongke Xintuo Instrument Co., Ltd.) with a sampling rate of 1000 Hz and 45 Hz low-pass filtering. At the beginning of the experiment, each subject completed the task according to the prompts in Figure 1. During the experiment, the subjects were required to avoid movements and blinking as much as possible during the visual imagery task. Figure 2 shows the schematic diagram of the electrode layout. e collected EEG data consisted of the following channels: Fp1 , Fp2, F7, F3, Fz, F4, F8, FT7, FC3,  FCz, FC4, FT8, T7, C3, Cz, C4, T8, TP7, TP8, CP3, CP4, CPz,  P7, P3, Pz, P4, P8, PO7, PO8, O1, O2, and Oz. e electrode layout included spatial positioning over the frontal lobe, parietal lobe, occipital lobe, and temporal lobe. e reference electrodes are denoted as A1 and A2.

Data Preparation.
Visual inspection was carried out on the collected EEG data to remove any data segments with substantial noise in the EEG signal. en, the EEG data of the 30 subjects under different task conditions were extracted.

Pretreatment.
First, baseline drift correction was performed on the extracted EEG signals to eliminate deviations of the EEG signal from the baseline [51]. An elliptic filter was used for 8-30 Hz band-pass filtering. en, independent component analysis (ICA) was used to remove electrooculogram artifacts and electromyographical artifacts. Additionally, alpha, beta, gamma, zeta, delta, and theta bands were extracted.

Time-Domain Analysis.
After preprocessing, the data were divided into three categories: MCI, IEG, and process of rest (REST). en, these three categories of data were stacked and averaged, and we compared the movementrelated cortical potentials (MRCPs) among them in the time domain and compared the MRCP average across channels and the grand average across subjects among these three categories of data.

Analysis of Brain Topography.
In order to determine the activation location of the brain region under the three tasks of MCI, IEG, and REST, the three categories of data were preprocessed and then averaged, after which the brain topographies during the three tasks were determined at 1, 100, 200, 300, 400, and 500 ms.

Power Spectrum Analysis.
In order to analyze the EEG power associated with the task and reduce the impact of absolute power between subjects and electrodes, the absolute powers, ABSP, of MCI and IEG were obtained by using the power of MCI minus the power of REST, as well as by using the power of IEG minus the power of REST. e quantitative method for determining the power spectral density was via a fast Fourier transform algorithm. ABSP was calculated as follows [52]: (1)

Brain Network Analysis.
In order to determine the connectivity of the brain function network under the three tasks of MCI, IEG, and REST, a multivariable autoregressive (MVAR) model was established to obtain a time-varying direct directed transfer function (dDTF). Because dDTF can eliminate the pseudocausal relationship derived from other leads [53], dDTF was used as a measure of causality, and Granger causality analysis was performed on the three types of preprocessed data to generate a dDTF weight matrix. e threshold value used in the present study was set to 0.2. en, we generated an adjacency matrix under the selected threshold. Finally, we obtained the brain function network under the three types of tasks.
Specifically, the original EEG data, X(t), of the N lead can be expressed as follows: Among them, X i (t) is the EEG time series of the ith lead. MVAR can be constructed according to X(t) [54]: Among them, A(k) is a matrix of model coefficients of N × N, which represents the dependence of the time delay, k; p is the model order, and E(t) is the random noise. Our present study used the Bayesian information criterion to determine the model order, p.
By Fourier transform of A(k), we obtain Furthermore, X(f) in the frequency domain can be obtained as follows: Among them, H(f) is the transfer matrix of the system. Furthermore, the dDTF can be obtained as [55] Among them, F 2 ij (f) is the full frequency DTF, and the calculation method is as follows: where C ij (f) is the partial coherence. e calculation method is as follows:

Computational Intelligence and Neuroscience
where M ij is a minor of spectral matrix (matrix of spectra and cross-spectra) with the i-th row and j-th column removed.

Measurement Analysis of the Brain Function Network.
In order to analyze the connectivity of the brain function network in detail, four kinds of network measures (characteristic path length, global efficiency, clustering coefficient, and local efficiency) were selected for comparative analysis of the three types of tasks. e characteristic path length, PL, and global efficiency, E glob , measure the global transmission capability of the network as follows [56]: In equations (9) and (10), the distance between two nodes, i and j, in the network is represented by d ij , which represents the edge number information of the shortest path, and N is the total number of nodes in the network. It can be determined from formulas (9) and (10) that the calculation method of global efficiency is similar to the characteristic path length. e difference is that the global efficiency, E glob , does not need to consider the isolated points in the network, which is suitable for a regular network, small world network, or a random network. e local efficiency E loc measures the local transmission capacity of the network. e local efficiency E loc of the network is the average of the global efficiency of the subgraphs corresponding to all nodes in the network: In equation (11), G i represents the subnetwork formed by all nodes directly connected with node i. e cluster coefficient CL a measures the degree of grouping of the network. e cluster coefficient CL a of a network is the mean value of the cluster coefficient cl i of all nodes in the network, where the cluster coefficient of node i is defined as follows: In equation (12), j,m a ij a im a jm is the number of edges included in the subnetwork composed of node i and is directly connected to nodes j and m. k i (k i − 1) is the maximum number of edges of the subnetwork, and k i is the degree of node i. erefore, the clustering coefficient CL a of a network is as follows: In equation (13), N represents the total number of nodes in the network.

Classification Verification.
In order to verify the difference between normal motor coordination imagery involving the lower limbs and that of poststroke hemiplegia gait imagery, the time-domain MRCP feature, power spectrum feature, time-frequency feature, network attribute feature, combined features of network attributes, and the adjacency matrix space were extracted. Since the SVM classifier has a good classification effect on small sample data [57], SVM was used for classification verification.
Specifically, the MRCP feature was used to extract the peak and phase of the 32 electrodes and to form a 32-dimensional feature matrix. e feature of the power spectrum was used to for the AR model to estimate the power spectrum and to obtain a 32-dimensional feature matrix. e time-frequency feature was to extract the time-frequency energy of the 32-lead EEGs by HHT and to form a 32-dimensional feature matrix. e network attribute feature was used to calculate the characteristic path length of the 32 electrodes to form a 32-dimensional characteristic matrix. e spatial feature of the adjacency matrix was used to extract the 16-dimensional feature matrix of the adjacency matrix using a spatial filter.

Statistical Analysis.
In order to test the significance of the task differences, the characteristic path lengths of 15 channels for the three types of tasks were compared in pairs by repeated-measures analysis of variance. e statistical significance level was set as p < 0.05.

Results
Figure 3(a) shows the average EEG waveforms across channels and the grand average across subjects during MCI. e colored line represents the average EEG data during MCI across all subjects, and the black line represents the grand average EEG data of MCI across the 32 channels. It can be seen from Figure 3(a) that a negative ERP was generated at about 200 ms. Figure 3(b) shows the average and grand average waveforms during IEG for all subjects. Figure 3(c) shows a comparison of the stack average of EEG data after preprocessing under MCI, IEG, and REST tasks. Figure 3(d) shows the phase difference of EEG data between MCI and REST, IEG and REST, and MCI and IEG. e EEG waveform trends of IEG and REST states were relatively consistent, while the EEG waveform trends of MCI and REST states were almost opposite to one another; additionally, the EEG waveform trends of MCI and IEG states were also opposite to one another (Figure 3(d)). Figure 4(a) shows the MRCP difference between the MCI and IEG on the C3 and C4 channels, and Figure 4(b) shows the MRCP phase difference between the MCI and IEG on the C3 and C4 channels. Figure 4 shows that, for the MCI MRCPs, the waveforms of the C3 and C4 channels were similar, but there was a significant difference at about 200 ms. For the IEG MRCPs, the waveforms of the C3 and C4 channels were nearly opposite to one another, and the waveforms of the C3 channel were half a cycle ahead of those of the C4 channel, while they reached a consistency at about 200 ms. Computational Intelligence and Neuroscience Figure 5 shows the brain topographic maps at 1, 100, 200, 300, 400, and 500 ms during MCI, IEG, and REST. At about 200 ms, during MCI, the brain topographic map was mainly activated in the right brain sensorimotor area (C4 channel); during IEG, the brain topographic map was mainly activated in the left brain sensorimotor area (C3 channel) and the parietal lobe (Pz channel), while during REST the brain topographic map was mainly activated in the left and right sensorimotor areas (C3 and C4).
In terms of timing, MCI first activated the right frontal lobe (F4 channel), left sensorimotor area (C3 channel), and parietal lobe (Pz channel) at 100 ms; the right sensorimotor area (C4 channel) was activated at 200 ms, whereas the left sensorimotor area (C3 channel) and parietal lobe (Pz channel) were activated at 300 ms. e central frontal area  Figure 6 shows the power spectrum during MCI, IEG, and REST tasks, as well as the absolute power spectrum of MCI and IEG. According to changes in the absolute power spectrum of MCI and IEG in Figures 6(d) and 6(e), the power distributions of MCI and IEG were in the θ wave (4-7 Hz), α wave (8-13 Hz), and β wave (14)(15)(16). e frequency band corresponding to the absolute power showed that the power related to MCI was mainly concentrated in the α wave, while the power related to IEG was mainly concentrated in the θ wave. Figure 7 shows the brain function networks of MCI, IEG, and REST at 1, 100, 200, and 300 ms. e main core nodes of the brain network were F7, FT7, T8, CP3, Pz, and PO8 at 1 ms for MCI, and the network connection took place from the left frontal lobe to the right occipital lobe (Figure 7). e main core nodes of the brain network were FP1, Fz, Oz, and PO8 at 100 ms, and the network connection took place from the forehead to the occipital lobe. e main core nodes of the brain network were F8, FT8, FCz, CP4, TP8, PO7, and O1 at 200 ms. e main core nodes of the brain network were FP2, F3, FT7, T8, TP8, PO8, and Oz at 300 ms, and the network connection occurred from the right forehead to the right occipital lobe. At 1 ms, the main core nodes during IEG were F7, FP2, F4, FCz, C4, P3, and O2, and the network connection took place in the left and right hemispheres; at 100 ms, the main core nodes of the brain network were F7, Fz, P4, and PO8, and the network connection took place from the left forehead to the right occipital lobe. At 200 ms, the main core nodes of the brain network were F4, FCZ, CP4, TP7, and Pz, and the network connection took place in the right hemisphere; at 300 ms, the main core nodes of the brain network were FCz, FT8, P7, Pz, and O1, and the network connection occurred from right forehead to left occipital lobe. When REST was 1 ms, the main core nodes of brain network were FP1, FP2, F7, FC3, T8, O1, and O2, and the network connection took place in the left prefrontal lobe and the right occipital lobe; when REST was 100 ms, the main core nodes of the brain network were PO7, CPz, CP4, P4, and TP8, and the network connection took place in the left and right occipital lobes. At 200 ms during REST, the main core nodes of the brain network were F7, T8, and O2, and the network connection took place in the left prefrontal lobe; when REST was 300 ms, the main core nodes of the brain network were FP2, CPZ, and PO8, and the network connection occurred from the right forehead to the occipital lobe. Figure 8 shows the comparison of brain network measures among MCI, IEG, and REST tasks. e four selected brain network measures were as follows: characteristic path length, global efficiency, clustering coefficient, and local efficiency. Among these four measures, the difference in the characteristic path lengths was the largest across MCI, IEG, and REST, while the other three measures exhibited only small differences across groups ( Figure 8). REST had the largest characteristic path length, followed by MCI and finally IEG. e characteristic path length showed that the speed of information transmission between nodes in IEG was the fastest, followed by MCI and REST.
In order to statistically analyze differences among MCI, IEG, and REST, the characteristic path lengths of the 15 channels during MCI, IEG, and REST were selected for repeated-measures single-factor analysis of variance. Table 1 shows the statistical results, where significant differences are marked in red. e difference in the characteristic path length between MCI and REST was significant, and only the O2 channel of MCI was not significantly different between MCI and REST (Table 1). ere was no significant difference between IEG and REST. Only the p values of P3, C3, CP4, and Pz channel were small during IEG, and there was a significant difference between MCI and IEG. Figure 9 shows the average classification results obtained by classifying MCI and IEG. e classifier adopted features of MRCPs, time-frequency analysis (Time-freq), power spectrum analysis (Power spec), network attribute analysis (Net-attr), a combination of the network attribute features, and the spatial features of the adjacency matrix (Net-attr + Adj-mat-space). e average classification accuracies of MRCPs, time-frequency analysis, power spectrum analysis, network attribute analysis, a combination of the network attribute features, and the spatial features of the adjacency matrix were 67.89% ± 9.14%, 78.12% ± 8.91%, 81.08% ± 8.99%, 89.12% ± 9.92%, and 92.96% ± 7.54%, respectively. According to the average classification results, the combination of the network attribute feature and the spatial feature of the adjacency matrix had the highest classification accuracies.   Figure 5: e brain topographic maps at 1, 100, 200, 300, 400, and 500 ms during MCI, IEG, and REST. e left, middle, and right columns show the brain topographic maps during MCI, IEG, and REST tasks, respectively. e brain topographic maps, from top to bottom, occurred at 1, 100, 200, 300, 400, and 500 ms during the corresponding task.
Computational Intelligence and Neuroscience

Discussion
Compared with the efficacy of traditional physical therapy, motor imagery therapy can induce brain plasticity, which can partially restore motor abilities of poststroke hemiplegic patients. At present, motor imagery therapy is mostly used for the recovery of upper-limb motor abilities but rarely for the recovery of lower-limb motor abilities. In addition, although there have been studies on the actual motor coordination of the lower limbs, few studies have investigated the neural mechanisms and decoding of imagined motor coordination of the lower limbs. erefore, the focus of our present study was distinct from that of traditional studies on motor coordination, such that we focused on motor coordination imagery involving the lower limbs.
First, our present study was different from the simple motor imagery used traditionally for a unilateral limb, such as studies that focused on the upper limbs [58][59][60][61][62][63][64], as well as lower limbs [65][66][67]. Prior to our present study, no related studies on EEG-based lower-limb motor coordination have been published. In order to determine a rehabilitation method for lower-limb motor coordination in poststroke hemiplegic patients, we used a paradigm of motor coordination imagery involving the lower limbs that consisted of both normal gait imagery and poststroke hemiplegic gait imagery. Simple motor imagery of a unilateral limb mainly activates the contralateral brain area [68], while motor coordination imagery involving both lower limbs involves coordination of the left and right lower limbs, which requires the coordinated participation of relevant brain regions on the left and right sides of the brain. is kind of motor coordination imagery training may be beneficial for inducing plasticity of brain network structure/ function or affected patients. In future research, we will design and improve our experimental research paradigm for poststroke hemiplegic patients with problems in motor coordination of the lower limbs. e EEG pattern induced by motor coordination imagery involving the lower limbs is closely related to the performance of motor imagery psychological activities of subjects [68]. In order to ensure the effective completion of the motor coordination imagery involving the lower limbs and the reliability of the corresponding data, before the experiment, subjects were selected according to their scores on a motor-vision imagery questionnaire; subjects with scores higher than 70 were selected for participation in our present study. After the experiment, for the designed imagery task, each subject was required to fill in a questionnaire about the implementation of psychological activities and the degree of difficulty of imagination during completion of the experiments in our present study. e subjects reported that they could complete MCI and IEG in a controllable way, but that for MCI was easier to imagine   Figure 7: Functional brain networks at 1, 100, 200, and 300 ms during MCI, IEG, and REST. e upper, middle, and lower rows show the brain function networks during MCI, IEG, and REST, respectively. From left to right, the brain function networks occur at 1, 100, 200, and 300 ms during the corresponding task, respectively. compared to that for IEG. In future research, we intend to use online neurofeedback to further improve subjects' motor coordination imagery performance of the lower limbs.
In terms of analysis content, previous studies have mainly focused on the control of redundant degrees of freedom, quantification of motor coordination, external environments, and evaluation of motor coordination, while our present study focused on the neural mechanisms and decoding of motor coordination imagery involving the lower limb. In the EEG analysis method, in addition to traditional sensorimotor rhythm analysis [58][59][60][61][62][63][64][65][66][67], we also used MRCPs and brain networks to study motor coordination imagery involving the lower limbs. Previous studies have shown that MRCPs can represent neural activities related to planning before each exercise, the execution of each exercise, and the end of each exercise [64]. Previous studies have also shown that features of brain networks can represent the coordination among brain regions during activities such as sensation, perception, representation, and cognition [53]. e motor coordination imagery in our present study involved planning, execution, and termination of motor imagery, as well as coordination between brain regions.
In the time domain, the EEG patterns of IEG and REST were more consistent, while the EEG patterns of MCI and REST were almost opposite to one another, and the EEG patterns of MCI and IEG were also nearly opposite to one another. ese findings may be due to the fact that the subjects were more familiar with resting states and walking 0.000 0.000 0.001 0.000 0.000 0.002 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 F4 0.000 0.002 0.007 0.000 0.002 0.014 0.009 0.000 0.000 0.001 0.000 0.001 0.002 0.000 0.000 FC3 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 FC4 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 C3 0.000 0.000 0.000 0.000 0.000 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 Cz 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 C4 0.000 0.000 0.000 0.000 0.000 0.001 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 CP3 0.  states but were not as familiar with poststroke hemiplegic gaits. Unusual imagery can trigger different EEG patterns [36]. In our present study, the MRCPs of MCI and IEG were significantly different in terms of their EEG signals in the C3 and C4 channels. e EEG waveforms of MCI in the C3 and C4 channels were similar, but there was a significant difference at about 200 ms. e EEG waveforms of IEG in the C3 and C4 channels were nearly opposite to one another. e waveform of the C3 channel was a half cycle ahead of the waveform of the C4 channel, whereas they reached a consistent value at about 200 ms. is result may be due to MCI mainly involving bilateral lower limbs, while IEG mainly involves unilateral lower limbs. MRCPs will trigger ERP waveforms around 200 ms [35], so the ERP waveform is formed in the contralateral brain.
As demonstrated by brain topography, motor coordination imagery involving the lower limbs is mainly activated in sensorimotor areas, which is due to motor coordination imagery of the lower limbs representing a form of motor imagery, the latter of which has been shown to activate sensorimotor areas [10]. In addition to activating sensorimotor areas, the visual area located in the occipital lobe also has corresponding activation, which may be because the cues provided to the subjects are processed through the visual pathway before being transmitted to the sensorimotor area [68].
In the present study, power spectrum analysis showed that the power distribution related to MCI was distributed in α waves, while the power related to IEG was distributed in θ waves. Simultaneously, the absolute power of the IEG was much higher than the absolute power of the MCI. is result may be due to MCI being easier to complete than IEG, resulting in MCI inducing lower energy.
Under the condition of motor coordination or imagery, external stimulation acts on sensory organs and causes nerve excitation.
e activation of brain-related areas and the coordinated control between brain areas transmit such nerve excitation to synergetic and antagonistic muscles, such that each muscle's activation is timed appropriately and ensures smooth completion of motor actions. Not only the activation of brain area but also the features of the brain network to represent the coordinated motion or imagination need to be calculated. In the present study, brain network analysis showed that brain functional connections during MCI, IEG, and REST mainly occurred in the right forehead, which may have been due to these psychological activities being related to the subconscious since human subconscious behavior is mainly related to the right forehead network [37,53,65,69]. Our present brain network analysis also showed that the characteristic path length was a good distinguishing feature for MCI and IEG and that the other three measures were not. In addition, because the characteristic path length and the global efficiency determined the efficiency of the entire information transmission of the network [53], the characteristic path length and the global efficiency of the three types of tasks were quite different. IEG had the smallest characteristic path length and the highest global efficiency, so its network-information transfer rate was the fastest, followed by MCI; the slowest was REST, which may have been due to IEG needing to recruit more neural resources.
In the present study, the features of MRCPs, time-frequency analysis, power spectrum analysis, network attribute analysis, a combination of the network attribute features, and the spatial features of the adjacency matrix in MCI and IEG were classified by SVM, and the average classification accuracies were 67.89 ± 9.14%, 78.12 ± 8.91%, 81.08 ± 8.99%, 89.12 ± 9.92%, and 92.96 ± 7.54%, respectively. e results showed that the average classification accuracy of the combination of the network attribute features and the spatial features of the adjacency matrix was 92.96%, confirming that these were the two best classification features. is result may be due to features representing collaborative information among related brain regions; also, the feature dimensionality was high, so the relevant features were extracted and retained comprehensively.
Statistical analysis showed that there was a significant difference in the characteristic path lengths between MCI and REST, for which only the O2 channel was not statistically significant.
ere was no significant difference between the characteristic path lengths between IEG and REST, whereas the p values of P3, C3, CP4, and Pz channels of IEG were small, but the difference between MCI and IEG was significant.
is result may be due to the fact that imagining MCI and REST is a relatively easy psychological activity, while imagining IEG is a relatively difficult psychological activity.
In Figure 5, there are some transition areas which are not clear by uncertainty. is is often more frequent in reality than one might think. We may segment each image obtained in order to delimit each activated area using a fuzzy procedure [70,71].

Conclusion
In order to determine a rehabilitation method for lowerlimb motor coordination in poststroke hemiplegic patients, a novel paradigm using lower-limb motor coordination imagery (normal gait imagery and poststroke hemiplegic gait imagery) was employed in the present study. Timedomain analysis, power spectrum analysis, time-frequency analysis, brain network analysis, and statistical analysis were used to investigate the neural mechanisms and decoding methods of motor coordination imaging of the lower limbs. e results showed that motor coordination imagery involving the lower limbs was mainly activated in sensorimotor areas and that the brain functional connection mainly occurred in the right forehead. e results also showed that the average classification accuracies of the combination feature of the network attributes and the spatial features of the adjacency matrix were 92.96% for normal gait imagery and poststroke hemiplegic gait imagery. Taken together, our findings may provide ideas for the rehabilitation of movement in poststroke hemiplegic patients based on motor coordination imagery involving the lower-limb braincomputer interface (BCI). Finally, our future directions are as follows: (1) EEG combined with fMRI to study the neural mechanisms of motor coordination imagery involving the lower limbs and (2) online verification of motor coordination imagery via BCI.

Computational Intelligence and Neuroscience
Data Availability e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.