一种基于MODWT的动态心电图实时心率估算方法与流程

本发明提出了一种基于modwt的动态心电图实时心率估算方法，涉及心电图智能诊断领域。

背景技术：

心电图检查是体检的常见项目，如果患者疑似有心律失常症状，通常也会去医院做个心电图，但是，在医院的几分钟或者十几分钟心电图，难以发现问题，这时候医生可能会给患者一个可以穿戴的心电图检测器，佩戴一天、两周或者更长时间，这会产生几十、几百小时的心电图，即所谓的动态心电图，这就会出现两个问题：其一，为了及时发现致命心脏疾患，心电图检测器最好能实时进行分析，给出重要心脏疾患指标；其二：医生要一秒一秒的检查长时间的心电图，非常耗费时间，为此，一条比较好的思路是，首先由心电图检测器实时的给出重要的心脏指标，以便及时预警，其次是由心电图检测器进行智能判断，找出可能的心律失常的时间位置，再由医生针对这些地方进行进一步诊断，显然，对动态心电图信号进行实时检测、给出辅助诊断结果的首要任务是检测心率的变化，即迅速、准确地识别出主波位置，本专利的发明初衷来源于此。

与常规心电图相比较，动态心电图更容易受到外界因素的干扰，信号中存在大量的工频干扰、基线漂移和肌电噪声，目前，对心电信号的噪声抑制一般采用低通滤波方法，其中主要的不足在于固定频带的滤波不能适应心电信号的随机性、非平稳性，小波变换具有良好的时-频局部分析特点，能较好的解决这些问题，但在强噪声背景下，进行小波变换后，不具有平移不变性，且幅度较小的q、s波处常出现gibbs振荡现象，造成信号失真,为此，我们提出一种modwt的心电信号自适应滤波方法，modwt是一种高度冗余的非正交小波变换，对样本容量没有要求，具有小波系数和尺度系数的平移不变性、无相位扭曲，从而保留了心电信号中有用的特征波形，避免了有用心电信息的损失。

就主波的定位而言，已有不少研究者从多个角度出发对心电图主波的识别和定位进行了研究，比如差分阈值法、小波分析法、模板匹配法、聚类方法、神经网络(深度学习等)等,因动态1秒心率更新的需要，与常规心电图主波定位相比，其要求的运行速度很快，为此，我们在进行modwt的基础上，结合心电图各特征波频带分布、幅值特点及不应期等相关心电规则，选择了主波所在的主要频带，给出了一种快速的主波定位算法，并估算心率。

该方法基于modwt及心电规则，对动态心电图实时心率进行了估算，易于理解，容易实现，精度较高。

技术实现要素：

本发明克服现有技术不足，设计一种基于modwt的动态心电图实时心率估算方法。

为了实现上述目的，本发明所采用的技术方案如下：

一种基于modwt的动态心电图实时心率估算方法，包括：a)预处理；b)基于modwt及各特征波形频带特点对动态心电图进行滤波处理；c)基于心电规则，确定动态心电图主波位置及实时心率，所述a)预处理方法为：

步骤1：采集心电图信号，采样频率为128hz，以5秒、10秒、20、30秒作为分析时长分别进行了分析测试，基于测试结果，确定动态心电图分析实时参数，对心电图信号进行分段处理；

步骤2：去除干扰程度较大的分段心电信号；

作为优选：所述b)基于modwt及各特征波形频带特点对动态心电图进行滤波处理方法为：

步骤1：基于采样频率及各特征波形频带特点对选定分段信号进行modwt处理；

步骤2：去除干扰频带及p波、t波所在频带，选定分析频带；

步骤3：对分析频带心电信号进行自适应滤波，所述自适应滤波方法为极值阈值滤波方法：

，当n＞32

xlim＝0，当n<＝32

其中，xlim为阈值，n为信号长度，σ为均方根，可由中位数绝对偏差(mad)计算获得：

mad＝median(|x-median(x)|

σ＝1.4826×mad

其中，x即为选定的某分割频带信号；

作为优选：所述c)基于心电规则，确定动态心电图主波位置及实时心率方法为：

步骤1：基于心电规则，确定分析频带主波位置：

(1)找出信号的所有极值点：对信号做一阶差分运算，差分值改变的位置即为信号极值(注：因主波位置可能为可能朝上为r波，也可能朝下为q波、s波或qs波，所以可能是极大值，也可能是极小值)；

(2)基于心室不应期，选择不应期范围内绝对幅值较大的极值点：心室有效不应期约为200-300ms，我们以250ms为界，初步选定其中绝对幅值最大的极值点为预选的的主波位置；

(3)去除绝对幅值过小的主波位置：当某处预选主波对应的幅值小于最大幅值的1/4时，去除该处主波；

(4)比较当前信号及前5秒信号得出的主波位置(注：初始5秒不进行校准)，对主波位置进行自校准：基本原则是当某处主波在越多处得到验证时，则越确定为主波位置，反之则相反；

步骤2：估算实时心率：

根据主波位置，计算主波间期，当某间期时长大于2秒时，排除该间期，对剩下间期求平均间期，基于采样频率由下式估算心率：

心率＝60×采样频率/平均主波间期；

附图说明

图1为1份截取的60秒动态心电数据心率估算结果；

具体实施方式

下面结合实例，对本发明的技术方案作进一步详细的说明；

一种基于modwt的动态心电图实时心率估算方法，包括以下步骤：

a)预处理：

步骤1：基于测试结果，确定动态心电图分析实时参数，对心电图信号进行分段处理；

采集心电图信号，采样频率为128hz，以5秒、10秒、20、30秒作为分析时长分别进行了分析测试，当以10秒作为分析时长时，正常主波个数10-15，有利于后期根据平均主波间期来估算心率，也能兼顾实时分析的要求，决定采用10秒作为分析时长，另考虑心率估算实时更新的要求，以1秒间隔作为更新时长；

步骤2：去除干扰程度较大的分段心电信号；

当某段10秒心电信号幅值特别大(比如大于5mv)时，显然为干扰，此时进行modwt滤波及定位是没有意义的，反倒容易造成有用特征波形的被掩盖，为此，舍弃该段信号，心率维持为前1秒心率；

b)基于modwt及各特征波形频带特点对动态心电图进行滤波处理：

步骤1：基于采样频率及各特征波形频带特点对选定分段信号进行modwt处理。

根据心电规则，正常心电信号频率范围约为0.01-100hz，其中90％的信号能量集中在0.5-45hz，qrs波群约为3-45hz，波峰位置约在8-16hz，p波、t波约为0.5-10hz，而就干扰信号而言，工频干扰为50hz，因为电极移动、人体呼吸等引起的基线漂移约为0.05-2hz，肌电干扰约为5-2000hz，显然，就心率估算及相对应的主波定位而言，需要尽可能的保留qrs波群中的主波所在频带，而去掉其他一切杂波；

与普通的离散小波变换类似，modwt也具有频带二分特点，最高分析频率为采样频率/2，所以，当采样频率已知，可逐次对较低频带(注：第一次分割的频带为0-采样频率/2)进行二分,考虑到后期主要目的是进行主波定位，可将频带分割到某层次较低频带的最高频率达到个位数(即<10hz)，此时,频带分割完毕；

步骤2：去除干扰频带及p波、t波所在频带，选定分析频带；

因最终层次较低频带主要为p波、t波、基线漂移或低频率的肌电干扰，30-40hz以上较高频带主要为工频干扰及高频率的肌电干扰，直接去除这些频带，此时剩余的分割频带即为待重构频带信号，也是主波主要所在频带；

步骤3：对分析频带心电信号进行自适应滤波；

在剩余频带中，可能还会跟主波频带相似的肌电干扰，但其幅值相对于主波而言会比较小，此时，在此基础上采用自适应滤波方法，为了尽可能的保留有用信号，我们选择较保守的极值阈值滤波方法：

，当n＞32

xlim＝0，当n<＝32

其中，xlim为阈值，n为信号长度，σ为均方根，可由中位数绝对偏差(mad)计算获得：

mad＝median(|x-median(x)|

σ＝1.4826×mad

其中，x即为选定的某分割频带信号；

在滤波之后，再进行逆modwt进行重构即可，从而获得处理过的分析信号。

c)基于心电规则，确定动态心电图主波位置及实时心率：

步骤1：基于心电规则，确定主波位置；

主要分为4步：

(1)找出信号的所有极值点：对信号做一阶差分运算，差分值改变的位置即为信号极值(注：因主波位置可能为可能朝上为r波，也可能朝下为q波、s波或qs波，所以可能是极大值，也可能是极小值)；

(2)基于心室不应期，选择不应期范围内绝对幅值较大的极值点：心室有效不应期约为200-300ms，我们以250ms为界，初步选定其中绝对幅值最大的极值点为预选的的主波位置；

(3)去除绝对幅值过小的主波位置：当某处预选主波对应的幅值小于最大幅值的1/4时，去除该处主波；

(4)比较当前信号及前5秒信号得出的主波位置(注：初始5秒不进行校准)，对主波位置进行自校准：基本原则是当某处主波在越多处得到验证时，则越确定为主波位置，反之则相反。

步骤2：估算实时心率；

根据主波位置，计算主波间期，当某间期时长大于2秒时，排除该间期，对剩下间期求平均间期，基于采样频率由下式估算心率：

心率＝60×采样频率/平均主波间期

实施例

实施例说明

为了检验方法的有效性，我们截取1份60秒动态心电数据进行了心率估算测试。采样频率：250hz；

如前所述，我们以10秒作为分析信号，1秒进行1次心率更新。modwt分到第4层，取d4频带信号作为分析频带，基于极值阈值滤波方法进行信号自适应滤波，并基于心室不应期、心电图特征波形幅值特点等心电规则，估算出实时心率。

计算结果

图1(a)为原始心电信号；

图1(b)为滤波后心电信号及主波位置。其中主波位置以红色“ο”标出；

图1(c)为心率估算趋势图(注意：从10秒位置开始)；

由图1可见，在干扰比较严重的情况下，主波位置定位准确，心率估算结果在55±1次/分之间,符合要求；

以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。