CN110720911A - Muscle movement unit extraction method - Google Patents

Muscle movement unit extraction method Download PDF

Info

Publication number
CN110720911A
CN110720911A CN201910967780.XA CN201910967780A CN110720911A CN 110720911 A CN110720911 A CN 110720911A CN 201910967780 A CN201910967780 A CN 201910967780A CN 110720911 A CN110720911 A CN 110720911A
Authority
CN
China
Prior art keywords
muscle movement
movement units
muscle
signals
extraction method
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910967780.XA
Other languages
Chinese (zh)
Other versions
CN110720911B (en
Inventor
何金保
胡庆波
骆再飞
周世官
李国君
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ningbo University of Technology
Original Assignee
Ningbo University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Ningbo University of Technology filed Critical Ningbo University of Technology
Priority to CN201910967780.XA priority Critical patent/CN110720911B/en
Publication of CN110720911A publication Critical patent/CN110720911A/en
Application granted granted Critical
Publication of CN110720911B publication Critical patent/CN110720911B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/389Electromyography [EMG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems

Abstract

The invention provides a muscle movement unit extraction method, which comprises the steps of firstly collecting multichannel surface electromyogram signals, then respectively extracting muscle movement units from two parts of surface electrode signals of odd-numbered channels and even-numbered channels, finally sorting and issuing results, confirming the movement units found by the two parts at the same time, and confirming the muscle movement units found by the odd-numbered channels or the even-numbered channels through differential signals. The muscle movement unit extraction method provided by the invention can effectively improve the correctness of the result by mutual verification through the respective extraction results of the odd and even part electrodes. The invention is simple to realize and meets the requirement of practical application.

Description

Muscle movement unit extraction method
Technical Field
The invention relates to a muscle movement unit extraction method.
Background
The issue of a Motor Unit (MU) is the basis for forming a surface electromyogram (surface EMG, sEMG), and obtaining detailed information of a Motor Unit through the sEMG is helpful for understanding the neuromuscular working mechanism. The method for extracting MU researched by scholars at home and abroad comprises the following steps: a K-means clustering algorithm, a template matching method, an Artificial Neural Network (ANN) algorithm, a real-time linear aliasing blind signal separation algorithm, an independent component analysis and reduction (ICA), a convolution kernel compensation algorithm and the like. The effects of the methods can only be verified through simulation, the effect of extracting MU from the real sEMG is difficult to verify, and in general, how to verify the extraction accuracy is one of the difficulties in the electromyography research field because the extraction of the myomotor unit by using the multichannel surface electromyography is still in the exploration stage.
Disclosure of Invention
In view of the above problems, an object of the present invention is to provide a muscle motion unit extraction method, including the steps of:
collecting multichannel surface electromyographic signals under muscle contraction, filtering the signals and weakening interference;
step two, extracting surface electromyographic signals of odd number channels, and calculating an inverse matrix of a signal cross-correlation matrix and a distribution sequence xid(n);
Extracting surface electromyographic signals of even number channels, and calculating an inverse matrix of a signal cross-correlation matrix and a distribution sequence xie(n);
Step four, the issuing sequence xid(n),ξe(n) comparing to confirm that the same dispensing sequence is a muscle movement unit;
step five, determining the positions of the muscle movement units of the rest issuing sequences by adopting differential signals, determining the muscle movement units with the positions outside the surface electrodes, and if not, not identifying the muscle movement units;
and step six, classifying and sorting all the movement units, removing unreasonable muscle movement units and optimizing the result.
Because the muscle movement unit extracted from the surface electromyogram signal has no gold standard, and whether the extraction result is correct or not is not checked, the muscle movement unit extraction method provided by the invention can effectively improve the correctness of the result, and mutual verification is further carried out through the respective extraction results of the odd-even two parts of electrodes. For muscle movement units that are not found in both parts at the same time, as muscle movement units that may be present outside the area covered by the surface electrodes, only some of them may be detected while the other part is not. Therefore, for this case, the differential signal is further employed to extract the muscle motion unit. The invention is simple to realize and greatly improves the accuracy of muscle movement unit extraction.
Drawings
FIG. 1 is a flow chart of the present invention.
Fig. 2 is a schematic surface electrode numbering diagram according to an embodiment of the invention.
Detailed Description
The present invention is further described in detail below with reference to the accompanying drawings, and can be easily implemented by those skilled in the art from the disclosure of the present specification.
Fig. 1 shows a flow chart of the present invention. The invention provides a muscle movement unit extraction method, which is characterized by comprising the following steps:
collecting multichannel surface electromyographic signals under muscle contraction, filtering the signals and weakening interference. Because the sEMG signal contains various interference signals, the preprocessing firstly needs to adopt a band-pass filter, reserves signals of a frequency range of 10 Hz-500 Hz, and then adopts a notch filter to filter 50Hz power frequency interference.
Step two, extracting surface electromyographic signals of odd number channels, and calculating an inverse matrix of a signal cross-correlation matrix and a distribution sequence xid(n) of (a). The parity order lane distribution is shown in fig. 2. Firstly, a cross-correlation matrix and a cross-correlation matrix inverse matrix of odd-number channel signals are calculated, wherein the cross-correlation matrix is expressed as:
C=E(S(n)ST(n))
where n is the sampling instant, S (n) is the multichannel signal at the nth sampling instant, ST(n) is the signal transposition for the nth sampling instant, E (-) is the order expectation, the inverse C of the cross correlation matrix is calculated-1I.e. by
C-1=[E(S(n)ST(n))]-1
The sampling time n is the median of sEMG signal energy, and the energy is calculated according to the following formula:
Δ=ST(n)C-1S(n)
the time n corresponding to the energy median is taken0. Obtaining a distribution sequence according to the signal correlation, and calculating a signal distribution sequence xi of an odd number channeld(n):
ξd(n)=ST(n0)C-1S(n)
Extracting surface electromyographic signals of even number channels, and calculating an inverse matrix of a signal cross-correlation matrix and a distribution sequence xieAnd (n) the method is the same as the second step.
Step four, the issuing sequence xid(n),ξe(n) comparing to confirm that the same dispensing sequence is a muscle motor unit. And confirming the muscle movement units which can be found in the surface electromyographic signals of the channels with the even serial numbers and the channels with the odd serial numbers, wherein in the comparison process, the error within 3 is allowed to exist at the moment of issuing the sequence.
And step five, determining the positions of the muscle movement units by adopting differential signals for the rest issuing sequences, and determining the muscle movement units with the positions outside the surface electrodes, otherwise, not identifying the muscle movement units. When the muscle movement units are distributed outside the electrodes, the muscle movement units can be found only on even-numbered channels or odd-numbered channels, so that the positions of the muscle movement units are determined through differential signals, and then the muscle movement units are determined.
And step six, classifying and sorting all the movement units, removing unreasonable muscle movement units and optimizing the result. For dispensing rejections with dispensing times less than 15ms, it is possible to supplement for times that are significantly missing.
In summary, the muscle movement unit extraction method provided by the invention can effectively improve the accuracy of the result, and further verify each other through the respective extraction results of the odd and even electrodes. For muscle movement units that may be present outside the surface electrode coverage area, only some of them may be detected while the other part is not. Therefore, for this case, the differential signal is further employed to extract the muscle motion unit. The invention is simple to realize and greatly improves the accuracy of muscle movement unit extraction.
The foregoing embodiments are merely illustrative of the principles and utilities of the present invention and are not intended to limit the invention. Any person skilled in the art can modify or change the above-mentioned embodiments without departing from the spirit and scope of the present invention. Accordingly, it is intended that all equivalent modifications or changes which can be made by those skilled in the art without departing from the spirit and technical spirit of the present invention be covered by the claims of the present invention.

Claims (1)

1. A muscle movement unit extraction method is characterized by comprising the following steps:
collecting multichannel surface electromyographic signals under muscle contraction, filtering the signals and weakening interference;
step two, extracting surface electromyographic signals of odd number channels, and calculating an inverse matrix of a signal cross-correlation matrix and a distribution sequence xid(n);
Extracting surface electromyographic signals of even number channels, and calculating an inverse matrix of a signal cross-correlation matrix and a distribution sequence xie(n);
Step four, the issuing sequence xid(n),ξe(n) comparing to confirm that the same dispensing sequence is a muscle movement unit;
step five, determining the positions of the muscle movement units of the rest issuing sequences by adopting differential signals, determining the muscle movement units with the positions outside the surface electrodes, and if not, not identifying the muscle movement units;
and step six, classifying and sorting all the movement units, removing unreasonable muscle movement units and optimizing the result.
CN201910967780.XA 2019-10-12 2019-10-12 Muscle movement unit extraction method Active CN110720911B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910967780.XA CN110720911B (en) 2019-10-12 2019-10-12 Muscle movement unit extraction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910967780.XA CN110720911B (en) 2019-10-12 2019-10-12 Muscle movement unit extraction method

Publications (2)

Publication Number Publication Date
CN110720911A true CN110720911A (en) 2020-01-24
CN110720911B CN110720911B (en) 2022-08-30

Family

ID=69220069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910967780.XA Active CN110720911B (en) 2019-10-12 2019-10-12 Muscle movement unit extraction method

Country Status (1)

Country Link
CN (1) CN110720911B (en)

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040225211A1 (en) * 2001-10-24 2004-11-11 Gozani Shai N. Devices and methods for the non-invasive detection of spontaneous myoelectrical activity
EP1854404A2 (en) * 2006-05-09 2007-11-14 GE Healthcare Finland Oy Monitoring of the state of the central nervous system of a subjet
WO2008106232A1 (en) * 2007-03-01 2008-09-04 Neurometrix, Inc. Estimation of f-wave times of arrival (toa) for use in the assessment of neuromuscular function
US20090281448A1 (en) * 2008-05-10 2009-11-12 Neural Signals, Inc. Wireless Skin Surface Potential Sensing System and Method
US20100292617A1 (en) * 2009-05-15 2010-11-18 Kin Fong Lei method and system for quantifying an intention of movement of a user
WO2013144866A1 (en) * 2012-03-27 2013-10-03 B10Nix S.R.L. System for the acquisition and analysis of muscle activity and operation method thereof
CA2939949A1 (en) * 2014-02-28 2015-09-03 Powell Mansfield, Inc. Systems, methods and devices for sensing emg activity
CN104997508A (en) * 2015-08-19 2015-10-28 宁波工程学院 Automatic decomposition method of array type sEMG (surface EMG) signal
CN105147251A (en) * 2015-08-19 2015-12-16 宁波工程学院 Muscle fatigue dynamic prediction method based on multi-channel sEMG
CN105326501A (en) * 2015-12-10 2016-02-17 宁波工程学院 Muscle disease monitoring method based on sEMG
JP2016063995A (en) * 2014-09-25 2016-04-28 国立大学法人九州工業大学 Evaluation system and evaluation method of action amount and action efficiency of muscle tissue system and nervous system, and program thereof
CN105975917A (en) * 2016-04-28 2016-09-28 宁波工程学院 Array type surface electromyogram signal decomposition method facing high interference
CN107518895A (en) * 2016-06-22 2017-12-29 宁波工程学院 Towards the same muscle detection method of rehabilitation training
CN108403115A (en) * 2018-02-11 2018-08-17 宁波工程学院 A kind of muscular movement element number method of estimation
JP2019030520A (en) * 2017-08-08 2019-02-28 日本電信電話株式会社 Myoelectricity measurement apparatus, method, and program

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040225211A1 (en) * 2001-10-24 2004-11-11 Gozani Shai N. Devices and methods for the non-invasive detection of spontaneous myoelectrical activity
EP1854404A2 (en) * 2006-05-09 2007-11-14 GE Healthcare Finland Oy Monitoring of the state of the central nervous system of a subjet
WO2008106232A1 (en) * 2007-03-01 2008-09-04 Neurometrix, Inc. Estimation of f-wave times of arrival (toa) for use in the assessment of neuromuscular function
US20090281448A1 (en) * 2008-05-10 2009-11-12 Neural Signals, Inc. Wireless Skin Surface Potential Sensing System and Method
US20100292617A1 (en) * 2009-05-15 2010-11-18 Kin Fong Lei method and system for quantifying an intention of movement of a user
WO2013144866A1 (en) * 2012-03-27 2013-10-03 B10Nix S.R.L. System for the acquisition and analysis of muscle activity and operation method thereof
CA2939949A1 (en) * 2014-02-28 2015-09-03 Powell Mansfield, Inc. Systems, methods and devices for sensing emg activity
JP2016063995A (en) * 2014-09-25 2016-04-28 国立大学法人九州工業大学 Evaluation system and evaluation method of action amount and action efficiency of muscle tissue system and nervous system, and program thereof
CN104997508A (en) * 2015-08-19 2015-10-28 宁波工程学院 Automatic decomposition method of array type sEMG (surface EMG) signal
CN105147251A (en) * 2015-08-19 2015-12-16 宁波工程学院 Muscle fatigue dynamic prediction method based on multi-channel sEMG
CN105326501A (en) * 2015-12-10 2016-02-17 宁波工程学院 Muscle disease monitoring method based on sEMG
CN105975917A (en) * 2016-04-28 2016-09-28 宁波工程学院 Array type surface electromyogram signal decomposition method facing high interference
CN107518895A (en) * 2016-06-22 2017-12-29 宁波工程学院 Towards the same muscle detection method of rehabilitation training
JP2019030520A (en) * 2017-08-08 2019-02-28 日本電信電話株式会社 Myoelectricity measurement apparatus, method, and program
CN108403115A (en) * 2018-02-11 2018-08-17 宁波工程学院 A kind of muscular movement element number method of estimation

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
何金保等: "基于运动单元的肌肉力估计新方法", 《航天医学与医学工程》 *
宁勇等: "基于FastICA和通道间相关的表面肌电信号分解研究", 《航天医学与医学工程》 *
邱青菊: "表面肌电信号的特征提取与模式分类研究", 《中国优秀硕士学位论文电子数据库》 *

Also Published As

Publication number Publication date
CN110720911B (en) 2022-08-30

Similar Documents

Publication Publication Date Title
CN108491077B (en) Surface electromyographic signal gesture recognition method based on multi-stream divide-and-conquer convolutional neural network
CN104997508B (en) Array sEMG signal automatic classifying methods
CN108403108B (en) Array type surface electromyographic signal decomposition method based on waveform optimization
CN107203692A (en) The implementation method of atrial fibrillation detection based on depth convolutional neural networks
CN110598676B (en) Deep learning gesture electromyographic signal identification method based on confidence score model
CN104166548A (en) Deep learning method based on motor imagery electroencephalogram data
CN109359619A (en) A kind of high density surface EMG Signal Decomposition Based method based on convolution blind source separating
CN107822629A (en) The detection method of extremity surface myoelectricity axle
CN112684891A (en) Electroencephalogram signal classification method based on multi-source manifold embedding migration
CN107518895B (en) Same muscle detection method for rehabilitation training
CN112633195A (en) Myocardial infarction identification and classification method based on frequency domain features and deep learning
CN105956547B (en) Based on the smooth decomposition method of array surface electromyogram signal
CN105326501A (en) Muscle disease monitoring method based on sEMG
CN113158964A (en) Sleep staging method based on residual learning and multi-granularity feature fusion
CN108388855B (en) Needle electrode electromyographic signal decomposition method
CN114732424A (en) Method for extracting complex network attribute of muscle fatigue state based on surface electromyographic signal
CN106264460B (en) The coding/decoding method and device of cerebration multidimensional time-series signal based on self study
CN110720911B (en) Muscle movement unit extraction method
Sultornsanee et al. Classification of electromyogram using recurrence quantification analysis
CN112698720B (en) Movement imagery identification method based on mutual information feature extraction and multi-person fusion
CN105975917B (en) Towards strongly disturbing array Decomposition Surface EMG method
CN103070681B (en) A kind of ventricular late potential separation method based on sparse component analysis
CN110738093B (en) Classification method based on improved small world echo state network electromyography
CN111783719A (en) Myoelectric control method and device
CN110558980A (en) Muscle dynamic contraction movement unit extraction method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant