Open Access
A masked least-squares smoothing procedure for artifact reduction in scanning-EMG recordings
(Springer, 2018) Corera Orzanco, Íñigo; Eciolaza Ferrando, Adrián; Rubio Zamora, Oliver; Malanda Trigueros, Armando; Rodríguez Falces, Javier; Navallas Irujo, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
Scanning-EMG is an electrophysiological technique in which the electrical activity of the motor unit is recorded at multiple points along a corridor crossing the motor unit territory. Correct analysis of the scanning-EMG signal requires prior elimination of interference from nearby motor units. Although the traditional processing based on the median filtering is effective in removing such interference, it distorts the physiological waveform of the scanning-EMG signal. In this study, we describe a new scanning-EMG signal processing algorithm that preserves the physiological signal waveform while effectively removing interference from other motor units. To obtain a cleaned-up version of the scanning signal, the masked least-squares smoothing (MLSS) algorithm recalculates and replaces each sample value of the signal using a least-squares smoothing in the spatial dimension, taking into account the information of only those samples that are not contaminated with activity of other motor units. The performance of the new algorithm with simulated scanning-EMG signals is studied and compared with the performance of the median algorithm and tested with real scanning signals. Results show that the MLSS algorithm distorts the waveform of the scanning-EMG signal much less than the median algorithm (approximately 3.5 dB gain), being at the same time very effective at removing interference components.
Open Access
Evaluation of the electromyography test for the analysis of the aerobic-anaerobic transition in elite cyclists during incremental exercise
(MDPI, 2019) Latasa Zudaire, Iban Alexander; Córdova Martínez, Alfredo; Quintana Ortí, Gregorio; Lavilla Oiz, Ana; Navallas Irujo, Javier; Rodríguez Falces, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
(1) Background: The aim of this study was to investigate the validity and reliability of surface electromyography (EMG) for automatic detection of the aerobic and anaerobic thresholds during an incremental continuous cycling test using 1 min exercise periods in elite cyclists. (2) Methods: Sixteen well-trained cyclists completed an incremental exercise test (25 W/1 min) to exhaustion. Surface bipolar EMG signals were recorded from the vastus lateralis, vastus medialis, biceps femoris, and gluteus maximus, and the root mean square (RMS) were assessed. The multi-segment linear regression method was used to calculate the first and second EMG thresholds (EMG(T1) and EMG(T2)). During the test, gas exchange data were collected to determine the first and second ventilatory thresholds (VT1 and VT2). (3) Results: Two breakpoints (thresholds) were identified in the RMS EMG vs. time curve for all muscles in 75% of participants. The two breakpoints, EMG(T1) and EMG(T2) , were detected at around 70%-80% and 90%-95% of VO2MAX, respectively. No significant differences were found between the means of VT(1 )and EMGT(1) for the vastii and biceps femoris muscles (p > 0.05). There were no significant differences between means of EMG(T2) and VT2 (p > 0.05). (4) Conclusions: It is concluded that the multi-segment linear regression algorithm is a valid non-invasive method for analyzing the aerobic-anaerobic transition during incremental tests with 1 min stage durations.
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
Modelling fibrillation potentials: analysis of time parameters in the muscle intracellular action potential
(IEEE, 2007) Rodríguez Falces, Javier; Malanda Trigueros, Armando; Gila Useros, Luis; Rodríguez Carreño, Ignacio; Navallas Irujo, Javier; Ingeniería Eléctrica y Electrónica; Ingeniaritza Elektrikoa eta Elektronikoa
A single fibre 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 T spl is a time parameter related to the duration of the IAP waveform. Our work 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 T spl 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 T spl values corresponding to a given SFAP from its measured MSD.
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
Validation of the filling factor index to study the filling process of the sEMG signal in the quadriceps
(Elsevier, 2023) Rodríguez Falces, Javier; Malanda Trigueros, Armando; Mariscal Aguilar, Cristina; Niazi, Imran Khan; Navallas Irujo, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
Introduction: The EMG filling factor is an index to quantify the degree to which an EMG signal has been filled. Here, we tested the validity of such index to analyse the EMG filling process as contraction force was slowly increased. Methods: Surface EMG signals were recorded from the quadriceps muscles of healthy subjects as force was gradually increased from 0 to 40% MVC. The sEMG filling process was analyzed by measuring the EMG filling factor (calculated from the non-central moments of the rectified sEMG). Results: (1) As force was gradually increased, one or two prominent abrupt jumps in sEMG amplitude appeared between 0 and 10% of MVC force in all the vastus lateralis and medialis. (2) The jumps in amplitude were originated when a few large-amplitude MUPs, clearly standing out from previous activity, appeared in the sEMG signal. (3) Every time an abrupt jump in sEMG amplitude occurred, a new stage of sEMG filling was initiated. (4) The sEMG was almost completely filled at 2–12% MVC. (5) The filling factor decreased significantly upon the occurrence of an sEMG amplitude jump, and increased as additional MUPs were added to the sEMG signal. (6) The filling factor curve was highly repeatable across repetitions. Conclusions: It has been validated that the filling factor is a useful, reliable tool to analyse the sEMG filling process. As force was gradually increased in the vastus muscles, the sEMG filling process occurred in one or two stages due to the presence of abrupt jumps in sEMG amplitude.
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
Long-range and high-resolution traffic monitoring based on pulse-compression DAS and advanced vehicle tracking algorithm
(Optica Publishing Group, 2022) Corera Orzanco, Íñigo; Piñeiro Ben, Enrique; Navallas Irujo, Javier; Sagüés García, Mikel; Loayssa Lara, Alayn; Ingeniería Eléctrica, Electrónica y de Comunicación; Institute of Smart Cities - ISC; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
We demonstrate traffic monitoring over tens of kilometres of road using an enhanced distributed acoustic sensing system based on optical pulse compression and a novel transformed-domain-based processing scheme with enhanced vehicle detection and tracking capabilities.
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
A deep learning image analysis method for renal perfusion estimation in pseudo-continuous arterial spin labelling MRI
(Elsevier, 2023) Oyarzun Domeño, Anne; Cia Alonso, Izaskun; Echeverría Chasco, Rebeca; Fernández Seara, María A.; Martín Moreno, Paloma L.; Bastarrika, Gorka; Navallas Irujo, Javier; Villanueva Larre, Arantxa; Ingeniería Eléctrica, Electrónica y de Comunicación; Institute of Smart Cities - ISC; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Accurate segmentation of renal tissues is an essential step for renal perfusion estimation and postoperative assessment of the allograft. Images are usually manually labeled, which is tedious and prone to human error. We present an image analysis method for the automatic estimation of renal perfusion based on perfusion magnetic resonance imaging. Specifically, non-contrasted pseudo-continuous arterial spin labeling (PCASL) images are used for kidney transplant evaluation and perfusion estimation, as a biomarker of the status of the allograft. The proposed method uses machine/deep learning tools for the segmentation and classification of renal cortical and medullary tissues and automates the estimation of perfusion values. Data from 16 transplant patients has been used for the experiments. The automatic analysis of differentiated tissues within the kidney, such as cortex and medulla, is performed by employing the time-intensity-curves of non-contrasted T1-weighted MRI series. Specifically, using the Dice similarity coefficient as a figure of merit, results above 93%, 92% and 82% are obtained for whole kidney, cortex, and medulla, respectively. Besides, estimated cortical and medullary perfusion values are considered to be within the acceptable ranges within clinical practice.