The Influence of the sEMG Amplitude Estimation Technique on the EMG–Force Relationship

Validation of the filling factor index to study the filling process of the sEMG signal in the quadriceps

INTRODUCTION

The EMG filling factor is an index to quantify the degree to which an EMG signal has been filled. Here, we tested the validity of such index to analyse the EMG filling process as contraction force was slowly increased.

METHODS

Surface EMG signals were recorded from the quadriceps muscles of healthy subjects as force was gradually increased from 0 to 40% MVC. The sEMG filling process was analyzed by measuring the EMG filling factor (calculated from the non-central moments of the rectified sEMG).

RESULTS

(1) As force was gradually increased, one or two prominent abrupt jumps in sEMG amplitude appeared between 0 and 10% of MVC force in all the vastus lateralis and medialis. (2) The jumps in amplitude were originated when a few large-amplitude MUPs, clearly standing out from previous activity, appeared in the sEMG signal. (3) Every time an abrupt jump in sEMG amplitude occurred, a new stage of sEMG filling was initiated. (4) The sEMG was almost completely filled at 2-12% MVC. (5) The filling factor decreased significantly upon the occurrence of an sEMG amplitude jump, and increased as additional MUPs were added to the sEMG signal. (6) The filling factor curve was highly repeatable across repetitions.

CONCLUSIONS

It has been validated that the filling factor is a useful, reliable tool to analyse the sEMG filling process. As force was gradually increased in the vastus muscles, the sEMG filling process occurred in one or two stages due to the presence of abrupt jumps in sEMG amplitude.

Tutorial. Surface electromyogram (sEMG) amplitude estimation: Best practices

This tutorial intends to provide insight, instructions and "best practices" for those who are novices-including clinicians, engineers and non-engineers-in extracting electromyogram (EMG) amplitude from the bipolar surface EMG (sEMG) signal of voluntary contractions. A brief discussion of sEMG amplitude extraction from high density sEMG (HDsEMG) arrays and feature extraction from electrically elicited contractions is also provided. This tutorial attempts to present its main concepts in a straightforward manner that is accessible to novices in the field not possessing a wide range of technical background (if any) in this area. Surface EMG amplitude, also referred to as the sEMG envelope [often implemented as root mean square (RMS) sEMG or average rectified value (ARV) sEMG], quantifies the voltage variation of the sEMG signal and is grossly related to the overall neural excitation of the muscle and to peripheral parameters. The tutorial briefly reviews the physiological origin of the voluntary sEMG signal and sEMG recording, including electrode configurations, sEMG signal transduction, electronic conditioning and conversion by an analog-to-digital converter. These topics have been covered in greater detail in prior tutorials in this series. In depth descriptions of state-of-the-art methods for computing sEMG amplitude are then provided, including guidance on signal pre-conditioning, absolute value vs. square-law detection, selection of appropriate sEMG amplitude smoothing filters and attenuation of measurement noise. The tutorial provides a detailed list of best practices for sEMG amplitude estimation.

A real-time and convex model for the estimation of muscle force from surface electromyographic signals in the upper and lower limbs

Surface electromyography (sEMG) is a signal consisting of different motor unit action potential trains and records from the surface of the muscles. One of the applications of sEMG is the estimation of muscle force. We proposed a new real-time convex and interpretable model for solving the sEMG-force estimation. We validated it on the upper limb during isometric voluntary flexions-extensions at 30%, 50%, and 70% Maximum Voluntary Contraction in five subjects, and lower limbs during standing tasks in thirty-three volunteers, without a history of neuromuscular disorders. Moreover, the performance of the proposed method was statistically compared with that of the state-of-the-art (13 methods, including linear-in-the-parameter models, Artificial Neural Networks and Supported Vector Machines, and non-linear models). The envelope of the sEMG signals was estimated, and the representative envelope of each muscle was used in our analysis. The convex form of an exponential EMG-force model was derived, and each muscle's coefficient was estimated using the Least Square method. The goodness-of-fit indices, the residual signal analysis (bias and Bland-Altman plot), and the running time analysis were provided. For the entire model, 30% of the data was used for estimation, while the remaining 20% and 50% were used for validation and testing, respectively. The average R-square (%) of the proposed method was 96.77 ± 1.67 [94.38, 98.06] for the test sets of the upper limb and 91.08 ± 6.84 [62.22, 96.62] for the lower-limb dataset (MEAN ± SD [min, max]). The proposed method was not significantly different from the recorded force signal (p-value = 0.610); that was not the case for the other tested models. The proposed method significantly outperformed the other methods (adj. p-value < 0.05). The average running time of each 250 ms signal of the training and testing of the proposed method was 25.7 ± 4.0 [22.3, 40.8] and 11.0 ± 2.9 [4.7, 17.8] in microseconds for the entire dataset. The proposed convex model is thus a promising method for estimating the force from the joints of the upper and lower limbs, with applications in load sharing, robotics, rehabilitation, and prosthesis control for the upper and lower limbs.

Use of Surface Electromyography to Estimate End-Point Force in Redundant Systems: Comparison between Linear Approaches

Estimation of the force exerted by muscles from their electromyographic (EMG) activity may be useful to control robotic devices. Approximating end-point forces as a linear combination of the activities of multiple muscles acting on a limb may lead to an inaccurate estimation because of the dependency between the EMG signals, i.e., multi-collinearity. This study compared the EMG-to-force mapping estimation performed with standard multiple linear regression and with three other algorithms designed to reduce different sources of the detrimental effects of multi-collinearity: Ridge Regression, which performs an L2 regularization through a penalty term; linear regression with constraints from foreknown anatomical boundaries, derived from a musculoskeletal model; linear regression of a reduced number of muscular degrees of freedom through the identification of muscle synergies. Two datasets, both collected during the exertion of submaximal isometric forces along multiple directions with the upper limb, were exploited. One included data collected across five sessions and the other during the simultaneous exertion of force and generation of different levels of co-contraction. The accuracy and consistency of the EMG-to-force mappings were assessed to determine the strengths and drawbacks of each algorithm. When applied to multiple sessions, Ridge Regression achieved higher accuracy (R2 = 0.70) but estimations based on muscle synergies were more consistent (differences between the pulling vectors of mappings extracted from different sessions: 67%). In contrast, the implementation of anatomical constraints was the best solution, both in terms of consistency (R2 = 0.64) and accuracy (74%), in the case of different co-contraction conditions. These results may be used for the selection of the mapping between EMG and force to be implemented in myoelectrically controlled robotic devices.

A Simulation Study to Assess the Factors of Influence on Mean and Median Frequency of sEMG Signals during Muscle Fatigue

Mean and Median frequency are typically used for detecting and monitoring muscle fatigue. These parameters are extracted from power spectral density whose estimate can be obtained by several techniques, each one characterized by advantages and disadvantages. Previous works studied how the implementation settings can influence the performance of these techniques; nevertheless, the estimation results have never been fully evaluated when the power density spectrum is in a low-frequency zone, as happens to the surface electromyography (sEMG) spectrum during muscle fatigue. The latter is therefore the objective of this study that has compared the Welch and the autoregressive parametric approaches on synthetic sEMG signals simulating severe muscle fatigue. Moreover, the sensitivity of both the approaches to the observation duration and to the level of noise has been analyzed. Results showed that the mean frequency greatly depends on the noise level, and that for Signal to Noise Ratio (SNR) less than 10dB the errors make the estimate unacceptable. On the other hand, the error in calculating the median frequency is always in the range 2-10 Hz, so this parameter should be preferred in the tracking of muscle fatigue. Results show that the autoregressive model always outperforms the Welch technique, and that the 3rd order continuously produced accurate and precise estimates; consequently, the latter should be used when analyzing severe fatiguing contraction.