QRS波群检测方法、心电检测装置及可读存储介质与流程

qrs波群检测方法、心电检测装置及可读存储介质
技术领域
1.本发明涉及一种医疗电子领域，尤其涉及一种qrs波群检测方法、心电检测装置及可读存储介质。


背景技术：

2.心电图(ecg:electrocardiograph)指心脏在每个心动周期中，由起搏点、心房、心室相继兴奋，伴随着生物电的变化，通过心电描记器从体表引出多种形式的电位变化的图形。通常使用多个电极片采集身体多个部位的电势差，再通过ad转换芯片生成连续信号。心电图是反映心脏兴奋的电活动过程，它对心脏基本功能及其病理研究方面，具有重要的参考价值。参阅图1所示，一个典型的ecg信号包括一个p波，一个qrs波群和一个t波，精确测量中可能还存在一个u波。因为r波普遍较为突出，所以在心电信号分析中，qrs波群的检测常作为关键节点，通常是后续各种生理分析的基石。
3.由于工频噪声、肌电信号、芯片底噪等基本都是全频段干扰信号，传统滤波方法如低通、高通、带通滤波器等均难以完全滤除。此外，心电信号中t波幅值也较大，在部分人群中甚至可能高于r波，容易造成误检。一些新兴的研究方法如神经网络算法、模版匹配算法、troika算法、隐马尔科夫模型和希尔伯特-黄变换等，由于需要整段信号同时处理，且计算量巨大，通常用于服务器端后台脱机处理，不适用于嵌入式等常用简易设备的实时分析。


技术实现要素：

4.鉴于以上内容，有必要提供一种qrs波群检测方法、心电检测装置及可读存储介质，可有效防止漏检及误检qrs波群，提高检测的灵敏度和准确性。
5.所述qrs波群检测方法包括：采集心电信号；利用预设的滤波器对所采集的心电信号进行滤波处理，获得经过滤波处理的心电信号；采用双斜率法处理所述经过滤波处理的心电信号，获得经过双斜率法处理的心电信号；从经过双斜率法处理的心电信号中检测r波峰；及基于所检测的r波峰输出qrs波群的位置。
6.优选地，所述预设的滤波器包括工频陷波器和二阶iir高通滤波器。
7.优选地，所述采用双斜率法处理所述经过滤波处理的心电信号包括：计算所述经过滤波处理的心电信号的任意一个信号采样点的左侧的最大平均斜率以及该任意一个信号采样点的右侧的最大平均斜率；将所述任意一个信号采样点的左侧的最大平均斜率和该任意一个信号采样点的右侧的最大平均斜率中的较大者作为该任意一个信号采样点的最大平均斜率。
8.优选地，所述任意一个信号采样点的左侧的最大平均斜率为s
l,max
，该任意一个信号采样点的右侧的最大平均斜率为s
r,max
；其中，
9.其中，n、m、k为预设的正整数，z
i
代表信
号采样点i的电压值，z
i-k
代表自该信号采样点i向左的第k个信号采样点的电压值，z
i+k
代表自该信号采样点i向右的第k个信号采样点的电压值。
10.优选地，所述从经过双斜率法处理的心电信号中检测r波峰包括：对每个信号采样点是否为波峰进行确认；及根据相邻两个波峰之间的间距确定该相邻两个波峰中的每个波峰是否为r波峰。
11.优选地，所述对每个信号采样点是否为波峰进行确认包括：当任意一个信号采样点的电压值大于第一阈值thr1且该任意一个信号采样点的电压值为极大值时，确定该任意一个信号采样点为波峰；所述根据相邻两个波峰之间的间距确定该相邻两个波峰中的每个波峰是否为r波峰包括：当相邻两个波峰之间的信号采样点的总数小于或者等于预设值时，确定该相邻两个波峰中对应更大的电压值的波峰为r波峰；及当相邻两个波峰之间的信号采样点的总数大于所述预设值时，确定该相邻两个波峰均为r波峰。
12.优选地，该方法还包括：基于下述公式根据当前所确定的所有r波峰的峰值中的最大峰值peak调整所述第一阈值thr1，以及调整第二阈值thr2，所述第二阈值thr2用于确定该第一阈值thr1；其中，
[0013][0014]
所述可读存储介质存储有至少一个指令，所述至少一个指令被处理器执行时实现所述qrs波群检测方法。
[0015]
所述心电检测装置包括：存储器；至少一个处理器；所述存储器中存储有多个模块，所述多个模块由所述至少一个处理器执行，该多个模块包括：采集模块，用于采集心电信号；执行模块，用于利用预设的滤波器对所采集的心电信号进行滤波处理，获得经过滤波处理的心电信号；所述执行模块，还用于采用双斜率法处理所述经过滤波处理的心电信号，获得经过双斜率法处理的心电信号；所述执行模块，还用于从经过双斜率法处理的心电信号中检测r波峰；及所述执行模块，还用于基于所检测的r波峰输出qrs波群的位置。
[0016]
相较于现有技术，所述qrs波群检测方法、心电检测装置及可读存储介质，计算复杂度较低，计算量较小，可有效防止漏检及误检qrs波群，提高检测的灵敏度和准确性。
附图说明
[0017]
图1示意了一个ecg信号。
[0018]
图2是本发明较佳实施例的心电检测装置的架构图。
[0019]
图3是本发明较佳实施例的检测系统的功能模块图。
[0020]
图4是本发明较佳实施例的qrs波群检测方法的流程图。
[0021]
图5a为利用滤波器处理过的心电信号的示意图。
[0022]
图5b为利用双斜率法对经过滤波处理的心电信号作了进一步的处理后的示意图。
[0023]
主要元件符号说明
[0024]
心电检测装置3
处理器32存储器31检测系统30采集模块301执行模块302
[0025]
如下具体实施方式将结合上述附图进一步说明本发明。
具体实施方式
[0026]
为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施例对本发明进行详细描述。需要说明的是，在不冲突的情况下，本发明的实施例及实施例中的特征可以相互组合。
[0027]
在下面的描述中阐述了很多具体细节以便于充分理解本发明，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0028]
除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。
[0029]
参阅图2所示，为本发明较佳实施例提供的心电检测装置的架构图。
[0030]
本实施例中，所述心电检测装置3包括互相之间电气连接的存储器31、至少一个处理器32、至少一条通信总线33。
[0031]
本领域技术人员应该了解，图2示出的心电检测装置3的结构并不构成本发明实施例的限定，所述心电检测装置3还可以包括比图2更多或更少的其他硬件或者软件，或者不同的部件布置。
[0032]
需要说明的是，所述心电检测装置3仅为举例，其他现有的或今后可能出现的电子产品如可适应于本发明，也应包含在本发明的保护范围以内，并以引用方式包含于此。
[0033]
在一些实施例中，所述存储器31可以用于存储计算机程序的程序代码和各种数据。例如，所述存储器31可以用于存储安装在所述心电检测装置3中的检测系统30，并在心电检测装置3的运行过程中实现高速、自动地完成程序或数据的存取。所述存储器31可以是包括只读存储器(read-only memory，rom)、可编程只读存储器(programmable read-only memory，prom)、可擦除可编程只读存储器(erasable programmable read-only memory，eprom)、一次可编程只读存储器(one-time programmable read-only memory，otprom)、电子擦除式可复写只读存储器(electrically-erasable programmable read-only memory，eeprom)、只读光盘(compact disc read-only memory，cd-rom)或其他光盘存储器、磁盘存储器、磁带存储器、或者任何其他能够用于携带或存储数据的非易失性的计算机可读的存储介质。
[0034]
在一些实施例中，所述至少一个处理器32可以由集成电路组成。例如，可以由单个封装的集成电路所组成，也可以是由多个相同功能或不同功能封装的集成电路所组成，包括一个或者多个中央处理器(central processing unit，cpu)、微处理器、数字处理芯片、图形处理器及各种控制芯片的组合等。所述至少一个处理器32是所述心电检测装置3的控
制核心(control unit)，利用各种接口和线路连接整个心电检测装置3的各个部件，通过执行存储在所述存储器31内的程序或者模块或者指令，以及调用存储在所述存储器31内的数据，以执行心电检测装置3的各种功能和处理数据，例如，对心电信号的qrs波群进行检测的功能(具体细节参后面对图4的介绍)。
[0035]
在一些实施例中，所述至少一条通信总线33被设置为实现所述存储器31与所述至少一个处理器32以及其他元件之间的连接通信。
[0036]
在本实施例中，检测系统30可以包括一个或多个模块，所述一个或多个模块存储在所述存储器31中，并由至少一个或多个处理器(本实施例为处理器32)执行，以实现对心电信号的qrs波群进行检测的功能(具体细节参后面对图4的介绍)。
[0037]
在本实施例中，参阅图3所示，所述检测系统30根据其所执行的功能，可以被划分为多个功能模块。所述多个功能模块包括采集模块301和执行模块302。本发明所称的模块是指一种能够被至少一个处理器(例如处理器32)所执行并且能够完成固定功能的一系列计算机可读的指令段，其存储在存储器(例如心电检测装置3的存储器31)中。在本实施例中，关于各模块的功能将在后续结合图4中详述。
[0038]
本实施例中，以软件功能模块的形式实现的集成的单元，可以存储在一个非易失性可读取存储介质中。上述软件功能模块包括若干计算机可读指令用以使得一台心电检测装置3或处理器(processor)执行本发明各个实施例的方法的部分，例如图4所示的对心电信号的qrs波群进行检测的方法。
[0039]
在进一步的实施例中，结合图3，所述至少一个处理器32可执行所述心电检测装置3的操作装置以及安装的各类应用程序(如所述的检测系统30)、程序代码等，例如，上述的各个模块。
[0040]
所述存储器31中存储有程序代码，且所述至少一个处理器32可调用所述存储器31中存储的程序代码以执行相关的功能。例如，图3中所述检测系统30的各个模块是存储在所述存储器31中的程序代码，并由所述至少一个处理器32所执行，从而实现所述各个模块的功能以达到对心电信号的qrs波群进行检测的目的(详见下文中对图4的描述)。
[0041]
在本发明的一个实施例中，所述存储器31存储多个计算机可读指令，所述多个计算机可读指令被所述至少一个处理器32所执行以实现对心电信号的qrs波群进行检测的目的。具体地，所述至少一个处理器32对上述计算机可读指令的具体实现方法详见下文中对图4的描述。
[0042]
需要说明的是，在其他实施例中，所述检测系统30也可以与存储器、处理器及其他必要的硬件或软件实现为一个嵌入式系统。
[0043]
图4是本发明较佳实施例提供的检测方法的流程图。
[0044]
在本实施例中，所述检测方法可以应用于心电检测装置3中，对于需要进行心电信号的qrs波群进行检测的心电检测装置3，可以直接在心电检测装置3上集成本发明的方法所提供的用于检测的功能，或者以软件开发工具包(software development kit，sdk)的形式运行在心电检测装置3上。
[0045]
如图4所示，所述检测方法具体包括以下步骤，根据不同的需求，该流程图中步骤的顺序可以改变，某些步骤可以省略。
[0046]
步骤s1、采集模块301采集心电信号。
[0047]
具体地，采集模块301可利用粘着在人体皮肤上的在心外部的区域中的电极来采集心电信号。
[0048]
步骤s2、执行模块302利用预设的滤波器对所采集的心电信号进行滤波处理，获得经过滤波处理的心电信号。
[0049]
本实施例中，所述预设的滤波器包括工频陷波器和二阶iir高通滤波器。
[0050]
在一个实施例中，所述工频陷波器的频率可以为50hz。所述执行模块302利用所述工频陷波器可以从所述心电信号中去除工频噪声干扰。在一个实施例中，所述二阶iir高通滤波器的截止频率可以为0.4hz，所述执行模块302利用所述二阶iir高通滤波器可以从所述心电信号中去除基线漂移。
[0051]
在其他实施例中，所述预设的滤波器也可以为其他型号的滤波器。
[0052]
步骤s3、执行模块302采用双斜率法处理所述经过滤波处理的心电信号，获得经过双斜率法处理的心电信号。
[0053]
在一个实施例中，所述采用双斜率法处理所述经过滤波处理的心电信号包括(a1)-(a2)：
[0054]
(a1)计算所述经过滤波处理的心电信号的任意一个信号采样点的左侧的最大平均斜率以及该任意一个信号采样点的右侧的最大平均斜率。
[0055]
在一个实施例中，所述任意一个信号采样点的左侧的最大平均斜率为s
l,max
，该任意一个信号采样点的右侧的最大平均斜率为s
r,max
。
[0056]
其中，
[0057]
其中，n、m、k为预设的正整数，z
i
代表信号采样点i的电压值，z
i-k
代表自该信号采样点i向左的第k个信号采样点的电压值，z
i+k
代表自该信号采样点i向右的第k个信号采样点的电压值。
[0058]
在一个实施例中，n和m的值根据所述心电信号的采样频率以及qrs波群的通常宽度来设置。
[0059]
举例而言，假设所述心电信号的采样频率为250hz，qrs波群的宽度为0.06-0.1s，那么n可以等于6，m可以等于15。
[0060]
(a2)将所述任意一个信号采样点的左侧的最大平均斜率和该任意一个信号采样点的右侧的最大平均斜率中的较大者作为该任意一个信号采样点的最大平均斜率。
[0061]
需要说明的是，在其他实施例中，当所述任意一个信号采样点的左侧的最大平均斜率等于该任意一个信号采样点的右侧的最大平均斜率时，则将所述任意一个信号采样点的左侧的最大平均斜率或者该任意一个信号采样点的右侧的最大平均斜率作为该任意一个信号采样点的最大平均斜率。
[0062]
采用本发明中的上述双斜率法对所述心电信号处理后，r波峰明显突出，更易于后续qrs波群峰的检测。效果可参图5a-图5b所示。
[0063]
具体地，图5a为利用所述预设的滤波器处理过的心电信号的示意图。从图5a中可以看出，经过滤波处理后心电信号的r波峰水平约为6800，噪声水平约为4300，噪声幅值比例约为63.24％。图5b为利用所述双斜率法对所述经过滤波处理的心电信号作进一步的处
理后的示意图。从图5b可以看出，经过双斜率法处理后的心电信号的r波峰水平约为1400，噪声水平约为400，噪声幅值比例约为28.57％。很明显，经过所述双斜率法处理后，心电信号中的t波及噪声干扰明显降低。此外，从图5a中可以看出横坐标8000-8500点处有较强干扰，从图5b中可以看出，经过双斜率法处理后横坐标8000-8500点处的干扰得到了完美滤除。
[0064]
步骤s4、执行模块302从经过双斜率法处理的心电信号中检测r波峰。
[0065]
本实施例中，所述从经过双斜率法处理的心电信号中检测r波峰包括(b1)-(b2)：
[0066]
(b1)对每个信号采样点是否为波峰进行确认。
[0067]
本实施例中，所述对每个信号采样点是否为波峰进行确认包括：当任意一个信号采样点的电压值大于第一阈值thr1且该任意一个信号采样点的电压值为极大值时，确定该任意一个信号采样点为波峰；及当任意一个信号采样点的电压值小于或者等于所述第一阈值thr1，及/或该任意一个信号采样点的电压值不为极大值时，确定该任意一个信号采样点不为波峰。
[0068]
本实施例中，任意一个信号采样点的电压值为极大值是指该任意一个信号采样点的电压值大于自该任意一个信号采样点向左的第一个信号采样点的电压值，且该任意一个信号采样点的电压值大于自该任意一个信号采样点向右的第一个信号采样点的电压值。
[0069]
(b2)根据相邻两个波峰之间的间距确定该相邻两个波峰中的每个波峰是否为r波峰。
[0070]
本实施例中，所述根据相邻两个波峰之间的间距确定该相邻两个波峰中的每个波峰是否为r波峰包括：
[0071]
当相邻两个波峰之间的信号采样点的总数小于或者等于预设值(例如60)时，确定该相邻两个波峰中对应更大的电压值的波峰为r波峰；及
[0072]
当相邻两个波峰之间的信号采样点的总数大于所述预设值时，确定该相邻两个波峰均为r波峰。
[0073]
步骤s5、执行模块302基于所检测的r波峰输出qrs波群的位置。
[0074]
具体地，qrs波群的位置可以用r波峰的位置来表示。r波峰的位置则是用对应r波峰的信号采样点所对应的时间来表示。
[0075]
在一个实施例中，执行模块302输出qrs波群的位置时还可以考虑滤波延时，即将r波峰的位置往后移动预设时间点后的位置作为qrs波群的位置。
[0076]
具体地，所述预设时间点可以根据所述预设的滤波器的阶数来确定。
[0077]
在其他实施例中，所述执行模块302于输出所述qrs波群的位置之前，还可以基于下述公式根据当前所确定的所有r波峰的峰值中的最大峰值peak调整所述第一阈值thr1，以及调整第二阈值thr2，所述第二阈值thr2用于确定该第一阈值thr1。具体地，
[0078][0079]
[0080]
需要说明的是，所述执行模块302根据当前所确定的所有r波峰的峰值中的最大峰值peak实时调整所述第一阈值thr1和第二阈值thr2，使得本发明可以实时适应不断变化的信号水平，进而可有效地防止漏检及误检r波峰，提高检测的灵敏度和准确性。
[0081]
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照以上较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或等同替换，而不脱离本发明技术方案的精神和范围。