Robust Muscle Activity Onset Detection Using an Unsupervised Electromyogram Learning Framework

Accurate muscle activity onset detection is an essential prerequisite for many applications of surface electromyogram (EMG). This study presents an unsupervised EMG learning framework based on a sequential Gaussian mixture model (GMM) to detect muscle activity onsets. The distribution of the logarithmic power of EMG signal was characterized by a two-component GMM in each frequency band, in which the two components respectively correspond to the posterior distribution of EMG burst and non-burst logarithmic powers. The parameter set of the GMM was sequentially estimated based on maximum likelihood, subject to constraints derived from the relationship between EMG burst and non-burst distributions. An optimal threshold for EMG burst/non-burst classification was determined using the GMM at each frequency band, and the final decision was obtained by a voting procedure. The proposed novel framework was applied to simulated and experimental surface EMG signals for muscle activity onset detection. Compared with conventional approaches, it demonstrated robust performance for low and changing signal to noise ratios in a dynamic environment. The framework is applicable for real-time implementation, and does not require the assumption of non EMG burst in the initial stage. Such features facilitate its practical application.     

Sample entropy analysis of surface EMG for improved muscle activity onset detection against spurious background spikes

Voluntary surface electromyogram (EMG) signal is sometimes contaminated by spurious background spikes of both physiological and extrinsic or accidental origins. A novel method of muscle activity onset detection against such spurious spikes was proposed in this study based primarily on the sample entropy (SampEn) analysis of the surface EMG. The method takes advantage of the nonlinear properties of the SampEn analysis to distinguish voluntary surface EMG signals from spurious background spikes in the complexity domain. To facilitate muscle activity onset detection, the SampEn analysis of surface EMG was first performed to highlight voluntary EMG activity while suppressing spurious background spikes. Then, a SampEn threshold was optimized for muscle activity onset detection. The performance of the proposed method was examined using both semi-synthetic and experimental surface EMG signals. The SampEn based methods effectively reduced the detection error induced by spurious background spikes and achieved improved performance over the methods relying on conventional amplitude thresholding or its extended version in the Teager Kaiser Energy domain.     

Use of surface EMG for studying motor unit recruitment during isometric linear force ramp.

The aim of this work was to demonstrate the rank order of motor unit (MU) recruitment by surface EMG based on a Laplacian detection technique and to document the MU features at their recruitment threshold. Surface EMG signals were recorded on the biceps brachii of 10 healthy subjects during linear force ramps. When achievable, the signals were decomposed into MU action potential (MUAP) trains. MU inter-pulse interval (IPI), conduction velocity (MUCV) and amplitude were estimated on the first 12 MUAPs of each detectable train in order to characterize the MU features at their firing onset. A strong correlation was found between MU recruitment threshold and IPI, MUCV, and amplitude, showing that the size principle can be demonstrated by a fully non-invasive EMG technique. However, signal decomposition was not possible on seven subjects due to the effects of the volume conductor when the skinfold thickness was too large. When requirements for an optimal detection of MUAP trains are met, surface EMG may be used to improve our understanding of MU activity.     

Tonic-to-phasic shift of lumbo-pelvic muscle activity during 8 weeks of bed rest and 6-months follow up.

Prior motor control studies in unloading have shown a tonic-to-phasic shift in muscle activation, particularly in the short extensors. Tonic muscle activity is considered critical for normal musculoskeletal function. The shift from tonic-to-phasic muscle activity has not been systematically studied in humans in unloading nor at the lumbo-pelvic (LP) region. Ten healthy young male subjects underwent 8 wk of bed rest with 6-mo follow up as part of the "Berlin Bed-Rest Study." A repetitive knee movement model performed in the prone position is used to stimulate tonic holding LP muscle activity, as measured by superficial EMG. Tonic and phasic activation patterns were quantified by relative height of burst vs. baseline electromyographic linear-envelope signal components. Statistical analysis shows a shift toward greater phasic activity during bed rest and follow up (P < 0.001) with a significant interaction across muscles (P < 0.001) specifically affecting the short lumbar extensors. These changes appear unrelated to skill acquisition over time (P all > or = 0.196). This change of a shift from tonic LP muscle activation to phasic is in line with prior research on the effects of reduced weight bearing on motor control.     

Long-term electromyographic activity in upper trapezius and low back muscles of women with moderate physical activity.

The habitual activity patterns of trapezius and postural back muscles (multifidus, iliocostalis, longissimus) of 23 female subjects with moderate physical activity were studied. Bilateral surface electromyographic (sEMG) recordings from start of work until bedtime were analyzed. The activity level was calibrated as percentage of root mean square-detected muscle activity at maximal voluntary contraction (EMG(max)). Sixty-six previous trapezius recordings of women with moderate physical activity were included in some analyses to pursue the full range of variation in trapezius activity. Twenty-six of these were recorded twice, separated by 16-28 mo. Median activity level and duration of periods with sEMG activity of <0.5% EMG(max) ("rest time"; only trapezius) and exceeding 2 ("burst time"), 10, 30, and 50% EMG(max) was determined. The trapezius median activity level ranged from 0.6 to 8.8% EMG(max), burst time from 9 to 84%, and rest time from 2 to 84%. The activity patterns of the back muscles showed similar large interindividual variation. Repeated trapezius recordings of the same subject showed high consistency; intraclass correlation coefficients ranged from 0.62 to 0.79 for different sEMG variables. Periods with high sEMG amplitude were of short duration; 7% of the trapezius recordings did not present time intervals (0.2-s duration) above 50% EMG(max). The activity patterns of the postural muscles, despite large interindividual variability, were distinctly different from activity patterns of upper and lower limb muscles reported by others (e.g., mean burst time 40-50 vs. 10-20%). We conclude that postural trunk muscles show idiosyncratic activity patterns with large interindividual variation. High-threshold motor units are activated to a very minor extent.     

Modeling nonlinear errors in surface electromyography due to baseline noise: A new methodology

The surface electromyographic (EMG) signal is often contaminated by some degree of baseline noise. It is customary for scientists to subtract baseline noise from the measured EMG signal prior to further analyses based on the assumption that baseline noise adds linearly to the observed EMG signal. The stochastic nature of both the baseline and EMG signal, however, may invalidate this assumption. Alternately, “true” EMG signals may be either minimally or nonlinearly affected by baseline noise. This information is particularly relevant at low contraction intensities when signal-to-noise ratios (SNR) may be lowest. Thus, the purpose of this simulation study was to investigate the influence of varying levels of baseline noise (approximately 2–40% maximum EMG amplitude) on mean EMG burst amplitude and to assess the best means to account for signal noise. The simulations indicated baseline noise had minimal effects on mean EMG activity for maximum contractions, but increased nonlinearly with increasing noise levels and decreasing signal amplitudes. Thus, the simple baseline noise subtraction resulted in substantial error when estimating mean activity during low intensity EMG bursts. Conversely, correcting EMG signal as a nonlinear function of both baseline and measured signal amplitude provided highly accurate estimates of EMG amplitude. This novel nonlinear error modeling approach has potential implications for EMG signal processing, particularly when assessing co-activation of antagonist muscles or small amplitude contractions where the SNR can be low.     

Influence of amplitude cancellation on the accuracy of determining the onset of muscle activity from the surface electromyogram

The purpose of the study was to quantify the influence of amplitude cancellation on the accuracy of detecting the onset of muscle activity based on an analysis of simulated surface electromyographic (EMG) signals. EMG activity of a generic lower limb muscle was simulated during the stance phase of human gait. Surface EMG signals were generated with and without amplitude cancellation by summing simulated motor unit potentials either before (cancellation EMG) or after (no-cancellation EMG) the potentials had been rectified. The two sets of EMG signals were compared at forces of 30% and 80% of maximum voluntary contraction (MVC) and with various low-pass filter cut-off frequencies. Onset time was determined both visually and by an algorithm that identified when the mean amplitude of the signal within a sliding window exceeded a specified standard deviation (SD) above the baseline mean. Onset error was greater for the no-cancellation conditions when determined automatically and by visual inspection. However, the differences in onset error between the two cancellation conditions appear to be clinically insignificant. Therefore, amplitude cancellation does not appear to limit the ability to detect the onset of muscle activity from the surface EMG.     

Simulations of the alpha motoneuron pool electromyogram reflex at different preactivation levels in man

The alpha motoneuron pool and the surface electromyogram (EMG) of the human soleus muscle are modelled, respectively, by an alpha motoneuron pool model generating the firing patterns in the motor units of e muscle and by a muscle model using these discharge patterns to simulate the surface EMG. In the alpha motoneuron pool model, we use a population of motoneurons in which cellular properties like cell size and membrane conductance are distributed according to experimentally observed data. By calculating the contribution from each motor unit, the muscle model predicts the EMG. Wave forms of the motor unit action potentials in the surface EMG are obtained from experimental data. Using the model, we are able to give a quantitative prediction of the motoneuron pool activity and the reflex EMG output at different preactivation levels. The simulated data are consistent with experimentally obtained results in healthy humans. During static isometric muscle preactivations, the simulations show that the reflex strength is highly dependent on the intrinsic threshold properties of the alpha motoneuron pool. Received: 27 April 1993/Accepted in revised form: 8 September 1993     

Digital filter design for peak detection of surface EMG.

A simple Low-Pass Differential (LPD) filter is often used for the Motor Unit Action Potential (MUAP) peak detection in needle EMG decomposition. Decomposition of surface EMG is much more difficult since surface EMG is more mixed than needle EMG. A successful peak detection of MUAPs is a first important step for EMG decomposition. We found that a simple LPD filter is not suitable for surface MUAP detection; instead the Weighted Low-Pass Differential (WLPD) filters are proposed. The filter performance is analyzed based on different window selection and varying MUAPs from simulated and recorded surface EMG. The sinusoidal WLPD filter is found to have better SNR improving factors, and to be more robust under the varying conditions.     

The influence of increasing steady-state walking speed on muscle activity in below-knee amputees

The goal of this study was to identify changes in muscle activity in below-knee amputees in response to increasing steady-state walking speeds. Bilateral electromyographic (EMG) data were collected from 14 amputee and 10 non-amputee subjects during four overground walking speeds from eight intact leg and five residual leg muscles. Using integrated EMG measures, we tested three hypotheses for each muscle: (1) there would be no difference in muscle activity between the residual and intact legs, (2) there would be no difference in muscle activity between the intact leg and non-amputee legs, and (3) muscle activity in the residual and intact legs would increase with speed. Most amputee EMG patterns were similar between legs and increased in magnitude with speed. Differences occurred in the residual leg biceps femoris long head, vastus lateralis and rectus femoris, which increased in magnitude during braking compared to the intact leg. These adaptations were consistent with the need for additional body support and forward propulsion in the absence of the plantar flexors. With the exception of the intact leg gluteus medius, all intact leg muscles exhibited similar EMG patterns compared to the control leg. Finally, the residual, intact and control leg EMG all had a significant speed effect that increased with speed with the exception of the gluteus medius.     

Time-frequency analysis of surface electromyographic signals during fatiguing isokinetic muscle actions

The purpose of this study was to use a wavelet analysis designed specifically for electromyography (EMG) signals in combination with a trend plot to examine changes in EMG intensity patterns during maximal, fatiguing isokinetic muscle actions. Eleven men (mean ± SD age = 22.4 ± 1.1 years) and 7 women (mean ± SD age = 22.7 ± 2.1 years) performed 50 consecutive maximal concentric isokinetic muscle actions of the dominant leg extensors at a velocity of 180°·s(-1). During each muscle action, a bipolar surface EMG signal was detected from the vastus lateralis. All signals were then processed with a wavelet analysis designed specifically for EMG signals, which resulted in EMG intensity patterns. The patterns for each subject were then analyzed with a trend plot, which provided information regarding the changes that occurred because of fatigue. The results indicated that for all the 18 subjects, the EMG intensity patterns moved in a predictable manner in pattern space, but the changes to the patterns were different for each subject. These findings reflect the complex changes that occur in the EMG signal during fatigue. These changes cannot be characterized fully with a single amplitude and center frequency parameter and can be useful for athletes and coaches who need to track the fatigue status of individual muscles.     

Comparing the lactate and EMG thresholds of recreational cyclists during incremental pedaling exercise

The purpose of this study was to determine the validity of using the electromyography (EMG) signal as a noninvasive method of estimating the lactate threshold (LT) power output in recreational cyclists. Using an electromagnetic bicycle ergometer and constant pedaling cadence of 80 rpm, 24 recreational cyclists performed an incremental exercise protocol that consisted of stepwise increases in power output of 25 W every 3 min until exhaustion. The EMG signal was recorded from the right vastus lateralis (VL) and right rectus femoris (RF) throughout the test. Blood samples were taken from the fingertip every 3 min. The LT was determined by examining the relation between the lactate concentration and the power output using a log-log transformation model. The root mean square (RMS) value from the EMG signal was calculated for every 1-second non-superimposing window. Sets of pairs of straight regression lines were plotted and the corresponding determination coefficients (R(2)) were calculated. The intersection point of the pair of lines with the highest R(2) product was chosen to represent the EMG threshold (EMGT). The results showed that the correlation coefficients (r) between EMGT and LT were significant (p < 0.01) and high for the VL (r = 0.826) and RF (r = 0.872). The RF and VL muscles showed similar behavior during the maximal incremental test and the EMGT and LT power output were equivalent for both muscles. The validity of using EMG to estimate the LT power output in recreational cyclists was confirmed.     

Relation between the defensive behavior of altricial nestlings and the parameters of an alarm signal. Neuroethologic approach

The acoustic organization of the defensive behaviour was studied in Pied Flycatcher nestlings. Electrodes were chronically implanted into the field L and the neck muscles for auditory EPs and EMG recordings. A well developed defensive reaction which under natural conditions is elicited by species-specific alarm call, could be artificially provoked by monotonal 4.5-6 kHz bursts with repetition frequency of 4-0,8/s. The tone frequency was shown to be the critical acoustic parameter permitting the altricial nestlings to identify alarm vocalization and differentiate it from alimentary signals. The role of the burst repetition frequency is restricted to maintaining the tonic defensive activation in the course of the alarm signal presentation. Certain correlation was found between the time of regeneration of excitability in auditory structures and the burst repetition frequency in the species specific alarm signal.     

Effect of age and weight on upper airway function in a mouse model

Defects in pharyngeal mechanical and neuromuscular control are required for the development of obstructive sleep apnea. Obesity and age are known sleep apnea risk factors, leading us to hypothesize that specific defects in upper airway neuromechanical control are associated with weight and age in a mouse model. In anesthetized, spontaneously breathing young and old wild-type C57BL/6J mice, genioglossus electromyographic activity (EMG(GG)) was monitored and upper airway pressure-flow dynamics were characterized during ramp decreases in nasal pressure (Pn, cmH?O). Specific body weights were targeted by controlling caloric intake. The passive critical pressure (Pcrit) was derived from pressure-flow relationships during expiration. The Pn threshold at which inspiratory flow limitation (IFL) developed and tonic and phasic EMG(GG) activity during IFL were quantified to assess the phasic modulation of pharyngeal patency. The passive Pcrit increased progressively with increasing body weight and increased more in the old than young mice. Tonic EMG(GG) decreased and phasic EMG(GG) increased significantly with obesity. During ramp decreases in Pn, IFL developed at a higher (less negative) Pn threshold in the obese than lean mice, although the frequency of IFL decreased with age and weight. The findings suggest that weight imposes mechanical loads on the upper airway that are greater in the old than young mice. The susceptibility to upper airway obstruction increases with age and weight as tonic neuromuscular activity falls. IFL can elicit phasic responses in normal mice that mitigate or eliminate the obstruction altogether.     

Experimental Evidence of the Tonic Vibration Reflex during Whole-Body Vibration of the Loaded and Unloaded Leg

Increased muscle activation during whole-body vibration (WBV) is mainly ascribed to a complex spinal and supraspinal neurophysiological mechanism termed the tonic vibration reflex (TVR). However, TVR has not been experimentally demonstrated during low-frequency WBV, therefore this investigation aimed to determine the expression of TVR during WBV.  Whilst seated, eight healthy males were exposed to either vertical WBV applied to the leg via the plantar-surface of the foot, or Achilles tendon vibration (ATV) at 25Hz and 50Hzfor 70s. Ankle plantar-flexion force, tri-axial accelerations at the shank and vibration source, and surface EMG activity of m. soleus (SOL) and m. tibialis anterior (TA) were recorded from the unloaded and passively loaded leg to simulate body mass supported during standing.  Plantar flexion force was similarly augmented by WBV and ATV and increased over time in a load- and frequency dependent fashion. SOL and TA EMG amplitudes increased over time in all conditions independently of vibration mode. 50Hz WBV and ATV resulted in greater muscle activation than 25Hz in SOL when the shank was loaded and in TA when the shank was unloaded despite the greater transmission of vertical acceleration from source to shank with 25Hz and WBV, especially during loading. Low-amplitude WBV of the unloaded and passively loaded leg produced slow tonic muscle contraction and plantar-flexion force increase of similar magnitudes to those induced by Achilles tendon vibration at the same frequencies. This study provides the first experimental evidence supporting the TVR as a plausible mechanism underlying the neuromuscular response to whole-body vibration.     

Development of coordinated movement in chicks: II. Temporal analysis of hindlimb muscle synergies at embryonic day 10 in embryos with spinal gap transections

Spinal neural circuits can recruit muscles to produce organized patterns of activity early in embryonic development. In a previous study, using multichannel electromyographic (EMG) recordings, we characterized burst parameters for these patterns in the legs of chick embryos during spontaneous motility in ovo at embryonic days (E) 9 and E10 (Bradley and Bekoff, 1990). Results of the study suggested both neural and biomechanical factors play an important role in the development of coordinated limb movements. In this study, to explore the contribution of descending neural inputs to the control of leg movements during motility, we applied similar methods to characterize motor patterns produced by the spinal cord in the absence of descending inputs. Thoracic spinal gap transections were performed at E2 and EMG patterns were recorded at E10. Several EMG features for chronic spinal embryos were similar to those for normal embryos and demonstrate that lumbar spinal circuits can be correctly assembled to control limb movements in the absence of connectivity with more rostral neural structures during early differentiation processes. However, certain aspects of the EMG patterns in chronic spinal embryos were different from patterns in normal embryos and provide support for conclusions drawn earlier by Oppenheim (1975). Specifically, our data support the view that propriospinal and/or supraspinal inputs function to regulate the timing of cyclic limb movements controlled by spinal neural circuits. Finally, we consider the possible long-term effects of chronic spinal gap transections as compared to acute spinal transections on the development of motility. © 1992 John Wiley & Sons, Inc.     

Development of coordinated movement in chicks: II. Temporal analysis of hindlimb muscle synergies at embryonic day 10 in embryos with spinal gap transections.

Spinal neural circuits can recruit muscles to produce organized patterns of activity early in embryonic development. In a previous study, using multichannel electromyographic (EMG) recordings, we characterized burst parameters for these patterns in the legs of chick embryos during spontaneous motility in ovo at embryonic days (E) 9 and E10 (Bradley and Bekoff, 1990). Results of the study suggested both neural and biomechanical factors play an important role in the development of coordinated limb movements. In this study, to explore the contribution of descending neural inputs to the control of leg movements during motility, we applied similar methods to characterize motor patterns produced by the spinal cord in the absence of descending inputs. Thoracic spinal gap transections were performed at E2 and EMG patterns were recorded at E10. Several EMG features for chronic spinal embryos were similar to those for normal embryos and demonstrate that lumbar spinal circuits can be correctly assembled to control limb movements in the absence of connectivity with more rostral neural structures during early differentiation processes. However, certain aspects of the EMG patterns in chronic spinal embryos were different from patterns in normal embryos and provide support for conclusions drawn earlier by Oppenheim (1975). Specifically, our data support the view that propriospinal and/or supraspinal inputs function to regulate the timing of cyclic limb movements controlled by spinal neural circuits. Finally, we consider the possible long-term effects of chronic spinal gap transections as compared to acute spinal transections on the development of motility.     

Cross-education after one session of unilateral surface electrical stimulation of the rectus femoris

Thirty-six adult men were randomly assigned to a remote stimulation group (RS; n = 18) or control group (CTL; n = 18). The RS group unilaterally performed a 10-minute surface electrical stimulation program (frequency 100 Hz, impulse 300 micros, 10 seconds on/10 seconds off) on the rectus femoris of the non-dominant leg. The subjects of the CTL group relaxed for 10 minutes without performing any training. Immediately before and after the surface electrical stimulation program, the isometric strength and the electromyographic (EMG) and mechanomyographic (MMG) response of the dominant leg was measured for all subjects. The dominant leg of the RS group showed a significant increase in the isometric force (5.11%; P < 0.001) and EMG activity of the agonist muscle (4.67%; P < 0.05), whereas a decrease in EMG activity of the antagonist muscles was observed (-10.27%; P < 0.05). The MMG activity did not show any alteration. No significant changes were observed for the CTL group. These results indicate that one unilateral surface electrical stimulation session on the rectus femoris improves the efficiency of the inactive leg. At a practical level, the results open a new way to rehabilitate muscle-skeletal injuries, especially weak members that cannot do any physical work. In this case, the muscle strength (and physical efficiency) can be improved by passive electrostimulation training on the healthy member.     

Electrophysiological and pharmacological analysis of L-dopa-induced dyskinesia and tardive dyskinesia (author's transl)

Electrophysiological and pharmacological analysis of L-Dopa-induced dyskinesia and tardive dyskinesia (L.DD) due to neuroleptics was performed on 12 patients with Parkinson's disease and on 12 others with psychotic diseases. This analysis included the examination of spinal reflexes, monosynaptic H reflex, polysynaptic cutaneous reflex of the lower limb, muscular responses to passive movement [stretch reflex and shortening reaction (SR)] and the study of the motor response to a dopaminergic stimulus (I.V. injection of Piribedil (PBD), a dopamine agonist). There was no difference in EMG activity between L.DD and TD. Three EMG patterns can be distinguished: anarchic discharge pattern (ADA), tonic grouping discharge pattern (AST) and rhythmic burst pattern (ABR). PBD effects indicate a possible relationship between the EMG patterns and the sensitivity level of the motor dopamine receptors. During L-Dopa dyskinesia and tardive dyskinesia, the same changes in spinal reflexes were observed. Muscle tone tested by muscular responses to passive movement (shortening and myotatic reaction) was normal. Monosynaptic excitability explored by H/M ratio was within the normal range. In contrast, the polysynaptic nociceptive reflex was increased in every case. In Parkinsonian patients with L-Dopa dyskinesia, this pattern of the spinal reflexes was significantly different in comparison to the rigid phase. Intravenous infusion of PBD suppressed tremor and provoked the occurrence of dyskinetic activity in Parkinsonian patients with L-Dopa dyskinesia during the rigid phase. During the dyskinetic phase, as in tardive dyskinesia, PBD increases these phenomena and changes EMG activity in rhythmic pattern. It is suggested that L-Dopa dyskinesia and tardive dyskinesia can be determined by testing EMG activity, spinal reflexes and dopaminergic reactivity. There is evidence to suggest that the various types of involuntary abnormal movement represent a single entity, and that dopamine receptor supersensitivity may be involved.     

Reaction time and electromyographic activity during a sprint start

Eight male sprinters were filmed running three maximal starts over 3 m on a long force platform. The subjects were divided into two groups (n = 4) according to the leg on which the electromyograph (EMG) electrodes were fixed. When in the set position one group had electrodes on the front leg (FLG) and the other group on the rear leg (RLG). The EMG activities of the gastrocnemius caput laterale muscle (GA), vastus lateralis muscle (VL), biceps femoris caput longum muscle (BF), rectus femoris muscle (RF) and gluteus maximus muscle (GM) were recorded telemetrically using surface electrodes. Total reaction time (TRT) was defined as the time from the gun signal until a horizontal force was produced with a value 10% above the base line. Pre-motor time was defined as the time from the gun signal until the onset of EMG activity and motor time (MT) as the time between the onset of EMG activity and that of force production. Reproducibility of the reaction time variables was satisfactory (r = 0.79-0.89; coefficient of variation = 8.8%-11.6%). The TRT was 0.121 s, SD 0.014 in FLG and 0.119 s, SD 0.011 in RLG. The MT ranged from 0.008 s, SD 0.009 (GM) to 0.057 s, SD 0.050 (GA) in FLG and from 0.018 s, SD 0.029 (GA) to 0.045 s, SD 0.009 (GM) in RLG. In some individual cases there were no MT values before horizontal force production.(ABSTRACT TRUNCATED AT 250 WORDS)