CN110720911A - Muscle movement unit extraction method - Google Patents
Muscle movement unit extraction method Download PDFInfo
- 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
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/389—Electromyography [EMG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification 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
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.
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)
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 |
-
2019
- 2019-10-12 CN CN201910967780.XA patent/CN110720911B/en active Active
Patent Citations (15)
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)
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 |