The Effect of Combining Stationary Wavelet Transform and Independent Component Analysis in the Multichannel SEMGs Hand Motion Identification System

Transcription

1 Journal of Medical and Biological Engineering, 6(): The Effect of Combining Stationary Wavelet Transform and Independent Component Analysis in the Multichannel SEMGs Hand Motion Identification System Yi-Chun Du, Liang-Yu Shyu Weichih Hu,* Department of Biomedical Engineering, Chung Yuan Christian University, Chung Li, Taiwan, 30 ROC Institute of Biomedical Engineering, Chung Kung University, Tainan, Taiwan,70 ROC Received Aug 005; Accepted 9 Jan 006 Abstract Independent component analysis (ICA) was studied to examine its effect on independent component extraction and data reduction for the hand motion identification system using multichannel surface electromyogram (SEMG) and stationary wavelet transform (SWT). The results indicated that features extracted from wavelet transformed SEMGs significantly increase the classification rate. However, these extra features extracted from different wavelet scales dramatically increase the amount of computation of the artificial neural network (A). When using the ICA to extract independent components, the number of training epochs that were required during the A training period increased. This study used the unsupervised FastICA to reduce the number of SEMG channels from 7 to 4 and achieved a data reduction rate of 43%. However the hand motion identification rates using the reduced channels are significantly lower (p < 0.05) when they are compared with the traditional method. The results demonstrate that ICA can effectively reduce the size of the neural network and in turn reduce the amount of required computation with the price of reduced identification rates. Keywords: Independent component analysis, Principal component analysis, Stationary wavelet transform, Surface Electromyogram Introduction It is known that activities of different muscle groups result in different surface electromyography (SEMG) patterns. Conversely, by analyzing these SEMG patterns, the corresponding muscle group action can be identified. Since 990, surface EMG has been used for prosthesis control [, ] and to construct a robotic manipulator control system [3] with high degree of success. In addition to prosthesis control, SEMG has also been used to identify the action of hands and legs during functional electrical stimulation (FES) [4]. These studies demonstrate that hand actions, such as extension of the forearm, hand opening and grasping or bending of the finger can be identified using SEMG signals from muscle groups that are not directly corresponding to that particular action. However, this research and most other studies required SEMG electrodes to be placed on designated locations for signals collection, which could not be easily performed outside the laboratory. Additionally, these studies identify only hand actions involving large muscle groups. o study has been conducted to identify finger actions using SEMG electrodes placed on the forearm. In 00, Shyu et al [5] proposed a hand action identification * Corresponding author: Weichih Hu Tel: ; Fax: weichih@ be.cycu.edu.tw system to discriminate hand actions including extension of the index finger, extension of the index and middle fingers, extension of the thumb and little fingers, extension of the thumb and index fingers, extension of all fingers except the thumb, and extension of all five fingers. The proposed system uses seven active electrodes placed around the forearm to collect SEMG signals from the major muscle groups of the forearm. Six features, including the integral of EMG (IEMG), waveform length (WL), variance (VAR), zero crossing (ZC), slope sign change (SSC) and Willison amplitude (WAMP), were extracted from each of the seven SEMG signals as the input to the artificial neural networks (A). The average discrimination rate was 93.% [5]. An average discrimination rate of 93.% appears to be an acceptable outcome. However, it was also true that for one of the subjects the discrimination rate was only 63.8% [6]. Additionally, for hand actions such as extension of the index and middle fingers and extension of all five fingers, the discrimination rates were low in most of the subjects. To make the system of any practical use, it is not only necessary to increase the precision of discrimination on average but also increase the precision for individual subject and individual hand action. Additional electrodes may be used to acquire more SEMGs to improve the hand action

2 0 J. Med. Biol. Eng., Vol. 6. o. 006 discrimination. However, this will not only require more hardware but also increase the complexity of the system. Alternately, there are signal processing methods, for example wavelet transform, that can reveal hidden features of the signal. Wavelet transformation can decompose an existing signal into multiple wavelet scales. Useful information can be extracted from different wavelet scales and be used successfully in other biomedical applications [7-0]. The potential drawback of wavelet transforming SEMGs is that it dramatically increases computation load. On the other hand, independent component analysis (ICA) has been used to separate the independent component from multi-channel EMGs [] and extract the independent component of motor unit action potentials of the recorded multichannel SEMG signals []. Thus, in this study, stationary wavelet transformation (SWT) is proposed to decompose the seven SEMG signals into multiple scales. The major advantage of SWT over discrete wavelet transformation is the preservation of time information of the original signal sequence at each level. Additionally, ICA is proposed to extract independent component from the transformed signals. Characteristic features can be extracted from these independent components and fed into the A identification system for hand action discrimination. The goal of this study is to establish a new and reliable method for hand action classification and, at the same time, without sacrificing the accuracy. Methodology A. stationary wavelet transform The wavelet transform of a function, f, with respect to a given mother wavelet, ψ, is defined as: x t ws f ( x) = f ψ s ( x) = f ( t) ψ dt () s s Where s is the scale factor. Let s = j (j Z, Z is the integral set), then the WT is called dyadic WT. The dyadic WT can be calculated as: j S j f ( n) = hk S j f ( n k) k Z j W j f ( n) = gks j f ( n k) k Z () (3) Where S j is a smoothing operator and W is the WT of j f () n digital signal f(n). On the other hand, hk s and g k s are the coefficients of corresponding low-pass and high-pass filters, respectively. The standard discrete wavelet transform decimates the wavelet coefficients at each level. Thus, the results of wavelet transform at each level are half the size of the original sequence. The stationary wavelet transform, on the other hand, pads the corresponding low-pass and high-pass filters with zeros and the two new sequences which each have the same length as the original sequence [3]. The major advantage of SWT is the preservation of time information of the original signal sequence at each level. On the other hand, the nd order Daubechies (db) wavelet transformation has been used successfully in other SEMG classifications [4]. Thus, in this study, the db stationary wavelet transform was used to decompose the seven SEMGs into six levels each. B. Independent component analysis ICA is a method for finding the linear transformation s = Wx to obtain s i that are as independent as possible in the sense of maximizing some independence measurement functions F(s,...,s m ), where x is the random vector. Using the approximation of negentropy as the index on statistical independence between different components and fixed-point algorithm for estimating the maximum likelihood, Hyvärinen developed the FastICA algorithm [5,6]. egentropy is defined as J(y) = H(y gauss )-H(y) (4) where y gauss is a Gaussian random variable of the same covariance matrix as y. This can be approximated as J(y) = [E{G(y)}-E{G(v)}] (5) where G is any non-quadratic function. More details about FastICA can be found at this reference [7]. A MatLab program for ICA computation was developed based on the FastICA code provided on the same website. C. EMG characteristic features Six EMG features, integral of EMG (IEMG), waveform length (WL), variance (VAR), zero crossing (ZC), slope sign change (SSC) and Willison amplitude (WAMP), were reported and used by several studies [8]. They are defined as: IEMG = k= WL = k= VAR = ZC = i= (6), where x = x x (7) k= [sgn( x x ) and x, sgn(x ) = 0, otherwise k k k k k+ k xk + 0.0] if x > 0 SSC = ( x ) ( x x ) for k k k+, k =,..., WAMP = f i= where f ( x) ( ) +, = 0,, where (8) (9) (0) () if x > 0. 3 otherwise They are either energy or frequency related. It was found that using two features, IEMG and WL, could achieve an averaged discrimination rate of 90% [9]. However, for some subjects and some particular hand motions, the classification

3 SEMG Discrimination Using SWT and ICA (a) Stationary wavelet transform (Extract features from 3 scales) (b) 7 channels SEMG ICA, select 4 most independent components Stationary wavelet transform (Extract features from 3 scales) neural network (c) Manually select 4 SEMG signals Stationary wavelet transform (Extract features from 3 scales) Figure. Block diagram of three experiment procedures. rates were lower than 40% [9]. Thus, modifications are proposed for these two features. Two new SEMG characteristic features that are modified from IEMG and WL are proposed to replace the six commonly used features. They are EMG envelope energy and EMG envelope and defined as: EMG _ Envelop _ Energy f ) () = k = Where f(x) = x, if x > threshold and f(x) = 0, otherwise. = EMG _ Envelop g ( f ( )) (3) T = k = Where f(x) =, if x >= threshold and f(x) = -, if x < threshold and Σf(x) > 0, f(x) = 0, otherwise. Whereas, g(x) = x, if x > threshold and g(x) = 0, otherwise. The EMG envelope energy is an improved version of the commonly used IEMG feature in that it only includes EMG signals where the amplitude is larger than the threshold in order to avoid the influence of noise. On the other hand, the EMG envelope is a combination of WL and slope sign change. All characteristic features were normalized before they were fed into the neural network. D. Signal Processing In order to correctly recognize hand actions, it is necessary to identify the beginning and end of action periods from continuously recorded SEMG signals. To avoid mis-identification of the action period due to noise and artifacts in the SEMG signals, we employed the multi-threshold method that was proposed in the previous study [5]. First, during the training period, while the subject was carrying out designated actions, thresholds and normalization factors for each SEMG channel were established based on the absolute value of the SEMG signal. Then, during the recognition period, SEMG signals from each channel were normalized and continuously compared with the corresponding thresholds. When two or more SEMG signals exceeded their corresponding thresholds, it was recognized as the beginning of one action. After the beginning of an action, when six or more SEMG signals were determined to be below the corresponding threshold, it was marked as the end of the action. ICA extracted four independent components from the original seven SEMG signals during each action period. Stationary wavelet decomposed these independent components into six wavelet signals. Characteristic features, EMG envelope and EMG envelope energy, were extracted and normalized prior to being fed into the A identification system. ( E. Experiment Design A total of hand motions were defined and studied in this study. Four of them were defined according to previous studies, including fist flexion up and down, pronation and supination [0]. The other seven hand motion were defined and used in a previous study including extension of the index finger, extension of index and middle fingers, extension of thumb and little fingers, extension of thumb and index fingers, extension of all fingers except the thumb, and extension of all five fingers [5]. Twelve subjects, nine males and three females, were recruited for this study, with an average age of 43 years. The new electrode system was placed on the right forearm five centimeters below the elbow. Experiments were conducted while subjects were standing. During the training period, each subject was asked to perform two hand actions, hand opening and closing, repeatedly. The SEMG signals acquired during this period were used to establish thresholds and normalization factors for each channel. After the initialization period, subjects were asked to carry out the hand actions 8 times in random order. The SEMGs of the first nine trials were used as the training set and the other nine were used as the testing set in the experiment. All SEMG signals were digitized using a 000Hz sampling rate and processed off-line using a MatLab program. To reduce the amount of the signal processing load, a pilot study was conducted to investigate the minimum and best combination of different wavelet scales that were required to maintain the discrimination rate for hand motion identification. Three different combinations, three scales (scales, 3 and 4), four scales (scales, 3, 4 and 5) or five scales (scales, 3, 4, 5 and 6), were tested. After the best combination of wavelet scales had been decided, three signal processing procedures were designed to test the effect of ICA. In the first procedure, no channel reduction was attempted, following SWT, features were extracted from seven SEMG signals as the input to the A (Fig. a). The second procedure used ICA to obtain four independent SEMG signals prior to the SWT and feature extraction (Fig. b). The third procedure manually selected four SEMG signals before the SWT and feature extraction (Fig. c). F. Performance Evaluation Due to the fact that the initial weights of A were selected randomly, the A may have different classification results following each successful training. Thus, the averaged

4 J. Med. Biol. Eng., Vol. 6. o. 006 Table. Results of traditional hand discrimination method. Where ()-() represent hand action flexion down, extension of thumb and little fingers, extension of index and middle fingers, extension of all five fingers, extension of thumb and index fingers, extension of the index finger, extension of all fingers except the thumb, pronation, supination and flexion up, respectively. Hand actions Discrimination Rate % 9.7.3% % % % % % % % 9.9.3% % % scales Table. Results of wavelet scale selection. Discrimination Rate Training Epochs Convergence Probability % # %* % # % % # %* * p < 0.05; # p < 0.05 vs. traditional method (a) (b) (c) (d) (e) (f) (g) (h) Figure. A typical result of SWT. Where (a) is the original SEMG signal, (b) (h) are signals of different wavelet scale. otice that the slow frequency components in (b) are not visible in the original signal. discrimination rate after five different trainings will be used for analysis in this study. To evaluate the performance of different procedures, the tested method not only needs to compare the accuracy of discrimination but also assess the feasibility of new procedures. In addition to the amount of computation required by A, the converging speed and ease of A convergence are good measures of system feasibility. The size of the A is directly proportional to the amount of computation required. In contrast, the number of training epoch is commonly used to evaluate the speed of convergence. To estimate the ease of A convergence, in this study, the convergence probability was defined as: convergence probability = 5 / (number of training to achieve five A convergences) * 00%. For example, for a training data set, if the A converged in five successive trainings, the convergence probability is 5 / 5 *00% = 00%. As opposed to, if the A achieve five convergences after ten trainings, the convergence probability will be 5 / 0 * 00% = 50%. Results First, for comparison, the collected SEMG signals were processed by the traditional hand identification procedure [] that uses six SEMG features from seven SEMG signals. On average, the discrimination rate of the traditional method is 90.9% (Table ) which is slightly lower than previously reported []. A typical result of SWT is illustrated in Fig.. It can be seen that, at different levels, the SWT reviews some hidden feature from the original SEMG signal. For example, the slow frequency components in scale six are not visible in the original signal. Thus, with features obtained from multiple wavelet scales, significant improvement in discrimination rate is expected. However, the size of A will increase due to the increased information from multiple wavelet scales, which in turn increase the computational load. The purpose of wavelet scale selections pilot study is to discover the best combination of wavelet scales that can both minimize the computational load and achieve an adequate discrimination rate. The results of wavelet scale selections study are illustrated in Table II. It is shown that using either three scales (scales, 3 and 4), four scales (scales, 3, 4 and 5) or five scales (scales, 3, 4, 5 and 6) provides significantly higher discrimination rates as compared to the traditional method. Among the three possible wavelet selections, there is no difference in the average discrimination rates. The number of training epochs increases as the number of selected scales decreases. In contrast, the convergence probability increases as the number of selected scales increases. There is a small but significant different in convergence probability between the three scales and five scales. The combination of scales, 3, and 4 uses the least amount of features in the process of classification and achieves comparable discrimination rates as compared to other combinations, thus, this combination was selected and used in the following study. After the determination of wavelet scales, the effect of channel reduction was studied. According to the experiment design, in this study, two channel reduction methods were investigated. First, two possible collections of four SEMG signals were manually selected from the seven input SEMGs. The second method selected four independent component signals according to the result of unsupervised FastICA. After extracting the independent components, the four signals that were considered the most independent were selected and used in the following study (Fig. 3).

5 SEMG Discrimination Using SWT and ICA 3 Table 3. Results of channel reduction. The discrimination rates after channel reduction were significantly lower than the procedure that incorporates no channel reduction. The 4.8% data reduction rates were achieved by reducing seven SEMG signals into four signals. Procedures Discrimination Rate Training Epochs Convergence Probability Data Reduction Without channel reduction % * * 83.7% 0 % Channel reduction using ICA %*,# *,# % 4.8 % Channel reduction by selecting channels, 3, 5, 7 Channel reduction by selecting channels, 4, 6, 7 * p < 0.05, # p < % * * 75.9.%* 4.8 % % # % # 7.73%* 4.8 % (a) (b) Figure 3. A typical result of ICA. The original seven SEMG signals are illustrated in (a). (b) are the results of ICA independent component extraction. The four signals, denoted as ICS-ICA4, are the four most independent components. Table 3 summarizes the system performances of three signal processing procedures including: ICA channel reduction, manual channel reduction and without channel reduction. It can be seen that, with channel reduction, performances significantly decreased as compared to those without channel reduction in terms of discrimination rate and the number of training epoch. However, both channel reduction methods, ICA and manual selection, reduced the seven SEMG signals to four and achieved a 4.8% data reduction. When compared with manual channel reduction, the discrimination rate and the number of training epoch following ICA reduction were significantly better. The convergence probability of ICA was better too, although the difference was not significant. This demonstrates that the ICA can reduce the amount of computation but cannot preserve the necessary information for discrimination. Discussion The significantly superior discrimination rates, when using SWT and two new SEMG features, indicate that the characteristic features extracted from different scales definitely provide valuable information in hand motion discrimination at the cost of larger A. At the same time, it demonstrates that the two new SEMG features can replace the six commonly used features, after SWT. The results also indicate that most of the useful SEMG features are in the mid-frequency range. On the other hand, the number of wavelet scales that were included in the process does not significantly affect the discrimination rates. An ideal A easily converges during the training phase and achieves a high classification rate during trial. However, not all A perform ideally. It is found that, for some input features combination, it is very difficult to train the A. However, once it is successfully trained, the A performs very well. On the other hand, for some other input features combination, the A may converge very easily during the training period but perform poorly during trial. For practical

6 4 J. Med. Biol. Eng., Vol. 6. o. 006 usage, an A that is relatively easy to train and performs well is much more desirable than an A that is difficult to train and has an excellent performance. In this study, the idea of convergence ratio is proposed. It measures how easily the A can be trained. The results indicate that by increasing the number of wavelet scales the convergence ratio increases, i.e. the A is more likely to converge. In spite of the significant different in convergence ratios between procedures using the combination of three wavelet scales and procedure using five wavelet scales, this study chose to use the procedure that uses three wavelet scales. This is due to the consideration that this difference is small and by using three wavelet scales requires much less computation, in the same time, without significantly reducing the discrimination rate. Moreover, it is found that when less then three wavelet scales were included for features extraction the A failed to converge. ICA, when it was used for the purposes of channel reduction, failed to maintain the level of accuracy that the system had prior to the channel reduction. However, when compared with manual channel reduction, the discrimination rate and the number of training epoch after ICA reduction were significantly better (Table 3). The convergence probability of ICA was better too, although the difference was not significant. The channel reduction through manually selecting four channels from the existing seven SEMG signals proved to be unreliable. In fact, the A failed to converge following manual channel selection for two of the twelve subjects. On the other hand, the inferior performance of ICA may be due to the fact that a fixed number of independent component signals were selected in this study. A further study of a reliable method to select appropriate independent signals after ICA will be conducted in the future. Conclusion This study proposes the use of SWT to extract features from multiple wavelet scales in order to increase the performance of the hand motion identification system using custom-made electrode. Additionally, two new SEMG features were also proposed. Extra features from different wavelet scales increase the performance of the system as is illustrated in this result. The results also indicate that not all wavelet scales are needed in the process of hand motion identification. Selecting only three wavelet scales achieves the same discrimination rate as using five wavelet scales. On the other hand, the results of this study indicate that the ICA is a reliable pre-process to reduce the amount of data and, in turn, reduce the amount of computation in hand motion classification. However, the performance of the system decreased significantly following ICA. References [] K. Kuribayashi, K. Okimura and T. Taniguchi; A Discrimination System Using eural etwork for EMG-controlled Prostheses, Proc. IEEE/RSJ Internal. Conference on Intelligent Robots and Systems, 63-88, 993 [] H.P. Huang and C.Y. Chiang, Development of a Myoelectric Discrimination System for a Multi-Degree Prosthetic Hand, Proc. IEEE Intenal. Conference on Robotics & Automation, 3: , 999 [3] O. Fukuda, T. Tsuji, H. Shigeyoshi and M. Kaneko, An EMG Controlled Human Supporting Robot Using eural etwork, Proc. IEEE/RSJ Internal Conference on Intelligent Robots and Systems, 3: , 999 [4] K.D. Strange, Restoration of use of paralyzed limb muscles using sensory nerve signals for state control of FES-assisted walking, IEEE Trans. Rehab. Eng., 7(3): , 999 [5] L.Y. Shyu, J.Y. Chen, R.W. Tatn, and W.C. Hu, A new electrode system for hand action discrimination, J. Med. d Biol. Eng., (4): -7, 00. [6] J.Y. Chen Multi-electrode hand motion identification system, Master thesis, Chung Yuan Christian University, Chung Li, Taiwan, 00 [7] B. Celler and P. de Chazal, Low computational cost classifiers for ECG diagnosis using neural networks, Proc. the 0 th Annual Internal Conference of IEEE Eng. Med. Biol. Society, 0(): , 998. [8] E. M. Lalitha, and L. Satish, Wavelet analysis for classification of multi-source PD patterns, IEEE Trans. Dielectrics and Electrical Insulation, 7(): 40-47, 000. [9] P.J. Sparto, M. Parnianpour, E.A. Barria, and J.M. Jagadeesh, Wavelet and short-time Fourier transform analysis of electromyography for detection of back muscle fatigue, IEEE Trans. Rehab. Eng., 8(3): , 000. [0] L.Y. Shyu, Y.S. Lin, C.P. Liu and W.C. Hu, The detection of impedance cardiogram characteristic points using wavelet transform, Computers in Biology and Medicine, 34: 65 75, 004. [] M.J. McKeown, R. Radtke, Phasic and tonic coupling between EEG and EMG demonstrated with independent component analysis, J. Clinical europhysioogy, 8(): 45 57, 00. [] H. akamura, M. Yoshida, M. Kotani, K. Akazawa, T. Moritani, The application of independent component analysis to the multichannel surface electromyographic signals for detecion of motor unit action potential trains: Part I measuring techniques, J. Electromyography Kinesiol., 4: 43-43, 004. [3] G.P. ason and B.W. Silverman, The stationary wavelet transform and some statistical application, Lecture otes in Statistics, 03: 8-99, 995. [4] P. Welling, and G. Moschytz, Analysis of wavelet features for myoelectric signal classification, Proc. IEEE Internal Conference Electronics, Circuits and Systems, 3: 09-, 998. [5] A. Hyvärinen, Survey on independent component analysis, eural Computing Surveys, : 94-8, 999. [6] A. Hyvärinen, E. Oja, A fast fixed-point algorithm for independent component analysis, eural Comput. 9: , 997. [7] A. Hyvärinen, The FastICA MATLAB package, Available: http: // /projects /ica /icademo/. [8] H.P. Huang and C.Y. Chiang, Development of a Myoelectric Discrimination System for a Multi-Degree Prosthetic Hand, Proc. IEEE Internal. Conf. on Robotics & Automation, 3: , 999. [9] Y.C. Du, Evaluation of stationary wavelet, PCA and ICA for hand motion EMG identification, M.S. thesis, Dep. Biomed. Eng., Chung-Yuan Christian Univ. Chung Li, Taipei, Taiwan, 003. [0] H.P. Huang and C.Y. Chiang, DSP-based controller for a multi-degree prosthetic hand, Proc. IEEE Internal. Conf. Robotics & Automation, : , 000.