一种基于多尺度斜率特征的可除颤节律识别方法与流程

1.本发明涉及心电信号处理技术领域和除颤设备领域，具体涉及一种基于多尺度斜率特征的可除颤节律识别方法。


背景技术：

2.心室颤动是严重的致命性心律失常，各部分心肌快速而不协调的颤动，心室丧失有效收缩功能，不能提供有效的血液循环。因此，对室颤及快速室性心动过速的检测显得尤为重要。
3.目前主要应用的可除颤节律检测方法是从心电信号中提取出有利于区分可除颤节律与不可除颤节律的特征，然后设计规则进行判断或者利用机器学习方法进行决策。特征提取的有效性是算法准确性的关键，常用的特征有时域特征（如幅值、斜率、信号分布、心率等）、频域特征（如频谱分布等）、复杂度特征（如复杂度、峭度、近似熵、样本熵等）。
4.现有的可除颤节律识别方法的准确度较差，对于环境复杂多变、要求算法识别速度快的场合中，目前的算法性能还有待提高。


技术实现要素：

5.本发明为了解决上述现有技术中的技术问题，提出一种基于多尺度斜率特征的可除颤节律识别方法。
6.本发明采用的技术方案是：本发明提出了一种基于多尺度斜率特征的可除颤节律识别方法，包括步骤：步骤1，获取预设长度的心电信号，将获取到的心电信号进行预处理；步骤2，将预处理后的心电信号以不同的差分间隔进行计算，得到不同尺度的斜率特征；步骤3，将不同尺度的斜率特征汇总得到特征向量；步骤4，根据特征向量判断是否为可除颤节律。
7.步骤2具体包括：步骤21，将预处理后的心电信号与不同的差分间隔代入公式进行计算，得到每个差分间隔对应的斜率信号；步骤22，将每个差分间隔对应的斜率信号归一化；步骤23，对每个归一化后的斜率信号中的斜率值进行统计，得到对应每一个斜率信号的多个斜率特征。
8.进一步的，步骤21的公式为：；为差分间隔为r的情况下第n个点的斜率值，k为预设的滑动平均点数，所述心电信号选取的待计算的点的总数为n，即每个所述斜率信号包括n个斜率值。
9.进一步的，步骤22具体包括：步骤221，将斜率信号平均划分为三个片段，选出每个片段中的最大斜率值；步骤222，将三个片段中的最大斜率值按大小依次排序，选出位于中间的最大斜率值作为mid_max_slope；步骤223，用mid_max_slope对斜率信号进行归一化，即：；为归一化后的斜率值。
10.进一步的，步骤23具体包括：步骤231，设置阈值th；步骤232，统计每个斜率信号中小于阈值th的斜率值所占的比例，作为特征，，n为斜率信号的斜率值总数；步骤233，统计斜率信号中跨过阈值th的次数，作为特征，。
11.优选地，所述步骤4具体为：将步骤3获取的特征向量，输入训练好的机器学习分类模型进行分类，根据分类结果决定该心电信号对应的心电片段是否为可除颤节律。
12.上述机器学习分类模型的训练方法为：步骤51，从预设心电数据库中，选取多个心电图片段，每个心电图片段的长度为步骤1中的预设长度；步骤52，对各个心电片段按照步骤1至3进行处理，获得对应的特征向量；步骤53，将步骤52得到的特征向量与预设的是否为可除颤节律信息输入机器学习分类模型进行训练，输出是否为可除颤节律的分类结果。
13.优选地，机器学习分类模型具体为神经网络、支持向量机或者随机森林。
14.进一步的，预处理的方式为：将预设长度的心电信号进行滤波处理得到滤波后的心电信号。
15.进一步的，所述不同的差分间隔为1，2，4，8...。
16.与现有技术比较，本发明具有以下创新点：1、采用多种差分间隔尺度的斜率特征分析，可以提取到不同类型心律波形的特征，有利于准确区分各种类型的可除颤心律和不可除颤心律。
17.2、采用分段最大斜率值的中值对斜率信号进行归一化，可以有效减轻尖峰等突变信号对特征提取的影响，提高算法在噪声干扰环境下的稳定性3、提供的算法计算简单、所需心电信号片段时间短，可以快速准确地识别心电信号是否为可除颤节律。
附图说明
18.为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些
实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
19.图1为本发明实施例中的流程图；图2为本发明实施例中归一化的流程图。
具体实施方式
20.为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
21.下面结合附图以及实施例对本发明的原理及结构进行详细说明。
22.如图1所示，本发明提出了一种基于多尺度斜率特征的可除颤节律识别方法，包括步骤：步骤1，获取预设长度的心电信号，将获取到的心电信号进行预处理；该预设长度可以是4秒，或者4秒左右的数值。
23.步骤2，将预处理后的心电信号以不同的差分间隔进行计算，得到不同尺度的斜率特征；通过将预处理后的心电信号采用多种间隔尺度进行计算，可以提取到不同类型心律波形的特征，从而有利于准确区分各种类型的可除颤心律和不可除颤心律。
24.步骤3，将多个尺度的斜率特征汇总得到特征向量；步骤4，根据特征向量判断是否为可除颤节律。
25.在具体的实施例中，预处理可以是进行30hz低通滤波以及1hz高通滤波，将心电信号重采样到100hz左右，得到滤波后的心电信号, 长度为l。
26.不同的间隔尺度具体可以选择为1、2、4、8。或者1、2、4、8、16...（数值到达），不同的差分间隔可以提取到不同类型心律信号的特征。例如正常窦性心律qrs波变化快，差分间隔小的斜率变化较灵敏；而室性心搏的波形较宽大，差分间隔大的斜率变化比较敏感。现以间隔尺度为1、2、4、8为例具体说明步骤2，步骤2具体包括：步骤21，将预处理后的心电信号与不同的差分间隔代入公式进行计算，得到每个差分间隔对应的斜率信号；具体为：在心电信号中一共选取n个点，计算每个点的斜率值，即差分间隔为1时能够计算出一组的斜率值，这组斜率值的集合为差分间隔为1时对应的斜率信号。即差分间隔为2时也能够计算出一组的斜率值，这组斜率值的集合为差分间隔为2时对应的斜率信号。具体公式如下：；为差分间隔为r的情况下第n个点的斜率值；k为预设的滑动平均点数，优选取100ms心电长度对应的点数。例如当采样率为100时，k = 100 * 0.1 = 10。
27.步骤22，将每个斜率信号归一化得到归一化后的斜率信号；如图2所示，步骤22具体包括：
步骤221，将斜率信号平均划分为三个片段，选出每个片段中的最大斜率值；即将一个斜率信号中所有的斜率值均分成三份，然后再选取每份中的最大斜率值；步骤222，将三个片段中的最大斜率值按大小依次排序，选出位于中间的最大斜率值作为mid_max_slope；步骤223，用mid_max_slope对斜率信号进行归一化，即对斜率信号中的每一个斜率值进行归一化，具体为：，为归一化后的斜率值；现有技术方案一般采用最大值进行归一化，但是最大值容易受到尖峰等突变干扰信号影响。本发明采用信号分段最大值的中值进行归一化，可以提高抗干扰性。
28.步骤23，对每个归一化后的斜率信号中的斜率值进行统计，得到对应每一个斜率信号的多个斜率特征。
29.其中步骤23具体包括：步骤231，设置阈值th；步骤232，统计每个斜率信号中小于阈值th的斜率值所占的比例，作为特征，，其中n为斜率信号的斜率值总数；步骤233，统计斜率信号中跨过阈值th的次数，作为特征，。
30.即一个斜率信号对应一个和一个。
31.步骤3，将多个尺度的斜率特征汇总得到特征向量具体为；；此外，也可以将其他现有技术所提取的特征纳入到该特征向量中。
32.步骤4具体为：将步骤3获取的特征向量，输入训练好的机器学习分类模型进行分类，根据分类结果决定该心电信号对应的心电片段是否为可除颤节律，也可以采用现有技术中的其他模型进行分类判断。其中机器学习分类模型的具体训练步骤包括：步骤51，从预设心电数据库中，选取多个心电图片段，每个心电图片段的长度为步骤1中的预设长度；步骤52，对各个心电片段按照步骤1至3进行处理，获得对应的特征向量；步骤53，将步骤42得到的特征向量与预设的是否为可除颤节律信息输入机器学习分类模型进行训练，输出是否为可除颤节律的分类结果。
33.机器学习分类模型可以选择神经网络、支持向量机、随机森林等，预设心电数据库具体为心电专家标注的心电数据库。
34.以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。