CN102727193B - Twa测量心电图仪、twa测量方法和twa测量系统 - Google Patents

Twa测量心电图仪、twa测量方法和twa测量系统 Download PDF

Info

Publication number
CN102727193B
CN102727193B CN201210090420.4A CN201210090420A CN102727193B CN 102727193 B CN102727193 B CN 102727193B CN 201210090420 A CN201210090420 A CN 201210090420A CN 102727193 B CN102727193 B CN 102727193B
Authority
CN
China
Prior art keywords
twa
coefficient
electrocardiograph
individual
patient
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
CN201210090420.4A
Other languages
English (en)
Other versions
CN102727193A (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.)
Nihon Kohden Corp
Original Assignee
Nihon Kohden Corp
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 Nihon Kohden Corp filed Critical Nihon Kohden Corp
Publication of CN102727193A publication Critical patent/CN102727193A/zh
Application granted granted Critical
Publication of CN102727193B publication Critical patent/CN102727193B/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/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/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/355—Detecting T-waves
    • 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/327—Generation of artificial ECG signals based on measured signals, e.g. to compensate for missing leads
    • 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/339—Displays specially adapted therefor
    • A61B5/341—Vectorcardiography [VCG]
    • 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/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • 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
    • A61B5/7235—Details of waveform analysis
    • A61B5/7253—Details of waveform analysis characterised by using transforms
    • A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms

Landscapes

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

Abstract

一种TWA测量心电图仪包括：变换系数存储部件，其构造成存储用于获得弗兰克矢量心电图的变换系数；心电图仪控制部件，其构造成从适于连接于受测者的测量电极的心电图信号产生标量心电图，并且构造成通过利用所述变换系数从所述标量心电图获得所述弗兰克矢量心电图；和TWA测量部件，其构造成基于获得的弗兰克矢量心电图来测量TWA的存在。

Description

TWA测量心电图仪、TWA测量方法和TWA测量系统
发明背景
本发明涉及一种能够测量TWA(T波交替)的存在的TWA测量心电图仪、TWA测量方法和TWA测量系统，尤其是涉及能够从标量心电图得到弗兰克矢量心电图，并且从弗兰克矢量心电图测量TWA的存在的TWA测量心电图仪、TWA测量方法和TWA测量系统。
TWA出现在诸如QT延长综合症、变异型心绞痛、急性心肌缺血、电解质异常、阵发性心动过速、心动过缓或心包膜液体积聚的病症发作。TWA是一种在心电图中出现的T波的幅度和极性交替变化的现象，是有效地预示突发性心脏死亡的指示。TWA不是总能由裸眼观测到的现象，因此其在临床上的应用受到限制。
从20世纪80年代以来，从而使细微的TWA(微伏特TWA:MTWA)能够被计算机检测到的技术得到发展。
当前提出的用于测量TWA的技术的实例是基于分别在美国专利6,668,189和美国专利5,935,082中公开的通用电气(GE)公司的MMA(ModifiedMovingAverage改进的移动平均)方法和CambridgeHeart(CH)有限公司的周期图的测量技术。
GE公司的测量技术指导了一种在时间域上分析时间波形的方法，据说对噪声有阻抗。然而，该技术作为测量技术并没有很长的时间，并且需要观察其临床反应。
相反，在频域上使用的CH有限公司的测量技术从20世纪80年代就已经使用，并且其临床效果已经被证实。因此，当前认为基于CH有限公司的周期图的测量技术在临床上比基于GE公司的MMA方法更有用。
关于基于CH有限公司的周期图的测量技术，在其公开之后，用于执行新处理的诸如测量电极的技术的各种技术被增加，并且现在增加的最新技术已经被使用。
然而，基于CH有限公司的周期图的测量技术存在几个问题。例如，必须使用消除噪声的专用的测量电极，并且必须获得通常不一定使用的弗兰克矢量心电图。
专用测量电极的使用或弗兰克矢量心电图的获得给测量人员增加负担，因为测量电极的设置等与通常要获得的12导联心电图的情况不同。而且，在弗兰克导联方法中，电极还必须被放置在背部。因此，也增加了受测者的负担。当普通测量人员使用不习惯的弗兰克导联方法时，与普遍使用的标准感应相比，看来可能产生由电极的连接位置的偏移导致的较大的测量误差。
发明内容
本发明的目的在于提供一种能够从诸如12导联心电图的标量心电图获得弗兰克矢量心电图，并且从弗兰克矢量心电图测量TWA的存在的TWA测量心电图仪、TWA测量方法和TWA测量系统，从而使测量人员能够以与测量人员获得标量心电图的情况相同的方式来测量TWA的存在。
本发明的目的还在于能够从通常使用的诸如12导联心电图仪、训练用心电图仪或霍尔特(Holter)心电图仪的心电图仪，通过利用基于临床有效性已经被证实的矢量心电图的TWA的分析方法来测量TWA的存在。
为了实现上述目的，根据本发明，提出一种TWA测量心电图仪，其包括：变换系数存储部件，其构造成存储用于获得弗兰克矢量心电图的变换系数；心电图仪控制部件，其构造成从适用于连接到受测者的测量电极的心电图信号产生标量心电图，并且构造成利用所述变换系数从所述标量心电图推导弗兰克矢量心电图；和TWA测量部件，其构造成基于所推导的弗兰克矢量心电图测量TWA的存在。
该变换系数存储部件可以包括：个人系数数据库，其构造成存储从受测者得到的，并且所述受测者专用的个人系数作为变换系数；和群系数数据库，其构造成存储作为变换系数的群系数，该群系数为从具有统计效果的人口的非特定数量的人得到的多个变换系数的平均值。
当个人系数数据库存储个人系数时，心电图仪控制部件可以通过利用个人系数作为变换系数来推导弗兰克矢量心电图，和当个人系数数据库没有存储个人系数时，心电图仪控制部件可以通过利用存储在群系数数据库中的群系数作为变换系数来推导弗兰克矢量心电图。
当个人系数数据库没有存储个人系数时，心电图仪控制部件使显示部件显示用于提示获得个人系数的消息。
心电图仪控制部件可以从测量电极的心电图信号计算个人系数，并且将所计算的个人系数存储在个人系数数据库中。
TWA测量部件可以基于从弗兰克矢量心电图获得的矢量大小的波形来测量TWA的存在。
TWA测量部件可以从矢量大小的波形来选择有助于测量TWA的存在的波形，并且从所选择的波形测量TWA的存在。
在选择的有助于测量TWA存在的波形中，TWA测量部件可以在奇数或偶数次心跳中检测对测量TWA的存在具有不利影响的异常值，并且校正包括所检测的异常值的心跳。
根据本发明，还提出一种TWA测量方法，其包括：将测量电极连接于受测者；从测量电极的心电图信号产生标量心电图；推导用于推导弗兰克矢量心电图的变换系数；通过利用所述变换系数从所述标量心电图推导弗兰克矢量心电图；和从所推导的弗兰克矢量心电图测量TWA的存在。
所述TWA测量方法还包括：在获得变换系数之前，将测量电极连接于受测者，并且产生标量心电图和弗兰克矢量心电图；和基于产生的标量心电图和弗兰克矢量心电图计算所述变换系数。
在推导弗兰克矢量心电图的过程中，当存在从受测者获得的并且所述受测者专用的个人系数时，该个人系数可以被用作变换系数，并且，当不存在个人系数时，群系数可以被用作变换系数，所述群系数是从具有统计效果的人口的非特定数目的人获得的多个变换系数的平均值。
在测量TWA存在的过程中，TWA的存在可以基于从弗兰克矢量心电图获得的矢量大小的波形来测量。
在测量TWA存在的过程中，有助于测量TWA的存在的波形可以从矢量大小的波形中选择，并且TWA存在可以从被选择的波形中测量。
在测量TWA的存在的过程中，在所选择的有助于测量TWA的存在的波形中，可以在奇数次心跳或偶数次心跳中检测对测量TWA的存在有不利影响的异常值，并且可以校正包括所检测的异常值的心跳。
根据本发明，还提供了一种TWA测量系统，其包括：TWA测量心电图仪；和心电图管理系统，其构造成管理与通过TWA测量心电图仪获得的弗兰克矢量心电图相关的信息，并包括个人系数数据库，其构造成提供个人系数作为对TWA测量心电图仪的变换系数。
TWA测量心电图仪的心电图仪控制部件可以从TWA测量心电图仪的个人系数数据库或从心电图管理系统的个人系数数据库得到个人系数。
附图说明
图1为实施例的TWA测量心电图仪的方框图。
图2为利用12导联心电图的实施例的TWA测量心电图仪的操作流程图。
图3为示出在实施例的TWA测量心电图仪中用于获得弗兰克矢量心电图的过程的流程图。
图4为利用波形图示出图3的操作流程图的过程的视图。
图5为示出在实施例的TWA测量心电图仪中用于测量TWA存在的过程的流程图。
图6为示出选择分析片断的过程的视图。
图7为示出选择分析片段的该过程的视图。
图8示出进行校正异常值的过程前后的波形图。
图9为示出计算周期图的过程的视图。
图10为示出计算周期图的该过程的视图。
图11为实施例的TWA测量心电图仪的操作流程图。
图12为实施例的TWA测量系统的方框图。
图13为在实施例中的心电图管理系统的方框图。
图14为在实施例中的病人ID管理部件的方框图。
图15为示出在实施例的TWA测量系统中当TWA测量心电图仪测量TWA存在期间所执行的过程的主要流程图。
图16为在图15的主要流程图中步骤S23(对照病人ID/病人姓名)的子程序流程图。
图17为在图15的主要流程图中步骤S25(搜索个人系数)的子程序流程图。
图18为在图15的主要流程图中步骤S28(搜索群系数)的子程序流程图。
图19为在计算病人的个人系数和所计算的个人系数被存储在个人系数数据库期间所执行的过程的主要流程图。
图20为示出在图19的主要流程图中步骤S33(计算个人系数)的过程的子程序流程图。
图21为示出在图19的主要流程图中步骤S34(在数据库中注册计算的个人系数)的过程的子程序流程图。
图22为示出在图21的子程序流程图中步骤S41(存储TWA测量心电图的个人系数数据库)的过程的子程序流程图。
图23为示出在图21的子程序流程图中步骤S42(存储心电图管理系统的个人系数数据库)的过程的子程序流程图。
实施例的详细描述
下文中，将参考附图通过举例说明利用12导联心电图的情况来详细描述本发明的TWA测量心电图仪、TWA测量方法和TWA测量系统的实施例。虽然在该实施例中将示范性地描述利用12导联心电图的情况，但是本发明也可以应用于利用包括所获得的导联的3导联标量心电图的情况。
在本发明的TWA测量心电图仪、TWA测量方法和TWA测量系统中，从测量电极的心电图信号产生12导联心电图。将要产生的12导联心电图目标在于通过任意一种测量方法产生的12导联心电图。例如，不仅通常的标准的12导联心电图，而且通过利用转换矩阵推导所获得12导联心电图都被作为目标。另外，就能够产生12导联心电图来说，导联数为二的霍尔特(Holter)心电图和训练心电图也被作为目标。
如上所述，在实施例的TWA测量心电图仪、TWA测量方法和TWA测量系统中，可以从通过任意一种测量方法产生的12导联心电图来获得弗兰克矢量心电图，并且在突发性心脏死亡的预后诊断中有效的TWA的存在可以从所获得的弗兰克矢量心电图来测量。
在实施例中，术语“病人”作为受测者的特定实例来使用。然而，受测者不仅包括在医院接收诊断的病人，还包括除了医院之外的公共设施，例如进行健康检查的体检中心或诊所，或通常的房子的使用者。
[TWA测量心电图的配置]
首先，将描述实施例的TWA测量心电图仪的配置。
图1为实施例的TWA测量心电图仪的方框图。
TWA测量心电图仪110包括病人信息输入部件111、病人信息存储部件112、变换系数存储部件114、测量电极115、显示部件116、心电图仪控制部件117、TWA测量部件118和通信部件119。变换系数存储部件114包括个人系数数据库113A和群系数数据库113B。
病人信息输入部件111用来通过测量人员的键盘操作输入病人信息。具体说，病人信息包含病人ID(专用ID和个人ID)(用于获得个人系数)、病人姓名(用于获得个人系数)、病人年龄(用于获得群系数)、和病人性别(用于获得该群系数)。
病人信息存储部件112存储通过病人信息输入部件111输入的病人信息。例如，病人的专用ID“C123”和个人ID“A123”、病人的姓名、病人的年龄和病人的性别作为病人ID被存储。
变换系数存储部件114存储用于从12导联心电图获得弗兰克矢量心电图的变换系数。
构成变换系数存储部件114的个人系数数据库113A存储先前从特定病人获得的并且对于该受测者是专用的个人系数作为变换系数。个人系数按照时间序列为每一个病人存储。在上述情况下，例如，病人的专用ID“C123”、病人的个人ID“A123”和病人的姓名作为个人系数的附加信息被添加，并且按照时间序列为每一个病人存储。
当获得弗兰克矢量心电图时，通常，利用一个最新的个人系数是足够的。即使这样，多个个人系数按照时间序列被存储，因为存在这样的情况，即，病症的进展或恢复状态可以通过观察个人系数的变化或通过利用过去的个人系数和当前的个人系数产生的弗兰克矢量心电图的变化来了解。
构成变换系数存储部件14的群系数数据库113B存储作为变换系数的群系数，该群系数是为了获得病人的弗兰克矢量心电图而从具有统计效果的人口的非特定数目的人获得的多个变换系数的平均值。群系数为每一种性别和年龄群而储存。
测量电极115是连接于病人的身体表面的电极。当要产生12导联心电图时或当要得到病人专用的个人系数时，利用测量电极115。实施例以由各种测量方法产生的12导联心电图为目标。因此，测量电极115分别连接于在所使用的测量方法中确定的病人的测量部分。
显示部件116在显示装置上显示通过病人信息输入部件111输入的病人信息，和由TWA测量部件118测量的TWA存在的测量结果，或者将它们打印出来。
心电图仪控制部件117从测量电极115的心电图信号产生12导联心电图，并且通过利用变换系数存储部件114的变换系数从该12导联心电图推导弗兰克矢量心电图。当从病人获得的个人系数存在于个人系数数据库113A中时，心电图仪控制部件117利用个人系数作为变换系数来推导弗兰克矢量心电图。这是因为，当利用病人的个人系数时，TWA的测量的精确性被提高。相反，当从病人获得的个人系数不存在于个人系数数据库113A中时，为了在不存在该病人专用的个人系数时能够进行TWA的测量，利用存在于群系数数据库113B中的群系数作为变换系数来推导弗兰克矢量心电图。
当在个人系数数据库113A中不存在该病人专用的个人系数时，心电图仪控制部件117使显示部件116显示用于提示获得病人的个人系数的消息。这是因为获得个人系数的提示有助于提高测量精确性。在要获得病人专用的个人系数的情况下，心电图仪控制部件117从连接到病人身体表面的测量电极115的心电图信号计算该病人专用的个人系数，并将计算的个人系数存储在个人系数数据库113A中。
在得到病人专用的个人系数的情况下，为了产生12导联心电图，首先，测量电极115分别连接于在每种测量方法中确定的病人的测量部分，并且从测量电极115的心电图信号产生12导联心电图。其次，测量电极115分别连接到为产生弗兰克矢量心电图而确定的病人的测量部分，并且从测量电极115的心电图信号产生弗兰克矢量心电图。测量电极115连接于用于产生12导联心电图的测量部分与测量电极115连接于用于产生弗兰克矢量心电图的测量部分不同。12导联心电图的产生与弗兰克矢量心电图的产生可以分开进行。然而，在实施例中，为了提高工作效率心电图的产生同时进行。心电图仪控制部件117计算病人专用的并且用来将12导联心电图与弗兰克矢量心电图匹配的个人系数。心电图仪控制部件117为每一个病人并按照时间序列将计算的个人系数存储在个人系数数据库113A中，同时将病人的专用ID、病人的个人ID和病人的姓名作为个人系数的附加信息添加到个人系数中。
心电图仪控制部件117通常控制上述各种操作，并且还控制TWA测量心电图仪110的所有操作。心电图仪控制部件117包括用于实现TWA测量心电图仪110的所有操作的程序。心电图仪控制部件117的操作将在稍后详细描述。
TWA测量部件118从由心电图仪控制部件117获得的弗兰克矢量心电图来测量TWA的存在。TWA测量部件118从形成弗兰克矢量心电图的波形中选择有助于测量TWA存在的波形，并且从所选择的波形测量TWA的存在。这是因为，当有助于测量TWA存在的波形被选择时，TWA测量的精确性被提高。TWA测量部件118的操作将在稍后详细描述。
通信部件119将病人的个人系数传输给心电图管理系统(EMS)140(参见图12)，并且相反地从EMS140接收由心电图仪控制部件117请求的病人的个人系数。在通信部件119中的个人系数的传输/接收可以通过医院内部网络160来执行。
[TWA测量心电图仪的操作]
接下来，将参考图2的流程图描述实施例的TWA测量心电图仪的操作。
在图2的流程图中，步骤S1到S3由TWA测量心电图仪的操作者(测量人员)来执行，而步骤S4到S8的操作由心电图仪控制部件117来执行。步骤S1到S8的操作还与实施例的TWA测量方法的过程对应。
步骤S1
为了产生12导联心电图，图1中所示的TWA测量心电图仪110的操作者将测量电极115连接到病人身体表面的预定部分。由于实施例的TWA测量心电图仪110以由各种测量方法产生的12导联心电图为目标，因此测量电极115分别连接到在使用的测量方法中确定的病人的测量部分。例如，在将要产生推导的所述类型的12导联心电图的情况下，其中V1、V3、V4和V6从4个导联I、II、V2和V5获得，测量电极115连接到总共六个位置，即，4个位置或左右胳膊(电极L、R)和左右下肢(电极LL、RL)以便获得导联I和II的电极信号，和两个位置或第四肋间空间的左下胸骨边缘，和左前腋窝一线与穿过第五肋间空间的水平线的交点以便获得两个胸导联(导联V2和导联V5)心电图信号。在标准12导联心电图将被产生的情况下，10个测量电极115连接到总共10个位置，即六个位置用于测量胸导联，和4个位置用于测量4个肢体导联。
步骤S2
接下来，操作者通过病人信息输入部件111输入病人信息。例如，输入专用ID“C123”和个人ID“A123”作为病人ID，然后输入病人的姓名、病人的年龄和病人的性别。
步骤S3
然后，操作者接通TWA测量心电图仪110的测量开关(未示出)。当测量开关接通时，测量TWA的存在就开始。
步骤S4
心电图仪控制部件117从在步骤S1连接到病人的六个测量电极115的心电图信号来产生12导联心电图。
在推导的12导联心电图将被产生的情况下，四导联，即导联I和II、和两个胸导联(导联V2和V5)的心电图信号代入下面的矩阵公式，并且将四个导联的心电图信号与转换矩阵相乘，从而产生实际上不被测量电极115检测的剩余四个胸导联，即导联V1、导联V3、导联V4和导联V6(推导的导联矢量)的推导的心电图。
V 1 V 3 V 4 V 6 = T 1 I T 1 II T 12 T 15 T 3 I T 3 II T 32 T 35 T 4 I T 4 II T 42 T 45 T 6 I T 6 II T 62 T 65 VI VII V 2 V 5
最后，产生12导联心电图。
在标准12导联心电图心电图将被产生的情况下，基于由十个测量电极检测的心电图信号，计算标准12导联的六个肢体导联波形(I、II、III、aVR、aVL和aVF)和标准12导联的六个胸导联波形(V1、V2、V3、V4、V5和V6)，并且最终产生12导联心电图。
步骤S5到S7
接下来，心电图控制部件117从变换系数存储部件114获得用来推导弗兰克矢量心电图的变换系数。
心电图仪控制部件117判断与在步骤S2中输入的个人信息一致的个人系数是否存在于个人系数数据库13A中。具体说，判断与从病人信息输入部件111输入的专用ID“C123”和个人ID“A123”一致的个人系数是否存在于个人系数数据库113A中。
如果与病人信息一致的个人系数在个人系数数据库113A中，那么心电图仪控制部件117从个人系数数据库113A获得病人的个人系数。个人系数为相应的病人并按照时间序列进行存储。在存在用于该病人的多个个人系数的情况下，获得最新的个人系数。
相反，如果与病人信息一致的个人系数不在个人系数数据库113A中，那么心电图仪控制部件117从群系数数据库113B中获得与作为病人信息输入的病人的年龄和性别一致的群系数。
步骤S8
心电图仪控制部件117通过利用在步骤S5到S7中获得的病人专用的个人系数或者群系数其中之一的变换系数，从在步骤S4中产生的12导联心电图来推导弗兰克矢量心电图。
推导弗兰克矢量心电图的过程将参考图3的流程图具体地描述。
步骤9
TWA测量部件118从在步骤S8中推导的弗兰克矢量心电图来测量TWA的存在。
测量TWA存在的过程将参考图5的操作流程图具体地描述。
图3的操作流程图示出获得弗兰克矢量心电图的过程。该流程图为图2中的步骤S8的子程序流程图。
步骤S8-1
心电图仪控制部件117获得在步骤S4中产生的12导联心电图，和在步骤S6中获得的个人系数或在步骤S7中获得的群系数其中任何一个，作为变换系数。
步骤S8-2
心电图仪控制部件117将12导联心电图与变换系数相乘以获得弗兰克矢量心电图的六个导联I、E、C、A、H和M。
在变换系数相乘之前，心电图仪控制部件117进行用于消除产生的12导联心电图的基线漂移预处理，以便12消除导联心电图的DC分量以与基线对齐。同时，进行用于消除12导联心电图的高频分量的预处理以便消除高频噪声，从而平滑12导联心电图的波形。
弗兰克矢量心电图可以由于以下的原理通过将12导联心电图与变换系数相乘来获得。
在12导联心电图L和弗兰克矢量心电图F之间，存在例如下面的公式1中所示的关系。在该公式中，α表示变换系数。
Lα＝F…公式1
具体说，在公式1中，L表示12导联心电图的排列(n×8)，α表示使用个人系数或群系数的变换系数，F表示弗兰克矢量心电图的排列(n×6)。在这个实例中，在12导联心电图将被测量的情况下，利用八个电极以便测量导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6。因此，12导联心电图构成(n×8)排列。
如公式1所示，为了从12导联心电图推导弗兰克矢量心电图的六个导联I、E、C、A、H和M，需要将12导联心电图与变换系数α相乘。
为了从12导联心电图获得弗兰克矢量心电图，必需先获得变换系数α。因为这个，实施例的TWA测量心电图仪110具有两种变换系数，即病人专用的个人系数，和其为从具有统计效果的人口的非特定数目的人中获得的多个变换系数的平均值的群系数。
当得到病人专用的个人系数时，例如，测量电极115连接于用于产生12导联心电图L(导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6)所确定的病人的特定部分，并且12导联心电图从测量电极115的心电图信号产生。接下来，测量电极115连接到用于产生弗兰克矢量心电图F(导联I、导联E、导联C、导联A、导联H和导联M)所确定的病人的特定部分，并且弗兰克矢量心电图从测量电极115的心电图信号产生。在该实例中，将所产生的12导联心电图L(导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6)和所产生的弗兰克矢量心电图F(导联I、导联E、导联C、导联A、导联H和导联M)代入下面的公式2中，并且获得作为变换系数α的病人专用的个人系数。用于获得作为变换系数的个人系数的过程将参考图11的流程图专门详细描述。
α＝(LTL)-1LTF…公式2
其中LT表示12导联心电图L的转置矩阵，而(LTL)-1表示(LTL)的逆矩阵。
如公式1所示，心电图仪控制部件117将12导联心电图L与获得的变换系数相乘，以获得弗兰克矢量心电图F的导联I、导联E、导联C、导联A、导联H和导联M。
步骤S8-3
当推导弗兰克矢量心电图F(导联I、导联E、导联C、导联A、导联H和导联M)时，于是心电图仪控制部件117将导联I、导联E、导联C、导联A、导联H和导联M的值代入下面的公式3中，从而获得弗兰克矢量心电图的导联VX、导联VY和导联VZ。
VX＝0.61A+0.171C-0.781I，
VY＝0.655F+0.345M-0.1H，
VZ＝0.133A+0.736M-0.264I-0.374E-0.231C。
以上公式表示公式3。
步骤S8-4
最后，心电图仪控制部件117将获得的导联VX、导联VY和导联VZ的值求平方，然后将平方值相加，并且计算平方值的总和的平方根，从而计算矢量大小VM的值。
上文中，已经描述了推导弗兰克矢量心电图时所进行的过程程。
下面将参考图4的波形图来被描述图3的操作流程图的过程。
当导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6的八个心电图波形作为在图4中的上面的波形图中所示的12导联心电图L被获得时，首先，对心电图波形施加消除基线漂移和高频分量的处理，并且变换系数α与最终得到的波形相乘。然后，获得导联I、导联E、导联C、导联A、导联H和导联M的值，并且获得弗兰克矢量心电图的导联VX、导联VY和导联VZ的心电图波形，如图4所示的中间的波形图所示。最后，如图4所示的下面的波形图所示，弗兰克矢量心电图的矢量大小VM的心电图波形从导联VX、导联VY和导联VZ的值获得。最终获得的TWA的存在将从矢量大小VM的心电图波形来获得。
图5的操作流程图示出用于测量TWA的存在的过程。该流程图为在图2中的步骤S9的子程序流程图。
步骤S9-1
TWA测量部件118从心电图仪控制部件117获得如上所述获得的弗兰克矢量心电图的矢量大小VM。TWA测量部件118从多个矢量大小VM的波形中选择似乎有助于测量TWA的存在的波形，如下所述。
TWA测量部件118从矢量大小VM的波形中选择分析片段。进行该分析片段的选择以便提高TWA存在测量的精确性。
如图6所示，通过从心电图仪控制部件117获得的矢量大小VM中选择128次心跳的ST-T片段来选择分析片段。与图6的情况不同，连续128次心跳的ST-T片段可能不被选择，而是仅仅可能从N次心跳的ST-T片段来被选择似乎容易产生TWA的128次心跳的ST-T片段。在连续的128次心跳的ST-T片段中，存在产生很多噪声并且TWA的存在很难被测量的心跳，因此测量精确性受到影响。因此，似乎容易产生TWA的128次心跳被不连续地选择。
步骤S9-2
TWA测量部件118进行校正异常值的处理。
在分组为奇数次心跳和偶数次心跳的同时，检测似乎对TWA存在的测量产生不利影响的异常值，并且对检测的心跳上进行校正。
具体说，对于奇数和偶数次心跳获得心电图波形的中间值，并且对于奇数和偶数心跳获得标准偏差。计算标准偏差作为参数使用的并且用于确定异常值的阈值，确定关于阈值的大小关系，并且确定异常值。确定为异常值的数值代入标准偏差作为参数使用的校正函数中，并且计算校正值。确定这样的阈值和校正函数的各种方法是已知的。在实施例中，可以利用任何已知的方法或函数。
图8示出进行校正异常值处理前后的波形图。在没有经历校正处理的原始信号的波形中，如图8的上面的波形图所示，中间部分存在电平急速下降的部分。相反，在经历了校正处理的校正信号的波形中，如图8的下面的波形图所示，在中间部分不存在电平急速下降的部分。可以发现校正异常值的处理达到了显著效果。
步骤S9-3
在异常值校正之后，TWA测量部件118对128次心跳的数据执行FFT处理，并且对处理结果求平均以计算周期图。
图9和10为示出周期图的视图。图9和10的上部波形图分别示出奇数和偶数次心跳的平均波形。图9和10的下部波形图示出在计算周期图之后获得的波形图。在计算周期图之后获得的图9波形示出，当频率(循环/心跳)为0.5时，矢量大小值为小。这意味着不存在TWA。相反，在计算周期图之后获得的图10的波形示出，当频率(循环/心跳)为0.5时，矢量大小值为大。这意味着TWA的值为大，并且存在TWA的可能性为大。
如上所述，当计算周期图时，TWA的存在可以根据矢量大小的值的大小在一定程度上被预测。
步骤S9-4
接下来，TWA测量部件118计算交替。在图9和10所示并且在计算周期图之后获得的波形图中，循环/心跳频率从0.44到0.49的区域被确定为噪声频带，并且获得该区域的平均值SNB和标准偏差σNB。当循环/心跳频率为0.5时获得的值表示为S0.5，并且计算下面的公式4，从而计算交替Valt。
Valt＝(S0.5-SNB)1/2…公式4
步骤S9-5
TWA测量部件118确定TWA的存在。通过利用在步骤S9-4中计算的标准偏差σNB和交替值Valt，来计算下面公式5，从而计算交替比k。
k＝Valt/σNB…公式5
然而，TWA的存在根据交替值Valt和交替比k来确定。用于确定TWA存在的条件为交替值Valt＞1.9μV，和交替率k＞3。当满足该确定条件时，可以判定TWA存在。
形时，能够确定TWA的存在，其中具有不同形状的T波在每个心跳(ABABAB…)交替出现。
接下来，在病人的个人系数没有被存储的情况下，实施例的TWA测量心电图仪110使显示部件116显示用于提示获得个人系数的消息。此刻，TWA测量心电图仪110的操作者获得病人的个人系数，并且将其存储在TWA测量心电图仪110的个人系数数据库113A中。这种情况的过程将参考图11的操作流程图详细描述。在图11的流程图中，步骤S10到S12的操作由TWA测量心电图仪110的操作者来执行，而步骤S13到S15的操作由心电图仪控制部件117来执行。步骤S13到步骤S15的操作与实施例的TWA测量方法的过程相对应。
步骤S10
为了获得12导联心电图和弗兰克矢量心电图，首先，操作者将测量电极连接到病人的特定部分。
具体说，测量电极115连接到确定的病人的特定部分以便产生12导联心电图(导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6)，并且连接到确定的病人的特定部分以便产生弗兰克矢量心电图(导联I、导联E、导联C、导联A、导联H和导联M)。
步骤S11
接下来，操作者通过病人信息输入部件111输入病人信息。例如，专用ID“C123”和个人ID“A123”作为病人ID被输入，然后输入病人的姓名和年龄。
步骤S12
然后，操作者接通TWA测量心电图仪110的测量开关(未显示)。
步骤S13
心电图仪控制部件117从连接到病人的测量电极获得测量部分的心电图信号。
步骤S14
心电图仪控制部件117从测量电极115的心电图信号产生12导联心电图，并且从测量电极115的心电图信号还产生弗兰克矢量心电图。在上述示例中，所产生的12导联心电图(导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6)和所产生的弗兰克矢量心电图(导联I、导联E、导联C、导联A、导联H和导联M)被代入下面的公式6中，并且计算病人专用的个人系数作为变换系数α。
α＝(LTL)-1LTF…公式6
其中LT表示12导联心电图L的转置矩阵，而(LTL)-1表示(LTL)的逆矩阵。
步骤S15
心电图仪控制部件117将计算的个人系数存储在个人系数数据库113A中。当个人系数将被存储时，病人的专用ID“C123”、病人的个人ID“A123”，和病人的姓名和年龄作为个人系数的附加信息被添加并按照时间序列存储。
在上文中，已经描述了实施例的TWA测量心电图仪和TWA测量方法的操作。
[TWA测量系统的配置]
接下来，将描述实施例的TWA测量系统的配置。
TWA测量系统500具有：设置在医院A的心电图管理设备100、设置在医院B的心电图管理设备200、在医院A和B中管理病人ID(专用ID和个人ID)的病人ID管理部件300、和连接线400。心电图管理设备100、200和病人ID管理部件300通过连接线400相互连接。
在该实施例中，将举例说明心电图管理设备100、200设置在医院中的情况。设置心电图管理设备100、200的机构不限于医院。例如，除了医院之外，它们可以设置在诸如进行健康检查的体检中心、诊所、学校和高级护理机构。尽管病人ID管理部件300设置在心电图管理设备100、200之外，但是部件可以设置在，例如，医院A的心电图管理设备100中。连接线400可以是通过心电图管理设备100、200连接到病人ID管理部件300的专用线路、或者进行安全性测量的有线或无线因特网线路。
医院A的心电图管理设备100包括TWA测量心电图仪110、120、130；心电图管理系统(EMS)140和显示装置150。TWA测量心电图仪110、120、130、心电图管理系统(EMS)140和显示装置150通过医院内部网160相互连接。医院B的心电图管理设备200包括TWA测量心电图仪210、220、230、心电图管理系统(EMS)240和显示装置250。TWA测量心电图仪210、220、230、心电图管理系统(EMS)240和显示装置250通过医院内部网260相互连接。医院内部网160、260可以是TWA心电图仪110、120、130、210、220、230和显示装置150、250通过它相互连接的专用线路，或有线或无线内部线路，或进行安全性测量的有线或无线因特网线路。
如上所述，TWA测量心电图仪110、120、130、210、220、230具有搜索病人的个人系数或群系数、计算并存储病人的个人系数、和测量TWA存在的功能。因此，用于完成这些功能的程序安装在TWA测量心电图仪110、120、130、210、220、230中。
心电图管理系统140通过医院内部网络160连接到TWA测量心电图仪110、120、130。心电图管理系统140传输并接收与TWA测量有关的信息，例如病人信息、产生的心电图数据以及来自和发往TWA测量心电图仪110、120、130的个人系数，并且管理这些信息。心电图管理系统240通过医院内部网络260连接到TWA测量心电图仪210、220、230。心电图管理系统240传输并接收与TWA测量有关的信息，例如病人信息、产生的心电图数据以及来自和发往TWA测量心电图仪210、220、230的个人系数，并且管理这些信息。
心电图管理系统140、240具有通过连接线400搜索TWA测量心电图仪110、120、130、210、220或230需要的病人的个人系数，并向需要系数的TWA测量心电图仪输出所搜索的个人系数的功能。因此，用于完成该功能的程序安装在心电图管理系统140、240中。
显示装置150接收由TWA测量心电图仪110、120、130产生的以及来自心电图管理系统140与病人的TWA测量有关的信息，并且显由示操作者指定的病人信息。而且，显示装置150显示由TWA测量心电图仪110、120、130获得的并且与TWA测量有关的信息。显示装置250接收由TWA测量心电图仪210、220、230产生的以及来自心电图管理系统240与病人TWA的测量有关的信息，并且显示操作者指定的病人信息。而且，显示装置显示由TWA测量心电图仪210、220、230获得的并且与TWA测量有关的信息。
病人ID管理部件300具有如下功能：当心电图管理系统140、240相互搜索特定病人的个人系数时，将由另一家医院提供的相同病人的个人ID转换为在该医院使用的个人ID。例如，对于相同的病人，假定有在医院A提供个人ID为“A123”，而在医院B提供个人ID为“B456”的情况。在这种情况下，个人ID“A123”转换为“B456”，而个人ID“B456”被转换为“A123”。在医院之间的个人ID的转换中，专用ID被使用。换句话说，专用ID是提供给病人的唯一ID。ID能够在全世界所有的医院通用，并且是该病人专用的。不能为同一个病人提供两个或两个以上的ID。病人ID管理部件300具有专用ID和个人ID对照表。对照表响应心电图管理系统140、240的请求不断被更新。
图13为在实施例中心电图管理系统的方框图。下文中，将描述心电图管理系统140的配置。心电图管理系统240以与心电图管理系统140相同的方式配置。
心电图管理系统140包括病人信息存储部件142、个人系数数据库143、个人系数管理部件147、和通信部件148。
个人信息存储部件142存储病人的信息，该信息与存储在TWA测量心电图仪110、120、130的病人信息存储部件(例如，在TWA测量心电图仪110中的病人信息存储部件112)中的信息相同。病人信息存储部件142集中地存储在连接于医院内部网160的TWA测量心电图仪110、120、130中所存储的所有信息。
个人系数数据库143存储个人系数，该系数与在存储在TWA测量心电图仪110、120、130的个人系数数据库(例如，在TWA测量心电图仪110中的个人系数数据库113A)中的个人系数相同。个人系数数据库143集中地存储在连接于医院内部网160的TWA测量心电图仪110、120、130中所存储的所有个人系数。
个人系数管理部件147通常控制心电图管理系统140的操作。个人系数管理部件147具有如下功能：将病人信息存储在病人信息存储部件142中、将病人的个人系数存储在个人系数数据库143中并且从个人系数数据库143获取病人的个人系数。个人系数管理部件147的操作将被详细描述。
通信部件148接收来自TWA测量心电图仪110、120、130的病人信息和个人系数，并且将个人系数传输给TWA测量心电图仪110、120、130。而且，通信部件148将病人ID传输给病人ID管理部件300，并接收转换后的病人ID。病人信息和个人系数的传输和接收通过医院内部网160来执行，并且病人ID的传输和接收通过连接线400来执行。
图14为实施例中的病人ID管理部件的方框图。
病人ID管理部件300包括病人ID数据库310、病人ID转换部件320和通信部件330。
病人ID数据库310对应地存储病人的专用和个人ID。在专用ID为“C123”，医院A的个人ID为“A123”，而在医院B的个人ID为“B456”的情况下，例如“C123”-“A123”和“C123”-“B456”被存储。
病人ID转换部件320在参考病人的专用ID时将病人ID转换为在各自医院的ID。
在医院A将要搜索医院B的个人系数的情况下，例如，参考专用ID“C123”进入病人ID数据库310，并且医院A中的个人ID“A123”被转换为医院B中的个人ID“B456”。
通信部件330接收从心电图管理系统140、240传输的病人的专用和个人ID，并且将在病人ID转换部件320中转换的个人ID传输给心电图管理系统140、240，
实施例的TWA测量系统具有上述配置。
[TWA测量系统的操作]
接下来，将参考图15到图23的操作流程图描述实施例的TWA测量系统的操作。将简要描述该操作，该操作与TWA测量心电图仪的操作相同。
图15为示出在TWA测量心电图仪110输出病人的TWA测量结果期间进行的过程的主要流程图。该主要流程由TWA测量心电图仪110的心电图仪控制部件117来执行。
步骤S21
心电图仪控制部件117获得从病人信息输入部件111提供的病人信息。病人信息包含病人ID(专用ID和个人ID)、病人的姓名、病人的年龄和病人的姓别。病人ID和病人的姓名被用于搜索病人的个人系数，而病人的年龄和性别用于搜索群系数。
步骤S22
心电图仪控制部件117从连接到病人的测量电极115获得测量部分的心电图信号。
步骤S23
心电图仪控制部件117将病人ID和病人姓名与存储在病人信息存储部件112中的病人ID和病人姓名进行对照。
步骤S24、S25
如果在病人信息存储部件112中存在对照的病人ID和对照的病人姓名，那么确定该对照为OK(步骤S24：是)，并且存储在个人系数数据库113A或143或心电图管理系统240的个人系数数据库的个人系数被搜索。
步骤S26、S27、S30
如果在个人系数数据库113A、143或心电图管理系统240的个人系数数据库的任何一个中存在病人的个人系数(步骤S26：是)，心电图仪控制部件117从个人系数数据库113A、143或心电图管理系统240的个人系数数据库之一获得病人的个人系数。心电图仪控制部件117通过利用获得的个人系数处理由多个测量电极115获得的心电图信号，并测量TWA测存在。
步骤S28
相反，如果在病人信息存储部件112中不存在对照的病人ID和对照的病人姓名(步骤S24：否)，或者如果在个人系数数据库113A、143或心电图管理系统240的个人系数数据库的任何一个中不存在个人系数(步骤S26：否)，那么心电图仪控制部件117搜索存储在变换系数存储部件114中的群系数。在检查病人的年龄和性别的同时，心电图仪控制部件117从变换系数存储部件114中获得对病人来说最优的群系数。
步骤S29
心电图仪控制部件117通过利用获得的群系数处理由多个测量电极检测到的心电图信号，并且测量TWA的存在。
步骤S30
心电图仪控制部件117将TWA存在的测量结果输出给显示部件116。
图16为在图15中的主要流程图的步骤S23(对照病人ID/病人姓名)的子程序流程图。该子程序流程图由心电图仪控制部件117来执行。
步骤S23-1
心电图仪控制部件117获得从病人信息输入部件111提供的病人信息中的病人ID和病人姓名。
步骤S23-2
心电图仪控制部件117搜索存储在病人信息存储部件112、142中的病人ID。
步骤S23-3、S23-4
如果在病人信息存储部件112或142任何一个中，存在与病人信息输入部件111提供的病人ID一致的病人ID(步骤S23-3：是)，那么心电图仪控制部件117搜索存储在病人信息存储部件112中的病人姓名。
步骤S23-5、S23-6
如果在病人信息存储部件112或142任何一个中，存在与病人信息输入部件111提供的病人姓名一致的病人姓名(步骤S23-5：是)，那么指示对照为OK的信号被输出。
步骤S23-3、S23-5、S23-7
如果在病人信息存储部件112或142任何一个中，不存在与病人信息输入部件111提供的病人ID一致的病人ID(步骤S23-3：否)，或者在病人信息存储部件112或142任何一个中，不存在与病人信息输入部件111提供的病人姓名一致的病人姓名(步骤S23-5：否)，那么指示对照为NG的信号被输出，并且心电图仪控制部件117使显示部件116显示消息“产生个人系数，请执行测量和注册”。
即，如果输入的病人ID和病人姓名与在TWA测量心电图仪110或心电图管理系统140中注册的一致，那么可以确定对照为OK，并且，如果病人ID和病人姓名之一不一致，那么可以确定对照为NG。在对照为NG的情况下，病人的个人系数不存储在个人系数数据库113A、143中，因而输出用于提示获得个人系数的消息。
图17为在图15的主流程图中的步骤S25(搜索个人系数)的子程序流程图。该子程序流程图由心电图仪控制部件117来执行。
步骤S25-1
心电图仪控制部件117在TWA测量心电图仪110的个人系数数据库113A中搜索与病人ID一致的个人系数。
步骤S25-2、S25-3
如果在TWA测量心电图仪110中存在病人的个人系数(步骤S25-2：是)，那么心电图仪控制部件117从个人系数数据库113A获得个人系数。
步骤S25-2、S25-4
相反，如果在TWA测量心电图仪110中不存在病人的个人系数(步骤S25-2：否)，那么心电图仪控制部件117在心电图管理系统140的个人系数数据库143搜索与病人ID一致的个人系数。即，在心电图管理系统140中搜索是否存在病人的个人系数。
步骤S25-5、S25-3
如果在心电图管理系统140中存在病人的个人系数(步骤S25-5：是)，那么心电图仪控制部件117从个人系数数据库143中获得个人系数。
步骤S25-5、S25-6
相反，如果在心电图管理系统140的个人系数数据库143中不存在病人的个人系数(步骤S25-5：否)，那么心电图仪控制部件117在心电图管理系统140、240中搜索与病人ID一致的个人系数。在该实施例中，从心电图管理系统240搜索与和病人ID和TWA测量结果的类型一致的个人系数。在心电图管理系统140、240中搜索个人系数通过在利用病人的专用ID的同时利用由病人ID转换部件320转换的个人ID来执行。
步骤S25-7、S25-3
如果在心电图管理系统240中存在病人的个人系数(步骤S25-7：是)，那么心电图仪控制部件117从心电图管理系统240的个人系数数据库获得个人系数。
步骤S25-8
如果在心电图管理系统240的个人系数数据库中不存在病人的个人系数，那么心电图仪控制部件117使显示部件显示“产生个人系数，请执行测量和注册”的消息。
图18为在图15的主要流程图中步骤S28(搜索群系数)的子程序流程图。该子程序流程图由心电图仪控制部件117来执行。
步骤S28-1
心电图仪控制部件117在TWA测量心电图仪110的群系数数据库113B中搜索与病人的年龄和性别一致的群系数。即，在TWA测量心电图仪110中搜索是否存在病人的群系数。
步骤S28-2
心电图仪控制部件117从TWA测量心电图仪110的群系数数据库113B获得群系数。
图19为在计算病人的个人系数和所计算的个人系数存储到个人系数数据库期间执行的过程的主要流程图。该主要流程图由心电图仪控制部件117来执行。
步骤S31
心电图仪控制部件117获得从病人信息输入部件111提供的病人信息。病人信息包含病人ID(专用ID和个人ID)、病人的姓名、病人的年龄和病人的性别。在病人信息中，病人ID和病人的姓名被用来在个人系数数据库中存储病人的病人系数。
步骤S32
心电图仪控制部件117从连接到病人的测量电极115获得任选测量部分的心电图信号。具体说，例如，测量电极115连接到确定的病人的特定部分以产生12导联心电图(导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6)，和连接到确定的病人的特定部分以产生弗兰克矢量心电图(导联I、导联E、导联C、导联A、导联H和导联M)。心电图仪控制部件117从测量电极115测量用于12导联心电图的导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6，和用于弗兰克矢量心电图的导联I、导联E、导联C、导联A、导联H和导联M。
步骤S33
心电图仪控制部件117将由测量电极115测量的导联代入上述公式2中以计算病人的个人系数α。
步骤S34
接下来，心电图仪控制部件117对每一个病人并按照时间序列将计算的个人系数存储在TWA测量系统的数据库中。
步骤S35
心电图仪控制部件117将在步骤S31中获得的病人信息存储在TWA测量心电图仪110的病人信息存储部件112中，同时将该信息存储在心电图管理系统140的病人信息存储部件142中。当病人信息被存储到病人信息存储部件142中时，心电图仪控制部件117通过医院内部网络160从通信部件148传输病人信息，并且个人系数管理部件147通过心电图管理系统140的通信部件148接收病人信息。接下来，由个人系数管理部件147接收的病人信息被存储在病人信息存储部件142中。
步骤S36
然后，心电图仪控制部件117将在步骤S31中获得的病人信息的病人ID(专用ID和个人ID)存储在病人ID管理部件300的病人ID数据库310中。当病人ID将被存储到病人ID数据库310中时，心电图仪控制部件117通过医院内部网160和连接线400从通信部件148传输病人ID，并且病人ID转换部件320通过病人ID管理部件300的通信部件330接收病人ID。接下来，由病人ID转换部件320接收的病人ID被存储到病人ID数据库310中。
图20为图19的主流程图的步骤S33(计算个人系数)的过程的子程序流程图。该流程图由心电图仪控制部件117来执行。
步骤S33-1
心电图仪控制部件117将由连接到病人的测量电极115测量的12导联心电图(导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6)，和弗兰克矢量心电图(导联I、导联E、导联C、导联A、导联H和导联M)代入上述公式2中。
步骤S33-2
接下来，心电图仪控制部件117从导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6的值，和导联I、导联E、导联C、导联A、导联H和导联M的值获得变换系数的值。
步骤S33-3
最后，心电图仪控制部件117将病人信息添加到所获得的变换系数以计算个人系数。具体说，被添加的病人信息包含病人ID(专用ID和个人ID)和病人姓名。在上述情况下，病人的专用ID“C123”、病人的个人ID“A123”和病人的姓名将作为附加信息被添加，并且计算个人系数。
图21为图19的主流程图中的步骤S34(在数据库中注册计算的个人系数)的过程的子程序流程图。该子程序流程图由心电图仪控制部件117和个人系数管理部件147来执行。
步骤S41
首先，心电图仪控制部件117在个人系数数据库113A中存储计算后的个人系数。
步骤S42
其次，个人系数管理部件147将从心电图仪控制部件117发送的个人系数存储在个人数据库143中。这些步骤的过程详细地示于在图22和图23的子程序流程图中。
图22为图21的子程序流程图中的步骤S41(在TWA测量心电图仪的个人系数数据库中存储)的过程的子程序流程图。该流程图由心电图仪控制部件117来执行。
步骤S41-1
首先，心电图仪控制部件117获得病人信息输入部件111提供的病人信息。病人信息包括病人ID(专用ID和个人ID)、病人的姓名、病人的年龄和病人的性别。
步骤S41-2
心电图仪控制部件117将该病人ID与存储在病人信息存储部件112中的ID进行对照。
步骤S41-3、S41-4
如果在病人信息存储部件112中存在对照的病人ID，那么可以确定对照为OK(步骤S41-3：是)，并且心电图仪控制部件117附加地将病人的个人系数地存储在个人系数数据库113A中。即，在相同的病人已经具有个人系数的情况下，新计算的个人系数按照时间序列被附加地存储。当个人系数按照时间序列被更新时，能够总是选择病人的最佳系数，以便能够产生TWA测量的精确结果。
步骤S41-3、S41-5
相反，如果在病人信息存储部件112中不存在对照的病人ID，那么可以确定对照为NG(步骤S41-3：否)，病人的个人系数不存在于个人系数数据库113A中，因此心电图仪控制部件117将个人系数最新地存储在个人系数数据库113A中。
图23为图21的子程序流程图中的步骤S42(存储在心电图管理系统的个人系数数据库中)的过程的子程序流程图。该流程图由个人系数管理部件147来执行。
步骤S42-1
首先心电图管理系统140的个人系数管理部件147从TWA测量心电图仪110获得病人信息。病人信息包含病人ID(专用ID和个人ID)、病人的姓名、病人的年龄和病人的姓别。
步骤S42-2
个人系数管理部件147将该病人ID与存储在病人信息存储部件142中的病人ID进行对照。
步骤S42-3、S42-4
如果在病人信息存储部件142中存在对照的病人ID，那么可以确定对照为OK(步骤S42-3：是)，个人系数管理部件147将病人的个人系数附加地存储在个人系数数据库143中。即，在相同的病人已经具有个人系数的情况下，最近计算的个人系数按照时间序列被附加地存储。
步骤S42-3、S42-5
相反，如果在病人信息存储部件142中不存在对照的病人ID，那么可以确定对照为NG(步骤S42-3：否)，病人的个人系数不存在于个人系数数据库143中，因此个人系数管理部件147将个人系数最新地存储在个人系数数据库143中。
如上所述，根据实施例的TWA测量心电图仪、TWA测量方法和TWA测量系统，变换系数的利用使得弗兰克矢量心电图可以从12导联心电图获得。因此，测量人员可以将测量电极连接到与获得12导联心电图的情况下的相同的位置。由于TWA的存在可以从获得的12导联心电图来测量，因此TWA的存在能够以与12导联心电图被获得的情况的相同的方式来测量。而且，病人专用的个人系数可以用来测量病人的TWA的存在，因此能够获得TWA测量的最佳结果，因此能够提高测量的精确性。
已经描述了个人系数对每一个病人并按照时间序列被存储的实施例。在TWA测量心电图仪的存储容量较小的情况下，个人系数可以不按照时间序列存储，而是最新的个人系数可以被更新存储。实施例如此构造：在存在特定病人的个人系数的情况下，例如，即使个人系数在两年或三年之前产生，该病人系数也可以被使用。然而，为了非常精确地检测病人的TWA存在，优选地，TWA测量的结果通过利用在时间上尽可能靠近当前时间所获得的个人系数来产生。因此，可选地，可以设置个人系数的可用周期(例如，从获得开始一年)。在可用周期到期的情况下，可以用与个人系数不存在的情况相似的方式来提示个人系数的获得，并且可以执行诸如在图20的流程图所示的过程。
在该实施例中，进行个人ID的转换以便在医院之间搜索并获得个人系数。个人ID的转换根据专用ID来执行。用于ID转换的各种技术是已知的。用于ID转换的技术不限于在实施例中例举的技术，并且任何已知的各种技术都可以使用。例如，称为OpenPIXPDQ的开放式软件应用程序就是通常所说的用于医院之间的ID转换的软件应用程序。ID转换可以利用这样的开放式软件应用程序来执行。
本发明不限于上述实施例。本发明的技术范围包括上述的实施例和由本领域技术人员能够修改的其他实施例。
根据本发明的一方面，心电图仪控制部件从测量电极的心电图信号产生标量心电图。心电图仪控制部件可以通过利用各种测量方法中的任何一种来产生标量心电图。连接到受测者的测量电极的位置和数量根据用来获得标量心电图的测量方法(例如，通常的标准12导联心电图、推导的12导联心电图、训练心电图、Holter心电图或监视心电图)的不同而不同。心电图仪控制部件通过将利用各种测量方法的任何一种所产生的标量心电图与存储在变换系数存储部件中的变换系数相乘来获得弗兰克矢量心电图。TWA测量部件基于由心电图仪控制部件获得的弗兰克矢量心电图的波形来测量TWA的存在。
根据本发明的一方面，TWA测量心电图仪包括：连接到受测者以便产生标量心电图的测量电极；存储用于获得弗兰克矢量心电图的变换系数的变换系数存储部件；从测量电极的心电图信号产生标量心电图并且通过利用变换系数从标量心电图获得弗兰克矢量心电图的心电图仪控制部件；和基于获得的弗兰克矢量心电图测量TWA存在的TWA测量部件。变换系数存储部件包括：存储作为变换系数从受测者获得的并受测者专用的个人系数的个人系数数据库；和存储作为变换系数的群系数的群系数数据库，该群系数是为了获得受测者的弗兰克矢量心电图从具有统计效果的人口的非特定数量的人获得的多个变换系数的平均值。
根据本发明的一方面，与设置在TWA测量心电图仪中的变换系数存储部件类似，心电图管理系统也包括个人系数数据库。因此，在受测者专用的个人系数不存在于TWA测量心电图仪的情况下，TWA测量心电图仪通过利用从心电图管理系统提供的个人系数作为变换系数来测量TWA的存在。
根据本发明的一方面，可以从标量心电图来获得弗兰克矢量心电图，因此测量的人可以将测量电极连接到与获得标量心电图的情况相同的位置。由于TWA的存在可以从获得的弗兰克矢量心电图测量，因此能够以与获得标量心电图的情况相同的方式来测量TWA存在。

Claims (14)

1.一种TWA测量心电图仪，包括：
变换系数存储部件，该变换系数存储部件构造成存储用于推导弗兰克矢量心电图的变换系数；
心电图仪控制部件，该心电图仪控制部件构造成从适于连接到受测者的测量电极的电信号产生标量心电图，并且构造成通过利用所述变换系数从所述标量心电图推导所述弗兰克矢量心电图；和
TWA测量部件，该TWA测量部件构造成基于所推导的弗兰克矢量心电图来测量TWA的存在，
其特征在于，所述变换系数存储部件包括：
个人系数数据库，该个人系数数据库构造成存储从所述受测者获取的并且对受测者特有的个人系数，作为所述变换系数；和
群系数数据库，该群系数数据库构造成存储群系数作为所述变换系数，所述群系数是从统计学上有效的总体的非特定数量的人获取的多个变换系数的平均值。
2.如权利要求1所述TWA测量心电图仪，其中
当所述个人系数数据库存储所述个人系数时，所述心电图仪控制部件通过利用所述个人系数作为所述变换系数来推导所述弗兰克矢量心电图，并且
当所述个人系数数据库未存储所述个人系数时，所述心电图仪控制部件通过利用存储在所述群系数数据库中的群系数作为所述变换系数来推导所述弗兰克矢量心电图。
3.如权利要求1所述TWA测量心电图仪，其中
当所述个人系数数据库未存储所述个人系数时，所述心电图仪控制部件使显示部件显示用于提示获取所述个人系数的消息。
4.如权利要求1所述TWA测量心电图仪，其中
所述心电图仪控制部件从所述测量电极的心电图信号计算所述个人系数，并且将所计算的个人系数存储在所述个人系数数据库中。
5.如权利要求1所述TWA测量心电图仪，其中
所述TWA测量部件基于从所述弗兰克矢量心电图获得的矢量大小的波形来测量TWA的存在。
6.如权利要求1所述TWA测量心电图仪，其中
所述TWA测量部件从矢量大小的波形选择能够有助于测量TWA的存在的波形，并且从所选择的波形测量TWA的存在。
7.如权利要求6所述TWA测量心电图仪，其中
在所选择的有助于测量TWA的存在的波形中，TWA测量部件检测在奇数次心跳或偶数次心跳中可能对TWA存在的测量具有不利影响的异常值，并且校正包括所检测的异常值的心跳。
8.一种TWA测量方法，包括：
将测量电极连接于受测者；
从所述测量电极的心电图信号产生标量心电图；
获取用于推导弗兰克矢量心电图的变换系数；
通过利用所述变换系数从所述标量心电图推导所述弗兰克矢量心电图；和
从所推导的弗兰克矢量心电图测量TWA的存在，
其特征在于，在推导所述弗兰克矢量心电图的过程中，
当存在从所述受测者获取的并且对所述受测者特有的个人系数时，将所述个人系数用作所述变换系数，并且
当不存在所述个人系数时，将群系数作为所述变换系数，所述群系数是从统计学上有效的总体的非特定数量的人获取的多个变换系数的平均值。
9.如权利要求8所述TWA测量方法，还包括：
在获取所述变换系数之前，将测量电极连接于受测者，并且产生标量心电图和弗兰克矢量心电图；和
基于所产生的标量心电图和弗兰克心电图计算所述变换系数。
10.如权利要求8所述TWA测量方法，其中，在测量TWA的存在的过程中，
基于从所述弗兰克矢量心电图获得的矢量大小的波形来测量所述TWA的存在。
11.如权利要求10所述TWA测量方法，其中，在测量TWA的存在的过程中，
从所述矢量大小的波形选择可有助于测量TWA的存在的波形，并且从所选择的波形测量TWA的存在。
12.如权利要求11所述TWA测量方法，其中，在测量TWA的存在的过程中，
在所选择的可有助于测量TWA的存在的波形中，在奇数次心跳或偶数次心跳中检测可能对测量TWA的存在具有不利影响的异常值，并且校正包括所检测的异常值的心跳。
13.一种TWA测量系统，其特征在于，包括：
根据权利要求1的TWA测量心电图仪；和
心电图管理系统，该心电图管理系统构造成管理与由所述TWA测量心电图仪推导的弗兰克矢量心电图有关的信息，并且包括个人系数数据库，该个人系数数据库构造成为TWA测量心电图仪提供作为变换系数的个人系数。
14.如权利要求13所述TWA测量系统，其中
所述TWA测量心电图仪的心电图仪控制部件从所述TWA测量心电图仪的所述个人系数数据库或所述心电图管理系统的个人系数数据库获取所述个人系数。
CN201210090420.4A 2011-03-29 2012-03-29 Twa测量心电图仪、twa测量方法和twa测量系统 Active CN102727193B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2011-073093 2011-03-29
JP2011073093A JP5631793B2 (ja) 2011-03-29 2011-03-29 Ｔｗａ測定心電計及びｔｗａ測定システム

Publications (2)

Publication Number Publication Date
CN102727193A CN102727193A (zh) 2012-10-17
CN102727193B true CN102727193B (zh) 2016-05-25

Family

ID=46044326

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210090420.4A Active CN102727193B (zh) 2011-03-29 2012-03-29 Twa测量心电图仪、twa测量方法和twa测量系统

Country Status (4)

Country Link
US (1) US9265437B2 (zh)
EP (1) EP2505135B1 (zh)
JP (1) JP5631793B2 (zh)
CN (1) CN102727193B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5926074B2 (ja) * 2012-03-06 2016-05-25 大名 魏 Ｔｗａ測定装置及びｔｗａ測定装置の作動方法
JP6002584B2 (ja) * 2013-01-10 2016-10-05 日本光電工業株式会社 Ｔｗａ測定装置
US10716483B2 (en) * 2013-02-08 2020-07-21 Ivana I. VRANIC Method and system for vector analysis of electrocardiograms
CN103431857A (zh) * 2013-09-09 2013-12-11 苏州百慧华业精密仪器有限公司 一种自动扫描动态心电图可疑twa阳性段的方法
CN107072582B (zh) * 2014-10-22 2020-11-27 株式会社心电技术研究所 脑电波测量装置、脑电波测量方法以及脑电波测量系统
CN110151165B (zh) * 2019-04-30 2022-02-01 杭州电子科技大学 一种基于非线性动力学特征的心电向量图分类方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5713367A (en) * 1994-01-26 1998-02-03 Cambridge Heart, Inc. Measuring and assessing cardiac electrical stability

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5935082A (en) 1995-01-26 1999-08-10 Cambridge Heart, Inc. Assessing cardiac electrical stability
US6804550B1 (en) * 1999-09-29 2004-10-12 Draeger Medical Systems, Inc. Method and apparatus for frank lead reconstruction from derived chest leads
US6453191B2 (en) * 2000-02-18 2002-09-17 Cambridge Heart, Inc. Automated interpretation of T-wave alternans results
EP1304958B1 (en) * 2000-08-03 2008-11-19 Draeger Medical Systems, Inc. An electrocardiogram system for synthesizing leads and providing an accuracy measure
AU2002232954A1 (en) * 2000-10-23 2002-05-06 Cornell Research Foundation Inc. Intracardiac detection and control of repolarization alternans
US6668189B2 (en) 2001-10-05 2003-12-23 Ge Medical Systems Information Technologies, Inc. Method and system for measuring T-wave alternans by alignment of alternating median beats to a cubic spline
US6901285B2 (en) 2002-05-17 2005-05-31 David M. Schreck System and method for synthesizing leads of an electrocardiogram
JP4777326B2 (ja) * 2003-02-26 2011-09-21 大名 魏 付加誘導機能を備えた心電計及び付加誘導心電図導出方法
JP4153950B2 (ja) 2003-02-26 2008-09-24 公立大学法人会津大学 付加誘導機能を備えた心電計及び付加誘導心電図導出方法
GB0317947D0 (en) * 2003-07-31 2003-09-03 Mar Reynolds Medical Del Ltd Reduced electrode electrocardiography system
RS49856B (sr) * 2004-01-16 2008-08-07 Boško Bojović Uređaj i postupak za vizuelnu trodimenzionalnu prezentaciju ecg podataka
US20070260151A1 (en) * 2006-05-03 2007-11-08 Clifford Gari D Method and device for filtering, segmenting, compressing and classifying oscillatory signals
JP5143375B2 (ja) * 2006-05-26 2013-02-13 フクダ電子株式会社 心電図解析装置
US8041416B2 (en) * 2007-09-18 2011-10-18 Harbinger Medical, Inc. Method and apparatus for determining susceptibility for arrhythmias using wedensky modulated electrocardiography and alternans analysis

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5713367A (en) * 1994-01-26 1998-02-03 Cambridge Heart, Inc. Measuring and assessing cardiac electrical stability

Also Published As

Publication number Publication date
US9265437B2 (en) 2016-02-23
JP2012205716A (ja) 2012-10-25
CN102727193A (zh) 2012-10-17
EP2505135B1 (en) 2015-01-07
EP2505135A1 (en) 2012-10-03
JP5631793B2 (ja) 2014-11-26
US20120253213A1 (en) 2012-10-04

Similar Documents

Publication Publication Date Title
Liu et al. An open access database for evaluating the algorithms of electrocardiogram rhythm and morphology abnormality detection
CN102727193B (zh) Twa测量心电图仪、twa测量方法和twa测量系统
Kligfield et al. Recommendations for the standardization and interpretation of the electrocardiogram: part I: the electrocardiogram and its technology: a scientific statement from the American Heart Association Electrocardiography and Arrhythmias Committee, Council on Clinical Cardiology; the American College of Cardiology Foundation; and the Heart Rhythm Society endorsed by the International Society for Computerized Electrocardiology
Anand et al. Explainable AI decision model for ECG data of cardiac disorders
US20130096447A1 (en) System and methods for serial analysis of electrocardiograms
JP5422513B2 (ja) 導出心電図生成システム及び導出心電図生成方法
Hnatkova et al. Clinical value of different QRS-T angle expressions
May et al. The WCT formula: a novel algorithm designed to automatically differentiate wide-complex tachycardias
Xia et al. Automatic detection of ECG electrode misplacement: a tale of two algorithms
Huebner et al. Electrocardiologic and related methods of non-invasive detection and risk stratification in myocardial ischemia: state of the art and perspectives
Kashou et al. The WCT formula II: an effective means to automatically differentiate wide complex tachycardias
US11963800B2 (en) ECG training and skill enhancement
Alkhodari et al. Investigating circadian heart rate variability in coronary artery disease patients with various degrees of left ventricle ejection fraction
AU2003231150A1 (en) System and method for predicting the onset of cardiac pathology using fractal analysis
US8843193B2 (en) TWA measuring apparatus and TWA measuring method
Madias A comparison of 2-lead, 6-lead, and 12-lead ECGs in patients with changing edematous states: implications for the employment of quantitative electrocardiography in research and clinical applications
CN104545882B (zh) 心电图测量装置和合成心电图生成方法
Zusterzeel et al. The 43rd International Society for Computerized Electrocardiology ECG initiative for the automated detection of strict left bundle branch block
US9149201B2 (en) TWA measuring apparatus and TWA measuring method
JP2013183801A5 (zh)
Al-Zaiti et al. The role of automated 12-lead ECG interpretation in the diagnosis and risk stratification of cardiovascular disease
Myerburg et al. Task force 2: training in electrocardiography, ambulatory electrocardiography, and exercise testing
Estes Jr et al. The electrocardiogram in left ventricular hypertrophy: past and future
Jabbour et al. A stepwise approach to reading ECGs using colour-coded electrical viewpoints
Warner et al. Optimized electrocardiographic criteria for prior inferior and anterior myocardial infarction

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20191203

Address after: Tokyo, Japan

Patentee after: Nihon Kohden Corp.

Address before: Chiba, Japan

Co-patentee before: Nihon Kohden Corp.

Patentee before: Wei Daming

TR01 Transfer of patent right