Stimulation artifact in surface EMG signal: effect of the stimulation waveform, detection system, and current amplitude using hybrid stimulation technique

Tutorial. Surface electromyogram (sEMG) amplitude estimation: Best practices

This tutorial intends to provide insight, instructions and "best practices" for those who are novices-including clinicians, engineers and non-engineers-in extracting electromyogram (EMG) amplitude from the bipolar surface EMG (sEMG) signal of voluntary contractions. A brief discussion of sEMG amplitude extraction from high density sEMG (HDsEMG) arrays and feature extraction from electrically elicited contractions is also provided. This tutorial attempts to present its main concepts in a straightforward manner that is accessible to novices in the field not possessing a wide range of technical background (if any) in this area. Surface EMG amplitude, also referred to as the sEMG envelope [often implemented as root mean square (RMS) sEMG or average rectified value (ARV) sEMG], quantifies the voltage variation of the sEMG signal and is grossly related to the overall neural excitation of the muscle and to peripheral parameters. The tutorial briefly reviews the physiological origin of the voluntary sEMG signal and sEMG recording, including electrode configurations, sEMG signal transduction, electronic conditioning and conversion by an analog-to-digital converter. These topics have been covered in greater detail in prior tutorials in this series. In depth descriptions of state-of-the-art methods for computing sEMG amplitude are then provided, including guidance on signal pre-conditioning, absolute value vs. square-law detection, selection of appropriate sEMG amplitude smoothing filters and attenuation of measurement noise. The tutorial provides a detailed list of best practices for sEMG amplitude estimation.

An Adaptive Spatial Filtering Method for Multi-Channel EMG Artifact Removal During Functional Electrical Stimulation With Time-Variant Parameters

Removing the stimulation artifacts evoked by the functional electrical stimulation (FES) in electromyogram (EMG) signals is a challenge. Previous researches on stimulation artifact removal have focused on FES modulation with time-constant parameters, which has limitations when there are time-variant parameters. Therefore, considering the synchronism of muscle activation induced by FES and the asynchronism of muscle activation induced by proprioceptive nerves, we proposed a novel adaptive spatial filtering method called G-S-G. It entails fusing the Gram-Schmidt orthogonalization (G-S) and Grubbs criterion (G) algorithms to remove the FES-evoked stimulation artifacts in multi-channel EMG signals. To verify this method, we constructed a series of simulation data by fusing the FES signal with time-variant parameters and the voluntary EMG (vEMG) signal, and applied the G-S-G method to remove any FES artifacts from the simulation data. After that, we calculated the root mean square (RMS) value for both preprocessed simulation data and the vEMG data, and then compared them. The simulation results showed that the G-S-G method was robust and effective at removing FES artifacts in simulated EMG signals, and the correlation coefficient between the preprocessed EMG data and the recorded vEMG data yielded a good performance, up to 0.87. Furthermore, we applied the proposed method to the experimental EMG data with FES-evoked stimulation artifact, and also achieved good performance with both the time-constant and time-variant parameters. This study provides a new and accessible approach to resolving the problem of removing FES-evoked stimulation artifacts.

Toward a wearable monitor of local muscle fatigue during electrical muscle stimulation using tissue Doppler imaging

Electrical muscle stimulation (EMS) is widely used in rehabilitation and athletic training to generate involuntary muscle contractions. However, EMS leads to rapid muscle fatigue, limiting the force a muscle can produce during prolonged use. Currently available methods to monitor localized muscle fatigue and recovery are generally not compatible with EMS. The purpose of this study was to examine whether Doppler ultrasound imaging can assess changes in stimulated muscle twitches that are related to muscle fatigue from electrical stimulation. We stimulated five isometric muscle twitches in the medial and lateral gastrocnemius of 13 healthy subjects before and after a fatiguing EMS protocol. Tissue Doppler imaging of the medial gastrocnemius recorded muscle tissue velocities during each twitch. Features of the average muscle tissue velocity waveforms changed immediately after the fatiguing stimulation protocol (peak velocity: -38%, p = .022; time-to-zero velocity: +8%, p = .050). As the fatigued muscle recovered, the features of the average tissue velocity waveforms showed a return towards their baseline values similar to that of the normalized ankle torque. We also found that features of the average tissue velocity waveform could significantly predict the ankle twitch torque for each participant (R2 = 0.255-0.849, p < .001). Our results provide evidence that Doppler ultrasound imaging can detect changes in muscle tissue during isometric muscle twitch that are related to muscle fatigue, fatigue recovery, and the generated joint torque. Tissue Doppler imaging may be a feasible method to monitor localized muscle fatigue during EMS in a wearable device.

Ultrasound Echogenicity as an Indicator of Muscle Fatigue during Functional Electrical Stimulation

Functional electrical stimulation (FES) is a potential neurorehabilitative intervention to enable functional movements in persons with neurological conditions that cause mobility impairments. However, the quick onset of muscle fatigue during FES is a significant challenge for sustaining the desired functional movements for more extended periods. Therefore, a considerable interest still exists in the development of sensing techniques that reliably measure FES-induced muscle fatigue. This study proposes to use ultrasound (US) imaging-derived echogenicity signal as an indicator of FES-induced muscle fatigue. We hypothesized that the US-derived echogenicity signal is sensitive to FES-induced muscle fatigue under isometric and dynamic muscle contraction conditions. Eight non-disabled participants participated in the experiments, where FES electrodes were applied on their tibialis anterior (TA) muscles. During a fatigue protocol under either isometric and dynamic ankle dorsiflexion conditions, we synchronously collected the isometric dorsiflexion torque or dynamic dorsiflexion angle on the ankle joint, US echogenicity signals from TA muscle, and the applied stimulation intensity. The experimental results showed an exponential reduction in the US echogenicity relative change (ERC) as the fatigue progressed under the isometric (R2=0.891±0.081) and dynamic (R2=0.858±0.065) conditions. The experimental results also implied a strong linear relationship between US ERC and TA muscle fatigue benchmark (dorsiflexion torque or angle amplitude), with R2 values of 0.840±0.054 and 0.794±0.065 under isometric and dynamic conditions, respectively. The findings in this study indicate that the US echogenicity signal is a computationally efficient signal that strongly represents FES-induced muscle fatigue. Its potential real-time implementation to detect fatigue can facilitate an FES closed-loop controller design that considers the FES-induced muscle fatigue.

Ultrasound Echogenicity-based Assessment of Muscle Fatigue During Functional Electrical Stimulation

The rapid onset of muscle fatigue during functional electrical stimulation (FES) is a major challenge when attempting to perform long-term periodic tasks such as walking. Surface electromyography (sEMG) is frequently used to detect muscle fatigue for both volitional and FES-evoked muscle contraction. However, sEMG contamination from both FES stimulation artifacts and residual M-wave signals requires sophisticated processing to get clean signals and evaluate the muscle fatigue level. The objective of this paper is to investigate the feasibility of computationally efficient ultrasound (US) echogenicity as a candidate indicator of FES-induced muscle fatigue. We conducted isometric and dynamic ankle dorsiflexion experiments with electrically stimulated tibialis anterior (TA) muscle on three human participants. During a fatigue protocol, we synchronously recorded isometric dorsiflexion force, dynamic dorsiflexion angle, US images, and stimulation intensity. The temporal US echogenicity from US images was calculated based on a gray-scaled analysis to assess the decrease in dorsiflexion force or motion range due to FES-induced TA muscle fatigue. The results showed a monotonic reduction in US echogenicity change along with the fatigue progression for both isometric (R² =0.870±0.026) and dynamic (R² =0.803±0.048) ankle dorsiflexion. These results implied a strong linear relationship between US echogenicity and TA muscle fatigue level. The findings indicate that US echogenicity may be a promising computationally efficient indicator for assessing FES-induced muscle fatigue and may aid in the design of muscle-in-the-loop FES controllers that consider the onset of muscle fatigue.

A Closed-Loop Neuromodulation Chipset With 2-Level Classification Achieving 1.5-Vpp CM Interference Tolerance, 35-dB Stimulation Artifact Rejection in 0.5ms and 97.8%-Sensitivity Seizure Detection

This work presents an 8-channel closed-loop neuromodulation chipset with 2-level seizure classification. The power-consuming fine classifier is only enabled when the coarse classifier in the frontend chip judges the patient's status as "suspected seizure". This scheme can reduce the overall power consumption extensively since seizure usually occurs with very low possibility. In the capacitive-coupled instrument amplifier (CCIA) of the front-end IC, a feedback based common-mode (CM) cancellation circuit is proposed to suppress large-scale CM interferences and the stimulation artifacts are suppressed by a mixed-signal loop with fast response. An auto-zero based pre- charge path is adopted to boost the input impedance, while the electrode DC offset is canceled by a DC servo loop with very-large and accurate time constant. The 2.32-mm² front-end chip and 3.51-mm² DSP chip implemented in 0.18 μm CMOS are applied in a deep-brain stimulation (DBS) neuromodulator. Measurement results show that the CCIA can suppress 1.5-Vpp CM interference, and achieve an accurate high-pass corner frequency as low as 0.1 Hz and an input impedance greater than 2.2 GΩ. The overall classifier achieves 97.8% sensitivity and consumes only 1.16-μW average power for the CHB-MIT database test. The chipset has been verified by in vivo measurement, showing that the stimulation artifact can be suppressed by 35 dB within 0.5 ms.

A Novel Technique to Reject Artifact Components for Surface EMG Signals Recorded During Walking With Transcutaneous Spinal Cord Stimulation: A Pilot Study

Transcutaneous spinal cord electrical stimulation (tSCS) is an emerging technology that targets to restore functionally integrated neuromuscular control of gait. The purpose of this study was to demonstrate a novel filtering method, Artifact Component Specific Rejection (ACSR), for removing artifacts induced by tSCS from surface electromyogram (sEMG) data for investigation of muscle response during walking when applying spinal stimulation. Both simulated and real tSCS contaminated sEMG data from six stroke survivors were processed using ACSR and notch filtering, respectively. The performance of the filters was evaluated with data collected in various conditions (e.g., simulated artifacts contaminating sEMG in multiple degrees, various tSCS intensities in five lower-limb muscles of six participants). In the simulation test, after applying the ACSR filter, the contaminated-signal was well matched with the original signal, showing a high correlation (r = 0.959) and low amplitude difference (normalized root means square error = 0.266) between them. In the real tSCS contaminated data, the ACSR filter showed superior performance on reducing the artifacts (96% decrease) over the notch filter (25% decrease). These results indicate that ACSR filtering is capable of eliminating artifacts from sEMG collected during tSCS application, improving the precision of quantitative analysis of muscle activity.

A compact system for simultaneous stimulation and recording for closed-loop myoelectric control

Background

Despite important advancements in control and mechatronics of myoelectric prostheses, the communication between the user and his/her bionic limb is still unidirectional, as these systems do not provide somatosensory feedback. Electrotactile stimulation is an attractive technology to close the control loop since it allows flexible modulation of multiple parameters and compact interface design via multi-pad electrodes. However, the stimulation interferes with the recording of myoelectric signals and this can be detrimental to control.

Methods

We present a novel compact solution for simultaneous recording and stimulation through dynamic blanking of stimulation artefacts. To test the system, a feedback coding scheme communicating wrist rotation and hand aperture was developed specifically to stress the myoelectric control while still providing meaningful information to the subjects. Ten subjects participated in an experiment, where the quality of closed-loop myoelectric control was assessed by controlling a cursor in a two degrees of freedom target-reaching task. The benchmark performance with visual feedback was compared to that achieved by combining visual feedback and electrotactile stimulation as well as by using electrotactile feedback only.

Results

There was no significant difference in performance between visual and combined feedback condition with regards to successfully reached targets, time to reach a target, path efficiency and the number of overshoots. Therefore, the quality of myoelectric control was preserved in spite of the stimulation. As expected, the tactile condition was significantly poorer in completion rate (100/4% and 78/25% for combined and tactile condition, respectively) and time to reach a target (9/2 s and 13/4 s for combined and tactile condition, respectively). However, the performance in the tactile condition was still good, with no significant difference in path efficiency (38/8%) and the number of overshoots (0.5/0.4 overshoots), indicating that the stimulation was meaningful for the subjects and useful for closed-loop control.

Conclusions

Overall, the results demonstrated that the developed system can provide robust closed-loop control using electrotactile stimulation. The system supports different encoding schemes and allows placing the recording and stimulation electrodes next to each other. This is an important step towards an integrated solution where the developed unit will be embedded into a prosthetic socket.

A hybrid method for real-time stimulation artefact removal during functional electrical stimulation with time-variant parameters

Objective. In this study, a hybrid method combining hardware and software architecture is proposed to remove stimulation artefacts (SAs) and extract the volitional surface electromyography (sEMG) in real time during functional electrical stimulations (FES) with time-variant parameters.Approach. First, an sEMG detection front-end (DFE) combining fast recovery, detector and stimulator isolation and blanking is developed and is capable of preventing DFE saturation with a blanking time of 7.6 ms. The fragment between the present stimulus and previous stimulus is set as an SA fragment. Second, an SA database is established to provide six high-similarity templates with the current SA fragment. The SA fragment will be de-artefacted by a 6th-order Gram-Schmidt (GS) algorithm, a template-subtracting method, using the provided templates, and this database-based GS algorithm is called DBGS. The provided templates are previously collected SA fragments with the same or a similar evoking FES intensity to that of the current SA fragment, and the lengths of the templates are longer than that of the current SA fragment. After denoising, the sEMG will be extracted, and the current SA fragment will be added to the SA database. The prototype system based on DBGS was tested on eight able-bodied volunteers and three individuals with stroke to verify its capacity for stimulation removal and sEMG extraction.Results.The average stimulus artefact attenuation factor, SA index and correlation coefficient between clean sEMG and extracted sEMG for 6th-order DBGS were 12.77 ± 0.85 dB, 1.82 ± 0.37 dB and 0.84 ± 0.33 dB, respectively, which were significantly higher than those for empirical mode decomposition combined with notch filters, pulse-triggered GS algorithm, 1st-order and 3rd-order DBGS. The sEMG-torque correlation coefficients were 0.78 ± 0.05 and 0.48 ± 0.11 for able-bodied volunteers and individuals with stroke, respectively.Significance.The proposed hybrid method can extract sEMG during dynamic FES in real time.

Use of Surface EMG in Clinical Rehabilitation of Individuals With SCI: Barriers and Future Considerations

Surface electromyography (sEMG) is a widely used technology in rehabilitation research and provides quantifiable information on the myoelectric output of a muscle. In this perspective, we discuss the barriers which have restricted the wide-spread use of sEMG in clinical rehabilitation of individuals with spinal cord injury (SCI). One of the major obstacles is integrating the time-consuming aspects of sEMG in the already demanding schedule of physical therapists, occupational therapists, and other clinicians. From the clinicians' perspective, the lack of confidence to use sEMG technology is also apparent due to their limited exposure to the sEMG technology and possibly limited mathematical foundation through educational and professional curricula. Several technical challenges include the limited technology-transfer of ever-evolving knowledge from sEMG research into the off-the-shelf EMG systems, lack of demand from the clinicians for systems with advanced features, lack of user-friendly intuitive interfaces, and the need for a multidisciplinary approach for accurate handling and interpretation of data. We also discuss the challenges in the application and interpretation of sEMG that are specific to SCI, which are characterized by non-standardized approaches in recording and interpretation of EMGs due to the physiological and structural state of the spinal cord. Addressing the current barriers will require a collaborative, interdisciplinary, and unified approach. The most relevant steps could include enhancing user-experience for students pursuing clinical education through revised curricula through sEMG-based case studies/projects, hands-on involvement in the research, and formation of a common platform for clinicians and technicians for self-education and knowledge share.

An Anti Stimulation Artifacts and M-waves Surface Electromyography Detector with a Short Blanking Time. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020;2020:4126-4129. [PMID: 33018906 DOI: 10.1109/embc44109.2020.9176373]

A surface electromyography (sEMG) detector, that not only removes stimulation artifacts entirely but also increases the recording time, has been developed in this paper. The sEMG detector consists of an sEMG detection circuit and a stimulation isolator. The sEMG detection circuit employs a stimulus isolate switch (SIS), a blanking (BLK) and non-linear feed-back (NFB) circuit to remove the artifacts and to increase the recording time. In the SIS, the connection between stimulator and stimulation electrodes, along with the stimulation electrodes and the ground are controlled by an opto-isolator, and the connection of instrument amplifier and the recording electrodes are controlled by CMOS-based switches. The mode switches of the BLK and the NFB circuit also employs CMOS-based switches. By an accurate timing adjustment, the voluntary EMG can be recorded during electrical stimulation. Two 6 able-bodied experiments have been performed to test the three anti-artifact sEMG detector: BLK, BLK&SIS, BLK&SIS&NFB. The results indicate that the BLK&SIS&NFB proposed in this work effectively removes stimulus artifacts and M-waves, and has a longer recording time compared with BLK and BLK&SIS circuits.

Restoring Finger-Specific Sensory Feedback for Transradial Amputees via Non-Invasive Evoked Tactile Sensation

Objective: This study assessed the feasibility to restore finger-specific sensory feedback in transradial amputees with electrical stimulation of evoked tactile sensation (ETS). Methods: Here we investigated primary somatosensory cortical (SI) responses of ETS using Magnetoencephalography. Results: SI activations revealed a causal correlation with peripheral stimulation of projected finger regions on the stump skin. Peak latency was accountable to neural transmission from periphery to SI. Peak intensity of SI response was proportional to the strength of peripheral stimulation, manifesting a direct neural pathway from skin receptors to SI neurons. Active regions in SI at the amputated side were consistent to the finger/hand map of homunculus, forming a mirror imaging to that of the contralateral hand. With sensory feedback, amputees can recognize a pressure at prosthetic fingers as that at the homonymous lost fingers. Conclusions: Results confirmed that the direct neural pathway from periphery to SI allows effective communication of finger-specific sensory information to these amputees.

A Method for Suppressing Electrical Stimulation Artifacts from Electromyography

When surface electromyography (EMG) signal is used in a real-time functional electrical stimulation (FES) system for feedback control, the artifact from electrical stimulation is a key challenge for EMG signal processing. To address this challenge, this study proposes a novel method to suppress stimulation artifacts in the EMG-driven closed-loop FES system. The proposed method is inspired by an experimental study that compares artifacts generated by electrical stimulations with different current intensities. It is found that (1) spikes of stimulation artifacts are susceptible to the current intensity and (2) tailing components are similar under different current intensities. Based on these observations, the proposed method combines the blanking and template subtracting strategies for suppressing stimulation artifact. The length of blanking window for suppressing the stimulation spike is adaptively determined by a spike detection algorithm and the first-order derivative analysis of signal. An autoregressive model is used to estimate the tailing part of stimulation artifact, which is an adaptive template for subtracting the artifact. The proposed method is evaluated on both semi-synthetic and experimental datasets. Verified on the semi-synthetic dataset, the proposed method achieves better performance than the classic blanking method. Validated on the experimental dataset, the proposed method substantially decreases the power of stimulation artifact in the EMG. These results indicate that the proposed method can effectively suppress the stimulation artifact while retains the useful EMG signal for an EMG-driven FES system.

Quantitative Assessment of Changes in Muscle Contractility Due to Fatigue During NMES: An Ultrasound Imaging Approach

OBJECTIVE

This paper investigates an ultrasound imaging-based non-invasive methodology to quantitatively assess changes in muscle contractility due to the fatigue induced by neuromuscular electrical stimulation (NMES).

METHODS

Knee extension experiments on human participants were conducted to record synchronized isometric knee force data and ultrasound images of the electrically stimulated quadriceps muscle. The data were first collected in a pre-fatigue stage and then in a post-fatigue stage. Ultrasound images were processed using a contraction rate adaptive speckle tracking algorithm. A two-dimensional strain measure field was constructed based on the muscle displacement tracking results to quantify muscle contractility.

RESULTS

Analysis of the strain images showed that, between the pre-fatigue and post-fatigue stages, there was a reduction in the strain peaks, a change in the strain peak distribution, and a decrease in an area occupied by the large positive strain.

CONCLUSION

The results indicate changes in muscle contractility due to the NMES-induced muscle fatigue.

SIGNIFICANCE

Ultrasound imaging with the proposed methodology is a promising tool for a direct NMES-induced fatigue assessment and facilitates new strategies to alleviate the effects of the NMES-induced fatigue.

Assessment of Skeletal Muscle Contractile Properties by Radial Displacement: The Case for Tensiomyography

Skeletal muscle operates as a near-constant volume system; as such muscle shortening during contraction is transversely linked to radial deformation. Therefore, to assess contractile properties of skeletal muscle, radial displacement can be evoked and measured. Mechanomyography measures muscle radial displacement and during the last 20 years, tensiomyography has become the most commonly used and widely reported technique among the various methodologies of mechanomyography. Tensiomyography has been demonstrated to reliably measure peak radial displacement during evoked muscle twitch, as well as muscle twitch speed. A number of parameters can be extracted from the tensiomyography displacement/time curve and the most commonly used and reliable appear to be peak radial displacement and contraction time. The latter has been described as a valid non-invasive means of characterising skeletal muscle, based on fibre-type composition. Over recent years, applications of tensiomyography measurement within sport and exercise have appeared, with applications relating to injury, recovery and performance. Within the present review, we evaluate the perceived strengths and weaknesses of tensiomyography with regard to its efficacy within applied sports medicine settings. We also highlight future tensiomyography areas that require further investigation. Therefore, the purpose of this review is to critically examine the existing evidence surrounding tensiomyography as a tool within the field of sports medicine.

A survey on the feasibility of surface EMG in facial pacing

A survey on the feasibility of surface electromyography (EMG) measurements in facial pacing is presented. Pacing for unilateral facial paralysis consists of the measurement of activity from the healthy side of the face and functional electrical stimulation to reanimate the paralyzed one. The goal of this study is to evaluate the feasibility of surface EMG as a measurement method to detect muscle activations and to determine their intensities. Prior work is discussed, and results from experiments where 12 participants carried out a set of facial movements are presented. EMG was registered from zygomaticus major (smile), orbicularis oris (lip pucker), orbicularis oculi (eye blink), corrugator supercilii (frown), and masseter (chew). Most important facial functions that are limited due to the paralysis are blinking, smiling, and puckering. With majority of the participants, crosstalk between the measured EMG channels was found to be acceptably small to be able to pace smiling and puckering based on detecting their contraction intensities from the healthy side. However, pacing blinking based on orbicularis oculi EMG measurement does not seem possible due to crosstalk from other muscles, but the electro-oculographic (EOG) signals that couple to the same measurement channel could help to detect eye blinks and trigger stimuli. Futhermore, masseter greatly disturbs EMG measurement of most facial muscles, which needs to be addressed in the pacing system to avoid falsely interpreting its activity as the activity of another muscle.

Surgical Model for Analysis of Signal Transfer Through Biologic and Synthetic Materials

BACKGROUND

High-fidelity volitional control of bioengineered prosthetic limbs with multiple degrees of freedom requires the implantation of multiple recording interfaces to detect independent control signals. However, interface utilization is complicated by interfering electrophysiological signals originating from surrounding muscles and nerves, leading to equivocal signal detection. We developed and validated a surgical model to characterize signal propagation through various biomaterials to identify insulating substrates for use in implantable interfaces. The identification of these insulating materials will facilitate the acquisition of noncontaminated prosthetic control signals, thus improving manipulation of advanced prosthetic limbs.

METHODS

Using a rat hindlimb model, 4 groups (n = 8/group) were tested. A medial gastrocnemius muscle flap was elevated, leaving the neurovascular pedicle intact. The flap was rotated into a chamber and secured to a silicone base. A stainless steel electrode was affixed to the surface of a muscle and encircled by 1-layer small intestinal submucosa (SIS), 4-layer SIS, silicone elastomer, or nothing (uninsulated). A superimposing electrode was attached, and an external silicone layer was wrapped around the construct and sutured in place. Electromyographic studies were then performed.

RESULTS

This model was found to correspond with expected signal isolation characteristics of the nonconductive silicone group, electrically inert single and multilayer SIS group, and the uninsulated group. Signal isolation of compound muscle action potential amplitude at stimulation threshold was significantly greater using silicone (51.4%) compared with the 1-layer SIS (-6.8%), 4-layer SIS (-3.3% ), or uninsulated groups (1.2%) (P = <0.001). Isolation of the maximum compound muscle action potential peak-to-peak amplitude was also greater with silicone (56.7%) versus the 1-layer SIS (1.5%), 4-layer SIS (1.1%), or uninsulated groups (-0.7%) (P = <0.001).

CONCLUSIONS

This study demonstrates and validates a novel surgical model to characterize in vivo signal propagation and subsequently identify insulating materials for use in implantable interface systems currently in development. Improved signal isolation through the utilization of these materials stands to greatly improve control fidelity of neuroprosthetic limbs.

High Density EMG investigation of H-reflex distribution over the soleus muscle

The spatial distribution of H-reflexes over soleus muscle was investigated through High-Density EMG in five healthy subjects. The posterior tibial nerve was stimulated with a staircase current envelope with 1mA steps. The regions where the incremental responses (incremental H-reflexes) occurred were identified for each stimulation step with a validated segmentation algorithm. The average centroid of the segmented areas was located over the Achilles tendon, 5 cm below the myo-tendinous junction of the medial gastrocnemius. The average dimension of these regions corresponded to 28% of the surface covered by the grid of electrodes. The amplitude of H-reflexes recorded in the segmented areas was higher than the average amplitude computed over the entire detection system as well as the H-reflex recorded by the electrode positioned according to SENIAM guidelines. These preliminary results suggest that: i) H-reflex detected from a specific soleus region unlikely reflects the whole muscle volume and ii) H-reflexes with greatest amplitude can be recorded over the Achilles tendon.

Suppression of stimulus artifact contaminating electrically evoked electromyography

BACKGROUND

Electrical stimulation of muscle or nerve is a very useful technique for understanding of muscle activity and its pathological changes for both diagnostic and therapeutic purposes. During electrical stimulation of a muscle, the recorded M wave is often contaminated by a stimulus artifact. The stimulus artifact must be removed for appropriate analysis and interpretation of M waves.

OBJECTIVES

The objective of this study was to develop a novel software based method to remove stimulus artifacts contaminating or superimposing with electrically evoked surface electromyography (EMG) or M wave signals.

METHODS

The multiple stage method uses a series of signal processing techniques, including highlighting and detection of stimulus artifacts using Savitzky-Golay filtering, estimation of the artifact contaminated region with Otsu thresholding, and reconstruction of such region using signal interpolation and smoothing. The developed method was tested using M wave signals recorded from biceps brachii muscles by a linear surface electrode array. To evaluate the performance, a series of semi-synthetic signals were constructed from clean M wave and stimulus artifact recordings with different degrees of overlap between them.

RESULTS

The effectiveness of the developed method was quantified by a significant increase in correlation coefficient and a significant decrease in root mean square error between the clean M wave and the reconstructed M wave, compared with those between the clean M wave and the originally contaminated signal. The validity of the developed method was also demonstrated when tested on each channel's M wave recording using a linear electrode array.

CONCLUSIONS

The developed method can suppress stimulus artifacts contaminating M wave recordings.

Closed-Loop Control of Myoelectric Prostheses With Electrotactile Feedback: Influence of Stimulation Artifact and Blanking

Electrocutaneous stimulation is a promising approach to provide sensory feedback to amputees, and thus close the loop in upper limb prosthetic systems. However, the stimulation introduces artifacts in the recorded electromyographic (EMG) signals, which may be detrimental for the control of myoelectric prostheses. In this study, artifact blanking with three data segmentation approaches was investigated as a simple method to restore the performance of pattern recognition in prosthesis control (eight motions) when EMG signals are corrupted by stimulation artifacts. The methods were tested over a range of stimulation conditions and using four feature sets, comprising both time and frequency domain features. The results demonstrated that when stimulation artifacts were present, the classification performance improved with blanking in all tested conditions. In some cases, the classification performance with blanking was at the level of the benchmark (artifact-free data). The greatest pulse duration and frequency that allowed a full performance recovery were 400 μs and 150 Hz, respectively. These results show that artifact blanking can be used as a practical solution to eliminate the negative influence of the stimulation artifact on EMG pattern classification in a broad range of conditions, thus allowing to close the loop in myoelectric prostheses using electrotactile feedback.

Muscle motor point identification is essential for optimizing neuromuscular electrical stimulation use

Transcutaneous neuromuscular electrical stimulation applied in clinical settings is currently characterized by a wide heterogeneity of stimulation protocols and modalities. Practitioners usually refer to anatomic charts (often provided with the user manuals of commercially available stimulators) for electrode positioning, which may lead to inconsistent outcomes, poor tolerance by the patients, and adverse reactions. Recent evidence has highlighted the crucial importance of stimulating over the muscle motor points to improve the effectiveness of neuromuscular electrical stimulation. Nevertheless, the correct electrophysiological definition of muscle motor point and its practical significance are not always fully comprehended by therapists and researchers in the field. The commentary describes a straightforward and quick electrophysiological procedure for muscle motor point identification. It consists in muscle surface mapping by using a stimulation pen-electrode and it is aimed at identifying the skin area above the muscle where the motor threshold is the lowest for a given electrical input, that is the skin area most responsive to electrical stimulation. After the motor point mapping procedure, a proper placement of the stimulation electrode(s) allows neuromuscular electrical stimulation to maximize the evoked tension, while minimizing the dose of the injected current and the level of discomfort. If routinely applied, we expect this procedure to improve both stimulation effectiveness and patient adherence to the treatment.

The aims of this clinical commentary are to present an optimized procedure for the application of neuromuscular electrical stimulation and to highlight the clinical implications related to its use.

A blink restoration system with contralateral EMG triggered stimulation and real-time artifact blanking

Patients suffering from facial paralysis are on the hazard of disfigurement and loss of vision due to loss of blink function. Functional-electrical stimulation (FES) is one possible way of restoring blink and other functions in these patients. A blink restoration system for uni-lateral facial paralyzed patients is described in this paper. The system achieves restoration of synchronized blink through processing the myoelectric signal of orbicularis oculi at the normal side in real-time as the trigger to stimulate the paralyzed eyelid. Design issues are discussed, including EMG processing, stimulating strategies and real-time artifact blanking. Two artifact removal approaches based on sample and hold and digital filtering technique are proposed and implemented. Finally, the whole system has been verified on rabbit models.

Finding physiological responses in vestibular evoked potentials

Vestibular prostheses are regarded as a promising tool to restore lost sensation in patients with vestibular disorders. These prostheses often electrically stimulate the vestibular nerve and stimulation efficacy is evaluated by measuring the vestibulo-ocular reflex (VOR). However, eye movement recording as intuitive metric of vestibular functionality is difficult to obtain outside the laboratory environment, and hence not available as an error signal in a closed-loop prosthesis. Recently we investigated vestibular evoked potentials (VEPs) by stimulating and recording in the same semicircular canal of a guinea pig. Here we studied the correlation between VOR and one region of VEP. We further analyzed a second portion of VEP, where vestibular nerve activity should occur using rectified bin integration (RBI). To this end, stimulation artifact was significantly reduced by hardware and software approaches. We found a high VEP-VOR correlation (R-squared=0.86), suggesting that VEP could substitute VOR as metric of vestibular function. Differences between below and above vestibular threshold stimulation were seen for the second portion of VEP. Further investigations are required to determine the specific parts of VEP that accurately represents vestibular function(s).

Mechanisms of cramp contractions: peripheral or central generation?

We analysed the cramp threshold (i.e. the minimum frequency of electrical stimulation capable of inducing a cramp) and the behaviour of individual motor units during cramps electrically elicited in the absence (intact condition) and presence (blocked condition) of a peripheral nerve block in eight healthy subjects. The cramp threshold was significantly greater in the blocked than in the intact condition (18 ± 3 Hz vs. 13 ± 3 Hz; P = 0.01). Cramp duration and peak EMG amplitude in the intact condition (55.6 ± 19.2 s and 47.5 ± 24.8 μV, respectively) were significantly greater compared to the blocked condition (2.6 ± 1.3 s and 13.9 ± 8.8 μV; P < 0.01). All motor units identified in the blocked condition (n = 38) had a shorter interval of activity and a greater discharge rate compared to the intact condition (n = 37) (respectively, 1.1 ± 1.0 s vs. 29.5 ± 21.8 s, P < 0.0001; 25.7 ± 11.6 pulses s(-1) vs. 20.0 ± 5.9 pulses s(-1); P < 0.05). The motor unit activity detected during the blocked condition corresponded to spontaneous discharges of the motor nerves, while in the intact condition the motor unit discharge patterns presented the typical characteristics of motor neuron discharges. These results indicate a spinal involvement at the origin of cramps and during their development.

Muscle fiber conduction slowing and decreased levels of circulating muscle proteins after short-term dexamethasone administration in healthy subjects

CONTEXT

Glucocorticoids are known to decrease protein synthesis and impair membrane excitability of muscle fibers. However, their short-term effects on muscle structure and function of healthy subjects remain poorly understood.

OBJECTIVE

Our objective was to investigate whether steroid administration could decrease the circulating levels of muscle proteins and modify myoelectric indexes of sarcolemmal excitability and fatigability.

DESIGN

We conducted a single-blind, placebo-controlled study in 20 men randomized to receive dexamethasone (8 mg/d) or placebo for 1 wk. Blood sampling, force measurements for knee extensors and elbow flexors, and electrophysiological tests for biceps brachii, vastus lateralis and medialis, and tibialis anterior muscles were performed before and after the intervention.

RESULTS

Dexamethasone administration improved force by 6.0 +/- 6.0% (P = 0.01) for elbow flexors and by 8.5 +/- 5.5% (P < 0.01) for knee extensors, decreased levels of creatine kinase by 50.5 +/- 30.0% (P < 0.01) and myoglobin by 41.8 +/- 17.5% (P < 0.01), and impaired sarcolemmal excitability, as shown by the decline of muscle fiber conduction velocity for the four muscles (range from -6 to -10.5%, P < 0.05). Moreover, significant reductions of the myoelectric manifestations of fatigue were observed for the four muscles; the decrease in the rate of change of the mean frequency of the electromyographic power spectrum ranged from -22.6 to -43.9% (P < 0.05). In contrast, no significant changes were observed in muscle excitability and fatigability in subjects who received the placebo.

CONCLUSIONS

The demonstration that glucocorticoid-induced muscle impairments can be unraveled by means of blood sampling and noninvasive electrophysiological tests has clinical implications for the early identification of subclinical or preclinical forms of myopathy in treated patients.

Investigation of motor unit recruitment during stimulated contractions of tibialis anterior muscle

This work investigated motor unit (MU) recruitment during transcutaneous electrical stimulation (TES) of the tibialis anterior (TA) muscle, using experimental and simulated data. Surface electromyogram (EMG) and torque were measured during electrically-elicited contractions at different current intensities, on eight healthy subjects. EMG detected during stimulation (M-wave) was simulated selecting the elicited MUs on the basis of: (a) the simulated current density distribution in the territory of each MU and (b) the excitation threshold characteristic of the MU. Exerted force was simulated by adding the contribution of each of the elicited MUs. The effects of different fat layer thickness (between 2 and 8mm), different distributions of excitation thresholds (random excitation threshold, higher threshold for larger MUs or smaller MUs), and different MU distributions within the muscle (random distribution, larger MU deeper in the muscle, smaller MU deeper) on EMG variables and torque were tested. Increase of the current intensity led to a first rapid increase of experimental M-wave amplitude, followed by a plateau. Further increases of the stimulation current determined an increase of the exerted force, without relevant changes of the M-wave. Similar results were obtained in simulations. Rate of change of conduction velocity (CV) and leading coefficient of the second order polynomial interpolating the force vs. stimulation level curve were estimated as a function of increasing current amplitudes. Experimental data showed an increase of estimated CV with increasing levels of the stimulation current (for all subjects) and a positive leading coefficient of force vs. stimulation current curve (for five of eight subjects). Simulations matched the experimental results only when larger MUs were preferably located deeper in the TA muscle (in line with a histochemical study). Marginal effect of MU excitation thresholds was observed, suggesting that MUs closer to the stimulation electrode are recruited first during TES regardless of their excitability.

Time and frequency domain analysis of surface myoelectric signals during electrically-elicited cramps

OBJECTIVES

To examine if different frequencies of electrical stimulation trigger different sized cramps in the abductor hallucis muscle and to analyze their surface electromyographic (EMG) behaviour in both time and frequency domains.

METHODS

Fifteen subjects were studied. Stimulation trains of 150 pulses were applied to the muscle motor point. Frequency was increased (starting from 4pps with 2-pps steps) until a cramp developed. Current intensity was 30% higher than that eliciting maximal M-waves. After the first cramp ("threshold cramp"), a 30-minute rest was provided before a second cramp ("above-threshold cramp") was elicited with a frequency increased by 50% with respect to that eliciting the first cramp.

RESULTS

We found greater EMG amplitude and a compression of the power spectrum for above-threshold cramps with respect to threshold cramps. M-wave changes (ranging between small decreases of M-wave amplitude to complete M-wave disappearance) occurred and progressively increased throughout stimulation trains. Significant positive correlations were found between estimates of EMG amplitude during cramps and estimated reductions of M-wave amplitude.

CONCLUSIONS

Varying frequencies of electrical stimulation triggered different sized cramps. Moreover, decreases in M-wave amplitude were observed during both threshold and above-threshold stimulations. The choice of the stimulation frequency has relevance for optimizing electrical stimulation protocols for the study of muscle cramps in both healthy and pathological subjects.

Pulse charge and not waveform affects M-wave properties during progressive motor unit activation

The aim of this study was to investigate changes in experimentally recorded M-waves with progressive motor unit (MU) activation induced by transcutaneous electrical stimulation with different pulse waveforms. In 10 subjects, surface electromyographic signals were detected with a linear electrode array during electrically elicited contractions of the biceps brachii muscle. Three different monophasic waveforms of 304-micros duration were applied to the stimulation electrode on the main muscle motor point: triangular, square, and sinusoidal. For each waveform, increasing stimulation current intensities were applied in 10 s (frequency: 20 Hz). It was found that: (a) the degree of MU activation, as indicated by M-wave average rectified value, was a function of the injected charge and not of the stimulation waveform, and (b) MUs tended to be recruited in order of increasing conduction velocity with increasing charge of transcutaneous stimulation. Moreover, the subjects reported lower discomfort during the contractions elicited by the triangular waveform with respect to the others. Since subject tolerance to the stimulation protocol must be considered as important as MU recruitment in determining the effectiveness of neuromuscular electrical stimulation (NMES), we suggest that both charge and waveform of the stimulation pulses should be considered relevant parameters for optimizing NMES protocols.

Distribution of electrical stimulation current in a planar multilayer anisotropic tissue

This study analytically addresses the problem of neuromuscular electrical stimulation for a planar, multilayer, anisotropic model of a physiological tissue (referred to as volume conductor). Both conductivity and permittivity of the volume conductor are considered, including dispersive properties. The analytical solution is obtained in the 2-D Fourier transform domain, transforming in the planes parallel to the volume conductor surface. The model is efficient in terms of computational cost, as the solution is analytical (only numerical Fourier inversion is needed). It provides the current distribution in a physiological tissue induced by an electrical current delivered at the skin surface. Three representative examples of application of the model are considered. 1) The simulation of stimulation artefact during transcutaneous electrical stimulation and EMG detection. Only the effect of the volume conductor is considered, neglecting the other sources of artefact (such as the capacitive coupling between the stimulating and recording electrodes). 2) The simulation of the electrical current distribution within the muscle and the low-pass filter effect of the volume conductor on sinusoidal stimulation currents with different stimulation frequencies. 3) The estimation of the amplitude modulated current distribution within the muscle for interferential stimulation. The model is devoted to the simulation of neuromuscular stimulation, but the same method could be applied in other fields in which the estimation of the electrical current distribution in a medium induced by the injection of a current from the boundary of the medium is of interest.

Surface EMG and spatial resolution analysis with estimation of electromyographic descriptors

Spatial filtering has become a common way to improve the resolution of surface electromyographic signals (SEMG) when used in connection with electrode arrays. The goal of this study is to observe the behavior of S-EMG amplitude and spectral descriptors when signals are submitted to a longitudinal quadruple differentiating spatial filter. Signals were acquired at 20% and 60% of the maximum voluntary contraction using a linear array of eight surface electrodes in order to understand the impact of the filtering technique in the S-EMG variables during fatiguing and non-fatiguing contractions. The final results show that the filtering procedure yields better selectivity, suggesting that single motor units can be better observed if spatial filters and measurement configurations with smaller pick-up areas are used. During fatiguing contractions, however, further analysis is needed.

Reliability of a novel neurostimulation method to study involuntary muscle phenomena

Experimental methods involving painful electrical stimulation of a peripheral nerve showed the existence of a minimum stimulation frequency capable of inducing cramp, termed "threshold frequency" (TF). Our aim was to test an alternative method to induce fasciculations and cramps electrically. Two daily sessions of electrical stimulation of the abductor hallucis muscle were performed in 19 volunteers on 3 days: stimulation trains of 150 monophasic square pulses (duration 152 micros) of increasing frequency (current intensity 30% higher than maximal; frequency of the first trial, 4 pps; recovery between trials, 1 min) were delivered to the main muscle motor point until a cramp developed. Once a cramp was induced the protocol was repeated after 30 min. To verify by electromyography that cramp occurred, a surface electrode array was placed between the motor point and the distal tendon. Ambient and skin temperature were kept constant in all sessions. Fasciculations and cramps were elicited in all subjects. We observed the following median (interquartile range) values of TF: day 1 (session 1), 13 (6) pps; day 1 (session 2), 16 (4) pps; day 2 (session 1), 16 (6) pps; day 2 (session 2), 18 (6) pps; day 3 (session 1), 17 (4) pps; day 3 (session 2), 18 (8) pps. TF intersession intraclass correlation coefficients were 0.82, 0.92, and 0.90 for days 1, 2, and 3, respectively. TF interday intraclass correlation coefficient was 0.85. The absence of pain due to the stimulation and the demonstration of TF reliability support the use of our method for the study of involuntary muscle phenomena.

Estimation of M-Wave Scale Factor During Sustained Contractions at High Stimulation Rate

In this paper, we propose a time-domain index to assess M-wave widening during high-frequency stimulation, as an objective parameter for quantifying muscle fatigue. At high stimulation frequencies, signal truncation, due to the delivery of the electrical stimulus before the M-wave generated by the previous stimulus extinguishes, biases the spectral frequency variables usually computed to estimate M-wave widening. Thus, we propose an estimator of the scale factor between two truncated M-waves. The estimator is derived from the Scale Transforms of the two signals, with an efficient implementation that avoids limits of resolution. The method was tested on both simulated and experimental signals. The simulations showed that the proposed technique is significantly less affected by signal truncation than previous approaches. The experimental recordings were collected from 11 subjects at stimulation frequencies of 20, 40, and 60 Hz. The scale factor estimation assessed M-wave widening in the three conditions, differentiating between the different rates of change of signal widening. The method proved to be significantly superior to M-wave spectral analysis. The technique can be applied to investigate myoelectric manifestations of muscle fatigue at stimulation rates that could not be analyzed in the past and, thus, opens new perspectives in the evaluation of electrical stimulation for training and rehabilitation protocols.

M-wave properties during progressive motor unit activation by transcutaneous stimulation

The aim of this study was to interpret changes in experimentally recorded M waves with progressive motor unit (MU) activation based on simulation of the surface electromyogram. Activation order during transcutaneous electrical stimulation was analyzed by investigating M-wave average rectified value, spectral properties, and conduction velocity (CV) during electrically elicited contractions. M-waves were detected from the biceps brachii muscle of 10 healthy male subjects by a linear adhesive array of eight electrodes. Electrical stimulation was delivered to the motor point at either constant current intensity (40, 60, 80, and 100% of the supramaximal stimulation current) or with linearly increasing current. A model of surface electromyogram generation that varied activation order based on MU size and location was used to interpret the experimental results. From the experimental and model analysis, it was found that 1) MUs tended to be activated from low to high CV and from the superficial to the deep muscle layers with increasing transcutaneous electrical stimulation of the biceps brachii muscle, and 2) characteristic spectral frequencies of the M-wave were affected by many factors other than average CV (such as the activation order by MU location or the spread of the MU innervation zones and CVs), thus decreasing with a concomitant increase in CV during progressive MU activation.