一种心电图低电压T波终点定位方法与流程

一种心电图低电压t波终点定位方法
技术领域
1.本发明涉及一种心电图低电压t波终点定位方法，具体的说，为一种基于波形重构和叠加波聚类技术的心电图低电压t波终点定位方法，属于心电图t波定位技术领域。


背景技术：

2.临床上很多心律失常呈阵发性、持续时间短等特点，就诊时心律失常可能已经终止，常规心电图检查往往难以捕获。而小型化、便携式的单导联心电设备可实现连续、长时程的心电监护，逐渐成为一种新的心电监护解决方案。随着近年来人工智能技术的快速发展，利用ai设计单导联心电算法帮助医生进行辅助诊断已经成为了必然趋势。心电图波形一个完整心搏周期主要是p波、qrs波群和t波组成，每种波都反映了心脏心房和心室中兴奋波的动态传导，以及异位搏动情况的丰富信息，准确检测各波形位置是心率计算和异常诊断的关键。
3.目前各领域专家学者通过不同的方法进行了心搏各波形位置检测，其中对t波检测的方法集中在切线法、斜率阈值法、小波变换法、数学形态法等方面。t波常见的形态有t波倒置、t波双向和t波低平等，由于t波是心电信号中的低幅低频波，易受基线漂移、肌电干扰等噪声影响，其终点检测的方法研究始终是心电图识别技术中的难点问题，特别是对变化比较敏感的低平t波（如图1），其斜率更不明显，容易淹没在噪声中，这就使得众多t波终点定位算法失去了对低平t波终点检测的准确性，进而直接影响qt间期、尼亚加拉状t波等的测量，导致对常见类型心电图心搏的识别造成困难。


技术实现要素：

4.针对上述现有技术的不足，本发明提供了一种心电图低电压t波终点定位方法，包括以下步骤：步骤一：单导联心电信号获取及预处理，根据预设的采样率，采用干电极获取t秒固定长度的单导联ecg样本，并对每个样本进行滤波去肌电、工频干扰和去基线漂移处理；步骤二：心电样本数据4个等级干扰的识别；步骤三：样本心搏叠加波构建，样本心搏叠加波构建，为后续步骤非线性变换重构的低平t波聚类分析初始模板建立基础，针对预处理后的心电样本，基于正则化最小二乘回归学习方法定位r波位置，以检测到的r波位置为中心，向前向后截取一定数量采样点的固定窗口对样本进行分割，共分割得到la个心搏，根据上一步骤的干扰识别，对只有1级和2级干扰级别的l个心搏进行聚类分析，共分为四个阶段，以获取心搏数量最多的两类心搏构建叠加波；步骤四：筛选心搏低平t波样本；并对其所有心搏划分t波终点搜索区域；步骤五：心搏t波终点搜索区域单独进行低通滤波消除干扰；针对低平t波样本叠加波及其所有单个心搏的t波终点搜索区域单独进行低通滤波，再次消除肌电干扰和工频干扰；
步骤六：低平t波样本所有心搏及叠加波t波终点搜索区域的重构；步骤七：重构后的心搏t波终点搜索区域聚类分析；步骤八：心搏t波终点搜索区域波形特征提取；样本所有重构后的心搏t波终点搜索区域进行特征提取，描述每个心搏搜索区域数据有意义的属性，该特征集包括了形态特征群、统计特征群、信息特征群;步骤九：模型训练；根据上述步骤八提取的心搏t波终点搜索区域的形态特征群、统计特征群、信息特征群，结合对应的t波终点到t波顶点的距离标签，将训练集中叠加波相似的心搏、不相似的心搏，分别输入到随机森林模型中进行训练，得到两种t波终点定位模型；步骤十：模型评估。
5.作为优选：所述步骤二中的4个等级干扰的识别方法如下：首先将预处理后的心电样本数据按照n秒钟长度切片，再进行间隔m个点的下采样，获得心电数据x，并对其进行标注得到y；其次，对心电数据x进行傅里叶变换f，并将其实部标准化norm之后和x相加，得到输出xf；再次，训练卷积神经网络，设计一个包含四层卷积层的一维卷积神经网络来提取特征，同时在每第二个卷积层之后加一个maxpooling层用来筛选特征，并规定每一层合适的通道数c，卷积核长度k和卷积步长strides，随后接一个averagepool层和两层全连接层来解码提取到的特征，全连接层的通道数分别为c和4，第二个全连接层的输出为xfc，xfc通过softmax层，得到4个类别概率的一维向量xout，设定学习率lr、损失函数loss，用adam优化器进行训练，直到模型收敛；最后，导入测试数据进行4个等级干扰的识别，对比最终输出的4个概率xp，取其中最高概率对应的等级作为n秒钟的干扰等级，公式如下：。
6.作为优选：所述步骤三中的四个阶段分别为：第一阶段，在对样本l个心搏进行聚类前，首先需要产生模板心搏，将目标心搏与模板心搏进行相似度距离计算，相似度距离满足一定阈值条件的目标心搏匹配归入相应模板类别，选取相似度距离阈值th，从样本的l个心搏中随机抽取ls个心搏，并选择第一个心搏作为第一个模板心搏类型，计算第二个心搏到的相似度距离，如果满足
＞th，则将匹配归入该类别，否则以为新建模板心搏，继续计算下一个心搏与已有的若干个模板心搏间相似度距离(j=0,1,2,
…
)，若满足，且＞th，则将该心搏匹配归于模板类别， 否则以为新建模板心搏类别，重复此步骤，直至完成随机ls个心搏的聚类，最后对随机的ls个心搏模板分类结果进行统计排序，选取心搏匹配个数最多的模板类别，将ls个心搏中属于的所有心搏数据做平均，得到一个平均模板心搏；第二阶段，对样本l个心搏第一次进行聚类，获取心搏匹配个数最多的一类，将上一阶段得到的平均模板心搏作为第一个心搏模板类型，重复第一阶段的心搏匹配模板的步骤，直到所有l个心搏匹配完毕，最后对全部l个模板心搏分类结果进行统计排序，得到心搏匹配个数最多的模板类别，获取其匹配的心搏集合及个数；第三阶段，样本中全部l个心搏去掉上一阶段最多一类个心搏后，再对剩余个的心搏第二次进行聚类，对个心搏重复第一阶段步骤，得到随机抽取心搏匹配个数最多的模板类别，接着对个心搏重复第二阶段步骤，得到匹配个数最多的模板类别，获取其匹配心搏集合及其个数，将第二阶段和第三阶段获得的心搏集合、合并为一个心搏集合；第四阶段，针对每个心电样本构建独特的心搏叠加波，将上述阶段获取的样本心搏集合用来构建叠加波，所有用于构建的心搏以r波位置为基准点向前向后对齐，并对每个心搏采样点对应的振幅幅值进行排序，截取第15百分位至第85百分位对应的心搏幅值求和取平均，作为当前采样点对应的叠加波振幅幅值，以此来构建叠加波，公式如下：
7.其中m是当前心搏的第m个采样点，和分别是每个心搏当前采样点对应的振幅幅值排序之后处于在15%和85%的位置，是叠加波第m个采样点的幅值，是第m个采样点处于第i百分位的心搏幅值。
8.作为优选：所述步骤四的具体方法为：针对样本所有心搏及上一步骤构建的样本叠加波，进行t波顶点的检测，选取t波顶点幅值阈值，筛选叠加波t波顶点绝对值最大值小于的心电样本，定义为低平t波样本，共筛选出份训练集和份测试集的低平t波单导联ecg样本，其中包括了窦性节律样本、室上性心律失常样本、室性心律失常样本、房扑房颤样本、传导阻滞样本等，然后以样本各心搏t波顶点为基准点，设定对应心搏低平t波终点的搜索区域为：。
9.作为优选：所述步骤五的具体方法为： 无限脉冲响应滤波器iir对低平t波终点搜
索区域低通滤波，具体设计了巴特沃斯数字低通滤波器，通带和阻带的截止频率为和，通带允许最大衰减是，阻带允许最大衰减是，iir滤波器差分方程的一般表达式如下：
10.其中x(n)为输入序列，y(n)为输出序列，和为滤波器系数，m为滤波器阶数；最后，将搜索区域的低平t波序列按顺序滤波，将所得结果在时域翻转后反向通过滤波器，再将所得结果在时域翻转后输出，使得相移为0，保证滤波后的波形位置精确地保留在滤波前序列中的位置。
11.作为优选：所述步骤六中对低平t波样本所有心搏及叠加波t波终点搜索区域的重构方法为： 首先将上一步骤处理后的样本心搏及其叠加波的t波终点搜索区域数据进行归一化处理，使其映射到合理区间：
12.其中x、为搜索区域数据映射前和映射后的值，分别为搜索区域数据最大值和最小值，a、b为映射区间的上限和下限；其次，通过e的指数运算将归一化后的数据进行非线性变换，从而使搜索区域的低平t波数据得到增强，最后，为使增强后的t波终点搜索区域具有更好的噪声抑制性能,对样本心搏及其叠加波增强后的t波终点搜索区域数据进行sg平滑滤波的重构，原理是基于搜索区域特征的多项式拟合，应用最小二乘法确定加权系数进行移动窗口加权平均滤波，具体公式如下：
13.其中x是原始搜索序列，是重构后的搜索序列，是滑动窗口的第i个时间序列的相关系数，2m+1是滑动窗口的大小。
14.作为优选：所述步骤七中重构后的心搏t波终点搜索区域聚类分析方法为：以步骤六重构得到的叠加波t波终点搜索区域为初始模板，将重构后的样本l个心搏t波终点搜索区域，重复步骤三第一阶段的心搏匹配模板的步骤，直到样本所有个心搏t波终点搜索区域匹配完毕，获取其匹配的心搏t波终点搜索区域集合及其补集；通过对低平t波的份训练集和份测试集样本区分叠加波相似的心搏和不相似的心搏，在后续步骤分开进行t波终点的建模和预测。
15.作为优选：所述步骤八中心搏t波终点搜索区域波形特征提取，具体提取的特征群为：1）形态特征群一阶差分：心搏t波终点搜索区域数据相邻两点之差；二阶差分：心搏t波终点搜索区域数据一阶差分相邻两点之差；心搏t波顶点幅值；rr间期；心搏r波顶点到t波顶点的距离；心搏r波顶点到t波顶点的距离和rr间期的比值；2）统计特征群心搏t波终点搜索区域数据平均值；心搏t波终点搜索区域数据标准差；心搏t波终点搜索区域数据峭度；心搏t波终点搜索区域数据偏度；心搏t波终点搜索区域数据平方的平均值；心搏t波终点搜索区域数据最大值和最小值的差值；一阶差分数据绝对值最小值位置距离t波顶点位置的距离二阶差分数据绝对值最小值位置距离t波顶点位置的距离3）信息特征群心搏t波终点搜索区域数据模糊熵；心搏t波终点搜索区域数据功率谱熵；心搏t波终点搜索区域数据奇异谱熵；心搏t波终点搜索区域数据排列熵。
16.作为优选：所述步骤十中的模型评估方法为：将训练好的两种低平t波定位模型应用在测试集的样本心搏数据上，测试集样本中的心搏分别应用其所属叠加波相似和不相似两种类型模型进行预测，将最后的t波终点位置输出结果与测试集样本心搏的t波终点位置标签进行对比，误差间隔是用来衡量波形定位的误差区间，当算法检测结果和医生标注位置位于误差区间范围内时，检测结果被认可，使用误差间隔内的样本占比评估模型的性能，各项评估指标的表达式如下：k：模型输出的低平t波终点在标注位置的
±
4ms范围内的样本数；l：模型输出的低平t波终点在标注位置的
±
5ms~
±
8ms范围内的样本数；m：模型输出的低平t波终点在标注位置的
±
9ms~
±
16ms范围内的样本数；n：模型输出的低平t波终点在标注位置的
±
16ms外的样本数；误差范围
±
4ms以内的样本占比：
17.误差范围
±
8ms以内的样本占比：
18.误差范围大于
±
16ms的样本占比：
19.其中误差间隔在本发明低平t波终点定位算法中与采样点关系如下：误差间隔4ms=1个采样点；误差间隔8ms=2个采样点；误差间隔16ms=4个采样点。
20.本发明在于对低平t波非线性变换的放大重构，增强了低平t波终点的细节特征，根据构建的样本心搏叠加波，结合样本心搏和叠加波的重构操作，聚类分析成两类心搏，增加的两种状态引入了医学信息，解决了两类心搏数据对应的提取特征高度相似问题，最后提取形态特征群、统计特征群、信息特征群分别建模预测，在实际应用场景下实现了对心搏低平t波终点的准确定位。
附图说明
21.图1是典型的低平t波ecg心搏示意图；图2是本发明的总体算法流程图；图3是本发明步骤二中4个等级干扰识别的算法流程图；图4是本发明步骤五中对t波终点搜索区域单独低通滤波消除干扰前后的心搏；图5是本发明步骤六中重构后的心搏和叠加波的t波终点搜索区域；图6是本发明步骤七中典型的两种聚类结果类型（与叠加波相似类型或与叠加波不相似类型）。
具体实施方式
22.下面将结合附图对本发明作详细的介绍：如图2所示，一种心电图低电压t波终点定位方法，所述方法包括如下步骤：步骤一：单导联心电信号获取及预处理。根据250hz的采样率，采用干电极获取30秒固定长度的单导联ecg样本，并对每个样本进行滤波去噪，去肌电、工频干扰和去基线漂移处理。
23.步骤二：心电样本数据4个等级干扰的识别。本步骤根据本公司专利一种可穿戴设备上的ecg信号的干扰波识别方法实现，具体算法流程如图3，其原理是首先将预处理后的心电样本数据按照3秒钟长度切片，再进行间隔2个点的下采样，获得心电数据x，并对其进行标注得到y；其次，对心电数据x进行傅里叶变换f，并将其实部标准化norm之后和x相加，得到输出xf；再次，训练卷积神经网络，设计一个包含四层卷积层的一维卷积神经网络来提取特征，同时在每第二个卷积层之后加一个maxpooling层用来筛选特征，其中卷积核长度k均为5，前两个卷积的步长strides为2，后两个为1，前三个卷积的通道数c为64，最后一个为128，随后接一个averagepool层和两层全连接层来解码提取到的特征，全连接层的通道数
分别为64和4，第二个全连接层的输出为xfc，xfc通过softmax层，得到4个类别概率的一维向量xout，设定学习率0.0001、损失函数loss，用adam优化器进行训练，直到模型收敛；最后，导入测试数据进行4个等级干扰的识别，对比最终输出的4个概率xp，取其中最高概率对应的等级作为此n秒钟的干扰等级。具体公式如下：
24.步骤三：样本心搏叠加波构建，为后续步骤非线性变换重构的低平t波聚类分析初始模板建立基础。针对预处理后的心电样本，基于正则化最小二乘回归学习方法定位r波位置，以检测到的r波位置为中心，向前取75ms采样点、向后取125ms采样点的固定窗口对样本进行分割，共分割得到la个心搏，根据上一步骤的干扰识别，对只有1级和2级干扰级别的l个心搏进行聚类分析，共分为四个阶段，以获取心搏数量最多的两类心搏构建叠加波。
25.第一阶段，在对样本l个心搏进行聚类前，首先需要产生模板心搏，将目标心搏与模板心搏进行余弦相似度计算，余弦相似度满足一定阈值条件的目标心搏匹配归入相应模板类别。选取余弦相似度阈值0.8，从样本的l个心搏中随机抽取10个心搏，并选择第一个心搏作为第一个模板心搏类型，计算第二个心搏到的余弦相似度距离，如果满足＞0.8，则将匹配归入该类别，否则以为新建模板心搏，继续计算下一个心搏与已有的若干个模板心搏间相似度距离(j=0,1,2,
…
)，若满足，且＞0.8，则将该心搏匹配归于模板类别，否则以为新建模板心搏类别，重复此步骤，直至完成随机10个心搏的聚类。最后对随机的10个心搏模板分类结果进行统计排序，选取心搏匹配个数最多的模板类别，将10个心搏中属于的所有心搏数据做平均，得到一个平均模板心搏。余弦相似度距离公式如下：
26.其中k表示当前样本的第k个心搏，j表示第j个模板心搏，表示当前样本第个心搏采样点的幅值，表示第个模板心搏采样点的幅值；第二阶段，对样本l个心搏第一次进行聚类，获取心搏匹配个数最多的一类。将上一阶段得到的平均模板心搏作为第一个心搏模板类型，重复第一阶段的心搏匹配模板的步骤，直到所有l个心搏匹配完毕。最后对全部l个模板心搏分类结果进行统计排序，得到心搏匹配个数最多的模板类别，获取其匹配的心搏集合及个数；第三阶段，样本中全部l个心搏去掉上一阶段最多一类个心搏后，再对剩余个的心搏第二次进行聚类。对个心搏重复第一阶段步骤，得到随机抽取心搏匹配个数最多的模板类别，接着对个心搏重复第二阶段步骤，得到匹配个数最多的模板类别，获取其匹配心搏集合及其个数。将第二阶段和第三阶段获得的心搏集合、合并为一个心搏集合。
27.第四阶段，针对每个心电样本构建独特的心搏叠加波。将上述阶段获取的样本心搏集合用来构建叠加波，所有用于构建的心搏以r波位置为基准点向前向后对齐，并对每个心搏采样点对应的振幅幅值进行排序，截取第15百分位至第85百分位对应的心搏幅值求和取平均，作为当前采样点对应的叠加波振幅幅值，以此来构建叠加波。公式如下：
28.其中m是当前心搏的第m个采样点，和分别是每个心搏当前采样点对应的振幅幅值排序之后处于在15%和85%的位置，是叠加波第m个采样点的幅值，是第m个采样点处于第i百分位的心搏幅值。
29.步骤四：筛选心搏低平t波样本,并对其所有心搏划分t波终点搜索区域。针对样本所有心搏及上一步骤构建的样本叠加波，进行t波顶点的检测，选取t波顶点幅值阈值为0.1mv,筛选叠加波t波顶点绝对值最大值小于0.1mv的心电样本，定义为低平t波样本，共筛选出8115份训练集和2074份测试集的低平t波单导联ecg样本，训练集中窦性节律样本2117份，室上性心律失常样本1618份，室性心律失常样本1516份，房扑房颤样本1395份，传导阻滞样本1469份，测试集中窦性节律样本451份，室上性心律失常样本406份，室性心律失常样本417份，房扑房颤样本402份，传导阻滞样本398份。然后以样本各心搏t波顶点为基准点，设定对应心搏低平t波终点的搜索区域为：。
30.步骤五：心搏t波终点搜索区域单独进行低通滤波消除干扰。根据资料表明，心电信号的整体频率范围为0.05~100hz，其能量主要集中在 0.5~45hz，t波频率范围为0.7~10hz，考虑到步骤一是对包括qrs波在内的整个心电波形的预处理，去除了心电波形上的大部分高频噪声，但本发明的低平t波终点定位算法需要使用到信号的细节特征，所以需要针
对低平t波样本叠加波及其所有单个心搏的t波终点搜索区域单独进行低通滤波，再次消除肌电干扰和工频干扰。处理后的心搏如图4所示，本发明采取的方法是无限脉冲响应滤波器iir对低平t波终点搜索区域低通滤波，具体设计了巴特沃斯数字低通滤波器，通带和阻带的截止频率为12hz和17hz，通带允许最大衰减是3db，阻带允许最大衰减是20db，iir滤波器差分方程的一般表达式如下：其中x(n)为输入序列，y(n)为输出序列，和为滤波器系数，m为滤波器阶数。
31.最后，将搜索区域的低平t波序列按顺序滤波，将所得结果在时域翻转后反向通过滤波器，再将所得结果在时域翻转后输出，使得相移为0，保证滤波后的波形位置精确地保留在滤波前序列中的位置。
32.步骤六：低平t波样本所有心搏及叠加波t波终点搜索区域的重构，重构后的结果如图5所示。首先将上一步骤处理后的样本心搏及其叠加波的t波终点搜索区域数据进行归一化处理，使其映射到[1,5]区间：其中x、为搜索区域数据映射前和映射后的值，分别为搜索区域数据最大值和最小值，a、b为映射区间的上限5和下限1。
[0033]
其次，通过e的指数运算将归一化后的数据进行非线性变换，从而使搜索区域的低平t波数据得到增强。最后，为使增强后的t波终点搜索区域具有更好的噪声抑制性能,对样本心搏及其叠加波增强后的t波终点搜索区域数据进行sg平滑滤波的重构，原理是基于搜索区域特征的多项式拟合，应用最小二乘法确定加权系数进行移动窗口加权平均滤波，具体公式如下：
[0034]
其中x是原始搜索序列，是重构后的搜索序列，是滑动窗口的第i个时间序列的相关系数，2m+1是滑动窗口的大小。
[0035]
步骤七：重构后的心搏t波终点搜索区域聚类分析。考虑到很多的搜索区域数据对应的提取特征高度相似，不存在一一对应关系，故对样本所有重构后的心搏t波终点搜索区域进行聚类分析，分成和叠加波相似和不相似的两种心搏。以上一步骤重构得到的叠加波t波终点搜索区域为初始模板，将重构后的样本l个心搏t波终点搜索区域，重复步骤三第一阶段的心搏匹配模板的步骤，直到样本所有l个心搏t波终点搜索区域匹配完毕，获取其匹
配的心搏t波终点搜索区域集合及其补集。通过对低平t波的份训练集和份测试集样本区分叠加波相似的心搏和不相似的心搏，在后续步骤分开进行t波终点的建模和预测。典型的两种聚类结果如图6所示。
[0036]
步骤八：心搏t波终点搜索区域波形特征提取。对样本所有重构后的心搏t波终点搜索区域进行特征提取，描述每个心搏搜索区域数据有意义的属性，该特征集包括了形态特征群、统计特征群、信息特征群：1）形态特征群一阶差分：心搏t波终点搜索区域数据相邻两点之差；二阶差分：心搏t波终点搜索区域数据一阶差分相邻两点之差；心搏t波顶点幅值；rr间期；心搏r波顶点到t波顶点的距离；心搏r波顶点到t波顶点的距离和rr间期的比值；2）统计特征群心搏t波终点搜索区域数据平均值；心搏t波终点搜索区域数据标准差；心搏t波终点搜索区域数据峭度；心搏t波终点搜索区域数据偏度；心搏t波终点搜索区域数据平方的平均值；心搏t波终点搜索区域数据最大值和最小值的差值；一阶差分数据绝对值最小值位置距离t波顶点位置的距离二阶差分数据绝对值最小值位置距离t波顶点位置的距离3）信息特征群心搏t波终点搜索区域数据模糊熵；心搏t波终点搜索区域数据功率谱熵；心搏t波终点搜索区域数据奇异谱熵；心搏t波终点搜索区域数据排列熵；步骤九：模型训练。根据上述步骤八提取的心搏t波终点搜索区域的形态特征群、统计特征群、信息特征群，结合对应的t波终点到t波顶点的距离标签，将训练集中叠加波相似的心搏、不相似的心搏，分别输入到随机森林模型中进行训练，得到两种t波终点定位模型。
[0037]
步骤十：模型评估。将训练好的两种低平t波定位模型应用在2074份低平t波测试集的样本心搏数据上，测试集样本中的心搏分别应用其所属叠加波相似和不相似两种类型模型进行预测，将最后的t波终点位置输出结果与测试集样本心搏的t波终点位置标签进行对比。误差间隔是用来衡量波形定位的误差区间，当算法检测结果和医生标注位置位于误差区间范围内时，检测结果被认可，使用误差间隔内的样本占比评估模型的性能，各项评估指标的表达式如下：k：模型输出的低平t波终点在标注位置的
±
4ms范围内的样本数；l：模型输出的低平t波终点在标注位置的
±
5ms~
±
8ms范围内的样本数；
m：模型输出的低平t波终点在标注位置的
±
9ms~
±
16ms范围内的样本数；n：模型输出的低平t波终点在标注位置的
±
16ms外的样本数；误差范围
±
4ms以内的样本占比：
[0038]
误差范围
±
8ms以内的样本占比：
[0039]
误差范围大于
±
16ms的样本占比：
[0040]
其中误差间隔在本发明低平t波终点定位算法中与采样点关系如下：误差间隔4ms=1个采样点误差间隔8ms=2个采样点误差间隔16ms=4个采样点。
[0041]
不同诊断类型的低平t波心搏的测试结果如表1。
[0042]

技术特征：
1.一种心电图低电压t波终点定位方法，其特征在于：所述方法包括以下步骤：步骤一：单导联心电信号获取及预处理，根据预设的采样率，采用干电极获取t秒固定长度的单导联ecg样本，并对每个样本进行滤波去肌电、工频干扰和去基线漂移处理；步骤二：心电样本数据4个等级干扰的识别；步骤三：样本心搏叠加波构建，样本心搏叠加波构建，为后续步骤非线性变换重构的低平t波聚类分析初始模板建立基础，针对预处理后的心电样本，基于正则化最小二乘回归学习方法定位r波位置，以检测到的r波位置为中心，向前向后截取一定数量采样点的固定窗口对样本进行分割，共分割得到la个心搏，根据上一步骤的干扰识别，对只有1级和2级干扰级别的l个心搏进行聚类分析，共分为四个阶段，以获取心搏数量最多的两类心搏构建叠加波；步骤四：筛选心搏低平t波样本；并对其所有心搏划分t波终点搜索区域；步骤五：心搏t波终点搜索区域单独进行低通滤波消除干扰；针对低平t波样本叠加波及其所有单个心搏的t波终点搜索区域单独进行低通滤波，再次消除肌电干扰和工频干扰；步骤六：低平t波样本所有心搏及叠加波t波终点搜索区域的重构；步骤七：重构后的心搏t波终点搜索区域聚类分析；步骤八：心搏t波终点搜索区域波形特征提取；样本所有重构后的心搏t波终点搜索区域进行特征提取，描述每个心搏搜索区域数据有意义的属性，该特征集包括了形态特征群、统计特征群、信息特征群;步骤九：模型训练；根据上述步骤八提取的心搏t波终点搜索区域的形态特征群、统计特征群、信息特征群，结合对应的t波终点到t波顶点的距离标签，将训练集中叠加波相似的心搏、不相似的心搏，分别输入到随机森林模型中进行训练，得到两种t波终点定位模型；步骤十：模型评估。2.根据权利要求1所述的心电图低电压t波终点定位方法，其特征在于：所述步骤二中的4个等级干扰的识别方法如下：首先将预处理后的心电样本数据按照n秒钟长度切片，再进行间隔m个点的下采样，获得心电数据x，并对其进行标注得到y；其次，对心电数据x进行傅里叶变换f，并将其实部标准化norm之后和x相加，得到输出xf；再次，训练卷积神经网络，设计一个包含四层卷积层的一维卷积神经网络来提取特征，同时在每第二个卷积层之后加一个maxpooling层用来筛选特征，并规定每一层合适的通道数c，卷积核长度k和卷积步长strides，随后接一个averagepool层和两层全连接层来解码提取到的特征，全连接层的通道数分别为c和4，第二个全连接层的输出为xfc，xfc通过softmax层，得到4个类别概率的一维向量xout，设定学习率lr、损失函数loss，用adam优化器进行训练，直到模型收敛；最后，导入测试数据进行4个等级干扰的识别，对比最终输出的4个概率xp，取其中最高概率对应的等级作为n秒钟的干扰等级，公式如下：
。3.根据权利要求1所述的心电图低电压t波终点定位方法，其特征在于：所述步骤三中的四个阶段分别为：第一阶段，在对样本l个心搏进行聚类前，首先需要产生模板心搏，将目标心搏与模板心搏进行相似度距离计算，相似度距离满足一定阈值条件的目标心搏匹配归入相应模板类别，选取相似度距离阈值th，从样本的l个心搏中随机抽取ls个心搏，并选择第一个心搏作为第一个模板心搏类型，计算第二个心搏到的相似度距离，如果满足＞th，则将匹配归入该类别，否则以为新建模板心搏，继续计算下一个心搏与已有的若干个模板心搏间相似度距离(j=0,1,2,
…
)，若满足，且＞th，则将该心搏匹配归于模板类别， 否则以为新建模板心搏类别，重复此步骤，直至完成随机ls个心搏的聚类，最后对随机的ls个心搏模板分类结果进行统计排序，选取心搏匹配个数最多的模板类别，将ls个心搏中属于的所有心搏数据做平均，得到一个平均模板心搏；第二阶段，对样本l个心搏第一次进行聚类，获取心搏匹配个数最多的一类，将上一阶段得到的平均模板心搏作为第一个心搏模板类型，重复第一阶段的心搏匹配模板的步骤，直到所有l个心搏匹配完毕，最后对全部l个模板心搏分类结果进行统计排序，得到心搏匹配个数最多的模板类别，获取其匹配的心搏集合及个数；第三阶段，样本中全部l个心搏去掉上一阶段最多一类个心搏后，再对剩余个的心搏第二次进行聚类，对个心搏重复第一阶段步骤，得到随机抽取心搏匹配个数最多的模板类别，接着对个心搏重复第二阶段步骤，得到匹配个数最多的模板类别，获取其匹配心搏集合及其个数，将第二阶段和第三阶段获得的心搏集合、合并为一个心搏集合；第四阶段，针对每个心电样本构建独特的心搏叠加波，将上述阶段获取的样本心搏集合用来构建叠加波，所有用于构建的心搏以r波位置为基准点向前向后对齐，并对每个
心搏采样点对应的振幅幅值进行排序，截取第15百分位至第85百分位对应的心搏幅值求和取平均，作为当前采样点对应的叠加波振幅幅值，以此来构建叠加波，公式如下：其中m是当前心搏的第m个采样点，和分别是每个心搏当前采样点对应的振幅幅值排序之后处于在15%和85%的位置，是叠加波第m个采样点的幅值，是第m个采样点处于第i百分位的心搏幅值。4.根据权利要求1所述的心电图低电压t波终点定位方法，其特征在于：所述步骤四的具体方法为：针对样本所有心搏及上一步骤构建的样本叠加波，进行t波顶点的检测，选取t波顶点幅值阈值 ，筛选叠加波t波顶点绝对值最大值小于的心电样本，定义为低平t波样本，共筛选出份训练集和份测试集的低平t波单导联ecg样本，其中包括了窦性节律样本、室上性心律失常样本、室性心律失常样本、房扑房颤样本、传导阻滞样本，然后以样本各心搏t波顶点为基准点，设定对应心搏低平t波终点的搜索区域为：。5.根据权利要求1所述的心电图低电压t波终点定位方法，其特征在于：所述步骤五的具体方法为： 无限脉冲响应滤波器iir对低平t波终点搜索区域低通滤波，具体设计了巴特沃斯数字低通滤波器，通带和阻带的截止频率为和，通带允许最大衰减是，阻带允许最大衰减是，iir滤波器差分方程的一般表达式如下：其中x(n)为输入序列，y(n)为输出序列，和为滤波器系数，m为滤波器阶数；最后，将搜索区域的低平t波序列按顺序滤波，将所得结果在时域翻转后反向通过滤波器，再将所得结果在时域翻转后输出，使得相移为0，保证滤波后的波形位置精确地保留在滤波前序列中的位置。6.根据权利要求1所述的心电图低电压t波终点定位方法，其特征在于：所述步骤六中对低平t波样本所有心搏及叠加波t波终点搜索区域的重构方法为： 首先将上一步骤处理后的样本心搏及其叠加波的t波终点搜索区域数据进行归一化处理，使其映射到合理区间：其中x、为搜索区域数据映射前和映射后的值，分别为搜索区域数据最大值和最小值，a、b为映射区间的上限和下限；其次，通过e的指数运算将归一化后的数据进行非线性变换，从而使搜索区域的低平t波数据得到增强，最后，为使增强后的t波终点搜索区域具有更好的噪声抑制性能,对样本
心搏及其叠加波增强后的t波终点搜索区域数据进行sg平滑滤波的重构，原理是基于搜索区域特征的多项式拟合，应用最小二乘法确定加权系数进行移动窗口加权平均滤波，具体公式如下：其中x是原始搜索序列，是重构后的搜索序列，是滑动窗口的第i个时间序列的相关系数，2m+1是滑动窗口的大小。7.根据权利要求1所述的心电图低电压t波终点定位方法，其特征在于：所述步骤七中重构后的心搏t波终点搜索区域聚类分析方法为：以步骤六重构得到的叠加波t波终点搜索区域为初始模板，将重构后的样本l个心搏t波终点搜索区域，重复步骤三第一阶段的心搏匹配模板的步骤，直到样本所有个心搏t波终点搜索区域匹配完毕，获取其匹配的心搏t波终点搜索区域集合及其补集；通过对低平t波的份训练集和份测试集样本区分叠加波相似的心搏和不相似的心搏，在后续步骤分开进行t波终点的建模和预测。8.根据权利要求1所述的心电图低电压t波终点定位方法，其特征在于所述步骤八中心搏t波终点搜索区域波形特征提取，具体提取的特征群为：1）形态特征群一阶差分：心搏t波终点搜索区域数据相邻两点之差；二阶差分：心搏t波终点搜索区域数据一阶差分相邻两点之差；心搏t波顶点幅值；rr间期；心搏r波顶点到t波顶点的距离；心搏r波顶点到t波顶点的距离和rr间期的比值；2）统计特征群心搏t波终点搜索区域数据平均值；心搏t波终点搜索区域数据标准差；心搏t波终点搜索区域数据峭度；心搏t波终点搜索区域数据偏度；心搏t波终点搜索区域数据平方的平均值；心搏t波终点搜索区域数据最大值和最小值的差值；一阶差分数据绝对值最小值位置距离t波顶点位置的距离二阶差分数据绝对值最小值位置距离t波顶点位置的距离3）信息特征群心搏t波终点搜索区域数据模糊熵；心搏t波终点搜索区域数据功率谱熵；心搏t波终点搜索区域数据奇异谱熵；心搏t波终点搜索区域数据排列熵。9.根据权利要求1所述的心电图低电压t波终点定位方法，其特征在于：所述步骤十中
的模型评估方法为：将训练好的两种低平t波定位模型应用在测试集的样本心搏数据上，测试集样本中的心搏分别应用其所属叠加波相似和不相似两种类型模型进行预测，将最后的t波终点位置输出结果与测试集样本心搏的t波终点位置标签进行对比，误差间隔是用来衡量波形定位的误差区间，当算法检测结果和医生标注位置位于误差区间范围内时，检测结果被认可，使用误差间隔内的样本占比评估模型的性能，各项评估指标的表达式如下：k：模型输出的低平t波终点在标注位置的
±
4ms范围内的样本数；l：模型输出的低平t波终点在标注位置的
±
5ms~
±
8ms范围内的样本数；m：模型输出的低平t波终点在标注位置的
±
9ms~
±
16ms范围内的样本数；n：模型输出的低平t波终点在标注位置的
±
16ms外的样本数；误差范围
±
4ms以内的样本占比：误差范围
±
8ms以内的样本占比：误差范围大于
±
16ms的样本占比：其中误差间隔在本发明低平t波终点定位算法中与采样点关系如下：误差间隔4ms=1个采样点误差间隔8ms=2个采样点误差间隔16ms=4个采样点。

技术总结
一种心电图低电压T波终点定位方法，所述方法包括以下步骤：步骤一：单导联心电信号获取及预处理，根据预设的采样率，采用干电极获取t秒固定长度的单导联ECG样本，并对每个样本进行滤波去肌电、工频干扰和去基线漂移处理；步骤二：心电样本数据4个等级干扰的识别；本发明在于对低平T波非线性变换的放大重构，增强了低平T波终点的细节特征，根据构建的样本心搏叠加波，结合样本心搏和叠加波的重构操作，聚类分析成两类心搏，增加的两种状态引入了医学信息，解决了两类心搏数据对应的提取特征高度相似问题，最后提取形态特征群、统计特征群、信息特征群分别建模预测，在实际应用场景下实现了对心搏低平T波终点的准确定位。现了对心搏低平T波终点的准确定位。现了对心搏低平T波终点的准确定位。


技术研发人员：许昆明 符灵建 刘凯歌 杨茂
受保护的技术使用者：浙江好络维医疗技术有限公司
技术研发日：2023.07.14
技术公布日：2023/8/14