| Literature DB >> 34229307 |
Aaron Fleming1,2, Nicole Stafford3, Stephanie Huang1,2, Xiaogang Hu1,2, Daniel P Ferris4, He Helen Huang1,2.
Abstract
Objective.Advanced robotic lower limb prostheses are mainly controlled autonomously. Although the existing control can assist cyclic movements during locomotion of amputee users, the function of these modern devices is still limited due to the lack of neuromuscular control (i.e. control based on human efferent neural signals from the central nervous system to peripheral muscles for movement production). Neuromuscular control signals can be recorded from muscles, called electromyographic (EMG) or myoelectric signals. In fact, using EMG signals for robotic lower limb prostheses control has been an emerging research topic in the field for the past decade to address novel prosthesis functionality and adaptability to different environments and task contexts. The objective of this paper is to review robotic lower limb Prosthesis control via EMG signals recorded from residual muscles in individuals with lower limb amputations.Approach.We performed a literature review on surgical techniques for enhanced EMG interfaces, EMG sensors, decoding algorithms, and control paradigms for robotic lower limb prostheses.Main results.This review highlights the promise of EMG control for enabling new functionalities in robotic lower limb prostheses, as well as the existing challenges, knowledge gaps, and opportunities on this research topic from human motor control and clinical practice perspectives.Significance.This review may guide the future collaborations among researchers in neuromechanics, neural engineering, assistive technologies, and amputee clinics in order to build and translate true bionic lower limbs to individuals with lower limb amputations for improved motor function. Creative Commons Attribution license.Entities:
Keywords: EMG; gait and balance; human motor control; neural–machine interface; robotic lower limb protheses
Mesh:
Year: 2021 PMID: 34229307 PMCID: PMC8694273 DOI: 10.1088/1741-2552/ac1176
Source DB: PubMed Journal: J Neural Eng ISSN: 1741-2552 Impact factor: 5.379
Supervisory EMG control review summary.
| Author[ | Muscles recorded[ | Classifier type[ | EMG features | Windowing/prediction time | Training | Control type | Amputation type | Activity |
|---|---|---|---|---|---|---|---|---|
| Au | TA, GAS, SOL | Neural network | Pre-processed EMG | Predictions at 0.5 Hz | Back propagation method with pre-processed EMG as input | None; offline classification only | Transtibial | Trajectory tracking sitting |
| Jin | ADDL, TFL, RF, VL, VM, BF, SM, ST | Distribution of different muscle features | Mean absolute value, waveform length, mean square value, zero tangent number, median power frequency | Entire gait cycle | 3–5 times along walkway for each condition | None | Transfemoral | Level-ground walking at varying speeds, stairs, and ramps |
| Huang | GME , GMA , SAR, RF, VL, VM, GRA, BFL, BFS, SEM, ADM | LDA and neural network | 3rd order autoregression coefficients, mean absolute value, zero crossings, waveform length, number of slope sign changes and root mean square | Range from 50 to 150 ms EMG window, 200 ms phase window | Ten walking trials of all activities with own prosthesis | None; offline classification only | Transfemoral | Level-ground walking, stepping over and obstacle, stairs, ipsilateral and contralateral turning, and standing |
| Au | GASL, GASM, TA | Neural network | Variance of EMG signals | 100 ms EMG window | Matching imagined ankle position to virtual ankle | Impedance control | Transtibial (bilateral) | Level-ground walking, stair descents |
| Ha | Quadriceps, hamstring | LDA, QDA | Not provided | Not provided | 100 s of knee flexion/extension visualizations | Online impedance control (prosthesis next to participant) | Transfemoral | Virtual tracking tasks (sitting) |
| Huang | GME, GMA, SAR, RF, VL, VM, ST, GRA, BFL, BFS, ADM | SVM, LDA | Zero crossings, signal direction change, mean absolute value, waveform length | 150 ms sliding EMG window prediction | 15 times each activity | None; online classification | Transfemoral | Walking, stair/ramp ascent and descent, stepping over obstacles (and transitions) |
| Simon | SM, SAR, TFL, ADM, GRA, RF, VL, BFL | Two LDA classifiers | Mean absolute value, zero crossings, waveform length, slope sign change number, five-count majority vote | 250 ms overlapping | Ten sit-to-stand trials | Impedance control | Transfemoral | Sitting down, walking and transition between them, repositioning prosthesis while sitting |
| Hargrove | Two reinnervated muscle segments, proximal BF, RF, VL, VM, SAR, GRA, ADM, TFL | DBN | Not provided | Not provided | Offline, 20 reps of locomotor circuit | Impedance control | Transfemoral | Level-ground walking, ramps, stairs and outside stairs |
| Young | ST, BF, TFL, VM, SAR, ADM, GRA | Two LDA models | Mean absolute value, waveform length, zero crossings, slope sign change, 1st two coefficients of 3rd order autoregressive model | 300 ms before heel contact and toe off | Offline, 20 reps of two locomotor circuits, train classifier on three subjects test 1 as novel user. Then add 5–10 min of level-ground walking for novel user to train classifier. | None; offline classification | Transfemoral | Level-ground walking, ramps, stairs |
| Zhang | RF, VL, VM, TFL, BFL, BFS, ST, ADM | SVM | Zero crossings, slope change number, mean absolute value, waveform length | EMG window 150 ms decision time between 45 and 28 ms | Offline training, at least five reps of 30 s of data per condition | None; online classification | Transfemoral | Ramp/stair ascent, descent, level-ground walking, sitting, standing |
| Tkach | TA, PL, GASL, GASM | LDA | Mean absolute value, zero crossings, slope sign change, and waveform length | 250 ms window | 12 trials of each ambulation mode (stairs and ramps), and 12 trials of level-ground walking | None; offline classification only | Transtibial | Level-ground walking, stairs, ramps, and transitions |
| Tkach | TA, PL, GASL, GASM, VM, VL, RF, BF | LDA | Six coefficients from 6th order autoregressive model, mean absolute value, zero crossings, slope sign change, waveform length | 250 ms prediction window 50 ms overlap | 18 s for data for each movement class (one-DOF up to three-DOF, including no movement flexion, rotation, in/eversion) | None; Offline classification only | Transtibial | Virtual tasks for varying DOF ankle movements (rotation, flexion, in/eversion) |
| Hargrove | ST, SAR, TFL, ADM, GRA, VM, RF, VL, BFL | LDA | Mean absolute value, zero crossings, slope sign change, waveform length | 250 ms window, decisions every 50 ms | EMG for virtual cases, four reps of 3 s for each movement, no feedback provided | Online impedance control and virtual avatar feedback | Transfemoral | Knee flexion/extension, ankle plantar flexion/dorsiflexion,internal/external tibial/femoral rotation, relaxation |
| Miller | TA, GASM, VL, BF | LDA and SVM | Mean absolute value, variance, wavelength, number slope sign changes, zero crossings | Three subwindows with varying times based on heel strike and toe off | Six trials each activity, additional trials for different walking speeds | None; offline classification only | Transtibial | Level-ground walking three speeds, ramp/stair ascent/descent |
| Du | SAR, RF, VL, VM, GRA, BFL, ST, BFS, | Adaptive algorithm: EBA and learning from test data | Absolute value, slope sign changes, waveform length, zero crossings | 160ms sliding window | 15–10 trials of each activity | None; offline classification only | Transfemoral | Level-ground walking, ramps, stairs and their transitions |
| Zhang | RF, VL, VM, BFL, SAR, ST, ADM | Not provided | Not provided | Not provided | Offline training with ten reps of walking level-ground and ramp course | Impedance control | Transfemoral | Standing, level-ground walking ramp ascent/descent |
| Young | ST, BF, TFL, RF, VL, VM, SAR, ADM, GRA | DBN | Mean absolute value, waveform length, zero crossings, slope sign change, two autoregressive coefficients | Tested 0–450 ms windows with 20 ms sliding window | Offline training with 20 reps of locomotor circuit | None; offline classification | Transfemoral | Level-ground walking, ramps, stairs |
| Hargrove | ST, BF, TFL, RF, VL, VM, SAR, ADM, GRA | LDA + DBN | Not provided | Not provided | Offline training, 20 reps of locomotor circuit, level-ground walking at variable speeds, stopping/starting/turning | Impedance control | Transfemoral | Walking, stairs, ramps |
| Spanias | ST, ADM, TFL, RF | DBN | Mean absolute value, waveform length, zero crossings, slope sign changes, six autoregressive coefficients | 300 ms window | 20 reps of locomotion circuit | None; offline classification only | Transfemoral | Level-ground walking, ramp/straight ascent/descent, turn around |
| Spanias | ST, BF, TFL, RF, VL, VM, SAR, ADM, GRA | LDA with log-likelihood | Mean absolute value, waveform length, zero crossings, slope sign change, 1st two coefficients of 3rd order autoregressive model | 300 ms before toe off and heel strike | Offline training with 20 reps of locomotor circuit | Impedance control | Transfemoral | Level-ground walking, ramps, stairs |
| Liu | RF, VL, VM, BFL, SM, TFL, BFS | EBA, TSVM | Absolute value, signal length, slope sign changes, zero crossings | 160 ms window avg processing time 45 ms | Four 1 min trials of each locomotion mode collected in separate training session | EBA online impedance control | Transfemoral | Level-ground walking, ramps, stairs, and their transitions |
| Spanias | Pairs of electrodes over RF, TFL, ST, ADM | Eight DBN classifiers + LDA | Mean absolute value, waveform length, zero crossings, slope sign changes, six autoregressive coefficients from 6th order model PCA and ULDA to features totals | Transition modes 300 ms (210 ms before gait event and 90 ms after gait event) | Offline training with 1st session data across locomotor modes. Online training 2nd session through forward prediction. | Impedance control | Transfemoral | Level-ground walking, ramps, stairs over multiple days |
| Hussain | GM, GASL, GASM, TFL, RF, VL, VM, BFL, SOL, TA | SVM and LDA | Bispectrum for high order frequency spectrum/non-Gaussian info, RMS, zero crossings, histogram, integrated EMG, sum of square integral, waveform length, mean/median frequency, autoregression, and reduced with PCA | Non-overlapping 150 ms-300 ms windows | Offline training five trials for each walking mode | None; offline classification only | Transtibial and transfemoral | Slow, steady state, and fast walking, ramps |
Authors are listed in order of year published.
Muscle abbreviations: adductor magnus (ADM), adductor longus (ADDL), biceps femoris (BF), biceps femoris long head (BFL), biceps femoris short head (BFS), extensor digitorum brevis (EDB), gastrocnemius (GAS), gastrocnemius medialis (GASM), gastrocnemius lateralis (GASL), gluteus maximus (GMA), gluteus medius (GME), gracilis (GRA), peroneus longus (PL), rectus femoris (RF), sartorius (SAR), semimembranosus (SM), semitendinosus (ST), soleus (SOL), tensor fasciae latae (TFL), tibialis anterior (TA), vastus lateralis (VL), vastus medialis (VM).
Classifier abbreviations: support vector machine (SVM), neural network (NN), linear discriminant analysis (LDA), quadratic descriminant analysis (QDA), dynamic Bayesian network (DBN), entropy-based adaption (EBA), transductive support vector machine (TSVM).
Direct EMG control review summary.
| Author[ | Muscle[ | EMG decoding method | Modulated control parameters | Amputation type | Activity |
|---|---|---|---|---|---|
| Ha | Quadriceps, hamstring (unspecified) | Envelope (2 Hz cutoff), 20% MVC threshold QDA | Reference angular velocity | Transfemoral | Virtual tracking task (sitting) |
| Hoover | VL, VM, RF, BF | Envelope (5–10 Hz cutoff). | Flexion/extension torque | Transfemoral | Level-ground walking |
| Hoover | BF, RF, ST, VL, VM | Envelope (2.5 Hz cutoff) | Flexion/extension torque | Transfemoral | Stair ascent |
| Dawley | Quadriceps, hamstring, (unspecified) | Initial processing not described, principal component analysis (flexion/extension) | Reference angular velocity, joint stiffness | Transfemoral | Level-ground walking |
| Wang | GAS (unspecified head) | Low-pass filtered (15 Hz), rectified envelope (200 ms moving average window) | Plantar flexor torque gain (push-off only) | Transtibial | Level-ground walking |
| Alcaide-Aguirre | TA | Envelope (10 Hz cutoff) | Virtual object acceleration | Transtibial | Virtual tracking task (sitting) |
| Chen | TA, GAS (unspecified head) | Envelope (2.5 Hz cutoff), PCA (flexion/extension) | Reference angular velocity, joint stiffness | Transtibial | Virtual target hitting (sitting) |
| Huang | GASL | Envelope (2 Hz cutoff) | Pneumatic artificial muscle force | Transtibial | Level-ground walking |
| Huang | GASM or GASL | Envelope (2 Hz cutoff) | Pneumatic artificial muscle force | Transtibial | Level-ground walking |
| Huang | TA, GASM or GASL | Envelope (2 Hz cutoff) | Virtual object position | Transtibial | Virtual ballistic target hitting (sitting) |
| Clites | TA, GASL, TP, PL | Envelope (100 ms moving average window) | Flexion/extension torque | Transtibial | Virtual target hitting, stair ascent/descent, obstacle walking |
| Fleming | TA, GASL | Envelope (2 Hz cutoff) | Virtual spring stiffness | Transtibial | Virtual balancing inverted pendulum (sitting) |
| Huang | TA, GASL | Envelope (2 Hz cutoff) | Virtual cursor position | Transtibial | Virtual control input space filling (sitting) |
| Dimitrov | TA, GASM, GASL | Envelope (5 Hz cutoff), non-negative matrix factorization (125 ms windows) | Equilibrium angle, joint stiffness | Transtibial | Target hitting (standing), walking (with passive device) |
| Fleming | TA, GASL | Envelope (2 Hz cutoff) | Pneumatic artificial muscle force | Transtibial | Quiet standing (vision, no vision, foam and firm surfaces), Sit-to-stand, load transfer. |
Authors are listed in order of year published.
Muscle abbreviations: biceps femoris (BF), gastrocnemius (GAS), gastrocnemius medialis (GASM), gastrocnemius lateralis (GASL), peroneus longus (PL), rectus femoris (RF), semitendinosus (ST), soleus (SOL), tibialis anterior (TA), vastus lateralis (VL), vastus medialis (VM), tibialis posterior (TP).
Figure 1.Supervisory EMG control paradigm for robotic lower-limb prosthesis. In supervisory EMG control, EMG signals and gait events are used to classify the user’s locomotion mode (such as level-ground walking, stair ascent/descent, ramp ascent/descent). The classifier’s decision determines transitions between the predefined finite-states and thus the specified low-level control (e.g. impedance control) for prosthesis operation associated with the identified locomotion mode.
Figure 2.Direct EMG control paradigm for robotic lower-limb prosthesis. In direct EMG control, the magnitude of of EMG signals recorded from antagonistic residual muscles directly and continuously modulate the prosthesis joint dynamics. Various control laws can be used to continuously map EMG activity to ankle control torque to drive prosthesis dynamics. For example, EMG magnitude of residual ankle antagonistic muscles (u1 and u2) can activate an EMG-driven musculoskeletal model to estimate intended ankle control torque.
Figure 3.A human motor control framework (adopted from the framework reported in [109] to guide future research in myoelectric control of robotic lower limb prostheses. The actual state of the distal limb is disrupted after limb amputation (red dashed lines). When motor commands (EMG signals of residual muscles) are used to drive a robotic prosthetic limb, humans need to adapt internal model control parameters (the inverse model and forward model) via repetitive motor practice to minimize errors between the desired state and the predicted state, between the desired state and the actual state, and between the predicted state and actual state. This review presents the framework as a means to facilitate future research improving an amputee’s capability to produce appropriate motor commands (residual muscle EMG activity) and control the robotic prosthesis.