CN110720912B - 一种基于波形相关性的肌肉运动单元提取方法 - Google Patents
一种基于波形相关性的肌肉运动单元提取方法 Download PDFInfo
- Publication number
- CN110720912B CN110720912B CN201910968169.9A CN201910968169A CN110720912B CN 110720912 B CN110720912 B CN 110720912B CN 201910968169 A CN201910968169 A CN 201910968169A CN 110720912 B CN110720912 B CN 110720912B
- Authority
- CN
- China
- Prior art keywords
- waveform
- muscle movement
- movement unit
- template
- muscle
- 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.)
- Expired - Fee Related
Links
- 210000003205 muscle Anatomy 0.000 title claims abstract description 40
- 238000000605 extraction Methods 0.000 title claims abstract description 10
- 238000012886 linear function Methods 0.000 claims abstract description 10
- 238000000034 method Methods 0.000 claims abstract description 5
- 210000001087 myotubule Anatomy 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 3
- 238000003064 k means clustering Methods 0.000 claims description 3
- 230000003313 weakening effect Effects 0.000 claims description 3
- 230000004118 muscle contraction Effects 0.000 abstract description 6
- 230000008602 contraction Effects 0.000 abstract description 3
- 101000739754 Homo sapiens Semenogelin-1 Proteins 0.000 description 2
- 102100037550 Semenogelin-1 Human genes 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000036982 action potential Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Fuzzy Systems (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明提供了一种基于波形相关性的肌肉运动单元提取方法,首先根据肌肉运动单元发放波形传播延时特征,精确提取各个电极上的发放波形,然后根据发放波形,确定发放时刻,最后根据发放时刻,采用聚类方法,提取肌肉运动单元。针对肌肉动态收缩,采用线性函数动态更新模板,大大提高了提取肌肉运动单元的准确性。本发明实现简单,为肌肉动态收缩提取运动单元提供了新的思路,满足实际应用的需要。
Description
技术领域
本发明涉及一种基于波形相关性的肌肉运动单元提取方法。
背景技术
表面肌电信号(Surface EMG,SEMG)是利用表面电极从人体体表检测肌电信号,它具有无创性,操作简单的特点,易于接受,因此应用前景广阔。SEMG是由多个运动单元(Motor Unit,MU)兴奋而发放出来的运动单元动作电位序列(Motor Unit ActionPotential Train,MUAPT)经由皮肤、脂肪等构成的容积导体滤波后,在皮肤表面检测电极处的时间和空间上综合叠加的结果。目前,对sEMG信号分析处理主要针对的是肌肉静态收缩的表面肌电信号,对肌肉动态收缩的研究非常少,本发明就是针对肌肉动态收缩,提出一种肌肉MU提取方法。
发明内容
鉴于上述问题,本发明的目的在于提供一种基于波形相关性的肌肉运动单元提取方法,其特征是包括以下步骤:
步骤一,采集肌肉在动态收缩时阵列式表面肌电信号,对信号滤波,削弱干扰;
步骤二,提取每个通道上的肌肉运动单元发放波形;
步骤三,以肌肉运动单元发放波形幅值最大的电极M为基准,沿肌纤维方向,根据电极间距和传播速度,计算各电极波形延时;
步骤四,以电极M发放波形的时刻和延时,计算沿肌纤维方向各电极发放波形时间域,并且提取时间域上的波形W1、W2、W3、…、Wn;
步骤五,计算W1、W2、W3、…、Wn相似性,满足相似性条件,即可确认为肌肉运动单元发放时刻;
步骤六,提取所有肌肉运动单元发放时刻的电极M上的波形,采用K均值聚类方法将所有的发放时刻分类,并且动态更新分类模板,输出分类结果;
步骤七,对所有发放序列归类整理,剔除重复以及不合理的肌肉运动单元,优化结果。
优化的措施包括:
步骤六中动态更新分类模板,根据波形幅值单调变化,采用线性函数更新模板。
与现有技术相比,本发明提供一种基于波形相关性的肌肉运动单元提取方法,采用波形传播特征,提取所有的发放时刻,并且考虑了传播延时特征,精确提取各电极上的发放波形,为后续的分类奠定基础。针对肌肉动态收缩,采用线性函数动态更新模板,大大提高了提取肌肉运动单元的准确性。本发明实现简单,为肌肉动态收缩提取运动单元提供了新的思路。
附图说明
图1是本发明流程图。
具体实施方式
以下结合附图实例对本发明作进一步详细描述,本领域技术人员可由本说明书所揭露的内容轻易地实现。
如图1所示为本发明的流程图。
本发明的目的在于提供一种基于波形相关性的肌肉运动单元提取方法,其特征是包括以下步骤:
步骤一,采集肌肉在动态收缩时阵列式表面肌电信号,对信号滤波,削弱干扰。由于sEMG信号中包含各种干扰信号,预处理首先需要采用带通滤波器,保留10Hz--500Hz频段信号,然后采用陷波滤波器,滤除50Hz工频干扰。
步骤二,提取每个通道上的肌肉运动单元发放波形。提取时,设置峰值阈值,高于阈值,提取峰值时刻前后40ms长度的波形,否则不提取。
步骤三,以肌肉运动单元发放波形幅值最大的电极M为基准,沿肌纤维方向,根据电极间距和传播速度,计算各电极波形延时。在各个通道上提取的肌肉运动单元发放波形,都会对应时刻,在这些时刻,找到幅值最大的电极M。计算延时的发放就是电极间距除以传播速度,电极间距可以从电极上测量得到,传播速度采取典型值4m/s。
步骤四,以电极M发放波形的时刻和延时,计算沿肌纤维方向各电极发放波形时间域,并且提取时间域上的波形W1、W2、W3、…、Wn。该步骤只需要确定波形的起始点就可以获得W1、W2、W3、…、Wn,起始点是由电极M上波形对应的时刻为基准,再由延时调整即可得到。
步骤五,计算W1、W2、W3、…、Wn相似性,满足相似性条件,即可确认为肌肉运动单元发放时刻。相似性si由公式计算得到,pi是动态输入变化的波形,zi是模板,N是波形长度。设置相似性阈值,小于阈值,即可确认为肌肉运动单元发放时刻。
步骤六,提取所有肌肉运动单元发放时刻的电极M上的波形,采用K均值聚类方法将所有的发放时刻分类,并且动态更新分类模板,输出分类结果。由于M是幅值最大的电极,最能反映肌肉运动单元特征,用M上波形聚类,准确性高。动态更新分类模板,根据波形幅值单调变化,采用线性函数更新模板。波形单调增,采用单增线性函数更新模板;波形单调减,采用单减线性函数更新模板。具体实施时,采用幅值与系数相乘更新模板。
步骤七,对所有发放序列归类整理,剔除重复以及不合理的肌肉运动单元,优化结果。不合理的肌肉运动单元指发放时刻小于15毫秒间隔的序列,需剔除。
综上所述,本发明提供一种基于波形相关性的肌肉运动单元提取方法,采用波形传播特征,提取所有的发放时刻,并且考虑了传播延时特征,精确提取各电极上的发放波形。针对肌肉动态收缩,采用线性函数动态更新模板,大大提高了提取肌肉运动单元的准确性。本发明实现简单,为肌肉动态收缩提取运动单元提供了新的思路,满足实际应用的需要。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (1)
1.一种基于波形相关性的肌肉运动单元提取方法,其特征是包括以下步骤:
步骤一,采集肌肉在动态收缩时阵列式表面肌电信号,对信号滤波,削弱干扰;
步骤二,提取每个通道上的肌肉运动单元发放波形;
步骤三,以肌肉运动单元发放波形幅值最大的电极M为基准,沿肌纤维方向,根据电极间距和传播速度,计算各电极波形延时;
步骤四,以电极M发放波形的时刻和延时,计算沿肌纤维方向各电极发放波形时间域,并且提取时间域上的波形W1、W2、W3、…、Wn;
步骤六,提取所有肌肉运动单元发放时刻的电极M上的波形,采用K均值聚类方法将所有的发放时刻分类,并且动态更新分类模板,输出分类结果,动态更新分类模板,根据波形幅值单调变化,采用线性函数更新模板,波形单调增,采用单增线性函数更新模板;波形单调减,采用单减线性函数更新模板;
步骤七,对所有发放序列归类整理,剔除重复以及不合理的肌肉运动单元,优化结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910968169.9A CN110720912B (zh) | 2019-10-12 | 2019-10-12 | 一种基于波形相关性的肌肉运动单元提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910968169.9A CN110720912B (zh) | 2019-10-12 | 2019-10-12 | 一种基于波形相关性的肌肉运动单元提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110720912A CN110720912A (zh) | 2020-01-24 |
CN110720912B true CN110720912B (zh) | 2022-04-26 |
Family
ID=69219907
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910968169.9A Expired - Fee Related CN110720912B (zh) | 2019-10-12 | 2019-10-12 | 一种基于波形相关性的肌肉运动单元提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110720912B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101317794A (zh) * | 2008-03-11 | 2008-12-10 | 清华大学 | 多指多自由度假手用肌电控制能力检测训练方法 |
CN101507606A (zh) * | 2008-12-09 | 2009-08-19 | 上海理工大学 | 脑诱发电位的采集方法 |
CN104586383A (zh) * | 2015-01-12 | 2015-05-06 | 北京世纪今科医疗器械有限公司 | 一种心电波形分类方法及装置 |
CN104997508A (zh) * | 2015-08-19 | 2015-10-28 | 宁波工程学院 | 阵列式sEMG信号自动分解方法 |
CN105326501A (zh) * | 2015-12-10 | 2016-02-17 | 宁波工程学院 | 基于sEMG的肌肉疾病监测方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7407485B2 (en) * | 2004-06-08 | 2008-08-05 | Instrumentarium Corporation | Monitoring pain-related responses of a patient |
WO2012040402A2 (en) * | 2010-09-21 | 2012-03-29 | Somaxis Incorporated | Methods for promoting fitness in connection with electrophysiology data |
CN102961203B (zh) * | 2012-12-10 | 2015-04-22 | 杭州电子科技大学 | 基于emd样本熵的表面肌电信号识别方法 |
US20150320575A1 (en) * | 2014-05-07 | 2015-11-12 | The Regents Of The University Of California | Intuitive prosthetic interface |
-
2019
- 2019-10-12 CN CN201910968169.9A patent/CN110720912B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101317794A (zh) * | 2008-03-11 | 2008-12-10 | 清华大学 | 多指多自由度假手用肌电控制能力检测训练方法 |
CN101507606A (zh) * | 2008-12-09 | 2009-08-19 | 上海理工大学 | 脑诱发电位的采集方法 |
CN104586383A (zh) * | 2015-01-12 | 2015-05-06 | 北京世纪今科医疗器械有限公司 | 一种心电波形分类方法及装置 |
CN104997508A (zh) * | 2015-08-19 | 2015-10-28 | 宁波工程学院 | 阵列式sEMG信号自动分解方法 |
CN105326501A (zh) * | 2015-12-10 | 2016-02-17 | 宁波工程学院 | 基于sEMG的肌肉疾病监测方法 |
Non-Patent Citations (5)
Title |
---|
The Decomposition of Multi-channel Surface;He Jinbao;《2015 Chinese Automation Congress (CAC)》;20160118;第979-982段 * |
The exaction of motor unit action potential from multi-channel SEMG signals;He Jinbao;《2015 4th International Conference on Computer Science and Network Technology (ICCSNT)》;20160616;第1424-1427页 * |
基于FastICA和通道间相关的表面肌电信号分解研究;宁勇等;《航天医学与医学工程》;20170615(第03期);第39-45 * |
基于表面肌电信号参数的人体上肢运动研究;刘涛然等;《天津科技大学学报》;20090425(第02期);第42-45页 * |
基于运动单元的肌肉力估计新方法;何金保等;《航天医学与医学工程》;20151015(第05期);第4-9页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110720912A (zh) | 2020-01-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104799854B (zh) | 一种表面肌电采集装置及其肌电信号处理方法 | |
CN101599127B (zh) | 眼电信号的特征提取与识别方法 | |
CN102426651B (zh) | 人体前臂表面肌电信号采集及模式识别系统 | |
Phinyomark et al. | Investigating long-term effects of feature extraction methods for continuous EMG pattern classification | |
CN104997508A (zh) | 阵列式sEMG信号自动分解方法 | |
Fattah et al. | Identification of motor neuron disease using wavelet domain features extracted from EMG signal | |
CN108403108A (zh) | 基于波形优化的阵列式表面肌电信号分解方法 | |
CN108388855B (zh) | 一种针电极肌电信号分解方法 | |
CN110720912B (zh) | 一种基于波形相关性的肌肉运动单元提取方法 | |
CN106845348B (zh) | 一种基于手臂表面肌电信号的手势识别方法 | |
CN105956547A (zh) | 基于阵列式表面肌电信号平滑的分解方法 | |
CN110558980B (zh) | 一种肌肉动态收缩的运动单元提取方法 | |
McCool et al. | 1-D local binary patterns for onset detection of myoelectric signals | |
CN108403113B (zh) | 一种半自动的阵列式表面肌电信号分解方法 | |
CN105975917B (zh) | 面向强干扰的阵列式表面肌电信号分解方法 | |
KR101446183B1 (ko) | 상관 특징 분석 기반 실시간 emg 패턴 인식 방법 | |
CN105769186B (zh) | 基于复杂网络的肌电信号采集位置选取方法 | |
Alam et al. | Real-time classification of multi-channel forearm emg to recognize hand movements using effective feature combination and lda classifier | |
CN110558979B (zh) | 一种基于神经网络的肌肉运动单元提取方法 | |
CN110507324B (zh) | 一种肌肉运动单元搜索方法 | |
Ling-Ling et al. | Electromyographic movement pattern recognition based on random forest algorithm | |
CN110720910B (zh) | 一种基于相关性的肌肉运动单元搜索方法 | |
CN114098768A (zh) | 基于动态阈值和EasyTL的跨个体表面肌电信号手势识别方法 | |
KR20220158462A (ko) | 근전도 신호 기반 인식 정보 추출 시스템과 이를 이용한 근전도 신호 기반 인식 정보 추출 방법 | |
CN111476295A (zh) | 一种基于gs-lssvm的多通道生理信号体感姿势识别方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220426 |