A novel channel selection method for multiple motion classification using high-density electromyography
© Geng et al.; licensee BioMed Central Ltd. 2014
Received: 6 March 2014
Accepted: 14 July 2014
Published: 25 July 2014
Selecting an appropriate number of surface electromyography (EMG) channels with desired classification performance and determining the optimal placement of EMG electrodes would be necessary and important in practical myoelectric control. In previous studies, several methods such as sequential forward selection (SFS) and Fisher-Markov selector (FMS) have been used to select the appropriate number of EMG channels for a control system. These exiting methods are dependent on either EMG features and/or classification algorithms, which means that when using different channel features or classification algorithm, the selected channels would be changed. In this study, a new method named multi-class common spatial pattern (MCCSP) was proposed for EMG selection in EMG pattern-recognition-based movement classification. Since MCCSP is independent on specific EMG features and classification algorithms, it would be more convenient for channel selection in developing an EMG control system than the exiting methods.
The performance of the proposed MCCSP method in selecting some optimal EMG channels (designated as a subset) was assessed with high-density EMG recordings from twelve mildly-impaired traumatic brain injury (TBI) patients. With the MCCSP method, a subset of EMG channels was selected and then used for motion classification with pattern recognition technique. In order to justify the performance of the MCCSP method against different electrode configurations, features and classification algorithms, two electrode configurations (unipolar and bipolar) as well as two EMG feature sets and two types of pattern recognition classifiers were considered in the study, respectively. And the performance of the proposed MCCSP method was compared with that of two exiting channel selection methods (SFS and FMS) in EMG control system.
The results showed that in comparison with the previously used SFS and FMS methods, the newly proposed MCCSP method had better motion classification performance. Moreover, a fixed combination of the selected EMG channels was obtained when using MCCSP.
The proposed MCCSP method would be a practicable means in channel selection and would facilitate the design of practical myoelectric control systems in the active rehabilitation of mildly-impaired TBI patients and in other rehabilitation applications such as the multifunctional myoelectric prostheses for limb amputees.
Surface electromyography (EMG) is an electrical manifestation of muscle contractions, which has been used as a convenient and proprioceptive control signals for the operations of limb prostheses , assistive exoskeletons , and rehabilitation robots . In recent years, one of the most popular applications of EMG signals would be to use EMG-based pattern recognition (EMG-PR) algorithms to control multifunctional powered prostheses [1, 4–7]. More recently, the EMG-PR also was adopted in the neural rehabilitation of the motor impaired patients who suffer from neurological injuries for development of an active rehabilitation system. The previous studies have proved that the active rehabilitation is promising in enhancing the therapeutic effect  and accelerating the brain plasticity . In light of these findings, the EMG-PR algorithm was proposed to identify the motor intentions of stroke survivors and patients with incomplete cervical spinal cord injury [10, 11]. The results of these previous studies have suggested that EMG-PR algorithms would be feasible and useful in active rehabilitation of patients suffering from neurological injuries.
High density EMG recordings have been increasingly applied in EMG-PR algorithms [12–14] with an attempt to capture more temporal and spatial information about the muscle activities and electrophysiology. With more channels of temporal information that can reflect the activities of some small forearm muscles, a better motion classification performance could be achieved when using high density EMG recordings [13, 14]. And the topographical maps of EMG amplitude can be used to examine the exact electrode locations where the strong myoelectric activity is experienced during a motion task [12–14]. However, the data processing of high-density EMG recordings from a large number of EMG electrodes is computationally expensive, making it impractical in real-time myoelectric control. Additionally, the high dimensional EMG recordings may cause the classifier to over-fit the training data due to the irrelevant or redundant information. Therefore, selecting an appropriate number of EMG channels with desired classification performance and determining the optimal placement locations of electrodes would be necessary in the practical myoelectric control for a potential user.
Currently, there are two commonly used ways to determine the number and the locations of EMG electrodes in the applications of EMG-based motion classification. One direct way is based on the physiologically known anatomical knowledge of skeletal muscles [15–17], in which the electrode sites would be chosen by viewing the intensity, repeatability and consistency of multi-channel EMG signal recordings. For the persons with the intact functions of skeletal muscles, using the muscle-physiology-based way we may quickly and easily determine the appropriate sites of EMG electrodes since their almost consistent muscle anatomy could ensure producing very similar muscle contractions when they actuate a limb activity. However, for the persons with some muscular issues such as the patients with a post-stroke, traumatic brain injury or limb amputation, when doing a limb activity, their muscle contractions would be different from healthy people and may be various between patients since the extents, the positions and the causes of their muscular damages might be inconsistent. Thus it would be hardly to use this kind of clinical method for determining the appropriate electrode sites in patients.
Another way to determine a subset of appropriate EMG channels is based on a certain optimizing criterion. A commonly used method is the sequential forward selection (SFS)-based channel selection algorithm [14, 16], in which the most informative EMG channels were determined by an optimization procedure with an optimizing objective of high motion classification accuracy. Since SFS method requires to repeatedly perform the searching procedure for the optimal channels until obtaining a given number of channels, this classification-accuracy-based method would take a long computation time to get a relatively optimal subset of EMG channels. Moreover, SFS method must rely on a specific set of EMG features and a specific EMG pattern-recognition algorithm to do the selection of the optimal channels for an EMG control application. When the EMG features and/or pattern-recognition algorithm would be changed for some reasons such as trying to use different features or algorithms for higher classification performance, the selected EMG channels might be not optimal yet. Thus the optimization procedure would be executed again to search the optimal channels accordingly. Some feature-dependent algorithms such as the quasi-optimal channel selection method based on partial Kullback-Leibler information  and coefficient-based channel selection method using principal component analysis  also have been proposed and used for channel selection. In addition, another two feature selection methods, Fisher-Markov selector (FMS) based on Fisher criterion  and the minimal-redundancy-maximal-relevance based on mutual information , also fall into the group of feature-dependent channel selection methods. It is obvious that these feature-dependent channel selection methods would rely on the EMG features. Similarly, when changing the set of EMG features such as increasing/decreasing a feature to/from it, the selected channels might not work well again, thus another channel selection procedure would be required to re-determine the appropriate channels accordingly.
Alternatively, the direct channel selection method (i.e., variable selection) which directly works on the raw EMG data also has been proposed to determine a set of appropriate channels for EMG recordings. The Monte Carlo method reported by Nagata et al. for hand motion classification is one of such kinds of the direct channel selection approaches . Unlike the classification-accuracy-based and the feature-dependent channel selection methods, this raw-data-based method is independent of EMG features and classification algorithms. Thus a unique set of optimal EMG channels might be determined from the multi-channel EMG recordings. As a result, when the features and/or classification algorithm are changed, the selected EMG channels would be retained. In this study, we proposed a novel direct channel selection approach named multi-class common spatial pattern (MCCSP) for channel selection in EMG pattern-recognition-based movement classification. Since MCCSP would be independent on EMG features and classification algorithms, it would be more convenient for channel selection in developing an EMG control system than the exiting methods. The performance of the proposed MCCSP method in selecting an optimal set of EMG channels was assessed with high-density EMG recordings from twelve mildly-impaired traumatic brain injured (TBI) patients. Generally, the TBI patients would suffer from the physical and behavioral disorders, so appropriately physical therapy should be employed for the rehabilitation of their physical functions. The EMG-control active rehabilitation systems would be a promising way for TBI patients’ neural function recovery . In addition, the performance of the proposed MCCSP method was compared with that of two exiting channel selection methods, SFS and FMS.
EMG data preprocessing
EMG channel selection
The proposed MCCSP-based channel selection algorithm as well as the SFS and FMS were applied on the two modes of high-density EMG recordings to select an optimal set of EMG channels (1-20), respectively.
MCCSP-based channel selection algorithm considers the EMG data from all the channels simultaneously. To clearly explain its principle and implementation, the principle of two-class common spatial pattern (CSP) algorithm, the CSP-based channel selection algorithm, and the MCCSP-based channel selection algorithm were successively introduced as follows.
where and . With the constraint condition D A + D B = I, the eigenvectors with the largest eigenvalues for D A have the smallest eigenvalues for D B and vice versa. Applying the filter matrix W to the raw signals X would give d output signals Y = W * X that were also called components. The variance of each component is indicated by its corresponding eigenvalue on the principal diagonal of D A , for class B of D B . In the binary classification problem of EEG signals, the first component that has the largest variance for class A and the smallest variance for class B and the last component that has the smallest variance for class A and the largest variance for class B were commonly used .
CSP-based channel selection
With W-1 (inverse matrix of W) that is called spatial patterns, the original signal matrix X can be reconstructed by X = W-1 * Y. The first and last columns of W-1 would be the most important spatial patterns that explain the largest variance of one motion and the smallest variance of the other. If considering the coefficients of spatial patterns as projection weights, the two channels corresponding to the maximal coefficients of spatial pattern vectors may be the channels that are most correlated with the motion specific sources. This method has been used by Wang et al. to select four EEG channels from 118 channels for classifying three motor imaginary tasks and achieved a high classification accuracy .
MCCSP-based channel selection
Accordingly, multiple spatial pattern matrices corresponding to all the included motion classes would be obtained, where c denotes the class label. Then we used these spatial pattern matrices to select the optimal channels for subsequent pattern recognition analysis. In the study, we adopted similar coefficients-based channel selection method as did in the two-class CSP. For each spatial pattern matrix corresponding to a specific motion class, the two channels corresponding to the maximal coefficients of the first and the last spatial pattern vectors were selected, and then all the selected two-channel from each spatial pattern matrix W C -1 were combined together, which were then used as the optimal channels in the subsequent pattern recognition analysis. Note that for the two channels selected from a specific spatial pattern matrix, e.g., W A -1, one or both of which may be overlapped with that obtained from another spatial pattern matrix W B -1. So the selected channels for all movements were ranked in the order of their occurrence frequency from high to low. The number of selected EMG channels might vary for different subjects.
SFS is an iterative searching procedure, in which one optimal channel that produces the highest classification accuracy was firstly selected among all the channels, and then another channel that can achieve the maximum classification accuracy with the combination with the selected channels was added [14, 16]. Suppose we require to choose an optimal subset with n channels from all the 56 channels, the classification computation procedure would be repeatedly implemented by (56 + 56 + 1 - n) * n/2 times.
FMS is essentially a feature selection algorithm, which was proposed to select the globally optimal subset of features from high-dimensional feature space in the spirit of Fisher’s class separation criterion . By using specific kernel functions and the Markov random field optimization techniques, the FMS is capable of selecting features and has efficient computational complexity. In this study, we firstly applied FMS to the EMG feature matrix to get the indices of all the features that were ranked in the order of their importance coefficients, and then the feature indexes were mapped to their corresponding channel indexes. In this way, the optimal EMG channels could be selected.
Feature extraction and motion classification
With a selected subset of EMG signals, a shifting analysis window with a time length of 150 ms and an increment of 100 ms (50 ms overlapping) was used to segment the EMG signals into a series of analysis windows. Then two commonly used time-domain feature sets were extracted from each analysis window, respectively. They were (1) four time-domain features (TD) , mean absolute value (MAV), number of zeros crossings (ZC), number of slope sign changes (SSC), and waveform length (WL), and (2) the six order autoregressive (AR) model coefficients  plus the root mean square (RMS) amplitude of EMG signals (TDAR). For each analysis window, the features extracted from all the selected EMG channels was concatenated to form a feature matrix, which was then fed into a classifier for motion classification.
Two classifier algorithms, the linear discriminant analysis (LDA) and the k-nearest neighbor (KNN), were used in this study. The LDA classifier was used for its merit of algorithm simplicity and low computational cost , and the KNN classifier was used for its tolerance to arbitrary data distribution . In order to find the best feature-classifier combination for each selected EMG subset, each of the two classifiers was separately combined with each of the two feature sets, yielding four different feature-classifier combinations.
The classification accuracy, defined as the percentage of the number of correct classification decisions over the total number of classification decisions, was used to assess the motion classification performance . Five-fold cross validation was performed in motion classification when using the EMG subset selected via MCCSP and FMS for the purpose of evaluating their performance in EMG-PR. The EMG feature matrix extracted from the selected EMG recordings was randomly divided into five subsets with equal length. Four of the five subsets were used as a training set and the remaining one subset was considered as a testing set in each classification process to compute the classification accuracy. The classification accuracies were averaged over the five validation results from the folds. For the EMG subset selected by SFS method, the first half of the feature matrix was used as the training set and its second half was used as the testing set. No cross validation was conducted herein considering the newly added EMG channel was determined according to the maximum classification accuracy during last round repetition.
To assess the statistical difference among the three channel selection methods, one-way ANOVA was performed in terms of motion classification accuracy with SPSS Statistical Modeling Software (SPSS 17.0 IBM Corp., Chicago, IL). In addition, we also analyzed the effect of feature set, classifier, and combination of feature set and classifier on the classification performance by performing two-way ANOVA. The interaction effect of feature set and classifier on classification performance was firstly assessed. If the effect was insignificant, the main effects of features set and classifiers were examined. The level of statistical significance was set to p < 0.05.
Number of selected EMG channels using MCCSP
Optimal feature-classifier combination
Comparison of channel selection methods and electrode configuration
Comparison of distribution of selected EMG channels
Generally speaking, using more EMG electrodes could capture more electrophysiological information that may improve the performance of movement classification in EMG-PR-based control system such as multifunctional myoelectric prostheses and active rehabilitation robots. However, this would simultaneously increase the complexity and cost of the EMG controlled systems. Thus it becomes necessary and important to find an appropriate number of electrodes and their locations for the high performance of a myoelectric control system before it is clinically viable. Theoretically, to get the optimal channels from high density channels, all the combinations of desired number of channels should be involved in the pattern recognition analysis. But it often be impractical to search a global optimal solution when a large number of channels are involved. For example, searching 6 from 56 channels would result in 32,468,436 different combinations. So some suboptimal methods have been proposed and used in the previous studies [14, 16, 18, 19]. Two commonly used methods for channel selection are SFS and FMS, which are dependent on either classification accuracy or EMG features. While these methods might be applicable for channel selection, they have some limitations, as mentioned above in Introduction section. In this study, a novel direct channel selection algorithm was proposed and its performance in selecting the appropriate channels has been evaluated by a high-density EMG recordings from twelve mildly-impaired TBI patients. Since the newly proposed method (MCCSP) for channel selection is independent on EMG features and classification algorithms, it should be more convenient for channel selection in comparison to the SFS and FMS, especially when changing the EMG features or the classification algorithms is needed.
By examining two different types of EMG features and two pattern-recognition algorithms, we found that the TD features combining with a KNN classifier would be a better configuration for the selected EMG subset via MCCSP, and TDAR features with a LDA classifier would be a better one for the determined channels by using SFS (Figure 5). In addition, the electrode configuration is also a factor that would affect the motion classification performance. When the number of selected electrodes was increasing, the motion classification accuracy in bipolar electrode configuration reached a plateau (Figures 5 and 7) more quickly than that in monopolar electrode configuration. Furthermore, using same number of EMG channels, the bipolar electrode configuration always achieved higher classification accuracy than the monopolar electrode configuration (Figure 8). These findings are similar as the conclusions obtained by Zhou and Huang et al. [14, 16]. Generally, a monopolar electrode can provide more electrophysiological information than a bipolar electrode. However, different EMG applications may need different information involved in EMG signals. For the EMG-based motion classification, the useful information should be these components of EMG signals that are characteristic of discrepancy for different movements, which could make a classifier to easily identify a movement. This may be one reason why the bipolar electrodes would outperform the monopolar electrodes in movement classification. Another possible reason underlying these findings would be that compared to a monopolar electrode, a bipolar electrode acts as a differential operator that can remove some common components such as direct-current part and common-mode noise in EMG signals. These common components in different channels might not provide any discriminate information for movement classification, even making the classification performance worse. In future study, we are interested in exploring the reasons.
As a commonly used channel selection method in EMG-PR based classification, the previous studies have shown that SFS would be an effective method in selecting the appropriate EMG channels [7, 13, 14, 16]. However, in comparison to the proposed MCCSP, SFS was consistently worse just with one exception of using only one EMG channel (Figure 7). With SFS method, the first chosen EMG channel would be the global optimum. However, when more channels were considered, it was impossible to get the globally optimal EMG subset by using the SFS method. Compared with FMS, the proposed MCCSP performed slightly better (Figure 7) in classifying different arm/hand movements. In addition, the computational cost should be also an interesting performance metric in channel selection. In this study, we found that the average cost time taken by searching 18 optimal monopolar EMG channels was about 0.98 s when using the proposed MCCSP, significantly lower than that when using SFS (with TD feature and LDA classifier) (about 2096 s), FMS (with TD feature) (10.56 s), or FMS (with TDAR feature) (18.67 s), respectively. Additionally, a fixed combination of the selected EMG channels could be achieved when using MCCSP (Figure 9). These outcomes of this study may suggest that as an alternative, the MCCSP might be an effective and practicable choice for channel selection in the design of a practical myoelectric control system. Note that the channels selected by the three different methods in the study were inconsistent, but the similar classification performance was achieved by these selected channels. This would be because the three methods depend on different principles in channel selection. As mentioned above, as a direct channel selection algorithm the proposed MCCSP uses the raw EMG data to select the appropriate channels, the FMS and the SFS rely on EMG features and classification accuracy for searching the optimal channels, respectively. Thus the different channel combinations would be obtained by using the three methods. However, this does not mean there is no need to carefully select the channels. In the EMG-PR movement classification, it is the patterns of EMG signals from the selected channels that are considered as a whole to train and test a classifier for movement classification. Although different channel combinations were determined by the three methods, their EMG patterns all could provide identifiable information to properly classify different movements. So a similar classification performance could be achieved by the three methods. This may indicate that the channel combination for similar classification performance would be not only one.
It should be noted that all the results of this study were conducted by means of post-processing and evaluated with offline classification accuracy which might not reflect the real-time performance directly . In real-time applications, the EMG signals would be easily affected by some factors like muscle contraction force , skin impedance , electrodes shift [30, 32], arm position variation [33, 34], and muscle fatigue [30, 35]. The stability of the proposed channel selection algorithm was not considered in current study. Note that Huang et al. investigated the temporal stability of the SFS-selected electrodes by validating the EMG data recorded from follow-up experiment in terms of classification performance . The repeatability of EMG signals recorded from the selected electrode sites and their inter-parameter agreement  would be also necessary to examine the stability of channel selection algorithm. With these important issues, we will examine the stability of the channels selected by the proposed MCCSP in our future work.
This work was jointly supported by the National Key Basic Research Program of China (#2013CB329505), the National Natural Science Foundation of China under Grants (#61135004, #51275101), Guangdong Innovation Research Team Fund for Low-cost Healthcare Technologies, and External Cooperation Program of the Chinese Academy of Sciences (GJHZ1212). The authors would like to thank members at Center for Neural Engineering, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences for their assistance in experiment, data preprocessing and language editing. And also we would like to thank the cooperators at Guangdong Provincial Industrial Injury Rehabilitation Center for their assistance in the recruitment of subjects.
- Dorcas D, Scott R: A three-state myoelectic controller. Med Biol Eng 1966, 4: 367–372. 10.1007/BF02476154View ArticleGoogle Scholar
- Kawamoto H, Lee S, Kanbe S, Sankai Y: Power assist method for HAL-3 using EMG-based feedback controller. IEEE International Conference on Systems, Man and Cybernetics (SMC 03) OCT 05–08, 2003 Washington, D.C., USA 2003, 1648–1653.Google Scholar
- Dipietro L, Ferraro M, Palazzolo JJ, Krebs HI, Volpe BT, Hogan N: Customized interactive robotic treatment for stroke: EMG-triggered therapy. IEEE Trans Neural Syst Rehabil Eng 2005, 13: 325–334. 10.1109/TNSRE.2005.850423View ArticleGoogle Scholar
- Kuiken T, Li G, Lock B, Lipschutz R, Miller L, Stubblefield K, Englehart K: Targeted muscle reinnervation for real-time myoelectric control of multifunction artificial arms. JAMA-J Am Med Assoc 2009, 301: 619–628. 10.1001/jama.2009.116View ArticleGoogle Scholar
- Englehart K, Hudgins B, Parker PA, Stevenson M: Classification of the myoelectric signal using time-frequency based representations. Med Eng Phys 1999, 21: 431–438. 10.1016/S1350-4533(99)00066-1View ArticleGoogle Scholar
- Hudgins B, Parker P, Scott R: A new strategy for multifunction myoelectric control. IEEE Trans Biomed Eng 1993, 40: 82–94. 10.1109/10.204774View ArticleGoogle Scholar
- Li G, Schultz A, Kuiken T: Quantifying pattern recognition-based myoelectric control of multifunctional transradial prostheses. IEEE Trans Neural Syst Rehabil Eng 2010, 18: 185–192.View ArticleGoogle Scholar
- Hu XL, Tong K-y, Song R, Zheng XJ, Leung WWF: A comparison between electromyography-driven robot and passive motion device on wrist rehabilitation for chronic stroke. Neurorehab Neural Re 2009, 23: 837–846. 10.1177/1545968309338191View ArticleGoogle Scholar
- Wolbrecht ET, Chan V, Reinkensmeyer DJ, Bobrow JE: Optimizing compliant, model-based robotic assistance to promote neurorehabilitation. IEEE Trans Neural Syst Rehabil Eng 2008, 16: 286–297.View ArticleGoogle Scholar
- Liu J, Zhou P: A novel myoelectric pattern recognition strategy for hand function restoration after incomplete cervical spinal cord injury. IEEE Trans Neural Syst Rehabil Eng 2013, 21: 96–103.View ArticleGoogle Scholar
- Zhang X, Zhou P: High-density myoelectric pattern recognition toward improved stroke rehabilitation. IEEE Trans Biomed Eng 2012, 59: 1649–1657.View ArticleGoogle Scholar
- Rojas-Martinez M, Mananas MA, Alonso JF, Merletti R: Identification of isometric contractions based on High Density EMG maps. J Electromyogr Kinesiol 2013, 23: 33–42. 10.1016/j.jelekin.2012.06.009View ArticleGoogle Scholar
- Daley H, Englehart K, Hargrove L, Kuruganti U: High density electromyography data of normally limbed and transradial amputee subjects for multifunction prosthetic control. J Electromyogr Kines 2012,22(3):478–484. 10.1016/j.jelekin.2011.12.012View ArticleGoogle Scholar
- Zhou P, Lowery M, Englehart K, Huang H, Li G, Hargrove L, Dewald J, Kuiken T: Decoding a new neural-machine interface for control of artificial limbs. J Neurophysiol 2007, 98: 2974–2982. 10.1152/jn.00178.2007View ArticleGoogle Scholar
- Kendell C, Lemaire ED, Losier Y, Wilson A, Chan A, Hudgins B: A novel approach to surface electromyography: an exploratory study of electrode-pair selection based on signal characteristics. J Neuroeng Rehabil 2012, 9: 24. 10.1186/1743-0003-9-24View ArticleGoogle Scholar
- Huang H, Zhou P, Li GL, Kuiken TA: An analysis of EMG electrode configuration for targeted muscle reinnervation based neural machine interface. IEEE Trans Neural Syst Rehabil Eng 2008, 16: 37–45.View ArticleGoogle Scholar
- Farrell TR, Weir RF: A comparison of the effects of electrode implantation and targeting on pattern classification accuracy for prosthesis control. IEEE Trans Biomed Eng 2008, 55: 2198–2211.View ArticleGoogle Scholar
- Shibanoki T, Shima K, Tsuji T, Otsuka A, Chin T: A quasi-optimal channel selection method for bioelectric signal classification using a partial kullback-leibler information measure. IEEE Trans Biomed Eng 2013, 60: 853–861.View ArticleGoogle Scholar
- Kvas G, Velik R: A Filter Approach for Myoelectric Channel Selection. 2008 6th IEEE International Conference on Industrial Informatics; 13–16 July 2008; Daejeon, Korea 2008, 1437–1440.Google Scholar
- Qiang C, Hongbo Z, Jie C: The fisher-markov selector: fast selecting maximally separable feature subset for multiclass classification with applications to high-dimensional data. IEEE Trans Pattern Anal Mach Intell 2011, 33: 1217–1233.View ArticleGoogle Scholar
- Peng H, Fulmi L, Ding C: Feature selection based on mutual information criteria of max-dependency, max-relevance, and min-redundancy. IEEE Trans Pattern Anal Mach Intell 2005, 27: 1226–1238.View ArticleGoogle Scholar
- Nagata K, Ando K, Magatani K, Yamada M: Annual International Conference of the IEEE Engineering in Medicine and Biology Society 2007. Lyon, France: Cité Internationale; 2007:5214–5217.View ArticleGoogle Scholar
- Geng YJ, Zhang LQ, Tang D, Zhang XF, Li GL: Pattern Recognition Based Forearm Motion Classification for Patients with Chronic Hemiparesis. 35th Annual International Conference of the IEEE EMBS; 3 - 7 July, 2013; Osaka, Japan. 2013, 5918–5921.Google Scholar
- Wang Y, Gao S, Gao X: Common spatial pattern method for channel selection in motor imagery based brain-computer interface. 27th Annual International Conference of the IEEE Engineering in Medicine and Biology Society; Shanghai, China 2005, 5392–5395.Google Scholar
- Ramoser H, Muller-Gerking J, Pfurtscheller G: Optimal spatial filtering of single trial EEG during imagined hand movement. IEEE Trans Rehabil Eng 2000, 8: 441–446. 10.1109/86.895946View ArticleGoogle Scholar
- Hahne JM, Graimann B, Muller K-R: Spatial filtering for robust myoelectric control. IEEE Trans Biomed Eng 2012, 59: 1436–1443.View ArticleGoogle Scholar
- Graupe D, Cline WK: Functional separation of EMG signals via ARMA identification methods for prosthesis control purposes. IEEE Trans Syst Man Cybern 1975, SMC-5: 252–259.View ArticleGoogle Scholar
- Nazarpour K, Sharafat AR, Firoozabadi SMP: Application of higher order statistics to surface electromyogram signal classification. IEEE Trans Biomed Eng 2007, 54: 1762–1769.View ArticleGoogle Scholar
- Lock B, Englehart K, Hudgins B: Real-time myoelectric control in a virtual environment to relate usability vs. accuracy. Proceedings of the 2005 Myoelectric Controls/Powered Prosthetics Symposium; August 17–19, 2005. Fredericton, New Brunswick, Canada 2005.Google Scholar
- Tkach D, Huang H, Kuiken T: Study of stability of time-domain features for electromyographic pattern recognition. J Neuroeng Rehabil 2010, 7: 21. 10.1186/1743-0003-7-21View ArticleGoogle Scholar
- Doheny EP, Caulfield BM, Minogue CM, Lowery MM: Effect of subcutaneous fat thickness and surface electrode configuration during neuromuscular electrical stimulation. Med Eng Phys 2010, 32: 468–474. 10.1016/j.medengphy.2010.03.004View ArticleGoogle Scholar
- Young A, Hargrove L, Kuiken T: Improving myoelectric pattern recognition robustness to electrode shift by changing interelectrode distance and electrode configuration. IEEE Trans Biomed Eng 2012, 59: 645–652.View ArticleGoogle Scholar
- Geng Y, Zhou P, Li G: Toward attenuating the impact of arm positions on electromyography pattern-recognition based motion classification in transradial amputees. J Neuroeng Rehabil 2012, 9: 74. 10.1186/1743-0003-9-74View ArticleGoogle Scholar
- Geng Y, Zhang F, Yang L, Zhang Y, Li G: Reduction of the effect of arm position variation on real-time performance of motion classification. 34th Annual International Conference of the IEEE Engineering in Medicine and Biology Society; August 28 - September 1, 2012; San Diego, CA, United states 2012, 2772–2775.View ArticleGoogle Scholar
- Dideriksen JL, Farina D, Enoka RM: Influence of fatigue on the simulated relation between the amplitude of the surface electromyogram and muscle force. Philos Trans R Soc A - Math Phys Eng Sci 2010, 368: 2765–2781.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.