[发明专利]一种心电信号T波及其终点自动检测方法

权利要求书

1.一种心电信号T波及其终点自动检测方法，其特征在于，该方法能获得待检测心信号序列的T峰和T波终点位置点集合；包括以下步骤：

步骤1、对待检测的心电信号序列进行预处理，获得去除噪声的预处理后的心电信号序列；

步骤2、从预处理后的心电信号序列中获取多个R峰位置点，以获得R峰位置点集合；

步骤3、根据步骤2获得的R峰位置点集合中每一个R峰位置点，利用极值原理获得每一个T峰位置点，以获得T峰位置点集合；

所述步骤3包括：

步骤3.1、在R峰位置点集合中计算相邻两个R峰位置点之间的水平距离值，获得RR间期；

步骤3.2、在R峰位置点集合中，以每一个R峰位置点为起点，根据T波位置特点和时间长度，向后进行加窗处理，获得预选T波片段序列；

步骤3.3、计算R峰至S峰的黎曼和S_RS、预选T波片段序列的黎曼和S_T以及R峰至S峰心电片段序列和预选T波片段序列的总黎曼和S；比较S与S_RS的大小，若S＞S_RS，则T波为正向波；若S＜S_RS，则T波为倒置波；

步骤3.4、在R峰位置点集合中，以每一个R峰位置点为起点，向后对预处理后的心电信号序列进行加窗处理，获得多个T峰搜索区间；若T波为正向波，则在每个搜索区间中利用极大值方法检测极大幅值对应的位置点，获得多个T峰位置点；若T波为倒置波，则在每个搜索区间中利用极小值方法检测极小幅值对应的位置点，获得多个T峰位置点；最终获得T峰位置点集合；

步骤4、根据步骤3获得的T峰位置点集合中每一个T峰位置点，通过拐点搜索法获得每一个T波终点位置点，以获得T波终点位置点集合。

2.如权利要求1所述的心电信号T波及其终点自动检测方法，其特征在于，所述步骤2包括：

步骤2.1、采用小波变换对预处理后的心电信号序列进行分解得到多个分解层；根据QRS波群所在的频带范围选择对应分解层的小波系数，并据此对该分解层信号进行重构，获得重构心电信号序列；

步骤2.2、选取重构心电信号序列的边缘波峰点，通过曲线拟合的方法对其进行拟合，获得拟合心电信号序列；在拟合心电信号序列中，用模极大值法计算拟合心电信号序列中所有局部极值点，获得局部极值点集合；

步骤2.3、在局部极值点集合中，以局部极值点为中心，对重构心电信号序列进行加窗处理，获得心电信号片段序列；

步骤2.4、在每一个心电信号片段序列中选取最大幅值对应的位置点，获得候选R峰位置点集合；

步骤2.5、根据候选R峰位置点集合，获得R峰位置点集合。

3.如权利要求2所述的心电信号T波及其终点自动检测方法，其特征在于，所述步骤2.5包括：

步骤2.5.1、在候选R峰位置点集合中计算相邻两个R峰位置点之间的水平距离值，获得多个RR间期，每个RR间期对应两个R峰位置点；

步骤2.5.2、在RR间期序列中找到大于第一阈值THR_rr1的RR间期，并在心电信号序列中找到其所对应的一段漏检心电信号子序列；重复执行本步骤，直至找到所有漏检心电信号子序列；其中，第一阈值THR_rr1＝1.0*fs，fs为心电信号采样率，单位为Hz；

步骤2.5.3、在每段漏检心电信号子序列中，将最大幅值对应的位置点作为新的R峰位置点，并将其添加至候选R峰位置点集合中，获得插入漏检点后新的候选R峰位置点序列；

步骤2.5.4、在新的候选R峰位置点集合中，计算相邻两个R峰位置点之间的水平距离值，获得多个新的RR间期，每个新的RR间期对应两个R峰位置点；

步骤2.5.5、在新的RR间期序列中找到小于第二阈值THR_rr2的新的RR间期，并在心电信号序列上找到其所对应的两个R峰位置点，获得多个初选误检R峰位置点；重复执行本步骤，直至找到所有初选误检R峰位置点；其中，第二阈值THR_rr2＝0.6*fs，fs为心电信号采样率，单位为Hz；

步骤2.5.6、比较小于第二阈值的新的RR间期对应的两个R峰位置点的幅值，删去幅值较小的点对应的误检R峰位置点，获得剔除误检点后新的R峰位置点集合，从而得到最终确定的R峰位置点集合。