一种动态心电图中的运动伪差识别方法及装置制造方法及图纸

本发明专利技术公开一种动态心电图中的运动伪差识别方法及装置。首先根据检测得到的心拍，提取每个心拍的特征参数；然后，计算心拍前后固定时间宽度内的噪声情况，计算得到心拍的信噪比以及置信度等指数；接着，根据设定阈值，将信噪比或置信度过低的心拍识别为运动伪差；最后，结合每个心拍的前后心拍的伪差情况，将连续伪差中未成功识别的心拍视为运动伪差，最终得到运动伪差识别结果。本发明专利技术方法解决现有无法快速准备识别心拍运动伪差问题，在已检测的心拍基础上，根据心拍特征参数提取出能反映噪声水平的信噪比以及置信度指数，并结合前后心拍噪声情况进一步识别运动伪差情况，最终自动剔除无效心拍。

A method and device for identifying motion pseudo difference in dynamic electrocardiogram

【技术实现步骤摘要】
一种动态心电图中的运动伪差识别方法及装置
本专利技术属于医学测试
，公开了一种动态心电图中的运动伪差识别方法及装置。
技术介绍
动态心电图是连续记录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毫秒)。本文档来自技高网...

【技术保护点】
1.一种动态心电图中的运动伪差识别方法，其特征在于该方法包含以下各步骤：/n(1)在已完成心拍检测的动态心电图基础上，提取每个心拍的特征参数；所述特征参数包括面积A、宽度W、R波高度H

【技术特征摘要】
1.一种动态心电图中的运动伪差识别方法，其特征在于该方法包含以下各步骤：
(1)在已完成心拍检测的动态心电图基础上，提取每个心拍的特征参数；所述特征参数包括面积A、宽度W、R波高度H1、Q波高度H2、S波高度H3等；
(2)根据已检测得到的心拍，获取每个心拍前后200毫秒宽度区域的平均面积AT；
(3)根据上述步骤(1)每个心拍的特征参数以及步骤(2)每个心拍前后200毫秒宽度区域的平均面积AT计算心拍信噪比以及置信度等指数，根据满足以下任一条件则认为该心拍为运动伪差：①信噪比小于阈值1.5，则认为运动伪差；②置信度小于阈值40，则认为运动伪差；
心拍信噪比SNR指数计算公式：



心拍置信度指数Rconf计算公式如下：
若心拍宽度W＜60毫秒或心拍宽度W＞200毫秒，则Rconf＝0；
若心拍高度H1＜(H2+H3)，则Rconf＝0；
若60毫秒≤心拍宽度W≤200毫秒，且心拍高度H1≥(H2+H3)，则



(4)从第一个心拍起，若当前心拍为非运动伪差，则计...

【专利技术属性】
技术研发人员：黄超，汤征，李顶立，
申请(专利权)人：苏州百慧华业精密仪器有限公司，
类型：发明
国别省市：江苏;32