Denoising of HD-sEMG signals using canonical correlation analysis.

High-density surface electromyography (HD-sEMG) is a recent technique that overcomes the limitations of monopolar and bipolar sEMG recordings and enables the collection of physiological and topographical informations concerning muscle activation. However, HD-sEMG channels are usually contaminated by noise in an heterogeneous manner. The sources of noise are mainly power line interference (PLI), white Gaussian noise (WGN) and motion artifacts (MA). The spectral components of these disruptive signals overlap with the sEMG spectrum which makes classical filtering techniques non effective, especially during low contraction level recordings. In this study, we propose to denoise HD-sEMG recordings at 20 % of the maximum voluntary contraction by using a second-order blind source separation technique, named canonical component analysis (CCA). For this purpose, a specific and automatic canonical component selection, using noise ratio thresholding, and a channel selection procedure for the selective version (sCCA) are proposed. Results obtained from the application of the proposed methods (CCA and sCCA) on realistic simulated data demonstrated the ability of the proposed approach to retrieve the original HD-sEMG signals, by suppressing the PLI and WGN components, with high accuracy (for five different simulated noise dispersions using the same anatomy). Afterward, the proposed algorithms are employed to denoise experimental HD-sEMG signals from five healthy subjects during biceps brachii contractions following an isometric protocol. Obtained results showed that PLI and WGN components could be successfully removed, which enhances considerably the SNR of the channels with low SNR and thereby increases the mean SNR value among the grid. Moreover, the MA component is often isolated on specific estimated sources but requires additional signal processing for a total removal. In addition, comparative study with independent component analysis, CCA-wavelet and CCA-empirical mode decomposition (EMD) proved a higher efficiency of the presented method over existing denoising techniques and demonstrated pointless a second filtering stage for denoising HD-sEMG recordings at this contraction level.

Evaluation of HD-sEMG Probability Density Function deformations in ramp exercise.

The aim of the present study is to propose a subject-specific screening approach of High Density surface EMG (HD-sEMG) Probability Density Function (PDF) shape evolution in experimental conditions following a ramp exercise from 0% to 50% of the Maximum Voluntary Contraction (MVC) during 25 seconds of isometric contractions of the Biceps Brachii from six healthy subjects. This method uses High Order Statistics (HOS), namely the kurtosis and the skewness for PDF shape screening examined on selectively positioned Laplacian sEMG channels obtained on an 8×8 HD-sEMG grid. For each subject, the position of the Laplacian channels was chosen based on the level of muscle activation obtained from the Signal to Noise Ratio (SNR) matrix computed for the 64 sEMG signals of the grid in order to obtain independent Laplacian configurations localized in areas with high SNRs indicating high muscle activation. Afterwards, we used the Principal Component Analysis (PCA) to obtain the principal trend of the kurtosis and the skewness computed from the selected Laplacian signals according to force level variation. The obtained results show a globally common increasing HOS trend according to force increase from 0% to 50% MVC for all the subjects regardless of the anatomical, instrumental and physiological variability that usually strongly influences these trends.

Automatic identification of motor unit action potential trains from electromyographic signals using fuzzy techniques.

A technique is proposed that allows automatic decomposition of electromyographic (EMG) signals into their constituent motor unit action potential trains (MUAPTs). A specific iterative algorithm with a classification method using fuzzy-logic techniques was developed. The proposed classification method takes into account imprecise information, such as waveform instability and irregular firing patterns, that is often encountered in EMG signals. Classification features were determined by the combining of time position and waveform information. Statistical analysis of inter-pulse intervals and spike amplitude provided an accurate estimation of features used in the classification step. Algorithm performance was evaluated using simulated EMG signals composed of up to six different discharging motor units corrupted with white noise. The algorithm was then applied to real signals recorded by a high spatial resolution surface EMG device based on a Laplacian spatial filter. On six groups of 20 simulated signals, the decomposition algorithm performed with a maximum and an average mean error rate of 2.13% and 1.37%, respectively. On real surface EMG signals recorded at different force levels (from 10% to 40% of the maximum voluntary contraction), the algorithm correctly identified 21 MUAPTs, compared with the 29 MUAPTs identified by an experienced neurophysiologist. The efficiency of the decomposition on surface EMG signals makes this method very attractive for non-invasive investigation of physiological muscle properties. However, it can also be used to decompose intramuscularly recorded EMG signals.

Effects of resistance training in humans on neck muscle performance, and electromyogram power spectrum changes.

The purpose of this study was to quantify the neuromuscular cervical adaptations to an 8 week strength training programme. Seven healthy men, with no pathological conditions of the neck, performed a lateral flexion isometric resistance-training programme three times a week. The training sessions consisted of one set of ten contractions, each of 6 s duration, at 60% of the predetermined maximal voluntary isometric torque (MVTim) (warm-up) and two sets of eight contractions, each of 6 s duration, at 80% MVTim. The training effects were evaluated in three ways: muscle size, strength and fatigability. The cross-sectional areas (CSA) of the trapezius (TRP) and sternocleidomastoideus (SCM) muscles were determined using a computerised tomographic scanner. Results showed an increase in the CSA of TRP and SCM muscles after training, 8.8% at C5 level and 6.4% at C7 level for SCM muscle and 12.2% at C7 level for TRP muscle. Strength increased significantly under both isometric and isokinetic conditions (35% and 20%, respectively). Muscle fatigability in lateral flexion was quantified during a sustained isometric contraction at 50% of MVTim. The shift of the mean power frequency of the electromyogram power spectrum density function of SCM muscle toward lower frequencies was less after training (14.6% compared to 6.8%). These results indicate the beneficial effect of a strength-training programme which increases neck muscle size and strength during lateral flexion, and decreases the fatigability of the superficial muscles of the neck.

Motor unit recruitment and EMG power spectra during ramp contractions of a bifunctional muscle.

Surface electromyograms (EMGs) were analysed on the short and long head of the biceps brachii (BBSH and BBLH) during single (F and S) or dual (F+S) flexion and supination tasks. It was confirmed, by the analysis of EMG root-mean-square (RMS) values, that the highest activations of BBSH and BBLH were obtained during a maximal dual task. This study was essentially concerned with the analysis of power spectra data obtained during progressive or ramp contractions (RCs). The shape of the power spectra established during the first second of the RCs differs between F, S and F+S tasks. Differences in mean power frequency (MPF) calculated during RCs would be representative of a recruitment of motor units (MUs) that is, at least partly, task-dependent. In order to compare MPF values calculated from RCs performed under different mechanical conditions (F, S and F+S), MPF-RMS(PSD) relationships have been established (RMS(PSD) being defined as the power spectrum density RMS). Both BBSH and BBLH exhibited initial MPF values higher in supination RC than in flexion RC. Because of plateau values reached at the same level of muscle activation whatever the task performed, the slope of the MPF-RMS(PSD) relationship was lower in S than in F. These results are in favour of MU recruitment that is, at least partly, different in F and in S conditions. Dual submaximal tasks seem to mix the activation of the F and S subpopulations of MUs as revealed by the spectral parameters obtained during F+S ramp contractions. This study could find some implication in the field of muscle rehabilitation or reinforcement.

Reproducibility of kinetics of electromyogram spectrum parameters during dynamic exercise.

During incremental exercise on a cycle ergometer a study was made of the reproducibility of changes in electromyographic activity (EMG) of human quadriceps muscles. Seven subjects performed three periods of incremental exercise either every 2 weeks (G1: four subjects) or 6 weeks (G2: three subjects). Each test was normalized with respect to the maximal aerobic power (MAP) of the subject. It consisted of a quasilinear increase in mechanical power (from 20% to 100% MAP) during 8 min. For rectus femoris muscle activity, changes in total power spectrum (PEMG) and in mean power frequency (MPF) were fitted by a 3rd order polynomial function (named profile) and normalized with respect to the maximal value for PEMG and to the mean value for MPF. A curvilinear increase was found for PEMG. The MPF kinetics varied from one subject to another. These changes were either a continuous increase, or a continuous decrease, or an increase followed by a decrease. Only for G2, was MAP determined before each test in order to update the limits of the test. A good reproducibility of PEMG was shown by its mean magnitude [95.88%, (SD 3.92)] and by computing a mean correlation coefficient between profiles two by two [r2 = 0.948 (SD 0.028) n = 21]. Intraclass coefficient correlation (ICC) calculated for each subject indicated a high level of reproducibility for five of the seven subjects (ICC > 0.80). No clear effect of MAP updating on PEMG and MPF profile reproducibility was observed. Thus it is suggested that MPF kinetics may characterize a subject at a given moment while PEMG kinetics may illustrate a normal profile, and they may both characterize EMG changes for a population during incremental exercise.

Surface electromyogram power spectrum in human quadriceps muscle during incremental exercise.

Spectral electromyographic (EMG) changes in human quadriceps muscles were studied to reinvestigate discrepant results concerning mean power frequency (MPF) changes during dynamic exercise. An incremental test consisting of a quasi-linear increase in mechanical power on a bicycle ergometer (for 20-100% of maximal aerobic power) was performed by forty subjects. During this test, surface EMGs from the quadriceps muscles showed that EMG total power (PEMG) increased with a curvilinear pattern for every subject, whereas MPF kinetics varied from one subject to another. PEMG changes had the same shape, which would lead to disappointing results in terms of discrimination between subjects. The ability of normalized MPF kinetics to define significant clusters of subjects was tested using a principal component analysis. This analysis led to the projection of all experiments onto a plane and revealed a relevant grouping of MPF profiles. Differences in MPF kinetics between clusters are interpreted in terms of various possibilities of balance between physiological events leading to an increase or a decrease in MPF.

Electromyogram power spectrum during dynamic contractions at different intensities of exercise.

During dynamic contractions performed on a cycle ergometer, we studied the influence of motor unit (MU) recruitment on the electromyographic (EMG) spectral content by exerting mechanical power of different intensities, which was chosen to remain below the maximal aerobic power (VO2max). The spectral parameters: EMG total power (PEMG), mean (MPF) and median (MED) power frequencies, which are the most representative of the EMG spectral content, were calculated according to the EMG activity of the vastus medialis muscle (VM) and soleus muscle (SOL) of the right leg. For VM and SOL, PEMG increased linearly with exerted power demonstrating an enhancement of MU recruitment. Moreover these relationships were less scattered when exerted power was expressed as a percentage of VO2max. Changes in MPF and MED with varying exercise intensities were different from one subject to another. For a set of subjects, MPF and MED were found to be independent of exerted power. Although VM and SOL muscles are different in fibre type composition, similar results were obtained for both EMG activities. We have concluded that for dynamic contractions performed at different intensities below VO2max, the recruitment of the MU has a poor effect on the EMG spectral content whatever the predominant type of fibre.

The fatigability of two agonistic muscles in human isometric voluntary submaximal contraction: an EMG study. I. Assessment of muscular fatigue by means of surface EMG.

During an external isometric constant torque (25% of the maximal voluntary contraction) maintained until the maximal endurance time (limit time), we analysed and compared the changes in electromyographic (EMG) activity illustrating muscular fatigue simultaneously with mechanical activity (the tangential acceleration theta") related to physiological tremor. The EMG activities recorded were of two agonistic flexors, the biceps brachii (BB) and the brachioradialis (BR) muscles and one of the main extensors, the triceps brachii (TB). The integrated EMG increase and the mean power frequency (MPF) of the power spectrum density function (PSDF) decrease were larger for BR than for BB activity. These two findings suggested a greater BR fatigability. However, it is shown that differences between BB and BR MPF changes could be related to differences in the PSDF upper frequency limit of the two muscles and also to the relative magnitude of their tremor component.

The fatigability of two agonistic muscles in human isometric voluntary submaximal contraction: an EMG study. II. Motor unit firing rate and recruitment.

The recruitment and firing rate of biceps brachii (BB) and brachioradialis (BR) motor units (MUs) were studied in the course of fatiguing isometric contractions at 20%-30% of maximal voluntary contraction (MVC). MU recruitment generally occurred throughout the maintained contraction and was similar for BB and BR muscles. Newly recruited MUs started to discharge in the form of bursts, the duration of which increased until a continuous rhythmical firing was achieved. Within each burst, the first interval between two consecutive discharges was usually the shortest. MU threshold was lowered just after the limit time of the maintained contraction. The MU's firing rate either increased or remained stable as a function of the elapsed time. It is concluded that (1) in fatiguing isometric contractions at 20%-30% MVC contractile failure is mainly compensated for by MU recruitment and a lowered MU threshold and (2) differences between in surface changes in the electromyogram of BB and BR muscles cannot easily be explained by related differences in MU firing rate and recruitment.