High-density surface EMG maps from upper-arm and forearm muscles
© Rojas-Martínez et al.; licensee BioMed Central Ltd. 2012
Received: 17 December 2011
Accepted: 29 November 2012
Published: 10 December 2012
sEMG signal has been widely used in different applications in kinesiology and rehabilitation as well as in the control of human-machine interfaces. In general, the signals are recorded with bipolar electrodes located in different muscles. However, such configuration may disregard some aspects of the spatial distribution of the potentials like location of innervation zones and the manifestation of inhomogineties in the control of the muscular fibers. On the other hand, the spatial distribution of motor unit action potentials has recently been assessed with activation maps obtained from High Density EMG signals (HD-EMG), these lasts recorded with arrays of closely spaced electrodes. The main objective of this work is to analyze patterns in the activation maps, associating them with four movement directions at the elbow joint and with different strengths of those tasks. Although the activation pattern can be assessed with bipolar electrodes, HD-EMG maps could enable the extraction of features that depend on the spatial distribution of the potentials and on the load-sharing between muscles, in order to have a better differentiation between tasks and effort levels.
An experimental protocol consisting of isometric contractions at three levels of effort during flexion, extension, supination and pronation at the elbow joint was designed and HD-EMG signals were recorded with 2D electrode arrays on different upper-limb muscles. Techniques for the identification and interpolation of artifacts are explained, as well as a method for the segmentation of the activation areas. In addition, variables related to the intensity and spatial distribution of the maps were obtained, as well as variables associated to signal power of traditional single bipolar recordings. Finally, statistical tests were applied in order to assess differences between information extracted from single bipolar signals or from HD-EMG maps and to analyze differences due to type of task and effort level.
Significant differences were observed between EMG signal power obtained from single bipolar configuration and HD-EMG and better results regarding the identification of tasks and effort levels were obtained with the latter. Additionally, average maps for a population of 12 subjects were obtained and differences in the co-activation pattern of muscles were found not only from variables related to the intensity of the maps but also to their spatial distribution.
Intensity and spatial distribution of HD-EMG maps could be useful in applications where the identification of movement intention and its strength is needed, for example in robotic-aided therapies or for devices like powered- prostheses or orthoses. Finally, additional data transformations or other features are necessary in order to improve the performance of tasks identification.
KeywordsHigh-Density surface electromyography 2D electrode arrays EMG pattern recognition Artifact detection Rehabilitation Robotics Prosthetics
The electromyographic (sEMG) signals detected on the skin above human muscles are important to infer motion intention and therefore could be used to control devices such as exoskeletons, biofeedback systems or assistive tools [1, 2]. Central to these goals is the extraction of information from myoelectric signal which is commonly detected with electrode pairs. Other purposes include the estimation of muscle force, biofeedback of the activity of the muscle and the analysis of myoelectric fatigue.
The main disadvantage of bipolar signals is that its amplitude depends on the distance between the active motor units and the recording electrodes . Due to the low spatial resolution of the bipolar signal, the standard surface EMG reflects the activity of a number of motor units (MU) within a delimited area of the muscle. However, amplitude variations are expected in both, the parallel and the perpendicular directions of propagation of the MU action potentials (MUAP): in the former, the amplitude of the signal varies with the proximity to innervation zones and tendons, while in the latter the amplitude is attenuated because of the propagation properties of the conductor volume . As pointed out by Zwarts et al. in , single channel approaches do not reflect the physical propagation of the potentials and therefore only the time-varying properties of the signals are usually analyzed, disregarding important spatial aspects of the propagation like extent and length of the muscle fibers, which are essential for the force-generating capacity of the muscle, and, if not well addressed can lead to incorrect conclusions. In recent years, on the other hand, the development and application of electrode arrays in one or two dimensions have allowed the study of the sEMG signal in the temporal and the spatial domain, opening new possibilities to the study of the neuromuscular system [5–7] and to the field of myoelectric control .
What is more, recent studies have demonstrated that the muscles do not activate homogeneously, that is, distinct regions of the muscle are activated differentially depending on the position of the joint  and the duration  and strength of the contraction . Such activation may be related to bundles of fiber types organized in different regions within the muscle, each of them following different recruitment strategies according to Henneman’s size principle .
Therefore, EMG amplitude information provided by a single bipolar channel is highly dependent on the location of the recording electrodes, even when they are well located away from innervation zones and tendons and it does not offer the possibility of tracking inhomogeneities in the activation of the muscles.
The recording of sEMG signals with 2D arrays in a wide area of the muscles and the processing of the signal in the space dimension [5, 12] can overcome some of the drawbacks of single-channel approaches, providing a quantification of the temporal and spatial properties of the electrical muscle activity [13, 14]. In this study, High Density surface EMG (HD-EMG) signals recorded with 2-dimensional (2D) arrays of closely spaced electrodes were used to calculate activation maps for the upper arm and forearm. These maps provided a larger amount of information related to the tracking of (task changing) skin surface areas where EMG amplitude is maximal and a better estimation of muscle activity by the proper selection of the most significant channels. The main objective is to analyze patterns in the activation maps associating them with four movement directions at the elbow joint and with different strengths of those movements. Although the activation pattern can be assessed with bipolar electrodes, HD-EMG signals from upper-arm and forearm muscles could enable the extraction of features that depend on the spatial distribution of the MUAPs and on the load-sharing between muscles in order to have a better differentiation between tasks and effort levels.
For this purpose, two additional automatic algorithms are proposed, one intended for the automatic detection of low-quality signals and the other, for the segmentation of global activation areas. In the first case, different algorithms based on the spatial spread of the potentials over the skin surface have been proposed in the literature [9, 15, 16]. Such algorithms considered only the information of the bulk of data for the detection of outliers in 2D multichannel recordings. The approach proposed here takes into account features extracted from channels in the close- proximity neighborhood to evaluate the quality of a given signal and to interpolate its value from neighboring channels when needed. Besides, the method was designed to compromise between precision and sensitivity of the detection in order to reduce the number of misclassifications of good-quality signals.
In the second case, watershed segmentation has been previously proposed for determining the regions of activity in HD-EMG recordings , however this method is prone to over-segmentation in the presence of multiple discontinuities , particularly, those introduced by dissimilarities in the electrode-skin impedance of the different channels. The segmentation proposed here offers an alternative, being less sensitive to local maxima. Results show that features extracted from HD-EMG maps could be useful in the identification of movement intention.
Three 2-D electrode arrays were fabricated in hydrophobic fabric for the recording of high density surface EMG signals (HD-EMG). The arrays were developed in our laboratory and consisted in silver-plated and gel-filled eyelets (external diameter of 5 mm) equally spaced by 10 mm in rows (y in the proximal-distal direction) and columns (x in the medial-lateral direction). Textile fabric adapts to the geometry of the muscle while maintaining inter-electrode distance and good adhesion to the skin provided by elastic straps. Array 1 (forearm) had 6 rows and a variable number of columns (between 17 and 19) depending on the dimensions of the limb of the subject whereas arrays 2 and 3 (upper arm) consisted of 120 electrodes distributed in 8 rows × 15 columns. About 350 channels were recorded for each subject and were sufficient to cover the muscles of interest.
Twelve healthy male volunteers (age, 28.3 ± 5.5 years; height: 177.8 ± 6.0 cm; body mass: 75.7± 8.7 kg) participated in the experiment. Subjects did not have any history of neuromuscular disorders, pain or regular training of the upper limb. All subjects gave informed consent to participate to the experimental procedure.
Five muscles associated with flexion, extension, pronation and supination of the forearm were selected for the study: Biceps and Triceps Brachii in the upper arm and Anconeus, Brachioradialis and Pronator Teres in the forearm. Array 1 (A1) was located on the forearm, with the most proximal electrode row 2 cm below the elbow crease, and was intended for covering the three forearm muscles whose edges were previously drawn over the surface of the skin according to Kendall et al. . The columns of the array laid along the axis oriented from lateral to medial direction in order to cover each selected muscle with at least three columns. Array 2 and array 3 (A2 and A3) were located in the distal and proximal regions of the upper arm, respectively, with their centers on the location recommended by SENIAM . Previous to the positioning of the electrode array, the skin was shaved and cleaned with abrasive paste. Conductive gel (20 μl) was inserted in each electrode of the array using a gel dispenser (Multipette Plus, Eppendorf, Germany). Lengths and circumferences of the upper arm and forearm segments were measured for each subject. The length of the ventral side of the upper-arm was measured from the acromion to the fossa cubit whereas the length of its dorsal side was measured from the posterior crista of the acromion to the olecranon. The length of the forearm was measured from the medial epicondyle to the apophysis of the radius. Circumferences of the arm segments were measured while contracting different muscles: the distal and proximal upper arm circumference were measured over the muscle belly of biceps and triceps respectively and the proximal forearm circumference was measured over the muscle belly of the Brachio Radialis (approximately 2 cm below the elbow crease). During the experiment, subjects sat in front of the mechanical brace with the back straight, the elbow joint flexed at 45°, shoulder abducted at 90° (arm in the sagittal plane), and forearm rotated 90°, midway between supination and pronation (Figure 1). Subject’s Maximal Voluntary Contraction (MVC) during flexion, extension, pronation and supination were obtained at the beginning of the experiment as the maximum of three trials for each task. Afterwards, subjects carried out a series of isometric contractions at 10%, 30% and 50% MVC. Contractions were performed in randomized order, each lasting for 10 seconds followed by 2 minutes rest in order to avoid cumulative fatigue. In addition, subjects were previously trained to maintain the hand and fingers at rest during signal recording.
Detection of low quality signals
When recording a large number of physiological signals it is very likely to observe some low quality channels affected by artifacts originated in deformations of the skin under the electrodes, movement of the cables or in bad contacts between each recording electrode and the skin, inducing capacitive couplings and enabling power line interference . Visual inspection of the outliers (channels with low-quality signals) is time-consuming and depends on the expertise of the operators. Thus, it becomes necessary to apply an automatic method to identify such signals (and perform adequate processing if necessary) before any kind of information extraction.
Relative power of low frequency components P l/t , from 0 to 12 Hz
Relative power of power-line components P line/t corresponding to 50Hz and its first four harmonics, and
Signal power calculated from the root mean square (RMS) of the signal
The Power Spectral Density of the signal was estimated with the FFT in non-overlapping signal epochs of 500 ms. Features (P l/t , P line/t and RMS) were computed for each channel as the mean of the values obtained from six epochs over segments of 3s where the exerted force remained constant.
Automatic algorithm for artifact removal
An expert system based on thresholds associated with the three features described previously was designed.
where IQR represents the Interquartile Range. The median was chosen instead of the mean because of its lower sensitivity to outliers since it considers the highest breakpoint (50%), that is, the smallest percentage of outliers that can cause an estimator to take arbitrary large values .
Sensitivity analysis for constant K line for the detection of power-line artifacts
k line (2.5)
Training and validation
where TP and TN is the number of channels correctly identified as low and good-quality signals respectively, FN is the number of low-quality signals not identified by the algorithm and FP is the number of good-quality signals identified as low-quality.
Signals were visually classified as low-quality or not by three experts based on the observation of similarity between different channels and on the examination of baseline fluctuations or periodicity patterns (related to movement artifacts and power-line interference respectively, see Figure 2a as an example). Two databases for training and validation were obtained, each composed of 20 signal sets S selected from different contractions, effort levels and arrays 1 to 3. Fleiss’ Kappa index  was used to measure agreement between experts, scoring 82.63% and 86.19% for the training and validation databases respectively and indicating an “almost perfect agreement”. The information of the three experts was combined by obtaining the majority vote in each case (i.e. the statistical mode of the three opinions) in order to obtain a binary classification-label as artifact/non artifacts for each single-channel in the set S.
Optimal values for k 1 and k 2 were tuned according to the following criteria: 1) by Receiver Operating Characteristic (ROC) curves (S vs. 1-SP), widely used in signal detection theory and clinical diagnostics  and, 2) by Precision-Recall representations (P vs. SP), which are commonly used in machine learning . Both methods assessed the accuracy of the prediction (or outcome) of the described method based on its Sensitivity, Specificity and Precision. The optimal classifier was found as a tradeoff between hit rates and false alarm rates. In the first case the optimal can be found as the minimum distance between the curve (S vs. 1-SP) and the point [0, 1], and in the second case, as the minimum distance between curve (P vs. SP) and the point [1, 1].
Segmentation of HD-EMG maps
where s is an EMG signal in the channel located at the position i, j of the 2D array (as explained before), N is the total number of samples in an epoch of 500 ms and the RMS value was averaged in M=6 non overlapped epochs corresponding to three seconds of signal. Prior to the calculation of RMS, the signals were filtered between 12 and 350 Hz with a 4th order Butterworth filter in forward and backward direction in order to correct for phase distortion following SENIAM recommendations for the processing of surface EMG signals . RMS values corresponding to signals previously identified as artifacts were replaced by triangle-based cubic interpolation based on Delaunay Triangulation for surface fitting proposed in .
Each channel in the map can be considered as a pixel located at the positions i and j of the electrode array whose intensity is given by I i,j . In other words, the map I can be thought of three dimensions where the intensity values represent elevations as in a topographical map. Intensity levels in the maps correspond to the activation level of a muscle (or a set of muscles) during the development of a specific task.
where b, the structuring element, is a disc of radius 1 and ⊖ and ⊕ are the operators for dilation and erosion respectively .
Average HD-EMG maps
Average maps for the 12 subjects who participated in the study were obtained by averaging individual segmented maps, in order to obtain useful information related to muscle co-activation pattern in terms of intensity values and its spatial distribution during different kind of tasks and levels of effort. Considering that upper-limb dimensions, specifically circumference and length, are different for every subject, it was necessary to normalize and interpolate the image in the 2-D space so that results could be comparable among subjects allowing the calculation of an average map for the population. In the case of arrays A2 and A3 (biceps and triceps), the zero of the coordinate system was defined to correspond to landmarks defined by SEMIAM project . In the case of forearm muscles, the origin of the x-axis laid in the intersection of the line connecting the origin and insertion of each muscle (Anconeus, Brachioradialis or Pronator Teres) and an arch traced around the forearm, 2 cm below the elbow crease which in turn was the zero of the y-axis. Values in the x dimension were normalized with respect to the circumferences of the different arm segments: proximal forearm for array A1, distal upper-arm for array A2 and proximal upper-arm for array A3, and values in the y dimension were normalized with respect to the muscle length as described in the protocol description. In order to obtain individual maps with the same (i, j) coordinates for the calculation of population’s average map, RMS values in the 2-D space were re-sampled by cubic splines interpolation in the x and y directions, considering units relative to the total length and total circumference of the arm segment as explained in the previous paragraph.
where Q k dim is the value of the projection of the maximum of the segmented map D’ h at the k th coordinate along the dimension dim= x or dim=y and μ dim corresponds to the median of the projection Q k dim . The median of the projections permitted to evaluate spatial shifts along the x and y- axes of the maps, both of them associated with changes in effort levels or even with different tasks.
On the other hand, data similar to that obtained with bipolar electrodes was extracted by selecting two monopolar channels for each muscle in the 2D arrays: the first electrode corresponded to the one located at the origin of the coordinate system (i.e. reference in SENIAM recommendations ) as described previously for each muscle, and the second one was located 10 mm apart in the direction of the muscular fibers. For each muscle, a single differential signal was obtained from these two channels and its corresponding RMSav-bip value was calculated at the same time-interval as in the case of HD-EMG maps. The variables RMSav-bip and RMSav-HD (from bipolar or HD-EMG configurations) were used in the identification of tasks and their performances were compared as later explained.
The statistical analysis was intended to assess differences between information extracted from single bipolar signals or from HD-EMG maps and also to analyze differences due to type of task and effort level. Such analysis was based in the variables determined in the previous section, that is, RMSav-bip and RMSav-HD both of them related to the signal power, and μx and μy associated with the spatial distribution of the maps.
Factors considered in the statistical analysis were: 1) the type of electrode (bipolar or HD-EMG), 2) The type of task, that is, flexion, extension, supination and pronation, and contraction level (10%, 30% and 50% MVC), and 3) Muscle (i.e. biceps, triceps, brachioradialis, anconeus and pronator teres).
Differences between RMSav-bip and RMSav-HD were evaluated with the non-parametric Friedman test in two-way layout. In order to avoid differences between the type of recording (bipolar and monopolar), both variables in each muscle were normalized with respect to the mean value obtained in all tasks and contraction levels. Data corresponding to different tasks and effort levels were pooled together by considering blocks with replication of cells of the factor muscle according to the procedure described in  and implemented in the statistical toolbox of MATLAB®.
On the other hand, the capacity of the extracted variables (RMSav-bip or RMSav-HD) for the identification of tasks and/or effort levels was evaluated by classifying the data into different groups based on linear discriminant analysis  and cross-validation with the Leaving One Out Method . Data was classified into four or twelve groups corresponding to type of task or to type of task and effort level respectively. In the former case, samples corresponding to the three levels of effort were pooled together for each type of task. The overall classification performance was obtained in terms of Accuracy, Sensitivity, Specificity and Precision as described in  and (Eq. 5). For this analysis TP were data samples well classified into any of the classes, FN corresponded to the number of missing samples in any of the classes, that is, samples belonging to a given class but that were classified in another, TN were not misclassified samples, and FP were samples misclassified in any of the classes.
In addition, changes in the spatial distribution of the maps due to type of task and effort level were analyzed with a non-parametric repeated measures design based on the Friedman test. In this case, variations of the variables μ x and μ y were evaluated in 12 different measures corresponding to four types of task by three levels of effort each. A Bonferroni correction was applied in order to take into account the multiple comparisons.
Finally, pair-wise comparisons were analyzed trough non-parametric Wilcoxon signed rank test.
Statistical significance was set to p=0.05. In the Friedman test, χ2 statistics was considered significant for χ2 (d.o.f=1)> 3.84 for one degree of freedom (d.o.f = 2 types of electrode −1) and for χ2 (d.o.f=11)> 27.28 for eleven degrees of freedom (d.o.f=12 measures −1) after the Bonferroni correction.
Detection of low-quality signals
Sensitivity analysis in the training database for constant k line in Eq. 3 is presented in Table 1. Higher values of k line increased the number of TP and FP affecting both, the sensitivity and the precision of the detection. As it can be observed on Table 1, k line between 2 and 3 is a good compromise between sensitivity and precision and it can be confirmed from the accuracy and specificity, where the highest values were obtained. For values lower than 2, precision decreases whereas for values higher than 3, sensitivity decreases.
Optimum values of k1 and k2 (Eq.2 and Eq.4) and their performance indexes
(k 1, k 2 )
Although high Sensitivity (S) and Specificity (SP) were desired, a lower number of FP became important when considering the next step of the analysis where RMS values of artifacts were interpolated from neighbor channels. If too many channels in the neighborhood were wrongly labeled as “artifacts” (i.e. too many FP), then the interpolation process was not possible. Because of the lower Precision of the algorithm with the ROC criterion (Table 2), constants k 1 and k 2 were finally selected according to optimization results obtained by Precision- Recall.
On the other hand, the algorithm for artifact detection had low computational complexity. The execution time of the algorithm, (mean and standard deviation for joint training and validation databases), was 201.2±12.85 ms, [min: 187.5 ms, max 234.4 ms] per signal set on a 2.13 GHz Intel® Core2™ processor. Each set had a total duration of 3s and was composed by a different number of channels between 109 and 120 channels.
The segmentation produced intensity maps for each muscle and subject, and allowed the calculation of average HD-EMG maps representative for the 12 subjects in the study.
Bipolar vs. High-density EMG signals
Significant differences were observed for the two-way Friedman test between RMSav-bip and RMSav-HD using the factor muscle as blocking factor (χ2 = 16.55, p<0.001). Thus, significant differences were observed when characterizing the different tasks and effort levels with information extracted from one or the other type of sensor.
Average HD-EMG Maps
Variability between the 12 individual maps associated with the average HD-EMG maps
Biceps (mean, [min max])
0.085, [0.07, 0.11]
0.18, [0.16, 0.2]
0.29, [0.23, 0.39]
Triceps (mean, [min max])
0.13, [0.084, 0.23]
0.15, [0.089, 0.2]
0.24, [0.11, 0.36]
Brachio Radialis (mean, [min max])
0.14, [0.063, 0.25]
0.16, [0.13, 0.21]
0.23, [0.17, 0.26]
Anconeus (mean, [min max])
0.13, [0.11, 0.14]
0.18, [0.12, 0.22]
0.2, [0.15, 0.23]
Pronator Teres (mean, [min max])
0.12, [0.063, 0.2]
0.16, [0.11, 0.2]
0.22, [0.15, 0.29]
In addition to the average maps at 50% MVC, their projections Q k x at 10%, 30% and 50% MVC for the two most active muscles in each type of task are presented in Figure 9. Each projection was normalized with respect to the maximum value reached at 50% MVC in order to observe changes of the intensity as a function of effort level and to compare among muscles and/or tasks. Changes observed are summarized as following:
Flexion: Intensity decreased with decreasing levels of effort (from 50 to 10% MVC) in Biceps and Brachioradialis (p<0.0005), preserving a similar proportion in both muscles (Figure 9a).
Extension: Intensity in Triceps decayed proportionally from 50% to 10% MVC (p<0.0005), but not in the Anconeus, where the maximum of the distribution was significantly different between 10% and 30% MVC (p<0.009) but not between 30% and 50% MVC (Figure 9b).
Supination: Intensity decreased from 50 to 10% MVC in Biceps and Anconeus (p<0.0005) but not in the same proportion in both muscles. Changes in the activation of the Anconeus were proportionally higher at 10% and 30% MVC when compared to changes in biceps, (p<0.003), although the former continued being the most active muscle in the contraction (Figure 9c).
Pronation: Intensity levels in Pronator Teres and Anconeus decreased similarly from 50 to 10% MVC. However, the normalized levels for Anconeus at 10%MVC were higher than those for Pronator Teres (p<0.001) showing that the contribution of the former was higher at lower levels of contraction. (Figure 9d).
Results for the repeated measured friedman test for the variables μ x and μ y extracted from HD_EMG
Furthermore, pair-wise comparisons between spatial-distribution variables obtained for different effort levels of the same task on a given muscle in Figure 9 were assessed by applying a Wilcoxon signed rank test to the medians. For example, μx in the Biceps, shifted to the left with increasing levels of contraction during flexion (p<0.02) and during supination (p<0.05). Such shifts were likely due to differences in the activation of the two heads of this muscle. Additionally, significant shifts were also found for the Pronator Teres during pronation (p<0.05).
Discussion and conclusions
The main objective of the study was to extract information from HD-EMG maps that could be associated with four tasks at the elbow joint (forearm pronation and supination, and elbow flexion and extension) at different effort levels. This objective has been reached using 2D arrays of electrodes on five muscles in the upper-arm and forearm. It was shown no only that the average signal power extracted from HD-EMG maps may improve the differentiation of tasks and effort levels but also that the spatial distribution of the maps differed between tasks. Variables related with the spatial distribution of the intensity may complement the information provided by the amplitude of the signals in the identification of motion intention. Additionally, average HD-EMG maps for the four types of tasks and for the group of 12 voluntary subjects were obtained.
Detection of low-quality signals
Due to the high number of sEMG channels, an algorithm for the automatic detection of low quality signals was developed. In this study, the artifactual channels were removed and interpolated based on RMS values of neighbor channels in the HD-EMG maps.
As the number of channels affected by artifacts is usually much lower than the number of non- artifact channels, the number of TN is usually high, so the accuracy and especially, the specificity, is in general very high (>99%). In this study, the percentage of low quality signals was between 0 and 13% of the total number of channels of each set in the training database. For this reason, the ROC method, which was intended to compromise between specificity and sensitivity, provided an overestimated value of the former (the latter was always very high as it was explained above) at the expense of reducing the precision of the algorithm. On the other hand, PR method maximized the sensitivity while preserving the precision. Considering that both, the specificity and the accuracy were always very high, the use of PR was more convenient in our case since we were dealing with offline detection and interpolation of artifacts. In addition, the accuracy was slightly higher with the PR method. Other applications, especially those intended for online detection, should consider the ROC approach in order to increase the sensitivity as much as possible.
In addition, the sensitivity analysis for the constants k line , k 1 and k 2 in Tables 1 and 2 showed that the selected values represented a good compromise between the Precision and Sensitivity of the detection.
Different approaches have also been suggested for the detection of low quality signals. In , a non-supervised method based on local distance-based outlier factor was proposed for HD-EMG signals recorded from the same muscles and electrodes arrays. That method did not required a training process and successfully detected low-quality signals with an average Accuracy of 91.9%, Sensitivity of 96.9% and Specificity of 96.4%. If the latter index is considered (refer to SP in Eq. 5), it is possible that such method was prone to the inclusion of FP, given that the number of TN is always very high as explained before. This condition was not desired in our study because intensity values associated with low quality signals were later replaced in the maps based on neighbor channels, and if those were wrongly identified as artifacts (FP), the interpolation for the replacement was not possible. Thus, another method intended to minimize the number of FP while preserving the accuracy of the detection was proposed in the present work. Moreover, the method in  did not took into account information provided by neighbor channels as in the algorithm described here.
A different study by Gronlund et al. in  assumed unimodal distribution of the amplitude of the signals and therefore it was not applicable to cases where the multichannel recording involved various regions with different levels of activity (as in Figure 2) or even with no activity at all. The latter corresponded to regions localized far away from the main activation areas or to muscles marginally active during the contraction (an example can be observed in Figure 6 bottom-left for the signals recorded in array A1-forearm during flexion). In contrast, the algorithm proposed in the present work labeled a signal as an artifact based on the amplitude information of 6 neighboring channels (Eq. 4), avoiding general assumptions on the distribution of the potentials recorded in the array as proposed in .
Finally, the algorithm reached very high values for the performance indexes in the validation dataset (see Table 2). In addition, the channels identified as artifacts were correctly interpolated from the information of neighboring non-artifact channels even at the edges of the map (see Figure 5), which provided a smooth surface for subsequent stages of the analysis. Therefore, both the proposed methodology and features showed a very good performance for the detection and off-line replacement of low-quality signals.
Segmentation of HD-EMG maps
An automatic algorithm for the segmentation of active zones was proposed. This segmentation allowed the calculation of average HD-EMG maps for the population of 12 subjects by determining the ranges of active zones in the x and y axes relative to upper-limb circumference and muscle length respectively, and referred to the electrode location recommended for sEMG recording on the analyzed muscles.
Furthermore, it is important to note that maps segmentation was not the aim of this study neither the determination of anatomical muscle regions but was an intermediate step in order to obtain average maps. It is known that potentials’ amplitude diffuses across skin surface, so the actual size of active muscle regions might be overestimated by the segmentation proposed. In spite of this, the segmentation permitted the extraction of surface areas corresponding to different muscles in individual maps and also permitted to focus on the regions of major activation, avoiding other active neighboring areas that could be more affected by the activity of nearby muscles (see Figure 6). By averaging individual segmented maps, it was possible to obtain representative HD-EMG maps associated with the activation of each muscle during the different tasks and effort levels. In this sense, the segmentation worked properly.
Average HD-EMG maps
Variability in the levels of intensity between subjects with respect to the average HD-EMG maps (Table 3) was found to be low enough to consider such maps as representative for the 12 subjects. The obtained differences were higher for higher levels of contraction (between 20% and 29% in the four tasks at 50% MVC) and being less than 12% for contractions at 10% MVC (see Table 3). In the average activation maps depicted in Figure 8, it is possible to observe changes in the activation pattern corresponding to different tasks. Such differences were confirmed with the classification by LDA shown in Figure 7. Consequently, results concerning the spatial distribution and the levels of intensity from average maps can be considered as globally associated with the muscle function and to the activation pattern of the subjects in the study.
On the other hand, the classification performance was higher for the average power extracted from HD-EMG maps (RMSav-HD) than from bipolar electrodes (RMSav-bip) either for four or twelve classes, especially when considering the precision and sensitivity of the classification. Therefore it is possible to conclude that information extracted from the amplitude of signals recorded in high-dimensional configuration has more power to differentiate between tasks and even effort levels than single bipolar signals. With this respect, a recent study by Tkach et al.  based on information extracted from bipolar signals showed that the classification accuracy for the identification of motor tasks worsens when considering distinct strengths of the same motion. Substantial drops were observed when training and testing the classifier with data of mixed high and low effort levels, obtaining a maximal accuracy of ~80% for the classification into the four tasks described in the present work. In our case, an overall classification accuracy of approximately 90% for 4 motor tasks with mixed data from very-low to medium-high effort levels (10%, 30% and 50% MVC) was obtained. What is more, classification accuracy for 12 classes was also in the order of 90%. However, in both cases the precision and sensitivity were not very high, thus additional data transformations or other features are necessary in order to improve the identification performance. In this sense, the Friedman test showed that variables related to the spatial distribution of the maps (μ x and μ y , see Table 4) may also assist in the discrimination between types of tasks and effort levels.
In addition, Tkach et al. also found that classification accuracy dropped when bipolar electrodes shifted by15 mm . Such shifts could be due to relative movements between the recording electrodes and the skin or because of errors when positioning the sensors, for example in different days. Thus another advantage of HD-EMG maps relies on the contact redundancy implied by the recording of a number of signals over a large surface of the muscle, as well as in the possibility of extracting features associated with spatial-changes induced by the central nervous system in the control of the muscles . All of this makes HD-EMG maps more robust to errors introduced by contact artifacts and by the relative location of the electrodes with respect to the origin of the potentials, especially in contractions involving joint movement or sensor repositioning.
When analyzing variables associated with the activation maps for the 12 subjects, it was possible to observe differences in the co-activation pattern of the muscles according to the kind of task and the effort level. For example during flexion at 50% MVC (see Figure 6), the Biceps, and the Brachioradialis were the most active muscles (as expected) but there was also an important activation of the Anconeus and of the Pronator Teres likely to stabilize the elbow joint and to compensate for the supination action of the biceps. The extension at 50% MVC was mainly produced by the Triceps and the Anconeus but the other two analyzed muscles in the forearm were also active. All selected muscles but Triceps were involved during Supination at 50% MVC with similar intensity levels among them. Finally, during Pronation at 50% MVC, naturally the most active muscle was the Pronator Teres but both the Anconeus and the Brachioradialis were also active while the muscles of the upper-arm appeared not to be active during the contraction. Additionally, when considering activation patterns for contractions at 30% and 10% MVC it was found that the levels of intensity did not decreased proportionally with effort level in all of the muscles, showing changes in the load-sharing of the involved muscles.
Differences in the average maps between 10%, 30% and 50% MVC were not only related to the levels of intensity but also to its spatial distribution. It was possible to observe differences in their projection over the x-axis (see Figure 9 and 10 and Table 4) for different effort levels depending on the task. These variations corresponded to shifts in the lateral to medial axis when the level of effort changed from 10% to 50% MVC. Furthermore, similar values of intensity were obtained for different kind of tasks in some muscles and differences were only found in the spatial distributions of the maps.
For all of this, HD-EMG maps instead of single bipolar signals, and variables related to maps intensity and spatial distribution might be useful in applications where identification of movement intention is needed: for example in robotic-aided therapies focused on the improvement of muscle coordination where interaction between patient and machine is involved and where the robot has to be able to sense patient’s intention . Additionally, other applications implying proportional control for devices like powered- prostheses or orthoses could benefit from information provided by HD-EMG maps regarding not only the task but also its strength.
Diven Right Leg
High Density Surface Electromyography
Maximal Voluntary Contraction
Power Spectral Density
Precision- Recall representation
Root Mean Square
Receiver Operating Characteristics
The experimental part of this work was carried out in the Laboratory for Engineering of the Neuromuscular System (LISiN), Politecnico di Torino, Italy, under the supervision of Prof R. Merletti. This work was partially supported by the MINECO of the Spanish Government (project DPI2011-22680), by the Grant BE-00194 from AGAUR and by the CIDEM (Generalitat de Catalunya) with the project Neurorehab 3E+D (RD08-2-0019). We are grateful to Dr. Roberto Merletti and Dr. Hamid Reza Marateb for reviewing a draft of this paper.
- Nishikawa D, Wenwei Y, Yokoi H, Kakazu Y: EMG prosthetic hand controller using real-time learning method. Proc IEEE Int Conf Syst Man Cybern 1999, 1: 153-158.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(3):325-334. 10.1109/TNSRE.2005.850423PubMed CentralView ArticlePubMedGoogle Scholar
- Drost G, Stegeman DF, Schillings ML, Horemans HLD, Janssen HMHA, Massa M, et al.: Motor unit characteristics in healthy subjects and those with postpoliomyelitis syndrome: A high-density surface EMG study. Muscle Nerve 2004,30(3):269-276. 10.1002/mus.20104View ArticlePubMedGoogle Scholar
- Kleine BU, van Dijk JP, Lapatki BG, Zwarts MJ, Stegeman DF: Using two-dimensional spatial information in decomposition of surface EMG signals. J Electromyogr Kinesiol 2007,17(5):535-548. 10.1016/j.jelekin.2006.05.003View ArticlePubMedGoogle Scholar
- Zwarts MJ, Stegeman DF: Multichannel surface EMG: Basic aspects and clinical utility. Muscle Nerve 2003,28(1):1-17. 10.1002/mus.10358View ArticlePubMedGoogle Scholar
- Holobar A, Minetto MA, Botter A, Negro F, Farina D: Experimental analysis of accuracy in the identification of motor unit spike trains from high-density surface EMG. IEEE Trans Neural Syst Rehabil Eng 2010,18(3):221-229.View ArticlePubMedGoogle Scholar
- Merletti R, Holobar A, Farina D: Analysis of motor units with high-density surface electromyography. J Electromyogr Kinesiol 2008,18(6):879-890. 10.1016/j.jelekin.2008.09.002View ArticlePubMedGoogle Scholar
- Farina D, Lorrain T, Negro F, Jiang N: High-density EMG E-textile systems for the control of active prostheses. Proc IEEE Int Conf Eng Med Biol Soc 2010, 3591-3593.Google Scholar
- Vieira TMM, Merletti R, Mesin L: Automatic segmentation of surface EMG images: Improving the estimation of neuromuscular activity. J Biomech 2010,43(11):2149-2158. 10.1016/j.jbiomech.2010.03.049View ArticlePubMedGoogle Scholar
- Tucker K, Falla D, Graven-Nielsen T, Farina D: Electromyographic mapping of the erector spinae muscle with varying load and during sustained contraction. J Electromyogr Kinesiol 2009,19(3):373-379. 10.1016/j.jelekin.2007.10.003View ArticlePubMedGoogle Scholar
- Holtermann A, Roeleveld K, Karlsson JS: Inhomogeneities in muscle activation reveal motor unit recruitment. J Electromyogr Kinesiol 2005,15(2):131-137. 10.1016/j.jelekin.2004.09.003View ArticlePubMedGoogle Scholar
- Staudenmann D, Kingma I, Daffertshofer A, Stegeman DF, van Dieën JH: Heterogeneity of muscle activation in relation to force direction: A multi-channel surface electromyography study on the triceps surae muscle. J Electromyogr Kinesiol 2009,19(5):882-895. 10.1016/j.jelekin.2008.04.013View ArticlePubMedGoogle Scholar
- Merletti R, Aventaggiato M, Botter A, Holobar A, Marateb H, Vieira TMM: Advances in surface EMG: Recent progress in detection and processing techniques. Crit Rev Biomed Eng 2010,38(4):305-345. 10.1615/CritRevBiomedEng.v38.i4.10View ArticlePubMedGoogle Scholar
- Merletti R, Botter A, Cescon C, Minetto MA, Vieira TMM: Advances in surface EMG: Recent progress in clinical research applications. Crit Rev Biomed Eng 2010,38(4):347-379. 10.1615/CritRevBiomedEng.v38.i4.20View ArticlePubMedGoogle Scholar
- Gronlund C, Roeleveld K, Holtermann A, Karlsson JS: On-line signal quality estimation of multichannel surface electromyograms. Med Biol Eng Comput 2005,43(3):357-364. 10.1007/BF02345813View ArticlePubMedGoogle Scholar
- Marateb HR, Rojas-Martinez M, Mansourian M, Merletti R, Mananas Villanueva MA: Outlier detection in high-density surface electromyographic signals. Med Biol Eng Comput 2012,50(1):79-89. 10.1007/s11517-011-0790-7View ArticlePubMedGoogle Scholar
- Vincent L: Morphological grayscale reconstruction in image analysis: Applications and efficient algorithms. IEEE Trans Image Process 1993,2(2):176-201. 10.1109/83.217222View ArticlePubMedGoogle Scholar
- Metting Van Rijn A, Peper A, Grimbergen C: High-quality recording of bioelectric events. Med Biol Eng Comput 1990,28(5):389-397. 10.1007/BF02441961View ArticlePubMedGoogle Scholar
- Kendall F, McCreary E, Provance P: Muscles: Testing and Function. 4th edition. Baltimore, Md: Williams & Wilkins; 1993.Google Scholar
- Freriks B, Hermens HJ: SENIAM9: European Recommendations for Surface ElectroMyoGraphy, Results of the SENIAM Project, ISBN: 90-75452-14-4 (CD-rom). The Netherlands: Roessingh Research and Development; 1999.Google Scholar
- Clancy EA, Morin EL, Merletti R: Sampling, noise-reduction and amplitude estimation issues in surface electromyography. J Electromyogr Kinesiol 2002,12(1):1-16. 10.1016/S1050-6411(01)00033-5View ArticlePubMedGoogle Scholar
- Hampel FR: A general qualitative definition of robustness. Ann Math Stat 1971,42(6):1887-1896. 10.1214/aoms/1177693054View ArticleGoogle Scholar
- Farina D, Colombo R, Merletti R, Baare Olsen H: Evaluation of intra-muscular EMG signal decomposition algorithms. J Electromyogr Kinesiol 2001,11(3):175-187. 10.1016/S1050-6411(00)00051-1View ArticlePubMedGoogle Scholar
- Fleiss JL: Measuring nominal scale agreement among many raters. Psychol Bull 1971,76(5):378-382.View ArticleGoogle Scholar
- Fawcett T: An introduction to ROC analysis. Pattern Recog Lett. 2006,27(8):861-874. 10.1016/j.patrec.2005.10.010View ArticleGoogle Scholar
- Landgrebe TC, Paclik P, Duin RP, Bradley AP: Precision-recall operating characteristic (P-ROC) curves in imprecise environments. In Proceedings of 18th International Conference on Pattern Recognition. Los Alamitos, California USA . IEEE Computer Society 2006, 4: 123-127.Google Scholar
- Hermens H, Freriks B, Disselhorst-Klug C, Rau G: Development of recommendations for sensors and sensor placement procedures. J Electromyogr Kinesiol 2000,10(5):361-374. 10.1016/S1050-6411(00)00027-4View ArticlePubMedGoogle Scholar
- Bradford Barber C, Dobkin DP, Huhdanpaa H: The quickhull algorithm for convex hulls. ACM Trans Math Softw 1996,22(4):469-483. 10.1145/235815.235821View ArticleGoogle Scholar
- Serra J: Image analysis and mathematical morphology. London: Academic; 1982.Google Scholar
- Zar JH: Biostatistical analysis. 5th edition. New Jersey: Prentice Hall; 2010.Google Scholar
- Krzanowski WJ: Principles of multivariate analysis: A user's perspective. New York: Oxford University Press; 1988.Google Scholar
- Kearns M, Ron D: Algorithmic stability and sanity-check bounds for leave-one-out cross-validation. Neural Comput 1999,11(6):1427-1453. 10.1162/089976699300016304View ArticlePubMedGoogle Scholar
- Tkach D, Huang H, Kuiken TA: Study of stability of time-domain features for electromyographic pattern recognition. J Neuroeng Rehabil 2010, 7: 21. 10.1186/1743-0003-7-21PubMed CentralView ArticlePubMedGoogle Scholar
- Hogan N, Krebs HI, Rohrer B, Palazzolo JJ, Dipietro L, Fasoli SE, et al.: Motions or muscles? some behavioral factors underlying robotic assistance of motor recovery. J Rehabil Res Dev 2006,43(5):605-618. 10.1682/JRRD.2005.06.0103View ArticlePubMedGoogle 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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.