Highly Accurate Real-Time Decomposition of Single Channel Intramuscular EMG.

OBJECTIVE: Real-time intramuscular electromyography (iEMG) decomposition, as an identification procedure of individual motor neuron (MN) discharge timings from a streaming iEMG recording, has the potential to be used in human-machine interfacing. However, for these applications, the decomposition accuracy and speed of current approaches need to be improved. METHODS: In our previous work, a real-time decomposition algorithm based on a Hidden Markov Model of EMG, using GPU-implemented Bayesian filter to estimate the spike trains of motor units (MU) and their action potentials (MUAPs), was proposed. In this paper, a substantially extended version of this algorithm that boosts the accuracy while maintaining real-time implementation, is introduced. Specifically, multiple heuristics that aim at resolving the problems leading to performance degradation, are applied to the original model. In addition, the recursive maximum likelihood (RML) estimator previously used to estimate the statistical parameters of the spike trains, is replaced by a linear regression (LR) estimator, which is computationally more efficient, in order to ensure real-time decomposition with the new heuristics. RESULTS: The algorithm was validated using twenty-one experimental iEMG signals acquired from the tibialis anterior muscle of five subjects by fine wire electrodes. All signals were decomposed in real time. The decomposition accuracy depended on the level of muscle activation and was when less than 10 MUs were identified, substantially exceeding previous real-time results. CONCLUSION: Single channel iEMG signals can be very accurately decomposed in real time with the proposed algorithm. SIGNIFICANCE: The proposed highly accurate algorithm for single-channel iEMG decomposition has the potential of providing neural information on motor tasks for human interfacing.

Recursive Decomposition of Electromyographic Signals With a Varying Number of Active Sources: Bayesian Modeling and Filtering.

OBJECTIVE: This paper describes a sequential decomposition algorithm for single-channel intramuscular electromyography (iEMG) generated by a varying number of active motor neurons. METHODS: As in previous work, we establish a hidden Markov model of iEMG, in which each motor neuron spike train is modeled as a renewal process with inter-spike intervals following a discrete Weibull law and motor unit action potentials are modeled as impulse responses of linear time-invariant systems with known prior. We then expand this model by introducing an activation vector associated with the state vector of the hidden Markov model. This activation vector represents recruitment/derecruitment of motor units and is estimated together with the state vector using Bayesian filtering. Non-stationarity of the model parameters is addressed by means of a sliding window approach, thus making the algorithm adaptive to variations in contraction force and motor unit action potential waveforms. RESULTS: The algorithm was validated using simulated and experimental iEMG signals with varying number of active motor units. The experimental signals were acquired from the tibialis anterior and abductor digiti minimi muscles by fine wire and needle electrodes. The decomposition accuracy in both simulated and experimental signals exceeded 90%. CONCLUSION: The recruitment/derecruitment was successfully tracked by the algorithm. Because of its parallel structure, this algorithm can be efficiently accelerated, which lays the basis for its real-time applications in human-machine interfaces. SIGNIFICANCE: The proposed method substantially broadens the domains of applicability of the algorithm.

On-Line Recursive Decomposition of Intramuscular EMG Signals Using GPU-Implemented Bayesian Filtering.

OBJECTIVE: Real-time intramuscular electromyography (iEMG) decomposition, which is needed in biofeedback studies and interfacing applications, is a complex procedure that involves identifying the motor neuron spike trains from a streaming iEMG recording. METHODS: We have previously proposed a sequential decomposition algorithm based on a Hidden Markov Model of EMG, which used Bayesian filter to estimate unknown parameters of motor unit (MU) spike trains, as well as their action potentials (MUAPs). Here, we present a modification of this original model in order to achieve a real-time performance of the algorithm as well as a parallel computation implementation of the algorithm on Graphics Processing Unit (GPU). Specifically, the Kalman filter previously used to estimate the MUAPs, is replaced by a least-mean-square filter. Additionally, we introduce a number of heuristics that help to omit the most improbable decomposition scenarios while searching for the best solution. Then, a GPU-implementation of the proposed algorithm is presented. RESULTS: Simulated iEMG signals containing up to 10 active MUs, as well as five experimental fine-wire iEMG signals acquired from the tibialis anterior muscle, were decomposed in real time. The accuracy of decompositions depended on the level of muscle activation, but in all cases exceeded 85 %. CONCLUSION: The proposed method and implementation provide an accurate, real-time interface with spinal motor neurons. SIGNIFICANCE: The presented real time implementation of the decomposition algorithm substantially broadens the domain of its application.