一种动态心电图中的运动伪差识别方法及装置与流程

本发明属于医学测试技术领域，公开了一种动态心电图中的运动伪差识别方法及装置。

背景技术：

动态心电图是连续记录24小时或更长时间内的体表心电图，其包含的信息量远远大于常规心电图，从而增加了心电图分析的工作量，这就使得心电图的自动检测与分析技术成为必要。并且在心电图采集过程中，容易受到日常活动的各种干扰，尤其是运动干扰的影响，导致实际采集的心电波形掺杂着各种运动伪差干扰，这些伪差多为无效信号，需要剔除分析，若仅靠手工剔除，无疑大大增加了医生的工作量。因此，快速准确的运动伪差干扰识别方法是提高医生工作效率的有效手段。

由于动态心电图的诊断均以心拍为单位进行，因而运动伪差也一般以心拍为基础进行识别。目前，市场上已有的产品对运动伪差的识别均以直方图或者心拍模板辅助手工删除为主，处理效率普遍偏低。

技术实现要素：

本发明的目的是针对现有无法快速、准确识别心电图中运动伪差信号的不足，提供一种动态心电图中的运动伪差识别方法，以达到使运动伪差识别更加精确、可靠，有效的完成对运动伪差的识别检测的目的。

本发明的动态心电图中的运动伪差识别方法，包括以下步骤：

(1)在已完成心拍检测的动态心电图基础上，提取每个心拍的特征参数；所述特征参数包括面积a、宽度w、r波高度h1、q波高度h2、s波高度h3等；

(2)根据已检测得到的心拍，获取每个心拍前后200毫秒宽度区域的平均面积at；

(3)根据上述步骤(1)每个心拍的特征参数以及步骤(2)每个心拍前后200毫秒宽度区域的平均面积at计算心拍信噪比以及置信度等指数，根据满足以下任一条件则认为该心拍为运动伪差：①信噪比小于阈值1.5，则认为运动伪差；②置信度小于阈值40，则认为运动伪差；

心拍信噪比指数snr计算公式如下：

若信噪比snr小于1.5，则视检测心拍为运动伪差；

心拍置信度指数rconf计算公式如下：

若心拍宽度w小于60毫秒或大于200毫秒，则rconf设为0；

若心拍高度h1小于(h2+h3)，则rconf设为0；

若置信度rconf小于40，则认为该心拍伪差干扰过大，将其判定为运动伪差；

(4)从第一个心拍起，若当前心拍为非运动伪差，则计算其前后n秒内的心拍伪差个数占比，若伪差比例大于一定阈值，则将当前心拍视为运动伪差；

(5)重复步骤(4)两次，完成运动伪差识别处理。

一种动态心电图中的运动伪差识别装置，包括：

特征提取单元，用于提取已完成心拍检测的动态心电图中每个心拍的特征参数；所述特征参数包括面积a、宽度w、r波高度h1、q波高度h2、s波高度h3等；

第一计算单元，用于获取每个心拍前后200毫秒宽度区域的平均面积at；at为t1点至q波以及s波至t2点的幅值积分除以宽度(400毫秒)；

第二计算单元，用于计算每个心拍的心拍信噪比和置信度；

第一运动伪差识别单元，用于初步判断当前心拍是否为运动伪差；具体是满足以下任一条件则认为该心拍为运动伪差：①信噪比小于阈值1.5，则认为运动伪差；②置信度小于阈值40，则认为运动伪差；

第二运动伪差识别单元，用于再次判断经第一运动伪差识别单元识别后的非运动伪差心拍是否为运动伪差；具体是通过计算当前心拍前后n秒内的心拍伪差个数占比，若伪差比例大于一定阈值，则将当前心拍视为运动伪差；

本发明方法解决现有无法快速准备识别心拍运动伪差问题，在已检测的心拍基础上，根据心拍特征参数提取出能反映噪声水平的信噪比以及置信度指数，并结合前后心拍噪声情况进一步识别运动伪差情况，最终自动剔除无效心拍。

附图说明

图1是典型的心电波形以及各特征点和特征参数标注，其中，横向虚线表示基线，w表示心拍宽度，h1/h2/h3分别表示r波/q波/s波高度，心拍面积a为q波至s波的幅值积分除以心拍宽度w，心拍前后200毫秒宽度的平均面积at为t1点至q波以及s波至t2点的幅值积分除以宽度(400毫秒)。

具体实施方式

下面结合附图对本发明作进一步说明：

动态心电图中的运动伪差识别方法，步骤如下：

(1)在已完成心拍检测的动态心电图基础上，提取每个心拍的特征参数，所述特征参数包括心拍面积a、心拍宽度w、r波高度h1、q波高度h2、s波高度h3等；

(2)根据已检测得到的心拍，计算每个心拍前后200毫秒宽度的平均面积at；

(3)根据计算的各项特征参数提取心拍信噪比以及置信度等指数，其中，心拍信噪比指数snr计算公式如下：

若信噪比snr小于1.5，则视检测心拍为运动伪差；

心拍置信度指数rconf计算公式如下：

进一步地，若心拍宽度w小于60毫秒或大于200毫秒，则rconf设为0；

进一步地，若心拍高度h1小于(h2+h3)，则rconf设为0；

若置信度rconf小于40，则将检测心拍判定为运动伪差；

(4)从第一个心拍起，若当前心拍为非运动伪差，则计算其前后n秒内的心拍伪差个数占比，若伪差比例大于70％，则将当前心拍视为运动伪差；

(5)重复步骤(4)两次，完成运动伪差识别处理。

一种动态心电图中的运动伪差识别装置，包括：

特征提取单元，用于提取已完成心拍检测的动态心电图中每个心拍的特征参数；所述特征参数包括面积a、宽度w、r波高度h1、q波高度h2、s波高度h3等；

第一计算单元，用于获取每个心拍前后200毫秒宽度区域的平均面积at；at为t1点至q波以及s波至t2点的幅值积分除以宽度(400毫秒)；

第二计算单元，用于计算每个心拍的心拍信噪比和置信度；

第一运动伪差识别单元，用于初步判断当前心拍是否为运动伪差；具体是满足以下任一条件则认为该心拍为运动伪差：①信噪比小于阈值1.5，则认为运动伪差；②置信度小于阈值40，则认为运动伪差；

第二运动伪差识别单元，用于再次判断经第一运动伪差识别单元识别后的非运动伪差心拍是否为运动伪差；具体是通过计算当前心拍前后n秒内的心拍伪差个数占比，若伪差比例大于一定阈值，则将当前心拍视为运动伪差；

图1是典型的心电波形以及各特征点和特征参数标注，其中，横向虚线表示基线，w表示心拍宽度，h1/h2/h3分别表示r波/q波/s波高度，心拍面积a为q波至s波的幅值积分除以心拍宽度w，心拍前后200毫秒宽度的平均面积at为t1点至q波以及s波至t2点的幅值积分除以宽度(400毫秒)。