Extracting Individual Muscle Drive and Activity From High-Density Surface Electromyography Signals Based on the Center of Gravity of Motor Unit

A motor unit action potential-based method for surface electromyography decomposition

OBJECTIVE

Surface electromyography (EMG) decomposition is crucial for identifying motor neuron activities by analyzing muscle-generated electrical signals. This study aims to develop and validate a novel motor unit action potential (MUAP)-based method for surface EMG decomposition, addressing the limitations of traditional blind source separation (BSS)-based techniques in computation complexity and motor unit (MU) tracking.

METHODS

Within the framework of the convolution kernel compensation algorithm, we developed a MUAP-based decomposition algorithm by reconstructing the MU filters from MUAPs and evaluated its performance using both simulated and experimental datasets. A systematic analysis was conducted on various factors affecting decomposition performance, including MU filter reconstruction methods, EMG covariance matrices, MUAP extraction techniques, and extending factors. The proposed method was subsequently compared to representative BSS-based techniques, such as convolution kernel compensation.

MAIN RESULTS

The MUAP-based method significantly outperformed traditional BSS-based techniques in identifying more MUs and achieving better accuracy, particularly under noisy conditions. It demonstrated superior performance with increased signal complexity and effectively tracked motor units consistently across decompositions. In addition, directly applying the MU filters reconstructed from MUAPs to decomposition exhibited marked computational efficiency.

CONCLUSION AND SIGNIFICANCE

The MUAP-based method enhances EMG decomposition accuracy, robustness, and efficiency, offering reliable motor unit tracking and real-time processing capabilities. These advancements highlight its potential for clinical diagnostics and neurorehabilitation, representing a promising step forward in non-invasive motor neuron analysis.

Wrist torque estimation by combining motor unit discharges with musculoskeletal model

OBJECTIVE

The application of electromyography (EMG) decomposition techniques in myoelectric control has gradually increased. However, most decomposition-based control schemes rely on machine learning, lacking interpretation of the biological mechanisms underlying movement generation and requiring large datasets for training. As neuromusculoskeletal modeling provides a promising alternative, this study proposes a decomposition-based musculoskeletal model for simultaneous and proportional myoelectric control.

METHODS

Sixteen able-bodied subjects participated in two experiments involving isometric wrist contractions in two degrees of freedom (DoF). High-density surface EMG signals and torques were recorded simultaneously. The EMG signals were decomposed into motor unit action potential trains (MUAPts). We proposed four clustering methods (two activation-based and two action potential-based) to group MUAPts, from which three neural features were extracted as neural excitations and input to the musculoskeletal model.

MAIN RESULTS

An activation-based clustering method with the twitch feature achieved a relatively high accuracy (R2 = 0.791 ± 0.101 and 0.622 ± 0.148 in the two experiments) with the highest smoothness (Roughness = 1.389 ± 0.211 and 1.140 ± 0.159).

CONCLUSION AND SIGNIFICANCE

The proposed MUAPt-based musculoskeletal model achieved promising accuracy in estimating continuous 2-DoF wrist torques, providing a novel approach for understanding the neuromechanical properties of multi-DoF movements and advancing the development of dexterous rehabilitation and robotic control.

Neuro-Musculoskeletal Modeling for Online Estimation of Continuous Wrist Movements from Motor Unit Activities

Decoding movement intentions from motor unit (MU) activities remains an ongoing challenge, which restricts our comprehension of the intricate transition mechanism from microscopic neural drive to macroscopic movements. This study presents an innovative neuro-musculoskeletal (NMS) model driven by MU activities for online estimation of continuous wrist movements. The proposed model employs a physiological and comprehensive utilization of MU firings and waveforms, thus facilitating the localization of MUs to muscle-tendon units (MTU) as well as the computation of MU-specific neural excitation. Subsequently, the MU-specific neural excitation was integrated to form the MTU-specific neural excitation, which were then inputted into a musculoskeletal model to accomplish the joint angle estimation. To assess the effectiveness of this model, high-density surface electromyography and angular data were collected from the forearms of eight subjects during their performance of wrist flexion-extension task. Two pieces of 8×8 electrode arrays and a motion capture system were employed for data acquisition. Following offline model calibration with a global optimization algorithm, online angle estimation results demonstrated a significant superiority of the proposed model over the state-of-the-art NMS models (p < 0.05), yielding the lowest normalized root mean square error ( 0.10 ± 0.02 ) and the highest determination coefficient ( 0.87 ± 0.06 ). This study provides a novel idea for the decoding of joint movements from MU activities. The research findings hold the potential to advance the development of NMS models towards the control of multiple degrees of freedom, with promising applications in the fields of motor control, biomechanics, and neuro-rehabilitation engineering.

Integration of Motor Unit Filters for Enhanced Surface Electromyogram Decomposition During Varying Force Isometric Contraction

Muscles generate varying levels of force by recruiting different numbers of motor units (MUs), and as the force increases, the number of recruited MUs gradually rises. However, current decoding methods encounter difficulties in maintaining a stable and consistent growth trend in MU numbers with increasing force. In some instances, an unexpected reduction in the number of MUs can even be observed as force intensifies. To address this issue, in this study, we propose an enhanced decoding method that adaptively reutilizes MU filters. Specifically, in addition to the normal decoding process, we introduced an additional procedure where MU filters are reused to initialize the algorithm. The MU filters are iterated and adapted to the new signals, aiming to decode motor units that were actually activated but cannot be identified due to heavy superimposition. We tested our method on both simulated and experimental surface electromyogram (sEMG) signals. We simulated isometric signals (10%-70%) with known MU firing patterns using experimentally recorded MU action potentials from forearm muscles and compared the decomposition results to two baseline approaches: convolution kernel compensation (CKC) and fast independent component analysis (fastICA). Our method increased the decoded MU number by a rate of 135.4% ± 62.5 % and 63.6% ± 20.2 % for CKC and fastICA, respectively, across different signal-to-noise ratios. The sensitivity and precision for MUs decomposed using the enhanced method remained at the same accuracy level (p <0.001) as those of normally decoded MUs. For the experimental signals, eight healthy subjects performed hand movements at five different force levels (10%-90%), during which sEMG signals were recorded and decomposed. The results indicate that the enhanced process increased the number of decoded MUs by 21.8% ± 10.9 % across all subjects. We discussed the possibility of fully capturing all activated motor units by appropriately reusing previously decoded MU filters and improving the balance of activated motor unit numbers across varying excitation levels.

Inferring position of motor units from high-density surface EMG

The spatial distribution of muscle fibre activity is of interest in guiding therapy and assessing recovery of motor function following injuries of the peripheral or central nervous system. This paper presents a new method for stable estimation of motor unit territory centres from high-density surface electromyography (HDsEMG). This completely automatic process applies principal component compression and a rotatable Gaussian surface fit to motor unit action potential (MUAP) distributions to map the spatial distribution of motor unit activity. Each estimated position corresponds to the signal centre of the motor unit territory. Two subjects were used to test the method on forearm muscles, using two different approaches. With the first dataset, motor units were identified by decomposition of intramuscular EMG and the centre position of each motor unit territory was estimated from synchronized HDsEMG data. These positions were compared to the positions of the intramuscular fine wire electrodes with depth measured using ultrasound. With the second dataset, decomposition and motor unit localization was done directly on HDsEMG data, during specific muscle contractions. From the first dataset, the mean estimated depth of the motor unit centres were 8.7, 11.6, and 9.1 mm, with standard deviations 0.5, 0.1, and 1.3 mm, and the respective depths of the fine wire electrodes were 8.4, 15.8, and 9.1 mm. The second dataset generated distinct spatial distributions of motor unit activity which were used to identify the regions of different muscles of the forearm, in a 3-dimensional and projected 2-dimensional view. In conclusion, a method is presented which estimates motor unit centre positions from HDsEMG. The study demonstrates the shifting spatial distribution of muscle fibre activity between different efforts, which could be used to assess individual muscles on a motor unit level.