A Study of Myoelectric Signal Processing

This dissertation of various aspects of electromyogram (EMG: muscle electrical activity) signal processing is comprised of two projects in which I was the lead investigator and two team projects in which I participated. The first investigator-led project was a study of reconstructing continuous EMG discharge rates from neural impulses. Related methods for calculating neural firing rates in other contexts were adapted and applied to the intramuscular motor unit action potential train firing rate. Statistical results based on simulation and clinical data suggest that performances of spline-based methods are superior to conventional filter-based methods in the absence of decomposition error, but they unacceptably degrade in the presence of even the smallest decomposition errors present in real EMG data, which is typically around 3-5%. Optimal parameters for each method are found, and with normal decomposition error rates, ranks of these methods with their optimal parameters are given. Overall, Hanning filtering and Berger methods exhibit consistent and significant advantages over other methods. In the second investigator-led project, the technique of signal whitening was applied prior to motion classification of upper limb surface EMG signals previously collected from the forearm muscles of intact and amputee subjects. The motions classified consisted of 11 hand and wrist actions pertaining to prosthesis control. Theoretical models and experimental data showed that whitening increased EMG signal bandwidth by 65-75% and the coefficients of variation of temporal features computed from the EMG were reduced. As a result, a consistent classification accuracy improvement of 3-5% was observed for all subjects at small analysis durations (\u3c 100 ms). In the first team-based project, advanced modeling methods of the constant posture EMG-torque relationship about the elbow were studied: whitened and multi-channel EMG signals, training set duration, regularized model parameter estimation and nonlinear models. Combined, these methods reduced error to less than a quarter of standard techniques. In the second team-based project, a study related biceps-triceps surface EMG to elbow torque at seven joint angles during constant-posture contractions. Models accounting for co-contraction estimated that individual flexion muscle torques were much higher than models that did not account for co-contraction

Experimental Investigations of EMG-Torque Modeling for the Human Upper Limb

The electrical activity of skeletal muscle—the electromyogram (EMG)—is of value to many different application areas, including ergonomics, clinical biomechanics and prosthesis control. For many applications, the EMG is related to muscular tension, joint torque and/or applied forces. In these cases, a goal is for an EMG-torque model to emulate the natural relationship between the central nervous system (as evidenced in the surface EMG) and peripheral joints and muscles. This thesis work concentrated on experimental investigations of EMG-torque modeling. My contributions include: 1) continuing to evaluate the advantage of advanced EMG amplitude estimators, 2) studying system identification techniques (regularizing the least squares fit and increasing training data duration) to improve EMG-torque model performance, and 3) investigating the influence of joint angle on EMG-torque modeling. Results show that the advanced EMG amplitude estimator reduced the model error by 21%—71% compared to conventional estimators. Use of the regularized least squares fit with 52 seconds of training data reduced the model error by 20% compared to the least squares fit without regulation when using 26 seconds of training data. It is also demonstrated that the influence of joint angle can be modeled as a multiplicative factor in slowly force-varying and force-varying contractions at various, fixed angles. The performance of the models that account for the joint angle are not statistically different from a model that was trained at each angle separately and thus does not interpolate across angles. The EMG-torque models that account for joint angle and utilize advanced EMG amplitude estimation and system identification techniques achieved an error of 4.06±1.2% MVCF90 (i.e., error referenced to maximum voluntary contraction at 90° flexion), while models without using these advanced techniques and only accounting for a joint angle of 90° generated an error of 19.15±11.2% MVCF90. This thesis also summarizes other collaborative research contributions performed as part of this thesis. (1) EMG-force modeling at the finger tips was studied with the purpose of assessing the ability to determine two or more independent, continuous degrees of freedom of control from the muscles of the forearm [with WPI and Sherbrooke University]. (2) Investigation of EMG bandwidth requirements for whitening for real-time applications of EMG whitening techniques [with WPI colleagues]. (3) Investigation of the ability of surface EMG to estimate joint torque at future times [with WPI colleagues]. (4) Decomposition of needle EMG data was performed as part of a study to characterize motor unit behavior in patients with amyotrophic lateral sclerosis (ALS) [with Spaulding Rehabilitation Hospital, Boston, MA]

Improving the Performance of Dynamic Electromyogram-to-Force Models for the Hand-Wrist and Multiple Fingers

Relating surface electromyogram (EMG) activity to force/torque models is used in many areas including: prosthesis control systems, to regulate direction and speed of movement in reaching and matching tasks; clinical biomechanics, to assess muscle deficiency and effort levels; and ergonomics analysis, to assess risk of work-related injury such as back pain, fatigue and skill tests. This thesis work concentrated on improving the performance of dynamic EMG-to-force models for the hand-wrist and multiple fingers. My contributions include: 1) rapid calibration of dynamic hand-wrist EMG-force models using a minimum number of electrodes, 2) efficiently training two degree of freedom (DoF) hand-wrist EMG-force models, and 3) estimating individual and combined fingertip forces from forearm EMG during constant-pose, force-varying tasks. My calibration approach for hand-wrist EMG-force models optimized three main factors for 1-DoF and 2-DoF tasks: training duration (14, 22, 30, 38, 44, 52, 60, 68, 76 s), number of electrodes (2 through 16), and model forms (subject-specific, DoF-specific, universal). The results show that training duration can be reduced from historical 76 s to 40–60 s without statistically affecting the average error for both 1-DoF and 2-DoF tasks. Reducing the number of electrodes depended on the number of DoFs. One-DoF models can be reduced to 2 electrodes with average test error range of 8.3–9.2% maximum voluntary contraction (MVC), depending on the DoF (e.g., flexion-extension, radial-ulnar deviation, pronation-supination, open-close). Additionally, 2-DoF models can be reduced to 6 electrodes with average error of 7.17–9.21 %MVC. Subject-specific models had the lowest error for 1-DoF tasks while DoF-specific and universal were the lowest for 2-DoF tasks. In the EMG-finger project, we studied independent contraction of one, two, three or four fingers (thumb excluded), as well as contraction of four fingers in unison. Using regression, we found that a pseudo-inverse tolerance (ratio of largest to smallest singular value) of 0.01 was optimal. Lower values produced erratic models and higher values produced models with higher errors. EMG-force errors using one finger ranged from 2.5–3.8 %MVC, using the optimal pseudoinverse tolerance. With additional fingers (two, three or four), the average error ranged from 5–8 %MVC. When four fingers contracted in unison, the average error was 4.3 %MVC. Additionally, I participated in two team projects—EMG-force dynamic models about the elbow and relating forearm muscle EMG to finger force during slowly force varying contractions. This work is also described herein

Studies of the relationship between the surface electromyogram, joint torque and impedance

This compendium-format dissertation (i.e., comprised mostly of published and in-process articles) primarily reports on system identification methods that relate the surface electromyogram (EMG)—the electrical activity of skeletal muscles—to mechanical kinetics. The methods focus on activities of the elbow and hand-wrist. The relationship between the surface EMG and joint impedance was initially studied. My work provided a complete second-order EMG-based impedance characterization of stiffness, viscosity and inertia over a complete range of nominal torques, from a single perturbation trial with slowly varied torque. A single perturbation trial provides a more convenient method for impedance evaluation. The RMS errors of the EMG-based method were 20.01% for stiffness and 7.05% for viscosity, compared with the traditional mechanical measurement. Three projects studied the relationship between EMG and force/torque, a topic that has been studied for a number of years. Optimal models use whitened EMG amplitude, combining multiple EMG channels and a polynomial equation to describe this relationship. First, we used three techniques to improve current models at the elbow joint. Three more features were extracted from the EMG (waveform length, slope sign change rate and zero crossing rate), in addition to EMG amplitude. Each EMG channel was used separately, compared to previous studies which combined multiple channels from biceps and, separately, from triceps muscles. Finally, an exponential power law model was used. Each of these improvement techniques showed better performance (P\u3c0.05 and ~0.7 percent maximum voluntary contraction (%MVC) error reduction from a nominal error of 5.5%MVC) than the current “optimal” model. However, the combination of pairs of these techniques did not further improve results. Second, traditional prostheses only control 1 degree of freedom (DoF) at a time. My work provided evidence for the feasibility of controlling 2-DoF wrist movements simultaneously, with a minimum number of electrodes. Results suggested that as few as four conventional electrodes, optimally located about the forearm, could provide 2-DoF simultaneous, independent and proportional control with error ranging from 9.0–10.4 %MVC, which is similar to the 1-DoF approach (error from 8.8–9.8 %MVC) currently used for commercial prosthesis control. The third project was similar to the second, except that this project studied controlling a 1-DoF wrist with one hand DoF simultaneously. It also demonstrated good performance with the error ranging from 7.8-8.7 %MVC, compared with 1-DoF control. Additionally, I participated in two team projects—EMG decomposition and static wrist EMG to torque—which are described herein

Influence of Electromyogram (EMG) Amplitude Processing in EMG-Torque Estimation

A number of studies have investigated the relationship between surface electromyogram (EMG) and torque exerted about a joint. The standard deviation of the recorded EMG signal is defined as the EMG amplitude. The EMG amplitude estimation technique varies with the study from conventional type of processing (i.e. rectification followed by low pass filtering) to further addition of different noise rejection and signal-to-noise ratio improvement stages. Advanced EMG amplitude processors developed recently that incorporate signal whitening and multiple-channel combination have been shown to significantly improve amplitude estimation. The main contribution of this research is a comparison of the performance of EMG-torque estimators with and without these advanced EMG amplitude processors. The experimental data are taken from fifteen subjects that produced constant-posture, non-fatiguing, force-varying contractions about the elbow while torque and biceps/triceps EMG were recorded. Utilizing system identification techniques, EMG amplitude was related to torque through a zeros-only (finite impulse response, FIR) model. The incorporation of whitening and multiple-channel combination separately reduced EMG-torque errors and their combination provided a cumulative improvement. A 15th-order linear FIR model provided an average estimation error of 6% of maximum voluntary contraction (or 90% of variance accounted for) when EMG amplitudes were obtained using a four-channel, whitened processor. The equivalent single-channel, unwhitened (conventional) processor produced an average error of 8% of maximum voluntary contraction (variance accounted for of 68%). This study also describes the occurrence of spurious peaks in estimated torque when the torque model is created from data with a sampling rate well above the bandwidth of the torque. This problem is anticipated when the torque data are sampled at the same rate as the EMG data. The problem is resolved by decimating the EMG amplitude prior to relating it to joint torque, in this case to an effective sampling rate of 40.96 Hz

ESTIMATION OF MULTI-DIRECTIONAL ANKLE IMPEDANCE AS A FUNCTION OF LOWER EXTREMITY MUSCLE ACTIVATION

The purpose of this research is to investigate the relationship between the mechanical impedance of the human ankle and the corresponding lower extremity muscle activity. Three experimental studies were performed to measure the ankle impedance about multiple degrees of freedom (DOF), while the ankle was subjected to different loading conditions and different levels of muscle activity. The first study determined the non-loaded ankle impedance in the sagittal, frontal, and transverse anatomical planes while the ankle was suspended above the ground. The subjects actively co-contracted their agonist and antagonistic muscles to various levels, measured using electromyography (EMG). An Artificial Neural Network (ANN) was implemented to characterize the relationship between the EMG and non-loaded ankle impedance in 3-DOF. The next two studies determined the ankle impedance and muscle activity during standing, while the foot and ankle were subjected to ground perturbations in the sagittal and frontal planes. These studies investigate the performance of subject-dependent models, aggregated models, and the feasibility of a generic, subject-independent model to predict ankle impedance based on the muscle activity of any person. Several regression models, including Least Square, Support Vector Machine, Gaussian Process Regression, and ANN, and EMG feature extraction techniques were explored. The resulting subject-dependent and aggregated models were able to predict ankle impedance with reasonable accuracy. Furthermore, preliminary efforts toward a subject-independent model showed promising results for the design of an EMG-impedance model that can predict ankle impedance using new subjects. This work contributes to understanding the relationship between the lower extremity muscles and the mechanical impedance of the ankle in multiple DOF. Applications of this work could be used to improve user intent recognition for the control of active ankle-foot prostheses

Neuromechanical Tuning for Arm Motor Control

Movement is a fundamental behavior that allows us to interact with the external world. Its importance to human health is most evident when it becomes impaired due to disease or injury. Physical and occupational rehabilitation remains the most common treatment for these types of disorders. Although therapeutic interventions may improve motor function, residual deficits are common for many pathologies, such as stroke. The development of novel therapeutics is dependent upon a better understanding of the underlying mechanisms that govern movement. Movement of the human body adheres to the principles of classic Newtonian mechanics. However, due to the inherent complexity of the body and the highly variable repertoire of environmental contexts in which it operates, the musculoskeletal system presents a challenging control problem and the onus is on the central nervous system to reliably solve this problem. The neural motor system is comprised of numerous efferent and afferent pathways with a hierarchical organization which create a complex arrangement of feedforward and feedback circuits. However, the strategy that the neural motor system employs to reliably control these complex mechanics is still unknown. This dissertation will investigate the neural control of mechanics employing a “bottom-up” approach. It is organized into three research chapters with an additional introductory chapter and a chapter addressing final conclusions. Chapter 1 provides a brief description of the anatomical and physiological principles of the human motor system and the challenges and strategies that may be employed to control it. Chapter 2 describes a computational study where we developed a musculoskeletal model of the upper limb to investigate the complex mechanical interactions due to muscle geometry. Muscle lengths and moment arms contribute to force and torque generation, but the inherent redundancy of these actuators create a high-dimensional control problem. By characterizing these relationships, we found mechanical coupling of muscle lengths which the nervous system could exploit. Chapter 3 describes a study of muscle spindle contribution to muscle coactivation using a computational model of primary afferent activity. We investigated whether these afferents could contribute to motoneuron recruitment during voluntary reaching tasks in humans and found that afferent activity was orthogonal to that of muscle activity. Chapter 4 describes a study of the role of the descending corticospinal tract in the compensation of limb dynamics during arm reaching movements. We found evidence that corticospinal excitability is modulated in proportion to muscle activity and that the coefficients of proportionality vary in the course of these movements. Finally, further questions and future directions for this work are discussed in the Chapter 5

Development and Biomechanical Analysis toward a Mechanically Passive Wearable Shoulder Exoskeleton

Shoulder disability is a prevalent health issue associated with various orthopedic and neurological conditions, like rotator cuff tear and peripheral nerve injury. Many individuals with shoulder disability experience mild to moderate impairment and struggle with elevating the shoulder or holding the arm against gravity. To address this clinical need, I have focused my research on developing wearable passive exoskeletons that provide continuous at-home movement assistance. Through a combination of experiments and computational tools, I aim to optimize the design of these exoskeletons. In pursuit of this goal, I have designed, fabricated, and preliminarily evaluated a wearable, passive, cam-driven shoulder exoskeleton prototype. Notably, the exoskeleton features a modular spring-cam-wheel module, allowing customizable assistive force to compensate for different proportions of the shoulder elevation moment due to gravity. The results of my research demonstrated that this exoskeleton, providing modest one-fourth gravity moment compensation at the shoulder, can effectively reduce muscle activity, including deltoid and rotator cuff muscles. One crucial aspect of passive shoulder exoskeleton design is determining the optimal anti-gravity assistance level. I have addressed this challenge using computational tools and found that an assistance level within the range of 20-30% of the maximum gravity torque at the shoulder joint yields superior performance for specific shoulder functional tasks. When facing a new task dynamic, such as wearing a passive shoulder exoskeleton, the human neuro-musculoskeletal system adapts and modulates limb impedance at the end-limb (i.e., hand) to enhance task stability. I have presented development and validation of a realistic neuromusculoskeletal model of the upper limb that can predict stiffness modulation and motor adaptation in response to newly introduced environments and force fields. Future studies will explore the model\u27s applicability in predicting stiffness modulation for 3D movements in novel environments, such as passive assistive devices\u27 force fields

Estimation of Upper Limb Joint Angle Using Surface EMG Signal

In the development of robot-assisted rehabilitation systems for upper limb rehabilitation therapy, human electromyogram (EMG) is widely used due to its ability to detect the user intended motion. EMG is one kind of biological signal that can be recorded to evaluate the performance of skeletal muscles by means of a sensor electrode. Based on recorded EMG signals, user intended motion could be extracted via estimation of joint torque, force or angle. Therefore, this estimation becomes one of the most important factors to achieve accurate user intended motion. In this paper, an upper limb joint angle estimation methodology is proposed. A back propagation neural network (BPNN) is developed to estimate the shoulder and elbow joint angles from the recorded EMG signals. A Virtual Human Model (VHM) is also developed and integrated with BPNN to perform the simulation of the estimated angle. The relationships between sEMG signals and upper limb movements are observed in this paper. The effectiveness of our developments is evaluated with four healthy subjects and a VHM simulation. The results show that the methodology can be used in the estimation of joint angles based on EMG

Relating forearm muscle electrical activity to finger forces

The electromyogram (EMG) signal is desired to be used as a control signal for applications such as multifunction prostheses, wheelchair navigation, gait generation, grasping control, virtual keyboards, and gesture-based interfaces [25]. Several research studies have attempted to relate the electromyogram (EMG) activity of the forearm muscles to the mechanical activity of the wrist, hand and/or fingers [41], [42], [43]. A primary interest is for EMG control of powered upper-limb prostheses and rehabilitation orthotics. Existing commercial EMG-controlled devices are limited to rudimentary control capabilities of either discrete states (e.g. hand close/open), or one degree of freedom proportional control [4], [36]. Classification schemes for discriminating between hand/wrist functions and individual finger movements have demonstrated accuracy up to 95% [38], [39], [29]. These methods may provide for increased amputee function, though continuous control of movement is not generally achieved. This thesis considered proportional control via EMG-based estimation of finger forces with the goal of identifying whether multiple degrees of freedom of proportional control information are available from the surface EMG of the forearm. Electromyogram (EMG) activity from the extensor and flexor muscles of the forearm was sensed with bipolar surface electrodes and related to the force produced at the four fingertips during constant-posture, slowly force-varying contractions from 20 healthy subjects. The contractions ranged between 30% maximum voluntary contractions (MVC) extension and 30% MVC flexion. EMG amplitude sampling rate, least squares regularization, linear vs. nonlinear models and number of electrodes used in the system identification were studied. Results are supportive that multiple degrees of freedom of proportional control information are available from the surface EMG of the forearm, at least in healthy subjects. An EMG amplitude sampling frequency of 4.096 Hz was found to produce models which allowed for good EMG amplitude estimates. Least squares regularization with a pseudo-inverse tolerance of 0.055 resulted in significant improvement in modeling results, with an average error of 4.69% MVC-6.59% MVC (maximum voluntary contraction). Increasing polynomial order did not significantly improve modeling results. Results from smaller electrode arrays remained fairly good with as few as six electrodes, with the average %MVC error ranging from 5.13%-7.01% across the four fingers. This study also identified challenges in the current experimental study design and subsequent system identification when EMG-force modeling is performed with four fingers simultaneously. Methods to compensate for these issues have been proposed in this thesis