<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Aging Neurosci.</journal-id>
<journal-title>Frontiers in Aging Neuroscience</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Aging Neurosci.</abbrev-journal-title>
<issn pub-type="epub">1663-4365</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fnagi.2022.892178</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Research on Rehabilitation Training Strategies Using Multimodal Virtual Scene Stimulation</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Xie</surname> <given-names>Ping</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/260414/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Wang</surname> <given-names>Zihao</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Li</surname> <given-names>Zengyong</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/343394/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Wang</surname> <given-names>Ying</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1779645/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Wang</surname> <given-names>Nianwen</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Liang</surname> <given-names>Zhenhu</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/175605/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Wang</surname> <given-names>Juan</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Chen</surname> <given-names>Xiaoling</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Key Laboratory of Measurement Technology and Instrumentation of Hebei Province, Institute of Electric Engineering, Yanshan University</institution>, <addr-line>Qinhuangdao</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>National Research Center for Rehabilitation Technical Aids</institution>, <addr-line>Beijing</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Peng Xu, University of Electronic Science and Technology of China, China</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Yunfa Fu, Kunming University of Science and Technology, China; Banghua Yang, Shanghai University, China</p></fn>
<corresp id="c001">&#x002A;Correspondence: Xiaoling Chen, <email>xlchen@ysu.edu.cn</email></corresp>
<fn fn-type="other" id="fn004"><p>This article was submitted to Neurocognitive Aging and Behavior, a section of the journal Frontiers in Aging Neuroscience</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>30</day>
<month>06</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>14</volume>
<elocation-id>892178</elocation-id>
<history>
<date date-type="received">
<day>08</day>
<month>03</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>21</day>
<month>04</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2022 Xie, Wang, Li, Wang, Wang, Liang, Wang and Chen.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Xie, Wang, Li, Wang, Wang, Liang, Wang and Chen</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>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.</p></license>
</permissions>
<abstract>
<p>It is difficult for stroke patients with flaccid paralysis to receive passive rehabilitation training. Therefore, virtual rehabilitation technology that integrates the motor imagery brain-computer interface and virtual reality technology has been applied to the field of stroke rehabilitation and has evolved into a physical rehabilitation training method. This virtual rehabilitation technology can enhance the initiative and adaptability of patient rehabilitation. To maximize the deep activation of the subjects motor nerves and accelerate the remodeling mechanism of motor nerve function, this study designed a brain-computer interface rehabilitation training strategy using different virtual scenes, including static scenes, dynamic scenes, and VR scenes. Including static scenes, dynamic scenes, and VR scenes. We compared and analyzed the degree of neural activation and the recognition rate of motor imagery in stroke patients after motor imagery training using stimulation of different virtual scenes, The results show that under the three scenarios, The order of degree of neural activation and the recognition rate of motor imagery from high to low is: VR scenes, dynamic scenes, static scenes. This paper provided the research basis for a virtual rehabilitation strategy that could integrate the motor imagery brain-computer interface and virtual reality technology.</p>
</abstract>
<kwd-group>
<kwd>brain-computer interface</kwd>
<kwd>motor imagery</kwd>
<kwd>virtual reality</kwd>
<kwd>neural activation</kwd>
<kwd>virtual rehabilitation</kwd>
</kwd-group>
<contract-num rid="cn001">U20A20192</contract-num>
<contract-num rid="cn002">21372001D</contract-num>
<contract-num rid="cn002">21372005D</contract-num>
<contract-sponsor id="cn001">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content></contract-sponsor>
<contract-sponsor id="cn002">Hebei Provincial Department of Bureau of Science and Technology<named-content content-type="fundref-id">10.13039/501100008238</named-content></contract-sponsor>
<counts>
<fig-count count="7"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="30"/>
<page-count count="10"/>
<word-count count="6222"/>
</counts>
</article-meta>
</front>
<body>
<sec id="S1" sec-type="intro">
<title>Introduction</title>
<p>The brain-computer interface (BCI) is a direct communication and control channel established between the human brain and computer or other electronic devices (<xref ref-type="bibr" rid="B26">Wolpaw et al., 2002</xref>). Through this channel, people can express ideas or manipulate equipment directly through the brain without any language or action, which can effectively enhance the ability of patients with severe physical disabilities to communicate with the outside world or to control the external environment to improve the quality of life of patients (<xref ref-type="bibr" rid="B19">Ren et al., 2004</xref>). Virtual reality (VR) technology creates a virtual world through the use of a new computer system, allowing the experiencer to integrate into the virtual environment and achieve mutual interaction. In addition to being used in the gaming industry, VR also can help stroke patients with flaccid paralysis perform purposeful training in a virtual environment, thereby improving the effect of rehabilitation training. Motor imagery (MI) is defined as the cognitive activity in which a subject imagines a movement without actually performing the movement (<xref ref-type="bibr" rid="B25">Vries and Mulder, 2007</xref>), and it is a common application paradigm in the field of brain-computer interface research (<xref ref-type="bibr" rid="B13">Mattia et al., 2012</xref>). This method realizes communication and control with external devices by imaging body movements (<xref ref-type="bibr" rid="B29">Xu et al., 2013</xref>).</p>
<p>Motor imagery brain-computer interface (MI-BCI) is a type of BCI that recognizes the patient&#x2019;s motor imagery intention by guiding the patient to perform motor imagery based on motor imagery therapy (<xref ref-type="bibr" rid="B6">Clark et al., 2019</xref>). It can be used to effectively remodel the central nervous system in patients with motor dysfunction (<xref ref-type="bibr" rid="B10">Guger et al., 2017</xref>). VR technology can provide patients with a more immersive training environment (<xref ref-type="bibr" rid="B18">Ren et al., 2020</xref>), help patients perform motor imagery more accurately (<xref ref-type="bibr" rid="B3">Bayliss and Ballard, 2000</xref>), and generate more easily identifiable electroencephalography (EEG) signals (<xref ref-type="bibr" rid="B11">Li et al., 2018</xref>). <xref ref-type="bibr" rid="B17">Remsik et al. (2016)</xref> reviewed 17 independent MI-BCI stroke rehabilitation studies, and they found that 16 produced significant treatment effects. Studies have shown that both motor imagery and actual movement can activate bilateral premotor areas (<xref ref-type="bibr" rid="B16">Pfurtscheller and Neuper, 2001</xref>), parietal lobes, basal ganglia, and cerebellum. Studies (<xref ref-type="bibr" rid="B15">Pfurtscheller and Neuper, 1997</xref>) also have shown that stroke patients can perform motor imagery to partially activate the damaged motor network (<xref ref-type="bibr" rid="B8">Fumanal-Idocin et al., 2021</xref>). In addition, studies have found that the rehabilitation treatment model combining BCI and VR technology is suitable for people suffering from stroke (<xref ref-type="bibr" rid="B21">Varsehi and Firoozabadi, 2021</xref>), depression, addiction and other diseases (<xref ref-type="bibr" rid="B20">Takenaka et al., 2021</xref>). VR-guided action offers the advantages of intuitive and specific actions and a strong sense of substitution (<xref ref-type="bibr" rid="B22">Velasquez-Martinez et al., 2020</xref>). Well-designed scene feedback can produce neural activation (<xref ref-type="bibr" rid="B24">Vourvopoulos et al., 2015</xref>).</p>
<p>Rehabilitation training strategies based on VR and MI-BCI have the following limitations (<xref ref-type="bibr" rid="B4">Benitez-Andonegui et al., 2020</xref>). First, the current rehabilitation strategy based on MI-BCI improves accuracy primarily by improving the algorithm without using a scene stimulation to improve the neural activation of the subject (<xref ref-type="bibr" rid="B27">Xiao and Fang, 2021</xref>), and thereby improving the quality of the EEG signals to improve accuracy (<xref ref-type="bibr" rid="B1">Bagarinao et al., 2020</xref>). Second, the virtual rehabilitation scene is singular (<xref ref-type="bibr" rid="B23">Vidaurre et al., 2020</xref>), the individual adaptability is poor, and few studies (<xref ref-type="bibr" rid="B12">Li et al., 2021</xref>) have compared the neural region activation and enhanced EEG signals feature mechanisms in different scenes (<xref ref-type="bibr" rid="B2">Barsotti et al., 2015</xref>). Third, most of the existing virtual rehabilitation training strategies are performed by observing virtual scenes on a computer screen (<xref ref-type="bibr" rid="B14">Moctezuma and Molinas, 2020</xref>). This training mode is not only less immersive, but also easily disturbed by the external environment, which introduces difficulties to the rehabilitation training (<xref ref-type="bibr" rid="B9">Gomez-Pilar et al., 2016</xref>). Finally, visual feedback training is lacking during rehabilitation. At the neural mechanism level, visual feedback training promotes brain plasticity changes and functional reorganization through the activation of the mirror neuron system, thereby promoting the recovery of motor function (<xref ref-type="bibr" rid="B5">Birbaumer et al., 2006</xref>).</p>
<p>In this study, to examine ways to improve the deep activation mechanism of motor nerves in stroke patients with flaccid paralysis, we designed different virtual scenes and training tasks to stimulate motor imagery from different angles. The virtual scenes include the following: limb MI in static scenes, such as text and pictures; limb MI in dynamic scenes of three-dimensional (3D) life and virtual games; and limb MI in the VR environment. This study compared and analyzed energy changes in the motor areas of the brain and the recognition rate of motor imagery before and after rehabilitation training using three-scene stimulation. This study verified the positive activation effect of rehabilitation training strategy on stroke patients and explored the brain activation mechanism using stimulation with different virtual scenes.</p>
</sec>
<sec id="S2" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec id="S2.SS1">
<title>Experimental Paradigm Under Virtual Scenes Stimulation</title>
<sec id="S2.SS1.SSS1">
<title>Design of Virtual Scenes and Training Tasks</title>
<p>As shown in <xref ref-type="fig" rid="F1">Figure 1</xref>, in this study, we designed a rehabilitation training strategy based on MI-BCI and VR. Multiple-evoked MI stimulation scenes (e.g., static scenes, dynamic scenes, and VR scenes) induced subjects to perform limb motor imagery. The EEG signals of the brain motor area were collected in real time and then the signals were subjected to preprocessing, intention feature extraction, and intention recognition. Finally, the results were output to the scene for interactive control.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p>Virtual rehabilitation training strategy map.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnagi-14-892178-g001.tif"/>
</fig>
<p>The virtual scenes were interesting and immersive, and the patients had a strong awareness of active training. In this study, we first built the scenes. Then, we added the image of the virtual left and right hands into the scenes and set the training tasks and audio-visual feedback. Finally, through human-computer interaction control, the subjects could feel that they were performing actual body movements. Then we guided the subjects to perform active limb motor imagery.</p>
<p>In addition, we observed the brain electrical activity mapping of the subject and analyzed the degree of neural activation during the subject&#x2019;s brain electrical activity mapping. We observed the activation changes in the brain areas through brain electrical activity mapping and analyzed the effect of training on the subjects. According to the current neural activity in the motor area of the brain, we adaptively adjusted the training scenes to ensure that the subjects could continuously achieve maximum activation of nerves and to accelerate the remodeling of nerve functions.</p>
<p>Different research protocols can lead to different promotion effects of the BCI rehabilitation training system on either the unaffected or affected side of the cerebral hemisphere (<xref ref-type="bibr" rid="B7">Dodd et al., 2017</xref>). Different subjects respond differently to different scenes, and the activation areas and intensity of the cerebral cortex and the EEG signal features also are different when the same subject performs different tasks. Considering these differences, this study used the Unity 3D platform to design the training scenes (static and dynamic scenes) in the computer screen as well as the training scenes in the VR environment. We also designed some life-skills training in the scene to feature rehabilitation tasks (e.g., holding goods, pouring water, and picking food).</p>
<p>As shown in <xref ref-type="fig" rid="F2">Figure 2</xref>, this study designed four different static and dynamic scenes without using VR. They included static text scenes, static picture scenes, dynamic 3D life scenes, and dynamic virtual game scenes.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p>Static and dynamic scenes without VR.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnagi-14-892178-g002.tif"/>
</fig>
<p>As shown in <xref ref-type="fig" rid="F3">Figure 3</xref>, we also designed scenes in the VR environment to enhance the immersion and realism of the subject&#x2019;s training in this study. The control method of the training tasks required the patient to perform motor imagery according to the prompt. In other words, when the patient&#x2019;s use of left (right) motor imagery was recognized, the corresponding limb in the scenes would perform the corresponding movement.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p>Training scenes in VR enviroment. <bold>(A)</bold> One scene for computer. <bold>(B)</bold> One scence for dinner.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnagi-14-892178-g003.tif"/>
</fig>
<p>The following describes the scene in detail.</p>
<p>Static scene: Mainly contains text, pictures, etc., the stimulation effect is poor, but the stimulation is more direct. Directly prompt the subjects to perform motor imagery through left, right, etc.</p>
<p>Dynamic scenes: Mainly include videos, games, etc., with general stimulation effects and strong interactivity. Realize motor imagery by playing videos or playing games.</p>
<p>VR scene: The main stimulation mode is the same as that under the computer monitor, but after VR rendering, the sense of immersion is strong, making the patient feel immersive, and the stimulation effect is the best.</p>
</sec>
<sec id="S2.SS1.SSS2">
<title>Experimental Design</title>
<sec id="S2.SS1.SSS2.Px1">
<title>Experimental Principle and Process</title>
<p>The purpose of this study was to verify the effect of the rehabilitation training strategy to improve the subjects&#x2019; limb control ability during motor imagery. At the same time, we studied the differences in brain activation when subjects were stimulated to perform motor imagery in different scenes. This enabled us to analyze the effect of neural activation when subjects performed motor imagery using stimulation with different virtual scenes.</p>
<p>In this study, we designed a controlled experiment to compare and analyze the changes in neural activity of the brain and to determine the recognition rate of motor imagery using stimulations with different virtual scenes to find the mechanism of neural deep activation. <xref ref-type="fig" rid="F4">Figure 4</xref> shows a schematic diagram of the experiment. We assessed the subjects&#x2019; motor imagery ability before training and then subjects performed multiple motor imagery training. After completing the training, participants finished a post-training assessment.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption><p>Experimental schematics.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnagi-14-892178-g004.tif"/>
</fig>
<p>We selected nine healthy college students in good mental condition as subjects (all male; average age: 24 &#x00B1; 2 years old). All subjects are right handedness. The experiment required that all the subjects had not completed similar experiments before and had no history of neurological diseases. All the subjects were informed of the research intention of the trial, the details of the study, and the potential dangers associated with the experiment. Additionally, all of the experiments were conducted within 3 h of the subject having eaten at noon, and each subject had to close his eyes and rest for 5 min before starting motor imagery. Doing so could relieve tension and anxiety and ensured that subjects were in a good mental state.</p>
<p>We grouped the nine subjects equally into three groups: static scene control group, called S1&#x2013;S3; dynamic scene experimental group, called S4&#x2013;S6; and VR scene experimental group, called S7&#x2013;S9. Except for the different virtual scenes for stimulating motor imagery training, all of the other conditions were the same. All the subjects performed a total of 17 days of experiments (14 days of motor imagery task training experiments and 3 days of motor imagery assessment experiments). The specific experimental process was as follows: Subjects completed three sets of enhanced motor imagery paradigm training per day, and each group included 40 trials. In other words, subjects had to complete a total of 120 limb motor imagery trials per day. In the 14-day motor imagery training using virtual scene stimulation, we obtained the correct rate of 1,680 limb motor images for each subject.</p>
<p>Subjects also performed three motor imagery assessment experiments. Three experiments were conducted on the day before the motor imagery training, the day after the training was completed for 7 days, and the day after the training was completed for 14 days. The subjects completed three sets of motor imagery tasks according to the assessment paradigm, and each group contained 30 motor imagery trials. We collected 19 channels of EEG data during the trials and obtained a total of 270 EEG data points for each subject in the motor imagery trials.</p>
</sec>
<sec id="S2.SS1.SSS2.Px2">
<title>Experimental Paradigm</title>
<p><xref ref-type="fig" rid="F5">Figure 5</xref> shows the experimental paradigm for the pre-assessment and post-assessment in the enhanced motor imagery training. A single trial contained three periods totaling 10 s. From 0 to 2 s, a red dot appeared in the center of the screen and then shrunk to remind the subjects to concentrate to start the next motor imagery. From 2 to 7 s, the subjects had to focus on the direction of the left/right movement of the red dot on the screen and performed left/right hand (or left/right limb) motor imagery. From 7 to 10 s, a plus sign (+) appeared in the center of the screen to remind the subjects that this trial had ended. During this period, subjects would rest for 3 s and then they would enter the next trial to repeat this experimental process.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption><p>Assessment paradigm before and after enhanced MI training.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnagi-14-892178-g005.tif"/>
</fig>
<p>The experimental paradigm of enhanced motor imagery training in different scenes is shown in <xref ref-type="fig" rid="F6">Figure 6</xref>. We divided the training process into three groups: the enhanced static scenes control group, the enhanced dynamic scenes experimental group, and the VR scenes experimental group. The motor imagery training process adopted an interactive mode of real-time motor imagery feedback. A single trial contained three periods totaling 10 s. From 0 to 1 s, a red dot appeared in the center of the screen to remind the subjects to concentrate on the start of the motor imagery. At the end of 1 s, the red dot disappeared, which reminded the subjects to start the limb motor imagery. From 1 to 5 s, subjects had to focus on the prompt information provided on the screen and performed left/right hand (or left/right limb) motor imagery. Different prompt information was given to different groups. From 5 to 7 s, the results of the classification model in the motor imagery appeared on the screen, and the virtual characters and limbs in the scenes were controlled to move accordingly. During this period, subjects took a short rest until the end of 7 s. Then, the subjects would enter the next trial and repeat this experimental process.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption><p>Experimental paradigm of enhanced MI training in different scenes.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnagi-14-892178-g006.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec id="S2.SS2">
<title>Electroencephalography Signals Acquisition and Preprocessing</title>
<p>The electrode distribution of EEG data acquisition adopted the international standard 10- to 20-electrode lead positioning. We set the reference electrode in the central area at the top of the head and set the sampling frequency to 1,000 Hz. The acquisition channels included 19 channels in the motor nerve-related area, namely, FC1, FC3, FC5, C1, C3, C5, CP1, CP3, CP5, Cz, FC2, FC4, FC6, C2, C4, C6, CP2, CP4, and CP6.</p>
<p>Most of the EEG acquisition devices were non-invasive dry electrodes. Although EEG acquisition was simpler and more convenient, it was unavoidable that various noises would appear in the acquisition process, such as EOG, ECG, and EMG. This seriously interfered with the subsequent analysis of EEG signals and affected the results of EEG signal analysis. Therefore, it was necessary to preprocess the collected raw EEG signals to remove noises.</p>
<p>The main noises included baseline drift caused by the device, 50 Hz power frequency interference, and EOG artifacts formed by blinking. We used the moving average method to remove baseline drift and used an adaptive 50 Hz filter to remove power frequency interference. We used an independent component analysis to remove EOG artifacts and used the sixth-order Butterworth filter as the band-pass filter. Finally, we acquired 8&#x2013;32 Hz EEG signals.</p>
</sec>
</sec>
<sec id="S3" sec-type="results">
<title>Results</title>
<sec id="S3.SS1">
<title>Analysis of Energy Changes in Brain Motor Areas</title>
<sec id="S3.SS1.SSS1">
<title>Brain Electrical Activity Mapping</title>
<p>Studies have shown that when subjects perform motor imagery, the energy features of EEG signals in the motor areas of the brain must change. In this study, we used brain electrical activity mapping to analyze the neural activation of multiple channels in the motor area, and visual neurofeedback was provided.</p>
<p>Brain electrical activity mapping (BEAM) is a commonly used method for multichannel EEG signal analysis. This method collects EEG signals from multiple channels of the head, performs fourth-order energy extraction, and then uses color bands and digitization to represent different gray levels. Finally, a dynamic BEAM can be drawn.</p>
</sec>
<sec id="S3.SS1.SSS2">
<title>Brain Motor Area Energy Analysis Based on Brain Electrical Activity Mapping</title>
<p>On the basis of these experiments, we obtained the EEG data of nine subjects when they conducted motor imagery assessment on the day before the motor imagery training and on the day after the training had been completed for 14 days. We obtained and normalized fourth-order cumulative energy value of each channel for a single trial for each subject. Then, we used the energy value of each channel for multiple trials of nine subjects, which then was used as the input parameter for drawing the BEAM. We obtained the spatial distribution map of the neural activation in different lead positions of the left and right limb motor imagery, as shown in <xref ref-type="fig" rid="F7">Figure 7</xref>.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption><p>Spatial distribution of the brain activation between the before and after enhanced MI training. <bold>(A)</bold> For control group. <bold>(B)</bold> For experimental group. <bold>(C)</bold> For experimental group II.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnagi-14-892178-g007.tif"/>
</fig>
<p><xref ref-type="fig" rid="F7">Figure 7</xref> shows the spatial distribution of neural activation in the brain regions of the subjects: S1, S4, and S8 before and after 10 trials of enhanced motor imagery training. <xref ref-type="fig" rid="F7">Figure 7A</xref> shows that after the subjects were trained in the rehabilitation strategy designed in this study, the activation range of motor nerves in the brain region was significantly expanded (i.e., the activation breadth of motor imagery nerves increased). <xref ref-type="fig" rid="F7">Figure 7B</xref> shows that after the subjects were trained in the rehabilitation strategy, the color of the activation area was darker (i.e., the ERD/ERS phenomenon was more obvious), which indicated that the activation depth of the motor imagery nerves increased. <xref ref-type="fig" rid="F7">Figure 7C</xref> shows that compared with static scenes and dynamic scenes, this strategy had a more obvious effect on the activation breadth and depth of motor imagery nerves in VR scenes.</p>
</sec>
</sec>
<sec id="S3.SS2">
<title>Electroencephalography Signals Feature Extraction, the Modeling of Electroencephalography Signals, and Analysis of the Recognition Rate of Motor Imagery</title>
<sec id="S3.SS2.SSS1">
<title>Electroencephalography Signal Feature Extraction</title>
<p>Feature extraction is the mathematical transformation or mapping of the input signal to obtain eigenvalues that are easier to observe and monitor. The widely used EEG signals that feature the extraction methods used during motor imagery mainly include the following: power spectrum estimation method, wavelet transform, independent component analysis, and common spatial pattern (<xref ref-type="bibr" rid="B28">Xin and Wang, 2017</xref>). In this study, we selected three feature combinations of mean square error, power spectral density, and common spatial pattern from the time domain, frequency domain, and spatial domain of EEG signal as the feature quantities for motor imagery intention recognition (<xref ref-type="bibr" rid="B30">Yang et al., 2012</xref>).</p>
<p>Furthermore, we used the existing particle swarm optimization (PSO)-support vector machine (SVM) as a classifier to perform motion imagery intention recognition.</p>
</sec>
<sec id="S3.SS2.SSS2">
<title>Analysis of the Recognition Rate of Motor Imagery</title>
<p>Studies have shown that there is a positive correlation between the degree of active participation, the activity degree of motor nerves, the performance of classification models, and the recognition rate of motor imagery when subjects perform motor imagery. The higher the recognition rate of motor imagery, the higher the active participation of the subjects in motor imagery and the stronger the control ability of motor nerves. The recognition rate of motor imagery can be used as an indicator to stimulate motor imagery ability in different scenes. In this study, the analysis of the recognition rate of motor imagery intention included the following two aspects: online MI training and offline MI assessment.</p>
<sec id="S3.SS2.SSS2.Px1">
<title>Online Motor Imagery Training Recognition Rate Analysis</title>
<p>We counted the recognition rate of the nine subjects who completed online MI training in virtual scenes within 14 days (<xref ref-type="table" rid="T1">Table 1</xref>). The recognition rate was the online recognition result obtained statistically using the experimental paradigm designed in this study. The feature quantities used in the classification model were mean square error, power spectral density, and common spatial pattern features. The classification model was PSO-SVM.</p>
<table-wrap position="float" id="T1">
<label>TABLE 1</label>
<caption><p>Statistical results of online recognition rate during MI training in different virtual scenes (%).</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">Group</td>
<td valign="top" align="center" colspan="3">Static scenes (control group)<hr/></td>
<td valign="top" align="center" colspan="3">Dynamic scenes (experimental group I)<hr/></td>
<td valign="top" align="center" colspan="3">VR scenes (experimental group II)<hr/></td>
</tr>
<tr>
<td valign="top" align="left">Training day</td>
<td valign="top" align="center">S1</td>
<td valign="top" align="center">S2</td>
<td valign="top" align="center">S3</td>
<td valign="top" align="center">S4</td>
<td valign="top" align="center">S5</td>
<td valign="top" align="center">S6</td>
<td valign="top" align="center">S7</td>
<td valign="top" align="center">S8</td>
<td valign="top" align="center">S9</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Day_01</td>
<td valign="top" align="center">60.9%</td>
<td valign="top" align="center">55.9%</td>
<td valign="top" align="center">57.6%</td>
<td valign="top" align="center">61.8%</td>
<td valign="top" align="center">58.5%</td>
<td valign="top" align="center">62.6%</td>
<td valign="top" align="center">62.6%</td>
<td valign="top" align="center">59.4%</td>
<td valign="top" align="center">60.3%</td>
</tr>
<tr>
<td valign="top" align="left">Day_02</td>
<td valign="top" align="center">65.7%</td>
<td valign="top" align="center">57.4%</td>
<td valign="top" align="center">60.7%</td>
<td valign="top" align="center">60.7%</td>
<td valign="top" align="center">64.8%</td>
<td valign="top" align="center">59.1%</td>
<td valign="top" align="center">64.1%</td>
<td valign="top" align="center">63.1%</td>
<td valign="top" align="center">64.7%</td>
</tr>
<tr>
<td valign="top" align="left">Day_03</td>
<td valign="top" align="center">63.4%</td>
<td valign="top" align="center">61.6%</td>
<td valign="top" align="center">65.2%</td>
<td valign="top" align="center">65.1%</td>
<td valign="top" align="center">67.4%</td>
<td valign="top" align="center">67.4%</td>
<td valign="top" align="center">67.4%</td>
<td valign="top" align="center">64.8%</td>
<td valign="top" align="center">71.1%</td>
</tr>
<tr>
<td valign="top" align="left">Day_04</td>
<td valign="top" align="center">69.1%</td>
<td valign="top" align="center">59.3%</td>
<td valign="top" align="center">61.5%</td>
<td valign="top" align="center">69.1%</td>
<td valign="top" align="center">59.3%</td>
<td valign="top" align="center">65.9%</td>
<td valign="top" align="center">71.8%</td>
<td valign="top" align="center">71.0%</td>
<td valign="top" align="center">73.9%</td>
</tr>
<tr>
<td valign="top" align="left">Day_05</td>
<td valign="top" align="center">67.4%</td>
<td valign="top" align="center">62.7%</td>
<td valign="top" align="center">70.7%</td>
<td valign="top" align="center">72.4%</td>
<td valign="top" align="center">62.6%</td>
<td valign="top" align="center">70.7%</td>
<td valign="top" align="center">65.7%</td>
<td valign="top" align="center">67.3%</td>
<td valign="top" align="center">74.5%</td>
</tr>
<tr>
<td valign="top" align="left">Day_06</td>
<td valign="top" align="center">65.9%</td>
<td valign="top" align="center">66.5%</td>
<td valign="top" align="center">65.9%</td>
<td valign="top" align="center">76.8%</td>
<td valign="top" align="center">65.7%</td>
<td valign="top" align="center">72.6%</td>
<td valign="top" align="center">76.8%</td>
<td valign="top" align="center">73.5%</td>
<td valign="top" align="center">72.2%</td>
</tr>
<tr>
<td valign="top" align="left">Day_07</td>
<td valign="top" align="center">56.5%</td>
<td valign="top" align="center">65.2%</td>
<td valign="top" align="center">67.6%</td>
<td valign="top" align="center">70.7%</td>
<td valign="top" align="center">70.7%</td>
<td valign="top" align="center">70.7%</td>
<td valign="top" align="center">72.6%</td>
<td valign="top" align="center">74.0%</td>
<td valign="top" align="center">70.7%</td>
</tr>
<tr>
<td valign="top" align="left">Day_08</td>
<td valign="top" align="center">67.7%</td>
<td valign="top" align="center">69.8%</td>
<td valign="top" align="center">56.6%</td>
<td valign="top" align="center">69.3%</td>
<td valign="top" align="center">69.3%</td>
<td valign="top" align="center">65.1%</td>
<td valign="top" align="center">73.2%</td>
<td valign="top" align="center">72.7%</td>
<td valign="top" align="center">75.9%</td>
</tr>
<tr>
<td valign="top" align="left">Day_09</td>
<td valign="top" align="center">69.3%</td>
<td valign="top" align="center">67.3%</td>
<td valign="top" align="center">63.5%</td>
<td valign="top" align="center">78.2%<xref ref-type="table-fn" rid="t1fns1">&#x002A;</xref></td>
<td valign="top" align="center">74.3%</td>
<td valign="top" align="center">74.3%</td>
<td valign="top" align="center">79.3%</td>
<td valign="top" align="center">79.2%<xref ref-type="table-fn" rid="t1fns1">&#x002A;</xref></td>
<td valign="top" align="center">80.8%</td>
</tr>
<tr>
<td valign="top" align="left">Day_10</td>
<td valign="top" align="center">72.4%</td>
<td valign="top" align="center">72.7%<xref ref-type="table-fn" rid="t1fns1">&#x002A;</xref></td>
<td valign="top" align="center">65.6%</td>
<td valign="top" align="center">72.6%</td>
<td valign="top" align="center">77.4%<xref ref-type="table-fn" rid="t1fns1">&#x002A;</xref></td>
<td valign="top" align="center">67.4%</td>
<td valign="top" align="center">76.6%</td>
<td valign="top" align="center">76.7%</td>
<td valign="top" align="center">74.2%</td>
</tr>
<tr>
<td valign="top" align="left">Day_11</td>
<td valign="top" align="center">70.1%</td>
<td valign="top" align="center">62.1%</td>
<td valign="top" align="center">69.3%</td>
<td valign="top" align="center">67.6%</td>
<td valign="top" align="center">69.1%</td>
<td valign="top" align="center">70.0%</td>
<td valign="top" align="center">82.4%<xref ref-type="table-fn" rid="t1fns1">&#x002A;</xref></td>
<td valign="top" align="center">72.1%</td>
<td valign="top" align="center">79.5%</td>
</tr>
<tr>
<td valign="top" align="left">Day_12</td>
<td valign="top" align="center">74.9%<xref ref-type="table-fn" rid="t1fns1">&#x002A;</xref></td>
<td valign="top" align="center">69.3%</td>
<td valign="top" align="center">70.7%</td>
<td valign="top" align="center">76.6%</td>
<td valign="top" align="center">68.4%</td>
<td valign="top" align="center">67.5%</td>
<td valign="top" align="center">73.4%</td>
<td valign="top" align="center">70.4%</td>
<td valign="top" align="center">83.0%</td>
</tr>
<tr>
<td valign="top" align="left">Day_13</td>
<td valign="top" align="center">68.2%</td>
<td valign="top" align="center">67.1%</td>
<td valign="top" align="center">71.6%<xref ref-type="table-fn" rid="t1fns1">&#x002A;</xref></td>
<td valign="top" align="center">72.5%</td>
<td valign="top" align="center">70.0%</td>
<td valign="top" align="center">73.5%</td>
<td valign="top" align="center">78.3%</td>
<td valign="top" align="center">76.8%</td>
<td valign="top" align="center">75.3%</td>
</tr>
<tr>
<td valign="top" align="left">Day_14</td>
<td valign="top" align="center">70.9%</td>
<td valign="top" align="center">67.6%</td>
<td valign="top" align="center">67.6%</td>
<td valign="top" align="center">74.2%</td>
<td valign="top" align="center">71.7%</td>
<td valign="top" align="center">74.0%<xref ref-type="table-fn" rid="t1fns1">&#x002A;</xref></td>
<td valign="top" align="center">80.0%</td>
<td valign="top" align="center">78.2%</td>
<td valign="top" align="center">88.8%<xref ref-type="table-fn" rid="t1fns1">&#x002A;</xref></td>
</tr>
<tr>
<td valign="top" align="left">Average recognition rate</td>
<td valign="top" align="center">67.3%</td>
<td valign="top" align="center">64.6%</td>
<td valign="top" align="center">65.3%</td>
<td valign="top" align="center">70.5%</td>
<td valign="top" align="center">67.8%</td>
<td valign="top" align="center">68.6%</td>
<td valign="top" align="center">73.2%</td>
<td valign="top" align="center">71.4%</td>
<td valign="top" align="center">74.6%</td>
</tr>
<tr>
<td valign="top" align="left">Standard deviation</td>
<td valign="top" align="center">0.045</td>
<td valign="top" align="center">0.042</td>
<td valign="top" align="center">0.046</td>
<td valign="top" align="center">0.052</td>
<td valign="top" align="center">0.051</td>
<td valign="top" align="center">0.043</td>
<td valign="top" align="center">0.060</td>
<td valign="top" align="center">0.057</td>
<td valign="top" align="center">0.070</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn id="t1fns1"><p><italic>&#x002A;denotes the highest MI recognition rate for each subject across the 14 days.</italic></p></fn>
</table-wrap-foot>
</table-wrap>
<p>Each subject performed 120 motor imagery trials per day. We counted the number of correct motor imagery classifications per day and calculated each subject&#x2019;s single-day recognition rate, average recognition rate, and standard deviation within 14 days. The data marked with an asterisk (&#x002A;) in the table indicate that the subject had the highest recognition rate on this day.</p>
<p>According to <xref ref-type="table" rid="T1">Table 1</xref>, by completing the enhanced motor imagery training task, the subjects&#x2019; motor imagery recognition rates improved to varying degrees. After training, the average recognition rate of motor imagery for three subjects (S1&#x2013;S3) in static scenes (control group) was between 64.5 and 67.5%, the average recognition rate of motor imagery for three subjects (S4&#x2013;S6) in dynamic scenes (experimental group I) was between 67.5 and 70.5%, and the average recognition rate of motor imagery for three subjects (S7&#x2013;S9) in VR scenes (experimental group II) was between 71.0 and 75.0%.</p>
<p>These data showed that the rehabilitation training strategy designed in this study could improve subjects&#x2019; limb motor imagery ability. The subjects&#x2019; motor imagination ability in these three scenarios is from weak to strong: static scene, dynamic scene, and VR scene.</p>
</sec>
<sec id="S3.SS2.SSS2.Px2">
<title>Offline Motor Imagery Assessment Recognition Rate Analysis</title>
<p>By analyzing the improvement range of the subjects&#x2019; motor imagery recognition rate under different training periods, we used an assessment method to evaluate whether the rehabilitation training strategy could improve the subjects&#x2019; motor imagery ability. We collected the EEG data of nine subjects when they conducted motor imagery assessments 1 day before the motor imagery training, the day after the training was completed for 7 days, and the day after the training was completed for 14 days. Each subject completed three sets of motor imagery tasks per day, and each set consisted of 30 random left and right hand (limb) motor imagery trials.</p>
<p>We subjected the collected data to offline recognition rate analysis. To obtain accurate results, we used a three-fold cross-validation method, as follows: First, we divided the 90 motor imagery datasets into three groups on average. One group was taken out each time as the test data, and the remaining data were used as the training data to complete the PSO-SVM model training. Then, we used the trained model to classify the test data. Finally, we averaged the obtained three groups of classification results as the final recognition rate. The statistical results are given in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<table-wrap position="float" id="T2">
<label>TABLE 2</label>
<caption><p>Offline evaluation of recognition rate on different virtual scenes before, during and after enhanced MI training (%).</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">Group</td>
<td valign="top" align="center">Subjects</td>
<td valign="top" align="center">Recognition rate before the training</td>
<td valign="top" align="center">Recognition rate after 7-day training</td>
<td valign="top" align="center">Increasing range (%)</td>
<td valign="top" align="center">Recognition rate after 14-day training</td>
<td valign="top" align="center">Increasing range (%)</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Static scenes</td>
<td valign="top" align="center">S1</td>
<td valign="top" align="center">63.6%</td>
<td valign="top" align="center">69.7%</td>
<td valign="top" align="center">6.1</td>
<td valign="top" align="center">72.8%</td>
<td valign="top" align="center">3.1</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">S2</td>
<td valign="top" align="center">61.6%</td>
<td valign="top" align="center">66.1%</td>
<td valign="top" align="center">4.5</td>
<td valign="top" align="center">68.9%</td>
<td valign="top" align="center">2.8</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">S3</td>
<td valign="top" align="center">64.1%</td>
<td valign="top" align="center">68.4%</td>
<td valign="top" align="center">4.3</td>
<td valign="top" align="center">72.1%</td>
<td valign="top" align="center">3.7</td>
</tr>
<tr>
<td valign="top" align="left">Dynamic scenes</td>
<td valign="top" align="center">S4</td>
<td valign="top" align="center">63.8%</td>
<td valign="top" align="center">69.6%</td>
<td valign="top" align="center">5.8</td>
<td valign="top" align="center">74.7%</td>
<td valign="top" align="center">5.1</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">S5</td>
<td valign="top" align="center">64.1%</td>
<td valign="top" align="center">71.4%</td>
<td valign="top" align="center">7.3</td>
<td valign="top" align="center">76.0%</td>
<td valign="top" align="center">4.6</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">S6</td>
<td valign="top" align="center">60.7%</td>
<td valign="top" align="center">66.4%</td>
<td valign="top" align="center">5.7</td>
<td valign="top" align="center">71.3%</td>
<td valign="top" align="center">4.9</td>
</tr>
<tr>
<td valign="top" align="left">VR scenes</td>
<td valign="top" align="center">S7</td>
<td valign="top" align="center">64.6%</td>
<td valign="top" align="center">74.4%</td>
<td valign="top" align="center">9.8</td>
<td valign="top" align="center">81.6%</td>
<td valign="top" align="center">7.2</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">S8</td>
<td valign="top" align="center">59.6%</td>
<td valign="top" align="center">72.5%</td>
<td valign="top" align="center">12.9</td>
<td valign="top" align="center">79.3%</td>
<td valign="top" align="center">6.8</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">S9</td>
<td valign="top" align="center">62.7%</td>
<td valign="top" align="center">73.4%</td>
<td valign="top" align="center">10.7</td>
<td valign="top" align="center">80.5%</td>
<td valign="top" align="center">7.1</td>
</tr>
</tbody>
</table>
</table-wrap>
<p><xref ref-type="table" rid="T2">Table 2</xref> shows that after training, the motor imagery recognition rate of the nine subjects was improved. In addition, the overall increase in the recognition rate in the static scenes (control group) and the dynamic scenes (experimental group I) was roughly the same, and the overall increase of the recognition rate in the VR scenes (experimental group II) was significantly higher than that of the other two groups. This result indicated that the rehabilitation training strategy designed in this study could improve the subjects&#x2019; EEG signal identification and ability to control EEG signals. Compared with static scenes and dynamic scenes, VR scenes had a more significant effect on improving subjects&#x2019; control of their motor nerves.</p>
</sec>
</sec>
</sec>
</sec>
<sec id="S4" sec-type="discussion">
<title>Discussion</title>
<p>The experiments designed in this study compared the degree of neural activation as well as the classification and recognition rate of motor imagery when subjects performed motor imagery using stimulation with different virtual scenes. The rehabilitation training strategy designed in this paper greatly improved the motor nerve activation of the subjects, accelerated the remodeling of the subjects&#x2019; neural functions, and improved the motor imagery ability of the subjects.</p>
<p>The experiments further investigated in which scene the subjects&#x2019; motor nerve activation layers were wider and deeper.</p>
<p>The virtual scenes designed in this study included static scenes (e.g., pictures, text), dynamic scenes (e.g., animation, games) displayed on a computer screen, and VR scenes. After completing the motor imagery training, we compared and analyzed the energy changes in the brain motor areas and the motor imagery recognition rate changes for the different subjects to explore the neural activation and motor imagery ability changes of subjects using stimulation of different virtual scenes.</p>
<p>First, we analyzed the energy changes in the brain motor areas of subjects under different scene simulation and analyzed the BEAM of one subject randomly selected from each of the three groups. By comparing the distribution of BEAM of each subject before and after training, we found that the range of energy distribution in the motor area after training expanded and the color depth deepened. This proved that after training with the rehabilitation strategy described in this paper, the neural activity intensity and energy in the motor areas of subjects increased, which could stimulate the remodeling of damaged motor nerves to a certain extent.</p>
<p>We further analyzed motor imagery recognition rates. First, the recognition rates of online training (<xref ref-type="table" rid="T1">Table 1</xref>) and offline assessment (<xref ref-type="table" rid="T2">Table 2</xref>) of the nine subjects under different virtual scene simulation remained at around 60%. The main reason that the recognition rate of the latter was slightly higher than that of the former was that the signal interference was less formed in the offline state. Second, from the data in <xref ref-type="table" rid="T1">Table 1</xref>, we found that with an increase in the number of training days, the motor imagery recognition rate of the nine subjects followed an overall upward trend. The fluctuation of the recognition rate in the short period was affected mainly by the mental state of the subjects during the experiment and the interference differences of the EEG signals in each acquisition. Third, by calculating the average recognition rate of 14 days, we found that the recognition rate of subjects in the VR scenes was significantly higher than that of the static scenes and dynamic scenes displayed on a computer screen. This result showed that the subjects could mobilize more nerve cells in the motor area to pursue regular physiological activities in the VR scenes, thereby improving the recognition of the subjects&#x2019; EEG signals. As shown in <xref ref-type="table" rid="T2">Table 2</xref>, to verify whether the improved recognition rate could be formed by the transient stimulation of the virtual scenes or by the change of the motor nerve activity mechanism, we adopted the same experimental paradigm of motor imagery for offline assessment. We found that the improvement of the recognition rate in the VR scenes was significantly better than that in the static and the dynamic scenes, which was consistent with previous results, and further indicated that the VR scenes had a better effect on improving the recognition of subjects&#x2019; EEG signals.</p>
</sec>
<sec id="S5" sec-type="conclusion">
<title>Conclusion</title>
<p>Today, the virtual rehabilitation that integrates MI-BCI and VR holds significant potential in the field of stroke rehabilitation, and this technology can greatly improve the rehabilitation effect for patients. To investigate how MI-BCI therapy can maximize the deep activation of subjects&#x2019; motor nerves and accelerate the remodeling mechanism of motor nerve function, we designed the following rehabilitation training strategy: we enhanced motor imagery training under different virtual scenes and compared and analyzed the degree of neural activation and the recognition rate of motor imagery in stroke patients after enhanced motor imagery training using stimulation with different virtual scenes. The experimental results showed that the motor imagery training using virtual scene stimulation could improve the motor nerve activation and motor imagery ability of the subjects. Compared with the static and dynamic scenes displayed on a computer screen, the VR scenes had a more significant effect in improving neural activation intensity and recognition rate.</p>
</sec>
<sec id="S6" sec-type="data-availability">
<title>Data Availability Statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="S7">
<title>Ethics Statement</title>
<p>The studies involving human participants were reviewed and approved by the Human Ethics and Administrative Council of Yanshan University. The patients/participants provided their written informed consent to participate in this study.</p>
</sec>
<sec id="S8">
<title>Author Contributions</title>
<p>PX, XC, and ZW designed this study. ZW wrote the manuscript. ZL, YW, NW, and ZL contributed to conception of the study. JW and XC wrote sections of the manuscript. All authors contributed to manuscript revision, and approved the submitted version.</p>
</sec>
<sec id="conf1" sec-type="COI-statement">
<title>Conflict of Interest</title>
<p>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.</p>
</sec>
<sec id="pudiscl1" sec-type="disclaimer">
<title>Publisher&#x2019;s Note</title>
<p>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.</p>
</sec>
</body>
<back>
<sec id="S9" sec-type="funding-information">
<title>Funding</title>
<p>This work was supported by National Natural Science Foundation of China (Grant Nos. U20A20192, 61901407, and 62076216) and the Key Research and Development Program of Hebei Province of China (Grant Nos. 21372001D and 21372005D).</p>
</sec>
<ack>
<p>We would like to thank editor and reviewers for helpful comments to improve the manuscript.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bagarinao</surname> <given-names>E.</given-names></name> <name><surname>Yoshida</surname> <given-names>A.</given-names></name> <name><surname>Terabe</surname> <given-names>K.</given-names></name> <name><surname>Kato</surname> <given-names>S.</given-names></name> <name><surname>Nakai</surname> <given-names>T.</given-names></name></person-group> (<year>2020</year>). <article-title>Improving Real-Time Brain State Classification of Motor Imagery Tasks During Neurofeedback Training.</article-title> <source><italic>Front. Neurosci.</italic></source> <volume>14</volume>:<issue>623</issue>. <pub-id pub-id-type="doi">10.3389/fnins.2020.00623</pub-id> <pub-id pub-id-type="pmid">32670011</pub-id></citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Barsotti</surname> <given-names>M.</given-names></name> <name><surname>Leonardis</surname> <given-names>D.</given-names></name> <name><surname>Loconsole</surname> <given-names>C.</given-names></name> <name><surname>Solazzi</surname> <given-names>M.</given-names></name> <name><surname>Sotgiu</surname> <given-names>E.</given-names></name> <name><surname>Procopio</surname> <given-names>C.</given-names></name><etal/></person-group> (<year>2015</year>). &#x201C;<article-title>A full upper limb robotic exoskeleton for reaching and grasping rehabilitation triggered by MI-BCI</article-title>,&#x201D; in <source><italic>IEEE International Conference on Rehabilitation Robotics</italic></source>, (<publisher-loc>Piscataway</publisher-loc>: <publisher-name>IEEE</publisher-name>), <pub-id pub-id-type="doi">10.1109/ICORR.2015.7281174</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bayliss</surname> <given-names>J. D.</given-names></name> <name><surname>Ballard</surname> <given-names>D. H.</given-names></name></person-group> (<year>2000</year>). <article-title>A virtual reality testbed for brain-computer interface research.</article-title> <source><italic>IEEE Trans. Rehabil. Eng.</italic></source> <volume>8</volume>:<issue>188</issue>. <pub-id pub-id-type="doi">10.1109/86.847811</pub-id> <pub-id pub-id-type="pmid">34823657</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Benitez-Andonegui</surname> <given-names>A.</given-names></name> <name><surname>Burden</surname> <given-names>R.</given-names></name> <name><surname>Benning</surname> <given-names>R.</given-names></name> <name><surname>Mckel</surname> <given-names>R.</given-names></name> <name><surname>Sorger</surname> <given-names>B.</given-names></name></person-group> (<year>2020</year>). <article-title>An Augmented-Reality fNIRS-Based Brain-Computer Interface: A Proof-of-Concept Study.</article-title> <source><italic>Front. Neurosci.</italic></source> <volume>14</volume>:<issue>346</issue>. <pub-id pub-id-type="doi">10.3389/fnins.2020.00346</pub-id> <pub-id pub-id-type="pmid">32410938</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Birbaumer</surname> <given-names>N.</given-names></name> <name><surname>Weber</surname> <given-names>C.</given-names></name> <name><surname>Neuper</surname> <given-names>C.</given-names></name> <name><surname>Buch</surname> <given-names>E.</given-names></name> <name><surname>Cohen</surname> <given-names>L.</given-names></name></person-group> (<year>2006</year>). <article-title>Physiological regulation of thinking: brain-computer interface (BCI) research.</article-title> <source><italic>Prog. Brain Res.</italic></source> <volume>159</volume> <fpage>369</fpage>&#x2013;<lpage>391</lpage>. <pub-id pub-id-type="doi">10.1016/S0079-6123(06)59024-7</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Clark</surname> <given-names>E.</given-names></name> <name><surname>Czaplewski</surname> <given-names>A.</given-names></name> <name><surname>Dourney</surname> <given-names>S.</given-names></name> <name><surname>Gadelha</surname> <given-names>A.</given-names></name> <name><surname>Korostenskaja</surname> <given-names>M.</given-names></name></person-group> (<year>2019</year>). &#x201C;<article-title>Brain-Computer Interface for Motor Rehabilitation</article-title>,&#x201D; in <source><italic>HCI International 2019 - Posters: 21st International Conference</italic></source>, (<publisher-loc>Orlando, FL</publisher-loc>: <publisher-name>HCII 2019</publisher-name>). <pub-id pub-id-type="doi">10.1088/1741-2552/aae4b8</pub-id> <pub-id pub-id-type="pmid">30260322</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dodd</surname> <given-names>K. C.</given-names></name> <name><surname>Nair</surname> <given-names>V. A.</given-names></name> <name><surname>Prabhakaran</surname> <given-names>V.</given-names></name></person-group> (<year>2017</year>). <article-title>Role of the Contralesional vs. Ipsilesional Hemisphere in Stroke Recovery.</article-title> <source><italic>Front. Hum. Neurosci</italic>.</source> <volume>11</volume>:<issue>469</issue>. <pub-id pub-id-type="doi">10.3389/fnhum.2017.00469</pub-id> <pub-id pub-id-type="pmid">28983244</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fumanal-Idocin</surname> <given-names>J.</given-names></name> <name><surname>Takac</surname> <given-names>Z.</given-names></name> <name><surname>Fernandez</surname> <given-names>J.</given-names></name> <name><surname>Sanz</surname> <given-names>J. A.</given-names></name> <name><surname>Bustince</surname> <given-names>H.</given-names></name></person-group> (<year>2021</year>). &#x201C;<article-title>Interval-valued aggregation functions based on Moderate deviations applied to Motor-Imagery-Based Brain Computer Interface</article-title>,&#x201D; in <source><italic>IEEE Transactions on Fuzzy Systems</italic></source>, (<publisher-loc>Piscataway</publisher-loc>: <publisher-name>IEEE</publisher-name>), <pub-id pub-id-type="doi">10.1109/TFUZZ.2021.3092824</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gomez-Pilar</surname> <given-names>J.</given-names></name> <name><surname>Corralejo</surname> <given-names>R.</given-names></name> <name><surname>Nicolas-Alonso</surname> <given-names>L. F.</given-names></name> <name><surname>&#x00C1;lvarez</surname> <given-names>D.</given-names></name> <name><surname>Hornero</surname> <given-names>R.</given-names></name></person-group> (<year>2016</year>). <article-title>Neurofeedback training with a motor imagery-based BCI: neurocognitive improvements and EEG changes in the elderly.</article-title> <source><italic>Med. Biol. Eng. Comput.</italic></source> <volume>54</volume> <fpage>1655</fpage>&#x2013;<lpage>1666</lpage>. <pub-id pub-id-type="doi">10.1007/s11517-016-1454-4</pub-id> <pub-id pub-id-type="pmid">26906278</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guger</surname> <given-names>C.</given-names></name> <name><surname>Allison</surname> <given-names>B.</given-names></name> <name><surname>Cao</surname> <given-names>F.</given-names></name> <name><surname>Edlinger</surname> <given-names>G.</given-names></name></person-group> (<year>2017</year>). <article-title>A Brain-Computer Interface for Motor Rehabilitation With Functional Electrical Stimulation and Virtual Reality.</article-title> <source><italic>Arch. Phys. Med. Rehabilit.</italic></source> <volume>98</volume>:<issue>e24</issue>. <pub-id pub-id-type="doi">10.1016/j.apmr.2017.08.074</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>F.</given-names></name> <name><surname>Fei</surname> <given-names>W.</given-names></name> <name><surname>Zhang</surname> <given-names>L.</given-names></name> <name><surname>Tao</surname> <given-names>Z.</given-names></name> <name><surname>Xu</surname> <given-names>P.</given-names></name></person-group> (<year>2018</year>). <article-title>The Dynamic Brain Networks of Motor Imagery: time-Varying Causality Analysis of Scalp EEG.</article-title> <source><italic>Int. J. Neural Syst.</italic></source> <volume>29</volume>:<issue>1850016</issue>. <pub-id pub-id-type="doi">10.1142/S0129065718500168</pub-id> <pub-id pub-id-type="pmid">29793372</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>Y.</given-names></name> <name><surname>Wei</surname> <given-names>Q.</given-names></name> <name><surname>Chen</surname> <given-names>Y.</given-names></name> <name><surname>Zhou</surname> <given-names>X.</given-names></name></person-group> (<year>2021</year>). <article-title>Transfer learning based on hybrid Riemannian and Euclidean space data alignment and subject selection in brain-computer interfaces.</article-title> <source><italic>IEEE Access</italic></source> <volume>9</volume> <fpage>6201</fpage>&#x2013;<lpage>6212</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2020.3048683</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mattia</surname> <given-names>D.</given-names></name> <name><surname>Pichiorri</surname> <given-names>F.</given-names></name> <name><surname>Molinari</surname> <given-names>M.</given-names></name> <name><surname>Rupp</surname> <given-names>R.</given-names></name></person-group> (<year>2012</year>). &#x201C;<article-title>Brain Computer Interface for Hand Motor Function Restoration and Rehabilitation</article-title>,&#x201D; in <source><italic>Towards Practical Brain-Computer Interfaces</italic></source>, <role>eds</role> <person-group person-group-type="editor"><name><surname>Brendan</surname> <given-names>Z. A.</given-names></name> <name><surname>Stephen</surname> <given-names>D.</given-names></name> <name><surname>Robert</surname> <given-names>L.</given-names></name> <name><surname>Jos&#x00E9;</surname> <given-names>D. R. M.</given-names></name> <name><surname>Anton</surname> <given-names>N.</given-names></name></person-group> (<publisher-loc>Berlin</publisher-loc>: <publisher-name>Springer Nature</publisher-name>), <pub-id pub-id-type="doi">10.1007/978-3-642-29746-5_7</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Moctezuma</surname> <given-names>L. A.</given-names></name> <name><surname>Molinas</surname> <given-names>M.</given-names></name></person-group> (<year>2020</year>). <article-title>EEG Channel-Selection Method for Epileptic-Seizure Classification Based on Multi-Objective Optimization.</article-title> <source><italic>Front. Neurosci.</italic></source> <volume>14</volume>:<issue>593</issue>. <pub-id pub-id-type="doi">10.3389/fnins.2020.00593</pub-id> <pub-id pub-id-type="pmid">32625054</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pfurtscheller</surname> <given-names>G.</given-names></name> <name><surname>Neuper</surname> <given-names>C.</given-names></name></person-group> (<year>1997</year>). <article-title>Motor imagery activates primary sensorimotor area in humans.</article-title> <source><italic>Neurosci. Lett.</italic></source> <volume>239</volume> <fpage>65</fpage>&#x2013;<lpage>68</lpage>. <pub-id pub-id-type="doi">10.1016/S0304-3940(97)00889-6</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pfurtscheller</surname> <given-names>G.</given-names></name> <name><surname>Neuper</surname> <given-names>C.</given-names></name></person-group> (<year>2001</year>). <article-title>Motor imagery and direct brain-computer communication.</article-title> <source><italic>Proc. IEEE</italic></source> <volume>89</volume> <fpage>1123</fpage>&#x2013;<lpage>1134</lpage>. <pub-id pub-id-type="doi">10.1109/5.939829</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Remsik</surname> <given-names>A.</given-names></name> <name><surname>Young</surname> <given-names>B.</given-names></name> <name><surname>Vermilyea</surname> <given-names>R.</given-names></name> <name><surname>Kiekoefer</surname> <given-names>L.</given-names></name> <name><surname>Abrams</surname> <given-names>J.</given-names></name> <name><surname>Evander Elmore</surname> <given-names>S.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>A review of the progression and future implications of brain-computer interface therapies for restoration of distal upper extremity motor function after stroke.</article-title> <source><italic>Expert Rev. Med. Devices</italic></source> <volume>13</volume> <fpage>445</fpage>&#x2013;<lpage>454</lpage>. <pub-id pub-id-type="doi">10.1080/17434440.2016.1174572</pub-id> <pub-id pub-id-type="pmid">27112213</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ren</surname> <given-names>S.</given-names></name> <name><surname>Wang</surname> <given-names>W.</given-names></name> <name><surname>Hou</surname> <given-names>Z. G.</given-names></name> <name><surname>Liang</surname> <given-names>X.</given-names></name> <name><surname>Shi</surname> <given-names>W.</given-names></name></person-group> (<year>2020</year>). &#x201C;<article-title>Enhanced Motor Imagery Based Brain-Computer Interface via FES and VR for Lower Limbs</article-title>,&#x201D; in <source><italic>IEEE Transactions on Neural Systems and Rehabilitation Engineering</italic></source>, (<publisher-loc>Piscataway</publisher-loc>: <publisher-name>IEEE</publisher-name>). <pub-id pub-id-type="doi">10.1109/TNSRE.2020.3001990</pub-id> <pub-id pub-id-type="pmid">32746291</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ren</surname> <given-names>Y.</given-names></name> <name><surname>Wang</surname> <given-names>G.</given-names></name> <name><surname>Cheng</surname> <given-names>M.</given-names></name></person-group> (<year>2004</year>). <article-title>The control on the artificial limb based on brain-computer interface.</article-title> <source><italic>CHIN. J. Rehabilit. Med</italic>.</source> <volume>19</volume> <fpage>330</fpage>&#x2013;<lpage>333</lpage>.</citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Takenaka</surname> <given-names>Y.</given-names></name> <name><surname>Suzuki</surname> <given-names>T.</given-names></name> <name><surname>Sugawara</surname> <given-names>K.</given-names></name></person-group> (<year>2021</year>). <article-title>Time course effect of corticospinal excitability for motor imagery.</article-title> <source><italic>Eur. J. Neurosci</italic>.</source> <volume>54</volume> <fpage>6123</fpage>&#x2013;<lpage>6134</lpage>. <pub-id pub-id-type="doi">10.1111/ejn.15404</pub-id> <pub-id pub-id-type="pmid">34328240</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Varsehi</surname> <given-names>H.</given-names></name> <name><surname>Firoozabadi</surname> <given-names>S.</given-names></name></person-group> (<year>2021</year>). <article-title>An EEG channel selection method for motor imagery based brain&#x2013;computer interface and neurofeedback using Granger causality.</article-title> <source><italic>Neural Netw.</italic></source> <volume>133</volume> <fpage>193</fpage>&#x2013;<lpage>206</lpage>. <pub-id pub-id-type="doi">10.1016/j.neunet.2020.11.002</pub-id> <pub-id pub-id-type="pmid">33220643</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Velasquez-Martinez</surname> <given-names>L. F.</given-names></name> <name><surname>Zapata-Castano</surname> <given-names>F.</given-names></name> <name><surname>Castellanos-Dominguez</surname> <given-names>G.</given-names></name></person-group> (<year>2020</year>). <article-title>Dynamic Modeling of Common Brain Neural Activity in Motor Imagery Tasks.</article-title> <source><italic>Front. Neurosci.</italic></source> <volume>14</volume>:<issue>714</issue>. <pub-id pub-id-type="doi">10.3389/fnins.2020.00714</pub-id> <pub-id pub-id-type="pmid">33328839</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vidaurre</surname> <given-names>C.</given-names></name> <name><surname>Haufe</surname> <given-names>S.</given-names></name> <name><surname>Jorajur&#x00ED;a</surname> <given-names>T.</given-names></name> <name><surname>M&#x00FC;ller</surname> <given-names>K.</given-names></name> <name><surname>Nikulin</surname> <given-names>V. V.</given-names></name></person-group> (<year>2020</year>). <article-title>Sensorimotor functional connectivity: a neurophysiological factor related to BCI performance.</article-title> <source><italic>Front. Neurosci.</italic></source> <volume>14</volume>:<issue>575081</issue>. <pub-id pub-id-type="doi">10.3389/fnins.2020.575081</pub-id> <pub-id pub-id-type="pmid">33390877</pub-id></citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vourvopoulos</surname> <given-names>A.</given-names></name> <name><surname>Cardona</surname> <given-names>J.</given-names></name> <name><surname>Badia</surname> <given-names>S.</given-names></name></person-group> (<year>2015</year>). &#x201C;<article-title>Optimizing Motor Imagery Neurofeedback through the Use of Multimodal Immersive Virtual Reality and Motor Priming</article-title>,&#x201D; in <source><italic>International Conference on Virtual Rehabilitation</italic></source>, (<publisher-loc>Valencia</publisher-loc>: <publisher-name>IEEE</publisher-name>), <pub-id pub-id-type="doi">10.1109/ICVR.2015.7358592</pub-id></citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vries</surname> <given-names>S. D.</given-names></name> <name><surname>Mulder</surname> <given-names>T.</given-names></name></person-group> (<year>2007</year>). <article-title>Motor imagery and stroke rehabilitation: a critical discussion.</article-title> <source><italic>J. Rehabilit. Med.</italic></source> <volume>39</volume> <fpage>5</fpage>&#x2013;<lpage>13</lpage>. <pub-id pub-id-type="doi">10.2340/16501977-0020</pub-id> <pub-id pub-id-type="pmid">17225031</pub-id></citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wolpaw</surname> <given-names>J. R.</given-names></name> <name><surname>Birbaumer</surname> <given-names>N.</given-names></name> <name><surname>McFarland</surname> <given-names>D. J.</given-names></name> <name><surname>Pfurtscheller</surname> <given-names>G.</given-names></name> <name><surname>Vaughan</surname> <given-names>T. M.</given-names></name></person-group> (<year>2002</year>). <article-title>Brain&#x2013;computer interfaces for communication and control.</article-title> <source><italic>Clin. Neurophysiol</italic>.</source> <volume>113</volume> <fpage>767</fpage>&#x2013;<lpage>791</lpage>. <pub-id pub-id-type="doi">10.1016/S1388-2457(02)00057-3</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xiao</surname> <given-names>X.</given-names></name> <name><surname>Fang</surname> <given-names>Y.</given-names></name></person-group> (<year>2021</year>). <article-title>Motor Imagery EEG Signal Recognition Using Deep Convolution Neural Network.</article-title> <source><italic>Front. Neurosci.</italic></source> <volume>15</volume>:<issue>655599</issue>. <pub-id pub-id-type="doi">10.3389/fnins.2021.655599</pub-id> <pub-id pub-id-type="pmid">33841094</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xin</surname> <given-names>X. U.</given-names></name> <name><surname>Wang</surname> <given-names>N.</given-names></name></person-group> (<year>2017</year>). <article-title>Feature extraction and classification of EEG signals in four kinds of motion imagination.</article-title> <source><italic>J. Nanjing Univ. Posts Telecommun.</italic></source> <volume>37</volume>, <fpage>18</fpage>&#x2013;<lpage>22</lpage>. <pub-id pub-id-type="doi">10.14132/j.cnki.1673-5439.2017.06.003</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xu</surname> <given-names>B. L.</given-names></name> <name><surname>Fu</surname> <given-names>Y. F.</given-names></name> <name><surname>Shi</surname> <given-names>G.</given-names></name> <name><surname>Yin</surname> <given-names>X. X.</given-names></name> <name><surname>Miao</surname> <given-names>L.</given-names></name> <name><surname>Wang</surname> <given-names>Z. D.</given-names></name><etal/></person-group> (<year>2013</year>). <article-title>Comparison of optical and concentration feature used for fNIRS-based BCI system using HMM.</article-title> <source><italic>Appl. Mech. Mat.</italic></source> <volume>385-386</volume> <fpage>1443</fpage>&#x2013;<lpage>1448</lpage>. <pub-id pub-id-type="doi">10.4028/www.scientific.net/amm.385-386.1443</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>B.</given-names></name> <name><surname>Lu</surname> <given-names>W.</given-names></name> <name><surname>He</surname> <given-names>M.</given-names></name> <name><surname>Liu</surname> <given-names>L.</given-names></name></person-group> (<year>2012</year>). <article-title>Novel feature extraction method for BCI based on WPD and CSP.</article-title> <source><italic>CHIN. J. Sci. Instr.</italic></source> <volume>33</volume> <fpage>2560</fpage>&#x2013;<lpage>2565</lpage>.</citation></ref>
</ref-list>
</back>
</article>
