Open Access
EMG probability density function: a new way to look at EMG signal filling from single motor unit potential to full interference pattern
(IEEE, 2023) Navallas Irujo, Javier; Eciolaza Ferrando, Adrián; Mariscal Aguilar, Cristina; Malanda Trigueros, Armando; Rodríguez Falces, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
An analytical derivation of the EMG signal's amplitude probability density function (EMG PDF) is presented and used to study how an EMG signal builds-up, or fills, as the degree of muscle contraction increases. The EMG PDF is found to change from a semi-degenerate distribution to a Laplacian-like distribution and finally to a Gaussian-like distribution. We present a measure, the EMG filling factor, to quantify the degree to which an EMG signal has been built-up. This factor is calculated from the ratio of two non-central moments of the rectified EMG signal. The curve of the EMG filling factor as a function of the mean rectified amplitude shows a progressive and mostly linear increase during early recruitment, and saturation is observed when the EMG signal distribution becomes approximately Gaussian. Having presented the analytical tools used to derive the EMG PDF, we demonstrate the usefulness of the EMG filling factor and curve in studies with both simulated signals and real signals obtained from the tibialis anterior muscle of 10 subjects. Both simulated and real EMG filling curves start within the 0.2 to 0.35 range and rapidly rise towards 0.5 (Laplacian) before stabilizing at around 0.637 (Gaussian). Filling curves for the real signals consistently followed this pattern (100% repeatability within trials in 100% of the subjects). The theory of EMG signal filling derived in this work provides (a) an analytically consistent derivation of the EMG PDF as a function of motor unit potentials and motor unit firing patterns; (b) an explanation of the change in the EMG PDF according to degree of muscle contraction; and (c) a way (the EMG filling factor) to quantify the degree to which an EMG signal has been built-up.
Open Access
A new muscle architecture model with non-uniform distribution of muscle fiber types
(World Academy of Science, Engineering and Technology, 2007) Navallas Irujo, Javier; Malanda Trigueros, Armando; Gila Useros, Luis; Rodríguez Falces, Javier; Rodríguez Carreño, Ignacio; Ingeniería Eléctrica y Electrónica; Ingeniaritza Elektrikoa eta Elektronikoa
According to previous studies, some muscles present a non-homogeneous spatial distribution of its muscle fiber types and motor unit types. However, available muscle models only deal with muscles with homogeneous distributions. In this paper, a new architecture muscle model is proposed to permit the construction of non-uniform distributions of muscle fibers within the muscle cross section. The idea behind is the use of a motor unit placement algorithm that controls the spatial overlapping of the motor unit territories of each motor unit type. Results show the capabilities of the new algorithm to reproduce arbitrary muscle fiber type distributions.
Open Access
Masked least-squares averaging in processing of scanning-EMG recordings with multiple-discharges
(Springer, 2020) Corera Orzanco, Íñigo; Malanda Trigueros, Armando; Rodríguez Falces, Javier; Navallas Irujo, Javier; Ingeniería Eléctrica y Electrónica; Ingeniaritza Elektrikoa eta Elektronikoa
Removing artifacts from nearby motor units is one of the main objectives when processing scanning-EMG recordings. Methods such as median filtering or masked least-squares smoothing (MLSS) can be used to eliminate artifacts in recordings with just one discharge of the motor unit potential (MUP) at each location. However, more effective artifact removal can be achieved if several discharges per position are recorded. In this case, processing usually involves averaging the discharges available at each position and then applying a median filter in the spatial dimension. The main drawback of this approach is that the median filter tends to distort the signal waveform. In this paper, we present a new algorithm that operates on multiple discharges simultaneously and in the spatial dimension. We refer to this algorithm as the multi masked least-squares smoothing (MMLSS) algorithm: an extension of the MLSS algorithm for the case of multiple discharges. The algorithm is tested using simulated scanning-EMG signals in different recording conditions, i.e., at different levels of muscle contraction and for different numbers of discharges per position. Results demonstrate that the algorithm eliminates artifacts more effectively than any previously available method and does so without distorting the waveform of the signal.
Open Access
Understanding EMG PDF changes with motor unit potential amplitudes, firing rates, and noise level through EMG filling curve analysis
(IEEE, 2024-08-30) Navallas Irujo, Javier; Mariscal Aguilar, Cristina; Malanda Trigueros, Armando; Rodríguez Falces, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoa eta Telekomunikazio Ingeniaritza
EMG filling curve characterizes the EMG filling process and EMG probability density function (PDF) shape change for the entire force range of a muscle.We aim to understand the relation between the physiological and recording variables, and the resulting EMG filling curves. We thereby present an analytical and simulation study to explain how the filling curve patterns relate to specific changes in the motor unit potential (MUP) waveforms and motor unit (MU) firing rates, the two main factors affecting the EMG PDF, but also to recording conditions in terms of noise level. We compare the analytical results with simulated cases verifying a perfect agreement with the analytical model. Finally, we present a set of real EMG filling curves with distinct patterns to explain the information about MUP amplitudes, MU firing rates, and noise level that these patterns provide in the light of the analytical study. Our findings reflect that the filling factor increases when firing rate increases or when newly recruited motor unit have potentials of smaller or equal amplitude than the former ones. On the other hand, the filling factor decreases when newly recruited potentials are larger in amplitude than the previous potentials. Filling curves are shown to be consistent under changes of the MUP waveform, and stretched under MUP amplitude scaling. Our findings also show how additive noise affects the filling curve and can even impede to obtain reliable information from the EMG PDF statistics.
Open Access
M-wave changes caused by brief voluntary and stimulated isometric contractions
(Springer, 2023) Rodríguez Falces, Javier; Malanda Trigueros, Armando; Navallas Irujo, Javier; Place, Nicolas; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Introduction Under isometric conditions, the increase in muscle force is accompanied by a reduction in the fbers’ length. The efects of muscle shortening on the compound muscle action potential (M wave) have so far been investigated only by computer simulation. This study was undertaken to assess experimentally the M-wave changes caused by brief voluntary and stimulated isometric contractions. Methods Two diferent methods of inducing muscle shortening under isometric condition were adopted: (1) applying a brief (1 s) tetanic contraction and (2) performing brief voluntary contractions of diferent intensities. In both methods, supramaximal stimulation was applied to the brachial plexus and femoral nerves to evoke M waves. In the frst method, electrical stimulation (20 Hz) was delivered with the muscle at rest, whereas in the second, stimulation was applied while participants performed 5-s stepwise isometric contractions at 10, 20, 30, 40, 50, 60, 70, and 100% MVC. The amplitude and duration of the frst and second M-wave phases were computed. Results The main fndings were: (1) on application of tetanic stimulation, the amplitude of the M-wave frst phase decreased (~10%, P<0.05), that of the second phase increased (~50%, P<0.05), and the M-wave duration decreased (~20%, P<0.05) across the frst fve M waves of the tetanic train and then plateaued for the subsequent responses; (2) when superimposing a single electrical stimulus on muscle contractions of increasing forces, the amplitude of the M-wave frst phase decreased (~20%, P<0.05), that of the second phase increased (~30%, P<0.05), and M-wave duration decreased (~30%, P<0.05) as force was raised from 0 to 60–70% MVC force. Conclusions The present results will help to identify the adjustments in the M-wave profle caused by muscle shortening and also contribute to diferentiate these adjustments from those caused by muscle fatigue and/or changes in Na+–K+ pump activity.
Open Access
Motor unit action potential duration, II: a new automatic measurement method based on the wavelet transform
(Lippincott, Williams & Wilkins, 2007) Rodríguez Carreño, Ignacio; Gila Useros, Luis; Malanda Trigueros, Armando; García Gurtubay, Ignacio; Mallor Giménez, Fermín; Gómez Elvira, Sagrario; Rodríguez Falces, Javier; Navallas Irujo, Javier; Ingeniería Eléctrica y Electrónica; Estadística e Investigación Operativa; Ingeniaritza Elektrikoa eta Elektronikoa; Estatistika eta Ikerketa Operatiboa
To present and evaluate a new algorithm, based on the wavelet transform, for the automatic measurement of motor unit action potential (MUAP) duration. A total of 240 MUAPs were studied. The waveform of each MUAP was wavelet-transformed, and the start and end points were estimated by regarding the maxima and minima points in a particular scale of the wavelet transform. The results of the new method were compared with the gold standard of duration marker positions obtained by manual measurement. The new method was also compared with a conventional algorithm, which we had found to be best in a previous comparative study. To evaluate the new method against manual measurements, the dispersion of automatic and manual duration markers were analyzed in a set of 19 repeatedly recorded MUAPs. The differences between the new algorithm’s marker positions and the gold standard of duration marker positions were smaller than those observed with the conventional method. The dispersion of the new algorithm’s marker positions was slightly less than that of the manual one. Our new method for automatic measurement of MUAP duration is more accurate than other available algorithms and more consistent than manual measurements.
Open Access
The probability density function of the surface electromyogram and its dependence on contraction force in the vastus lateralis
(BMC, 2024-10-26) Rodríguez Falces, Javier; Malanda Trigueros, Armando; Mariscal Aguilar, Cristina; Recalde Villamayor, Silvia; Navallas Irujo, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoa eta Telekomunikazio Ingeniaritza; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Introduction: the probability density function (PDF) of the surface electromyogram (sEMG) depends on contraction force. This dependence, however, has so far been investigated by having the subject generate force at a few fixed percentages of MVC. Here, we examined how the shape of the sEMG PDF changes with contraction force when this force was gradually increased from zero. Methods: voluntary surface EMG signals were recorded from the vastus lateralis of healthy subjects as force was increased in a continuous manner vs. in a step-wise fashion. The sEMG filling process was examined by measuring the EMG filling factor, computed from the non-central moments of the rectified sEMG signal. Results: in 84% of the subjects, as contraction force increased from 0 to 10% MVC, the sEMG PDF shape oscillated back and forth between the semi-degenerate and the Gaussian distribution; the PDF–force relation varied greatly among subjects for forces between 0 and ~ 10% MVC, but this variability was largely reduced for forces above 10% MVC; the pooled analysis showed that, as contraction force gradually increased, the sEMG PDF evolved rapidly from the semi-degenerate towards the Laplacian distribution from 0 to 5% MVC, and then more slowly from the Laplacian towards the Gaussian distribution for higher forces. Conclusions: the study demonstrated that the dependence of the sEMG PDF shape on contraction force can only be reliably assessed by gradually increasing force from zero, and not by performing a few constant-force contractions. The study also showed that the PDF–force relation differed greatly among individuals for contraction forces below 10% MVC, but this variability was largely reduced when force increased above 10% MVC.
Open Access
EMG modeling
(InTechOpen, 2012) Rodríguez Falces, Javier; Navallas Irujo, Javier; Malanda Trigueros, Armando; Ingeniería Eléctrica y Electrónica; Ingeniaritza Elektrikoa eta Elektronikoa
The aim of this chapter is to describe the approaches used for modelling electromyographic (EMG) signals as well as the principles of electrical conduction within the muscle. Sections are organized into a progressive, step-by-step EMG modeling of structures of increasing complexity. First, the basis of the electrical conduction that allows for the propagation of the EMG signals within the muscle is presented. Second, the models used for describing the electrical activity generated by a single fibre described. The third section is devoted to modeling the organization of the motor unit and the generation of motor unit potentials. Based on models of the architectural organization of motor units and their activation and firing mechanisms, the last section focuses on modeling the electrical activity of a complete muscle as recorded at the surface.
Open Access
Métodos de procesamiento y análisis de señales electromiográficas
(Gobierno de Navarra, 2009) Gila Useros, Luis; Malanda Trigueros, Armando; Rodríguez Carreño, Ignacio; Rodríguez Falces, Javier; Navallas Irujo, Javier; Ingeniería Eléctrica y Electrónica; Ingeniaritza Elektrikoa eta Elektronikoa
La electromiografía clínica es una metodología de registro y análisis de la actividad bioeléctrica del músculo esquelético orientada al diagnóstico de las enfermedades neuromusculares. Las posibilidades de aplicación y el rendimiento diagnóstico de la electromiografía han evolucionado paralelamente al conocimiento de las propiedades de la energía eléctrica y al desarrollo de la tecnología eléctrica y electrónica. A mediados del siglo XX se introdujo el primer equipo comercial de electromiografía para uso médico basado en circuitos electrónicos analógicos. El desarrollo posterior de la tecnología digital ha permitido disponer de sistemas controlados por microprocesadores cada vez más fiables y potentes para captar, representar, almacenar, analizar y clasificar las señales mioeléctricas. Es esperable que el avance de las nuevas tecnologías de la información y la comunicación pueda conducir en un futuro próximo a la aplicación de desarrollos de inteligencia artificial que faciliten la clasificación automática de señales así como sistemas expertos de apoyo al diagnóstico electromiográfico.
Open Access
Sliding window averaging in normal and pathological motor unit action potential trains
(Elsevier, 2018) Malanda Trigueros, Armando; Navallas Irujo, Javier; Rodríguez Falces, Javier; Porta Cuéllar, Sonia; Fernández Martínez, Miguel; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
Objective: To evaluate the performance of a recently proposed motor unit action potential (MUAP) averaging method based on a sliding window, and compare it with relevant published methods in normal and pathological muscles. Methods: Three versions of the method (with different window lengths) were compared to three relevant published methods in terms of signal analysis-based merit figures and MUAP waveform parameters used in the clinical practice. 218 MUAP trains recorded from normal, myopathic, subacute neurogenic and chronic neurogenic muscles were analysed. Percentage scores of the cases in which the methods obtained the best performance or a performance not significantly worse than the best were computed. Results: For signal processing figures of merit, the three versions of the new method performed better (with scores of 100, 86.6 and 66.7%) than the other three methods (66.7, 25 and 0%, respectively). In terms of MUAP waveform parameters, the new method also performed better (100, 95.8 and 91.7%) than the other methods (83.3, 37.5 and 25%). Conclusions: For the types of normal and pathological muscle studied, the sliding window approach extracted more accurate and reliable MUAP curves than other existing methods. Significance: The new method can be of service in quantitative EMG.