Simulation Techniques in Electromyography

A motor unit action potential (MUAP) recorded in clinical electromyography (EMG) is the spatial and temporal summation of the action potentials (AP's) from all muscle fibers in a motor unit (MU). An important determinant of MUAP waveform characteristics is the size of the recording electrode. In this paper, we have described the use of a modified line source model of single muscle fiber action potentials to simulate MUAP's as recorded by single fiber (SF) EMG, concentric needle (CN) EMG, and macro-EMG electrodes. Results indicate that SFEMG recordings from a normal MU contain mainly the AP's of the closest one to three muscle fibers of the MU. The amplitude, area, and duration of the simulated CNEMG MUAP's are determined mainly by the number and size of muscle fibers within a semicircular territory of 0.5, 1.5, and 2.5 mm, respectively, around the tip of the electrode. The amplitude and area of simulated macro-EMG MUAP's increase with the number of muscle fibers in the MU.     

A technique to track individual motor unit action potentials in surface EMG by monitoring their conduction velocities and amplitudes

The speed of propagation of an action potential along a muscle fiber, its conduction velocity (CV), can be used as an indication of the physiological or pathological state of the muscle fiber membrane. The motor unit action potential (MUAP), the waveform resulting from the spatial and temporal summation of the individual muscle fiber action potentials of that motor unit (MU), propagates with a speed referred to as the motor unit conduction velocity (MUCV). This paper introduces a new algorithm, the MU tracking algorithm, which estimates MUCVs and MUAP amplitudes for individual MUs in a localized MU population using SEMG signals. By tracking these values across time, the electrical activity of the localized MU pool can be monitored. An assessment of the performance of the algorithm has been achieved using simulated SEMG signals. It is concluded that this analysis technique enhances the suitability of SEMG for clinical applications and points toward a future of noninvasive diagnosis and assessment of neuromuscular disorders.     

A model for the generation of synthetic intramuscular EMG signals to test decomposition algorithms

As more and more intramuscular electromyogram (EMG) decomposition programs are being developed, there is a growing need for evaluating and comparing their performances. One way to achieve this goal is to generate synthetic EMG signals having known features. Features of interest are: the number of channels acquired (number of detection surfaces), the number of detected motor unit action potential (MUAP) trains, their time-varying firing rates, the degree of shape similarity among MUAPs belonging to the same motor unit (MU) or to different MUs, the degree of MUAP superposition, the MU activation intervals, the amount and type of additive noise. A model is proposed to generate one or more channels of intramuscular EMG starting from a library of real MUAPs represented in a 16-dimensional space using their Associated Hermite expansion. The MUAP shapes, regularity of repetition rate, degree of superposition, activation intervals, etc. may be time variable and are described quantitatively by a number of parameters which define a stochastic process (the model) with known statistical features. The desired amount of noise may be added to the synthetic signal which may then be processed by the decomposition algorithm under test to evaluate its capability of recovering the signal features.     

Selectivity of spatial filters for surface EMG detection from the tibialis anterior muscle

Many spatial filters have been proposed for surface electromyographic (EMG) signal detection. Although theoretical and modeling predictions on spatial selectivity are available, there are no extensive experimental validations of these techniques based on single motor unit (MU) activity detection. The aim of this study was to compare spatial selectivity of one- and two-dimensional (1-D and 2-D) spatial filters for EMG signal detection. Intramuscular and surface EMG signals were recorded from the tibialis anterior muscle of ten subjects. The simultaneous use of intramuscular wire and surface recordings (with the spike triggered averaging technique) allowed investigation of the activity of single MUs at the skin surface. The surface EMG signals were recorded with a grid of point electrodes (3 x 3 electrodes) and a ring electrode system at 15 locations over the muscle, with the wires detecting signals from the same intramuscular location. For most subjects, it was possible to classify, from the intramuscular recordings, the activity of the same MUs for all the contractions. The surface EMG signals were averaged with the intramuscularly detected MU action potentials as triggers. In this way, eight spatial filters--longitudinal and transversal, single and double differential (LSD, TSD, LDD, TDD), Laplacian (NDD), inverse binomial filter of the second order (IB2), inverse rectangle filter (IR), and differential ring system (C1)--could be compared on the basis of their spatial selectivity. The distance from the source (transversal with respect to the muscle fiber orientation) after which the surface detected potential did not exceed +/- 5% of the maximal peak-to-peak amplitude (detection distance) was statistically smaller for the 2-D systems and TDD than for the other filters. The MU action potential duration was significantly shorter with LDD and with the 2-D systems than with the other filters. The 2-D filters investigated (including C1) showed very similar performance and were, thus, considered equivalent from the point of view of spatial selectivity.     

A software package for the decomposition of long-term multichannel EMG signals using wavelet coefficients

This paper presents a method to decompose multichannel long-term intramuscular electromyogram (EMG) signals. In contrast to existing decomposition methods which only support short registration periods or single-channel recordings of signals of constant muscle effort, the decomposition software EMG-LODEC (ElectroMyoGram LOng-term DEComposition) is especially designed for multichannel long-term recordings of signals of slight muscle movements. A wavelet-based, hierarchical cluster analysis algorithm estimates the number of classes [motor units (MUs)], distinguishes single MUAPs from superpositions, and sets up the shape of the template for each class. Using three channels and a weighted averaging method to track action potential (AP) shape changes improve the analysis. In the last step, nonclassified segments, i.e., segments containing superimposed APs, are decomposed into their units using class-mean signals. Based on experiments on simulated and long-term recorded EMG signals, our software is capable of providing reliable decompositions with satisfying accuracy. EMG-LODEC is suitable for the study of MU discharge patterns and recruitment order in healthy subjects and patients during long-term measurements.     

Blind separation of linear instantaneous mixtures of nonstationary surface myoelectric signals

Electromyographic (EMG) recordings detected over the skin may be mixtures of signals generated by different active muscles due to the phenomena related to volume conduction. Separation of the sources is necessary when single muscle activity has to be detected. Signals generated by different muscles may be considered uncorrelated but in general overlap in time and frequency. Under certain assumptions, mixtures of surface EMG signals can be considered as linear instantaneous but no a priori information about the mixing matrix is available when different muscles are active. In this study, we applied blind source separation (BSS) methods to separate the signals generated by two active muscles during a force-varying task. As the signals are non stationary, an algorithm based on spatial time-frequency distributions was applied on simulated and experimental EMG signals. The experimental signals were collected from the flexor carpi radialis and the pronator teres muscles which could be activated selectively for wrist flexion and rotation, respectively. From the simulations, correlation coefficients between the reference and reconstructed sources were higher than 0.85 for signals largely overlapping both in time and frequency and for signal-to-noise ratios as low as 5 dB. The Choi-Williams and Bessel kernels, in this case, performed better than the Wigner-Ville one. Moreover, the selection of time-frequency points for the procedure of joint diagonalization used in the BSS algorithm significantly influenced the results. For the experimental signals, the interference of the other source in each reconstructed source was significantly attenuated by the application of the BSS method. The ratio between root-mean-square values of the signals from the two sources detected over one of the muscles increased from (mean +/- standard deviation) 2.33 +/- 1.04 to 4.51 +/- 1.37 and from 1.55 +/- 0.46 to 2.72 +/- 0.65 for wrist flexion and rotation, respectively. This increment was statistically significant. It was concluded that the BSS approach applied is promising for the separation of surface EMG signals, with applications ranging from muscle assessment to detection of muscle activation intervals, and to the control of myoelectric prostheses.     

Investigation of optimum electrode locations by using an automatized surface electromyography analysis technique.

Identification of the innervation zone is widely used to optimize the accuracy and precision of noninvasive surface electromyography (EMG) signals because the EMG signal is strongly influenced by innervation zones. However, simply structured fusiform muscle, such as biceps brachii muscle, has been employed mainly due to the simplicity with which the propagation from raw EMG signals can be observed. In this study, the optimum electrode location (OEL), free from innervational influence, was investigated by the propagation pattern of action potentials for brachii muscles and more complicated deltoid muscle structures using an automatized signal analysis technique. The technique employed newly developed computer software with additional clinical uses and minimized subjective differences. EMG signals were recorded using surface array electrodes during voluntary isometric contractions obtained from 12 healthy male subjects. Peaks in EMG signals were detected and averaged for each muscle. The propagation patterns and OEL were examined from biceps brachii muscles for all subjects and from deltoid muscles for seven subjects. The estimated locations were partially confirmed by comparing the root mean squares of the EMG signals. These results show that propagation patterns and OEL could be estimated simply and automatically even from the surface EMG signals of deltoid muscles.     

Simulation of surface EMG signals for a multilayer volume conductor with a superficial bone or blood vessel

This study analytically describes surface electromyogram (EMG) signals generated by a planar multilayer volume conductor constituted by different subdomains modeling muscle, bone (or blood vessel), fat, and skin tissues. The bone is cylindrical in shape, with a semicircular section. The flat portion of the boundary of the bone subdomain is interfaced with the fat layer tissue, the remaining part of the boundary is in contact with the muscle layer. The volume conductor is a model of physiological tissues in which the bone is superficial, as in the case of the tibia bone, backbone, and bones of the forearm. The muscle fibers are considered parallel to the axes of the bone, so that the model is space invariant in the direction of propagation of the action potential. The proposed model, being analytical, allows faster simulations of surface EMG with respect to previously developed models including bone or blood vessels based on the finite-element method. Surface EMG signals are studied by simulating a library of single-fiber action potentials (SFAP) of fibers in different locations within the muscle domain, simulating the generation, propagation, and extinction of the action potential. The decay of the amplitude of the SFAPs in the direction transversal to the fibers is assessed. The decay in the direction of the bone has a lower rate with respect to the opposite direction. Similar results are obtained by simulating motor unit action potentials (MUAPs) constituted by 100 fibers with territory 5 mm2. M waves and interference EMG signals are also simulated based on the library of SFAPs. Again, the decay of the amplitude of the simulated interference EMG signals is lower approaching the bone with respect to going farther from it. The findings of this study indicate the effect of a superficial bone in enhancing the EMG signals in the transversal direction with respect to the fibers of the considered muscle. This increases the effect of crosstalk. The same mathematical method used to simulate a superficial bone can be applied to simulate other physiological tissues. For example, superficial blood vessels (e.g., basilic vein, brachial artery) can influence the recorded EMG signals. As the electrical conductivity of blood is high (it is of the same order as the longitudinal conductivity in the muscle), the effect on EMG signals is opposite compared to the effect of a superficial bone.     

Estimation of muscle fiber conduction velocity with a spectral multidip approach

We propose a novel method for estimation of muscle fiber conduction velocity from surface electromyographic (EMG) signals. The method is based on the regression analysis between spatial and temporal frequencies of multiple dips introduced in the EMG power spectrum through the application of a set of spatial filters. This approach leads to a closed analytical expression of conduction velocity as a function of the auto- and cross-spectra of monopolar signals detected along the direction of muscle fibers. The performance of the algorithm was compared with respect to that of the classic single dip approach on simulated and experimental EMG signals. The standard deviation of conduction velocity estimates from simulated single motor unit action potentials was reduced from 1.51 m/s [10 dB signal-to-noise ratio (SNR)] and 1.06 m/s (20 dB SNR) with the single dip approach to 0.51 m/s (10 dB) and 0.23 m/s (20 dB) with the proposed method using 65 dips. When 200 active motor units were simulated in an interference EMG signal, standard deviation of conduction velocity decreased from 0.95 m/s (10 dB SNR) and 0.60 m/s (20 dB SNR) with a single dip to 0.21 m/s (10 dB) and 0.11 m/s (20 dB) with 65 dips. In experimental signals detected from the abductor pollicis brevis muscle, standard deviation of estimation decreased from (mean +/- SD over 5 subjects) 1.25 +/- 0.62 m/s with one dip to 0.10 +/- 0.03 m/s with 100 dips. The proposed method does not imply limitation in resolution of the estimated conduction velocity and does not require any iterative procedure for the estimate since it is based on a closed analytical formulation.     

Modeling of surface myoelectric signals. I. Model implementation

The relationships between the parameters of active motor units (MU's) and the features of surface electromyography (EMG) signals have been investigated using a mathematical model that represents the surface EMG as a summation of contributions from the single muscle fibers. Each MU has parallel fibers uniformly scattered within a cylindrical volume of specified radius embedded in an anisotropic medium. Two action potentials, each modeled as a current tripole, are generated at the neuromuscular junction, propagate in opposite directions and extinguish at the fiber-tendon endings. The neuromuscular junctions and fiber-tendon endings are uniformly scattered within regions of specified width. Muscle fiber conduction velocity and average fiber length to the right and left of the center of the innervation zone are also specified. The signal produced by MU's with different geometries and conduction velocities are superimposed. Monopolar, single differential and double differential signals are computed from electrodes placed in equally spaced locations on the surface of the muscle and are displayed as functions of any of the model's parameters. Spectral and amplitude variables and conduction velocity are estimated from the surface signals and displayed as functions of any of the model's parameters. The influence of fiber-end effects, electrode misalignment, tissue anisotropy, MU's location and geometry are discussed. Part II of this paper will focus on the simulation and interpretation of experimental signals.     

Fatigue compensation of the electromyographic signal for prostheticcontrol and force estimation

During a sustained muscle contraction, the amplitude of electromyographic (EMG) signals increases and the spectrum of the EMG signal shifts toward lower frequencies. These effects are due to muscular fatigue and can cause problems in the control of myoelectric prostheses and in the estimation of contraction level from the EMG signal. It has been well known that the fatigue effects can be explained by the conduction velocity changes during the fatigue process and by the idea that the conduction velocity is linearly proportional to the median frequency of EMG signals. Hence the fatigue process can be monitored by measuring the median frequency. A fatigue compensation preprocessor has been developed. It uses the widely accepted power spectrum density model of EMG signals that contains the conduction velocity as a measure of fatigue. It was verified that the preprocessor scales down the amplitude of the fatigued EMG signal and decompresses the spectrum. Hence, the preprocessor eliminates the increase in amplitude and the shift in frequency and enables consistent EMG signals to be used to control prostheses     

Evaluation of the Effectiveness of EMG Parameters in the Study of Neurogenic Diseases-A Statistical Approach Using Clinical and Simulated Data

Of interest here is the problem of determining to what extent combinations of parameters derived from the EMG signal allow 1) discriminating two subclasses of neurogenic myopathies, and 2) recognizing different morphologies of the motor unit action potential underlying a measured EMG signal. EMG signals measured on clinical subjects and computer-simulated EMG signals were collected in a database and used cooperatively in this study. Suitable statistical models were developed which allow testing hypotheses on the role of accepted EMG parameters for the two purposes named above, and deriving new suitable combinations of EMG parameters. Results support the hypothesis that frequency-domain parameters are very clearly related to the morphology of the motor unit action potential. However, the attempt to use them in order to discriminate the two pathologic subclasses considered appears to be jeopardized by the fact that the signal may be measured in territories which do not reflect the morphology of the motor unit action potential dominant in such subclasses. On the basis of time-domain parameters, a significant discrimination was obtained between the two subclasses, and such discrimination is related mainly to a time-domain parameter which has already proved successful in the discrimination between myopathic and normal subjects. Data corroborate the hypothesis that the diagnostic yield improves when time-domain EMG parameters are measured at recruitment.     

Design of perturbation signals for the estimation of proprioceptive reflexes

This study aimed to identify the functional contribution of reflexes to human motor control during posture maintenance. Continuous random force disturbances were applied at the hand while the subjects were instructed to minimize the deviation resulting from the force disturbances. The results were analyzed in the frequency domain with frequency response functions (FRFs). Two FRFs were evaluated: 1) the mechanical admittance and 2) the reflexive impedance, expressing the dynamic relation between position and muscle activation (assessed via electromyography, EMG). The reflexive impedance is a direct measure of the proprioceptive reflexes. To record all relevant dynamical characteristics of the arm, wide bandwidth signals were used as force disturbance. Distributing the power of the signal over fewer frequencies within the bandwidth improved the signal-to-noise-ratio SNR of the EMG recordings, facilitating reliable estimation of the reflexive impedance. The coherence indicated that the relation between force disturbance and EMG is linear under the given conditions and improved with the SNR. The method of designing disturbance signals and the estimation of the reflexive impedance are useful for studies aiming to quantify proprioceptive reflexes and to investigate its functionality.     

Spatial Filtering of Noninvasive Multielectrode EMG: Part I?Introduction to Measuring Technique and Applications

Complementary to its conventional applications, surface EMG is also suited to gain more detailed information on the functional state of a muscle, when measurement configurations with smaller pickup areas are used. A new category of suitable measurement configurations is obtained by application of the spatial filtering principle to electromyography. In a spatial filter unit, the signals of several recording electrodes are combined to form one output signal channel. The filter characteristic is determined by the weighting factors used and by the geometrical arrangement of the electrodes. Extended multielectrode arrays and multichannel recording make possible the detection of correlated excitations at different sites of the muscle. Even in high levels of muscle contraction, single motor unit impulses that are suitably shaped by filtering can be repeatedly recognized in the surface EMG signal. In clinical studies, pathologically shaped impulses have been identified indicating multiple innervation zones. The initiation and the propagation of excitation within single motor units can be detected with improved accuracy even from very small muscles.     

Autoregressive Modeling of Surface EMG and Its Spectrum with Application to Fatigue

The following is an investigation of the ability of the autoregressive (AR) model to describe the spectrum of the processes underlying the recorded surface EMG. Surface EMG (SEMG) spectrum is influenced by two major factors; one attributed to the motor units (MU) firing rate and the second, the higher frequency one, to the morphology of the action potentials (AP) traveling along the muscle fiber. In the present paper, SEMG measurements were carried out on the biceps brachii muscle with fixed surface electrodes arrangement and isotonic conditions. Sufficient averaging of 0.5 s segments enabled the identification of the low-frequency peak related to the firing rates of the MU's.     

Differences in the activity of the muscles in the forearm of individuals with a congenital absence of the hand

The spectral content of the myoelectric signals from the muscles of the remnant forearms of three persons with congenital absences (CA) of their forearms was compared with signals from their intact contra-lateral limbs, similar muscles in three persons with acquired losses (AL) and seven persons without absences [no loss (NL)]. The observed bandwidth for the CA subjects was broader with peak energy between 200 and 300 Hz. While the signals from the contra-lateral limbs and the AL and NL subjects was in the 100-150 Hz range. The mean skew of the signals from the AL subjects was 46.3 +/- 6.7 and those with NL of 45.4 +/- 8.7, while the signals from those with CAs had a skew of 11.0 +/- 11. The structure of the muscles of one CA subject was observed ultrasonically. The muscle showed greater disruption than normally developed muscles. It is speculated that the myographic signal reflects the structure of the muscle which has developed in a more disorganized manner as a result of the muscle not being stretched by other muscles across the missing distal joint, even in the muscles that are used regularly to control arm prostheses.     

Short time Fourier analysis of the electromyogram: fast movements and constant contraction

We applied short-time Fourier analysis to surface electromyograms (EMG) recorded during rapid movements, and during isometric contractions at constant forces. We selected a portion of the data to be transformed by multiplying the signal by a Hamming window, then computed the discrete Fourier transform. Shifting the window along the data record, we computed a new spectrum each 10 ms. We displayed the transformed data in spectograms or "voiceprints." This short-time technique allowed us to see time-dependencies in the EMG that are normally averaged in the Fourier analysis of these signals. Spectra of EMG's during isometric contractions at constant force vary in the short (10-20 ms) term. Moments of the spectral distribution show this variability. Short-time spectra from EMG's recorded during rapid movements were much less variable. The windowing technique picked out the typical "three-burst pattern" in EMG's from both wrist and head movements. Spectra during the bursts were more consistent than those during isometric contractions. Furthermore, there was a consistent shift in spectral statistics in the course of the three bursts. Both the center frequency and the variance of the spectral energy distribution grew from the first burst to the second burst in the same muscle. We discuss this pattern with respect to the origin of the EMG bursts in rapid movement. We also extend the analogy between electromyograms and speech signals to argue for future applicability of short-time spectral analysis of EMG.     

A fast and reliable technique for muscle activity detection from surface EMG signals

The estimation of on-off timing of human skeletal muscles during movement is an important issue in surface electromyography (EMG) signal processing with relevant clinical applications. In this paper, a novel approach to address this issue is proposed. The method is based on the identification of single motor unit action potentials from the surface EMG signal with the use of the continuous wavelet transform. A manifestation variable is computed as the maximum of the outputs of a bank of matched filters at different scales. A threshold is applied to the manifestation variable to detect EMG activity. A model, based on the physical structure of the muscle, is used to test the proposed technique on synthetic signals with known features. The resultant bias of the onset estimate is lower than 40 ms and the standard deviation lower than 30 ms in case of additive colored Gaussian noise with signal-to-noise ratio as low as 2 dB. Comparison with previously developed methods was performed, and representative applications to experimental signals are presented. The method is designed for a complete real-time implementation and, thus, may be applied in clinical routine activity.     

Spatial filtering for robust myoelectric control

Pattern recognition techniques have been applied to extract information from electromyographic (EMG) signals that can be used to control electrical powered hand prostheses. In this paper, optimized spatial filters that enhance separation properties of EMG signals are investigated. In particular, different multiclass extensions of the common spatial patterns algorithm are applied to high-density surface EMG signals acquired from the forearms of ten healthy subjects. Visualization of the obtained filter coefficients provides insight into the physiology of the muscles related to the performed contractions. The CSP methods are compared with a commonly used pattern recognition approach in a six-class classification task. Cross-validation results show a significant improvement in performance and a higher robustness against noise than commonly used pattern recognition methods.     

EMG-based speech recognition using hidden markov models with global control variables

It is well known that a strong relationship exists between human voices and the movement of articulatory facial muscles. In this paper, we utilize this knowledge to implement an automatic speech recognition scheme which uses solely surface electromyogram (EMG) signals. The sequence of EMG signals for each word is modelled by a hidden Markov model (HMM) framework. The main objective of the work involves building a model for state observation density when multichannel observation sequences are given. The proposed model reflects the dependencies between each of the EMG signals, which are described by introducing a global control variable. We also develop an efficient model training method, based on a maximum likelihood criterion. In a preliminary study, 60 isolated words were used as recognition variables. EMG signals were acquired from three articulatory facial muscles. The findings indicate that such a system may have the capacity to recognize speech signals with an accuracy of up to 87.07%, which is superior to the independent probabilistic model.