相关文章推荐
另类的墨镜  ·  实战总结 Vue ...·  5 月前    · 
乐观的四季豆  ·  NLP ...·  1 年前    · 
The .gov means it’s official. Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site. The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely. As a library, NLM provides access to scientific literature. Inclusion in an NLM database does not imply endorsement of, or agreement with, the contents by NLM or the National Institutes of Health.
Learn more: PMC Disclaimer
Sheng Wu Yi Xue Gong Cheng Xue Za Zhi. 2018 Dec; 35(6): 953–958.
PMCID: PMC9935191

Language: Chinese | English

一种心电干扰下的表面肌电起止点检测方法

An endpoint-detection algorithm of surface electromyography insensitive to electrocardiogram interference

伟钊 徐

华南理工大学 自动化科学与工程学院(广州 510641), School of Automation Science and Engineering, South China University of Technology, Guangzhou 510641, P.R.China

Find articles by 伟钊 徐

鸿强 莫

华南理工大学 自动化科学与工程学院(广州 510641), School of Automation Science and Engineering, South China University of Technology, Guangzhou 510641, P.R.China

Find articles by 鸿强 莫

联房 田

华南理工大学 自动化科学与工程学院(广州 510641), School of Automation Science and Engineering, South China University of Technology, Guangzhou 510641, P.R.China

Find articles by 联房 田

德苗 欧

华南理工大学 自动化科学与工程学院(广州 510641), School of Automation Science and Engineering, South China University of Technology, Guangzhou 510641, P.R.China 华南理工大学 自动化科学与工程学院(广州 510641), School of Automation Science and Engineering, South China University of Technology, Guangzhou 510641, P.R.China

莫鸿强,Email: nc.ude.tucs@omgnaiqh

不管是 sEMG 的时域分析还是频域分析,肌电起止点(即动作起止时刻)的准确判断均是保证分析结果正确性的关键前提。当数据量不大时,可以由经验丰富的分析员人工选定 sEMG 的起止点,但数据量大时这种方法的效率很低。sEMG 与语音信号类似,都属于短时平稳信号 [ 2 - 3 ] ,可以借鉴语音信号的端点检测算法。语音识别中常见的端点检测方法是基于短时能量的端点检测方法 [ 2 ] ,将信号按固定帧长和帧移分帧,根据单帧信号均方根值是否超过阈值来判断信号端点。直接运用这一算法来检测 sEMG 端点,当被观察肌电幅值较强时,可以有效地检测到动作的起止点;但若在测量位置处被观察肌电幅值较弱而心电(electrocardiogram,ECG)干扰幅值相对较大时,其容易因 ECG 干扰而误判动作起止点。

克服 ECG 干扰影响的最直接的方法是抑制所测肌电信号中的 ECG 成分,即根据 ECG 和 sEMG 信号的频谱差别设计合适的滤波器以滤除 ECG。决定 ECG 幅值的主要能量集中于 30 Hz 以下的频段 [ 4 - 5 ] ,因此,作截止频率为 20 Hz 或 30 Hz 的高通滤波可以大大削弱 ECG 干扰 [ 6 - 8 ] 。不过这种高通滤波方法势必也滤去一部分 sEMG 有效信号,影响肌电分析结果。为此,研究人员提出了很多新的方法。例如:Mak 等 [ 9 ] 提出了利用独立成分分析法(independent component analysis,ICA)减少 ECG 对 sEMG 的干扰;Costa 等 [ 10 ] 提出了单通道 ICA 分离的方法;Nougarou 等 [ 11 ] 提出离散小波变换(discrete wavelet transform,DWT)与 ICA 混合算法;周权等 [ 12 ] 则提出结合小波变换(wavelet transform,WT)和能量熵的混合算法。这些新方法在很大程度上能有效地抑制 ECG 波形,保证了 sEMG 分析结果的准确性,但计算量往往较大,在穿戴式设备上实现有一定的难度。当前移动医疗健康监测设备的应用日益广泛,受成本、体积、功耗等因素限制,移动医疗设备通常应采用计算量较小的算法。

鉴于此,本文借鉴语音识别的端点检测方法,建立一种计算量小的 sEMG 起止点检测的算法,其综合利用 sEMG 短时过零率高于 ECG 以及 sEMG 强度高于背景噪声的特点,确定 sEMG 的起点和终点。腹直肌的 sEMG 测量部位靠近心脏,其 sEMG 信号受 ECG 干扰比较严重,因此本文研究对象以腹直肌的 sEMG 为主,但不局限于腹直肌的 sEMG。

1. 算法原理及步骤

1.1. 传统的均方根阈值法及其缺点

记第 n 时刻 sEMG 值为 为所需分析的数据的点数。取帧长为 ,根据式(1)计算第 k 帧的均方根值 作四舍五入的取整运算。

由于 sEMG 信号幅值显著高于背景噪声(无 sEMG 时)的幅值,故可根据背景噪声的统计参数设定阈值以判断 sEMG 的起止点。选取一段足够长(例如 32 帧)的背景噪声信号,根据式(1)计算每一帧噪声的均方根值,再求所有帧噪声均方根值的均值 。通常假设背景噪声满足高斯分布,故可根据其统计特性设置阈值,例如可取阈值为 ,当第 k 帧的 sEMG 的均方根值大于该阈值,则认为该帧信号非背景噪声而为 sEMG 信号。

图 1a 显示的是一段采样频率为 2 kHz,并经过通带为[10 Hz, 700 Hz]的二阶带通滤波器和一个二阶工频陷波器滤波后的腹直肌 sEMG 数据,其中椭圆圈中所标示的为 ECG 片段。根据上述均方根阈值方法对其进行分析,所得结果如 图 1b 所示,竖线表示检测到的 sEMG 起止点位置。可以发现,sEMG 信号起止点的检测结果并不准确,将一部分 ECG 片段当成 sEMG 检测出来。由于该段信号中 ECG 幅值较高,即便调大阈值,检测结果亦不理想。增大帧长至远大于单次 ECG 持续时间,则可减小误判的概率,但是过度增大帧长又会加大 sEMG 起止点的定位误差。

sEMG and endpoints detected by means of a classical root-mean-square method

sEMG 信号及均方根阈值起止点检测结果

a. sEMG (ECGs marked out within the ovals);b. endpoints found (marked out by vertical lines)

a. sEMG 信号(椭圆圈中为 ECG 片段);b. sEMG 信号起止点检测结果(竖线表示)

1.2. 基于短时能量和短时过零率的 sEMG 信号起止点检测方法

语音信号处理中,通常联合短时能量和短时过零率以克服高幅值噪声的干扰,实现准确的端点检测 [ 2 ] 。与语音信号一样,sEMG 亦具有短时平稳特性,此外,sEMG 短时过零率高于 ECG 而强度高于背景噪声,故此,本文借鉴语音信号的端点检测方法实现 sEMG 的起止点检测,其主要计算过程包括短时能量包络和短时过零率计算及阈值判断两部分。

1.2.1. 基于短时能量的检测

考虑到均方根计算量偏大且绝对值能起到与均方根相近的效果,语音处理中通常也以求绝对值之和的方式来确定单帧的短时能量。本文方法亦以绝对值方式计算单帧的短时能量。记第 j 帧的短时能量值为 ,其值根据式(2)计算。 。 图 2 描述了本文所采取的短时能量计算过程。

Determination of short-time energy

分帧计算短时能量示意图

帧长 需根据采样频率确定。提高起止点时间定位的精度应选择较小的帧长,但瞬时尖峰噪声往往会使个别相邻两帧的短时能量值出现很大的差别,加大了确定短时能量阈值的难度。本文采用先对帧短时能量 阶低通滤波得短时能量包络 作滑动平均处理得 的方法,以尽可能去除 曲线中存在的高幅值的峰、谷点。所采用的 阶低通滤波的脉冲传递函数如式(3)所示。 帧的能量包络值。应说明的是,式(4)中取 2 的负整数幂次为滤波器参数,故可通过移位实现乘法运算,有利于降低计算量。但这些参数的选择并不唯一,只要低通滤波截止频率合适,均可较好地求得包络线,消除短时能量曲线中可能出现的尖峰。 作如式(5)所示的滑动平均处理,得平滑包络线幅值 的选取将会在 1.2.3 节介绍。

1.2.2. 基于短时过零率的检测

本文所述的短时过零率定义为一窗口(即连续数帧)内肌电信号穿越零值的次数。过零率决定于信号主要分量的频率。健康成人 ECG 的主频带虽落于 0.05~100 Hz 之间 [ 13 ] ,但由于其以 QRS 波群为主,频谱分布较为集中,ECG 能量集中在 30 Hz 以下频段 [ 4 - 5 ] ;而 sEMG 信号的频段为 5~450 Hz [ 14 ] ,且能量的频段分布不如 ECG 那么集中。因此,选取长度合适的窗口来观察 sEMG 和 ECG 波形,当采样频率相同时,两种信号在单个窗口内穿越零值的次数将有显著差别,ECG 穿越零值次数比 sEMG 少,可利用这一特点基于短时过零率检测 sEMG 的起止点。

计算短时过零率选取的窗长要合适,窗长过小难以确定合适的短时过零率阈值,过大则会增大计算量,且起止点定位精度低。分别记 ECG 和 sEMG 主要分量对应的平均频率为 。记窗长为 ,则短时过零次数 及采样频率 的关系满足式(6)。 0.0040.240.80.0080.481.60.0160.963.20.0321.685.60.0402.408.00.0563.3611.2

为确认 表 1 结果,可进一步做如下统计:测量背景噪声幅值相对较小的 ECG 与 sEMG 信号,计算在不同 值下相应的短时过零次数 ;运用 sEMG 采集仪对十位处于静止状态的健康成年测试者(九位男性一位女性,年龄在 20~25 岁之间)采集 sEMG 信号,每位测试者采集三次,每次五组数据。求出 的平均值和标准差,画出折线图如 图 4 所示,其中折线代表的是 取不同值时过零次数 的平均值,垂直线段表示根据过零次数标准差确定的过零次数上下限。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-35-6-953-4.jpg

for different values of length of window/sampling frequency (

窗长/采样频率( )取不同值时的

结合 表 1 图 4 可知, 时,根据短时过零次数 即可较好地区分 ECG 和 sEMG。但窗长 越大,计算量越高,故折中考虑选取 。此时,如果 =2 kHz,则短时过零次数 ,短时过零率阈值 可选定为 4。对于不同测试者,当 时,过零次数均值波动不大,如 图 5 所示,因此所选阈值参数具有一定的普遍适用性。具体应用时,若 取 4 依然容易因 ECG 干扰导致误判,则可选择更大的短时过零率阈值;不过阈值一般不超过 8,以免增加计算量。

Comparison of among ten adults (

十位测试者 的对比 (

1.2.3. 判断阈值的确定

如 1.2.1 和 1.2.2 节所述,算法实现时需先确定短时能量检测的阈值 和短时过零率的阈值 ,其中短时过零率的阈值 由 1.2.2 节确定,还需要确定短时能量检测的阈值 。由于阈值 用于区分背景噪声和 EMG 信号,且前者幅值远小于后者,因此 可根据实测信号包络线幅值的平均值来确定,如式(7)所示。

2. 实验结果

应用与广州雪利昂公司合作开发的 sEMG 信号采集仪(技术指标满足中国医药行业标准 YY/T 1095-2015 对肌电生物反馈仪的要求)采集十例处于静止状态健康成年测试者的腹直肌收缩-舒张动作的 sEMG 信号,其中九例男性,一例女性,年龄在 20~25 岁之间。每例测试者采集三组数据,每组至少包含三次收缩-舒张动作。采集环境条件保持一致,且每次测量完成后让受测试者休息一段时间。对采集所得的三组数据检测 sEMG 的起止点,记录处理过程的中间结果:短时过零率、短时能量 、短时能量包络 以及平滑包络线幅值 ,并记录起止点的最终检测结果。将检测结果中端点检测不受 ECG 干扰影响的记为“正确判断”,反之记为“有误判断”,据此统计得检测结果的样本数为 86 个,准确率为 95.6%。 图 6 给出了其中三位测试者的起止点端点的测试结果。如 图 6 所示,ECG 片段对应的短时过零率与 sEMG 的短时过零率差别很大;另外,从代表短时能量的曲线可以发现,背景噪声的短时能量很低,ECG 片段和 sEMG 的短时能量较大。 图 6 结果表明,该算法基本克服了 ECG 干扰的影响,可较准确地检测出 sEMG 的起止点(图中竖线标注)。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-35-6-953-6.jpg

Endpoint-detection results of rectus abdominis muscle’s sEMG

腹直肌 sEMG 的端点检测效果

a. the first group; b. the second group; c. the third group. Sub-graphs presented from top to bottom are origin signal, short-time zero-crossing rate, short-time energy ( ), the envelope of short-time energy ( ), the average envelope of short-time energy ( ), and the endpoint-detection result of sEMG

a.第一个测试者;b.第二个测试者;c.第三个测试者。各子图中曲线由上至下依次是原始信号、短时过零率、短时能量 、短时能量包络 、平滑包络线幅值 、sEMG 起止点检测结果

仔细观察端点检测结果出现错误的样本,发现出现这些错误时,ECG 干扰均恰好处于 sEMG 波形的边缘;如 图 6c 所示,sEMG 波形的终点附近出现一个 ECG 波形,终点检测结果受此干扰判断错误。但通常 ECG 波形的持续时间显著小于 sEMG,即使端点检测所得 sEMG 起止点位置与真实位置之间存在一定偏差,但该偏差相对于 sEMG 的持续时间而言很小,对 sEMG 的后续分析影响不大。

如算法所描述,对于一段长度为 L D 个采样点的肌电实测信号,算法实现的程序指令时间约为 T L D )≈ L D L F + L A +5)/ M +2 L D ,即其算法复杂度为 O L D ),复杂度低,可满足快速运算的要求。

3. 结论

本文借鉴语音信号端点检测的方法,建立了一种可克服 ECG 干扰影响的 sEMG 起止点检测算法。该算法根据 sEMG 与 ECG 短时过零率的差异以及 sEMG 与背景噪声短时能量的差异确定算法参数,从而实现 sEMG 的起止点检测。本文以腹直肌 sEMG 信号为典型应用,介绍了确定该算法的短时能量和短时过零率判断阈值等各参数的方法,并给出了实验结果。实验结果表明,该方法不易受 ECG 干扰而可实现 sEMG 起止点的准确检测,检测准确率达到 95.6%。

Funding Statement

广东省科技计划资助项目(2017B020210002)

References

1. Jeffrey R C, Glenn S K. Introduction to surface electromyography. Maryland: An Aspen Publication, 1998.
2. Jalil M, Butt F A, Malik A. Short-time energy, magnitude, zero crossing rate and autocorrelation measurement for discriminating voiced and unvoiced segments of speech signals//2013 International Conference on Technological Advances in Electrical, Electronics and Computer Engineering (TAEECE). Konya, Turkey: IEEE, 2013: 208-212.
3. 阳国清. 连续语音流中咳嗽信号的识别. 广州: 华南理工大学, 2010.
4. Hu Y, Li X H, Xie X B, et al. Applying independent component analysis on ECG cancellation technique for the surface recording of trunk electromyography//27th Annual International Conference of the Engineering in Medicine and Biology Society, 2005. IEEE-EMBS 2005. Shanghai, China: IEEE, 2006: 3647-3649.
5. Hu Y, Mak J, Liu H, et al. ECG cancellation for surface electromyography measurement using independent component analysis//IEEE International Symposium on Circuits and Systems. New Orleans, LA, USA: IEEE, 2007: 3235-3238.
6. Kadefors R, Broman H, de Luca C J. Units, terms and standards in the reporting of EMG research. Boston: the Ad Hoc Committee of the International Society of Electrophysiogical Kinesiology, 1980.
7. Redfern M, Hughes R, Chaffin D High-pass filtering to remove electrocardiographic interference from torso EMG recordings. Clin Biomech (Bristol, Avon) 1993; 8 (1):44–48. doi: 10.1016/S0268-0033(05)80009-9. [ PubMed ] [ CrossRef ] [ Google Scholar ]
8. Drake J D M, Callaghan J P Elimination of electrocardiogram contamination from electromyogram signals: An evaluation of currently used removal techniques. Journal of Electromyography & Kinesiology. 2006; 16 (2):175–187. [ PubMed ] [ Google Scholar ]
9. Mak J N F, Hu Y, Luk K D K ICA-based ECG removal from surface electromyography and its effect on low back pain assessment//International IEEE/EMBS Conference on Neural Engineering. Kohala Coast, HI, USA: IEEE. 2007:646–649. [ Google Scholar ]
10. Costa J D, Seixas J M, de Sa A M F L M. Independent component analysis to remove ECG interference from single channel electromyographic signal//2013 ISSNIP Biosignals and Biorobotics Conference: Biosignals and Robotics for Better and Safer Living (BRC). Rio de Janerio, Brazil: IEEE, 2013: 1-5.
11. Nougarou F, Massicotte D, Descarreaux M. Efficient combination of DWT and ICA to localize and remove ECG from surface electromyography measurement//201318th International Conference on Digital Signal Processing (DSP). Fira, Greece: IEEE, 2013: 1-6.
12. 周权, 杨智, 范正平, 等 结合小波变换和能量熵在膈肌肌电中去除心电干扰的应用 生物医学工程学杂志 2013; 30 (1):16–21. [ Google Scholar ]
13. Sano T, Shiozawa T, Uenosono N, et al Recommendations for standardization of leads and of specifications for instruments in electrocardiography and vectorcardiography (report of Committee on Electrocardiography, American Heart Association) Journal of Japanese Medical Instruments. 1970; 40 :242–258. [ Google Scholar ]
14. Merletti R Standards for Reporting EMG Data. Journal of Electromyography & Kinesiology. 2014; 24 (2):I–II. [ Google Scholar ]

Articles from Sheng Wu Yi Xue Gong Cheng Xue Za Zhi = Journal of Biomedical Engineering are provided here courtesy of West China Hospital of Sichuan University