A Novel HD-sEMG Preprocessing Method Integrating Muscle Activation Heterogeneity Analysis and Kurtosis-Guided Filtering for High-Accuracy Joint Force Estimation

This study proposes a novel preprocessing method integrating muscle activation heterogeneity analysis and kurtosis-guided filtering to realize high-accuracy surface electromyogr-aphy (sEMG)-based force estimation. A total of 10 subjects were recruited. Each subject performed isometric elbow flexion tasks at 20%, 40%, and 60% maximum voluntary contraction (MVC) target force levels, and the joint force and high-density sEMG (HD-sEMG) signals from biceps brachii and brachialis were collected synchronously. The force estimation model was built using three-order polynomial fitting technique. The input signal extraction of the force model, also named as the preprocessing of HD-sEMG signal, was carried out in the following procedures: first, HD-sEMG signals were decomposed by principal component analysis into principal components and weight vectors; second, the first several weight maps were segmented to obtain heterogeneity information by the Otsu and Moore-Neighbor tracing methods, and the principal component covering the most activated areas was selected; and last, a kurtosis-guided filter was designed to process the selected principal component to obtain the input signal. For the sake of comparison, the joint force estimation experiments based ON five preprocessing methods were conducted. The experimental results demonstrated that the proposed method obtained 52%, 53%, and 59% reduction in the mean root mean square difference at 20% MVC, 40% MVC, and 60% MVC force-level tasks, respectively, compared to the preprocessing method with the first principal component plus fixed parameter filtering. This proposed HD-sEMG pre-processing method has reliable neuromuscular electro-physiological foundation, and has good application value for realizing high-accuracy muscle/joint force estimation in the fields of rehabilitation engineering, sports biomechanics, and muscle disease diagnosis etc.

[1]  Miguel Angel Mañanas,et al.  Spatial distribution of HD-EMG improves identification of task and force in patients with incomplete spinal cord injury , 2016, Journal of NeuroEngineering and Rehabilitation.

[2]  N. Hogan,et al.  Relating agonist-antagonist electromyograms to joint torque during isometric, quasi-isotonic, nonfatiguing contractions , 1997, IEEE Transactions on Biomedical Engineering.

[3]  Dacheng Xiu,et al.  Using Principal Component Analysis to Estimate a High Dimensional Factor Model with High-Frequency Data , 2016 .

[4]  Stephen H. M. Brown,et al.  Less is more: high pass filtering, to remove up to 99% of the surface EMG signal power, improves EMG-based biceps brachii muscle force estimates. , 2004, Journal of electromyography and kinesiology : official journal of the International Society of Electrophysiological Kinesiology.

[5]  Fuchun Sun,et al.  sEMG-Based Joint Force Control for an Upper-Limb Power-Assist Exoskeleton Robot , 2014, IEEE Journal of Biomedical and Health Informatics.

[6]  Xun Chen,et al.  HD-sEMG-based research on activation heterogeneity of skeletal muscles and the joint force estimation during elbow flexion , 2018, Journal of neural engineering.

[7]  Changmok Choi,et al.  A Study on Estimation of Joint Force Through Isometric Index Finger Abduction With the Help of SEMG Peaks for Biomedical Applications , 2016, IEEE Transactions on Cybernetics.

[8]  P. Cavanagh,et al.  Electromechanical delay in human skeletal muscle under concentric and eccentric contractions , 1979, European Journal of Applied Physiology and Occupational Physiology.

[9]  Catherine Marque,et al.  Evaluation of muscle force classification using shape analysis of the sEMG probability density function: a simulation study , 2014, Medical & Biological Engineering & Computing.

[10]  Jaap H van Dieën,et al.  Increase in heterogeneity of biceps brachii activation during isometric submaximal fatiguing contractions: a multichannel surface EMG study. , 2014, Journal of neurophysiology.

[11]  Sofiane Boudaoud,et al.  Analysis of the sEMG/force relationship using HD-sEMG technique and data fusion: A simulation study , 2017, Comput. Biol. Medicine.

[12]  R. Merletti,et al.  Estimation of shape characteristics of surface muscle signal spectra from time domain data , 1995, IEEE Transactions on Biomedical Engineering.

[13]  Naomi S. Altman,et al.  Points of Significance: Principal component analysis , 2017, Nature Methods.

[14]  Andreas Daffertshofer,et al.  Improving EMG-based muscle force estimation by using a high-density EMG grid and principal component analysis , 2006, IEEE Transactions on Biomedical Engineering.

[15]  Claudio Castellini,et al.  Evaluating subsampling strategies for sEMG-based prediction of voluntary muscle contractions , 2013, 2013 IEEE 13th International Conference on Rehabilitation Robotics (ICORR).

[16]  N. Otsu A threshold selection method from gray level histograms , 1979 .

[17]  B. Qiu,et al.  An isometric muscle force estimation framework based on a high-density surface EMG array and an NMF algorithm , 2017, Journal of neural engineering.

[18]  Xu Zhang,et al.  Muscle-tendon units localization and activation level analysis based on high-density surface EMG array and NMF algorithm , 2016, Journal of neural engineering.

[19]  R. A. Ekstrom,et al.  Electromyographic analysis of core trunk, hip, and thigh muscles during 9 rehabilitation exercises. , 2007, The Journal of orthopaedic and sports physical therapy.

[20]  S. Boudaoud,et al.  Muscle force estimation using data fusion from high-density SEMG grid , 2013, 2013 2nd International Conference on Advances in Biomedical Engineering.

[21]  Toshio Tsuji,et al.  A Variance Distribution Model of Surface EMG Signals Based on Inverse Gamma Distribution , 2017, IEEE Transactions on Biomedical Engineering.

[22]  Evelyn Morin,et al.  Force Modelling of Upper Limb Biomechanics Using Ensemble Fast Orthogonal Search on High-Density Electromyography , 2016, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[23]  Xiang Chen,et al.  ICA-based muscle–tendon units localization and activation analysis during dynamic motion tasks , 2018, Medical and Biological Engineering and Computing.

[24]  S Boudaoud,et al.  On the benefits of using HD-sEMG technique for estimating muscle force , 2015, Computer methods in biomechanics and biomedical engineering.

[25]  IMPROVED SEARCH OF PRINCIPAL COMPONENT ANALYSIS DATABASES FOR SPECTRO-POLARIMETRIC INVERSION , 2013, 1307.0061.

[26]  R. Scott,et al.  A Nonstationary Model for the Electromyogram , 1977, IEEE Transactions on Biomedical Engineering.

[27]  Philip E. Martin,et al.  The relationship between cadence and lower extremity EMG in cyclists and noncyclists. , 1995, Medicine and science in sports and exercise.

[28]  D. Stegeman,et al.  Redundancy or heterogeneity in the electric activity of the biceps brachii muscle? Added value of PCA-processed multi-channel EMG muscle activation estimates in a parallel-fibered muscle. , 2013, Journal of electromyography and kinesiology : official journal of the International Society of Electrophysiological Kinesiology.

[29]  Mohan P. Pradhan,et al.  Contour Line Tracing Algorithm for Digital Topographic Maps , 2010 .

[30]  Aydin Akan,et al.  Prediction of externally applied forces to human hands using frequency content of surface EMG signals , 2010, Comput. Methods Programs Biomed..

[31]  Parvin Mousavi,et al.  EMG-force modeling using parallel cascade identification. , 2012, Journal of electromyography and kinesiology : official journal of the International Society of Electrophysiological Kinesiology.

[32]  V. Medved,et al.  Locomotion diagnostics: Some neuromuscular and robotic aspects , 1991, IEEE Engineering in Medicine and Biology Magazine.

[33]  I. Kingma,et al.  Towards optimal multi-channel EMG electrode configurations in muscle force estimation: a high density EMG study. , 2005, Journal of electromyography and kinesiology : official journal of the International Society of Electrophysiological Kinesiology.