申请/专利权人：北京康博众联电子科技有限公司

申请日：2017-10-27

公开（公告）日：2020-06-26

公开（公告）号：CN107997759B

主分类号：A61B5/0402(20060101)

分类号：A61B5/0402(20060101);A61B5/0452(20060101);A61B5/0456(20060101)

优先权：

专利状态码：有效-授权

法律状态：2020.06.26#授权;2018.06.01#实质审查的生效;2018.05.08#公开

摘要：本发明公开了一种心电信号中QRS波群的检测方法，对于从心电信号中提取的极值点集合Ex‑a，如果一组连续或非连续极值具有相似的高度和两边斜率，且这组极值中相邻的极值间隔基本相同，则称之为一个相似极值序列；对于每一个相似极值序列，在这个序列的覆盖的范围内，去除所有不属于这个极值序列的极值；这样处理完每一个相似极值序列，形成优选极值点集合Ex‑b；针对优选极值点集合标记QRS波群。本发明还公开了一种承载上述方法的计算机可读存储介质和计算机设备。使用本发明能够自动检测心电信号中的QRS波群位置，速度快，准确率高。

主权项：1.一种心电信号中QRS波群的检测方法，其特征在于，包括：对于从心电信号中提取的极值点集合Ex-a，如果一组连续或非连续极值具有相似的高度和两边斜率，且这组极值中相邻的极值间隔基本相同，则称之为一个相似极值序列；对于每一个相似极值序列，在这个序列的覆盖的范围内，去除所有不属于这个极值序列的极值；按同样的方式处理极值点集合Ex-a中的每一个相似极值序列，形成优选极值点集合Ex-b；利用优选极值点集合标记QRS波群；进行相似极值序列的处理之前，进一步去除小短边极值，在相似极值序列的处理之后，进一步去除小长边极值；所述去除小短边极值为：对于每个极值点，其短边幅值为exs，该极值点左右两个相邻极值点的短边幅值分别为exs1和exs2；如果exs小于exs1和exs2，且exs与exs1和exs2的长度差均超过设定比例，则认为该极值点短边过小，从极值点集合中去除；所述去除小长边极值为：对于每个极值点，其长边幅值为exl，该极值点左右相邻两个极值点的长边幅值分别为exl1和exl2；如果exl小于exl1和exl2，且exl与exl1和exl2的长度差均超过设定比例，则认为该极值点长边过小，从极值点集合中去除。

全文数据：心电信号中QRS波群的检测方法、存储介质和计算机设备技术领域[0001]本发明属于心电信号的信号处理技术领域，具体涉及一种心电信号中QRS波群的检测方法、一种计算机可读存储介质和一种计算机设备。背景技术[0002]心电信号的自动分析与处理一直是生物医学工程领域中的重要研究课题，其中，心电信号的自动分析与处理关键问题是QRS波群的检测，因为可靠的QRS波群检测不仅是诊断心律失常的重要依据，而且也只有在QRS波群确定之后才有可能计算心率并进行心率变异分析和分析心电信号的其它细节信息。[0003]QRS波群检测学术界有很多方法。但是对于实际检测，会发现心电信号具有特有的噪声特点。如图1a的心电信号，其特点是噪声密集、噪声幅度与信号幅度相当，噪声区内QRS波群被覆盖，对于这种情况不应该检测噪声区的QRS波群;如图1⑹所示的心电信号，其特点为小幅度密集噪声，噪声区内QRS波群仍然可见，对于这种情况应该设法检测出噪声区的QRS波群，避免噪声区的QRS误检;如图1c所示的心电信号，其特点为单点噪声，应该避免QRS误检。而对于图1d所示的心电信号，其特点为稀疏大幅度噪声，对于这种情况也不应该检测噪声段内的QRS波群。[0004]现有的QRS波群检测没有针对上述情况进行具有针对性的设计，对各种形态的噪声均具有较高的敏感度，尤其针对图1⑹和图1c的情况，现有算法会产生误报，也就是将噪声检测为QRS波群，导致QRS波群检测率不是很高。发明内容[0005]有鉴于此，本发明提供了一种心电信号中QRS波群的检测方法，能够提高QRS波群检测的准确率，而且速度快。[0006]为解决上述技术问题，本发明具体方法如下：[0007]一种心电信号中QRS波群的检测方法，包括：[0008]对于从心电信号中提取的极值点集合Ex-a，如果一组连续或非连续极值具有相似的高度和两边斜率，且这组极值中相邻的极值间隔基本相同，则称之为一个相似极值序列；对于每一个相似极值序列，在这个序列的覆盖的范围内，去除所有不属于这个极值序列的极值;按同样的方式处理极值点集合Ex-a中的每一个相似极值序列，形成优选极值点集合Ex-b;利用优选极值点集合标记QRS波群。^[0009]优选地，进行相似极值序列的处理之前，进一步去除小短边极值，在相似极值序列的处理之后，进一步去除小长边极值；[0010]所述去除小短边极值为:对于每个极值点，其短边幅值为exs，该极值点左右两个相邻极值点的短边幅值分别为exsl和exs2;如果exs小于exsl和exs2,且exs与exsl和exs2的长度差均超过设定比例，则认为该极值点短边过小，从极值点集合中去除；[0011]所述去除小长边极值为:对于每个极值点，其长边幅值为exl，该极值点左右相邻两个极值点的长边幅值分别为exll和exl2;如果exl小于exll和exl2,且exl与exll和exl2的长度差均超过设定比例，则认为该极值点长边过小，从极值点集合中去除。[0012]优选地，所述exs与exsl和exs2的长度差均超过设定比例的确定方式为:如果exsexslX0•4并且exsexs2X0•4，则认为这个极值点短边过小；[0013]所述exl与exll和exl2的长度差均超过设定比例的确定方式为:如果exlexllX0.4并且exl〈exl2X0•4,则认为这个极值点长边过小。[0014]优选地，所述去除小短边极值之后，相似极值序列处理之前，进一步从极值点集合中去除间距小于设定小值的一对极值中短边小的一个极值点。[0015]优选地，所述设定小值为0.17s。[0016]优选地，提取所述极值点集合Ex-a的方式为:获取去除基线信号的心电信号的极值点，组成集合Exl;从集合Exl中去除不属于QRS波群的极值，得到集合Ex-a。[0017]优选地，所述去除基线信号后的心电信号的获取方式为:采用对心电信号F进行分段中值滤波的方式得到基线信号B，进而获得所述去除基线信号B后的心电信号。[0018]优选地，所述去除不属于QRS波群的极值为:在心电信号F中找出所有时间长度为0.4的片段，且片段内峰值大于1毫伏的极值数量达到8个，则在极值点集合中去除这些片段内的极值；[0019]并且，从极值点集合中去除斜率过小的极值点。[0020]本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现上述任意一种检测方法。[0021]本发明还提供了一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述任意一种检测方法。[00221有益效果[0G23]⑴针对图1a〜图1⑹出现的各种情况，现有算法对这些噪声没有特别处理，所以只能处理其中一部分情况。尤其是图2和图3,现有算法会产生误报，也就是将噪声检测为QRS波群。本发明的“去除相似极值序列内的其他极值”的操作，正是针对图1⑹和图1c的情况设计的步骤，可以去除覆盖在QRS波群中的噪声，从而去除噪声带来的误报，提高QRS波群检测率。[0024]⑵通过去除小短边极值、去除临近小极值和去除小长边极值，从而去除不可能为心跳的噪声。其中，将去除小长边极值的步骤放在“去除相似极值序列内的其他极值”之后，可以避免误删除，误删除会导致寻找相似极值序列时失败。[0025]3通过去除噪声段内的极值、去除小斜率极值的操作，去除确定为非QRS波群的部分，提高了后续操作的准确度。[0026]⑷该方法简单，易实现，不采用复杂的公式，因此检测速度快。附图说明[0027]图1a〜图id为四种噪声的示意图。[0028]图2为本发明流程图。[0029]图3为去除小斜率极值的示意图。[0030]图4为去除小短边极值的示意图。[0031]图5为相似极值序列和夹杂在其中的噪声示意图。[0032]图6为标记QRS波群的示意图。[0033]图7为用来实现本申请实施例的检测方法的计算机设备5〇〇的结构示意图。具体实施方式[0034]下面将参照附图来说明本发明的实施例。[0035]本发明提供了一种心电信号中QRS波群的检测方法，其基本思想是:对于从心电信号中提取的极值点集合Ex-a，如果一组连续或非连续极值具有相似的高度和两边斜率，且这组极值中相邻的极值间隔基本相同，则称之为一个相似极值序列;对于每一个相似极值序列，在这个序列的覆盖的范围内，去除所有不属于这个极值序列的极值;这样处理完每一个相似极值序列，形成优选极值点集合Ex-b;最后，根据优选极值点集合标记QRS波群。[0036]本发明“去除相似极值序列内的其他极值”的操作，正是针对图1⑹和图1c的情况设计的步骤，可以去除覆盖在QRS波群中的噪声，从而去除噪声带来的误报，提高QRS波群检测率。[0037]基于上述基本思想，为了进一步降低噪声，还设计了多个噪声去除步骤。结合噪声去除步骤，本发明实施例的具体流程参见图2,该方法包括如下步骤：[0038]步骤一、去除基线。[0039]本步骤对心电信号F进行分段中值滤波，得到基线信号B，则去除基线信号B后的心电信号F’=F-B。[0040]步骤二、检测极值。[0041]本步骤检测心电信号F的所有极值点，包括最小值和最大值，将所有极值点组成的集合记为Ex。[0042]步骤三、去除噪声段内的极值。[0043]本步骤中，在心电信号F中，找出所有时间长度为〇.4秒的片段，且片段内峰值大于1毫伏的极值数量达到8个。在集合Ex中去除这些片段内的极值，得到集合Ex1。本步骤的目的是去除图1a中的幅度与信号相当的密集噪声。[0044]步骤四、去除小斜率极值。[0045]本步骤中，在集合Exl中去除斜率过小小于设定值的极值点，得到集合Ex2。[0046]如图3所示，箭头处是一个极值点，其左边设定区域本实施例为0.2秒区域）内斜率最大值定义为左边斜率，右边斜率同样可得。对于一个极值点，如果这两个斜率中至少有一个小于设定值本实施例中为2.7毫伏每秒），则去除该极值点。箭头处的极值点就是一个斜率过于小的例子，将被去掉。这是一个T波，不是QRS波群。该点左边那个尖锐的波形是QRS波群。[0047]步骤五、去除小短边极值。[0048]每个极值点两边的高斜率的片段称为“边”，其中将极值点左右两个边中的较长者称为长边，较短者称为短边。本步骤在集合2中去除短边小的极值，形成集合Ex3。[0049]本步骤中，判断小短边的标准是这样的:如图4所示，对于一个极值点ex，其短边幅值为exs〇ex左右相邻的两个极值点为exl，ex2,其短边幅值分别为exsl，exs2。如果exs小于exsl和exs2,且exs与exsl和exs2的长度差均超过设定比例，则认为这个极值点短边过小，会被去除。[0050]本实施例中，exs与exsl和exs2的长度差均超过设定比例的判断标准为:exsexslX0•4并且exsexs2X0.4。[0051]通常情况下，QRS波群的短边不会太小。所以短边小说明QRS波群形态不大符合标准。长边小则说明QRS波群整体幅度小。那么短边小可以确信是噪声，而长边小大部分情况是噪声，还有很小的可能性是幅度小的QRS波群。这里先去除短边小的，不去除长边小的，以避免误删除。误删除会导致步骤七在寻找相似极值序列时失败。所以去除长边小的极值的步骤放在步骤七之后。[0052]步骤六、去除临近小极值。[0053]在集合Ex3中，相邻很近的两个极值，去除一个高度小的，形成集合Ex4。[0054]通过上述的步骤以后，还会存在0.17秒以内有两个极值的情况，而心跳是不会相邻这么近的，所以去除一个。其中，去除短边较小的。[0055]步骤七、去除相似极值序列内的其他极值。[0056]如果一组极值，这组极值不必是连续的，但都具有相似的高度和两边斜率，且这组极值中相邻极值间隔基本相同，则称之为一个相似极值序列。在集合Ex4中，对于每一个相似极值序列，在这个序列的覆盖的范围内，去除所有不属于这个极值序列的极值。这样处理完每一个相似极值序列，就形成集合Ex5。[0057]这里所述的相似和基本相同可以通过判断差值是否小于设定阈值为来判定。如图5所示，下部箭头的极值点构成了相似极值序列。那么在这个序列覆盖的这段范围内，还存在两个极值，即上部箭头处。这两个极值将被当作噪声去除。[0058]步骤八、去除小长边极值。[0059]本步骤可参考步骤五的“去除小短边极值”。在集合Ex5中去除长边小的极值，形成集合Ex6。[0060]所述去除小长边极值为:对于每个极值点，其长边幅值为exl，该极值点左右相邻两个极值点的长边幅值分别为exll和exl2;exl小于exll和exl2,且exl与exll和exl2的长度差均超过设定比例，则认为该极值点长边过小，从极值点集合中去除。本实施例中，所述exl与exll和exl2的长度差均超过设定比例的确定方式为:如果exKexllXO.4并且exlex12X0.4，则认为这个极值点长边过小。[0061]步骤九、标记QRS波群。[0062]在集合Ex6中，每一个极值记录为一个R波，极值左边的高斜率片段的端点为Q波，极值右边的高斜率片段的端点为S波。所谓高斜率片段就是指的前面提到的长边和短边。如图6所示箭头所指分别为Q、R、S点，QR段和RS段就是所谓的高斜率片段。[0063]至此，本流程结束。[0064]参考图7,其示出了适于用来实现本申请实施例的检测方法的计算机设备5〇〇的结构示意图。[0065]如图7所示，计算机设备500包括中央处理单元CPU5〇1，其可以根据存储在只读存储器ROM502中的程序或者从存储部分508加载到随机访问存储器RAM5〇3中的程序而执行各种适当的动作和处理。在RAM503中，还存储有系统500操作所需的各种程序和数据。CPU501、ROM502以及RAM503通过总线504彼此相连。输入输出（IO接口5〇5也连接至总线504。^6]以:部件^接至10接口505:包括键盘、鼠标等的输入部分5〇6;包括诸如阴极射，官c=、液晶显示器LCD等以及扬声器等的输出部分5〇7;包括硬盘等的存储部分5〇8;以及包括诸如LAN卡、调制解调器等的网络接口卡的通信部分s〇9。通信部分5〇9经由诸如因特网的网络执行通信处理。驱动器510也根据需要连接至1〇接口5〇5。可拆卸介质5U，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器51〇上，以便于从其上读出的计算机程序根据需要被安装入存储部分508。该计算机程序实现本申请的QRS波群检测方法。[0067]本申请还提供了一种计算机可读存储介质，该计算机可读存储介质可以是上述方案所述的计算机设备中所包含的计算机可读存储介质；也可以是单独存在，未装配入设备中的计算机可读存储介质。计算机可读存储介质存储有一个或者一个以上程序，所述程序被一个或者一个以上的处理器用来执行描述于本申请的QRS波群检测方法。[0068]综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

权利要求：1.一种心电信号中QRS波群的检测方法，其特征在于，包括：对于从心电信号中提取的极值点集合Ex-a，如果一组连续或非连续极值具有相似的高度和两边斜率，且这组极值中相邻的极值间隔基本相同，则称之为一个相似极值序列;对于每一个相似极值序列，在这个序列的覆盖的范围内，去除所有不属于这个极值序列的极值；按同样的方式处理极值点集合Ex-a中的每一个相似极值序列，形成优选极值点集合Ex-b;利用优选极值点集合标记QRS波群。2.如权利要求1所述的方法，其特征在于，进行相似极值序列的处理之前，进一步去除小短边极值，在相似极值序列的处理之后，进一步去除小长边极值；所述去除小短边极值为:对于每个极值点，其短边幅值为exs，该极值点左右两个相邻极值点的短边幅值分别为exsl和exs2;如果exs小于exsl和exs2,且exs与exsl和exs2的长度差均超过设定比例，则认为该极值点短边过小，从极值点集合中去除；所述去除小长边极值为:对于每个极值点，其长边幅值为exl，该极值点左右相邻两个极值点的长边幅值分别为exll和exl2;如果exl小于exll和exl2,且exl与exll和exl2的长度差均超过设定比例，则认为该极值点长边过小，从极值点集合中去除。3.如权利要求2所述的方法，其特征在于，所述exs与exsl和exs2的长度差均超过设定比例的确定方式为:如果exsexslX0.4并且exsexs2X0.4，则认为这个极值点短边过小；所述exl与exll和exl2的长度差均超过设定比例的确定方式为:如果exlexlIX0.4并且exlexl2X0.4,则认为这个极值点长边过小。4.如权利要求2所述的方法，其特征在于，所述去除小短边极值之后，相似极值序列处理之前，进一步从极值点集合中去除间距小于设定小值的一对极值中短边小的一个极值点。5.如权利要求4所述的方法，其特征在于，所述设定小值为0.17s。6.如权利要求1所述的方法，其特征在于，提取所述极值点集合Ex-a的方式为：获取去除基线信号的心电信号的极值点，组成集合Exl;从集合Exl中去除不属于QRS波群的极值，得到集合Ex-a。7.如权利要求6所述的方法，其特征在于，所述去除基线信号后的心电信号的获取方式为:采用对心电信号F进行分段中值滤波的方式得到基线信号B，进而获得所述去除基线信号B后的心电信号。8.如权利要求6所述的方法，其特征在于，所述去除不属于QRS波群的极值为:在心电信号F中找出所有时间长度为0.4的片段，且片段内峰值大于1毫伏的极值数量达到8个，则在极值点集合中去除这些片段内的极值；并且，从极值点集合中去除斜率过小的极值点。9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1〜8任意一项所述的方法。10.—种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求丨〜8任意一项所述的方法。

百度查询： 北京康博众联电子科技有限公司 心电信号中QRS波群的检测方法、存储介质和计算机设备

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD