一种分析含伪差心电图的方法和装置与流程

本发明属于心电图技术领域，具体地涉及一种通过衍生导联分析含伪差心电图的方法和装置。

背景技术：

心电图主要反映心脏激动的电学活动，因此对各种心率失常和传导障碍的诊断分析具有肯定价值，到目前为止，尚没有任何其它方法能替代心电图在这方面的作用。心电图记录的是心脏激动时所产生于体表的电位差，不管是采用3导联、5导联还是12导联采集心电信号，都可能包含多种不是由心脏激动而发生于心电图上的改变的干扰，称为伪差。伪差通常是由强磁场、静电干扰、患者运动造成的导联线拉扯或肌电干扰等因素产生。尤其是动态心电图，因佩戴后患者需要正常活动，产生伪差几乎是不可避免的。心电图伪差不仅没有利用价值，还会对医生诊断造成影响，降低心电监护准确性。

现有不少技术用来识别伪差，如公开专利：CN104644160A，但其将伪差段数据排除在分析范围外，另外也有一些技术用来从多个导联中选择干扰最小的导联作为分析导联，如公开专利：CN105411567A，以提高分析精度。但是如果伪差发生时间较长，或者肢体导联影响了Wilson中心电端造成每个导联伪差都比较大，就会严重影响其分析结果的准确性。

技术实现要素：

本发明的目的在于为解决上述问题而提供一种通过衍生导联分析含伪差心电图的方法和装置，提高了含伪差心电图的分析准确性，且计算简单，鲁棒性和适应性好。

为实现上述目的，本发明采用的技术方案为：一种分析含伪差心电图的方法，包括以下步骤：

S1，按时间段识别出心电图中含显著伪差部分和不含显著伪差部分；

S2，筛选出受干扰小的心电电极集合SET；

S3，将含显著伪差部分根据S2步骤所筛选出的心电电极集合SET生成衍生导联L_d；

S4，从衍生导联L_d中选取质量最好的一条或者多条衍生导联L_d进行进一步处理。

进一步的，所述步骤S2中，是通过对不含显著伪差部分的分析处理中筛选出受干扰小的心电电极组合SET。

更进一步的，所述步骤S2具体为：从心电电极中选取所有不重复的任意两个电极来组成第一过渡心电电极对集合SET′，将不含显著伪差部分按第一过渡心电电极对集合SET′生成第一过渡衍生导联L′_d1；并根据该第一过渡衍生导联L′_d1设立一个阈值Threshold1，选取第一过渡衍生导联L′_d1幅度大于阈值Threshold1的所有电极对，形成心电电极集合SET。

进一步的，所述步骤S2中，是通过对含显著伪差部分和不含显著伪差部分相配合的分析处理中筛选出受干扰小的心电电极组合SET。

更进一步的，所述步骤S2具体为：从心电电极中选取所有不重复的任意两个电极来组成第一过渡心电电极对集合SET′，将不含显著伪差部分按第一过渡心电电极对集合SET′生成第一过渡衍生导联L′_d1；并根据该第一过渡衍生导联L′_d1设立一个阈值Threshold1，选取第一过渡衍生导联L′_d1幅度小于阈值Threshold1的所有电极对，形成第二过渡心电电极集合SET′1；

将含显著伪差部分按第二过渡心电电极集合SET′1来生成第二过渡衍生导联L′_d2；选取第二过渡导联L′_d2幅度仍小于阈值Threshold1的电极对，然后从第一过渡心电电极对集合SET′中找出与这对电极相关的所有元素，构成心电电极集合SET。

进一步的，所述步骤S4具体为选取质量最好的一条或者多条衍生导联L_d进行显示处理。

更进一步的，所述显示处理还包括：同时显示该衍生导联L_d反映的是哪两个电极之间的电位差。

进一步的，将含显著伪差部分根据S2步骤所筛选出的心电电极集合SET生成衍生导联L_d具体方法为：计算心电电极集合SET中的每个元素组间的电位差产生衍生导联L_d。

进一步的，将不含显著伪差部分按第一过渡心电电极对集合SET′生成第一过渡衍生导联L′_d1具体方法为：计算第一过渡心电电极对集合SET′中的每个元素组间的电位差产生第一过渡衍生导联L′_d1。

进一步的，所述步骤S2中的阈值Threshold1采用自适应方法设立。

更进一步的，所述步骤S2中的阈值Threshold1采用自适应方法设立具体为：计算心搏数据在第一过渡衍生导联L′_d1上的幅度，选所有幅度的中位数或者平均值作为阈值threshold1。

进一步的，所述步骤S4具体为：计算各条衍生导联L_d心电信号的能量，并按心电信号能量占总能量的比重排序；选取心电信号能量比重最大的一条或者多条衍生导联L_d进行进一步处理。

进一步的，还包括如下步骤；将不含显著伪差部分的心电图按步骤S3和S4的方式生成衍生导联数据，以利于对比判别。

本发明还公开了一种分析含伪差心电图的装置，包括

识别模块，用于按时间段识别出心电图中含显著伪差部分和不含显著伪差部分；

筛选模块，用于筛选出受干扰小的心电电极集合SET；

衍生导联生成模块，用于将含显著伪差部分根据筛选模块所筛选出的心电电极集合SET生成衍生导联L_d；

处理模块，用于从衍生导联L_d中选取质量最好的一条或者多条衍生导联L_d进行进一步处理。

本发明的有益技术效果：

本发明根据心电图采集使用的导联体系，筛选出受干扰小的心电电极组合SET，将含显著伪差部分根据心电电极组合SET生成衍生导联L_d，并将质量最好的一条或多条衍生导联L_d通过显示界面呈现给阅图者，或者输出给自动分析程序做进一步分析，实现了对含伪差信号的心电图分析，提高了其分析的准确性，适用于多种导联体系(如Wilson，EASI，15导，18导等)，且计算简单，鲁棒性和适应性好。

附图说明

图1为本发明实施例的方法流程图；

图2为本发明实施例的装置结构示意图；

图3为本发明实施例一的方法流程图；

图4为本发明实施例二的方法流程图；

图5为本发明实施例二的含显著伪差部分的心电图；

图6为本发明实施例二的处理后的心电图。

具体实施方式

现结合附图和具体实施方式对本发明进一步说明。

如图1所示，一种分析含伪差心电图的方法，包括以下步骤：

S1，按时间段识别出心电图中含显著伪差部分和不含显著伪差部分；

S2，筛选出受干扰小的心电电极集合SET；

S3，将含显著伪差部分根据S2所筛选出的心电电极集合SET生成衍生导联L_d；

S4，从衍生导联L_d中选取质量最好的一条或者多条衍生导联L_d进行进一步处理。

本发明还公开了一种分析含伪差心电图的装置，包括

识别模块1，用于按时间段识别出心电图中含显著伪差部分和不含显著伪差部分；

筛选模块2，用于筛选出受干扰小的心电电极集合SET；

衍生导联生成模块3，用于将含显著伪差部分根据筛选模块所筛选出的心电电极集合SET生成衍生导联L_d；

处理模块4，用于从衍生导联L_d中选取质量最好的一条或者多条衍生导联L_d进行进一步处理。

下面将以具体实施例来说明本发明。

心电图机采用Wilson标准12导联就可以拾取心脏电活动的全貌，不仅能从额面，也能从横面观察心脏电信息活动，已得到国际上的认可而在世界各国通用。

根据Wilson标准12导心电图连接方式，考虑到肢体导联的其他数据都可以通过II，III导联推导得出，因此定义导联向量L＝(Ⅱ Ⅲ V1 V2 V3 V4 V5 V6)^T，即独立导联数为8，电极电势向量E＝(E_F E_RA E_LA E_C1 E_C2 E_C3 E_C4 E_C5 E_C6)^T，根据Wilson导联的计算方式：

E_WCT＝1/3*(E_F+E_RA+E_LA)

II＝E_F-E_RA

III＝E_F-E_LA

V_i＝E_Ci-E_WCT

可以推导出导联向量转换为电极电压向量的转换矩阵T_L2E，如下：

其中E_WCT为Wilson中心电端，可见三个肢体导联中只要任何一个肢体电极产生干扰，胸导都会受到不同程度的影响，而如果两个肢体电极遭遇干扰，则整个12导联都会受到严重影响，导致心电图基本不可读。

下面就以Wilson标准12导联心电图为例介绍本发明，但不以此为限。

实施例一

如图3所示，一种分析含伪差心电图的方法，包括如下步骤：

A1，识别并按时间段区分心电图中含显著伪差部分和不含显著伪差部分。

具体的，采用上述的Wilson标准12导联采集心电图信号，然后识别并按时间段区分心电图中含显著伪差部分(简称为伪差段)和不含显著伪差部分(简称为参考段)。本实施例中，采用现有技术区分出伪差段，例如峰值峰谷对法等方法，也可以采用公开专利：CN104644160A的方法进行识别，具体可以参照现有技术，此不再细说。

A2，将心电图导联数据转换成心电电极的电压值。

本具体实施例中，将心电图导联数据转换成心电电极的电压值具体方法为：E＝T_L2E*L，其中，本实施例中计算转换矩阵T_L2E时以左下肢F为零电势点，故上述的转换矩阵T_L2E中第一行全为0。在其它实施例中也可以取其它导联为零电势点，获得不同的转换矩阵T_L2E，方法相同，不再一一说明。

A3，从心电电极中选取所有不重复的任意两个电极来组成第一过渡心电电极对集合SET′，第一过渡心电电极对集合SET′中的每个元素为一个心电电极对。

本具体实施例中，从9个电极任意选取2个电极为一对，形成第一过渡心电电极对集合SET′，第一过渡心电电极对集合SET′中有即36个元素。计算每个元素组间的电位差可以产生36种衍生导联。例如选择第一过渡心电电极对集合SET′中的元素(C1，C5)得到一个没有肢体导联参与的新衍生导联V_C1-C5≡E_C1-E_C5。

A4，将参考段的电极电压值，按第一过渡心电电极对集合SET′中的每个元素生成第一过渡衍生导联L′_d1；并根据第一过渡衍生导联L′_d1设立一个阈值Threshold1，选取第一过渡衍生导联L′_d1幅度大于阈值Threshold1的所有元素，形成心电电极集合SET。

具体的，针对不同的患者，有些电极之间的电位差比较小，不利于识别，需要先排除。为此，本具体实施例通过将参考段的心电数据生成36种第一过渡衍生导联L′_d1，并计算心搏数据在各条第一过渡衍生导联L′_d1上的幅度。选取幅度大于阈值threshold1的所有第一过渡衍生导联L′_d1的所有元素，形成心电电极集合SET。阈值threshold1可以选所有幅度的中位数或者平均值。

A5，将伪差段按心电电极集合SET中的元素生成衍生导联L_d，衍生导联L_d的数目和心电电极集合SET中的元素个数相同，生成方法为：计算心电电极集合SET中每个元素组间的电位差。

A6，计算伪差段生成的各条衍生导联L_d的心电信号的质量。

本具体实施例中，计算伪差段生成的各条衍生导联L_d的心电信号的能量，并按心电信号能量占总能量的比重排序，心电信号能量比重越大，则说明衍生导联L_d的质量越好。计算方法采用现有的计算方法，此是本领域技术人员可以轻易实现的，不再细说。

A7，选取质量最好的一条或者多条衍生导联L_d进行处理分析。

具体的，选取质量最好(心电信号能量比重最大)的一条或者多条衍生导联L_d显示。同时显示该衍生导联L_d反映的是哪两个电极之间的电位差，方便阅图者识别。

进一步的，为了方便对比，本发明中，还可以将参考段的心电图按按步骤A5和A6的同样的组合生成衍生导联数据，以利于对比判别，提高分析准确性。

实施例二

如图4所示，一种分析含伪差心电图的方法，包括如下步骤：

A1，识别并按时间段区分心电图中含显著伪差部分和不含显著伪差部分。

具体的，采用上述的Wilson标准12导联采集心电图信号，然后识别并按时间段区分心电图中含显著伪差部分(简称为伪差段)和不含显著伪差部分(简称为参考段)。本实施例中，采用现有技术区分出伪差段，例如峰值峰谷对法等方法，也可以采用公开专利：CN104644160A的方法进行识别，具体可以参照现有技术，此不再细说。

A2，将心电图导联数据转换成心电电极的电压值。

本具体实施例中，将心电图导联数据转换成心电电极的电压值具体方法为：E＝T_L2E*L，其中，本实施例中计算转换矩阵T_L2E时以左下肢F为零电势点，故上述的转换矩阵T_L2E中第一行全为0。在其它实施例中也可以取其它导联为零电势点，获得不同的转换矩阵T_L2E，方法相同，不再一一说明。

A3，从心电电极中选取所有不重复的任意两个电极来组成第一过渡心电电极对集合SET′，第一过渡心电电极对集合SET′中的每个元素为一个心电电极对。

本具体实施例中，从9个电极任意选取2个电极为一对，形成第一过渡心电电极对集合SET′，第一过渡心电电极对集合SET′中有即36个元素。计算每个元素组间的电位差可以产生36种衍生导联。例如选择第一过渡心电电极对集合SET′中的元素(C1，C5)得到一个没有肢体导联参与的新衍生导联V_C1-C5≡E_C1-E_C5。

A4，将参考段的的电极电压值，按第一过渡心电电极对集合SET′中的每个元素生成第一过渡衍生导联L′_d1；并根据该第一过渡衍生导联L′_d1设立一个阈值Threshold1，选取第一过渡衍生导联L′_d1幅度小于阈值Threshold1的所有元素，形成第二过渡心电电极集合SET′1。

本具体实施例通过将参考段的心电数据生成36种第一过渡衍生导联L′_d1，并计算心搏数据在各条第一过渡衍生导联L′_d1上的幅度。选取幅度小于阈值threshold1的所有第一过渡衍生导联L′_d1的元素，形成第二过渡心电电极集合SET′1。阈值threshold1可以选所有幅度的中位数或者平均值。由于第二过渡心电电极集合SET′中元素生成的衍生导联幅度小，便于判断该导联是否受干扰。

A5，将伪差段按第二过渡心电电极集合SET′1中的元素来生成第二过渡衍生导联L′_d2；第二过渡衍生导联L′_d2的数目和第二过渡心电电极集合SET′1中的元素个数相同，生成方法为：计算第二过渡心电电极集合SET′1中每个元素组间的电位差。选取第二过渡衍生导联L′_d2幅度仍小于阈值Threshold1的电极对，然后从第一过渡心电电极对集合SET′中找出与这对电极相关的所有元素，构成心电电极集合SET。例如，从第二过渡心电电极集合SET′1中选中电极对(C1，C5)，则心电电极集合SET包含第一过渡心电电极对集合SET′中具有C1或C5的所有元素。

A6，将伪差段按照心电电极集合SET中的元素来生成衍生成衍生导联L_d，计算衍生导联L_d的心电信号的质量。

本具体实施例中，计算伪差段生成的各条衍生导联L_d的心电信号的能量，并按心电信号能量占总能量的比重排序，心电信号能量比重越大，则说明衍生导联L_d的质量越好。计算方法采用现有的计算方法，此是本领域技术人员可以轻易实现的，不再细说。

A7，选取质量最好的一条或者多条导联进行处理分析。

具体的，选取质量最好(心电信号能量比重最大)的一条或者多条衍生导联L_d显示。同时显示该衍生导联L_d反映的是哪两个电极之间的电位差，方便阅图者识别。

进一步的，为了方便对比，本发明中，还可以将参考段的心电图按按步骤A5和A6的同样的组合生成衍生导联数据，以利于对比判别，提高分析准确性。

图4为伪差段的心电图，按照上述方法分析处理后如图5所示，图5表明C1电极和C5电极受干扰较小，同时此显示结果也便于医生和心电分析软件分析，达到了衍生导联辅助分析含伪差段心电图的效果。

本发明根据心电图采集使用的导联体系，筛选出受干扰小的心电电极组合，将伪差段按照这些心电电极组合生成衍生导联L_d，并将质量最好的一条或多条衍生导联L_d通过显示界面呈现给阅图者，或者输出给自动分析程序做进一步分析，实现了对含伪差信号的心电图分析，提高了其分析的准确性，适用于多种导联体系，且计算简单，鲁棒性和适应性好。

当然，上述方法也可以用于其它导联体系，如EASI，15导，18导等，只需对转换矩阵、第一过渡心电电极对集合SET′等做相应的简单调整即可实现。此是本领域技术人员可以轻易实现的，不再详细说明。

上述实施例给出了两种筛选出受干扰小的心电电极集合SET的方法，当然，也可以采用现有的其它方法来实现，此是本领域技术人员可以轻易实现的，不再详细说明。

尽管结合优选实施方案具体展示和介绍了本发明，但所属领域的技术人员应该明白，在不脱离所附权利要求书所限定的本发明的精神和范围内，在形式上和细节上可以对本发明做出各种变化，均为本发明的保护范围。