Decomposition of Multi-Channel Intramuscular EMG Signals by Cyclostationary-Based Blind Source Separation

We propose a novel decomposition method for electromyographic signals based on blind source separation. Using the cyclostationary properties of motor unit action potential trains (MUAPt), it is shown that the MUAPt can be decomposed by joint diagonalization of the cyclic spatial correlation matrix of the observations. After modeling the source signals, we provide the proof of orthogonality of the sources and of their delayed versions in a cyclostationary context. We tested the proposed method on simulated signals and showed that it can decompose up to six sources with a probability of correct detection and classification >95%, using only eight recording sites. Moreover, we tested the method on experimental multi-channel signals recorded with thin-film intramuscular electrodes, with a total of 32 recording sites. The rate of agreement of the decomposed MUAPt with those obtained by an expert using a validated tool for decomposition was >93%.

[1]  Philippe Ravier,et al.  Cyclostationary analysis of electromyographic signals , 2013, 21st European Signal Processing Conference (EUSIPCO 2013).

[2]  Dario Farina,et al.  Covariance and Time-Scale Methods for Blind Separation of Delayed Sources , 2011, IEEE Transactions on Biomedical Engineering.

[3]  Christian Faye,et al.  Real time detector for cyclostationary RFI in radio astronomy , 1998, 9th European Signal Processing Conference (EUSIPCO 1998).

[4]  George S. Moschytz,et al.  High-precision EMG signal decomposition using communication techniques , 2000, IEEE Trans. Signal Process..

[5]  B Conrad,et al.  A neuroelectric signal recognition system. , 1972, Electroencephalography and clinical neurophysiology.

[6]  P. A. Karthick,et al.  Analysis of Muscle Fatigue Progression using Cyclostationary Property of Surface Electromyography Signals , 2015, Journal of Medical Systems.

[7]  J. L. Coatrieux,et al.  Presentation de methodes de reconnaissance des signaux d'electromyographie , 1981 .

[8]  Damjan Zazula,et al.  Gradient Convolution Kernel Compensation Applied to Surface Electromyograms , 2007, ICA.

[9]  Christine Servière,et al.  New second order cyclostationary analysis and application to the detection and characterization of a runner's fatigue , 2014, Signal Process..

[10]  Philippe Ravier,et al.  A new cyclostationarity-based blind approach for motor unit's firing rate automated detection in electromyographic signals , 2014, IEEE-EMBS International Conference on Biomedical and Health Informatics (BHI).

[11]  Ping Zhou,et al.  A Novel Framework Based on FastICA for High Density Surface EMG Decomposition , 2016, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[12]  Georgios B. Giannakis,et al.  Bibliography on cyclostationarity , 2005, Signal Process..

[13]  Damjan Zazula,et al.  Surface EMG Decomposition Using a Novel Approach for Blind Source Separation , 2003 .

[14]  Damjan Zazula,et al.  Real-Time Motor Unit Identification From High-Density Surface EMG , 2013, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[15]  Kevin C. McGill,et al.  Optimal resolution of superimposed action potentials , 2002, IEEE Transactions on Biomedical Engineering.

[16]  Valentina Agostini,et al.  An Algorithm for the Estimation of the Signal-To-Noise Ratio in Surface Myoelectric Signals Generated During Cyclic Movements , 2012, IEEE Transactions on Biomedical Engineering.

[17]  A Holobar,et al.  Non-invasive characterization of motor unit behaviour in pathological tremor , 2012, Journal of neural engineering.

[18]  C. D. De Luca,et al.  Control scheme governing concurrently active human motor units during voluntary contractions , 1982, The Journal of physiology.

[19]  H. Clamann.  Statistical analysis of motor unit firing patterns in a human skeletal muscle. , 1969, Biophysical journal.

[20]  Ronald S. Lefever,et al.  A Procedure for Decomposing the Myoelectric Signal Into Its Constituent Action Potentials-Part II: Execution and Test for Accuracy , 1982, IEEE Transactions on Biomedical Engineering.

[21]  Christian Jutten,et al.  Detection de grandeurs primitives dans un message composite par une architecture de calcul neuromime , 1985 .

[22]  V J Prochazka,et al.  On-line multi-unit sorting with resolution of superposition potentials. , 1973, Electroencephalography and clinical neurophysiology.

[23]  Hossein Parsaei,et al.  EMG Signal Decomposition Using Motor Unit Potential Train Validity , 2013, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[24]  Pierre Comon,et al.  Independent component analysis, A new concept? , 1994, Signal Process..

[25]  H. Broman.  Knowledge-based signal processing in the decomposition of myoelectric signals , 1988, IEEE Engineering in Medicine and Biology Magazine.

[26]  Kiran Marri,et al.  Analysis of fatigue conditions in triceps brachii muscle using sEMG signals and spectral correlation density function , 2014, 2014 International Conference on Informatics, Electronics & Vision (ICIEV).

[27]  Kevin C. McGill,et al.  EMGLAB: An interactive EMG decomposition program , 2005, Journal of Neuroscience Methods.

[28]  Dario Farina,et al.  Analysis of intramuscular electromyogram signals , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[29]  Antoine Souloumiac,et al.  Jacobi Angles for Simultaneous Diagonalization , 1996, SIAM J. Matrix Anal. Appl..

[30]  Ronald S. Lefever,et al.  A Procedure for Decomposing the Myoelectric Signal Into Its Constituent Action Potentials - Part I: Technique, Theory, and Implementation , 1982, IEEE Transactions on Biomedical Engineering.

[31]  Dario Farina,et al.  Time-varying delay estimators for measuring muscle fiber conduction velocity from the surface electromyogram , 2015, Biomed. Signal Process. Control..

[32]  Kenzo Akazawa,et al.  Decomposition of Synthetic Multi-channel Surface-Electromyogram Using Independent Component Analysis , 2004, ICA.

[33]  Damjan Zazula,et al.  Multichannel Blind Source Separation Using Convolution Kernel Compensation , 2007, IEEE Transactions on Signal Processing.

[34]  Eric Moulines,et al.  A blind source separation technique using second-order statistics , 1997, IEEE Trans. Signal Process..

[35]  Damjan Zazula,et al.  A novel approach to convolutive blind separation of close-to-orthogonal pulse sources using second-order statistics , 2004, 2004 12th European Signal Processing Conference.

[36]  C. D. De Luca,et al.  Behaviour of human motor units in different muscles during linearly varying contractions , 1982, The Journal of physiology.

[37]  Robert B. Randall,et al.  THE RELATIONSHIP BETWEEN SPECTRAL CORRELATION AND ENVELOPE ANALYSIS IN THE DIAGNOSTICS OF BEARING FAULTS AND OTHER CYCLOSTATIONARY MACHINE SIGNALS , 2001 .

[38]  C. Disselhorst-Klug,et al.  Principles of high-spatial-resolution surface EMG (HSR-EMG): single motor unit detection and application in the diagnosis of neuromuscular disorders. , 1997, Journal of electromyography and kinesiology : official journal of the International Society of Electrophysiological Kinesiology.

[39]  Damjan Zazula,et al.  Blind Deconvolution of Close-to-Orthogonal Pulse Sources Applied to Surface Electromyograms , 2004, ICA.

[40]  S. Gandevia,et al.  Accurate and representative decoding of the neural drive to muscles in humans with multi‐channel intramuscular thin‐film electrodes , 2015, The Journal of physiology.

[41]  Qiang Li,et al.  The Decomposition of Surface EMG Signals Based on Blind Source Separation of Convolved Mixtures* , 2005, 2005 IEEE Engineering in Medicine and Biology 27th Annual Conference.

[42]  Kevin C. McGill,et al.  High-Resolution Alignment of Sampled Waveforms , 1984, IEEE Transactions on Biomedical Engineering.

[43]  Dario Farina,et al.  Unsupervised Bayesian Decomposition of Multiunit EMG Recordings Using Tabu Search , 2010, IEEE Transactions on Biomedical Engineering.

[44]  C. J. De Luca,et al.  Probability Distribution Function of the Inter-Pulse Intervals of Single Motor Unit Action Potentials during Isometric Contractions , 1973 .

[45]  Ganesh R. Naik,et al.  Single-Channel EMG Classification With Ensemble-Empirical-Mode-Decomposition-Based ICA for Diagnosing Neuromuscular Disorders , 2016, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[46]  Antonio Napolitano,et al.  Cyclostationarity: Half a century of research , 2006, Signal Process..

[47]  R. Stein,et al.  The relation between the surface electromyogram and muscular force. , 1975, The Journal of physiology.

[48]  B Mambrito,et al.  A technique for the detection, decomposition and analysis of the EMG signal. , 1984, Electroencephalography and clinical neurophysiology.

[49]  Antonio Napolitano,et al.  Cyclostationarity: New trends and applications , 2016, Signal Process..

[50]  C. Sherrington,et al.  Recruitment and some other features of reflex inhibition , 1925 .

[51]  Dario Farina,et al.  Characterization of Human Motor Units From Surface EMG Decomposition , 2016, Proceedings of the IEEE.

[52]  Frédéric Bonnardot,et al.  Comparaison entre les analyses angulaire et temporelle des signaux vibratoires de machines tournantes : étude du concept de cyclostationnarité floue , 2004 .

[53]  J. Vibert,et al.  Spike separation in multiunit records: a multivariate analysis of spike descriptive parameters. , 1979, Electroencephalography and clinical neurophysiology.

[54]  Dario Farina,et al.  The extraction of neural strategies from the surface EMG: an update. , 2014, Journal of applied physiology.

[55]  Kevin C. McGill,et al.  Automatic Decomposition of the Clinical Electromyogram , 1985, IEEE Transactions on Biomedical Engineering.

[56]  D. Farina,et al.  Multi-channel intramuscular and surface EMG decomposition by convolutive blind source separation , 2016, Journal of neural engineering.

[57]  D.W. Stashuk,et al.  Probabilistic inference-based classification applied to myoelectric signal decomposition , 1992, IEEE Transactions on Biomedical Engineering.