CN113208599B - 一种基于心电单极导联v1、v5r推衍导联v3r、v4r的方法及设备 - Google Patents

一种基于心电单极导联v1、v5r推衍导联v3r、v4r的方法及设备 Download PDF

Info

Publication number
CN113208599B
CN113208599B CN202110532230.2A CN202110532230A CN113208599B CN 113208599 B CN113208599 B CN 113208599B CN 202110532230 A CN202110532230 A CN 202110532230A CN 113208599 B CN113208599 B CN 113208599B
Authority
CN
China
Prior art keywords
lead
data
clinical data
electrocardio
linear regression
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110532230.2A
Other languages
English (en)
Other versions
CN113208599A (zh
Inventor
张斌
周一彬
段扬
张永红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Pengyang Fengye Technology Co ltd
Original Assignee
Beijing Pengyang Fengye Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Pengyang Fengye Technology Co ltd filed Critical Beijing Pengyang Fengye Technology Co ltd
Priority to CN202110532230.2A priority Critical patent/CN113208599B/zh
Publication of CN113208599A publication Critical patent/CN113208599A/zh
Application granted granted Critical
Publication of CN113208599B publication Critical patent/CN113208599B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Abstract

本发明公开一种基于心电单极导联V1、V5R推衍导联V3R、V4R的方法及设备，属于心电导联技术领域。该方法包括以下步骤：获取心电单极导联V1的临床数据V1和心电单极导联V5R的临床数据V5R；对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行基于临近单极导联线性相关的推衍处理，得到导联V3R的推衍数据V3 ′ R和导联V4R的推衍数据V4 ′ R；对所述导联V3R的推衍数据V3 ′ R和所述导联V4R的推衍数据V4 ′ R分别进行幅值转化处理，得到所述导联V3R的合成数据V3R和所述导联V4R的合成数据V4R，以便利用所述导联V3R的合成数据V3R和所述导联V4R的合成数据V4R进行心电图分析。

Description

一种基于心电单极导联V1、V5R推衍导联V3R、V4R的方法及 设备
技术领域
本发明涉及心电导联技术领域，尤其涉及一种基于心电单极导联V1、V5R推衍导联V3R、V4R的方法及设备。
背景技术
动态心电图是通过动态心电图仪在患者日常生活状态下连续24小时或更长时间记录其心电活动的全过程，并借助计算机进行分析处理，以发现在常规体表心电图检查时不易发现的心律失常和心肌缺血等，为临床诊断、治疗及判断疗效提供重要的客观依据。动态心电图仪又称Holter心电图，已由单导、双导发展为12导联、18导联全记录。
随着导联增加，心电数据更加全面，利于心脏诊断。然而，导联增加，导联电极数量增加，且患者需要佩戴24小时或更长时间，不利于患者活动，且舒适性降低。因此，在现有的心电算法中，一般采用推衍导联的方式，以减少动态心电图实际所需的导联电极数量。
导联电极及安放位置如下：V1，电极安放于胸骨右缘第四肋间；V2，电极安放于胸骨左缘第四肋间；V3，电极安放于V2～V4连接线的中间位置；V4，电极安放于左锁骨中线与第五肋间相交处；V5，电极安放于左腋前线和V4同一水平的位置；V6，电极安放于左腋中线和V4、V5同一水平的位置；V7，电极安放于左腋后线和V5～V6同一水平的位置；V8，电极安放于左肩胛线和V5～V6同一水平的位置；V9，电极安放于左脊柱旁线和V5～V6同一水平的位置；V3R：电极安放于与V3相对应的右侧胸壁处；V4R：电极安放于与V4相对应的右侧胸壁处；V5R：电极安放于与V5相对应的右侧胸壁处。
目前，已有心电算法较少能推衍或准确推衍V3R和V4R。
发明内容
本发明实施例的主要目的在于提供一种基于心电单极导联V1、V5R推衍导联V3R、V4R的方法及设备，旨在推衍特别是准确推衍V3R和V4R，以进一步减少导联电极。
本发明实施例提供了一种基于心电单极导联V1、V5R推衍导联V3R、V4R的方法，所述方法包括以下步骤：获取心电单极导联V1的临床数据V1和心电单极导联V5R的临床数据V5R；对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行基于临近单极导联线性相关的推衍处理，得到导联V3R的推衍数据V′3R和导联V4R的推衍数据V′4R；对所述导联V3R的推衍数据V′3R和所述导联V4R的推衍数据V′4R分别进行幅值转化处理，得到所述导联V3R的合成数据V3R和所述导联V4R的合成数据V4R，以便利用所述导联V3R的合成数据V3R和所述导联V4R的合成数据V4R进行心电图分析。
优选地，所述对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行基于临近单极导联线性相关的推衍处理，得到所述导联V3R的推衍数据V′3R和所述导联V4R的推衍数据V′4R包括：利用用于推衍所述导联V3R的第一线性回归方程，对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行推衍处理，得到所述导联V3R的推衍数据V′3R，其中，所述第一线性回归方程为V′3R＝β3R1*V1+β3R2*V5R+β3R0，所述β3R1，β3R2，β3R0是所述第一线性回归方程的参数；利用用于推衍所述导联V4R的第二线性回归方程，对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行推衍处理，得到所述导联V4R的推衍数据V′4R，其中，所述第二线性回归方程为V4R＝β4R1*V1+β4R2*V5R+β4R0，其中，β4R1，β4R2，β4R0是所述第二线性回归方程的参数。
优选地，所述第一线性回归方程的参数β3R1、β3R2、β3R0的参数值通过以下步骤得到：获取第一时间段内所述心电单极导联V1的已有临床数据、所述心电单极导联V3R的已有临床数据和所述心电单极导联V5R的已有临床数据；利用所述第一时间段内的所述心电单极导联V1的已有临床数据、所述心电单极导联V3R的已有临床数据和所述心电单极导联V5R的已有临床数据，确定所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值；对所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β3R1、β3R2、β3R0的最终的参数值。
优选地，所述对所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β3R1、β3R2、β3R0的最终的参数值包括：根据所述第一时间段内所述心电单极导联V1的已有临床数据，确定所述心电单极导联V1的已有临床数据平均值V1amp；根据所述第一时间段内所述心电单极导联V5R的已有临床数据，确定所述心电单极导联V5R的已有临床数据平均值V5Ramp；利用所述心电单极导联V1的已有临床数据平均值V1amp和所述心电单极导联V5R的已有临床数据平均值V5Ramp，对所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β3R1、β3R2、β3R0的最终的参数值。
优选地，所述第二线性回归方程的参数β4R1、β4R2、β4R0的参数值通过以下步骤得到：获取第二时间段内所述心电单极导联V1的已有临床数据、所述心电单极导联V4R的已有临床数据和所述心电单极导联V5R的已有临床数据；利用所述第二时间段内的所述心电单极导联V1的已有临床数据、所述心电单极导联V4R的已有临床数据和所述心电单极导联V5R的已有临床数据，确定所述第二线性回归方程的参数β4R1、β4R2、β4R0的初始的参数值；对所述第二线性回归方程的参数β4R1、β4R2、β4R0的初始参的数值进行归一化处理，得到所述第二线性回归方程的参数β4R1、β4R2、β4R0的最终的参数值。
优选地，所述对所述第一线性回归方程的参数β4R1、β4R2、β4R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β4R1、β4R2、β4R0的最终的参数值包括：根据所述第二时间段内所述心电单极导联V1的已有临床数据，确定所述心电单极导联V1的已有临床数据平均值V1amp；根据所述第二时间段内所述心电单极导联V5R的已有临床数据，确定所述心电单极导联V5R的已有临床数据平均值V5Ramp；利用所述心电单极导联V1的已有临床数据平均值V1amp和所述心电单极导联V5R的已有临床数据平均值V5Ramp，对所述第一线性回归方程的参数β4R1、β4R2、β4R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β4R1、β4R2、β4R0的最终的参数值。
优选地，所述对所述导联V3R的推衍数据V′3R和所述导联V4R的推衍数据V′4R分别进行幅值转化处理，得到所述导联V3R的合成数据V3R和所述导联V4R的合成数据V4R包括：获取推衍时间段内所述心电单极导联V1的临床数据，并根据所述推衍时间段内的所述心电单极导联V1的临床数据，确定所述推衍时间段内所述心电单极导联V1的临床数据平均值V1AMP；获取所述推衍时间段内所述心电单极导联V5R的临床数据，并根据所述推衍时间段内的所述心电单极导联V5R的临床数据，确定所述推衍时间段内所述心电单极导联V5R的临床数据平均值V5RAMP；利用所述推衍时间段内所述心电单极导联V1的临床数据平均值V1AMP和所述心电单极导联V5R的临床数据平均值V5RAMP，对所述推衍时间段内所述导联V3R的推衍数据V′3R进行幅值转化处理，得到所述导联V3R的合成数据V3R；利用所述推衍时间段内所述心电单极导联V1的临床数据平均值V1AMP和所述心电单极导联V5R的临床数据平均值V5RAMP，对所述推衍时间段内所述导联V4R的推衍数据V′4R进行幅值转化处理，得到所述导联V4R的合成数据V4R。
本发明实施例还提供了一种基于心电单极导联V1、V5R推衍导联V3R、V4R的设备，所述设备包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的程序，所述程序被所述处理器执行时实现上述的基于心电单极导联V1、V5R推衍导联V3R、V4R的方法的步骤。
本发明提供的一种基于心电单极导联V1、V5R推衍导联V3R、V4R的方法及设备，其通过心电单极导联V1的临床数据V1和心电单极导联V5R的临床数据V5R，得到导联V3R、V4R的合成数据，减少了导联电极数量，且合成数据的准确性高。
附图说明
图1是本发明提供的基于心电单极导联V1、V5R推衍导联V3R、V4R的方法的流程图；
图2是本发明提供的基于心电单极导联V1、V5R推衍导联V3R、V4R的详细流程图；
图3a和图3b分别是本发明提供的V1原始数据、V5R原始数据对应的波形图以及V3R原始数据、V3R合成数据对应的波形图；
图4a和图4b分别是本发明提供的V1原始数据、V5R原始数据对应的波形图以及V4R原始数据、V4R合成数据对应的波形图；
图5是本发明提供的基于心电单极导联V1、V5R推衍导联V3R、V4R的设备的结构框图。
具体实施方式
应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
在后续的描述中，使用用于表示元件的诸如“模块”、“部件”或“单元”的后缀仅为了有利于本发明的说明，其本身没有特有的意义。因此，“模块”、“部件”或“单元”可以混合地使用。
图1是本发明提供的基于心电单极导联V1、V5R推衍导联V3R、V4R的方法的流程图，如图1所示，该方法可以包括以下步骤：
步骤S101：获取心电单极导联V1的临床数据V1和心电单极导联V5R的临床数据V5R。
心电单极导联V1的临床数据V1是安放于胸骨右缘第四肋间的电极采集的心电数据，心电单极导联V5R的临床数据V5R是安放于与V5相对应的右侧胸壁处的电极采集的心电数据，其中，V5对应的电极安放于左腋前线和V4同一水平的位置，V4对应的电极安放于左锁骨中线与第五肋间相交处。
本发明旨在利用心电单极导联V1的临床数据V1和心电单极导联V5R的临床数据V5R。准确地推衍出导联V3R的数据和导联V4R的数据，以下通过步骤S102和步骤S103进行具体说明。
步骤S102：对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行基于临近单极导联线性相关的推衍处理，得到导联V3R的推衍数据V′3R和导联V4R的推衍数据V′4R。
具体地，通过以下步骤得到所述导联V3R的推衍数据V′3R：利用用于推衍所述导联V3R的第一线性回归方程，对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行推衍处理，得到所述导联V3R的推衍数据V′3R，其中，所述第一线性回归方程为V′3R＝β3R1*V1+β3R2*V5R+β3R0，所述β3R1，β3R2，β3R0是所述第一线性回归方程的参数。
上述β3R1，β3R2，β3R0的参数值可以通过以下步骤得到：获取第一时间段内所述心电单极导联V1的已有临床数据、所述心电单极导联V3R的已有临床数据和所述心电单极导联V5R的已有临床数据；利用所述第一时间段内的所述心电单极导联V1的已有临床数据、所述心电单极导联V3R的已有临床数据和所述心电单极导联V5R的已有临床数据，确定所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值；对所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β3R1、β3R2、β3R0的最终的参数值。具体地说，根据所述第一时间段内所述心电单极导联V1的已有临床数据，确定所述心电单极导联V1的已有临床数据平均值V1amp，根据所述第一时间段内所述心电单极导联V5R的已有临床数据，确定所述心电单极导联V5R的已有临床数据平均值V5Ramp，然后利用所述心电单极导联V1的已有临床数据平均值V1amp和所述心电单极导联V5R的已有临床数据平均值V5Ramp，对所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β3R1、β3R2、β3R0的最终的参数值，例如，先计算V1amp与V5Ramp之和的一半，得到调整系数，将参数β3R1、β3R2、β3R0的初始的参数值除以该调整系数，得到参数β3R1、β3R2、β3R0的最终的参数值。
同样地，通过以下步骤得到所述导联V4R的推衍数据V′4R：利用用于推衍所述导联V4R的第二线性回归方程，对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行推衍处理，得到所述导联V4R的推衍数据V′4R，其中，所述第二线性回归方程为V′4R＝β4R1*V1+β4R2*V5R+β4R0，其中，β4R1，β4R2，β4R0是所述第二线性回归方程的参数。
上述β4R1，β4R2，β4R0的参数值可以通过以下步骤得到：获取第二时间段内所述心电单极导联V1的已有临床数据、所述心电单极导联V4R的已有临床数据和所述心电单极导联V5R的已有临床数据；利用所述第二时间段内的所述心电单极导联V1的已有临床数据、所述心电单极导联V4R的已有临床数据和所述心电单极导联V5R的已有临床数据，确定所述第二线性回归方程的参数β4R1、β4R2、β4R0的初始的参数值；对所述第二线性回归方程的参数β4R1、β4R2、β4R0的初始参的数值进行归一化处理，得到所述第二线性回归方程的参数β4R1、β4R2、β4R0的最终的参数值。具体地说，根据所述第二时间段内所述心电单极导联V1的已有临床数据，确定所述心电单极导联V1的已有临床数据平均值V1amp，根据所述第二时间段内所述心电单极导联V5R的已有临床数据，确定所述心电单极导联V5R的已有临床数据平均值V5Ramp，然后利用所述心电单极导联V1的已有临床数据平均值V1amp和所述心电单极导联V5R的已有临床数据平均值V5Ramp，对所述第一线性回归方程的参数β4R1、β4R2、β4R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β4R1、β4R2、β4R0的最终的参数值。例如，先计算V1amp与V5Ramp之和的一半，得到调整系数，将参数β4R1、β4R2、β4R0的初始的参数值除以该调整系数，得到参数β3R1、β3R2、β3R0的最终的参数值。
需要指出的是，上述第一时间段和第二时间段可以是同一时间段，也可以是不同时间段。上述心电单极导联V1的已有临床数据、心电单极导联V3R的已有临床数据、心电单极导联V4R的已有临床数据和心电单极导联V5R的已有临床数据为历史数据，用于确定上述各个参数的参数值。确定参数值时，可以利用上述第一/第二时间段的全部数据，确定一组参数值，也可以将上述第一/第二时间段划分成若干个子时间段，然后利用各个子时间段的数据，确定回归效果较好的一组或多组参数值，当存在回归效果较好的多组参数值时，可以选取回归效果最好的一组参数值作为最终参数值，也可以对回归效果较好的多组参数值进行求平均处理，并将参数平均值作为最终参数值。这样，在确定上述各个参数的最终参数值后，基于步骤S101获取的临床数据，通过步骤S102和步骤S103的处理即可得到导联V3R的数据和导联V4R的数据。
步骤S103：对所述导联V3R的推衍数据V′3R和所述导联V4R的推衍数据V′4R分别进行幅值转化处理，得到所述导联V3R的合成数据V3R和所述导联V4R的合成数据V4R，以便利用所述导联V3R的合成数据V3R和所述导联V4R的合成数据V4R进行心电图分析。
具体地，步骤S103包括：获取推衍时间段内所述心电单极导联V1的临床数据，并根据所述推衍时间段内的所述心电单极导联V1的临床数据，确定所述推衍时间段内所述心电单极导联V1的临床数据平均值V1AMP；获取所述推衍时间段内所述心电单极导联V5R的临床数据，并根据所述推衍时间段内的所述心电单极导联V5R的临床数据，确定所述推衍时间段内所述心电单极导联V5R的临床数据平均值V5RAMP；利用所述推衍时间段内所述心电单极导联V1的临床数据平均值V1AMP和所述心电单极导联V5R的临床数据平均值V5RAMP，对所述导联V3R的推衍数据V′3R进行幅值转化处理，得到所述导联V3R的合成数据V3R；利用所述推衍时间段内所述心电单极导联V1的临床数据平均值V1AMP和所述心电单极导联V5R的临床数据平均值V5RAMP，对所述导联V4R的推衍数据V′4R进行幅值转化处理，得到所述导联V4R的合成数据V4R。其中，推衍时间段是需要进行推衍合成的导联V3R、V4R的数据的时间段。
本发明提供的基于心电单极导联V1、V5R推衍导联V3R、V4R的方法，通过心电单极导联V1的临床数据V1和心电单极导联V5R的临床数据V5R，得到导联V3R、V4R的合成数据，减少了导联电极数量，且合成数据的准确性高。
在本领域中，还可以提供一种计算机可读存储介质，该计算机可读存储介质存储有程序，该程序可被处理器执行，以实现上述基于心电单极导联V1、V5R推衍导联V3R、V4R的方法的各个步骤。
图2是本发明提供的基于心电单极导联V1、V5R推衍导联V3R、V4R的详细流程图，如图2所示，推衍导联V3R的步骤可以包括：
步骤S201：设线性回归方程。
对心电单极导联V1、V3R、V5R的数据V1、V″3R、V5R进行多元线性回归分析，设方程为：
V″3R＝β1*V1+β2*V5R+β0
用最小二乘法计算出公式中对应的β1、β2、β0三个参数；
其中，n是数据的长度；V11，V12…V1n是心电单极导联V1的n个数据；V3R1，V3R2…V3Rn是心电单极导联V3R的n个数据；V5R1，V5R2…V5Rn是心电单极导联V5R的n个数据。
步骤S202：对已有的临床数据进行测试，得到多组β1、β2、β0的值，选取其中回归效果较好的几组进行分析。
将利用一组参数值推导出来的心电单极导联V3R的数据与实际测试得到的心电单极导联V3R的数据比较，若两者相似度高于预设相似度阈值，则认为该组参数值的回归效果较好，可以选取该组参数。实施时，可以将最高相似度对应的一组阈值应用于后续计算步骤，也可以对回归效果较好的几组参数值进行求平均处理，并将平均值应用于后续计算步骤。
步骤S203：因为不同的心电数据的幅值不同，为了计算统一需要对参数进行归一化处理：
β3R1＝β1/((V1amp+V5Ramp)*0.5)
β3R2＝β2/((V1amp+V5Ramp)*0.5)
β3R0＝β0/((V1amp+V5Ramp)*0.5)
V1amp是求参数时对应心电数据V1导联幅值的平均值，V5Ramp是求参数时对应心电数据V5R导联幅值的平均值。
将已测得的参数和数据进行综合分析，得出最终的参数值，例如，得到的最终参数值分别为：β3R1＝-0.01；β3R2＝0.445；β3R0＝0.261。
步骤S204：测试数据时，在归一化的方程上乘以一定的倍数。
对待测数据进行幅值转化，即在归一化后的方程上乘以V1和V5R导联幅值的平均值相加的一半，具体如下：
V3R＝(V1AMP+V5RAMP)*0.5*(β3R1*V1+β3R2*V5R+β3R0)
其中，V1AMP是待测心电数据中V1导联幅值的平均值，V5RAMP是待测心电数据中V5R导联幅值的平均值。
步骤S205：计算出导联V3R的最终数据。
即，获取V1和V5R，并根据V1和V5R，计算V′3R＝β3R1*V1+β3R2*V5R+β3R0，然后乘以(V1AMP+V5RAMP)*0.5，得到导联V3R的合成数据V3R。
如图3a和图3b所示，心电单极导联V3R的原始波形与由心电单极导联V1和V5R合成的导联V3R的合成波形相似度高，两者基本一致。
推衍导联V4R与推衍导联V3R的方法基本相同，如图2所示，可以包括如下步骤：
步骤S201：设线性回归方程。
1.对心电单极导联V1、V4R、V5R的数据V1、V″4R、V5R进行多元线性回归分析，设方程为：
V″4R＝β1*V1+β2*V5R+β0
用最小二乘法计算出公式中对应的β1、β2、β0三个参数；
其中，n是数据的长度；V11，V12…V1n是心电单极导联V1的n个数据；V4R1，V4R2…V4Rn是心电单极导联V4R的n个数据；V5R1，V5R2…V5Rn是心电单极导联V5R的n个数据。
步骤S202：对已有的临床数据进行测试，得到多组β1、β2、β0的值，选取其中回归效果较好的几组进行分析。
将利用一组参数值推导出来的心电单极导联V4R的数据与实际测试得到的心电单极导联V4R的数据比较，若两者相似度高于预设相似度阈值，则认为该组参数值的回归效果较好，可以选取该组参数。实施时，可以将最高相似度对应的一组阈值应用于后续计算步骤，也可以对回归效果较好的几组参数值进行求平均处理，并将平均值应用于后续计算步骤。
步骤S203：因为不同的心电数据的幅值不同，为了计算统一需要对参数进行归一化处理：
β4R1＝β1/((V1amp+V5Ramp)*0.5)
β4R2＝β2/((V1amp+V5Ramp)*0.5)
β4R0＝β0/((V1amp+V5Ramp)*0.5)
V1amp是求参数时对应心电数据V1导联幅值的平均值，V5Ramp是求参数时对应心电数据V5R导联幅值的平均值。
将已测得的参数和数据进行综合分析，得出最终的参数值，例如，得到的最终参数值分别为：β4R1＝-0.046；β4R2＝0.425；β4R0＝0.225。
步骤S204：测试数据时，在归一化的方程上乘以一定的倍数。
对待测数据进行幅值转化，即在归一化后的方程上乘以V1和V5R导联幅值的平均值相加的一半，具体如下：
V4R＝(V1AMP+V5RAMP)*0.5*(β4R1*V1+β4R2*V5R+β4R0)
其中，V1AMP是待测心电数据中V1导联幅值的平均值，V5RAMP是待测心电数据中V5R导联幅值的平均值。
步骤S205：计算出导联V4R的最终数据。
即，获取V1和V5R，并根据V1和V5R，计算V′4R＝β4R1*V1+β4R2*V5R+β4R0，然后乘以(V1AMP+V5RAMP)*0.5，得到导联V4R的合成数据V4R。
如图4a和图4b所示，心电单极导联V4R的原始波形与由心电单极导联V1和V5R合成的导联V4R的合成波形基本一致。
在现有的心电算法中，所有的推衍导联都是基于双极的，如基于Wilson推导的肢体导联，或基于Frank的正交心电图的导联，本发明提出了一种基于心电单极导联以及利用临近心电单极导联之间线性相关推衍出对应导联的方法，提高了人体在扭转运动等不同姿态时导联重建的准确性。另外，传统心电图的判别标准都是基于Wilson导联体系的，本发明的基于心电单极导联的推衍方法使重建后的波形更加接近Wilson导联，便于心电图的判读，同时减少了电极数量。
图5是本发明提供的基于心电单级导联V1、V5R推衍导联V3R、V4R的设备的结构框图，如图5所示，所述设备20包括存储器21、处理器22以及存储在所述存储器21上并可在所述处理器22上运行的程序，所述程序被所述处理器22执行时实现上述的基于心电单级导联V1、V5R推衍导联V3R、V4R的方法的步骤。
本发明提供的基于心电单极导联V1、V5R推衍导联V3R、V4R的设备，通过心电单极导联V1的临床数据V1和心电单极导联V5R的临床数据V5R，得到导联V3R、V4R的合成数据，减少了导联电极数量，且合成数据的准确性高。
以上参照附图说明了本发明的优选实施例，并非因此局限本发明的权利范围。本领域技术人员不脱离本发明的范围和实质内所作的任何修改、等同替换和改进，均应在本发明的权利范围之内。

Claims (6)

1.一种基于心电单极导联V1、V5R推衍导联V3R、V4R的方法，其特征在于，所述方法包括以下步骤：
获取心电单极导联V1的临床数据V1和心电单极导联V5R的临床数据V5R；
仅对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行基于临近单极导联线性相关的推衍处理，得到导联V3R的推衍数据V'3R和导联V4R的推衍数据V'4R；
对所述导联V3R的推衍数据V'3R和所述导联V4R的推衍数据V'4R分别进行幅值转化处理，得到所述导联V3R的合成数据V3R和所述导联V4R的合成数据V4R，以便利用所述导联V3R的合成数据V3R和所述导联V4R的合成数据V4R进行动态心电图分析；
其中，所述仅对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行基于临近单极导联线性相关的推衍处理，得到所述导联V3R的推衍数据V'3R和所述导联V4R的推衍数据V'4R包括：
利用用于推衍所述导联V3R的第一线性回归方程，对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行推衍处理，得到所述导联V3R的推衍数据V'3R，其中，所述第一线性回归方程为V'3R＝β3R1*V1+β3R2*V5R+β3R0，所述β3R1，β3R2，β3R0是所述第一线性回归方程的参数；
利用用于推衍所述导联V4R的第二线性回归方程，对所述心电单极导联V1的临床数据V1和所述心电单极导联V5R的临床数据V5R进行推衍处理，得到所述导联V4R的推衍数据V'4R，其中，所述第二线性回归方程为V'4R＝β4R1*V1+β4R2*V5R+β4R0，其中，β4R1，β4R2，β4R0是所述第二线性回归方程的参数；
其中，所述对所述导联V3R的推衍数据V'3R和所述导联V4R的推衍数据V'4R分别进行幅值转化处理，得到所述导联V3R的合成数据V3R和所述导联V4R的合成数据V4R包括：
获取推衍时间段内所述心电单极导联V1的临床数据，并根据所述推衍时间段内的所述心电单极导联V1的临床数据，确定所述推衍时间段内所述心电单极导联V1的临床数据平均值V1AMP；
获取所述推衍时间段内所述心电单极导联V5R的临床数据，并根据所述推衍时间段内的所述心电单极导联V5R的临床数据，确定所述推衍时间段内所述心电单极导联V5R的临床数据平均值V5RAMP；
利用所述推衍时间段内所述心电单极导联V1的临床数据平均值V1AMP和所述心电单极导联V5R的临床数据平均值V5RAMP，对所述推衍时间段内所述导联V3R的推衍数据V'3R进行幅值转化处理，得到所述导联V3R的合成数据V3R；
利用所述推衍时间段内所述心电单极导联V1的临床数据平均值V1AMP和所述心电单极导联V5R的临床数据平均值V5RAMP，对所述推衍时间段内所述导联V4R的推衍数据V'4R进行幅值转化处理，得到所述导联V4R的合成数据V4R。
2.根据权利要求1所述的方法，其特征在于，所述第一线性回归方程的参数β3R1、β3R2、β3R0的参数值通过以下步骤得到：
获取第一时间段内所述心电单极导联V1的已有临床数据、所述心电单极导联V3R的已有临床数据和所述心电单极导联V5R的已有临床数据；
利用所述第一时间段内的所述心电单极导联V1的已有临床数据、所述心电单极导联V3R的已有临床数据和所述心电单极导联V5R的已有临床数据，确定所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值；
对所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β3R1、β3R2、β3R0的最终的参数值。
3.根据权利要求2所述的方法，其特征在于，所述对所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β3R1、β3R2、β3R0的最终的参数值包括：
根据所述第一时间段内所述心电单极导联V1的已有临床数据，确定所述心电单极导联V1的已有临床数据平均值V1amp；
根据所述第一时间段内所述心电单极导联V5R的已有临床数据，确定所述心电单极导联V5R的已有临床数据平均值V5Ramp；
利用所述心电单极导联V1的已有临床数据平均值V1amp和所述心电单极导联V5R的已有临床数据平均值V5Ramp，对所述第一线性回归方程的参数β3R1、β3R2、β3R0的初始的参数值进行归一化处理，得到所述第一线性回归方程的参数β3R1、β3R2、β3R0的最终的参数值。
4.根据权利要求1所述的方法，其特征在于，所述第二线性回归方程的参数β4R1、β4R2、β4R0的参数值通过以下步骤得到：
获取第二时间段内所述心电单极导联V1的已有临床数据、所述心电单极导联V4R的已有临床数据和所述心电单极导联V5R的已有临床数据；
利用所述第二时间段内的所述心电单极导联V1的已有临床数据、所述心电单极导联V4R的已有临床数据和所述心电单极导联V5R的已有临床数据，确定所述第二线性回归方程的参数β4R1、β4R2、β4R0的初始的参数值；
对所述第二线性回归方程的参数β4R1、β4R2、β4R0的初始参的数值进行归一化处理，得到所述第二线性回归方程的参数β4R1、β4R2、β4R0的最终的参数值。
5.根据权利要求4所述的方法，其特征在于，所述对所述第二线性回归方程的参数β4R1、β4R2、β4R0的初始的参数值进行归一化处理，得到所述第二线性回归方程的参数β4R1、β4R2、β4R0的最终的参数值包括：
根据所述第二时间段内所述心电单极导联V1的已有临床数据，确定所述心电单极导联V1的已有临床数据平均值V1amp；
根据所述第二时间段内所述心电单极导联V5R的已有临床数据，确定所述心电单极导联V5R的已有临床数据平均值V5Ramp；
利用所述心电单极导联V1的已有临床数据平均值V1amp和所述心电单极导联V5R的已有临床数据平均值V5Ramp，对所述第二线性回归方程的参数β4R1、β4R2、β4R0的初始的参数值进行归一化处理，得到所述第二线性回归方程的参数β4R1、β4R2、β4R0的最终的参数值。
6.一种基于心电单极导联V1、V5R推衍导联V3R、V4R的设备，其特征在于，所述设备包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的程序，所述程序被所述处理器执行时实现如权利要求1-5任一项所述的基于心电单极导联V1、V5R推衍导联V3R、V4R的方法的步骤。
CN202110532230.2A 2021-05-17 2021-05-17 一种基于心电单极导联v1、v5r推衍导联v3r、v4r的方法及设备 Active CN113208599B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110532230.2A CN113208599B (zh) 2021-05-17 2021-05-17 一种基于心电单极导联v1、v5r推衍导联v3r、v4r的方法及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110532230.2A CN113208599B (zh) 2021-05-17 2021-05-17 一种基于心电单极导联v1、v5r推衍导联v3r、v4r的方法及设备

Publications (2)

Publication Number Publication Date
CN113208599A CN113208599A (zh) 2021-08-06
CN113208599B true CN113208599B (zh) 2024-02-27

Family

ID=77092163

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110532230.2A Active CN113208599B (zh) 2021-05-17 2021-05-17 一种基于心电单极导联v1、v5r推衍导联v3r、v4r的方法及设备

Country Status (1)

Country Link
CN (1) CN113208599B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4850370A (en) * 1987-07-22 1989-07-25 Dower Gordon E Method and apparatus for sensing and analyzing electrical activity of the human heart
CN1531902A (zh) * 2003-03-21 2004-09-29 比顿（北京）医用设备有限公司 以校正的正交心电图方式转换多种同步心电图导联方法
CN106073760A (zh) * 2016-04-29 2016-11-09 北京蓬阳丰业医疗设备有限公司 基于威尔逊导联的18导联动态心电图分析方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4587008B2 (ja) * 2000-07-24 2010-11-24 大名 魏 標準１２誘導心電図の構築方法および心電図検査装置
US6636761B2 (en) * 2000-12-29 2003-10-21 Ge Medical Systems Information Technologies, Inc. Method and apparatus for generating a twelve-lead ECG from fewer than ten electrodes

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4850370A (en) * 1987-07-22 1989-07-25 Dower Gordon E Method and apparatus for sensing and analyzing electrical activity of the human heart
CN1531902A (zh) * 2003-03-21 2004-09-29 比顿（北京）医用设备有限公司 以校正的正交心电图方式转换多种同步心电图导联方法
CN106073760A (zh) * 2016-04-29 2016-11-09 北京蓬阳丰业医疗设备有限公司 基于威尔逊导联的18导联动态心电图分析方法及系统

Also Published As

Publication number Publication date
CN113208599A (zh) 2021-08-06

Similar Documents

Publication Publication Date Title
US7142907B2 (en) Method and apparatus for algorithm fusion of high-resolution electrocardiograms
Feild et al. Improved EASI coefficients: their derivation, values, and performance
Meo et al. Spatial variability of the 12-lead surface ECG as a tool for noninvasive prediction of catheter ablation outcome in persistent atrial fibrillation
Castells et al. Estimation of atrial fibrillatory wave from single-lead atrial fibrillation electrocardiograms using principal component analysis concepts
CN112237431A (zh) 一种基于深度学习的心电参数计算方法
Llamedo et al. Cross-database evaluation of a multilead heartbeat classifier
EP2505135B1 (en) TWA measuring electrocardiograph, TWA measuring method, and TWA measurement system
Besrour et al. ECG beat classifier using support vector machine
CN113208599B (zh) 一种基于心电单极导联v1、v5r推衍导联v3r、v4r的方法及设备
CN113208600B (zh) 一种基于心电单极导联v6、v9推衍导联v7、v8的方法及设备
Plesinger et al. Multichannel QRS Morphology Clustering-Data Preprocessing for Ultra-High-Frequency ECG Analysis
CN113223702B (zh) 一种基于心电单极导联v2、v5推衍导联v3、v4的方法及设备
Sahoo et al. Adaptive thresholding based EMD for delineation of QRS complex in ECG signal analysis
Vozda et al. Novel Method for Deriving Vectorcardiographic Leads Based on Artificial Neural Networks
Nallikuzhy et al. Enhancement of the spatial resolution of ECG using multi-scale Linear Regression
Nasario-Junior et al. Beat-to-beat T-peak T-end interval duration variability assessed by RR-interval histogram analysis in health sedentary and Athlete
Schmidt et al. Evaluation of Ventricular Repolarization Variability in Patients With Nonischemic Dilated Cardiomyopathy From Vectorcardiography
Padhy et al. Synthesis of 12-lead ECG from a reduced lead set using singular value decomposition
Bear et al. Application of an inverse-forward approach to derive the 12-lead ECG from body surface potential maps
Kavya et al. Abnormality diagnosis in ECG signal using Daubechies wavelet
Mirmohamadsadeghi et al. Respiratory rate estimation from multi-lead ECGs using an adaptive frequency tracking algorithm
Loewe et al. Maximization of left atrial information through the optimization of ECG lead systems
Maan et al. A vector cardiographic based method to determine the culprit artery in acute coronary syndrome
Sueaseenak Denoising of electrocardiogram measurement system based on statistical signal processing
Guldenring et al. Influence of the Training Set Composition on the Estimation Performance of Linear ECG-Lead Transformations

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant