Simulation of Surface EMG Signals for a Multilayer Volume Conductor With Triangular Model of the Muscle Tissue (original) (raw)
Related papers
IEEE Transactions on Biomedical Engineering, 2008
This study analytically describes surface electromyogram (sEMG) signals generated by a model of a triangular muscle, i.e., a muscle with fibres arranged in a fan shape. Examples of triangular muscles in the human body are the deltoid, the pectoralis major, the trapezius, the adductor pollicis. A model of triangular muscle is proposed. It is a sector of a cylindrical volume conductor (with the fibres directed along the radial coordinate) bounded at the muscle/fat interface. The muscle conductivity tensor reflects the fan anisotropy. Edge effects have been neglected. A solution of the non space invariant problem for a triangular muscle is provided in the Fourier domain. An approximate analytical solution for a two plane layer volume conductor model is obtained by introducing a homogeneous layer (modelling the fat) over the triangular muscle. The results are implemented in a complete sEMG generation model (including the finite length of the fibres), simulating single fibre action potentials. The model is not space invariant due to the changes of the volume conductor along the direction of action potential propagation. Thus the detected potentials at the skin surface change shape as they propagate. This determines problems in the extraction and interpretation of parameters. As a representative example of application of the simulation model, the influence of the inhomogeneity of the volume conductor in CV estimation is addressed (for two channels; maximum likelihood and reference point methods). Different fibre depths, electrode placements and small misalignments of the detection system with respect to the fibre have been simulated. The error in CV estimation is large when the depth of the fibre increases, when the detection system is not aligned with the fibre and close to the innervation point and to the tendons.
Estimation of average muscle fiber conduction velocity from simulated surface EMG in pinnate muscles
Journal of Neuroscience Methods, 2007
The aim of this simulation study was to assess the bias in estimating muscle fiber conduction velocity (CV) from surface electromyographic (EMG) signals in muscles with one and two pinnation angles. The volume conductor was a layered medium simulating anisotropic muscle tissue and isotropic homogeneous subcutaneous tissue. The muscle tissue was homogeneous for one pinnation angle and inhomogeneous for bipinnate muscles (two fiber directions). Interference EMG signals were obtained by simulating recruitment thresholds and discharge patterns of a set of 100 and 200 motor units for the pinnate and bipinnate muscle, respectively (15° pinnation in both cases). Without subcutaneous layer and muscle fibers with CV 4 m/s, average CV estimates from the pinnate (bipinnate) muscle were 4.81 0.18 m/s (4.80 0.18 m/s) for bipolar, 4.71 0.19 m/s (4.71 0.12 m/s) for double differential, and 4.78 0.16 m/s (4.79 0.15 m/s) for Laplacian recordings. When subcutaneous layer was added (thickness 1 mm) in the same conditions, estimated CV values were 4.93 0.25 m/s (5.16 0.41 m/s), 4.70 0.21 m/s (4.83 0.33 m/s), and 4.89 0.21 m/s (4.99 0.39 m/s), for the three recording systems, respectively. The main factor biasing CV estimates was the propagation of action potentials in the two directions which influenced the recording due to the scatter of the projection of end-plate and tendon locations along the fiber direction, as a consequence of pinnation. The same problem arises in muscles with the line of innervation zone locations not perpendicular to fiber direction. These results indicate an important limitation in reliability of CV estimates from the interference EMG when the innervation zone and tendon locations are not distributed perpendicular to fiber direction.
Medical & Biological Engineering & Computing, 2004
The study of surface electromyographic (EMG) signals under dynamic contractions is becoming increasingly important. However, knowledge of the methodological issues that may affect such analysis is still limited. The aim of the study was to analyse the effect of fibre shortening on estimates of conduction velocity (CV) and mean power spectral frequency (MNF) from surface EMG signals. Single fibre action potentials were simulated, as detected by commonly used spatial filters, for different fibre lengths. No physiological modifications were included with changes in fibre length, and thus only geometrical artifacts related to fibre shortening were investigated. The simulation results showed that the dependence of CV and MNF on fibre shortening is affected by the fibre location, electrode position and the spatial filter applied. With shortening of up to 50% for a fibre of 50 mm semi-length, the variations in CV and MNF estimates with shortening in bipolar recordings were 0.5% (CV) and 0.7% (MNF) for superficial fibres, and 3.6% and 5.1% for deeper fibres. Using the longitudinal double differential filter, under the same conditions, the percent variation was 0% and 0.2%, and 24.7% and 15.8%, respectively. The main conclusions were, first, muscle fibre shortening can significantly affect estimates of CV and MNF, especially for short fibre lengths. However, for long (semi-length > 50 mm) and superficial fibres, this effect is limited for shortenings of up to 50% of the initial fibre length. Secondly, CV and MNF are almost equally affected by changes in muscle length; and, thirdly, sensitivity to fibre shortening depends on the spatial filter applied for signal detection.
Medical & Biological Engineering & Computing, 2004
Surface electromyographic (EMG) signal modelling is important for signal interpretation, testing of processing algorithms, detection system design and didactic purposes. Various surface EMG signal models have been proposed in the literature. This study focuses on the proposal of a method for modelling surface EMG signals, using either analytical or numerical descriptions of the volume conductor for space-conductor by numerical approaches, accurately describing the volume conductor geometry and the conductivity, as mainly done in the past, but also the conductivity tensor of the muscle tissue. For volume conductors that are space-invariant in the direction of source propagation, the surface potentials generated by any source can be computed by one-dimensional convolutions, once the volume conductor transfer function has been derived (analytically or numerically). Conversely, more complex volume conductors require a complete numerical approach. In a numerical approach, the conductivity tensor of the muscle tissue should be matched with the fibre orientation. In some cases (e.g. multi-pinnate muscles), accurate description of the conductivity tensor can be very complex. A method for relating the conductivity tensor of the muscle tissue, to be used in a numerical approach, to the curve describing the muscle fibres is presented and applied to investigate representatively a bi-pinnate muscle with rectilinear and curvilinear fibres. The study thus proposes an approach for surface EMG signal simulation in space invariant systems, as well as new models of the volume conductor using numerical methods.
Annals of Biomedical Engineering, 1993
From regular one-channel registrations of single muscle fiber action potential no measures on the recording configuration can be derived. When multichannel recordings are made, experimental parameters such as the distance between muscle fiber and needle electrode can be estimated. With the help of a volume conductor model, the single fiber activity at each of the electrodes can be predicted as a function of the recording conditions. Within known physical and physiological constraints such a model approach can be inverted (the inverse model) and used to estimate basic experimental conditions. From eight simultaneous single fiber action potential recordings we estimated (a) the distance between fiber and needle, (b) the axial position of the needle with respect to the muscle fiber, (c) a factor related to the muscle tissue anisotropy, and (d) a factor combining the muscle fiber diameter and the effective muscle tissue conductivity. With the help of a model describing the influence of the needle shaft it is made plausible that the needle inhomogeneity influences the results of the proposed procedure.
Journal of Electromyography and Kinesiology, 2005
Changes in muscle fibre length and surface electrode position with respect to the muscle fibres affect the amplitude and frequency characteristics of surface electromyography (SEMG) in different ways. Knowledge of changes in muscle fibre length would help towards a better interpretation of the signals. The possibility of estimating the length through SEMG during voluntary contractions was checked in this study. The fibresÕ semi-length was estimated from the product of the conduction velocity and conduction time during which the wave of excitation propagated from the end-plate region to the ends of the fibres. Short (10 s), moderate (30% of maximum voluntary contraction) isometric contractions were performed by 10 subjects at different elbow joint angles (80-140°in steps of 20°). Monopolar signals were detected non-invasively, using a two-dimensional electrode array. High spatial resolution EMG and a decomposition technique were utilised to extract single motor unit activities for triggered averaging and to estimate conduction velocity. A significant increase with joint angle was found in conduction time and estimated fibre semi-length. Changes in conduction velocity with joint angle were found to be not significant. The methodology described allows the relative changes in fibresÕ semi-length to be estimated from SEMG data.
Volume conductor models in surface electromyography: Computational techniques
Computers in Biology and Medicine, 2013
Simulation of surface electromyogram (EMG) provided support for the development and test of algorithms and for the interpretation of experimental results. This is the second part of a pair of papers: the first explains the methods for the development of structure based models of surface EMG; this one provides a summary of interesting applications for the interpretation of signals and test of algorithms. Simulations from some models already introduced in the literature are compared. Moreover, estimation of indexes from an innovative layered volume conductor including a subcutaneous tissue with variable thickness in space is studied.
IEEE Transactions on Biomedical Engineering, 2004
Surface EMG signal modeling has important applications in the interpretation of experimental EMG data. Most models of surface EMG generation considered volume conductors homogeneous in the direction of propagation of the action potentials. However, this may not be the case in practice due to local tissue in-homogeneities or to the fact that there may be groups of muscle fibers with different orientations. This study addresses the issue of analytically describing surface EMG signals generated by bi-pinnate muscles, i.e., muscles which have two groups of fibers with two orientations. The approach will also be adapted to the case of a muscle with fibers inclined in the depth direction. Such muscle anatomies are in-homogeneous in the direction of propagation of the action potentials with the consequence that the system can not be described as space invariant in the direction of source propagation. In these conditions, the potentials detected at the skin surface do not travel without shape changes. This determines numerical issues in the implementation of the model which are addressed in this work. The study provides the solution of the non-homogenous, anisotropic problem, proposes an implementation of the results in complete surface EMG generation models (including finite length fibers), and shows representative results of the application of the models proposed.
Volume Conduction in an Anatomically Based Surface EMG Model
IEEE Transactions on Biomedical Engineering, 2004
A finite-element model to simulate surface electromyography (EMG) in a realistic human upper arm is presented. The model is used to explore the effect of limb geometry on surface-detected muscle fiber action potentials. The model was based on magnetic resonance images of the subject's upper arm and includes both resistive and capacitive material properties. To validate the model geometry, experimental and simulated potentials were compared at different electrode sites during the application of a subthreshold sinusoidal current source to the skin surface. Of the material properties examined, the closest approximation to the experimental data yielded a mean root-mean-square (rms) error of the normalized surface potential of 18% or 27%, depending on the site of the applied source. Surface-detected action potentials simulated using the realistic volume conductor model and an idealized cylindrical model based on the same limb geometry were then compared. Variation in the simulated limb geometry had a considerable effect on action potential shape. However, the rate of decay of the action potential amplitude with increasing distance from the fiber was similar in both models. Inclusion of capacitive material properties resulted in temporal low-pass filtering of the surface action potentials. This effect was most pronounced in the end-effect components of action potentials detected at locations far from the active fiber. It is concluded that accurate modeling of the limb geometry, asymmetry, tissue capacitance and fiber curvature is important when the specific action potential shapes are of interest. However, if the objective is to examine more qualitative features of the surface EMG signal, then an idealized volume conductor model with appropriate tissue thicknesses provides a close approximation.