心电图信号特征提取方法、装置、设备、系统和存储介质
阅读：368发布：2020-05-11
专利汇可以提供心电图信号特征提取方法、装置、设备、系统和存储介质专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种 心电图 信号 特征提取方法、装置、设备、系统和存储介质，其中方法包括如下步骤：提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；在平滑心电图中，提取P波；在平滑心电图中，提取T波；依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断 阈值 的±N％ 波动 范围内，则将该生理参数标记为 临界状态 ，N的取值范围为10～15。本发明提供的心电图信号特征提取方法、装置、设备、系统和存储介质，用于甄别临界状态的心电图信号，充分解读心电图中的信息。，下面是心电图信号特征提取方法、装置、设备、系统和存储介质专利的具体信息内容。
1.一种心电图信号特征提取方法，其特征在于，包括如下步骤：
提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；
在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
2.如权利要求1所述的心电图信号特征提取方法，其特征在于，利用小波分解提取心电图中QRS波群中的R波，在[tR-100ms，tR+100ms]范围内，tR时间之前第一个极小值数据点为Q波，tR时间之后第一个极小值数据点为S波。
3.如权利要求1所述的心电图信号特征提取方法，其特征在于，所述预定阈值为平滑心电图中最大值数据点电位幅值的0.4～0.6倍。
4.如权利要求1所述的心电图信号特征提取方法，其特征在于，所述Dmin的取值范围为
0.13～0.2，所述Dmax的取值范围为0.5～0.9。
5.一种心电图信号特征提取装置，其特征在于，包括：
Q波、R波和S波提取模块，用于提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
P波提取模块，用于在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；
T波提取模块，用于在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
标记模块，用于依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为
10～15。
6.一种计算机设备，包括：存储器和处理器，所述存储器内存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时，实现如权利要求1～4任一项所述的心电图信号特征提取方法。
7.一种心电图信号特征提取系统，包括终端以及服务器，所述服务器包括存储器和处理器，所述存储器内存储有计算机程序，其特征在于，所述服务器从终端获取心电图；所述处理器执行所述计算机程序时，实现如权利要求1～4任一项所述的心电图信号特征提取方法。
8.一种计算机可读存储介质，所述计算机可读存储介质内存储有计算机程序，其特征在于，所述计算机程序被计算机处理器执行时实现如权利要求1～4任一项所述的心电图信号特征提取方法。
心电图信号特征提取方法、装置、设备、系统和存储介质
技术领域
背景技术
[0002] 心电图表征了人体的心脏电活动，由于心脏疾病种类众多，同一种疾病在不同阶段表现在心电图上的特征也往往不同，对于临床中通常使用的十二导联的心电图而言，存在U波、P波、T波、以及QRS波群，根据各特征波计算得到的心脏参数存在临界状态，现有技术中没有有效方法对临界状态进行甄别，使心电图中的信息并未充分利用，可能发生遗漏。
发明内容
[0003] 本发明提供了一种心电图信号特征提取方法、装置、设备、系统和存储介质，用于甄别临界状态的心电图信号，充分解读心电图中的信息。
[0004] 一种心电图信号特征提取方法，包括如下步骤：
[0005] 提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
[0006] 在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
[0007] 若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
[0008] 若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
[0009] 若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
[0010] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
[0011] 依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
[0012] 以下还提供了若干可选方式，但并不作为对上述总体方案的额外限定，仅仅是进一步的增补或优选，在没有技术或逻辑矛盾的前提下，各可选方式可单独针对上述总体方案进行组合，还可以是多个可选方式之间进行组合。
[0013] 可选地，利用小波分解提取心电图中QRS波群中的R波，在[tR-100ms，tR+100ms]范围内，tR时间之前第一个极小值数据点为Q波，tR时间之后第一个极小值数据点为S波。
[0014] 可选地，所述预定阈值为平滑心电图中最大值数据点电位幅值的0.4～0.6倍。
[0015] 可选地，所述Dmin的取值范围为0.13～0.2，所述Dmax的取值范围为0.5～0.9。
[0016] 一种心电图信号特征提取装置，包括：
[0017] Q波、R波和S波提取模块，用于提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
[0018] P波提取模块，用于在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
[0019] 若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
[0020] 若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；
[0021] T波提取模块，用于在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
[0022] 若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
[0023] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
[0024] 标记模块，用于依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
[0026] 提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
[0027] 在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
[0028] 若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
[0029] 若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；
[0030] 在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
[0031] 若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
[0032] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
[0033] 依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
[0034] 一种心电图信号特征提取系统，包括终端以及服务器，所述服务器包括存储器和处理器，所述存储器内存储有计算机程序，所述服务器从终端获取心电图；所述处理器执行所述计算机程序时，实现如下步骤：
[0035] 提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
[0036] 在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
[0037] 若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
[0038] 若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；
[0039] 在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
[0040] 若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
[0041] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
[0042] 依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
[0043] 一种计算机可读存储介质，所述计算机可读存储介质内存储有计算机程序，所述计算机程序被计算机处理器执行时实现如下步骤：
[0044] 提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
[0045] 在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
[0046] 若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
[0047] 若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；
[0048] 在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
[0049] 若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
[0050] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
[0051] 依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
[0054] 图2为一个实施例中的心电图；
[0055] 图3为一个实施例中心电图及其基线的示意图；
[0056] 图4为一个实施例中去除基线漂移后的心电图；
[0057] 图5为一个实施例中心电图小波变换的示意图；
[0058] 图6为一个实施例中提取出的R波的示意图；
[0059] 图7为一个实施例中提取出的Q波的示意图；
[0060] 图8为一个实施例中提取出的S波的示意图；
[0061] 图9为一个实施例中削平QRS波群后的心电图；
[0062] 图10为一个实施例中削平QRS波群且平滑后的平滑心电图；
[0063] 图11为一个实施例中提取出的P波的示意图；
[0064] 图12为一个实施例中提取出的T波的示意图；
[0065] 图13为一个实施例中R波、T波、P波的位置关系示意图；
[0066] 图14为一个实施例中计算机设备的示意图。
具体实施方式
[0067] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0068] 如图1所示，一种心电图信号特征提取方法，包括如下步骤：
[0069] 步骤100，提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图。
[0070] 心电图优选采用十二导联心电图，也可以采用其他导联方式得到的心电图。在心电图中截取Ⅱ导联中任意一段时间范围内的图像，该时间范围内至少包括三个心动周期，针对该段时间范围内的图像进行Q波、R波、S波的提取。
[0071] 心电图上的横轴为时间，单位为ms，纵轴为电位，单位为mV。
[0072] 如果针对图片形式的心电图进行信息提取，需将图片上的心电图图像转化为心电图数据信息。心电图图像为二值化图像，即白色背景上绘制黑色心电图图形，将心电图图像上的各黑色像素的位置信息转化为包含时间和电位的心电图数据信息，并根据心电图数据信息重新绘制得到如图2所示的心电图，针对图2所示的心电图进行后续特征波提取。
[0074] 去除基线漂移可采用现有技术，例如可采用多项式拟合的方法去除基线漂移，拟合方程为p(x)＝p1x6+p2x5+p3x4+p4x3+p5x2+p6x+p7，利用最小二乘法曲线拟合原理，得到最佳的拟合系数p1～p7，该拟合曲线即为心电图信号的基线，基线如图3中的平滑曲线所示。
[0075] 将原始心电图信息减去相应时间点处的基线电位，即去除基线漂移，去除基线漂移后的心电图如图4所示。
[0076] Q波、R波和S波的提取可采用现有技术，在一个实施例中，利用小波分解提取心电图中QRS波群中的R波，具体包括：
[0077] 步骤110，如图5所示，对心电图进行N层分解，每层小波的高频系数为d1,d2,d3,…dN；
[0078] 步骤120，遍历d4小波信号，找出d4小波信号的最大值dmax，将信号值大于0.15×dmax的信号组成信号序列；
[0079] 步骤130，i，j分别为信号序列中两个信号点对应的的时间指针，如果|i-j|<100ms，则i，j对应的两个信号点属于同一QRS波群；
[0080] 步骤140，提取同一QRS波群内极大值信号点对应的时间t1，在心电图中，在(t1-t，t1+t)时间段内寻找极大值数据点，如图6所示，该极大值数据点即为R波。极大值数据点的搜索过程为动态过程，t值以5～10ms为步长逐渐增加，直至在(t1-t，t1+t)时间段内寻找到一个极大值数据点。
[0081] 成年人QRS波群时间通常为40ms～100ms，婴幼儿QRS波群时间通常为40ms～80ms，以R波的时间点tR为起始点，分别向前向后追溯一段时间，即可得到Q波和S波，向前搜索得到的第一个负向极值为Q波，向后搜索得到的第一个负向极值为S波。
[0082] 如图7所示，在[tR-100ms，tR+100ms]范围内，tR时间之前第一个极小值数据点为Q波，如图8所示，tR时间之后第一个极小值数据点为S波。
[0083] 在心电图中，P波和T波为电位幅值较小的正向波，容易受到高频噪声影响，为了便于识别P波和T波，首先削平QRS波群，削平QRS波群的含义为，将Q波数据点和S波数据点直接连线，削平QRS波群后的心电图如图9所示，对削平QRS波群后的心电图做平滑处理(平滑处理可以采用现有技术)后得到的平滑心电图如图10所示。
[0084] 步骤200，在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找
[0085] [tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
[0086] 若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
[0087] 若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点。P波位置如图11所示。
[0088] 第一极大值数据点同时满足两个条件，其一，同一心动周期内的极大值，其二，大于预定阈值。
[0089] 在一个实施例中，所述预定阈值为平滑心电图中最大值数据点电位幅值的0.5倍。
[0090] 步骤300，在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
[0091] 若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期(识别R波后通过计算得到R波间期)；
[0092] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波。T波位置如图12所示；
[0093] 若[Dmin×RR,Dmax×RR]时间范围内不存在第二极大值数据点，则没有P波出现。
[0094] 在一个实施例中，所述Dmin的取值为0.17，所述Dmax的取值为0.7。
[0095] T波为S波之后的第一个正向波，RTmin＝Dmin×RR，RTmax＝Dmax×RR，如图13所示，RTmin为第一个R波波峰与T波的最小距离，RTmax为第一个R波波峰与T波的最大距离。
[0096] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则除了识别出T波，还同时识别出U波或其他噪声，选择离前一个R波波峰最近的极大值为T波。
[0097] 提取得到Q波、R波、S波、P波和T波，同时得到特征波的电位信息和特征波对应的时间点。
[0098] 步骤400，依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
[0099] 在一个实施例中N取值为15。
[0100] 例如，心电图检查中，P波频率>100次/分，判断为窦性心动过速。若某一心电图中，计算得到的P波频率在100在性心动次/分内，则判定为临界状态，需要筛选出来供医生判断。
[0101] 又如，P波频率<60次/分，判定为窦性心动过缓。若某一心电图中，计算得到的P波频率在60±15％次/分内，则判定为临界状态，需要筛选出来供医生判断。
[0102] 又如，PP间期差距>120ms，PR间期在120～200ms之间，判定为窦房结内游走性心律不齐。PP间距差距120±20距ms，PR间期在120±15％ms内为临界状态。
[0103] 又如，P波高尖，II导联振幅≥0.25mV，判定为右心房异常。若某一心电图中，P波电位在0.25±15％mV内为临界状态。
[0104] 应该理解的是，虽然图1的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
[0105] 在一个实施例中，提供了一种心电图信号特征提取装置，包括：
[0106] Q波、R波和S波提取模块，用于提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
[0107] P波提取模块，用于在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
[0108] 若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
[0109] 若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；
[0110] T波提取模块，用于在平滑心电图中，针对同一心动周期，在
[0111] [Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
[0112] 若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
[0113] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
[0114] 标记模块，用于依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
[0115] 关于各模块中功能限定可参见上文中对于心电图信号特征提取方法的限定，在此不再赘述。上述心电图信号特征提取装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
[0116] 本实施例提供的心电图信号特征提取装置可配置在远端，通过与之相连的远程终端获取心电图信号，还可以是本实施例装置本身就配置在终端(例如用户使用的计算机或医用检测设备)，直接通过心电图采集装置获取心电图信号。
[0117] 在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图14所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现心电图信号特征提取方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
[0118] 本领域技术人员可以理解，图14中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0119] 在一个实施例中，提供了一种计算机设备，包括：存储器和处理器，所述存储器内存储有计算机程序，所述处理器执行所述计算机程序时，实现如下步骤：
[0120] 提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
[0121] 在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
[0122] 若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
[0123] 若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；
[0124] 在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
[0125] 若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
[0126] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
[0127] 依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
[0128] 在一个实施例中，提供了一种心电图信号特征提取系统，包括终端以及服务器，所述服务器包括存储器和处理器，所述存储器内存储有计算机程序，所述服务器从终端获取心电图；所述处理器执行所述计算机程序时，实现如下步骤：
[0129] 提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
[0130] 在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
[0131] 若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
[0132] 若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；
[0133] 在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
[0134] 若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
[0135] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
[0136] 依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
[0137] 在一个实施例中，提供了一种计算机可读存储介质，所述计算机可读存储介质内存储有计算机程序，所述计算机程序被计算机处理器执行时实现如下步骤：
[0138] 提取心电图中QRS波群的Q波、R波和S波，在心电图中削平QRS波群，并对削平QRS波群后的心电图做平滑处理，得到平滑心电图；
[0139] 在平滑心电图中，针对同一心动周期，筛选出电位幅值大于预定阈值的若干第一极大值数据点，若[tR-200ms,tR]时间范围外存在所述第一极大值数据点，则舍弃这些第一极大值数据点，查找[tR-200ms,tR]时间范围内是否存在所述第一极大值数据点，并进行如下判断：
[0140] 若[tR-200ms,tR]时间范围内存在至少一个第一极大值数据点，则距离Q波最近的第一极大值数据点为P波；
[0141] 若[tR-200ms,tR]时间范围内不存在所述第一极大值数据点，则取[tR-200ms,tR]时间范围内的极大值数据点为P波，tR为R波的时间点；
[0142] 在平滑心电图中，针对同一心动周期，在[Dmin×RR,Dmax×RR]时间范围内筛选第二极大值数据点，并进行如下判断：
[0143] 若[Dmin×RR,Dmax×RR]时间范围内仅存在一个第二极大值数据点，则该第二极大值数据点为T波，Dmin、Dmax为系数，RR为R波间期；
[0144] 若[Dmin×RR,Dmax×RR]时间范围内存在至少两个第二极大值数据点，则离时间在前的R波波峰最近的第二极大值数据点为T波；
[0145] 依据Q波、R波、S波、P波和T波计算生理参数，若某一生理参数位于某一心脏状态判断阈值的±N％波动范围内，则将该生理参数标记为临界状态，N的取值范围为10～15。
[0146] 本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种基于磁信号的Cr-Ni-Fe合金蠕变损伤早期及失效临界状态的判断与预警方法 | 2020-05-08 | 697 |
| 机械发电/电力发电系统 | 2020-05-11 | 758 |
| 一种球蠕灰梯度石墨铸铁玻璃模具材料的制备方法 | 2020-05-11 | 220 |
| 一种梯级水电站优化调度系统、方法及安全电网系统 | 2020-05-11 | 971 |
| 一种高值耗材柜的自动控制方法 | 2020-05-12 | 867 |
| 一种考虑校准箱气体质量变化的文氏管流量系数标定方法 | 2020-05-12 | 912 |
| 判断起辉状态的方法和装置、气相沉积方法和设备 | 2020-05-13 | 423 |
| 桥梁桥面的行车稳定性分析方法及装置 | 2020-05-08 | 45 |
| 滤光片切换方法和滤光片切换装置 | 2020-05-12 | 328 |
| 一种基于视觉伺服的磁流体液滴运动控制系统 | 2020-05-12 | 507 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
临界状态热门专利
-
6 电池单元
QQ群二维码
意见反馈
意见反馈