Simulation of electromyographic signals

A computational model of surface electromyography signal alterations after spinal cord injury

Objective. Spinal cord injury (SCI) can cause significant impairment and disability with an impact on the quality of life for individuals with SCI and their caregivers. Surface electromyography (sEMG) is a sensitive and non-invasive technique to measure muscle activity and has demonstrated great potential in capturing neuromuscular changes resulting from SCI. The mechanisms of the sEMG signal characteristic changes due to SCI are multi-faceted and difficult to studyin vivo. In this study, we utilized well-established computational models to characterize changes in sEMG signal after SCI and identify sEMG features that are sensitive and specific to different aspects of the SCI.Approach. Starting from existing models for motor neuron pool organization and motor unit action potential generation for healthy neuromuscular systems, we implemented scenarios to model damages to upper motor neurons, lower motor neurons, and the number of muscle fibers within each motor unit. After simulating sEMG signals from each scenario, we extracted time and frequency domain features and investigated the impact of SCI disruptions on sEMG features using the Kendall Rank Correlation analysis.Main results. The commonly used amplitude-based sEMG features (such as mean absolute values and root mean square) cannot differentiate between injury scenarios, but a broader set of features (including autoregression and cepstrum coefficients) provides greater specificity to the type of damage present.Significance. We introduce a novel approach to mechanistically relate sEMG features (often underused in SCI research) to different types of neuromuscular alterations that may occur after SCI. This work contributes to the further understanding and utilization of sEMG in clinical applications, which will ultimately improve patient outcomes after SCI.

Methodology of surface electromyography in gait analysis: review of the literature

Gait analysis is a significant diagnostic procedure for the clinicians who manage musculoskeletal disorders. Surface electromyography (sEMG) combined with kinematic and kinetic data is a useful tool for decision making of the appropriate method needed to treat such patients. sEMG has been used for decades to evaluate neuromuscular responses during a range of activities and develop rehabilitation protocols. The sEMG methodology followed by researchers assessed the issues of noise control, wave frequency, cross talk, low signal reception, muscle co-contraction, electrode placement protocol and procedure as well as EMG signal timing, intensity and normalisation so as to collect accurate, adequate and meaningful data. Further research should be done to provide more information related to the muscle activity recorded by sEMG and the force produced by the corresponding muscle during gait analysis.

Generation of synthetic surface electromyography signals under fatigue conditions for varying force inputs using feedback control algorithm

Surface electromyography is a non-invasive technique used for recording the electrical activity of neuromuscular systems. These signals are random, complex and multi-component. There are several techniques to extract information about the force exerted by muscles during any activity. This work attempts to generate surface electromyography signals for various magnitudes of force under isometric non-fatigue and fatigue conditions using a feedback model. The model is based on existing current distribution, volume conductor relations, the feedback control algorithm for rate coding and generation of firing pattern. The result shows that synthetic surface electromyography signals are highly complex in both non-fatigue and fatigue conditions. Furthermore, surface electromyography signals have higher amplitude and lower frequency under fatigue condition. This model can be used to study the influence of various signal parameters under fatigue and non-fatigue conditions.

Experimental and modelling investigation of surface EMG spike analysis

A pattern classification method based on five measures extracted from the surface electromyographic (sEMG) signal is used to provide a unique characterization of the interference pattern for different motor unit behaviours. This study investigated the sensitivity of the five sEMG measures during the force gradation process. Tissue and electrode filtering effects were further evaluated using a sEMG model. Subjects (N=8) performed isometric elbow flexion contractions from 0 to 100% MVC. The sEMG signals from the biceps brachii were recorded simultaneously with force. The basic building block of the sEMG model was the detection of single fibre action potentials (SFAPs) through a homogeneous, equivalent isotropic, infinite volume conduction medium. The SFAPs were summed to generate single motor unit action potentials. The physiologic properties from a well-known muscle model and motor unit recruitment and firing rate schemes were combined to generate synthetic sEMG signals. The following pattern classification measures were calculated: mean spike amplitude, mean spike frequency, mean spike slope, mean spike duration, and the mean number of peaks per spike. Root-mean-square amplitude and mean power frequency were also calculated. Taken together, the experimental data and modelling analysis showed that below 50% MVC, the pattern classification measures were more sensitive to changes in force than traditional time and frequency measures. However, there are additional limitations associated with electrode distance from the source that must be explored further. Future experimental work should ensure that the inter-electrode distance is no greater than 1cm to mitigate the effects of tissue filtering.

A muscle architecture model offering control over motor unit fiber density distributions

The aim of this study was to develop a muscle architecture model able to account for the observed distributions of innervation ratios and fiber densities of different types of motor units in a muscle. A model algorithm is proposed and mathematically analyzed in order to obtain an inverse procedure that allows, by modification of input parameters, control over the output distributions of motor unit fiber densities. The model's performance was tested with independent data from a glycogen depletion study of the medial gastrocnemius of the rat. Results show that the model accurately reproduces the observed physiological distributions of innervation ratios and fiber densities and their relationships. The reliability and accuracy of the new muscle architecture model developed here can provide more accurate models for the simulation of different electromyographic signals.

A wavelet-based adaptive filter for removing ECG interference in EMGdi signals

Diaphragmatic electromyogram (EMGdi) signals convey important information on respiratory diseases. In this paper, an adaptive filter for removing the electrocardiographic (ECG) interference in EMGdi signals based on wavelet theory is proposed. Power spectrum analysis was performed to evaluate the proposed method. Simulation results show that the power spectral density (PSD) of the extracted EMGdi signal from an ECG corrupted signal is within 1.92% average error relative to the original EMGdi signal. Testing on clinical EMGdi data confirm that this method is also efficient in removing ECG artifacts from the corrupted clinical EMGdi signal.

A simulation study for a surface EMG sensor that detects distinguishable motor unit action potentials

An advanced volume conductor model was used to simulate the surface-detected motor unit action potentials (MUAPs) due to current sources located at different depths within the muscle tissue of the biceps brachii. Seven different spatial filters were investigated by linear summation of the monopolarly detected surface MUAPs on a square array of nine electrodes. The criterion of the relative energy-of-difference (EOD) between the MUAPs was used to rank spatial filters for their ability to distinguish two motor units located at different depths. Using the same criterion pair wise combinations of spatial filters were ranked for their ability to generate different MUAP shape representations of the same motor unit. In both analyses, the bi-transversal double-differential (BiTDD) configurations and pair wise combinations involving a BiTDD configuration consistently ranked highest. Varying electrode spacing did not change the results in a relevant way. Based on the EOD calculations, a four-channel detection system using all available electrodes of the array is proposed. The implications of using only six electrodes, effectively reducing contact area of the sensor in half, are discussed.

Simulation analysis of interference EMG during fatiguing voluntary contractions. Part I: What do the intramuscular spike amplitude–frequency histograms reflect?

Decline in amplitude of EMG signals and in the rate of counts of intramuscularly recorded spikes during fatigue is often attributed to a progressive reduction of the neural drive only. As a rule, alterations in intracellular action potential (IAP) are not taken into account. To test correctness of the hypothesis, the effect of various discharge frequency patterns as well as changes in IAP shape and muscle fibre propagation velocity (MFPV) on the spike amplitude-frequency histogram of intramuscular interference EMG signals were simulated and analyzed. It was assumed that muscle was composed of four types of motor units (MUs): slow-twitch fatigue resistant, fast-twitch fatigue resistant, fast intermediate, and fast fatigable. MFPV and IAP duration at initial stage before fatigue as well as their changes differed for individual MU types. Fatigability of individual MU types in normal conditions as well as in the case of ischaemic or low oxygen conditions due to restricted blood flow was also taken into account. It was found that spike amplitude-frequency histogram is poorly sensitive to MU firing frequency, while it is highly sensitive to IAP profile lengthening. It is concluded that spike amplitude-frequency analysis can hardly provide a correct measure of MU rate-coding pattern during fatigue.

Experimental and modeling investigation of spectral compression of biceps brachii SEMG activity with increasing force levels

This study investigated possible motor unit (MU) firing patterns underlying changes in biceps brachii (BB) surface electromyographic (SEMG) activity in 96 participants who performed isometric actions of the elbow flexors at 40%, 60%, 80%, and 100% of maximum voluntary contraction (MVC). We also conducted a modeling investigation to determine the extent to which a model would fit the experimental results. Experimentally, there was a linear increase (277%; p<0.01) in root-mean-square (RMS) amplitude with increasing force. The mean power frequency (MNF) remained stable from 40% to 80% of MVC, but there was a decrease (8.2%; p<0.01) between 80% and 100% of MVC. A modeling approach was taken wherein well-known recruitment and rate-coding schemes activated MUs whose basic building block was the muscle fibre action potential. Two conditions were investigated: (1) an increase in firing rate (rate-coding) and (2) synchronization. The levels of rate-coding and synchronization were selected to produce a linear RMS-force relationship as observed in the experimental data. Then, the impact of these two strategies on changes in MNF was assessed. The MNF remained stable from 40% to 80% of maximum excitation for both the rate-coding and synchronization conditions. There was a decrease in MNF between 80% and 100% of maximum excitation for both modeling conditions, similar to that observed for the experimental data. Thus, at these high forces at which experimental data are technically difficult to obtain, the model supports the idea that both rate-coding and synchronization are responsible for the changes observed in surface EMG amplitude and frequency characteristics.

Modelling fibrillation potentials--analysis of time parameters in the muscle intracellular action potential

A single fiber action potential (SFAP) can be modelled as the convolution of a biolectrical source and a filter impulse response. In the Dimitrov-Dimitrova (D-D) convolutional model, the first temporal derivative of the intracellular action potential (IAP) is used as the source, and Tspl is a time parameter related to the duration of the IAP waveform. This paper is centred on the relation between Tspl and the main spike duration (MSD), defined as the time interval between the first and third phases of the SFAP. We show that Tspl essentially determines the MSD parameter. As experimental data, we used fibrillation potentials (FPs) of two different muscles to study the D-D model. We found that Tspl should have a certain statistical variability in order to explain the variability in the MSD of our FPs. In addition, we present a method to estimate the Tspl values corresponding to a given SFAP from its measured MSD.

Modeling Fibrillation Potentials—a New Analytical Description for the Muscle Intracellular Action Potential

The single-fiber action potential (SFAP) can be modeled as a convolution of a biolectrical source (the excitation) and a transfer function, representing the electrical volume conduction. In the Dimitrov-Dimitrova (D-D) convolutional model, the first temporal derivative of the intracellular action potential (IAP) is used as the source. In this model, the ratio between the amplitudes of the second and first phases of the SFAP (which we call the PPR, after peak-to-peak ratio) increases invariably with radial distance. This is not the case of real recorded fibrillation potentials (FPs). Moreover, FPs show a wider PPR range than that which the D-D model can provide. These discrepancies suggest that the D-D model should be revised. Since the volume conduction parameters seem to have no apparent effects on the PPR, we assume that the origin of this difference lies in the excitation source. This paper presents a new analytical description of the IAP based on that expressed in the D-D model. The new approximation is shown to model FPs with a range of PPRs comparable to that observed in a set of real FPs which we used as our test signals.

A simulation model of the surface EMG signal for analysis of muscle activity during the gait cycle

This work describes a model able to synthetize the surface EMG (electromyography) signal acquired from tibialis anterior and gastrocnemious medialis muscles during walking of asymptomatic adult subjects. The model assumes a muscle structure where the volume conductor is represented by multiple layers of anisotropic media. This model originates from analysis of the single fiber action potential characterized by the conduction velocity. The surface EMG of voluntary contraction is calculated by gathering motor unit action potentials estimated by the summation of all activities of muscle fibers assumed to have a uniformly parallel distribution. The parameters related to the gait cycle, such as onset and cessation timings of muscle activation, amplitude of muscle contraction, periods and sequences of motor units' recruitment, are included in the model presented. In addition, the relative positions of the electrodes during gait can also be specified in order to adapt the simulation to the different acquisition settings.

Physiologically based simulation of clinical EMG signals

An algorithm that generates electromyographic (EMG) signals consistent with those acquired in a clinical setting is described. Signals are generated using a model constructed to closely resemble the physiology and morphology of skeletal muscle, combined with line source models of commonly used needle electrodes positioned in a way consistent with clinical studies. The validity of the simulation routines is demonstrated by comparing values of statistics calculated from simulated signals with those from clinical EMG studies of normal subjects. The simulated EMG signals may be used to explore the relationships between muscle structure and activation and clinically acquired EMG signals. The effects of motor unit (MU) morphology, activation, and neuromuscular junction activity on acquired signals can be analyzed at the fiber, MU and muscle level. Relationships between quantitative features of EMG signals and muscle structure and activation are discussed.

Surface electromyogram signal modelling

The paper reviews the fundamental components of stochastic and motor-unit-based models of the surface electromyogram (SEMG). Stochastic models used in ergonomics and kinesiology consider the SEMG to be a stochastic process whose amplitude is related to the level of muscle activation and whose power spectral density reflects muscle conduction velocity. Motor-unit-based models for describing the spatio-temporal distribution of individual motor-unit action potentials throughout the limb are quite robust, making it possible to extract precise information about motor-unit architecture from SEMG signals recorded by multi-electrode arrays. Motor-unit-based models have not yet been proven as successful, however, for extracting information about recruitment and firing rates throughout the full range of contraction. The relationship between SEMG and force during natural dynamic movements is much too complex to model in terms of single motor units.

A fast and reliable technique for muscle activity detection from surface EMG signals

The estimation of on-off timing of human skeletal muscles during movement is an important issue in surface electromyography (EMG) signal processing with relevant clinical applications. In this paper, a novel approach to address this issue is proposed. The method is based on the identification of single motor unit action potentials from the surface EMG signal with the use of the continuous wavelet transform. A manifestation variable is computed as the maximum of the outputs of a bank of matched filters at different scales. A threshold is applied to the manifestation variable to detect EMG activity. A model, based on the physical structure of the muscle, is used to test the proposed technique on synthetic signals with known features. The resultant bias of the onset estimate is lower than 40 ms and the standard deviation lower than 30 ms in case of additive colored Gaussian noise with signal-to-noise ratio as low as 2 dB. Comparison with previously developed methods was performed, and representative applications to experimental signals are presented. The method is designed for a complete real-time implementation and, thus, may be applied in clinical routine activity.

Motor unit recruitment strategies investigated by surface EMG variables

During isometric contractions of increasing strength, motor units (MUs) are recruited by the central nervous system in an orderly manner starting with the smallest, with muscle fibers that usually show the lowest conduction velocity (CV). Theory predicts that the higher the velocity of propagation of the action potential, the higher the power at high frequencies of the detected surface signal. These considerations suggest that the power spectral density of the surface detected electromyogram (EMG) signal may give indications about the MU recruitment process. The purpose of this paper is to investigate the potential and limitations of spectral analysis of the surface EMG signal as a technique for the investigation of muscle force control. The study is based on a simulation approach and on an experimental investigation of the properties of surface EMG signals detected from the biceps brachii during isometric linearly increasing torque contractions. Both simulation and experimental data indicate that volume conductor properties play an important role as confounding factors that may mask any relation between EMG spectral variables and estimated CV as a size principle parameter during ramp contractions. The correlation between spectral variables and CV is thus significantly lower when the MU pool is not stable than during constant-torque isometric contractions. Our results do not support the establishment of a general relationship between spectral EMG variables and torque or recruitment strategy.

A novel approach for precise simulation of the EMG signal detected by surface electrodes

We propose a new electromyogram generation and detection model. The volume conductor is described as a nonhomogeneous (layered) and anisotropic medium constituted by muscle, fat and skin tissues. The surface potential detected in space domain is obtained from the application of a two-dimensional spatial filter to the input current density source. The effects of electrode configuration, electrode size and inclination of the fibers with respect to the detection system are included in the transfer function of the filter. Computation of the signal in space domain is performed by applying the Radon transform; this permits to draw considerations about spectral dips and clear misunderstandings in previous theoretical derivations. The effects of generation and extinction of the action potentials at the fiber end plate and at the tendons are included by modeling the source current, without any approximation of its shape, as a function of space and time and by using again the Radon transform. The approach, based on the separation of the temporal and spatial properties of the muscle fiber action potential and of the volume conductor, includes the capacitive tissue properties.

EMG signal decomposition: how can it be accomplished and used?

Electromyographic (EMG) signals are composed of the superposition of the activity of individual motor units. Techniques exist for the decomposition of an EMG signal into its constituent components. Following is a review and explanation of the techniques that have been used to decompose EMG signals. Before describing the decomposition techniques, the fundamental composition of EMG signals is explained and after, potential sources of information from and various uses of decomposed EMG signals are described.

Surface EMG models: properties and applications

After a general introduction on the kind of models and the use of models in the natural sciences, the main body of this paper reviews potential properties of structure based surface EMG (sEMG) models. The specific peculiarities of the categories (i) source description, (ii) motor unit structure, (iii) volume conduction, (iv) recording configurations and (v) recruitment and firing behaviour are discussed. For a specific goal, not all aspects conceivable have to be part of a model description. Therefore, finally an attempt is made to integrate the 'question level' and the 'model property level' in a matrix providing direction to the development and application of sEMG models with different characteristics and varying complexity. From this overview it appears that the least complex are models describing how the morphological muscle features are reflected in multi-channel EMG measurements. The most challenging questions in terms of model complexity are related to supporting the diagnosis of neuromuscular disorders.

New aspects to event-synchronous cancellation of ECG interference: an application of the method in diaphragmatic EMG signals

An "event-synchronous interference canceller" (ESC) for cancellation of electrocardiographic (ECG) interference in diaphragmatic electromyographic (EMGdi) signals is addressed in this paper. ESC pursues the concept of the "event synchronous adaptive interference canceller" (ESAIC), which was proposed in [1] as a specific application of the well known "adaptive noise canceller" (ANC) paradigm, but ESC uses a simple adaptive gain control (AGC) instead of the complex adaptive filter of the ANC. The proposed ESC method is evaluated using both computer simulations and real EMGdi data, and its efficiency in interference cancellation is compared to that of ESAIC. Of particular interest is the result that the ESC can replace the ESAIC providing better performance as well as a considerable reduction of computational costs.

A model of EMG generation

Simulation models are unavoidable in experimental research when the point is to develop new processing algorithms to be applied on real signals in order to extract specific parameter values. Such algorithms have generally to be optimized by comparing true parameter values to those deduced from the algorithm. Only a simulation model can allow the user to access and control the actual process parameter values. This constraint is especially true when dealing with biomedical signals like surface electromyogram (SEMG). This work is an attempt to produce an efficient SEMG simulation model as a help for assessing algorithms related to SEMG features description. It takes into account the most important parameters which could influence these characteristics. This model includes all transformations from intracellular potential to surface recordings as well as a fast implementation of the extracellular potential computation. In addition, this model allows multiple graphically-programmable electrode-set configurations and SEMG simulation in both voluntary and elicited contractions.

Detecting single fiber contributions to motor unit action potentials

The ability to detect muscle fiber action potential (MFAP) contributions to motor unit action potentials (MUAPs) measured using single fiber (SF) and concentric needle (CN) electrodes was studied using simulated MUAPs. Various MFAP-acceleration thresholds were used to define significant fiber contributions. Attempts to detect the significant MFAP contributions, by locating peaks in filtered MUAPs or MUAP accelerations using various MUAP-based thresholds, were then made. Considering filtered MUAPs and a significant contribution threshold of 7.5 kV/s2, and using fiber-density peak-detection criteria, at best 46% and 50% of significant MFAP contributions were detected for the SF and CN MUAPs, respectively. Considering MUAP accelerations and a significant contribution threshold of 7.5 kV/s2, 80% and 84% of significant MFAP contributions could be detected, respectively. Most significant contributions were created from fibers located within approximately 350 microm of the electrode. The results suggest that significant peaks, defined using MUAP-based thresholds, within the acceleration of CN MUAPs can strongly correspond to individual fiber activity and may be useful for measuring fiber density and neuromuscular jitter.