Hostname: page-component-586b7cd67f-t7czq Total loading time: 0 Render date: 2024-11-29T04:31:34.500Z Has data issue: false hasContentIssue false

Enhancing prosthetic hand control: A synergistic multi-channel electroencephalogram

Published online by Cambridge University Press:  28 November 2024

Pooya Chanu Maibam*
Affiliation:
Embedded Systems and Robotics Lab, Tezpur University, Tezpur, Assam, India
Dingyi Pei
Affiliation:
Vinjamuri Lab, University of Maryland, Baltimore County, Baltimore, MD, USA
Parthan Olikkal
Affiliation:
Vinjamuri Lab, University of Maryland, Baltimore County, Baltimore, MD, USA
Ramana Kumar Vinjamuri
Affiliation:
Vinjamuri Lab, University of Maryland, Baltimore County, Baltimore, MD, USA
Nayan M. Kakoty
Affiliation:
Embedded Systems and Robotics Lab, Tezpur University, Tezpur, Assam, India
*
Corresponding author: Maibam Chanu; Email: [email protected]

Abstract

Electromyogram (EMG) has been a fundamental approach for prosthetic hand control. However it is limited by the functionality of residual muscles and muscle fatigue. Currently, exploring temporal shifts in brain networks and accurately classifying noninvasive electroencephalogram (EEG) for prosthetic hand control remains challenging. In this manuscript, it is hypothesized that the coordinated and synchronized temporal patterns within the brain network, termed as brain synergy, contain valuable information to decode hand movements. 32-channel EEGs were acquired from 10 healthy participants during hand grasp and open. Synergistic spatial distribution pattern and power spectra of brain activity were investigated using independent component analysis of EEG. Out of 32 EEG channels, 15 channels spanning the frontal, central and parietal regions were strategically selected based on the synergy of spatial distribution pattern and power spectrum of independent components. Time-domain and synergistic features were extracted from the selected 15 EEG channels. These features were employed to train a Bayesian optimizer-based support vector machine (SVM). The optimized SVM classifier could achieve an average testing accuracy of 94.39 $ \pm $ .84% using synergistic features. The paired t-test showed that synergistic features yielded significantly higher area under curve values (p < .05) compared to time-domain features in classifying hand movements. The output of the classifier was employed for the control of the prosthetic hand. This synergistic approach for analyzing temporal activities in motor control and control of prosthetic hands have potential contributions to future research. It addresses the limitations of EMG-based approaches and emphasizes the effectiveness of synergy-based control for prostheses.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2024. Published by Cambridge University Press

1. Background and summary

The field of prosthetic hand control has witnessed remarkable advancements in recent years, aiming to enhance the functionality and usability of prosthetic hands by individuals with hand amputation. In research on motor control, electromyogram (EMG) has been a cornerstone enabling users to manipulate prosthetic hands through the detection of residual muscle activity (Campbell et al., Reference Campbell, Phinyomark and Scheme2020). However, EMG-controlled prostheses heavily rely on the availability and functionality of residual muscles. For individuals with amputations or muscle impairments, this reliance poses challenges in achieving consistent and reliable control (Fleming et al., Reference Fleming, Stafford, Huang, Hu, Ferris and Huang2021). Additionally, the overall versatility of prosthetic hand movements in EMG-controlled systems is constrained by the number of accessible muscle and their associated contractions. The complexity of muscle-based interfaces hinders the naturalness and speed of interaction. The users of EMG-controlled prostheses often face a learning curve associated with specific muscle contractions for precise control (Kumar et al., Reference Kumar, Jelfs, Sui and Arjunan2019).

With the advancement of neuroscientific understanding and technology, brain–computer interface (BCI) system with an electroencephalogram (EEG) serves as a means of enhancing the control precision and usability of prosthetic devices (Abiri et al., Reference Abiri, Borhani, Sellers, Jiang and Zhao2019). This will establish a direct channel for communication and control between the human brain and external devices, eliminating the need for reliance on peripheral nerves and muscles. Myokinetic interfaces were introduced for controlling prosthetic hands that rely on detecting muscle movements using surgically implanted sensors in or near the muscles (Gherardini et al., Reference Gherardini, Masiero, Ianniciello and Cipriani2023). This technique presents an alternative to using scalp EEG for prosthetic hand control. However, this interface involves surgically implanting sensors to detect muscle activity posing several challenges. These challenges include the invasiveness of the surgery, risks of infection, rejection, and long-term biocompatibility issues (Salminger et al., Reference Salminger, Sturma, Hofer, Evangelista, Perrin, Bergmeister and Aszmann2019). Additionally, the implanted sensors can cause discomfort for users and require periodic maintenance or replacement resulting in higher costs. However, EEG provides comprehensive motor and cognitive activity insights without the associated risks and complexities of myokinetic interfaces. Bypassing the need for intact muscles, EEG offers an alternative for individuals with nonfunctional or absent muscles, broadening the scope of prosthesis accessibility. EEG allows for the decoding of more intricate neural patterns, potentially offering a greater range of distinguishable hand movements (Hazrati & Erfanian, Reference Hazrati and Erfanian2010). Furthermore, EEG-based control strategies are more anthropomorphic enabling the interaction between the user and the prosthetic hand intuitively. Previous research has highlighted distinctions in EEG patterns associated with various hand movements such as opening, closing, supination, and pronation (Zou et al., Reference Zou, Zhao, He, Zhou and Yin2023). Additionally, differences have been observed among different types of grasps, including pincer, power, and intermediate grasps (Su et al., Reference Su, Wang, Wang, Wang and Hou2023). However, these control strategies using EEG have limitations and challenges in detecting complex hand movements in higher dimensions (Zou et al., Reference Zou, Zhao, He, Zhou and Yin2023; Su et al., Reference Su, Wang, Wang, Wang and Hou2023). EEG acquisition and processing require further engineering to integrate smoothly with prosthetic hands. To address these challenges, synergies are being explored that mimic the intelligent approach employed by the human brain for motor control (Santello et al., Reference Santello, Baud-Bovy and Jorntell2013). Studies have demonstrated the effective identification of synergies in hand movement (Antuvan et al., Reference Antuvan, Bisio, Marini, Yen, Cambria and Masia2016).

Moreover, the complex functional structure of cortical motor areas and the reasons for their numerous interaction during voluntary movement remain not fully elucidated. Neurons within these regions possess distinctive characteristics and engage in diverse stages of movement, from planning to execution through functional integration (Andres & Gerloff, Reference Andres and Gerloff1999; Tononi et al., Reference Tononi, Edelman and Sporns1998). Hand synergy correlation in kinematic movement (Pei et al., Reference Pei, Burns, Chandramouli and Vinjamuri2020), hand movement classification (Erdoĝan et al., Reference Erdoĝan, Özsarfati, Dilek, Kadak, Hanoĝlu and Akın2019), effects of beta rebound, and alpha-coherence during execution of hand movements (Wang et al., Reference Wang, Chen, Yang and Sawan2022; Formaggio et al., Reference Formaggio, Storti, Boscolo Galazzo, Gandolfi, Geroin, Smania and Manganotti2015) have been reported involving EEG channels spanning over central, frontal, occipital, and parietal regions. Studying the synchronization of temporal patterns of the parallel pathways connecting frontal, central, parietal, and occipital regions of the brain provides an understanding of their dynamic roles in voluntary hand movement. Previous studies have shown that the analysis of muscle synergy is employed to investigate the coordination of signals associated with muscle activity (Santello et al., Reference Santello, Bianchi, Gabiccini, Ricciardi, Salvietti, Prattichizzo and Bicchi2016; Li et al., Reference Li, Li, Li, Liu and Qiao2021). The muscle synergy theory is grounded in the concept that the highly redundant musculoskeletal system requires a mechanism to streamline the degrees of freedom in motor control. This mechanism enables the execution of movements. The relevance of muscle synergy analysis becomes particularly evident in studies focusing on motor impairment. Currently, there is no method that ensures consistent categorization of distinct hand movements based on synergy in the brain using noninvasive EEG.

In this study, it is hypothesised that coordinating temporal patterns within the different regions of the brain, referred to as brain synergy, contain sufficient information for decoding various hand movements for prosthetic hand control. 32-channels EEG were recorded from 10 healthy participants during hand grasp and open tasks. Movement-related synergistic brain activity was analyzed using coherence of spatial power distribution pattern and power spectral density (PSD) of independent components of EEG. From the 32 channels, 15 channels located in the frontal, central, and parietal regions were selected using independent component analysis (ICA) based spatial power distribution and power spectral analysis to accurately decode the neural activity associated with hand movements. Time-domain and synergistic features (coherence of spatial power distribution and power spectral) were extracted from the selected 15 channels. These features were used to train a support vector machine (SVM) classifier. The SVM classifier was optimized using a Bayesian optimizer. The optimized SVM classifier achieved an average testing accuracy of 94.39 $ \pm $ .84% across 10 participants using synergistic features. As a proof of concept, the output of the classifier was utilized to control a prosthetic hand. This synergistic method for analyzing temporal neural activities in motor control and controlling prosthetic devices offers promising contributions to future research in synergy-based prosthetic control.

2. Materials and methods

2.1. Human participants

Studies have shown that motor intentions and signals can still be detected after limb loss (Reilly et al., Reference Reilly, Mercier, Schieber and Sirigu2006). Individuals with hand amputation sometimes experience sensations where they perceive the presence of the lost limb (Wijk & Carlsson, Reference Wijk and Carlsson2015). These sensations are generated by the same motor cortical areas that were active before limb loss and thereby indicate the capability to generate motor signals. Furthermore, studies have shown that the motor intentions and signals remain the same for both healthy and individuals with limb amputations (Bruurmijn et al., Reference Bruurmijn, Pereboom, Vansteensel, Raemaekers and Ramsey2017; Chen et al., Reference Chen, Yao, Kuiken and Dewald2013). Following these, our study involved 10 right-handed healthy participants, consisting of four males and six females, with an average age of $ 23.0\pm 3.1 $ years. All procedures were conducted in compliance with relevant standards and regulations. The experimental protocol received approval from the Institutional Review Board at the Stevens Institute of Technology under Protocol No. 2015-022(20-AR6). Prior to the experiment, all participants provided informed written consent, demonstrating their understanding and agreement to participate in the study.

2.2. Task performance and data acquisition

In the experimental setup for data acquisition, an EEG cap (g.GAMMA cap from g.tec, Schieldberg, Austria), an EEG electrode set known as g.Ladybird, and BCI2000 EEG acquisition system (Schalk et al., Reference Schalk, McFarland, Hinterberger, Birbaumer and Wolpaw2004) were employed. Additionally, a computational unit and a cylindrical water bottle were integral parts of the setup. The schematic representation of the experimental methodology employed in this study is illustrated in Figure 1. To evaluate the SVM classifier used for the classification of hand movement, a prosthetic hand (Kakoty et al., Reference Kakoty, Gohain, Saikia, Kalita and Borah2022) controlled using the classification results was employed. The EEG was acquired with 10/20 configuration 32-active high-density electrodes covering central, frontal, parietal, and occipital areas and additional intermediate positions (In1, In2, In3, In4, In5, In6, In7, and In8). The reference electrode was positioned on the earlobe, either at the left (A1) or right (A2), with the ground electrode on the nasion (NZ).

Figure 1. The schematic of the proposed method consists of data acquisition, computational unit, and interfacing with the prosthetic hand for the control of hand open and close tasks.

Following the EEG cap setup, participants were seated in a relaxed position placing their dominant hand palm-down on the table. A water bottle was placed 40 cm away from the body midline of the participant to grasp comfortably as presented in Figure 2(a). The experimental timeline, illustrated in Figure 2(b), guided the participants through the task. Upon hearing an initiation signal from the computational unit, the participants were directed to grasp the cylindrical water bottle and release it when a subsequent signal was heard. The recording duration was set at 4 s, and each participant performed the experiment 30 times. Subjects were advised to minimize blinking and swallowing to mitigate potential artifacts.

Figure 2. Experimental setup and timeline. (a) The illustration of a participant grasping an object during EEG acquisition (image has been adapted from (Pei et al., Reference Pei, Patel, Burns, Chandramouli and Vinjamuri2019)). (b) Experimental timeline for participants to perform the grasping and opening tasks.

2.3. Data and statistical analysis

2.3.1. Data preprocessing

The EEG was recorded continuously with electrode impedance maintained below 10 kOhms and 256 Hz sampling rate, capturing the temporal dynamics of brain activity at a high resolution. To enhance data quality, trials containing artifacts due to eye blinking or swallowing were rejected. Furthermore, the raw EEG data underwent a fourth order Butterworth bandpass filter from a frequency range between .53 and 60 Hz. The utilization of a Butterworth filter (Subasi & Gursoy, Reference Subasi and Gursoy2010), helped to achieve a smooth and effective filtering of the EEG signals, contributing to the precision of subsequent analyses. This specific frequency range was chosen to focus on the relevant neural activity while attenuating unwanted noise and potential artifacts. ICA was applied to the 32 channels of filtered EEG data to separate them into independent components (ICs) in EEGLAB v 2022.1 (Lee et al., Reference Lee, Girolami and Sejnowski1999; Delorme & Makeig, Reference Delorme and Makeig2004). The spatial power distribution pattern associated with each ICs are examined. Furthermore, channels exhibiting task-relevant patterns of PSD across the ICs are also identified. Based on the coherence of both the spatial power distribution patterns and the PSD of ICs, 15 channels corresponding to these ICs are selected. The acquired 32-channel EEG were narrowed down to a subset of 15 channels for the classification of hand grasp and open tasks. These selected 15 channels span the frontal, central, and parietal regions of the human brain. The spatial distribution of the selected EEG channels across the scalp is depicted in Figure 3. This strategic placement of electrodes ensures a representative sampling of neural activity, effectively covering key regions that are directly influenced during the initiation and execution of hand movement indicating brain synergy (Erdoĝan et al., Reference Erdoĝan, Özsarfati, Dilek, Kadak, Hanoĝlu and Akın2019; Teplan, Reference Teplan2002; Wang et al., Reference Wang, Chen, Yang and Sawan2022; Pei et al., Reference Pei, Burns, Chandramouli and Vinjamuri2020).

Figure 3. The spatial distribution of the selected 15 EEG channels across the scalp for this study. The specific locations of these channels are critical for the accurate analysis of EEG, ensuring a representative sampling of neural activity from frontal, central, parietal, and occipital areas.

2.3.2. Feature extraction in the time domain

Eleven time-domain features were extracted from the filtered EEG using a window size of .5 s with an overlap of .3 s (Meier et al., Reference Meier, Dittrich, Schulze-Bonhage and Aertsen2008). The extracted time-domain features are minimum value (min), maximum value (max), median amplitude (med), mean amplitude ( $ \mu $ ), standard deviation ( $ \sigma $ ), variance (var), waveform length (WL), mean absolute value (MAV), root mean square ( $ {X}_{rms} $ ), skewness (skew), and kurtosis (kur). These features encompass crucial aspects of EEG signals, including amplitude, regularity, and synchronization (Emigdio et al., Reference Emigdio, Trujillo, Sotelo, Legrand and Coria2016; Olcay & Karaçalı, Reference Olcay and Karaçalı2019). A 9,000 $ \times $ 11 feature matrix was created to classify hand movements. Each row in the matrix corresponds to 10 individuals conducting hand open and grasp operations. Each subject performed the task 30 times, resulting in 600 trials. The matrix includes data from 15 EEG channels, resulting in 9,000 data points. Each column represents the 11 extracted features across 15 selected EEG channels. The details of the extracted time domain features were reported in (Pooya Chanu et al., Reference Pooya Chanu, Pei, Olikkal, Vinjamuri and Kakoty2023).

2.3.3. Computation of power spectral density

EEG signals collected during grasp and open tasks were segmented into specific time windows. These time windows were .5–1 s, 1.5–2 s, 2.5–3 s, and 3.5–4 s. ICA was employed on the segmented data to decompose the EEG data into ICs where each component represented a distinct neural source. After extracting ICs using ICA, the PSD of each IC was calculated by applying the Fourier transform using Welch’s method (Zhao & He, Reference Zhao and He2013). PSD values were computed for each frequency band, including Broadband (4–100 Hz), theta (3–8 Hz), alpha (8–12 Hz), beta (13–35 Hz), and gamma (36–100 Hz). The frequency range for each band was determined based on previous literature cited (Ali et al., Reference Ali, Afridi, Soomro, Khan, Khan and Chowdhry2022; Van Albada & Robinson, Reference Van Albada and Robinson2013). The coherence of spatial power distribution pattern and PSD during hand grasp and opening were estimated from the ICs. The analysis resulted in 54,000 synergistic features per window, considering all frequency bands across 15 EEG channels.

2.4. Classification

The recorded EEG data results in a high-dimensional feature space. The versatility of kernel selection in SVM allows for nonlinear mapping of input features to higher dimensional spaces, making it suitable for EEG data. The SVM operates by transforming the training dataset through mapping to a nonlinear vector space of high dimensions, accommodating the complexity of EEG data (Lotte et al., Reference Lotte, Congedo, Lécuyer, Lamarche and Arnaldi2007). Additionally, SVM is known for its robustness to overfitting which is a crucial consideration in EEG analysis (Hosseini et al., Reference Hosseini, Aminitabar and Shalchyan2023; Thanigaivelu et al., Reference Thanigaivelu, Sridhar and Sulthana2023; Martínez-Ramón et al., Reference Martínez-Ramón, Rojo-Alvarez, Camps-Valls, Muñoz-Marí, Soria-Olivas and Figueiras-Vidal2006). Based on these, the SVM classifier was chosen as the model for classification in this study.

For the classification training phase, two sets of features were utilized: time-domain features and synergistic features. Initially, the SVM classifier was trained using a combined matrix of time-domain features, which encompassed class labels for hand grasp and open states. The trials of each participant were divided into 24 folds to ensure robust training and evaluation. Subsequently, the synergistic features were employed for classification training using the SVM classifier. This involved a similar methodology to the time-domain feature training incorporating synergistic information.

To further enhance the performance of the SVM classifier, hyperparameter optimization was performed. Given the significant impact of internal parameters on SVM effectiveness, Bayesian optimizer (Luo et al., Reference Luo, Gao, Zhu, Wang, Lu and Wang2020) was employed for optimization of SVM. The objective function for training the model was defined as the mean square error. The Bayesian optimizer iteratively adjusted the hyperparameters within specified search ranges to minimize error in cross-validation. This approach ensured that the SVM classifier was fine-tuned for optimal performance in classifying EEG data.

Figure 4 depicts the average training classification accuracy of linear, cubic, Gaussian, and quadratic kernel functions, both before and after optimization with standard deviation. It utilizes time-domain and synergistic features from 10 participants. Averaging the training set, the minimized classification error was computed considering the input-to-output potential mappings. The adaptation of optimal values for internal parameters underscores the significance of hyper-parameters in minimizing prediction error. The highest accuracy of classification during training was achieved using gaussian kernel with both time-domain and synergistic features.

Figure 4. Comparative plot of classification accuracy using different kernel functions pre-optimization and post-optimization across 10 participants.

3. Results and discussion

3.1. SVM optimization

Figure 5(a and b) depicts the SVM hyperparameter optimization using Bayesian optimization for time domain synergistic features-based classification. The green surface illustrates the performance at different hyperparameter points, while the sampled parameter space is in blue points. The estimated objective function value corresponds to the misclassification rate. The hyperparameter values in the optimized SVM that resulted in the lowest misclassification error during the training phase were identified. The optimized values of the SVM classifier hyperparameter were achieved with Gaussian kernel, kernel scale of 341.21 and box constrain level of 5.93. The optimized values for synergistic features with Gaussian kernel are a kernel scale of .87 and a box constraint level of 271.80. The optimized misclassification rate with time-domain and synergistic features are $ 2.2\times {10}^{-1} $ and $ 3.3\times {10}^{-1} $ respectively during training of the optimized SVM. This optimized configuration enhances the performance of the SVM classifier, considering both time-domain and synergistic features.

Figure 5. The SVM hyper-parameters optimization using a Bayesian optimization approach. (A) The optimized hyper-parameters were obtained for time domain feature-based classification. (B) The optimized hyper-parameters were obtained for synergistic feature-based classification.

3.2. Spatial distribution of power

Studies have shown that investigation of the spatial distribution of power in different regions of the brain demonstrates the specific neural networks engaged during hand movements (Andres & Gerloff, Reference Andres and Gerloff1999; Tononi et al., Reference Tononi, Edelman and Sporns1998). The identification of task-specific activation patterns during hand grasp and open tasks contributes to characterizing distinct motor control processes. The spatial distribution of power, expressed in decibels (dB), over time (in s) during hand grasp and open across 10 participants is depicted in Figure 6. A dynamic changes in the power levels and their spatial patterns was observed during hand grasp and opening. However, no distinct difference was observed during the initiation of hand grasp and execution of hand opening. During the initiation of hand opening at 2.5 s, the increase in power in F4, FC4, FC2, CP1, and P3 shows the formulation of motor commands necessary to perform the hand open task and continuous sensory feedback in refining the movement.

Figure 6. Spatial distribution of power (in decibels (dB) with time (in s) during hand grasp and open states across 10 participants.

3.3. Temporal neural coordination

The graphs depicting inter-trial coherence (ITC) with a significance level of .01 during hand grasp and open for the EEG channels located in the frontal, central, and parietal regions are shown in Figures 79. The x-axis represents time in ms, while the y-axis denotes the coherence values, reflecting the degree of synchronization among the brain regions. The scale on the right of the ITC plots ranging from green to red indicates the magnitude of coherence and the blue squares present in the graph indicate no significant coherence. Furthermore, Figures 79 show the event-related potentials (ERP) and frequency power changes over time during hand grasp and open. The ITC analysis illustrates the temporal dynamics of neural synchronization during the hand grasp and open task providing insights into the coordinated activity across different brain regions. This significance level indicates a high level of confidence in the coherence values in which the observed synchronization patterns are unlikely due to random chance. These results highlight the temporal dynamics of synergy between different EEG channels during the motor tasks for visual comparison of coherence patterns. In Figure 7, the ITC values indicate consistent engagement in synchronized neural activity in the frontal region during hand grasp. Additionally, distinct peaks in the ERP waveforms show active neural processing in the frontal region for planning and decision-making involved in initiating and executing hand movements. In Figure 8, the ITC values for the central region show strong synchronization and neural activation during hand movements. The high ITC values observed in Figure 8 reflect significant engagement of the motor cortex, with the ERP waveforms indicating the direct involvement motor cortex in controlling and executing grasp tasks in a time-locked manner. In Figure 9, the ITC values for the parietal region exhibit significant synchronization patterns with high coherence. The ERP alterations in this region indicate its role in guiding motor responses during hand grasp and open tasks with sensory feedback. These findings collectively illustrate the coordinated neural activity known as synergy across different brain regions during hand movements with ITC and ERP analyses.

Figure 7. ITC during hand grasping and opening in the EEG channel located at the frontal region with the significance values of .01.

Figure 8. ITC during hand grasping and opening in the EEG channel located at the central region with the significance values of .01.

Figure 9. ITC during hand grasping and opening in the EEG channel located at the parietal region with the significance values of .01.

The frequency versus time plot for EEG during hand grasp and open tasks is presented in Figure 10. This result shows the temporal variations in neural activity, highlighting a surge in power within 3–8 Hz frequency range during the hand grasping phase. This elevation in power in the theta frequency band signifies intensified neural engagement associated with motor planning, execution, and sensorimotor integration of voluntary movements. The observed increase in power in theta range (3–8 Hz) reflects the complex neural processes involved in hand grasping. The frequency versus time plot not only captures variations in power but also provides insights into the temporal dynamics and functional connectivity patterns contributing to the successful execution of the hand grasp task. However, during the hand-opening phase, distinct neural activation patterns are not observed with no power surge. The neural activation diminishes after 2 s of hand grasp. The distinct surge in power within the theta frequency range (3–8 Hz) during hand grasp and no power surge during hand opening shows the more intensive neural engagement due to the need for precise motor control and sensorimotor feedback integration. This neural activity can be a critical feature for differentiating hand grasp and open task classification. These results underscore the importance of considering both the temporal dynamics and frequency-specific power changes in neural activity for effective synergistic classification of hand movements. Studies reported that during hand grasp, the complex network of cortical areas is involved in planning and execution leading to a significant rise in theta activity (Goldenkoff et al., Reference Goldenkoff, Deluisi, Destiny, Lee, Michon, Brissenden and Vesia2023; Iturrate et al., Reference Iturrate, Chavarriaga, Pereira, Zhang, Corbet, Leeb and Millán2018). However, during hand opening, the less complex motor planning and execution were involved compared to hand grasp (Goldenkoff et al., Reference Goldenkoff, Deluisi, Destiny, Lee, Michon, Brissenden and Vesia2023). Therefore, there may not be a significant rise in theta activity as the demands on the cortical grasping network and associated neural populations are less.

Figure 10. Spatial distribution of power (in decibels (dB) with time (in s)) during hand grasp and open across 10 participants.

3.4. Evaluation of classifier performance

The analysis of the area under the curve (AUC) for the classifier trained using the time domain and synergistic features presented in Figure 11 provides valuable insights into the SVM classifier in the context of distinguishing hand movement classes. The AUC value for synergistic features is .945 which is higher than the AUC for time-domain features which .926. The resulted values of AUC value indicates that the classifier is optimal for classifying the hand movement dataset as reported in (Metz, Reference Metz1978). This distinction suggests that synergistic features play a more effective role in enhancing the accuracy of a classifier to discern subtle variations in hand movement patterns.

Figure 11. AUC values indicate the performance of 10 participants using time-domain and synergistic features with p < .005. Grey dots = individual performance AUC values, thick dotted line = group mean.

Z-score normalization was applied to both the synergistic and time-domain features to standardize each feature. A paired t-test was performed to compare the performance of classifiers trained with the normalized time-domain features versus synergistic features. This test evaluates the significance of the difference in AUC values between the two classifiers. The results revealed a p-value < .05, indicating that the difference in AUC values is statistically significant. This finding shows that synergistic features offer a more accurate and reliable basis for classifying hand movement patterns compared to time-domain features. The reliance of the SVM classifier on synergistic features highlights the significance of frequency-specific information in capturing the neural signals associated with distinct motor activities. This finding has practical implications for the design and optimization of classifiers in applications for synergy-based prosthetic control.

For the classification of hand grasp and open using the SVM classifier, 80% of the dataset was randomly selected for training, while the remaining 20% was used for testing. This division of the dataset ensures a comprehensive evaluation of the performance of the classifier. This leveraged a significant portion of the data for training while reserving a sufficient sample for testing to validate the accuracy of the classifier. The classification accuracy using time-domain and synergistic features of 10 participants is presented in Figure 12. During the training, the classifier exhibits an accuracy of 96.8 $ \pm $ .98% with time-domain features. During testing of the classifier, an average accuracy of 93.4 $ \pm $ 1.16% was obtained across 10 participants. The lowest accuracy of classification was 92.1% in the S10 and the highest accuracy of 95.4% was observed in S8.

Figure 12. Classification accuracy of 10 participants with time domain features and synergistic features.

In the second set of classifications with synergistic features, the classifier maintained a high accuracy of 97.1 $ \pm $ .9134% during the training and achieved an accuracy of 94.39 $ \pm $ .8434% during testing. The classification accuracy varied among subjects, with the highest accuracy recorded for S10 at 96.47% and the lowest accuracy observed in S2 at 91.24%. These findings underscore the robustness of the classifier across different sets of data, indicating its adaptability and effectiveness in consistently distinguishing between hand grasp and open task.

3.5. Proof-of-concept

The optimized SVM classifier used for the classification of hand grasp and open was deployed to control prosthetic hand grasping operations. A seamless interface was established between the output of the SVM classifier and a prosthetic hand for the visualisation of classification performance. The output of the SVM classifier framed in transistor-transistor logic was directed to an 8-bit microcontroller. The controller translates the output of the SVM classifier into the corresponding actuating signal. The actuating signal was transmitted to the actuators on the prosthetic hand through a current buffer circuit. The control circuit was powered by a standalone 9-volt lithium polymer battery.

This controlled operation effectively demonstrated accurate classification performance by the optimized SVM. The integration of selected EEG channels, time-domain features and synergistic features underscored the versatility and adaptability of the SVM classifier. Figure 13 visually depicts the grasp and open task of the prosthetic hand, showcasing the successful implementation of the optimized SVM classifier in real-time prosthetic hand control. This innovative proof of concept not only validates the efficacy of the synergistic features-based SVM classification but also highlights its potential for practical applications in neuroprosthetics, offering a promising avenue for enhancing the lives of individuals with limb loss. However, the integration of EEG with prosthetic hand control presents significant challenges during both data acquisition and control of prosthetic devices. Noise from muscle activity and environmental factors distort neural signals, necessitating the use of artifact removal techniques. Variability in classification accuracy is influenced by fluctuations in EEG signals and the user, suggests increasing the number of participants for improved classification outcomes. While the integration of SVM classifier outputs with prosthetic hand actuators has shown promise in translating EEG into hand movement commands, residual noise in EEG and classification variability remain as a limiting factor. To overcome these challenges, advanced signal processing techniques, algorithm optimization, and the development of compact EEG systems will be essential in fully realizing the potential of EEG-based prosthetic control in practical applications.

Figure 13. Prosthetic hand executing grasping and opening task with the optimized SVM classifier using synergistic features.

4. Conclusions

The study demonstrated the complex neural dynamics associated with hand grasp and open tasks using EEG-based analysis. The spatial distribution of power across different brain regions during these tasks observed dynamic changes in power levels and spatial patterns. Inter-trial coherence analysis further reveals consistent synchronization patterns in the coordination of neural activity. Examining ERP and frequency dynamics during hand grasp and open tasks demonstrated temporal variations in neural activity. Notably, the alteration in power within the theta frequency range during hand grasp indicates neural engagement associated with motor planning and execution. The temporal dynamics and functional connectivity patterns revealed through these analyses contribute to our understanding of the synergy between different EEG channels during motor tasks.

The classifier performance analysis underscores the robustness of an optimized SVM classifier in distinguishing hand movement classes. Both time-domain and synergistic features contribute to high accuracy during training and robust generalization during testing. The higher AUC for synergistic features highlights their effectiveness in enhancing classifier accuracy, emphasizing the significance of frequency-specific information in capturing neural signals for hand movement. The high accuracy during training and robust generalization during testing underscores the potential of this approach for real-time applications in neuroprosthetics. These findings contribute to advancing the field, offering a promising avenue for enhancing the quality of life for individuals with limb loss through brain synergy-based control of prosthesis.

Data availability statement

Not applicable.

Authorship contribution

Conceptualization, M.P.C., D.P., P.O., R.V., N.M.K.; Data curation, M.P.C., D.P.; Data acquisition, D.P.; Formal analysis, M.P.C.; Funding acquisition, R.V., N.M.K.; Project administration, R.V., N.M.K.; Resources, R.V., N.M.K.; Supervision, R.V., N.M.K.; Methodology, M.P.C., R.V., N.M.K.; Investigation, M.P.C.; Software M.P.C.; Visualization M.P.C.; Writing-original draft M.P.C.; Writing-review and editing M.P.C., D.P., P.O., R.V., N.M.K. All authors have reviewed and approved the final version of the manuscript for publication.

Funding statement

The research received funding from the Innovation Hub Foundation for Cobotics, IIT Delhi, Department of Science and Technology (DST), Government of India (project number: GP/2021/RR/017), BITS BioCyTiH Foundation project number BBF/BITS/FY2023–24/ART-103 and the National Science Foundation (NSF) CAREER Award, grant number HCC–2053498, and NSF Planning IUCRC Award, grant number 20422.

Competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Ethical standard

The research was carried out following the approved Institutional Review Board protocol at the Stevens Institute of Technology under Protocol No. 2015-022(20-AR6), with approval granted in September 2020. All participants provided informed consent before participating in the study.

References

Abiri, R., Borhani, S., Sellers, E. W., Jiang, Y., and Zhao, X. (2019) A comprehensive review of EEG-based brain–computer interface paradigms. Journal of Neural Engineering, 16(1) 011001.CrossRefGoogle ScholarPubMed
Ali, A., Afridi, R., Soomro, T. A., Khan, S. A., Khan, M. Y. A., and Chowdhry, B. S. (2022) A single-channel wireless EEG headset enabled neural activities analysis for mental healthcare applications. Wireless Personal Communications, 125(4) 36993713.CrossRefGoogle ScholarPubMed
Andres, F. G., and Gerloff, C. (1999) Coherence of sequential movements and motor learning. Journal of Clinical Neurophysiology, 16(6), 520.CrossRefGoogle ScholarPubMed
Antuvan, C.W., Bisio, F., Marini, F., Yen, S.-C., Cambria, E. and Masia, L. (2016) Role of muscle synergies in real-time classification of upper limb motions using extreme learning machines. Journal of Neuroengineering and Rehabilitation, 13(1), 115.CrossRefGoogle ScholarPubMed
Bruurmijn, M. L., Pereboom, I. P., Vansteensel, M. J., Raemaekers, M. A., and Ramsey, N. F. (2017) Preservation of hand movement representation in the sensorimotor areas of amputees. Brain, 140(12), 31663178.CrossRefGoogle ScholarPubMed
Campbell, E. and Phinyomark, A. and Scheme, E. (2020) Current trends and confounding factors in myoelectric control: limb position and contraction intensity. Sensors 20(6), 1613.CrossRefGoogle ScholarPubMed
Chen, A., Yao, J., Kuiken, T., and Dewald, J. P. (2013) Cortical motor activity and reorganization following upper-limb amputation and subsequent targeted reinnervation. NeuroImage: Clinical, 3, 498506.CrossRefGoogle ScholarPubMed
Delorme, A., and Makeig, S. (2004) EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. Journal of Neuroscience Methods, 134(1), 921.CrossRefGoogle ScholarPubMed
Emigdio, Z., Trujillo, L., Sotelo, A., Legrand, P., and Coria, L. N. (2016) Regularity and Matching Pursuit feature extraction for the detection of epileptic seizures. Journal of Neuroscience Methods, 266, 107125.Google Scholar
Erdoĝan, S. B., Özsarfati, E., Dilek, B., Kadak, K. S., Hanoĝlu, L., and Akın, A. (2019) Classification of motor imagery and execution signals with population-level feature sets: implications for probe design in fNIRS based BCI. Journal of Neural Engineering, 16(2), 026029.CrossRefGoogle ScholarPubMed
Fleming, A., Stafford, N., Huang, S., Hu, X., Ferris, D. P., and Huang, H. H. (2021) Myoelectric control of robotic lower limb prostheses: a review of electromyography interfaces, control paradigms, challenges and future directions. Journal of Neural Engineering, 18 (4), 041004.CrossRefGoogle ScholarPubMed
Formaggio, E., Storti, S. F., Boscolo Galazzo, I., Gandolfi, M., Geroin, C., Smania, N., … and Manganotti, P (2015) Time–frequency modulation of ERD and EEG coherence in robot-assisted hand performance. Brain Topography, 28, 352363.CrossRefGoogle ScholarPubMed
Gherardini, M., Masiero, F., Ianniciello, V., and Cipriani, C. (2023) The myokinetic interface: Implanting permanent magnets to restore the sensory-motor control loop in amputees. Current Opinion in Biomedical Engineering, 27, 100460.CrossRefGoogle Scholar
Goldenkoff, E. R., Deluisi, J. A., Destiny, D. P., Lee, T. G., Michon, K. J., Brissenden, J. A., … and Vesia, M. (2023) The behavioral and neural effects of parietal theta burst stimulation on the grasp network are stronger during a grasping task than at rest. Frontiers in Neuroscience, 17, 1198222.CrossRefGoogle ScholarPubMed
Hazrati, M. K., and Erfanian, A. (2010) An online EEG-based brain–computer interface for controlling hand grasp using an adaptive probabilistic neural network. Medical Engineering and Physics, 32(7), 730739.CrossRefGoogle ScholarPubMed
Hosseini, S. M., Aminitabar, A. H., and Shalchyan, V. (2023) Investigating the application of graph theory features in hand movement directions decoding using EEG signals. Neuroscience Research, 194, 2435.CrossRefGoogle ScholarPubMed
Iturrate, I., Chavarriaga, R., Pereira, M., Zhang, H., Corbet, T., Leeb, R., and Millán, J. D. R. (2018) Human EEG reveals distinct neural correlates of power and precision grasping types. NeuroImage, 181, 635644.CrossRefGoogle ScholarPubMed
Kakoty, N. M., Gohain, L., Saikia, J. B., Kalita, A. J., and Borah, S. (2022) Real-time EMG based prosthetic hand controller realizing neuromuscular constraint International Journal of Intelligent Robotics and Applications, 6(3), 530542.CrossRefGoogle Scholar
Kumar, D. K., Jelfs, B., Sui, X., and Arjunan, S. P. (2019) Prosthetic hand control: a multidisciplinary review to identify strengths, shortcomings, and the future. Biomedical Signal Processing and Control, 53, 101588.CrossRefGoogle Scholar
Lee, T. W., Girolami, M., and Sejnowski, T. J. (1999) Independent component analysis using an extended infomax algorithm for mixed subgaussian and supergaussian sources. Neural Computation, 11(2), 417441.CrossRefGoogle ScholarPubMed
Li, G., Li, Z., Li, J., Liu, Y., and Qiao, H. (2021) Muscle-synergy-based planning and neural-adaptive control for a prosthetic arm. IEEE Transactions on Artificial Intelligence, 2(5), 424436.CrossRefGoogle Scholar
Lotte, F., Congedo, M., Lécuyer, A., Lamarche, F., and Arnaldi, B. (2007) A review of classification algorithms for EEG-based brain–computer interfaces. Journal of Neural Engineering, 4(2), R1.CrossRefGoogle ScholarPubMed
Luo, J., Gao, X., Zhu, X., Wang, B., Lu, N., and Wang, J. (2020) Motor imagery EEG classification based on ensemble support vector learning. Computer Methods and Programs in Biomedicine, 193, 105464.CrossRefGoogle ScholarPubMed
Martínez-Ramón, M., Rojo-Alvarez, J. L., Camps-Valls, G., Muñoz-Marí, J., Soria-Olivas, E., and Figueiras-Vidal, A. R. (2006) Support vector machines for nonlinear kernel ARMA system identification. IEEE Transactions on Neural Networks, 17(6), 16171622.CrossRefGoogle ScholarPubMed
Meier, R., Dittrich, H., Schulze-Bonhage, A., and Aertsen, A.(2008) Detecting epileptic seizures in long-term human EEG: a new approach to automatic online and real-time detection and classification of polymorphic seizure patterns. Journal of Clinical Neurophysiology, 25(3), 119131.CrossRefGoogle Scholar
Metz, C. E. (1978) Basic principles of ROC analysis Seminars in Nuclear Medicine, 8(4), 283298.CrossRefGoogle ScholarPubMed
Olcay, B. O., and Karaçalı, B. (2019) Evaluation of synchronization measures for capturing the lagged synchronization between EEG channels: a cognitive task recognition approach. Computers in Biology and Medicine, 114, 103441.CrossRefGoogle ScholarPubMed
Pei, D., Burns, M., Chandramouli, R., and Vinjamuri, R. (2020) Neural decoding of upper limb movements using electroencephalography. Brain–Computer Interface Research: A State-of-the-Art Summary, 8, 2533.CrossRefGoogle Scholar
Pei, D., Patel, V., Burns, M., Chandramouli, R. and Vinjamuri, R. (2019) Neural decoding of synergy-based hand movements using electroencephalography. IEEE Access, 7, 1815518163.CrossRefGoogle Scholar
Pooya Chanu, M., Pei, D., Olikkal, P., Vinjamuri, R. K., and Kakoty, N. M. (2023) Electroencephalogram based Control of Prosthetic Hand using Optimizable Support Vector Machine. In Proceedings of the 6th International Conference on Advances in Robotics, 16.Google Scholar
Reilly, K. T., Mercier, C., Schieber, M. H., and Sirigu, A. (2006) Persistent hand motor commands in the amputees’ brain. Brain, 129(8), 22112223.CrossRefGoogle ScholarPubMed
Salminger, S., Sturma, A., Hofer, C., Evangelista, M., Perrin, M., Bergmeister, K. D., … and Aszmann, O. C. (2019) Long-term implant of intramuscular sensors and nerve transfers for wireless control of robotic arms in above-elbow amputees. Science Robotics, 4(32), eaaw6306.CrossRefGoogle ScholarPubMed
Santello, M., Baud-Bovy, G. and Jorntell, H. (2013) Neural bases of hand synergies. Frontiers in Computational Neuroscience, 7, 23.CrossRefGoogle ScholarPubMed
Santello, M., Bianchi, M., Gabiccini, M., Ricciardi, E., Salvietti, G., Prattichizzo, D., … and Bicchi, A. (2016) Hand synergies: integration of robotics and neuroscience for understanding the control of biological and artificial hands. Physics of Life Reviews, 17, 123.CrossRefGoogle ScholarPubMed
Schalk, G., McFarland, D. J., Hinterberger, T., Birbaumer, N., and Wolpaw, J. R.(2004) BCI2000: a general-purpose brain-computer interface (BCI) system IEEE Transactions on Biomedical Engineering, 51(6), 10341043.CrossRefGoogle ScholarPubMed
Su, J., Wang, J., Wang, W., Wang, Y. and Hou, Zeng-Guang (2023) Enhanced motor imagery based brain–computer interface via vibration stimulation and robotic glove for post-stroke rehabilitation. International Conference on Neural Information Processing, 1963, 326337.Google Scholar
Subasi, A. and Gursoy, M. I. (2010) EEG signal classification using PCA, ICA, LDA and support vector machines. Expert Systems with Applications, 37(12), 86598666.CrossRefGoogle Scholar
Teplan, M. (2002)Fundamentals of EEG measurement. Measurement Science Review, 2(2), 111.Google Scholar
Thanigaivelu, P. S., Sridhar, S. S., and Sulthana, S. F. (2023) OISVM: optimal incremental support vector machine-based EEG classification for brain–computer interface model. Cognitive Computation, 15(3), 888903.CrossRefGoogle Scholar
Tononi, G., Edelman, G. M., and Sporns, O. (1998) Complexity and coherency: Integrating information in the brain. Trends in Cognitive Sciences, 2(12), 474484.CrossRefGoogle ScholarPubMed
Van Albada, S. J. and Robinson, P. A. (2013) Relationships between electroencephalographic spectral peaks across frequency bands. Frontiers Media SA, 7, 56.Google ScholarPubMed
Wang, J., Chen, Y. H., Yang, J., and Sawan, M. (2022) Intelligent classification technique of hand motor imagery using EEG Beta rebound follow-up pattern. Biosensors, 12(6), 384.CrossRefGoogle ScholarPubMed
Wijk, U., and Carlsson, I. (2015) Forearm amputees’ views of Prosthesis use and Sensory Feedback. Journal of Hand Therapy, 28(3), 269278.CrossRefGoogle ScholarPubMed
Zhao, L., and He, Y. (2013) Power spectrum estimation of the Welch method based on imagery EEG. Applied Mechanics and Materials, 278, 12601264.CrossRefGoogle Scholar
Zou, R., Zhao, L., He, S., Zhou, X. and Yin, X. (2023) Effect of the period of EEG signals on the decoding of motor information. Physical and Engineering Sciences in Medicine, 47, 112.Google Scholar
Figure 0

Figure 1. The schematic of the proposed method consists of data acquisition, computational unit, and interfacing with the prosthetic hand for the control of hand open and close tasks.

Figure 1

Figure 2. Experimental setup and timeline. (a) The illustration of a participant grasping an object during EEG acquisition (image has been adapted from (Pei et al., 2019)). (b) Experimental timeline for participants to perform the grasping and opening tasks.

Figure 2

Figure 3. The spatial distribution of the selected 15 EEG channels across the scalp for this study. The specific locations of these channels are critical for the accurate analysis of EEG, ensuring a representative sampling of neural activity from frontal, central, parietal, and occipital areas.

Figure 3

Figure 4. Comparative plot of classification accuracy using different kernel functions pre-optimization and post-optimization across 10 participants.

Figure 4

Figure 5. The SVM hyper-parameters optimization using a Bayesian optimization approach. (A) The optimized hyper-parameters were obtained for time domain feature-based classification. (B) The optimized hyper-parameters were obtained for synergistic feature-based classification.

Figure 5

Figure 6. Spatial distribution of power (in decibels (dB) with time (in s) during hand grasp and open states across 10 participants.

Figure 6

Figure 7. ITC during hand grasping and opening in the EEG channel located at the frontal region with the significance values of .01.

Figure 7

Figure 8. ITC during hand grasping and opening in the EEG channel located at the central region with the significance values of .01.

Figure 8

Figure 9. ITC during hand grasping and opening in the EEG channel located at the parietal region with the significance values of .01.

Figure 9

Figure 10. Spatial distribution of power (in decibels (dB) with time (in s)) during hand grasp and open across 10 participants.

Figure 10

Figure 11. AUC values indicate the performance of 10 participants using time-domain and synergistic features with p < .005. Grey dots = individual performance AUC values, thick dotted line = group mean.

Figure 11

Figure 12. Classification accuracy of 10 participants with time domain features and synergistic features.

Figure 12

Figure 13. Prosthetic hand executing grasping and opening task with the optimized SVM classifier using synergistic features.