Optimal spatio-temporal filter for the reduction of crosstalk in surface electromyogram

Nonlinear spatio-temporal filter to reduce crosstalk in bipolar electromyogram

Objective.The wide detection volume of surface electromyogram (EMG) makes it prone to crosstalk, i.e. the signal from other muscles than the target one. Removing this perturbation from bipolar recordings is an important open problem for many applications.Approach.An innovative nonlinear spatio-temporal filter is developed to estimate the EMG generated by the target muscle by processing noisy signals from two bipolar channels, placed over the target and the crosstalk muscle, respectively. The filter is trained on some calibration data and then can be applied on new signals. Tests are provided in simulations (considering different thicknesses of the subcutaneous tissue, inter-electrode distances, locations of the EMG channels, force levels) and experiments (from pronator teres and flexor carpi radialis of 8 healthy subjects).Main results.The proposed filter allows to reduce the effect of crosstalk in all investigated conditions, with a statistically significant reduction of its root mean squared of about 20%, both in simulated and experimental data. Its performances are also superior to those of a blind source separation method applied to the same data.Significance.The proposed filter is simple to be applied and feasible in applications in which single bipolar channels are placed over the muscles of interest. It can be useful in many fields, such as in gait analysis, tests of myoelectric fatigue, rehabilitation with EMG biofeedback, clinical studies, prosthesis control.

Comparison of Blind Source Separation Methods to Surface Electromyogram for Extensor Muscles of the Index and Little Fingers

Crosstalk is the result of the propagation of muscle electrical signals on surface electromyogram channels simultaneously. The objective of this paper is to study the behavior of three blind source separation (BSS) methods for crosstalk reduction during finger extensor muscle contractions: FastICA, joint diagonalization of covariance matrices and optimal filtering. These methods have been tested on artificial mixtures defined by a temporal sum of the real signals from isolated contraction of two independent biomechanical muscles for the extension of the index and little finger. Artificial mixtures display a ground truth for comparison between the methods. The separation was better using the optimal filtering compared to the other two methods. The optimal filtering have then be tested on real mixtures recorded during a simultaneous contraction of the two muscles. The results are less satisfactory but open doors to new perspectives.

Analytical Modelling of Surface EMG Signals Generated by Curvilinear Fibers with Approximate Conductivity Tensor

OBJECTIVE

Mathematical modelling of surface electromyographic (EMG) signals has been proven a valuable tool to interpret experimental data and to validate signal processing techniques. Most analytical EMG models only consider muscle fibers with specific and fixed arrangements. However, the fiber orientation and curvature may change along the fiber paths and may differ from fiber to fiber. Here we propose a subject-specific EMG model that simulates the fiber trajectories in muscles of the upper arm and analytically derives the action potentials assuming an approximate conductivity tensor.

METHODS

Magnetic Resonance (MR) images were acquired to identify and generate muscle fiber paths and to determine the muscle locations in a cylindrical volume conductor. While the propagation of the action potentials followed the identified curvilinear fiber paths, the conductivity tensor was not adapted to the fiber direction but approximated along the longitudinal axis of the cylindrical volume conductor. Single fiber action potentials (SFAPs) were computed by simulating the generation, propagation, and extinction of membrane current sources. To validate the assumption of the approximate conductivity tensor, two numerical models were implemented for comparison with the analytical solution. The first numerical model reproduced the analytical model and therefore included an approximation for the conductivity tensor. The second numerical model included the exact conductivity tensor derived from the fiber curvatures.

RESULTS

The motor unit action potentials generated by the proposed analytical model and the two numerical models were highly similar (cross-correlation >0.98, normalized root mean square error, nRMSE 0.04, relative error in the median frequency of the simulated waveforms of approximately 3%). The proposed analytical model was also evaluated by comparing simulated and experimentally recorded compound muscle action potentials (CMAPs). The CMAPs simulated with the proposed model better matched the experimental data (cross-correlation >0.90 and nRMSE <0.25 for the majority of the channels) than a model with straight fibers. Finally, the proposed model was representatively used to test the accuracy of an EMG decomposition algorithm, providing a realistic benchmark.

CONCLUSIONS AND SIGNIFICANCE

The proposed analytical model generates action potentials that reflect the spatial distributions of muscle fibers with curvilinear paths. The simulated signals are more realistic than signals generated by analytical models with straight fibers and can therefore be applied for testing EMG processing algorithms with a trade-off between simulation accuracy and computational speed.

Inverse modelling to reduce crosstalk in high density surface electromyogram

Surface electromyogram (EMG) has a relatively large detection volume, so that it could include contributions both from the target muscle of interest and from nearby regions (i.e., crosstalk). This interference can prevent a correct interpretation of the activity of the target muscle, limiting the use of surface EMG in many fields. To counteract the problem, selective spatial filters have been proposed, but they reduce the representativeness of the data from the target muscle. A better solution would be to discard only crosstalk from the signal recorded in monopolar configuration (thus, keeping most information on the target muscle). An inverse modelling approach is here proposed to estimate the contributions of different muscles, in order to focus on the one of interest. The method is tested with simulated monopolar EMGs from superficial nearby muscles contracted at different force levels (either including or not model perturbations and noise), showing statistically significant improvements in information extraction from the data. The median over the entire dataset of the mean squared error in representing the EMG of the muscle under the detection electrode was reduced from 11.2% to 4.4% of the signal energy (5.3% if noisy data were processed); the median bias in conduction velocity estimation (from 3 monopolar channels aligned to the muscle fibres) was decreased from 2.12 to 0.72 m/s (1.1 m/s if noisy data were processed); the median absolute error in the estimation of median frequency was reduced from 1.02 to 0.67 Hz in noise free conditions and from 1.52 to 1.45 Hz considering noisy data.

Pathway-specific modulatory effects of neuromuscular electrical stimulation during pedaling in chronic stroke survivors

Background

Neuromuscular electrical stimulation (NMES) is extensively used in stroke motor rehabilitation. How it promotes motor recovery remains only partially understood. NMES could change muscular properties, produce altered sensory inputs, and modulate fluctuations of cortical activities; but the potential contribution from cortico-muscular couplings during NMES synchronized with dynamic movement has rarely been discussed.

Method

We investigated cortico-muscular interactions during passive, active, and NMES rhythmic pedaling in healthy subjects and chronic stroke survivors. EEG (128 channels), EMG (4 unilateral lower limb muscles) and movement parameters were measured during 3 sessions of constant-speed pedaling. Sensory-level NMES (20 mA) was applied to the muscles, and cyclic stimulation patterns were synchronized with the EMG during pedaling cycles. Adaptive mixture independent component analysis was utilized to determine the movement-related electro-cortical sources and the source dipole clusters. A directed cortico-muscular coupling analysis was conducted between representative source clusters and the EMGs using generalized partial directed coherence (GPDC). The bidirectional GPDC was compared across muscles and pedaling sessions for post-stroke and healthy subjects.

Results

Directed cortico-muscular coupling of NMES cycling was more similar to that of active pedaling than to that of passive pedaling for the tested muscles. For healthy subjects, sensory-level NMES could modulate GPDC of both ascending and descending pathways. Whereas for stroke survivors, NMES could modulate GPDC of only the ascending pathways.

Conclusions

By clarifying how NMES influences neuromuscular control during pedaling in healthy and post-stroke subjects, our results indicate the potential limitation of sensory-level NMES in promoting sensorimotor recovery in chronic stroke survivors.

Association of anthropometric parameters with amplitude and crosstalk of mechanomyographic signals during forearm flexion, pronation and supination torque tasks

This study aimed to quantify the association of four anthropometric parameters of the human arm, namely, the arm circumference (CA), arm length (LA), skinfold thickness (ST) and inter-sensor distance (ISD), with amplitude (RMS) and crosstalk (CT) of mechanomyography (MMG) signals. Twenty-five young, healthy, male participants were recruited to perform forearm flexion, pronation and supination torque tasks. Three accelerometers were employed to record the MMG signals from the biceps brachii (BB), brachialis (BRA) and brachioradialis (BRD) at 80% maximal voluntary contraction (MVC). Signal RMS was used to quantify the amplitude of the MMG signals from a muscle, and cross-correlation coefficients were used to quantify the magnitude of the CT among muscle pairs (BB & BRA, BRA & BRD, and BB & BRD). For all investigated muscles and pairs, RMS and CT showed negligible to low negative correlations with CA, LA and ISD (r = −0.0001–−0.4611), and negligible to moderate positive correlations with ST (r = 0.004–0.511). However, almost all of these correlations were statistically insignificant (p > 0.05). These findings suggest that RMS and CT values for the elbow flexor muscles recorded and quantified using accelerometers appear invariant to anthropometric parameters.

Tutorial. Surface EMG detection in space and time: Best practices

This tutorial is aimed to non-engineers using, or planning to use, surface electromyography (sEMG) as an assessment tool in the prevention, monitoring and rehabilitation fields. Its first purpose is to address the issues related to the origin and nature of the signal and to its detection (electrode size, distance, location) by one-dimensional (bipolar and linear arrays) and two-dimensional (grids) electrode systems while avoiding advanced mathematical, physical or physiological issues. Its second purpose is to outline best practices and provide general guidelines for proper signal detection. Issues related to the electrode-skin interface, signal conditioning and interpretation will be discussed in subsequent tutorials.

A review on crosstalk in myographic signals

PURPOSE

Crosstalk in myographic signals is a major hindrance to the understanding of local information related to individual muscle function. This review aims to analyse the problem of crosstalk in electromyography and mechanomyography.

METHODS

An initial search of the SCOPUS database using an appropriate set of keywords yielded 290 studies, and 59 potential studies were selected after all the records were screened using the eligibility criteria. This review on crosstalk revealed that signal contamination due to crosstalk remains a major challenge in the application of surface myography techniques. Various methods have been employed in previous studies to identify, quantify and reduce crosstalk in surface myographic signals.

RESULTS

Although correlation-based methods for crosstalk quantification are easy to use, there is a possibility that co-contraction could be interpreted as crosstalk. High-definition EMG has emerged as a new technique that has been successfully applied to reduce crosstalk.

CONCLUSIONS

The phenomenon of crosstalk needs to be investigated carefully because it depends on many factors related to muscle task and physiology. This review article not only provides a good summary of the literature on crosstalk in myographic signals but also discusses new directions related to techniques for crosstalk identification, quantification and reduction. The review also provides insights into muscle-related issues that impact crosstalk in myographic signals.