Single site electromyograph amplitude estimation

Previous investigators have experimentally demonstrated and/or analytically predicted that temporal whitening of the surface electromyograph (EMG) waveform prior to demodulation improves the EMG amplitude estimate. However, no systematic study of the influence of various whitening filters upon amplitude estimate performance has been reported. The authors describe a phenomenological mathematical model of a single site of the surface EMG waveform and reports on experimental studies which examined the performance of several temporal whitening filters. Surface EMG waveforms were sampled during nonfatiguing, constant-force, isometric contractions of the biceps or triceps muscles, over the range of 10-75% maximum voluntary contraction. A signal-to-noise ratio (SNR) was computed from each amplitude estimate (deviations about the mean value of the estimate were considered as noise). A moving average root mean square estimator (245 ms window) provided an average±standard deviation (A±SD) SNR of 10.7±3.3 for the individual recordings. Temporal whitening with one fourth-order whitening filter designed per site improved the A±SD SNR to 17.6±6.0     

Influence of smoothing window length on electromyogram amplitudeestimates

A systematic, experimental study of the influence of smoothing window length on the signal-to-noise ratio (SNR) of electromyogram (EMG) amplitude estimates is described. Surface EMG waveforms were sampled during nonfatiguing, constant-force, constant-angle contractions of the biceps or triceps muscles, over the range of 10%-75% maximum voluntary contraction. EMG amplitude estimates were computed with eight different EMG processor schemes using smoothing length durations spanning 2.45-500 ms. An SNR was computed from each amplitude estimate (deviations about the mean value of the estimate were considered as noise). Over these window lengths, average ± standard deviation SNR's ranged from 1.4±0.28 to 16.2±5.4 for unwhitened single-channel EMG processing and from 3.2±0.7 to 37.3±14.2 for whitened, multiple-channel EMG processing (results pooled across contraction level). It was found that SNR increased with window length in a square root fashion. The shape of this relationship was consistent with classic theoretical predictions, however none of the processors achieved the absolute performance level predicted by the theory. These results are useful in selecting the length of the smoothing window in traditional surface EMG studies. In addition, this study should contribute to the development of EMG processors which dynamically tune the smoothing window length when the EMG amplitude is time varying     

Adaptive whitening of the electromyogram to improve amplitude estimation

Previous research showed that whitening the surface electromyogram (EMG) can improve EMG amplitude estimation (where EMG amplitude is defined as the time-varying standard deviation of the EMG). However, conventional whitening via a linear filter seems to fail at low EMG amplitude levels, perhaps due to additive background noise in the measured EMG. This paper describes an adaptive whitening technique that overcomes this problem by cascading a nonadaptive whitening filter, an adaptive Wiener filter, and an adaptive gain correction. These stages can be calibrated from two, five second duration, constant-angle, constant-force contractions, one at a reference level [e.g., 50% maximum voluntary contraction (MVC)] and one at 0% MVC. In experimental studies, subjects used real-time EMG amplitude estimates to track a uniform-density, band-limited random target. With a 0.25-Hz bandwidth target, either adaptive whitening or multiple-channel processing reduced the tracking error roughly half-way to the error achieved using the dynamometer signal as the feedback. At the 1.00-Hz bandwidth, all of the EMG processors had errors equivalent to that of the dynamometer signal, reflecting that errors in this task were dominated by subjects' inability to track targets at this bandwidth. Increases in the additive noise level, smoothing window length, and tracking bandwidth diminish the advantages of whitening.     

Nonlinear Spatial Filtering of Multichannel Surface Electromyogram Signals During Low Force Contractions

This study introduces the application of nonlinear spatial filters to help identify single motor unit discharge from multiple channel surface electromyogram (EMG) signals during low force contractions. The nonlinear spatial filters simultaneously take into account the instantaneous amplitude and frequency information of a signal. This property was used to enhance motor unit action potentials (MUAPs) in the surface EMG record. The advantages of nonlinear spatial filtering for surface MUAP enhancement were investigated using both simulation and experimental approaches. The simulation results indicate that when compared with various linear spatial filters, nonlinear spatial filtering achieved higher SNR and higher kurtosis of the surface EMG distribution. Over a broad range of SNR and kurtosis levels for the input signal, nonlinear spatial filters achieved at least 32 times greater SNR and 11% higher kurtosis for correlated noise, and at least 15 times greater SNR and 1.7 times higher kurtosis for independent noise, across electrode array channels. The improvements offered by nonlinear spatial filters were further documented by applying them to experimental surface EMG array recordings. Compared with linear spatial filters, nonlinear spatial filters achieved at least nine times greater SNR and 25% higher kurtosis. It follows that nonlinear spatial filters represent a potentially useful supplement to linear spatial filters for detection of motor unit activity in surface EMG at low force contractions.     

Influence of joint angle on the calibration and performance of EMGamplitude estimators

Multiple-channel electromyogram (EMG) amplitude estimators incorporating temporal whitening filters and/or spatial uncorrelation filters contain a characterization of the EMG waveform (specifically, auto- and cross-correlation information) which may vary with joint angle. This paper reports on an experimental study which investigated the influence of joint angle on these EMG amplitude estimators. It was found that little or no relative improvement in estimator performance resulted from altering either temporal whitening or spatial uncorrelation filters as a function of joint angle. Also, the absolute performance level of these estimators did not vary with joint angle     

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.     

Probability density of the surface electromyogram and its relationto amplitude detectors

When the surface electromyogram (EMG) generated from constant-force, constant-angle, nonfatiguing contractions is modeled as a random process, its density is typically assumed to be Gaussian. This assumption leads to root-mean-square (RMS) processing as the maximum likelihood estimator of the EMG amplitude (where EMG amplitude is defined as the standard deviation of the random process). Contrary to this theoretical formulation, experimental work has found the signal-to-noise-ratio [(SNR), defined as the mean of the amplitude estimate divided by its standard deviation] using mean-absolute-value (MAV) processing to be superior to RMS. This paper reviews RMS processing with the Gaussian model and then derives the expected (inferior) SNR performance of MAV processing with the Gaussian model. Next, a new model for the surface EMG signal, using a Laplacian density, is presented. It is shown that the MAV processor is the maximum likelihood estimator of the EMG amplitude for the Laplacian model. SNR performance based on a Laplacian model is predicted to be inferior to that of the Gaussian model by approximately 32%. Thus, minor variations in the probability distribution of the EMG may result in large decrements in SNR performance. Lastly, experimental data from constant-force, constant-angle, nonfatiguing contractions were examined. The experimentally observed densities fell in between the theoretic Gaussian and Laplacian densities. On average, the Gaussian density best fit the experimental data, although results varied with subject. For amplitude estimation, MAV processing had a slightly higher SNR than RMS processing.     

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.     

Adaptive whitening in electromyogram amplitude estimation for epoch-based applications

Epoch-based electromyogram (EMG) amplitude estimates have not incorporated signal whitening, even though whitening has demonstrated significant improvements for stream-based estimates. This paper presents new epoch-based algorithms, for both single- and multiple-channel EMG, which include a whitening stage. The best multiple-channel whitening processor provided a 21.4%-22.5% improvement over single-channel unwhitened estimation in an EMG-to-torque application.     

Relating agonist-antagonist electromyograms to joint torque duringisometric, quasi-isotonic, nonfatiguing contractions

Describes an experimental study which relates simultaneous elbow flexor-extensor electromyogram (EMG) amplitude to joint torque. Investigation was limited to the case of isometric, quasi-isotonic (slowly force-varying), nonfatiguing contractions. For each of the flexor and extensor muscle groups, the model relationship between muscle group torque contribution and EMG amplitude was constrained to be a sum of basis functions which had a linear dependence on a set of fit parameters. With these constraints, the problem of identifying the EMG-to-torque relationship was reduced to a linear least squares problem. Surface EMGs from elbow flexors and extensors, and joint torque were simultaneously recorded for nonfatiguing, quasi-isotonic, isometric contractions spanning 0-50% maximum voluntary contraction. Single-/multiple-channel unwhitened/whitened/adaptively-whitened EMG amplitude processors were used to identify an EMG-to-torque relation, and then estimate joint torque based on this relation. Each unwhitened multiple-channel EMG-to-torque estimator had a standard error (SE) approximately 70% of its respective single-channel estimator. The adaptively whitened multiple-channel joint torque estimator had an SE approximately 90% of the unwhitened multiple-channel estimator, providing an estimation error ≈3% of the combined flexion/extension torque range. The experimental studies demonstrated that higher fidelity EMG amplitude processing led to improved joint torque estimation     

Analysis and simulation of changes in EMG amplitude during high-level fatiguing contractions

Changes in surface electromyographic (EMG) amplitude during sustained, fatiguing contractions are commonly attributed to variations in muscle fiber conduction velocity (MFCV), motor unit firing rates, transmembrane action potentials and the synchronization or recruitment of motor units. However, the relative contribution of each factor remains unclear. Analytical relationships relating changes in MFCV and mean motor unit firing rates to the root mean square (RMS) and average rectified (AR) value of the surface EMG signal are derived. The relationships are then confirmed using model simulation. The simulations and analysis illustrate the different behaviors of the surface EMG RMS and AR value with changing MFCV and firing rate, as the level of motor unit superposition varies. Levels of firing rate modulation and short-term synchronization that, combined with variations in MFCV, could cause changes in EMG amplitude similar to those observed during sustained isometric contraction of the brachioradialis at 80% of maximum voluntary contraction were estimated. While it is not possible to draw conclusions about changes in neural control without further information about the underlying motor unit activation patterns, the examples presented illustrate how a combined analytical and simulation approach may provide insight into the manner in which different factors affect EMG amplitude during sustained isometric contractions.     

Comments on "Is accurate recording of the ECG surface Laplacian feasible?"

The authors comment on the work of D.B. Geselowitz and J. Ferrara (ibid., vol. 47, p. 1138, 2000) who concluded that the signal-to-noise ratio (SNR) of the ECG surface Laplacian (SL) is very low. Averaged over all recording channels on the anterolateral chest and averaged over all six subjects, when the noise level is estimated during the window of the P-R segment, our experimental results indicated that the potential ECG, regular Laplacian ECG, and diagonal Laplacian ECG have SNR of 105.70±22.49, 33.73±9.86, and 46.83±12.84, respectively. This is in contrast with the magnitude of SNR reported in the text of the above paper, which was too low becaase of their choice of a residue window containing noise as well as signal. Geselowitz and Ferrara reply to these comments     

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.92 +/- 0.07 with wavelet optimization, 0.74 +/- 0.09 with the wavelet leading to the poorest performance, 0.85 +/- 0.07 with Wigner-Ville distribution, 0.86 +/- 0.07 with Choi-Williams distribution, and 0.73 +/- 0.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.2 +/- 10.0% before BSS and was reduced to 15.2 +/- 6.3% with wavelet optimization, 30.1 +/- 15.0% with the worst wavelet, 28.3 +/- 12.3% with Wigner-Ville distribution, 26.2 +/- 12.0% with Choi-Williams distribution, and 35.1 +/- 15.5% with second-order statistics. In conclusion, the proposed approach resulted in better performance than previous methods for the separation of nonstationary myoelectric signals.     

Modeling of surface myoelectric signals. II. Model-based signalinterpretation

Experimental electromyogram (EMG) data from the human biceps brachii were simulated using the model described in [10] of this work. A multichannel linear electrode array, spanning the length of the biceps, was used to detect monopolar and bipolar signals, from which double differential signals were computed, during either voluntary or electrically elicited isometric contractions. For relatively low-level voluntary contractions (10%-30% of maximum force) individual firings of three to four-different motor units were identified and their waveforms were closely approximated by the model. Motor unit parameters such as depth, size, fiber orientation and length, location of innervation and tendonous zones, propagation velocity, and source width were estimated using the model. Two applications of the model are described. The first analyzes the effects of electrode rotation with respect to the muscle fiber direction and shows the possibility of conduction velocity (CV) over- and under-estimation. The second focuses on the myoelectric manifestations of fatigue during a sustained electrically elicited contraction and the interrelationship between muscle fiber CV, spectral and amplitude variables, and the length of the depolarization zone. It is concluded that a) surface EMG detection using an electrode array, when combined with a model of signal propagation, provides a useful method for understanding the physiological and anatomical determinants of EMG waveform characteristics and b) the model provides a way for the interpretation of fatigue plots.     

Entropy-based optimization of wavelet spatial filters

A new class of spatial filters for surface electromyographic (EMG) signal detection is proposed. These filters are based on the 2-D spatial wavelet decomposition of the surface EMG recorded with a grid of electrodes and inverse transformation after zeroing a subset of the transformation coefficients. The filter transfer function depends on the selected mother wavelet in the two spatial directions. Wavelet parameterization is proposed with the aim of signal-based optimization of the transfer function of the spatial filter. The optimization criterion was the minimization of the entropy of the time samples of the output signal. The optimized spatial filter is linear and space invariant. In simulated and experimental recordings, the optimized wavelet filter showed increased selectivity with respect to previously proposed filters. For example, in simulation, the ratio between the peak-to-peak amplitude of action potentials generated by motor units 20 degrees apart in the transversal direction was 8.58% (with monopolar recording), 2.47% (double differential), 2.59% (normal double differential), and 0.47% (optimized wavelet filter). In experimental recordings, the duration of the detected action potentials decreased from (mean +/- SD) 6.9 +/- 0.3 ms (monopolar recording), to 4.5 +/- 0.2 ms (normal double differential), 3.7 +/- 0.2 (double differential), and 3.0 +/- 0.1 ms (optimized wavelet filter). In conclusion, the new class of spatial filters with the proposed signal-based optimization of the transfer function allows better discrimination of individual motor unit activities in surface EMG recordings than it was previously possible.     

Maximum-likelihood synchronization of a single user forcode-division multiple-access communication systems

Code-division multiple access (CDMA) has emerged as an access protocol well-suited for voice and data transmission. One significant limitation of the conventional CDMA system is the near-far problem where strong signals interfere with the detection of a weak signal. Multiuser detectors assume knowledge of all of the modulation waveforms and channel parameters, and exploit this information to eliminate multiple-access interference (MAI) and to achieve near-far resistance. A major problem in practical application of multiuser detection is the estimation of the signal and channel parameters in a near-far limited system. We consider maximum-likelihood estimation of users delay, amplitude, and phase in a CDMA communication system. We present an approach for decomposing this multiuser estimation problem into a series of single-user problems. In this method the interfering users are treated as colored non-Gaussian noise. The observation vectors are preprocessed to be able to apply a Gaussian model for the MAI. The maximum-likelihood estimate (MLE) of each user's parameters based on the processed observation vectors becomes tractable. The estimator includes a whitening filter derived from the sample covariance matrix which is used to suppress the MAI, thus yielding a near-far resistant estimator     

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     

The identification of time series models of lower extremity EMG forthe control of prostheses using Box-Jenkins criteria

The Box-Jenkins method of model identification and diagnostic checking was used on the EMG of the posterio-lateral thigh while contraction level, electrode position, and limb function were varied. Models remained exclusively autoregressive (AR) of order typically less than five for isometric contractions ranging from 25 to 50% of the maximum voluntary knee flexions and extensions, and for recording sites located within 30 mm of the original electrode position. Diagnostics performed on the residuals from AR(4) models indicate that a fourth-order model is adequate for myoelectric control applications under many circumstances     

Motor Unit Firing Rate During Static Contraction Indicated by the Surface EMG Power Spectrum

In this study, power spectral density functions (PSDF's) were computed of interference EMG of various facial and jaw-elevator muscles during nonfatiguing submaximal static contractions, recorded with surface electrodes. A distinct peak was found in the PSDF's in the frequency region below 40 Hz. It was shown that the peak was due to genuine EMG activity and that it could not be considered as an artifact, which was caused by electrode displacements during contraction. An increase of contraction strength resulted in a shift of the peak to higher frequencies and a decrease of peak amplitude relative to the power spectral estimates above 40 Hz, which were shown to be determined by the shape of the motor unit (MU) action potentials. In accordance with mathematical models of the EMG PSDF, it was demonstrated that the peak indicates the dominant firing rate of the sampled MU's. Our results suggest that this can be defined as the firing rate of the first recruited low-threshold MU's, which may be expected to dominate the interference EMG signal because of their preponderance in number. The data further suggest that the peak can be more readily observed in PSDF's of facial and jaw-elevator muscles than in PSDF's of limb muscles. This might be related to differences in MU firing statistics.     

Efficient Complex Sphere Decoding for MC-CDMA Systems

Maximum likelihood (ML) joint detection of multi-carrier code division multiple access (MC-CDMA) systems can be efficiently implemented with a sphere decoding (SD) algorithm. In this paper, we examine the application of complex instead of real SD to detect MC-CDMA, which solves many problems in a more elegant manner and extends SD adaptability to any constellation. We first propose a new complex SD algorithm whose efficiency is based on not requiring an estimate of the initial search radius but selecting the Babai point as the initial sphere radius instead; also, efficient strategies regarding sorting the list of possible lattice points are applied. Indeed, complex SD allows complex matrix operations which are faster than real counterparts in double dimension. Next, a novel lattice representation for the MC-CDMA system is introduced, which allows optimum multiuser detection directly from the received signal. This avoids noise whitening operation, and also despreading and equalization procedures are not required further at the receiver side