Simulation of intramuscular EMG signals detected using implantable myoelectric sensors (IMES)

The purpose of this study was to test the feasibility of recording independent electromyographic (EMG) signals from the forearm using implantable myoelectric sensors (IMES), for myoelectric prosthetic control. Action potentials were simulated using two different volume conductor models: a finite-element (FE) model that was used to explore the influence of the electrical properties of the surrounding inhomogeneous tissues and an analytical infinite volume conductor model that was used to estimate the approximate detection volume of the implanted sensors. Action potential amplitude increased progressively as conducting electrodes, the ceramic electrode casing and high resistivity encapsulation tissue were added to the model. For the muscle fiber locations examined, the mean increase in EMG root mean square amplitude when the full range of material properties was included in the model was 18.2% (+/-8.1%). Changing the orientation of the electrode with respect to the fiber direction altered the shape of the electrode detection volume and reduced the electrode selectivity. The estimated detection radius of the IMES electrode, assuming a cylindrical muscle cross section, was 4.8, 6.2, and 7.5 mm for electrode orientations of 0 degree, 22.5 degrees, and 45 degrees with respect to the muscle fiber direction.     

基于灰度理论模型的截肢受试者手势分类方法研究

针对截肢者手势动作特征提取复杂、动作识别率较低的问题,该文提出一种基于灰度模型的特征提取方法。首先对预处理后的肌电信号与加速度信号经滑动窗信号截取。然后提取表面肌电信号均值、灰度模型的驱动项系数和加速度信号的绝对值均值构成特征向量,最后对滑动窗截取信号特征进行连续的识别。该文采用NinaPro(Non invasive adaptive Prosthetics)公开数据集对提出的方法进行验证,实验表明该文算法能够有效提取肌电和加速度信号的特征,对9名截肢受试者的17类手势动作的平均识别率达到91.14%,提高了17类手势的识别准确率,为仿生假肢人机交互控制算法提供了一种新的思路。     

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.     

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.     

Assessment of average muscle fiber conduction velocity from surface EMG signals during fatiguing dynamic contractions

In this paper, we propose techniques of surface electromyographic (EMG) signal detection and processing for the assessment of muscle fiber conduction velocity (CV) during dynamic contractions involving fast movements. The main objectives of the study are: 1) to present multielectrode EMG detection systems specifically designed for dynamic conditions (in particular, for CV estimation); 2) to propose a novel multichannel CV estimation method for application to short EMG signal bursts; and 3) to validate on experimental signals different choices of the processing parameters. Linear adhesive arrays of electrodes are presented for multichannel surface EMG detection during movement. A new multichannel CV estimation algorithm is proposed. The algorithm provides maximum likelihood estimation of CV from a set of surface EMG signals with a window limiting the time interval in which the mean square error (mse) between aligned signals is minimized. The minimization of the windowed mse function is performed in the frequency domain, without limitation in time resolution and with an iterative computationally efficient procedure. The method proposed is applied to signals detected from the vastus laterialis and vastus medialis muscles during cycling at 60 cycles/min. Ten subjects were investigated during a 4-min cycling task. The method provided reliable assessment of muscle fatigue for these subjects during dynamic contractions.     

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 multiple-layer finite-element model of the surface EMG signal

The effect of skin, muscle, fat, and bone tissue on simulated surface electromyographic (EMG) signals was examined using a finite-element model. The amplitude and frequency content of the surface potential were observed to increase when the outer layer of a homogeneous muscle model was replaced with highly resistive skin or fat tissue. The rate at which the surface potential decreased as the fiber was moved deeper within the muscle also increased. Similarly, the rate at which the surface potential decayed around the surface of the model, for a constant fiber depth, increased. When layers of subcutaneous fat of increasing thickness were then added to the model, EMG amplitude, frequency content, and the rate of decay of the surface EMG signal around the limb decreased, due to the increased distance between the electrodes and the active fiber. The influence of bone on the surface potential was observed to vary considerably, depending on its location. When located close to the surface of the volume conductor, the surface EMG signal between the bone and the source and directly over the bone increased, accompanied by a slight decrease on the side of the bone distal to the active fiber. The results emphasize the importance of distinguishing between the effects of material properties and the distance between source and electrode when considering the influence of subcutaneous tissue, and suggest possible distortions in the surface EMG signal in regions where a bone is located close to the skin surface.     

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     

Singularity characteristics of needle EMG IP signals

Clinical electromyography (EMG) interference pattern (IP) signals can reveal more diagnostic information than their constituents, the motor unit action potentials (MUAPs). Singularities and irregular structures typically characterize the mathematically defined content of information in signals. In this paper, a wavelet transform method is used to detect and quantify the singularity characteristics of EMG IP signals using the Lipschitz exponent (LE) and measures derived from it. The performance of the method is assessed in terms of its ability to discriminate healthy, myopathic and neuropathic subjects and how it compares with traditionally used Turns Analysis (TA) methods and a method recently developed by the authors, interscale wavelet maximum (ISWM). Highly significant intergroup differences were found using the LE method. Most of the singularity measures have a performance similar to that of ISWM and considerably better than that of TA. Some measures such as the ratio of the mean LE value to the number of singular points in the signal have considerably superior performance to both methods. These findings add weight to the view that wavelet analysis methods offer an effective way forward in the quantitative analysis of EMG IP signal to assist the clinician in the diagnosis of neuromuscular disorders.     

Robotic index finger prosthesis using stackable double 4-BAR mechanisms

This paper suggests a robotic index finger prosthesis realized to be one degree-of-freedom by using stackable double 4-bar mechanisms. Also its control method makes use of two electromyographic (EMG) signals measured on skin surfaces of flexor digitorum superficialis (FDS) and extensor indicis (EI) in a lower arm. In this paper, we assume that EMG signals have some relations with velocity of muscle movement by neglecting finger dynamics due to its negligible small mass. In order to obtain desired position and velocity of robotic index finger, the measured raw EMG signals are processed by sequential procedures such as root mean squaring, applying threshold operation to extract the initial burst part, subtracting antagonistic EMG signal, and integrating by every 2 millisecond. Finally the effectiveness of the suggested mechanism and control method is verified through experiments.     

A comparative study of simultaneous vibromyography andelectromyography with active human quadriceps

Vibromyographic (VMG) signals, which are low-frequency vibration signals generated during muscle contraction, were studied in comparison with electromyographic (EMG) signals recorded simultaneously during isometric contraction of the human quadriceps muscles. The comparison was accomplished by evaluating the averaged root mean squared (rms) value, mean frequency (MF), and peak frequency (PF) of the VMG and EMG signals for four muscle contraction levels at joint angles of 30 degrees, 60 degrees, and 90 degrees. The four contraction levels, namely 20, 40, 60, and 80% of maximum voluntary contraction (MVC), were estimated and controlled by the torque readings of a Cybex II dynamometer. It was found that the VMG and EMG under the same conditions on the same muscle group are in general equally sensitive to the levels of muscle contraction. Results show that the rms value of the VMG signal increases linearly, in a manner similar to the EMG rms/%MVC relationship, with increasing muscle contraction levels. Furthermore, the study indicates that the averaged MF (6-24 Hz) and PF (9-19 Hz) of the VMG signals are much lower than the MF (75-109 Hz) and PF (40-80 Hz) of the EMG signals. The slopes of MF/%MVC curves for the VMG and EMG are approximately the same for 60 degrees and 90 degrees joint angles (approximately 3.1 Hz per 20% MVC for VMG and approximately 2.6 Hz per 20% MVC for EMG).(ABSTRACT TRUNCATED AT 250 WORDS)     

Concentric-ring electrode systems for noninvasive detection ofsingle motor unit activity

New recording techniques for detecting surface electromyographic (EMG) signals based on concentric-ring electrodes are proposed in this paper. A theoretical study of the two-dimensional (2-D) spatial transfer function of these recording systems is developed both in case of rings with a physical dimension and in case of line rings. Design criteria for the proposed systems are presented in relation to spatial selectivity. It is shown that, given the radii of the rings, the weights of the spatial filter can be selected in order to improve the rejection of low spatial frequencies, thus increasing spatial selectivity. The theoretical transfer functions of concentric systems are obtained and compared with those of other detection systems. Signals detected with the ring electrodes and with traditional one-dimensional and 2-D systems are compared. The concentric-ring systems show higher spatial selectivity with respect to the traditional detection systems and reduce the problem of electrode location since they are invariant to rotations. The results shown are very promising for the noninvasive detection of single motor unit (MU) activities and decomposition of the surface EMG signal into the constituent MU action potential trains.     

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.     

Estimation of motor unit conduction velocity from surface EMG recordings by signal-based selection of the spatial filters

Muscle fiber conduction velocity (CV) can be estimated by the application of a pair of spatial filters to surface electromagnetic (EMG) signals and compensation of the spatial filter transfer function with equivalent temporal filters. This method integrates the selection of the spatial filters for signal detection to the estimation of CV. Using this approach, in this paper, we propose a novel technique for signal-based selection of the spatial filter pair that minimizes the effect of nonpropagating signal components (end-of-fiber effects) on CV estimates (optimal filters). The technique is applicable to signals with one propagating and one nonpropagating component, such as single motor unit action potentials. It is shown that the determination of the optimal filters also allows the identification of the propagating and nonpropagating signal components. The new method was applied to simulated and experimental EMG signals. Simulated signals were generated by a cylindrical, layered volume conductor model. Experimental signals were recorded from the abductor pollicis brevis with a linear array of 16 electrodes. In the simulations, the proposed approach provided CV estimates with lower bias due to nonpropagating signal components than previously proposed methods based on the entire signal waveform. In the experimental signals, the technique separated propagating and nonpropagating signal components with an average reconstruction error of 2.9 +/- 0.9% of the signal energy. The technique may find application in single motor unit studies for decreasing the variability and bias of CV estimates due to the presence and different weights of the nonpropagating components.     

基于频域块LMS算法的实时虚拟暗室测试方法

针对电磁辐射现场测试被测设备信号和干扰未知的情况,提出了一种基于频域块最小均方算法的实时虚拟暗室测试方法。该方法采用双通道接收机,根据测试通道和背景通道中干扰信号的相关性设计自适应滤波器,在频域对背景通道信号滤波以趋近测试通道中的干扰分量,采用瞬时双通道信号迭代更新滤波器系数,滤波器系数收敛后系统输出中只有被测设备信号。仿真与分析表明,该方法在背景通道有无被测设备信号泄露的情况下都能有效抑制干扰,与基于时域最小均方算法的方法相比,在滤波器长度相同的情况下其计算复杂度更低,适用于实时现场测试。     

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     

Continuous monitoring of sonomyography, electromyography and torque generated by normal upper arm muscles during isometric contraction: sonomyography assessment for arm muscles

The aim of this study is to demonstrate the feasibility of using the continuous signals about the thickness and pennation angle changes of muscles detected in real-time from ultrasound images, named as sonomyography (SMG), to characterize muscles under isometric contraction, along with synchronized surface electromyography (EMG) and generated torque signals. The right biceps brachii muscles of seven normal young adult subjects were tested.We observed that exponential functions could well represent the relationships between the normalized EMG root-mean-square (RMS) and the torque, the RMS and the muscle deformation SMG, and the RMS and the pennation angle SMG for the data of the contraction phase, with exponent coefficients of 0.0341 +/- 0.0148 (Mean SD), 0.0619 +/- 0.0273, and 0.0266 +/- 0.0076, respectively. In addition, the preliminary results also demonstrated linear relationships between the normalized torque and the muscle deformation as well as the pennation angle with the ratios of 9 .79 +/- 3.01 and 2.02 +/- 0.53, respectively. The overall mean R2 for the regressions was approximately 0.9 and the overall mean relative root mean square error (RRMSE) smaller than 15%. The potential values of SMG together with EMG to provide a more comprehensive assessment for the muscle functions should be further investigated with more subjects and more muscle groups.     

Feature Layer Fusion of Linear Features and Empirical Mode Decomposition of Human EMG Signal

To explore the influence of the fusion of different features on recognition, this paper took the electromyography(EMG) signals of rectus femoris under different motions(walk, step, ramp, squat, and sitting) as samples, linear features(time-domain features(variance(VAR) and root mean square(RMS)), frequency-domain features(mean frequency(MF) and mean power frequency(MPF)), and nonlinear features(empirical mode decomposition(EMD)) of the samples were extracted. Two feature fusion algorithms, the s...     

Upper Extremity Limb Function Discrimination Using EMG Signal Analysis

A signal analysis technique is developed for discriminating a set of lower arm and wrist functions using surface EMG signals. Data wete obtained from four electrodes placed around the proximal forearm. The functions analyzed included wrist flexion/extension, wrist abduction/adduction, and forearm pronation/supination. Multivariate autoregression models were derived for each function; discrimination was performed using a multiple-model hypothesis detection technique. This approach extends the work of Graupe and Cline [1] by including spatial correlations and by using a more generalized detection philosophy, based on analysis of the time history of all limb function probabilities. These probabilities are the sufficient statistics for the problem if the EMG data are stationary Gauss-Markov processes. Experimental results on-normal subjects are presented which demonstrate the advantages of using the spatial and time correlation of the signals. This technique should be useful in generating control signals for prosthetic devices.     

Independent component analysis of high-density electromyography in muscle force estimation

Accurate force prediction from surface electromyography (EMG) forms an important methodological challenge in biomechanics and kinesiology. In a previous study (Staudenmann et al., 2006), we illustrated force estimates based on analyses lent from multivariate statistics. In particular, we showed the advantages of principal component analysis (PCA) on monopolar high-density EMG (HD-EMG) over conventional electrode configurations. In the present study, we further improve force estimates by exploiting the correlation structure of the HD-EMG via independent component analysis (ICA). HD-EMG from the triceps brachii muscle and the extension force of the elbow were measured in 11 subjects. The root mean square difference (RMSD) and correlation coefficients between predicted and measured force were determined. Relative to using the monopolar EMG data, PCA yielded a 40% reduction in RMSD. ICA yielded a significant further reduction of up to 13% RMSD. Since ICA improved the PCA-based estimates, the independent structure of EMG signals appears to contain relevant additional information for the prediction of muscle force from surface HD-EMG.