一种心室除极关键点检测方法

摘要

本发明属于计算机技术与心电图技术的交叉领域，公开一种心室除极关键点检测方法。基于阈值法和局部距离变换法思想，对预处理过的ECG信号，在已知QRS波群主峰位置的条件下，对对应心室除极的多种形态QRS波群准确、快速地定位其起点、Q点、R点、S点和J点，并且计算简单、易于实现，为多种ECG处理问题提供一种新的解决途径，比如QRS波群分类、心电轴计算等。相对于小波变换方法，本方法速度更快、适用范围更广、更适用于商业使用；相对于传统的斜率比值法，本方法的抗干扰性更强、适用范围更广。

说明书

技术领域

本发明属于计算机技术与心电图技术的交叉领域，是一种基于阈值法和局部距离变换法的心室除极关键点检测方法。

背景技术

心电图(Electrocardiogram，ECG)记录了心脏在每个心动周期中心脏活动伴随着的电信号的变化，记录了心脏兴奋的发生、传播及恢复的过程，其形态变化反映了人体心脏功能，是判断心脏状态的参考依据，准确识别心电信号的变化特征显得至关重要。

心电信号的幅度范围大约在10μV-5mV，频率范围大约在0.05-100Hz。典型的心电信号波形如图2所示：心电信号由许多特征部分组成，分别为P波、Q波、R波、S波、T波，其中Q波、R波、S波构成QRS波群。QRS波群是心室除极在ECG上的表现。QRS波群起点对应于Q波峰之前的第一个斜率突变点，QRS波群的结束点对应于S波峰之后的第一个斜率突变点，这个点也叫做J点。QRS波群的起点、Q波峰值点、R波峰值点、S波峰值点与QRS波群结束点共同组成心室除极关键点。

对于ECG自动诊断系统来说，QRS波群起点、Q点、R点、S点、J点的准确定位极其重要。一些自动诊断算法以QRS波群起点的水平线作为基线，是计算其他各个波幅值的基础。QRS波群形态是诊断多种疾病的重要依据。其中以QRS波群的时限、类型(主要包括qRs，qR，Rs，QS，qS，rSR’，等等)、幅值、起点、J点等参数最为关键。

目前有公开的寻R算法有很多，但是在实际使用中并不理想，因其没有考虑QRS波群的多种形态，例如寻R算法对QS型、rSR’型无力；对病态的QRS波群检测的准确率低下。在寻QRS波群起点和结束点的算法上，目前还没有公认成熟的算法。目前亟需一种实用的准确的寻找心室除极关键点的算法。

发明内容

本发明为了克服现有算法存在的不足，提供一种适用于各种心室除极形态的通用算法，用以寻找QRS波群起点、Q点、R点、S点、QRS波群终点，其中QRS波群终点又称为J点。对采集到的人体心电数字信号标记为V(t)，其中V是电压，单位mV，t是时间，单位s。心电信号为离散值，每连续两个采样点之间的时间间隔是相等的，与采样率相关。该类问题首先需要进行信号质量评估，在评估结果的基础上选择去基线漂移算法及滤波算法，以得到相对理想的心电信号y(x)。后续操作均基于处理后的心电信号。

心室除极关键点检测步骤如下：

Step1，对经过去基线漂移及滤波处理过的心电信号y(x)，使用寻QRS波群主峰算法，寻得QRS波群主峰MP位置；该位置可能是R波主峰也可能是Q、S波主峰。在MP点左右0.07s内扫描所有候选波，候选波包括正向波和负向波；所述候选波定义为时限大于0.04s、高度大于0.15mV的波；由于干扰不可避免的存在和滤波产生的微弱形变，该标准不同于明尼苏达编码。

Step2，在候选波中，正向波峰值点描述为CR，负向波峰值点描述为CS；根据CR和CS分为以下四种情况：

1)当CR数量为0且CS数量等于1时，确定此情况下没有R波，则CS点为Q点与S点的重合点，此时QRS波群是QS型，CS点为CQ点；

2)当CR数量为1时，在CR点左0.06s内最低点，接近R者优先，作为CQ点；

3)当CR数量为2时，在第一个CR左0.06s内最低点，接近R者优先，作为CQ点；

4)当不属于1)、2)和3)的其他情况时，则QRS波群寻找失败，程序终止；

Step3，若不属于情况4)，根据局部距离变换法，过CQ(x_CQ,y_CQ)点做一条与过CQ点的垂线夹角为10°的直线LQ，此时直线LQ的表达式通过点斜式获取；在CQ点左侧，如果心电图散点中存在一点P(x_P,y_P)，使得y(x_p)>y_LQ(x_p)且y(x_p-1)<y_LQ(x_p-1)，则计算心电图散点在x<x_CQ且x>x_P时，与直线LQ的距离，且距离最大的心电图散点是QRS波群起点，CQ点为Q点；如果不存在上述点P，则QRS波群无Q波，且CQ点为QRS波群起点，无Q点；

至此QRS起点、Q点都可以确定；以过QRS波群起点的水平线作为基线，根据CR点是否高于基线判断该点是否是R点；R点为高于基线的CR点；

Step4，如果存在R点，连接R(x_R,y_R)点与之后0.2s的点M(x_M,y_M)，用两点式得到该直线RM的表达式，如果有两个R点选择右面的R点进行连接，计算心电图散点在x>x_R且x<x_M时，与直线RM的距离，距离最大值的点命名为CS1(x_CS1,y_CS1)；若不存在R点，则QRS波群为QS型，把CS点重命名为CS1点；

过CS1点做直线LJ与过CS1点的垂线夹角15°，如果存在一点P1(x_p1,y_p1),x_p1>x_CS1，使得y(x_p1)>y_LJ(x_p1)且y(x_p1+1)<y_LJ(x_p1+1)，则CS1点是S点；进一步计算心电图散点在x>x_CS1且x<x_P1时，与直线LJ的距离，距离最大值的ECG散点是J点；若不存在上述点P1，则QRS波群R波右面没有S波，且CS1点是J点；

至此，QRS波群起点、Q点、R点、S点、J点查找完毕。

本发明的有益效果为：基于阈值法和局部距离变换法思想，对预处理过的ECG信号，在已知QRS波群主峰位置的条件下，对对应心室除极的多种形态QRS波群准确、快速地定位其起点、Q点、R点、S点和J点，并且计算简单、易于实现，为多种ECG处理问题提供一种新的解决途径，比如QRS波群分类、心电轴计算、疾病诊断等。相对于小波变换方法，本方法速度更快、适用范围更广、更适用于商业使用；相对于传统的斜率比值法，本方法的抗干扰性更强、适用范围更广。

附图说明

图1为整体程序流程图。

图2为多种心电图QRS波群形态。其中子图(a)、(b)、(c)、(d)、(e)、(f)分别是qRs、Rs、R、qR、QS、rSR’型QRS波群示意图。

图3为候选区域和候选波示意图。其中子图(a)、(b)、(c)分别是qRs、rSR’、QS型时的候选区域和候选波示意图。

图4为QRS波群起点寻找示意图情况1。其中子图(a)是存在Q波时的示意图，子图(b)是关键部分的放大图。

图5为QRS波群起点寻找示意图情况2，是不存在Q波时的示意图。

图6为S点寻找示意图。其中子图(a)、(b)、(c)是在三种典型QRS形态时的示意图。

图7为J点寻找示意图。其中子图(a)、(b)、(c)、(d)是在四种典型QRS形态时的示意图。

具体实施方式

心电信号反映了心肌细胞有规律的除极复极过程，采用心电检测仪采集出人体心电信号V(t)，首先经过去工频干扰，然后通过信号质量评估分析噪声等级并针对不同等级噪声采用不同方式的滤波。此过程去除了基线漂移、工频干扰、肌电干扰等检测时参杂的干扰信号，并获得较为纯净的心电信号，以保证信号质量的可靠性。

QRS波群可以有多种形态，如附图2，不同形态之间的差异度很大，心室除极的关键点中的Q点、R点、S点很多时候不会规律出现，有时R点、S点会出现两次。但是在变之中还存在着不变，即QRS波群起点、J点必然存在；Q、R、S点必然不会同时都不存在；若Q波存在，必定在R波左面一定区域；若S波存在，一定在R波右面一定区域。利用这些不变性质可以计算所有关键点。

Step1，把ECG信号放在二维坐标系中，电压做Y轴(单位：毫伏mV)，时间做X轴(单位：秒s)。首先对处理过的心电信号y(x)，使用现存的寻QRS波群主峰算法寻得QRS波群主峰MP位置。MP点必然在QRS波群之中。在MP点左右0.07s之内寻找候选波，0.14s可以在很大程度上覆盖了QRS波群R点出现的区域。根据明尼苏达编码，结合滤波对波形的损失，规定时限大于0.04s且高度大于0.15mV的波属于候选波。

Step2，候选波分正向峰和负向峰，峰值点分别为CR和CS。接下来根据正向峰的数量做出如下判断：

当CR数量为0且CS数量为1时，可判断当前QRS波群QS型，且负向波CS为CQ点和S点(特殊情况，Q和S是同一点)，无R波。

当CR数量为1时，在CR点左0.06s内的最小值点定义为CQ点。

当CR数量为2时，在第一个CR点左0.06s内的最小值点定义为CQ点。

否则QRS波群检测有误，程序终止。

Step3，检测到CQ(x_CQ，y_CQ)点之后，通过CQ点做一条与过CQ点的垂线VQ夹角为10°的直线LQ，夹角在第二象限。

在已知斜率和过某一点的情况下可得到直线LQ的公式：

y_LQ＝k₁(x_LQ-x_CQ)+y_CQ

其中已知(x_CQ，y_CQ)，k₁＝tan(10°)。

如果存在一点P(x_P，y_P)，x_P＜x_CQ，使得y(x_P)＞y_LQ(x_P)且y(x_P-1)＜y_LQ(x_P-1)则：

直线LQ与ECG在CQ点左存在交点，且交点为P。通过点到直线的距离公式可得到任意点到直线LQ的表达式，其中x是自变量：

带入所有ECG上的点(x，y)，x＜x_CQ且x＞x_P，则QRS波群起点为使得D₁(x)达到最大值的点。

如果不存在上述点P，则QRS波群无Q波，且CQ点为QRS波群起点。

以过QRS波群起点的水平线作为基线，可根据候选CR点是否高于基线判断该点是否是R点。目前为止，QRS起点、Q点、R点都可以确定，接下来寻找S点和J点。

Step4，对于寻找S点，采用局部距离变换法。具体过程与寻QRS波群起点相似：如果存在R，连接R点(x_R，y_R)(如果存在两个R则选择右面R)与之后0.2s的点M(x_M，y_M)，该直线的表达式：

其中自变量是x_LS，因变量是y_LS。

公式可以简写成：

y_LS＝k₂x_LS+b

并计算这0.2s内距离这条直线最远的点CS1(x_CS1，y_Cs1)。

如果不存在R，此时QRS波群是QS型或qS型，把CS点命名为CS1。

已知，过QRS波群的S点与过S点的垂线的夹角在没有外界干扰的情况下会小于15度。因此过CS1点作该直线LJ时，如果存在S波，该直线会与ECG相交于一点P1(x_P1，y_P1)，交点之前y_LJ(x_p1)＜y(x_p1)，之后y_LJ(x_p1+1)＞y(x_P1+1)；如果无S波时，该直线与ECG散点在ST段无交点P1。综上，过CS1点做直线LJ与过CS1点的垂线夹角15°，该直线表达式用点斜式：

y_LJ＝k₃(x_LJ-x_CS1)+y_CS1

其中已知(x_CS1，y_CS1)，k₃＝tan(15°)。

设动点(x，y)，x＞x_CS1，如果存在P1点，则CS1点是S(x_S，y_S)点；否则，该QRS波群R右面没有S波，且CS1点是J点。

若存在S(x_S，y_S)点，再次利用上述的直线LJ，利用局部距离变换法计算在S点与P1点之间距离直线LJ最远的点。过程：

计算ECG上的点(x，y)，在x＞x_S且x＜x_P1时，使得D₂(x)达到最大值的点即为J点。

综上，QRS波群起点、Q点、R点、S点、J点全部寻找完毕。