Optimized Wavelets for Blind Separation of Nonstationary Surface Myoelectric Signals

Surface electromyography (EMG) signals detected over the skin surface may be mixtures of signals generated by many active muscles due to poor spatial selectivity of the recording. In this paper, we propose a new method for blind source separation (BSS) of nonstationary signals modeled as linear instantaneous mixtures. The method is based on whitening of the observations and rotation of the whitened observations. The rotation is performed by joint diagonalization of a set of spatial wavelet distributions (SWDs). The SWDs depend on the selection of the mother wavelet which can be defined by unconstrained parameters via the lattice parameterization within the multiresolution analysis framework. As the sources are classically supposed to be mutually uncorrelated, the design parameters of the mother wavelet can be blindly optimized by minimizing the average (over time lags) cross correlation between the estimated sources. The method was tested on simulated and experimental surface EMG signals and results were compared with those obtained with spatial time-frequency distributions and with second-order statistics (only spectral information). On a set of simulated signals, for 10-dB signal-to-noise ratio (SNR), the cross-correlation coefficient between original and estimated sources was 0.92plusmn0.07 with wavelet optimization, 0.74plusmn0.09 with the wavelet leading to the poorest performance, 0.85plusmn0.07 with Wigner-Ville distribution, 0.86plusmn0.07 with Choi-Williams distribution, and 0.73plusmn0.05 with second-order statistics. In experimental conditions, when the flexor carpi radialis and pronator teres were concomitantly active for 50% of the time, crosstalk was 55.2plusmn10.0% before BSS and was reduced to 15.2plusmn6.3% with wavelet optimization, 30.1plusmn15.0% with the worst wavelet, 28.3plusmn12.3% with Wigner-Ville distribution, 26.2plusmn12.0% with Choi-Williams distribution, and 35.1plusmn15.5% with second-order statistics. In conclusion, the proposed approach resulted in better performance than previous methods for the separation of nonstationary myoelectric signals.

[1]  Abdelhak M. Zoubir,et al.  Joint anti-diagonalization for blind source separation , 2001, 2001 IEEE International Conference on Acoustics, Speech, and Signal Processing. Proceedings (Cat. No.01CH37221).

[2]  Dario Farina,et al.  Blind separation of linear instantaneous mixtures of nonstationary surface myoelectric signals , 2004, IEEE Transactions on Biomedical Engineering.

[3]  R N Lemon,et al.  A novel algorithm to remove electrical cross‐talk between surface EMG recordings and its application to the measurement of short‐term synchronisation in humans , 2002, The Journal of physiology.

[4]  E. Henneman.  Relation between size of neurons and their susceptibility to discharge. , 1957, Science.

[5]  Marie-Françoise Lucas,et al.  Signal-dependent wavelets for electromyogram classification , 2006, Medical and Biological Engineering and Computing.

[6]  Damjan Zazula,et al.  Single autoterms selection for blind source separation in time-frequency plane , 2002, 2002 11th European Signal Processing Conference.

[7]  W. Lawton.  Tight frames of compactly supported affine wavelets , 1990 .

[8]  Moeness G. Amin,et al.  Blind source separation based on time-frequency signal representations , 1998, IEEE Trans. Signal Process..

[9]  Lang Tong,et al.  Indeterminacy and identifiability of blind identification , 1991 .

[10]  C. Torrence,et al.  A Practical Guide to Wavelet Analysis. , 1998 .

[11]  C. Burrus,et al.  Introduction to Wavelets and Wavelet Transforms: A Primer , 1997 .

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

[13]  Cédric Févotte,et al.  Two contributions to blind source separation using time-frequency distributions , 2004, IEEE Signal Processing Letters.

[14]  P. Vaidyanathan.  Multirate Systems And Filter Banks , 1992 .

[15]  L.H. Lindstrom,et al.  Interpretation of myoelectric power spectra: A model and its applications , 1977, Proceedings of the IEEE.

[16]  Marie-Françoise Lucas,et al.  Optimal wavelets for biomedical signal compression , 2006, Medical and Biological Engineering and Computing.

[17]  Wayne Lawton,et al.  Multiresolution properties of the wavelet Galerkin operator , 1991 .

[18]  R. Cantello,et al.  Conduction velocity along human muscle fibers in situ , 1983, Neurology.

[19]  Dario Farina,et al.  A surface EMG generation model with multilayer cylindrical description of the volume conductor , 2004, IEEE Transactions on Biomedical Engineering.

[20]  W. Lawton.  Necessary and sufficient conditions for constructing orthonormal wavelet bases , 1991 .

[21]  Cédric Févotte,et al.  Approche temps-fréquence pour la séparation aveugle de sources non-stationnaires , 2003 .

[22]  J A Stephens,et al.  Correlation between the discharges of motor units recorded from the same and from different finger muscles in man. , 1991, The Journal of physiology.

[23]  D. Winter,et al.  Models of recruitment and rate coding organization in motor-unit pools. , 1993, Journal of neurophysiology.

[24]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[25]  Robert Bregovic,et al.  Multirate Systems and Filter Banks , 2002 .

[26]  C. D. De Luca,et al.  Surface myoelectric signal cross-talk among muscles of the leg. , 1988, Electroencephalography and clinical neurophysiology.