CN107440707A - 基于wilson心电图导联体系的18导/21导动态心电图方法 - Google Patents

基于wilson心电图导联体系的18导/21导动态心电图方法 Download PDF

Info

Publication number
CN107440707A
CN107440707A CN201710370553.XA CN201710370553A CN107440707A CN 107440707 A CN107440707 A CN 107440707A CN 201710370553 A CN201710370553 A CN 201710370553A CN 107440707 A CN107440707 A CN 107440707A
Authority
CN
China
Prior art keywords
lead
mtd
leads
mtr
mrow
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.)
Pending
Application number
CN201710370553.XA
Other languages
English (en)
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 MEDSUN MEDICAL EQUIPMENT Co Ltd
Original Assignee
BEIJING MEDSUN MEDICAL EQUIPMENT 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 MEDSUN MEDICAL EQUIPMENT Co Ltd filed Critical BEIJING MEDSUN MEDICAL EQUIPMENT Co Ltd
Priority to CN201710370553.XA priority Critical patent/CN107440707A/zh
Publication of CN107440707A publication Critical patent/CN107440707A/zh
Pending legal-status Critical Current

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/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/319—Circuits for simulating ECG signals
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Cardiology (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明提出了一种基于wilson心电图导联体系的18导/21导动态心电图方法，包括：基于Wilson导联体系计算胸前12导联至15导联同步心电图，包括：设置多个胸前的固定位置导联电极，根据胸前导联电极建立X轴、Y轴和Z轴；计算出心脏在胸廓内位置参考值为0电位；根据胸前的多个导联相对体表位置，以心脏0电位作为参考角度，推算出胸前各个导联轴角度和放大系数；基于Wilson导联体系推算肢体6导联体系；将胸前12导联同步心电图联合肢体6导联同步心电图，重建18同步导联至21同步导联心电图体系。本发明可以减少使用的导联电极数量，通过较少的电极数量即可实现18导联至21导联同步动态心电图，重建效果和精度高。

Description

基于wilson心电图导联体系的18导/21导动态心电图方法
技术领域
本发明涉及心电导联技术领域，特别涉及一种基于wilson心电图导联体系的18导/21导动态心电图方法。
背景技术
现有技术中的心电图系统中，对于胸导联采用如下设计：胸前Ⅴ1至V6导联采用一对一电极摆放方式，胸前导联采用Frank导联体系横面向量投影方式。即用X和Z轴以双极导联方式推导出胸前包含Ⅴ1、V2…V7R、V9等12个导联抽角度和设置一定放大系数来完成心电图同步描记。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此，本发明的目的在于提出一种基于wilson心电图导联体系的18导/21导动态心电图方法，可以大大减少了使用的导联电极数量，通过较少的电极数量即可实现18导联同步动态心电图，并且重建效果和精度高。
为了实现上述目的，本发明的实施例提供一种基于wilson心电图导联体系的18导/21导动态心电图方法，包括如下步骤:
步骤S1，基于Wilson导联体系计算胸前12导联至15导联同步心电图，包括如下步骤：
步骤S11，设置多个胸前的固定电极导联位置，根据所述胸前导联电极建立X轴、Y轴和Z轴，并根据所述XYZ三轴构建三维线性空间，所述构建三维线性空间，包括如下步骤：
以右肩窝电极位置及左肋下联电极位置联为Y轴，以V5电极位置和V5R电极位置联为X轴，以V2电极位置和V9电极位置联为Z轴，以上述XYZ轴构建三维线性空间；
现有技术是随机选取的三个互相垂直的导联作为XYZ轴，不同个体对此空间坐标标准的符合程度不同，而且实际测定中电极位置的放置存在差别，导致重建结果误差。以右肩窝电极位置及左肋下联电极位置联为Y轴，以V5电极位置和V5R电极位置联为X轴，以V2电极位置和V9电极位置联为Z轴的XYZ轴是固定的三维线性空间，而且V5、V5R、V2、V9是四个相对稳定且互相垂直的电极，所有个体均以此为空间坐标，重建结果更准确。
步骤S12，利用所述三维线性空间中的三轴位置，计算出心脏在胸廓内相对三维空间位置，采用数学方式计算得出心脏在胸廓内对各个角度的参照值均为0的位置，记为心脏0电位；
步骤S13，设置4个胸前的固定位置导联电极，分别对应胸前导联Ⅴ2、V9、V5R、Ⅴ5，以胸前导联V2至Ⅴ9为一轴，以Ⅴ5至V5R为另一轴，以心脏0电位为参考点，作为固有参照电极位置，推算出胸前各个导联轴角度和放大系数，采样上述两条相对垂直轴向电极形成两条胸前同步心电采样通道，将该两条胸前同步导联心电采样通道和两条肢导采样通道组成四条同步心电图通道；
步骤S14，利用上述四条同步心电图通道同步采集胸前导联和肢体导联的全信息心电信号，从而根据胸前的固定位置电极导联轴推算出的胸前虚拟电极导联轴，模拟出胸前12导联至15导联全信息，生成12至15胸前同步导联心电图电极；
其中，胸前的固定位置电极的数量在同步18导联心电图不少于12个，胸前的固定位置电极的数量在同步21导联心电图不少于15个；
步骤S2，基于Wilson导联体系推算肢体6导联体系；
步骤S3，将步骤S1中的胸前12导联至15导联同步心电图联合步骤S2中的肢体6导联同步心电图，重建18同步导联至21同步导联心电图体系。
进一步，在所述步骤S2中，所述基于Wilson导联体系推算肢体6导联体系，采用以下两种方式之一：
(1)分别在左右肩窝和左侧肋缘下放置电极，左右肩窝联联为I导联，右肩至左肋缘下联为II导联，推导出III、avR、avL、avF导联，构成标准肢体导联系统；
(2)采用步骤S1推算出的胸前各个导联轴角度和放大系数，进而推算出肢体导联信息。
进一步，所述导联方法使用电极数量位于4至10之间。
进一步，所述导联方法使用电极数量为9个，其中，胸前导联7个电极，包括胸前固定电极4个；肢体导联2个电极。
进一步，所述推算出胸前各个导联轴角度和放大系数，包括如下步骤：
根据胸前各导联相对体表位置，胸前其他导联分别以多个固定电极位置角度为参照角度向一定方向推算出各虚拟导联轴相应角度和放大系数，包括：
胸前导联V3在固定电极V2基础上向左推20至25度，为+60至65度，放大系数为1，其中，V2为+80度；
胸前导联V4R在固定电极V5R约+150至160度基础上向左推20至25度，为+140至130度，放大系数为1。
胸前导联V4在固定电极V5基础上向右推10至15度，为+40至45度，放大系数为1，其中，V5为+30度；
胸前导联V1在固定电极V5R约+150至160度基础上向左推50至55度，为+100至110度，放大系数为1。
胸前导联V6在固定电极V5基础上向左推10至15度，为+15至20度，放大系数为1，其中，V5为+30度；
胸前导联V3R在固定电极V5R约+150至160度基础上向左推40至45度，为+110至120度，放大系数为1。
进一步，在所述步骤S13之后，还包括如下步骤:利用病例数据库中的多个病种的实际病例对得到的各个导联轴角度和放大系数进行调整，直至符合临床认可标准。
进一步，在所述步骤S3中，所述重建18同步导联心电图体系，包括如下步骤：
选取m个导联作为初始的用于重建的导联，这个重建导联组S可表示为：
S＝{L(1),...,L(i),...,L(m)}
18导联心电图ECG系统表示为：
E＝{Ι,П,Ш,aVR,aVL,aVF,V1,V2,V3,V4,V5,V6,V7,V8,V9,V3R,V4R,V5R}
当有n组样本(S1,E1),...,(Sn,En),(n＞m)时，E与S之间的线性模型可表示为：
E＝M·β+ε
其中，E是样本Ei,(i＝1,...,n)，β是转换矩阵，ε是一个误差矩阵。而M定义为：
分析标准18导联采样点的数据，用样本导联的采样点数据去拟合剩下的导联数据，拟合过程采用最小二乘法来评估，确定估计矩阵b，
则重建后的18导联矩阵可以表示为：
根据本发明实施例的基于wilson心电图导联体系的18导/21导动态心电图方法，采用标准Wilson方式，并用硬件和软件相合实现的一种全数字方式全新导联体系。本发明的心电信号釆集全部运用Wilson导联体系为基本原理，胸部导联采用Wilson体系单极导联方式为基础设计，运用虚拟(模拟)和半虚拟(半模拟)方式，在极少量固有参照电极基础上运用数学计算推导方式，以胸部少量体表电极导出多条同步导联心电图。即，以少量固定位置电极以WiIson方式推导出12至15个包含V1至V9及V3R至V5R同步胸前心电图导联，再联合肢体导联6导联，组成18至21同步导联心电图体系。本发明可以大大减少了使用的导联电极数量，通过较少的电极数量即可实现18导联至21导联同步动态心电图，并且重建效果和精度高。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为根据本发明实施例的基于wilson心电图导联体系的18导/21导动态心电图方法的流程图；
图2a和图2b为根据本发明实施例的以心脏0电位为参考点的导联轴向角度的示意图。
具体实施方式
下面详细描述本发明的实施例，实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
本发明提出一种基于心电图的Wilson导联体系的虚拟或半虚拟单极胸导联方法，通过本发明可以重建18同步导联至21同步导联心电图体系。
在介绍本发明之前，首先对同步18导联动态心电图系统DCG(DynamicElectroCardioGram，动态心电图)进行说明。
18导联动态DCG可以有效地覆盖全心脏各个部位，从而可以作到同步采集心脏各部位长程心电图的改变，由于覆盖面全，因此称为全信息心电图。并且由于是同步长程监控，这对临床来讲，可以向一线医务人员客观地、全面地监控病人心肌缺血情况和病情严重程度提供了有效手段。
18导联DCG不同于传统的12导联DCG，除了12导DCG监控的前壁、下壁、侧壁、间壁心电情况，同时可同步监控左心室正后壁和右心室心肌缺血性疾病的动态变化情况，尤其是可及时、有效观察到这两个部位心肌缺血、梗死等病情发展和变化的全过程。
具体来说，这两个部位心肌缺血和梗死在临床上有发病时间短、常和其它部位并发的特点，这个在12导联无法检测到的。利用18导联DCG可以向心内科一线临床提供了更全面、更及时监控心肌缺血性疾病的有效手段。因此，可以理解为12导联DCG取代了3导联DCG，实现了DCG仅仅是用于心律失常监控的桎梏，而18导联DCG则在12导联DCG基础上更全面地监控心脏各部位、各种类心肌缺血性疾病的全面动态情况，是心电图应用技术的一个重要发展。
如图1所示，本发明实施例的基于心电图Wilson的导联体系的虚拟或半虚拟单极胸导联方法，包括如下步骤：
步骤S1，基于Wilson导联体系计算胸前12至15导联同步心电图，包括如下步骤：
步骤S11，设置多个胸前的固定电极导联位置，根据胸前导联电极建立X轴、Y轴和Z轴，并根据XYZ三轴构建三维线性空间。
需要说明的是，本发明提供的导联方法使用电极数量位于4至10之间。
优选的，导联方法使用电极数量为8个，其中，胸前导联4个电极，肢体导联4个电极。
在本发明的一个实施例中，设置4个胸前的固定位置导联电极，分别对应胸前导联Ⅴ2、V9、V5R、Ⅴ5。
步骤S12，利用三维线性空间中的三轴位置，相对三维空间位置，采用数学方式计算得出心脏在胸廓内对各个角度的参照值均为0的位置，记为心脏0电位。
下面对数学计算的过程进行说明。
Wilson中心电端(WCT)由RA、LA、LL电极电位的平均值，即WCT＝(RA+LA+LL)/3，Kirchhoff定律不要求WCT电位在心动周期中是0或保持恒定。WCT与RA、LA、LL之间的电位差被分别定义为新的额面肢体导联VR、VL和VF。Wilson称这些电极对为“单极”导联。
具体地，构建三维线性空间，包括如下步骤：
以右肩窝电极位置及左肋下联电极位置联为Y轴，以V5电极位置和V5R电极位置联为X轴，以V2电极位置和V9电极位置联为Z轴，以上述XYZ轴构建三维线性空间。利用三轴组成空间位置，用数学方式实现计算设定心脏在胸腔内相对位置(第四至五肋间水平)为0。
在本发明的一个实施例中，在步骤S12和S13之间，还包括如下步骤：以胸前导联V2至Ⅴ9为一轴，以Ⅴ5至V5R为另一轴，采样上述两条相对垂直轴向电极形成两条胸前同步心电采样通道，将该两条胸前同步导联心电采样通道和两条肢导采样通道组成四条同步心电图通道，利用上述四条同步心电图通道同步采集胸前导联和肢体导联的全信息心电信号。
步骤S13，利用多个胸前的固定位置导联电极，以心脏0电位为参考点，作为固有参照电极位置，推算出胸前各个导联轴角度和放大系数，从而根据胸前的固定位置电极导联轴推算出的胸前虚拟电极导联轴，模拟出胸前12导联至15导联全信息，生成12至15胸前同步导联心电图电极。其中，12导联包括：12个V1至V9及V3R至V5R同步胸前心电图导联，即V1、V2、V3、V4、V5、V5、V6、V7、V8、V9、V3R、V4R、V5R。
需要说明的是，胸前的固定位置电极的数量在同步18导联心电图均不少于12个，胸前的固定位置电极的数量在同步21导联心电图均不少于15个。
例如，真实存在的胸前固定电极数量为4个，由该4个真实电极推算出8个虚拟导联电极，共同组成胸前12导联电极；
真实存在的胸前固定电极数量为4个，由该4个真实电极推算出11个虚拟导联电极，共同组成胸前15导联电极。
在本发明的一个实施例中，推算出胸前各个导联轴角度和放大系数，包括如下步骤：
参考图2a和图2b，根据胸前各导联相对体表位置，胸前其他导联分别以多个固定电极位置(Ⅴ2、V9、V5R、Ⅴ5)角度为参照角度向一定方向推算出各导联轴相应角度和放大系数，包括：
胸前导联V3在固定电极V2基础上向左推20至25度，为+60至65度，放大系数为1，其中，V2为+80度；
胸前导联V4R在固定电极V5R约+150至160度基础上向左推20至25度，为+140至130度，放大系数为1。通过采用上述方式可以推算出胸前12个同步心电图导联轴并形成同步胸导心电图。
需要说明的是，上述角度不是算出来的，而是根据wilson体系标准的电极安放位置转换成角度数据便于运算。
此外，本发明还包括如下步骤：利用病例数据库中的足够数量的多个病种的实际病例对得到的各个导联轴角度和放大系数进行调整，直至符合临床认可标准。
步骤S2，基于Wilson导联体系推算肢体6导联体系。
在本发明的一个实施例中，基于Wilson导联体系推算肢体6导联体系，采用以下两种方式之一：
(1)分别在左右肩窝和左侧肋缘下放置电极，左右肩窝联联为I导联，右肩至左肋缘下联为II导联，推导出III、avR、avL、avF导联，构成标准肢体导联系统。
具体地，在步骤S2中，基于Wilson导联体系，分别在左右肩窝(锁骨下缘部位)和左侧肋缘下(腋前线)放置电极，左右肩窝联联为I导联，右肩至左肋缘下联为II导联，用计算方式推导出III、avR、avL、avF导联，构成标准肢体导联系统。
(2)采用步骤S1推算出的胸前各个导联轴角度和放大系数，进而推算出肢体导联信息。
具体地，通过从中心电端电位中除去单个探查电极的电位，Goldberger“加压单极”肢体导联，之所以如此称谓是因为从数学角度上说，记录的波幅与利用Wilson中心电端记录的相比增大50％。加压肢体导联的Goldberger中心电端通过以下方法计算获得，即aVR导联为(LA+LL)/2，aVL导联为(RA+LL)/2，aVF为(RA+LA)/2。因此，aVL导联代表了左上肢与改良后的Goldberger中心电端的电位差，其计算方法为LA-(RA+LL)/2，其可简化为(Ⅰ导联-Ⅲ导联)/2。同样的，aVR导联为RA-(LA+LL)/2，可简化为-(Ⅰ导联+Ⅱ导联)/2；aVF导联为LL-(LA+RA)/2，可简化为(Ⅱ导联+Ⅲ导联)/2。这些衍生的导联提供了额面内的新向量视角。
步骤S3，将步骤S1中的胸前12导联至15导联同步心电图联合步骤S2中的肢体6导联同步心电图，重建18同步导联至21同步导联心电图体系。
首先，对导联之间的关系进行说明：
标准12导联ECG信号由肢体导联(Ⅰ、Ⅱ、Ⅲ、aVR、aVL和aVF)和胸导联(V1、V2、V3、V4、V5和V6)组成，18导联是在标准12导联的基础上又向右胸增加V3R、V4R和V5R和向左胸增加V7、V8和V9导联。其中，各导联之间的关系如下：
Ⅰ-Ⅱ-Ⅲ＝0；
aVR＝-(Ⅰ+Ⅱ)/2；
aVL＝Ⅰ-Ⅱ/2；
aVF＝Ⅱ-Ⅰ/2。
下面对心电图的重建过程进行说明。
根据重建算法的不同，导联重建方法主要可以分为线性重建法和非线性重建方法。下面以线性导联重建方法为例进行说明。绝大多数的导联重建系统都用线性变换的方法来进行ECG导联的重建。
假设选取m个导联作为初始的用于重建的导联，这个重建导联组S可表示为：
S＝{L(1),...,L(i),...,L(m)}
18导联ECG系统表示为：
E＝{Ι,П,Ш,aVR,aVL,aVF,V1,V2,V3,V4,V5,V6,V7,V8,V9,V3R,V4R,V5R}
当有n组样本(S1,E1),...,(Sn,En),(n＞m)时，E与S之间的线性模型可表示为：
E＝M·β+ε
其中，E是样本Ei,(i＝1,...,n)，β是转换矩阵，ε是一个误差矩阵。而M定义为：
因此，只需将转换矩阵β的估计矩阵b算出来就能将18导联重建出来。
本发明利用统计的方法来计算估计矩阵b，通过分析标准18导联采样点的数据，用样本导联的采样点数据去拟合剩下的导联数据，拟合过程采用最小二乘法来评估，最终确定估计矩阵b，实际上，在矩阵运算中，b＝(MTM)-1MTE即为线性回归运算后的结果。
在转换矩阵β的估计矩阵b确定后，则重建后的18导联矩阵可以表示为：
由此实现了导联的重建，而理论上统计的方法可以最少用2个导联来重建出标准的18导联，但随着用于重建的原始导联数量的减少，重建效果也会越差，因此需要选取合适的导联数量及导联种类来进行重建。
下面以9导联为例，对重建18导联心电图的重建算法进行说明：
心向量图与体表的心电数据之间存在着投影的关系，按照导联向量的概念，每个导联有一个向量L与之相对应，导联向量L与心向量H的点积，就是该导联的电压值V，即V＝L·H。
设心向量为(l1,l2,l3,l4,l5,l6,l7,l8,l9)T时，记录的九导联数据为(a1,a2,a3,a4,a5,a6,a7,a8,a9)T,
则(a1,a2,a3,a4,a5,a6,a7,a8,a9)T＝A·(l1,l2,l3,l4,l5,l6,l7,l8,l9)T，其中A为投影矩阵。
由于九导联不在同一平面，线性无关，矩阵A可逆，
则(l1,l2,l3,l4,l5,l6,l7,l8,l9)T＝A-1·(a1,a2,a3,a4,a5,a6,a7,a8,a9)T可由此求出心电向量的分量，即在L1至L9九个导联轴上的投影。对这九个导联进行正交化变换，计算出导联向量逆投影的值，即得出相应心向量9个导联。
导联向量是各导联在九个轴上的投影系数：
L＝(a1,a2,a3,a4,a5,a6,a7,a8,a9)·(l1,l2,l3,l4,l5,l6,l7,l8,l9)T，因此
L＝(a1,a2,a3,a4,a5,a6,a7,a8,a9)·A-1·(l1,l2,l3,l4,l5,l6,l7,l8,l9)T。
下面令
(a1,a2,a3,a4,a5,a6,a7,a8,a9)·A-1＝(a1',a2',a3',a4',a5',a6',a7',a8',a9')，
则(a1',a2',a3',a4',a5',a6',a7',a8',a9')为各导联在L1～L9这9个原始导联上的投影系数，进而可以重建18导联的心电图数据。
本发明中，从线性代数的角度看，心电向量是三维线性空间中的向量，XYZ轴是这个线性空间中的一个基，导联轴也是这个线性空间中的向量。当九个导联轴的向量线性无关时，它们就构成了线性空间的一个基，从而可以求出各导联在这九个原始导联上的投影系数并作为重建系数，进而可以计算出十八导联的心电图数据。
本申请采用改良的Wilson威尔逊导联体系，采用较少的体表电极少、具有图形还原好、干扰少(运用新的导联体系，采用新的电极贴法，可有效在一定程度上减少干扰的产生)的鲜明持点。本发明仅运用仅用4条心电采样通道，相比真18导为16个体表电极、14条心电通道，并采用全Wilson方式，而非正交方式的方案，本发明体现了心电信号采样同步性好、通道少、电极片少、效率高、图形还原率高等特性。
本发明在面向医院和病人角度，均具有一定有益效果：
1)对医院来讲，可减少耗材，体现了节约精神，少消耗获得大信息量；
2)对病人来讲，由于电板少配带明显较为舒适，易被病人接受；
3)对于发展看，由于通道数明显减少，可有效控制信息量，充分扩大信息利用空间，极有利于远程同步心电信号传输工作的延伸。
根据本发明实施例的基于wilson心电图导联体系的18导/21导动态心电图方法，采用标准Wilson方式，并用硬件和软件相合实现的一种全数字方式全新导联体系。本发明的心电信号釆集全部运用Wilson导联体系为基本原理，胸部导联采用Wilson体系单极导联方式为基础设计，运用虚拟(模拟)和半虚拟(半模拟)方式，在极少量固有参照电极基础上运用数学计算推导方式，以胸部少量体表电极导出多条同步导联心电图。即，以少量固定位置电极以WiIson方式推导出12个至15个包含V1至V9及V3R至V5R同步胸前心电图导联，再联合肢体导联6导联，组成18至21同步导联心电图体系。本发明可以大大减少了使用的导联电极数量，通过较少的电极数量即可实现18导联至21导联同步动态心电图，并且重建效果和精度高。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。

Claims (7)

1.一种基于wilson心电图导联体系的18导/21导动态心电图方法，其特征在于，包括如下步骤:
步骤S1，基于Wilson导联体系计算胸前12导联至15导联同步心电图，包括如下步骤：
步骤S11，设置多个胸前的固定电极导联位置，根据所述胸前导联电极建立X轴、Y轴和Z轴，并根据所述XYZ三轴构建三维线性空间，所述构建三维线性空间，包括如下步骤：
以右肩窝电极位置及左肋下联电极位置联为Y轴，以V5电极位置和V5R电极位置联为X轴，以V2电极位置和V9电极位置联为Z轴，以上述XYZ轴构建三维线性空间；
步骤S12，利用所述三维线性空间中的三轴位置，计算出心脏在胸廓内相对三维空间位置，采用数学方式计算得出心脏在胸廓内对各个角度的参照值均为0的位置，记为心脏0电位；
步骤S13，设置4个胸前的固定位置导联电极，分别对应胸前导联Ⅴ2、V9、V5R、Ⅴ5，以胸前导联V2至Ⅴ9为一轴，以Ⅴ5至V5R为另一轴，以心脏0电位为参考点，作为固有参照电极位置，推算出胸前各个导联轴角度和放大系数，采样上述两条相对垂直轴向电极形成两条胸前同步心电采样通道，将该两条胸前同步导联心电采样通道和两条肢导采样通道组成四条同步心电图通道；
步骤S14，利用上述四条同步心电图通道同步采集胸前导联和肢体导联的全信息心电信号，从而根据胸前的固定位置电极导联轴推算出的胸前虚拟电极导联轴，模拟出胸前12导联至15导联全信息，生成12至15胸前同步导联心电图电极；
其中，胸前的固定位置电极的数量在同步18导联心电图不少于12个，胸前的固定位置电极的数量在同步21导联心电图不少于15个；
步骤S2，基于Wilson导联体系推算肢体6导联体系；
步骤S3，将步骤S1中的胸前12导联至15导联同步心电图联合步骤S2中的肢体6导联同步心电图，重建18同步导联至21同步导联心电图体系。
2.如权利要求1所述的基于wilson心电图导联体系的18导/21导动态心电图方法，其特征在于，在所述步骤S2中，所述基于Wilson导联体系推算肢体6导联体系，采用以下两种方式之一：
(1)分别在左右肩窝和左侧肋缘下放置电极，左右肩窝联联为I导联，右肩至左肋缘下联为II导联，推导出III、avR、avL、avF导联，构成标准肢体导联系统；
(2)采用步骤S1推算出的胸前各个导联轴角度和放大系数，进而推算出肢体导联信息。
3.如权利要求1所述的基于wilson心电图导联体系的18导/21导动态心电图方法，其特征在于，所述导联方法使用电极数量位于4至10之间。
4.如权利要求4所述的基于心电图Wilson联体系的虚拟或半虚拟单极胸导联方法，其特征在于，所述导联方法使用电极数量为9个，其中，胸前导联7个电极，包括胸前固定电极4个；肢体导联2个电极。
5.如权利要求1所述的基于wilson心电图导联体系的18导/21导动态心电图方法，其特征在于，所述推算出胸前各个导联轴角度和放大系数，包括如下步骤：
根据胸前各导联相对体表位置，胸前其他导联分别以多个固定电极位置角度为参照角度向一定方向推算出各虚拟导联轴相应角度和放大系数，包括：
胸前导联V3在固定电极V2基础上向左推20至25度，为+60至65度，放大系数为1，其中，V2为+80度；
胸前导联V4R在固定电极V5R约+150至160度基础上向左推20至25度，为+140至130度，放大系数为1。
6.如权利要求1所述的基于wilson心电图导联体系的18导/21导动态心电图方法，其特征在于，在所述步骤S13之后，还包括如下步骤:利用病例数据库中的多个病种的实际病例对得到的各个导联轴角度和放大系数进行调整，直至符合临床认可标准。
7.如权利要求1所述的基于wilson心电图导联体系的18导/21导动态心电图方法，其特征在于，在所述步骤S3中，所述重建18同步导联心电图体系，包括如下步骤：
选取m个导联作为初始的用于重建的导联，这个重建导联组S可表示为：
S＝{L(1),...,L(i),...,L(m)}
18导联心电图ECG系统表示为：
E＝{Ι,П,Ш,aVR,aVL,aVF,V1,V2,V3,V4,V5,V6,V7,V8,V9,V3R,V4R,V5R}
当有n组样本(S1,E1),...,(Sn,En),(n＞m)时，E与S之间的线性模型可表示为：
E＝M·β+ε
其中，E是样本Ei,(i＝1,...,n)，β是转换矩阵，ε是一个误差矩阵。而M定义为：
<mrow> <mi>M</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mi>L</mi> <msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <mrow> <mi>L</mi> <msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> <mtd> <mtable> <mtr> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> </mtr> </mtable> </mtd> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mi>L</mi> <msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>n</mi> </msub> </mrow> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <mrow> <mi>L</mi> <msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mi>n</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
分析标准18导联采样点的数据，用样本导联的采样点数据去拟合剩下的导联数据，拟合过程采用最小二乘法来评估，确定估计矩阵b，
则重建后的18导联矩阵可以表示为：
CN201710370553.XA 2017-05-23 2017-05-23 基于wilson心电图导联体系的18导/21导动态心电图方法 Pending CN107440707A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710370553.XA CN107440707A (zh) 2017-05-23 2017-05-23 基于wilson心电图导联体系的18导/21导动态心电图方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710370553.XA CN107440707A (zh) 2017-05-23 2017-05-23 基于wilson心电图导联体系的18导/21导动态心电图方法

Publications (1)

Publication Number Publication Date
CN107440707A true CN107440707A (zh) 2017-12-08

Family

ID=60486955

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710370553.XA Pending CN107440707A (zh) 2017-05-23 2017-05-23 基于wilson心电图导联体系的18导/21导动态心电图方法

Country Status (1)

Country Link
CN (1) CN107440707A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111543979A (zh) * 2020-05-13 2020-08-18 许祥林 一种常规导联输出心电向量图的方法
WO2020259725A1 (en) * 2019-06-27 2020-12-30 Vladimir Kranz A pocket ecg for heart attack diagnostics
CN112353398A (zh) * 2020-12-01 2021-02-12 山东众阳健康科技集团有限公司 一种基于无监督学习的心电向量重建方法
CN113223702A (zh) * 2021-05-17 2021-08-06 北京蓬阳丰业科技有限公司 一种基于心电单极导联v2、v5推衍导联v3、v4的方法及设备
CN117100283A (zh) * 2023-07-26 2023-11-24 兴化市人民医院 一种十二导联心电图的生成方法和装置
CN117442211A (zh) * 2023-12-26 2024-01-26 纳龙健康科技股份有限公司 一种多导联心电图同步采集方法、终端设备及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020087088A1 (en) * 2000-12-29 2002-07-04 Brodnick Donald Eugene Method and apparatus for generating a twelve-lead ECG from fewer than ten electrodes
CN1531902A (zh) * 2003-03-21 2004-09-29 比顿（北京）医用设备有限公司 以校正的正交心电图方式转换多种同步心电图导联方法
CN101011242A (zh) * 2007-02-06 2007-08-08 赵峰 一种向量心电图仪及实现方法
US20110028821A1 (en) * 2009-07-31 2011-02-03 Newcardio, Inc. Electrocardiographic Monitoring System and Method Using Orthogonal Electrode Pattern
CN102934994A (zh) * 2012-10-26 2013-02-20 杭州师范大学 心电轴及钟向转位测量盘
CN106073760A (zh) * 2016-04-29 2016-11-09 北京蓬阳丰业医疗设备有限公司 基于威尔逊导联的18导联动态心电图分析方法及系统
CN106073753A (zh) * 2016-04-29 2016-11-09 北京蓬阳丰业医疗设备有限公司 18导联全息动态及静态心电图分析方法及系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020087088A1 (en) * 2000-12-29 2002-07-04 Brodnick Donald Eugene Method and apparatus for generating a twelve-lead ECG from fewer than ten electrodes
CN1531902A (zh) * 2003-03-21 2004-09-29 比顿（北京）医用设备有限公司 以校正的正交心电图方式转换多种同步心电图导联方法
CN101011242A (zh) * 2007-02-06 2007-08-08 赵峰 一种向量心电图仪及实现方法
US20110028821A1 (en) * 2009-07-31 2011-02-03 Newcardio, Inc. Electrocardiographic Monitoring System and Method Using Orthogonal Electrode Pattern
CN102934994A (zh) * 2012-10-26 2013-02-20 杭州师范大学 心电轴及钟向转位测量盘
CN106073760A (zh) * 2016-04-29 2016-11-09 北京蓬阳丰业医疗设备有限公司 基于威尔逊导联的18导联动态心电图分析方法及系统
CN106073753A (zh) * 2016-04-29 2016-11-09 北京蓬阳丰业医疗设备有限公司 18导联全息动态及静态心电图分析方法及系统

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020259725A1 (en) * 2019-06-27 2020-12-30 Vladimir Kranz A pocket ecg for heart attack diagnostics
CN111543979A (zh) * 2020-05-13 2020-08-18 许祥林 一种常规导联输出心电向量图的方法
CN111543979B (zh) * 2020-05-13 2024-02-06 许祥林 一种常规导联输出心电向量图的方法
CN112353398A (zh) * 2020-12-01 2021-02-12 山东众阳健康科技集团有限公司 一种基于无监督学习的心电向量重建方法
CN113223702A (zh) * 2021-05-17 2021-08-06 北京蓬阳丰业科技有限公司 一种基于心电单极导联v2、v5推衍导联v3、v4的方法及设备
CN113223702B (zh) * 2021-05-17 2022-05-03 北京蓬阳丰业科技有限公司 一种基于心电单极导联v2、v5推衍导联v3、v4的方法及设备
CN117100283A (zh) * 2023-07-26 2023-11-24 兴化市人民医院 一种十二导联心电图的生成方法和装置
CN117100283B (zh) * 2023-07-26 2024-05-03 兴化市人民医院 一种十二导联心电图的生成方法和装置
CN117442211A (zh) * 2023-12-26 2024-01-26 纳龙健康科技股份有限公司 一种多导联心电图同步采集方法、终端设备及存储介质
CN117442211B (zh) * 2023-12-26 2024-03-15 纳龙健康科技股份有限公司 一种多导联心电图同步采集方法、终端设备及存储介质

Similar Documents

Publication Publication Date Title
CN107440707A (zh) 基于wilson心电图导联体系的18导/21导动态心电图方法
US7136693B2 (en) Method and apparatus for generating a twelve-lead ECG from fewer than ten electrodes
US5711304A (en) Signal processing apparatus and method for adding additional chest leads to the 12-lead electrocardiogram without additional electrodes
Feild et al. Improved EASI coefficients: their derivation, values, and performance
US6052615A (en) Method and apparatus for sensing and analyzing electrical activity of the human heart using a four electrode arrangement
CN106073753A (zh) 18导联全息动态及静态心电图分析方法及系统
Ihara et al. Atrial repolarization as observable during the PQ interval
CN103829941A (zh) 一种多维心电信号成像系统及方法
CN1377244A (zh) 采用胸导联重建Frank导联的方法和装置
US6496720B1 (en) Process for sensing and analyzing electrical activity of the human heart utilizing one lead system with an egg monitor designed for use with another lead system
CN107440710A (zh) 基于wilson心电图导联体系的4电极至7电极同步采集方法
CN203138474U (zh) 一种18导联心电工作站
CN107440711A (zh) 基于wilson心电图导联体系的8电极和10电极同步采集方法
US10085662B2 (en) Device and method for digital-to-analog transformations and reconstructions of multi-channel electrocardiograms
CN107440706A (zh) 基于wilson心电图导联体系中的9电极同步采集方法
Hoekema et al. Lead system transformation of body surface map data
Mulyadi et al. Improving accuracy of derived 12-lead electrocardiography by waveform segmentation
CN111543979B (zh) 一种常规导联输出心电向量图的方法
Bear et al. Application of an inverse-forward approach to derive the 12-lead ECG from body surface potential maps
CN104027106A (zh) 一种心电断层扫描成像系统及方法
Cluitmans et al. Spatiotemporal activation time estimation improves noninvasive localization of cardiac electrical activity
Sindreu et al. Standard 12-lead ECG synthesis from homecare wearable measures
CN112353398B (zh) 一种基于无监督学习的心电向量重建方法
Taccardi Present and future of body surface electrocardiographic mapping
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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 100010 room 206, South Main Building, 9-17 Zhongli, yongdingmennei East Street, Dongcheng District, Beijing

Applicant after: Beijing Pengyang Fengye Technology Co.,Ltd.

Address before: 100050 South Main Building 206, Zhongli 9-17, mennei East Street, Yongding, Dongcheng District, Beijing

Applicant before: BEIJING PENGYANG FENGYE MEDICAL EQUIPMENT Co.,Ltd.

RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20171208