Comparison of Feature Vector Compositions to Enhance the Performance of NIRS-BCI-Triggered Robotic Hand Orthosis for Post-Stroke Motor Recovery

: Recently, brain–computer interfaces, combined with feedback systems and goal-oriented training, have been investigated for their capacity to promote functional recovery after stroke. Accordingly, we developed a brain–computer interface-triggered robotic hand orthosis that assists hand-closing and hand-opening for post-stroke patients without su ﬃ cient motor output. In this system, near-infrared spectroscopy is used to monitor the a ﬀ ected motor cortex, and a linear discriminant analysis-based binary classiﬁer estimates hand posture. The estimated posture then wirelessly triggers the robotic hand orthosis. For better performance of the brain–computer interface, we tested feature windows of di ﬀ erent lengths and varying feature vector compositions with motor execution data from seven neurologically intact participants. The interaction between a feature window and a delay in the hemodynamic response signiﬁcantly a ﬀ ected both classiﬁcation accuracy (Matthew Correlation Coe ﬃ cient) and detection latency. The ‘preserving channels’ feature vector was able to increase accuracy by 13.14% and decrease latency by 29.48%, relative to averaging. Oxyhemoglobin combined with deoxyhemoglobin improved accuracy by 3.71% and decreased latency by 6.01% relative to oxyhemoglobin alone. Thus, the best classiﬁcation performance resulted in an accuracy of 0.7154 and a latency of 2.8515 s. The hand rehabilitation system was successfully implemented using this feature vector composition, which yielded better classiﬁcation performance.


Introduction
Clinical demand is growing for a new neurorehabilitation strategy in which post-stroke patients with insufficient or no remaining hand motor function can participate. Currently, these patients cannot benefit from Active Movement Therapy (AMT) [1], one of the widely applied rehabilitation therapies, because the level of remaining motor output in their affected hands is insufficient. Although passive movement therapy is an alternative option, it can only bring minimal motor function recovery [2]. Therefore, roughly one-third of post-stroke patients with motor impairments [3] cannot expect much functional recovery from current rehabilitation strategies.
Incorporating a brain-computer interface (BCI) [4] into neurorehabilitation has been recognized and investigated as a promising solution for post-stroke patients with complete or almost complete our previously tested prototype, and some of the data has been presented in [13]. In the current study, we have improved upon the prototype system by evaluating the different window lengths and feature vector compositions. This has substantially improved classification performance in terms of accuracy and detection latency.
Thus, we here demonstrate the implementation of our NIRS-BCI-triggered hand rehabilitation system that integrates a NIRS with our previously developed robotic hand orthosis [14]. We have evaluated how the following items have affected the performance (accuracy and detection latency) of the NIRS-BCI module: (1) different hemodynamic delays and lengths of the feature window, and (2) different feature vector compositions. The feature vector compositions that we tested were: (a) 'preserving channels' vs. 'averaging', (b) using [HbO] together with [HbR] vs. [HbO] or [HbR] separately, and (c) using all channels vs. criterion-selected channels. The purpose of this study was to describe the system implementation and to discuss the results of the evaluation as we develop our hand rehabilitation system.

System and Elements
The hand rehabilitation system consists of a robotic hand orthosis, the NIRS system, and operating personal computer (PC) with a monitor, as shown in Figure 1a. A user, sitting on a comfortable chair in a relaxed manner, performs (or imagines) a hand motion based on the instructions provided on the monitor, located at a 50-60 cm distance from the user's eyes. The NIRS system measures local cortical activation in real-time and transmits the data to the operating PC. The operating program on the PC processes the data and classifies it into one of two hand motion states: hand-closing and hand-opening. Then, the classification output is wirelessly transmitted to the orthosis. The operation program is coded and implemented with MATLAB ® . A system demonstration is available in Video S1 (please see Supplementary Materials).
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 13 [13]. In the current study, we have improved upon the prototype system by evaluating the different window lengths and feature vector compositions. This has substantially improved classification performance in terms of accuracy and detection latency. Thus, we here demonstrate the implementation of our NIRS-BCI-triggered hand rehabilitation system that integrates a NIRS with our previously developed robotic hand orthosis [14]. We have evaluated how the following items have affected the performance (accuracy and detection latency) of the NIRS-BCI module: (1) different hemodynamic delays and lengths of the feature window, and (2) different feature vector compositions. The feature vector compositions that we tested were: (a) 'preserving channels' vs. 'averaging', (b) using [HbO] together with [HbR] vs. [HbO] or [HbR] separately, and (c) using all channels vs. criterion-selected channels. The purpose of this study was to describe the system implementation and to discuss the results of the evaluation as we develop our hand rehabilitation system.

System and Elements
The hand rehabilitation system consists of a robotic hand orthosis, the NIRS system, and operating personal computer (PC) with a monitor, as shown in Figure 1a. A user, sitting on a comfortable chair in a relaxed manner, performs (or imagines) a hand motion based on the instructions provided on the monitor, located at a 50-60 cm distance from the user's eyes. The NIRS system measures local cortical activation in real-time and transmits the data to the operating PC. The operating program on the PC processes the data and classifies it into one of two hand motion states: hand-closing and hand-opening. Then, the classification output is wirelessly transmitted to the orthosis. The operation program is coded and implemented with MATLAB ® . A system demonstration is available in Video S1 (please see Supplementary Materials). consists of three basic components: a near-infrared spectroscopy (NIRS) system, a robotic hand orthosis, and an operating personal computer (PC). The NIRS system measures changes in local [Hb] in the motor cortex and sends them to the PC. The PC detects hand motor intentions from the NIRS data and wirelessly triggers the robotic orthosis; (b) Experimental Design. One session had a resting state period and a task period comprising 15 repeated unit trials. A unit trial was a sequence of ready, hand-closing, and hand-opening periods. The duration of the whole session was 810 s. consists of three basic components: a near-infrared spectroscopy (NIRS) system, a robotic hand orthosis, and an operating personal computer (PC). The NIRS system measures changes in local [Hb] in the motor cortex and sends them to the PC. The PC detects hand motor intentions from the NIRS data and wirelessly triggers the robotic orthosis; (b) Experimental Design. One session had a resting state period and a task period comprising 15 repeated unit trials. A unit trial was a sequence of ready, hand-closing, and hand-opening periods. The duration of the whole session was 810 s.

NIRS Measurement
We used the LABNIRS system (Shimadzu Corporation, Kyoto, Japan) to measure variation in [HbO] and [HbR] from local cortical areas. The LABNIRS system applies three wavelengths of near-infrared light (780 nm, 805 nm, and 830 nm) to measure hemoglobin concentrations. The system offers both measured light intensities and hemoglobin concentrations converted simultaneously in real time. We used the [HbO] and [HbR] directly offered by the system.
We used four pairs of NIRS emitting and detecting optodes, resulting in nine measurement channels. The distance between an emitter and a detector was 30 mm, enabling a measurement depth 20 mm from the scalp. The optodes were positioned on the scalp surface around the motor cortex (C3) according to the international 10-20 system ( Figure 2). In particular, emitting optode 4 was on the C3. We combined a Shimadzu FLASH holder kit (Flexible Adjustable Surface Holder) with a swimming cap for more stable fixation of the NIRS optodes and for prevention of the potential effect of ambient light. We set the sampling interval to 120 ms, resulting in the sampling rate of 8.33 Hz. This sampling interval was sufficient, because the hemodynamic responses induced by the motor task occurred at around 0.1 Hz [10].
Appl. Sci. 2019, 9, x FOR PEER REVIEW 4 of 13 We used the LABNIRS system (Shimadzu Corporation, Kyoto, Japan) to measure variation in [HbO] and [HbR] from local cortical areas. The LABNIRS system applies three wavelengths of nearinfrared light (780 nm, 805 nm, and 830 nm) to measure hemoglobin concentrations. The system offers both measured light intensities and hemoglobin concentrations converted simultaneously in real time. We used the [HbO] and [HbR] directly offered by the system.
We used four pairs of NIRS emitting and detecting optodes, resulting in nine measurement channels. The distance between an emitter and a detector was 30 mm, enabling a measurement depth ~20 mm from the scalp. The optodes were positioned on the scalp surface around the motor cortex (C3) according to the international 10-20 system ( Figure 2). In particular, emitting optode 4 was on the C3. We combined a Shimadzu FLASH holder kit (Flexible Adjustable Surface Holder) with a swimming cap for more stable fixation of the NIRS optodes and for prevention of the potential effect of ambient light. We set the sampling interval to 120 ms, resulting in the sampling rate of 8.33 Hz. This sampling interval was sufficient, because the hemodynamic responses induced by the motor task occurred at around 0.1 Hz [10].

Robotic Hand Orthosis
For assisting hand motion, we used our robotic hand orthosis [14]. The exoskeleton-type orthosis assists with hand-closing and hand-opening motions, with a force of 10 N and 30 N, respectively. The orthosis consists of the wearable, exoskeletal hand part and an external control box ( Figure 2). The exoskeleton is attached to the dorsal side of the hand and covers both the dorsum and the fingertips. A linear actuator on the dorsum actuates all finger components simultaneously in the same direction, which thus enables flexion/extension of all fingers. The control box wirelessly controls the actuator.

Robotic Hand Orthosis
For assisting hand motion, we used our robotic hand orthosis [14]. The exoskeleton-type orthosis assists with hand-closing and hand-opening motions, with a force of 10 N and 30 N, respectively. The orthosis consists of the wearable, exoskeletal hand part and an external control box ( Figure 2). The exoskeleton is attached to the dorsal side of the hand and covers both the dorsum and the fingertips. A linear actuator on the dorsum actuates all finger components simultaneously in the same direction, which thus enables flexion/extension of all fingers. The control box wirelessly controls the actuator.

Wireless Transmission
A pair of XBee ® S1 802.15.4 RF modules offers wireless transmission of control commands from the control box to the wearable portion of the orthosis. There are only two control commands that correspond to the two hand motions (open and close). The default posture is hand-open.

Participants
We recruited seven neurologically intact participants for this study (all right-handed males, aged 29.9 ± 5.0 years). No participant had a history of neurological or psychiatric disorders, and none had ever participated in any BCI-related experiment. We obtained informed, written consent from all participants before the experiment. The study was approved by the institutional review boards and ethics committees of Kyushu University Hospital. All procedures were conducted in accordance with the latest version of the Declaration of Helsinki.

Experimental Protocol
All participants completed one session of the motor execution task. In this task, participants closed their right hands during a given period (default posture: open hand). A session consisted of a resting state (60 s) followed by 15 unit trials. A unit trial was a sequence of ready (5 s), hand-closing (15 s), and hand-opening (30 s) periods. Thus, the total length of the session was 810 s (Figure 1b). During the hand-closing period, participants loosely closed their right hands and kept them closed until the next instruction. During resting, ready, and hand-opening periods, they kept their hands opened in a relaxed manner. There was no intended finger hyperextension.

NIRS Channel Selection
We used the Contrast-to-Noise Ratio (CNR) of [Hb] as an indicator of task-induced hemodynamic response. Thus, each channel had two CNR values, one corresponding to [HbO] and the other to [HbR]. The CNR compares the amplitude contrast between hand-closing and hand-opening states. A task-induced response has a CNR value greater than zero for [HbO] and lower for [HbR] as in (1). A larger absolute value means a greater task-induced response. CNRs have been used for selecting NIRS channels containing task-evoked hemodynamic responses [15,16].
The entire hand-closing period (15 s) and the final 15 s of the hand-opening period were used for calculating the CNRs for the respective periods. Here, the CNR values for a channel are the mean CNRs of 15 trials after bandpass filtering.

Preprocessing
For real-time bandpass filtering, we applied a MACD (moving average convergence/divergence) filter (2). The MACD filter is the difference between two EMA (exponential moving average) filters (3) of different parameters. Several NIRS-BCI studies have tested [17] and successfully used MACD for their systems [16,18].
EMA αF is for fast components and EMA αS for slow ones. The applied passband was between 0.01 and 0.2 Hz to remove global ascending trend and the effect of physiological activities, such as respiration, heart-beating, and so on [19]. The corresponding parameters of α F and α S were calculated by (4) [20] (F C = f C / f S , f C: the cut-off frequency, f S : the sampling frequency).

Binary Classification
Linear discrimination analysis (LDA) was used to solve the binary classification of the two hand motions. An LDA-classifier finds a line that maximizes the difference between the mean values of two data and minimizes the variance within each individual class belonging to a training data. Then, the classifier applies the line to classify a test datum [21]. LDA has been one of the successfully applied binary classification algorithms for the NIRS-BCI systems [10]. Additionally, our previous work [13] demonstrated that LDA exhibited better classification performance than a support vector machine (SVM). We used MATLAB ® Machine Learning Toolbox™ to train and test the LDA-based classifier.

Different Lengths of the Feature Window
The length of the feature window determines the time interval in which the features, mean and slope, are calculated with MATLAB functions, mean and polyfit (set 'Degree of polynomial fit' to 1, meaning 'linear fitting') respectively. The length of the feature window can affect classification performance. Thus, we investigated how classification results varied depending on the length of the feature window (L FW , in s). We tested nine different feature window lengths (0.6, 1, 3, 5, 7, 9, 11, 13, and 15 s). The longest length, 15 s, corresponds to the length of a hand-closing period (task period). The shorter lengths, except for 0.6 s, were determined by narrowing down the maximum length with a decrement of 2 s. The shortest one, 0.6 s, matches with five data points. Additionally, delay times in the task-induced hemodynamic response (D HR ) from 0 to 10 s were introduced, ranging in increments of 1 s. These delays are time gaps between the onset of hand-closing and the visible change in hemoglobin concentrations.

Feature Vector Composition
We assessed the classification performance of 12 feature vector compositions (FCs) ( Table 1)

Classification Evaluation
We performed leave-5-out cross-validation after excluding the resting period. We partitioned the trials into training and test data sets, and the order of the trials in each set was preserved. This was for evaluating the detection latency.
We used two indices: (1) classification accuracy and (2) detection latency. As measures of classification accuracy, we used Balanced Accuracy (BACC) (5) and Matthews Correlation Coefficient (MCC) (6) because our data were class-imbalanced due to the experimental design. A detection latency is the time difference (in s) between the start time of a trial and the classifier-detected start of the motion. We labeled a trial as 'no-detection' if its detection latency was longer than the hand-closing period of 15 s. Then, we excluded 'no-detection' trials when averaging the detection latencies over all test trials.

Classifier Training with Different Delays and Window Lengths
Applying different D HR and L FW affected classification performance. Figure 3 shows the MCC RS and detection latency maps. Depending on the applied D HR and L FW , classification accuracy rose from 0.6489 up to 0.7880 (a 21.42 % increase) in the MCC RS . The D HR and L FW were positively related to each other, as evidenced by the antidiagonal direction of the map. Longer D HR with longer L FW tended to generate better classification accuracy. However, after D HR reached 9 s and L FW reached 13 s, classification accuracy decreased. The combination of long D HR and short L FW , or vice versa resulted in a relatively poor classification accuracy of around 0.65.
Detection latency also depended on the D HR and L FW , ranging between 0.1640 s and 4.7761 s. The detection latency map in Figure 3 was obtained based on the introduced D HR . To obtain the true detection latency value, the corresponding D HR must be added to the value of each pixel.

Classifier Training with Different Delays and Window Lengths
Applying different DHR and LFW affected classification performance. Figure 3 shows the MCCRS and detection latency maps. Depending on the applied DHR and LFW, classification accuracy rose from 0.6489 up to 0.7880 (a 21.42 % increase) in the MCCRS. The DHR and LFW were positively related to each other, as evidenced by the antidiagonal direction of the map. Longer DHR with longer LFW tended to generate better classification accuracy. However, after DHR reached 9 s and LFW reached 13 s, classification accuracy decreased. The combination of long DHR and short LFW, or vice versa resulted in a relatively poor classification accuracy of around 0.65.
Detection latency also depended on the DHR and LFW, ranging between 0.1640 s and 4.7761 s. The detection latency map in Figure 3 was obtained based on the introduced DHR. To obtain the true detection latency value, the corresponding DHR must be added to the value of each pixel.  The pixel values of the detection latency map were calculated on the basis of assumed delay in hemodynamic response. Thus, for true detection latency values, the corresponding D HR must be added to each pixel value. Note that each pixel is the average from the seven participants. The pixel values for both maps are given in Table S1 (Supplementary Materials).

All-and CNR-Selected Channels
The compositions with all channels (FC-1 to 3 and FC-7 to 9) performed better than the ones with CNR-selected channels (FC-4 to 6 and FC-10 to 12). Specifically, the all-channels compositions increased MCC RS by 1.86% (p = 0.11) and curtailed detection latency by 23.05% (p < 0.05) ( Table 2). Finally, we carried out one-way ANOVA tests to examine if the 12 different feature vector compositions affected classification accuracy and detection latency, respectively. The results indicate that the compositions significantly improved both classification accuracy (for MCC RS , F = 4.69; p < 0.00005) and detection latency (F = 2.5; p < 0.05). results are displayed after being grouped. ChPrsv means 'channel-preserving, and ChAvg means 'channel-averaging'. All indicates 'using all channels', whereas Contrast-to-Noise Ratio (CNR) means 'CNR-selected'. The values used for the bar graph are accessible in Table S1 (Supplementary Materials).

All-and CNR-selected Channels
The compositions with all channels (FC-1 to 3 and FC-7 to 9) performed better than the ones with CNR-selected channels (FC-4 to 6 and FC-10 to 12). Specifically, the all-channels compositions increased MCCRS by 1.86% (p = 0.11) and curtailed detection latency by 23.05% (p < 0.05) ( Table 2). Finally, we carried out one-way ANOVA tests to examine if the 12 different feature vector compositions affected classification accuracy and detection latency, respectively. The results indicate that the compositions significantly improved both classification accuracy (for MCCRS, F = 4.69; p < 0.00005) and detection latency (F = 2.5; p < 0.05).

Discussion
The final goal of the rehabilitation system presented in this study is to bring the functional recovery of the impaired hands of post-stroke patients. The essential element of the system is the NIRS-BCI module, and the most important characteristic is how well NIRS-BCI can detect motor intention. We therefore tried several approaches and compared the results to obtain the best performance of the BCI module. To determine the best accuracy, we compared the classification results of differing (1) delays in hemodynamic responses and lengths of the feature window, (2) channel composition, (3) the NIRS signal, and (4) channel selection.
Accuracy and latency depended on D HR and L FW selections, which turned out to be related to each other. This suggests that we must carefully select D HR when L FW depending on the purpose of the NIRS-based system. Although a D HR between 8 s and 9 s and an L FW between 11 s and 13 s guaranteed the maximum classification accuracy, detection latency was distributed between 9 s and 11 s. These combinations are acceptable and could be useful if the system is for offline analysis, where classification accuracy is of importance and detection latency can be ignored. However, these selections are not suitable for real-time use. When a D HR is between 0 s and 2 s and an L FW is around 1 s, some classification accuracy is lost, but detection latency is shortened to between 2 s and 3 s. These selections are preferred for a real-time NIRS system, such as what we have presented here.
Channel-preserved feature vector compositions generated better performance than channel-averaged compositions in terms of classification accuracy and detection latency. Classification after averaging all channels have reportedly produced accuracies of 78% [22] and 65-75% [23]. This approach might be effective and useful for a BCI system with fewer channels that is only concerned with one cortical area. This is indeed the case for our system of only nine channels that are focused on the motor cortex, especially when all channels record similar temporal patterns. However, our comparisons revealed that the channel-preserving composition improved classification accuracy and shortened detection latency. Cui et al. (2010) also reported that the feature space taking account of spatial information (i.e., channel-preserved) increased accuracy by 7.7% and shortened detection latency by 2.4 s. Thus, even though channel-averaging might guarantee an acceptable accuracy and is computationally effective, by reducing the dimensions of feature vectors, we recommend the channel-preserving feature vector for an accurate BCI with short detection latency.
Our results also indicated that CNR is not a perfect criterion for selecting informative NIRS channels for classification. CNR-selected channels could not achieve the levels of both classification accuracy and detection latency that were obtained using all channels ( Table 2). These results might indicate that there were some channels whose [HbO] or [HbR] correlated with the imagined hand motions but did not meet the CNR criteria. Thus, it might suggest that additional criteria are necessary. The coefficient of variance or the point-biserial correlation coefficient can be a good candidate to investigate. However, applying CNR can diminish the number of channels that are incorporated into a feature vector, thereby reducing computational time and cost. CNR could be more useful and effective if a system adopts many more NIRS channels, possibly producing more uninformative or noisy channels. This needs to be examined in future studies.
The concentration of HbR helps contribute to classification. The contribution of [HbR] alone is limited in the sense that classification improved when [HbO] was included in the feature vector. The feature compositions that included only [HbR] led to poorer classification performance, which dropped to the level of random guessing ( Table 2). Based on our analysis, the worst approach is incorporating only [HbR] into a feature vector after averaging. In contrast, the compositions with both [HbO] and [HbR] improved both classification accuracy and detection latency over those with only [HbO] through all of the tested compositions, across all participants. Most NIRS-BCI systems [10] have not favored [HbR] because of its unreliability, relating to its small dynamic range and its signal-to-noise ratio (SNR). However, [HbR] is conditionally and marginally useful for better classification.
The approach and results of this study with the neurologically intact participants could be translated into clinical cases. We consider that the final user of our system will be patients who suffer from hand motor impairments owing to subcortical stroke. Because of their impairments, the patients must perform motor imagery, not motor execution. However, because of the lesion location being subcortical, the ipsilesional cortical areas are intact and can produce motor imagery-induced activations, which are detectable. The MI-evoked hemodynamic responses by a subcortical stroke patient could be similar to that of a neurologically intact participant. Studies of post-stroke patients that use hemodynamics-related modalities, such NIRS, have revealed that trained patients could produce clear motor imagery-induced cortical activity similar to motor execution [12]. In addition, the clinical studies of the BCI systems based on electrophysiological responses, being neurovascularly coupled with hemodynamic responses, demonstrate the existence of classifiable, MI-evoked responses. In the magnetoencephalography (MEG)-and sensorimotor rhythm (SMR)-BCI studies, post-stroke patients having subcortical lesions could learn how to control a BCI system [24,25]. The successful cases of EEG-BCI-based rehabilitation indirectly support a post-stroke patient's ability to induce and control a task-related cortical activation [26,27]. Moreover, the pattern of a task-evoked hemodynamic response after a subcortical stroke seems preserved. The fMRI study of patients with cerebral ischemia [28] demonstrated an altered pattern of task-induced activation, i.e., an increase of [HbR] during a task, but this kind of alteration in a patient with a subcortical stroke has not been reported to our knowledge. Even though the altered activation can take place with a subcortical stroke, it implies that there is a task-related activation that exists and that is differentiable from a baseline. Thus, it can be detectable by a machine learning-based classifier, such as LDA.
This study has three main limitations. First, the sample size was small, and evaluation only used one technique (LDA) and two features (mean and slope). The inclusion of non-tested features, such as variance, skewness, and so on, and the testing of different classification techniques, such as SVM, neural networks, hidden Markov model, and naïve Bayesian with sets of tested and non-tested features could bring different results. Second, we must note that the two features here were calculated within a feature window of the same length. Each specific feature might contribute more to classification with a window of its own optimized length. Further research should be conducted for clarifying these problems. Finally, the influence of task-evoked extracerebral activation on our NIRS-BCI performance was not excluded due to the limitation of our NIRS system setup. The hemodynamics in the extracerebral layers of the head, such as scalp, can be affected by task-evoked changes in scalp blood flow and/or volume [29] and by those in the autonomic nervous system activity [30], which could contaminate NIRS signals. Although our single source-detector separation (SDS) of 3 cm could not filter out those extracerebral influences, it is unlikely that the inclusion of those extracerebral influences might compromise our study result. This is because the BCI study with a high-density multi-distance NIRS [23] showed that the use of multi-distance SDS NIRS improved accuracy by 5.2% compared to that of single SDS NIRS. Accordingly, it is inferred that our study results could be improved with a technique against the extracerebral contamination. Our further study will adopt the use of an additional short SDS of 0.5 cm [31], which is a proper option to filter out those extracerebral contamination in real time.

Conclusions
We have successfully developed a hand rehabilitation system integrating a reliable, commercial NIRS system with our previously developed robotic hand orthosis, in order to help post-stroke patients participate in effective rehabilitation therapy, such as AMT. Moreover, a series of evaluations was conducted to improve the classification performance of the BCI module, the most fundamental element of the system, in aspects of classification accuracy and detection latency. The evaluation results indicated that a shorter feature window is appropriate for a real-time application because it diminishes detection latency into between 2 s and 3 s, although it can preserve only 91.46% of the maximum classification accuracy that is attainable by a longer window. In addition, the results recommend a feature vector composition that preserves channel information and incorporates both [HbO] and [HbR] for better classification performance. Finally, we have to mention that the testing of our system on patients with subcortical stroke is ongoing, showing a positive translatability of the approach and results of this study into clinical cases.
Supplementary Materials: The following are available online at http://www.mdpi.com/2076-3417/9/18/3845/s1, Video S1: System demonstration with a neurologically intact participant performing motor imagery,  . The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.