Physiological modules for generating discrete and rhythmic movements: Component analysis of EMG signals

A central question in Neuroscience is that of how the nervous system generates the spatiotemporal commands needed to realize complex gestures, such as handwriting. A key postulate is that the central nervous system (CNS) builds up complex movements from a set of simpler motor primitives or control modules. In this study we examined the control modules underlying the generation of muscle activations when performing different types of movement: discrete, point-to-point movements in eight different directions and continuous figure-eight movements in both the normal, upright orientation and rotated 90 degrees. To test for the effects of biomechanical constraints, movements were performed in the frontal-parallel or sagittal planes, corresponding to two different nominal flexion/abduction postures of the shoulder. In all cases we measured limb kinematics and surface electromyographic activity (EMB) signals for seven different muscles acting around the shoulder. We first performed principal component analysis (PCA) of the EMG signals on a movement-by-movement basis. We found a surprisingly consistent pattern of muscle groupings across movement types and movement planes, although we could detect systematic differences between the PCs derived from movements performed in each sholder posture and between the principal components associated with the different orientations of the figure. Unexpectedly we found no systematic differences between the figute eights and the point-to-point movements. The first three principal components could be associated with a general co-contraction of all seven muscles plus two patterns of reciprocal activatoin. From these results, we surmise that both "discrete-rhythmic movements" such as the figure eight, and discrete point-to-point movement may be constructed from three different fundamental modules, one regulating the impedance of the limb over the time span of the movement and two others operating to generate movement, one aligned with the vertical and the other aligned with the horizontal.

A simulation of the surface EMG for analysis of muscle activity during whole body vibratory stimulation

Many studies have accounted for whole body vibration effects in the fields of exercise physiology, sport and rehabilitation medicine. Generally, surface EMG is utilized to assess muscular activity during the treatment; however, large motion artifacts appear superimposed to the raw signal, making sEMG recording not suitable before any artifact filtering. Sharp notch filters, centered at vibration frequency and at its superior harmonics, have been used in previous studies, to remove the artifacts. [6, 10] However, to get rid of those artifacts some true EMG signal is lost. The purpose of this study was to reproduce the effect of motor-unit synchronization on a simulated surface EMG during vibratory stimulation. In addition, authors mean to evaluate the EMG power percentage in those bands in which are also typically located motion artifact components. Model characteristics were defined to take into account two main aspect: the muscle MUs discharge behavior and the triggering effects that appear during local vibratory stimulation. [7] Inter-pulse-interval, was characterized by a polimodal distribution related to the MU discharge frequency (IPI 55-80ms, σ=12ms) and to the correlation with the vibration period within the range of ±2 ms due to vibration stimulus. [1, 7] The signals were simulated using different stimulation frequencies from 30 to 70 Hz. The percentage of the total simulated EMG power within narrow bands centered at the stimulation frequency and its superior harmonics (± 1 Hz) resulted on average about 8% (± 2.85) of the total EMG power. However, the artifact in those bands may contain more than 40% of the total power of the total signal. [6] Our preliminary results suggest that the analysis of the muscular activity of muscle based on raw sEMG recordings and RMS evaluation, if not processed during vibratory stimulation may lead to a serious overestimation of muscular response.

Teleoperation of a 7 DOF humanoid robot arm using human arm accelerations and EMG signals

We present our work on tele-operating a complex humanoid robot with the help of bio-signals collected from the operator. The frameworks (for robot vision, collision avoidance and machine learning), developed in our lab, allow for a safe interaction with the environment, when combined. This even works with noisy control signals, such as, the operator’s hand acceleration and their electromyography (EMG) signals. These bio-signals are used to execute equivalent actions (such as, reaching and grasping of objects) on the 7 DOF arm.

Short- and long-term changes in joint co-contraction associated with motor learning as revealed from surface EMG.

In the field of motor control, two hypotheses have been controversial: whether the brain acquires internal models that generate accurate motor commands, or whether the brain avoids this by using the viscoelasticity of musculoskeletal system. Recent observations on relatively low stiffness during trained movements support the existence of internal models. However, no study has revealed the decrease in viscoelasticity associated with learning that would imply improvement of internal models as well as synergy between the two hypothetical mechanisms. Previously observed decreases in electromyogram (EMG) might have other explanations, such as trajectory modifications that reduce joint torques. To circumvent such complications, we required strict trajectory control and examined only successful trials having identical trajectory and torque profiles. Subjects were asked to perform a hand movement in unison with a target moving along a specified and unusual trajectory, with shoulder and elbow in the horizontal plane at the shoulder level. To evaluate joint viscoelasticity during the learning of this movement, we proposed an index of muscle co-contraction around the joint (IMCJ). The IMCJ was defined as the summation of the absolute values of antagonistic muscle torques around the joint and computed from the linear relation between surface EMG and joint torque. The IMCJ during isometric contraction, as well as during movements, was confirmed to correlate well with joint stiffness estimated using the conventional method, i.e., applying mechanical perturbations. Accordingly, the IMCJ during the learning of the movement was computed for each joint of each trial using estimated EMG-torque relationship. At the same time, the performance error for each trial was specified as the root mean square of the distance between the target and hand at each time step over the entire trajectory. The time-series data of IMCJ and performance error were decomposed into long-term components that showed decreases in IMCJ in accordance with learning with little change in the trajectory and short-term interactions between the IMCJ and performance error. A cross-correlation analysis and impulse responses both suggested that higher IMCJs follow poor performances, and lower IMCJs follow good performances within a few successive trials. Our results support the hypothesis that viscoelasticity contributes more when internal models are inaccurate, while internal models contribute more after the completion of learning. It is demonstrated that the CNS regulates viscoelasticity on a short- and long-term basis depending on performance error and finally acquires smooth and accurate movements while maintaining stability during the entire learning process.

Accurate real-time feedback of surface EMG during fMRI.

Real-time acquisition of EMG during functional MRI (fMRI) provides a novel method of controlling motor experiments in the scanner using feedback of EMG. Because of the redundancy in the human muscle system, this is not possible from recordings of joint torque and kinematics alone, because these provide no information about individual muscle activation. This is particularly critical during brain imaging because brain activations are not only related to joint torques and kinematics but are also related to individual muscle activation. However, EMG collected during imaging is corrupted by large artifacts induced by the varying magnetic fields and radio frequency (RF) pulses in the scanner. Methods proposed in literature for artifact removal are complex, computationally expensive, and difficult to implement for real-time noise removal. We describe an acquisition system and algorithm that enables real-time acquisition for the first time. The algorithm removes particular frequencies from the EMG spectrum in which the noise is concentrated. Although this decreases the power content of the EMG, this method provides excellent estimates of EMG with good resolution. Comparisons show that the cleaned EMG obtained with the algorithm is, like actual EMG, very well correlated with joint torque and can thus be used for real-time visual feedback during functional studies.

EMG Biofeedback Videogame System for the Gait Rehabilitation of Hemiparetic Individuals

Gemstone Team CHIP

Decomposition of surface electromyographic signal using Hidden Markov Model

The detection of physiological signals from the motor system (electromyographic signals) is being utilized in the practice clinic to guide the therapist in a more precise and accurate diagnosis of motor disorders. In this context, the process of decomposition of EMG (electromyographic) signals that includes the identification and classification of MUAP (Motor Unit Action Potential) of a EMG signal, is very important to help the therapist in the evaluation of motor disorders. The EMG decomposition is a complex task due to EMG features depend on the electrode type (needle or surface), its placement related to the muscle, the contraction level and the health of the Neuromuscular System. To date, the majority of researches on EMG decomposition utilize EMG signals acquired by needle electrodes, due to their advantages in processing this type of signal. However, relatively few researches have been conducted using surface EMG signals. Thus, this article aims to contribute to the clinical practice by presenting a technique that permit the decomposition of surface EMG signal via the use of Hidden Markov Models. This process is supported by the use of differential evolution and spectral clustering techniques. The developed system presented coherent results in: (1) identification of the number of Motor Units actives in the EMG signal; (2) presentation of the morphological patterns of MUAPs in the EMG signal; (3) identification of the firing sequence of the Motor Units. The model proposed in this work is an advance in the research area of decomposition of surface EMG signals.

Reliability of surface EMG measurements of the quadriceps during maximal isometric contractions following water immersion

Background: The influence of water immersion on neuromuscular function is of importance to a number of disciplines; however, the reliability of surface electromyography (SEMG) following water immersion is not known. This study examined the reliability of SEMG amplitude during maximal voluntary isometric contractions (MVICs) of the vastus lateralis following water immersion.

Methods: Using a Biodex isokinetic dynamometer and in a randomized order, 12 healthy male subjects performed four MVICs at 60° knee flexion on both the dominant and nondominant kicking legs, and the SEMG was recorded. Each subject's dominant and nondominant kicking leg was then randomly assigned to have SEMG electrodes removed or covered during 15 min of water immersion (20°C–25°C). Following water immersion, subjects performed a further four MVICs.

Results: Intraclass correlation coefficient (ICC) and the relative standard error of measurement (%SEM) of SEMG amplitude showed moderate to high trial-to-trial reliability when electrodes were covered (0.93% and 2.79%) and removed (0.95% and 2.10%, respectively).

Conclusions: The results of the this study indicate that SEMG amplitude of the vastus lateralis may be accurately determined during maximal voluntary contractions following water immersion if electrodes are either removed or covered with water-resistive tape during the immersion.

Azimuthal asymmetry in the risetime of the surface detector signals of the Pierre Auger Observatory.

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, (sec theta)(max), sensitive to the mass composition of cosmic rays above 3 x 10(18) eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modeling that must be resolved before the mass composition can be inferred from (sec theta)(max).

Simulation of surface EMG for the analysis of muscle activity during whole body vibratory stimulation

This study aims to reproduce the effect of motor-unit synchronization on surface EMG recordings during vibratory stimulation to highlight vibration evoked muscle activity. The authors intended to evaluate, through numerical simulations, the changes in surface EMG spectrum in muscles undergoing whole body vibration stimulation. In some specific bands, in fact, vibration induced motion artifacts are also typically present. In addition, authors meant to compare the simulated EMGs with respect to real recordings in order to discriminate the effect of synchronization of motor units discharges with vibration frequencies from motion artifacts. Computations were performed using a model derived from previous studies and modified to consider the effect of vibratory stimulus, the motor unit synchronization and the endplates-electrodes relative position on the EMG signal. Results revealed that, in particular conditions, synchronization of MUs' discharge generates visible peaks at stimulation frequency and its harmonics. However, only a part of the total power of surface EMGs might be enclosed within artifacts related bands (±1. Hz centered at the stimulation frequency and its superior harmonics) even in case of strong synchronization of motor units discharges with the vibratory stimulus. © 2013 Elsevier Ireland Ltd.

Back extensor muscle fatigue at submaximal workloads assessed using frequency banding of the electromyographic signal

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Muscle fiber conduction velocity of the upper trapezius muscle during dynamic contraction of the upper limb in patients with chronic neck pain

The purpose of this study was to compare average muscle fiber conduction velocity (CV) and its changes over time in the upper trapezius muscle during a repetitive upper limb task in people with chronic neck pain and in healthy controls. Surface EMG signals were detected bilaterally from the upper trapezius muscle of 19 patients and nine healthy controls using linear adhesive arrays of four electrodes. Subjects were asked to tap their hands in a cyclic manner between targets positioned mid-thigh and 120 degrees of shoulder flexion, to the beat of a metronome set at 88 beats/min for up to 5 min. Muscle fiber CV and instantaneous mean power spectral frequency were estimated for each cycle at the time instant corresponding to 90 degrees of shoulder flexion. Average muscle fiber CV of the upper trapezius muscle was higher in people with chronic neck pain (mean +/- SE, 4.8 +/- 0.1 m/s) than in control subjects (4.4 +/- 0.1 m/s; P

Diagnostic accuracy of the EMG parameters associated with anterior knee pain in the diagnosis of patellofemoral pain syndrome

Objective To assess the diagnostic accuracy of the surface electromyography (sEMG) parameters associated with referred anterior knee pain in diagnosing patellofemoral pain syndrome (PFPS). Design Sensitivity and specificity analysis. Setting Physical rehabilitation center and laboratory of biomechanics and motor control. Participants Pain-free subjects (n=29) and participants with PFPS (n=22) selected by convenience. Interventions Not applicable. Main Outcome Measure The diagnostic accuracy was calculated for sEMG parameters’ reliability, precision, and ability to differentiate participants with and without PFPS. The selected sEMG parameter associated with anterior knee pain was considered as an index test and was compared with the reference standard for the diagnosis of PFPS. Intraclass correlation coefficient, SEM, independent t tests, sensitivity, specificity, negative and positive likelihood ratios, and negative and positive predictive values were used for the statistical analysis. Results The medium-frequency band (B2) parameter was reliable (intraclass correlation coefficient=.80–.90), precise (SEM=2.71–3.87 normalized unit), and able to differentiate participants with and without PFPS (P<.05). The association of B2 with anterior knee pain showed positive diagnostic accuracy values (specificity, .87; sensitivity, .70; negative likelihood ratio, .33; positive likelihood ratio, 5.63; negative predictive value, .72; and positive predictive value, .86). Conclusions The results provide evidence to support the use of EMG signals (B2 – frequency band of 45–96Hz) of the vastus lateralis and vastus medialis muscles with referred anterior knee pain in the diagnosis of PFPS.