CN103300847A - Twa测量装置和twa测量方法 - Google Patents
Twa测量装置和twa测量方法 Download PDFInfo
- Publication number
- CN103300847A CN103300847A CN2013100854515A CN201310085451A CN103300847A CN 103300847 A CN103300847 A CN 103300847A CN 2013100854515 A CN2013100854515 A CN 2013100854515A CN 201310085451 A CN201310085451 A CN 201310085451A CN 103300847 A CN103300847 A CN 103300847A
- Authority
- CN
- China
- Prior art keywords
- twa
- patient
- electrocardiogram
- coefficient
- measurement
- 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
Links
Images
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
-
- 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/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]
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
一种TWA测量装置包括:心电图仪控制部,该心电图仪控制部构造成从被检对象的心电信号产生心电图;以及TWA测量部,该TWA测量部构造成从心电图选择有助于TWA的测量的至少两个波形,并且该TWA测量部构造成通过利用所选择的至少两个波形测量TWA的存在。
Description
相关申请的交叉引用
本申请基于并要求2012年3月6日提交的在先日本专利申请号2012-049525的优先权,其全部内容通过引用并入此处。
技术领域
当前公开的主题涉及能测量TWA(T波交替)的存在的TWA测量装置和TWA测量方法。
背景技术
TWA出现在诸如QT延长综合症、变异型心绞痛、急性心肌缺血、电解质异常、阵发性心动过速、心动过缓或心包积液的疾病的开始。TWA是其中出现在心电图中的T波的振幅和极性交替改变的现象,并且是有效预测突然心脏死的指标。TWA不是总能用肉眼观察到的现象,并因此TWA在临床中的应用是有限的。
因此,从20世纪80年代起,已研制用于使微小TWA(微伏TWA:MTWA)能够由计算机测量的技术。
用于测量TWA的当前提出的技术的示例是基于通用电气(GE)公司的MMA(修正移动平均)方法的测量技术和基于剑桥心脏(CH)公司的周期图的测量技术,上述测量技术分别在美国专利号6,668,189和美国专利号5,935,082中公开。
GE公司的测量技术涉及一种分析时域中的时间波形的方法,并被认为具有对噪声的抵抗性。然而,该技术作为测量技术没有长的历史,并需要监视其临床效果。
相反,从20世纪80年代起已使用作为频域中的技术的、CH公司的测量技术,由此已证明其在临床中的有效性。因此,现今人们认为基于CH公司的周期图的测量技术在临床上比基于GE公司的MMA方法的测量技术更有用。
相对于基于CH公司的周期图的测量技术,在该测量技术的宣布之后,追加了诸如测量电极的技术的用于进行新过程的各种技术,并且目前仍使用追加的最新技术。
在美国专利No.5,935,082中公开的基于周期图测量TWA的技术基于的假设是满足各种条件,诸如心率在运动试验中必须保持恒定,并且在诸如额外收缩的异常情况下的心率必须被抑制成特定值或较小。
然而,取决于被检对象,存在心率在运动试验中不能保持恒定的情况。存在的另外的情况是即使当心率能大致保持恒定,噪声可能在多次搏动中进入。此外,被检对象不能在诸如额外收缩的异常情况下将心率抑制成特定值或较小。
在美国专利号No.5,935,082中公开的用于测量TWA的技术中,即使试图在不满足以上条件的情况下进行TWA的测量,但并不开始用于测量TWA的数据收集。即使在开始数据收集的情况下,当以上条件在测量中变得不满足时,数据收集复位。
此外,由于在满足以上条件的同时进行测量,所以检查时间倾向于延长,由此将额外的负担强加于测量人和被检对象。除非满足以上条件,否则TWA测量不结束。因此,在时间受限制的情况下,不能完成TWA测量。
发明内容
当前公开的主题可提供TWA测量装置和TWA测量方法,该TWA测量心电图仪和TWA测量方法能减轻对测量人和被检对象的负担,并且能从通过各种测量方法中的任何一种获取的心电图测量TWA的存在。
TWA测量装置可包括:心电图仪控制部,该心电图仪控制部构造成从被检对象的心电信号产生心电图;以及TWA测量部,该TWA测量部构造成从心电图中选择有助于TWA的测量的至少两个波形,并且该TWA测量部构造成通过利用所选择的所述至少两个波形测量TWA的存在。
TWA测量装置还可包括构造成存储变换系数的变换系数存储部。心电图可以是Frank矢量心电图,心电图仪控制部可通过利用变换系数从心电信号合成Frank矢量心电图,以及TWA测量部可从由Frank矢量心电图计算的矢量幅度的波形来选择所述到少两个波形。
变换系数存储部可包括:个人系数数据库,该个人系数数据库构造成存储对于被检对象特定的并且从被检对象获取的个人系数,作为所述变换系数;以及群组系数数据库,该群组系数数据库构造成存储群组系数作为所述变换系数,该群组系数是从统计有效的人口的非特定数量的人所获取的多个变换系数的平均数,以便合成被检对象的Frank矢量心电图。
当从被检对象获取的个人系数存在于个人系数数据库中时,心电图仪控制部可通过利用该个人系数作为变换系数而合成Frank矢量心电图。当在个人系数数据库中不存在个人系数时,心电图仪控制部可通过利用在群组系数数据库中存在的群组系数作为变换系数来合成Frank矢量心电图。
当在个人系数数据库中不存在被检对象的个人系数时,心电图仪控制部可使显示部显示用于提示获取个人系数的消息。
心电图仪控制部可由心电信号计算被检对象的个人系数,并将计算的个人系数存储在个人系数数据库中。
心电图可以是Frank矢量心电图,心电图仪控制部可从心电信号直接产生Frank矢量心电图,并且TWA测量部可从由Frank矢量心电图计算的矢量幅度的波形中选择所述至少两个波形。
心电图可以是标量心电图,心电图仪控制部可以从心电信号合成或者直接产生标量心电图,并且TWA测量部可从标量心电图的波形中选择所述至少两个波形。
标量心电图可以是12导联心电图、Holter心电图、事件心电图、运动心电图和监测心电图中的一种。
TWA测量部可从多次心搏的波形中选择有助于TWA的测量的至少两个波形。
在有助于TWA的测量的所述至少两个波形中包括异常值的情况下,TWA测量部可校正异常值。
TWA测量方法可包括:从被检对象的心电信号产生心电图;从所产生的心电图选择有助于TWA的测量的至少两个波形;以及通过利用所选择的所述至少两个波形测量TWA的存在。
心电图可以是Frank矢量心电图,可通过利用变换系数从心电信号合成Frank矢量心电图,并且可从由Frank矢量心电图计算的矢量幅度的波形中选择所述至少两个波形。
心电图可以是Frank矢量心电图,可从心电信号直接产生Frank矢量心电图,并且可从由Frank矢量心电图计算的矢量幅度的波形中选择所述至少两个波形。
心电图可以是标量心电图,所述标量心电图可以从心电信号合成或直接产生,并且可从标量心电图的波形中选择所述至少两个波形。
标量心电图可以是12导联心电图、Holter心电图、事件心电图、运动心电图和监测心电图中的一种。
可从多次心搏的波形中选择有助于TWA的测量的至少两个波形。
在有助于TWA的测量的波形中包括异常值的情况下,可校正异常值。
附图说明
图1是实施例1的TWA测量心电图仪的方框图。
图2是实施例的TWA测量心电图仪的操作流程图。
图3是示出用于在实施例1的TWA测量心电图仪中导出Frank矢量心电图的过程的流程图。
图4示出图示图3的操作流程图的过程的波形图。
图5是示出用于在实施例1的TWA测量心电图仪中测量TWA的存在的过程的流程图。
图6是图示选择分析段的处理的视图。
图7是图示选择分析段的处理的视图。
图8示出在进行校正离群值的处理之前和之后的波形图。
图9是图示计算周期图的处理的视图。
图10是图示计算周期图的所述处理的视图。
图11是实施例1的TWA测量心电图仪的操作流程图。
图12是实施例1中的TWA测量系统的方框图。
图13是实施例1中的心电图管理系统的方框图。
图14是实施例1中的患者ID管理部的方框图。
图15是示出直到当TWA测量心电图仪在实施例1中的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(存储心电图管理系统的个人系数数据库)的处理的子过程流程图。
图24是实施例2的TWA测量心电图仪的方框图。
图25是实施例2的TWA测量心电图仪的操作流程图。
图26是示出用于测量实施例2的TWA测量心电图仪中的TWA的存在的过程的流程图。
图27示出图示图26的操作流程图的处理的波形图。
图28是实施例3的TWA测量心电图仪的流程图。
图29是示出用于测量实施例3的TWA测量心电图仪中的TWA的存在的过程的流程图。
具体实施方式
在当前公开的主题的TWA测量装置和TWA测量方法中,可从通过各种测量方法中的任何一种测量方法获取的心电图来测量TWA的存在。
在下文中,将分别在实施例1至实施例3中描述当前公开的主题的TWA测量心电图仪和TWA测量方法的各实施例,所述TWA测量心电图仪是TWA测量装置的一个实例。
在实施例1中,从诸如12导联心电图的标量心电图合成Frank矢量心电图,并且从Frank矢量心电图测量TWA的存在。由于从标量心电图合成Frank矢量心电图,所以测量人能够以与获取标量心电图的相同的方式通过简单地将测量电极装接到被检对象来测量TWA的存在。
在下文中,将通过以将12导联心电图用作为标量心电图为例的情况描述实施例1的TWA测量心电图仪和TWA测量方法。尽管在实施例1中将示例性地描述使用12导联心电图的情况,但当前公开的主题还可适用于使用包括合成的导联的至少3导联的标量心电图的情况。至少3导联的标量心电图的示例是12导联心电图、Holter心电图、事件心电图、运动心电图和监测心电图。
(实施例1)
在实施例1的TWA测量心电图仪和TWA测量方法中,由测量电极的心电信号产生12导联心电图。如上所述,在实施例1的TWA测量心电图仪和TWA测量方法中,可从12导联心电图合成Frank矢量心电图,并且可从合成的Frank矢量心电图测量在突然心脏死的预测诊断中有效的TWA的存在。
在实施例中,术语“患者”用作为被检对象的特定示例。然而,患者不仅包括在医院中接受诊断的患者,而且包括诸如进行健康检查的医疗检查中心或诊所或者普通住宅的、除医院以外的机构的用户。
(TWA测量心电图仪的构造)
首先,将描述实施例1的TWA测量心电图仪的构造。图1是实施例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导联心电图导出Frank矢量心电图的变换系数。
构成变换系数存储部114的个人系数数据库113A将先前从特定患者获取并且对于被检对象特定的个人系数存储为变换系数。对每个患者并且以时间序列存储个人系数。在上述情况下,例如,将患者的特定ID“C123”、患者的单独ID“A123”和患者的姓名添加为个人系数的附加信息,并且对于每个患者并以时间序列存储。
当要合成Frank矢量心电图时,通常,使用一个最新的个人系数是足够的。即便如此,由于存在通过观察个人系数的变化或通过利用过去的个人系数和当前的个人系数所合成的Frank矢量心电图的变化能知道症状的发展或恢复状态的情况,所以多个个人系数也以时间序列存储。
构成变换系数存储部114的群组系数数据库113B存储群组系数作为变换系数,该群组系数是从统计有效的人口的非特定数量的人获取的多个变换系数的平均数,以便合成患者的Frank矢量心电图。对于每个性别和年龄群组存储群组系数。
测量电极115是将要附接到患者的体表的电极。当要产生12导联心电图时,或者当要获取对于患者特定的个人系数时,使用测量电极115。
显示部116在显示装置上显示或打印出通过患者信息输入部111输入的患者信息和由TWA测量部118测量的TWA的存在的测量结果。
心电图仪控制部117由测量电极的心电信号产生12导联心电图,并通过利用变换系数存储部114的变换系数从12导联心电图合成Frank矢量心电图。当从患者获取的个人系数存在于个人系数数据库113A中时,心电图仪控制部117在将个人系数用作变换系数的同时合成Frank矢量心电图。这是因为当使用患者的个人系数时,提高TWA的测量的精度。相反,当在个人系数数据库113A中不存在从患者获取的个人系数时,在将存在于群组系数数据库113B中的群组系数用作变换系数的同时合成Frank矢量心电图,以便使得即使当不存在对于患者特定的个人系数时也能够进行TWA的测量。
当在个人系数数据库113A中不存在对于患者特定的个人系数时,心电图仪控制部117使显示部116显示用于提示获取患者的个人系数的消息。这是因为个人系数的获取的提示有助于测量精度的提高。在要获取对于患者特定的个人系数的情况下,心电图仪控制部117从附接到患者的体表的测量电极115的心电信号计算对患者特定的个人系数,并将计算的个人系数存储在个人系数数据库113A中。
在要获取对于患者特定的个人系数的情况下,为了产生12导联心电图,首先,分别将测量电极115附接到在每种测量方法中确定的患者的测量部分,并从测量电极115的心电信号生成12导联心电图。
接下来,分别将测量电极115附接到确定用于产生Frank矢量心电图的患者的测量部分,并且从测量电极115的心电信号产生Frank矢量心电图。测量电极115为了产生12导联心电图而附接到的部分不同于该测量电极115为了产生Frank矢量心电图而附接到的部分。
12导联心电图的产生可与Frank矢量心电图的产生分开进行。然而,在实施例1中,为了提高工作效率,同时进行心电图的产生。
心电图仪控制部117计算个人系数,该个人系数对于患者是特定的并用于使12导联心电图与Frank矢量心电图匹配。在将患者的特定I D、患者的单独ID和患者的姓名添加到个人系数来作为个人系数的附加信息的同时,心电图仪控制部117对每个患者并以时间序列将计算的个人系数存储在个人系数数据库113A中。
心电图仪控制部117总体地控制上述各种操作,以及TWA测量心电图仪110的所有操作。心电图仪控制部117包括用于实现TWA测量心电图仪110的所有操作的程序。后面将详细地描述心电图仪控制部117的操作。
TWA测量部118从由心电图仪控制部117合成的Frank矢量心电图测量TWA的存在。在TWA的存在的测量中,TWA测量部118从由心电图仪控制部117合成的多个Frank矢量心电图的波形中选择能有助于TWA的存在的测量的波形,并从所选择的波形中测量TWA的存在。这是因为当选择能有助于TWA的存在的测量的波形时,提高了TWA的测量的精度。后面将详细地描述TWA测量部118的操作。
通信部119将患者的个人系数发送到心电图管理系统(EMS)140(见图12),并反相地从EMS140接收由心电图仪控制部117请求的患者的个人系数。可通过医院内部网络160进行通信部119中的个人系数的发送/接收。
(TWA测量心电图仪的操作)
接下来,将参考图2的操作流程图描述实施例1的TWA测量心电图仪的操作。
在图2的操作流程图中,由TWA测量心电图仪的操作者(测量人)进行步骤S1到S3的操作,并且由心电图仪控制部117进行步骤S4至S8的操作。步骤S1至S9的操作还对应于实施例1的TWA测量方法的过程。由TWA测量部118进行步骤S9的操作。
<步骤S1>
为了产生12导联心电图,图1所示的TWA测量心电图仪110的操作者将测量电极115附接到患者的体表的预定部分。由于实施例1的TWA测量心电图仪110以由各种测量方法产生的12导联心电图为目标,所以测量电极115分别附接到在所采用的测量方法中确定的患者的测量部分。在要产生其中从4导联I、II、V2和V5合成V1、V3、V4和V6的类型的合成的12导联心电图的情况下,例如,将测量电极115附接到总共六个位置,即为了获取导联I和II的心电信号的四个位置或者右臂与左臂(电极L、R)和右下肢与左下肢(电极LL、RL);以及为了获取两个胸导联(导联V2和导联V5)的心电信号的两个位置或者第四肋间隙的左下胸骨边缘位置和左前腋线与横切第五肋间隙的水平线的交点位置。在要产生标准12导联心电图的情况下,将十个测量电极115附接到总共十个位置,即,用于测量胸导联的六个位置和用于测量四个肢体导联的四个位置。
<步骤S2>
接下来,操作者通过患者信息输入部111输入患者信息。例如,输入作为患者ID的特定ID“C123”和单独ID“A123”,然后输入患者的姓名、患者的年龄和患者的性别。
<步骤S3>
然后,操作者接通TWA测量心电图仪110的测量开关(未示出)。当测量开关接通时,开始TWA的存在的测量。
<步骤S4>
心电图仪控制部117从在步骤S1中附接到患者的六个测量电极115的心电信号产生12导联心电图。
在要产生合成的12导联心电图的情况下,在以下的矩阵公式中代入四个导联、即导联I与II和两个胸导联(导联V2和导联V5)的心电信号(测量导联矢量),并且将该测量导联矢量乘以转换矩阵,从而产生了实际上并不通过测量电极115检测的其余四个胸导联、即导联V1、导联V3、导联V4和导联V6(合成导联矢量)的合成心电图。
其中 是合成导联矢量, 是变换矩阵,并且 是测量导联矢量。
最后,产生12导联心电图。
在要产生标准12导联心电图的情况下,基于由十个测量电极检测到的心电信号,计算标准12导联的六个肢体导联波形(I、II、III、aVR、aVL和aVF)和标准12导联的六个胸导联波形(V1、V2、V3、V4、V5和V6),并最终产生12导联心电图。
<步骤S5至S7>
接下来,心电图仪控制部117从变换系数存储部114获取用于导出Frank矢量心电图的变换系数。
心电图仪控制部117判定在个人系数数据库113A中是否存在与在步骤S2中输入的患者信息一致的个人系数。具体地,判定在个人系数数据库113A中是否存在与从患者信息输入部111输入的特定ID“C123”和单独ID“A123”一致的个人系数。
如果与患者信息一致的个人系数存在于个人系数数据库113A中,则心电图仪控制部117从个人系数数据库113A获取患者的个人系数。对各个患者并且以时间序列存储个人系数。因此,在对于患者存在多个个人系数的情况下,获取最新的个人系数。
相反,如果与患者信息一致的个人系数不在个人系数数据库113A中,则心电图仪控制部117从群组系数数据库113B获取与作为患者信息输入的患者的年龄和性别一致的群组系数。
<步骤S8>
心电图仪控制部117通过利用作为在步骤S5至S7中获取的对于患者特定的个人系数或者群组系数中的任意一者的变换系数,从在步骤S4中产生的12导联心电图合成Frank矢量心电图。
将参考图3的操作流程图具体地描述导出Frank矢量心电图的处理。
<步骤S9>
TWA测量部118从在步骤S8中合成的Frank矢量心电图测量TWA的存在。
如上所述,已知用于测量TWA的存在的多种技术。各技术中的一种是在美国专利号4,802,491中公开的技术,该技术基于频谱分析,并且认为该技术对临床现场最有用。在基于频谱分析的技术中,需要对于其中在负荷施加到被检对象的周期期间HR等于或小于110的128次连续搏动的矢量幅度波形。
因此,当波形受噪声的影响或者不满足HR=110或更低的条件时,不开始搏动计数。考虑到ST与T之间的时间周期取决于HR而变化,所以希望的是尽可能保持HR恒定。然而,在实践中,HR很难保持恒定。即使在HR能保持恒定的情况下,在128次搏动中也可能包括异常值。
在实施例1中,为了排除在基于频谱分析的所述技术中的这种问题,首先,从固定区间中的搏动中选择其中HR为从105到110的128次以上搏动的波形、例如150次搏动的波形。基本上,所选择的150次搏动是连续的。在搏动不连续的情况下,进行选择,以便维持交替。例如,通过将夹在所选择的搏动之间的未选择的搏动的数量设定成偶数或者利用例如所选择的搏动的形状相似性,来进行其中维持交替的选择。
当通过这样的技术提取对于150次搏动的矢量幅度波形时,排除了在施加负荷期间保持HR恒定的必要性,因此能容易地测量TWA的存在。
接下来,从所选择的对于150次搏动的矢量幅度波形中选择对于其中波形之间的相关性等于或大于预定阈值的128次搏动的波形。
此外,在对于128次搏动的矢量幅度波形中包括异常值(离群值)的情况下,通过利用另一值校正该异常值。
作为这样的处理的结果,能尽可能减轻异常值(离群值)的影响,并且提高最后获得的TWA的存在的测量的精度。
将参考图5的操作流程图具体地描述测量TWA的存在的处理。
图3的操作流程图示出了用于导出Frank矢量心电图的过程。流程图是图2中的步骤S8的子过程流程图。
<步骤S8-1>
心电图仪控制部117获取在步骤S4中产生的12导联心电图以及在步骤S6中获取的个人系数或在步骤S7中获取的群组系数中的任何一者作为变换系数。
<步骤S8-2>
心电图仪控制部117将12导联心电图乘以变换系数,以合成Frank矢量心电图的六个导联I、E、C、A、H和M。
在变换系数的相乘之前,心电图仪控制部117进行用于去除所产生的12导联心电图的基线变动的预处理,以便去除12导联心电图的DC分量,以校准基线。同时,进行用于去除12导联心电图的高频分量的预处理,以便去除高频噪声,从而使12导联心电图的波形平滑。
由于以下原理,可通过将12导联心电图乘以变换系数获取Frank矢量心电图。
在12导联心电图L与Frank矢量心电图F之间,存在诸如以下公式1所示的关系。在公式中,α表示变换系数。
Lα=F ......公式1
具体地,在公式1中,L指示12导联心电图的排列(n×8),α指示其中使用个人系数或群组系数的变换系数,并且F指示Frank的矢量心电图的排列(n×6)。
在示例中,在要测量12导联心电图的情况下,使用八个电极,以便测量导联1和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6。因此,12导联心电图构成(n×8)的排列。
如公式1所示,为了从12导联心电图合成Frank矢量心电图的六个导联I、E、C、A、H和M,要求将12导联心电图乘以变换系数α。
为了从12导联心电图合成Frank矢量心电图,需要预先获得变换系数α。正因如此,实施例1的TWA测量心电图仪110设置有两种变换系数,即对于患者特定的个人系数和作为从统计有效的人口的非特定数量的人获取的多个变换系数的平均数的群组系数。
在要获取对患者特定的个人系数的情况下,例如,将测量电极115附接到为产生12导联心电图L(导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6)而确定的患者的特定部分,并且从测量电极115的心电信号产生12导联心电图。接下来,将测量电极115附接到为产生Frank矢量心电图F(导联I、导联E、导联C、导联A、导联H和导联M)而确定的患者的特定部分,并且从测量电极115的心电信号产生Frank矢量心电图。在示例中,在以下的公式2中代入所产生的12导联心电图L(导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6)和所产生的Frank矢量心电图F(导联I、导联E、导联C、导联A、导联H和导联M),并获得对于患者特定的个人系数来作为变换系数α。将参考图11的流程图具体地详细描述用于获取个人系数作为变换系数的处理。
α=(LTL)-1LTF ......公式2
其中LT表示12导联心电图L的转置矩阵,并且(LTL)-1表示(LTL)的逆矩阵。
如公式1所示,心电图仪控制部117将12导联心电图L乘以如此获得的变换系数,以合成Frank矢量心电图F的导联I、导联正、导联C、导联A、导联H和导联M。
<步骤S8-3>
当合成Frank矢量心电图F(导联I、导联E、导联C、导联A、导联H和导联M)时,心电图仪控制部117于是在以下的公式3中代入导联I、导联E、导联C、导联A、导联H和导联M的值,从而获得Frank矢量心电图的导联VX、导联VY和导联VZ。
VX=0.61A+0.171C-0.781I,
VY=0.655F+0.345M-1.0H,
VZ=0.133A+0.736M-0.264I-0.374E-0.231C.
以上的公式由公式3表示。
<步骤S8-4>
最后,心电图仪控制部117将所获得的导联VX、导联VY和导联VZ的值的求平方,将平方值加到一起,并计算平方值的和的平方根,从而计算矢量幅度VM的值。
在上文中,已描述了当要合成Frank矢量心电图时进行的处理。
将参考图4的波形图描述图3的操作流程图的处理。
当如图4的上波形图所示获得导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6的八个心电图波形来作为12导联心电图L时,首先,对心电图波形应用去除基线变动和高频分量的处理,并且使变换系数α乘以所得到的波形。
然后,获得导联I、导联E、导联C、导联A、导联H和导联M的值,并如图4的中间波形图所示获得Frank矢量心电图的导联VX、导联VY和导联VZ的心电图波形。
最后,如图4的下方波形图所示,从导联VX、导联VY和导联VZ的值获得Frank矢量心电图的矢量幅度VM的心电图波形。将从矢量幅度VM的心电图波形获得最终要获得的TWA的存在。
图5的操作流程图示出了用于测量TWA的存在的过程。流程图是图2中的步骤S9的子过程流程图。
<步骤S9-1>
TWA测量部118从心电图仪控制部117获取如上所述获得的Frank矢量心电图的矢量幅度VM的心电图波形。为了提高TWA的存在的测量的精度,TWA测量部118从所获取的矢量幅度VM的心电图波形去除异常值等,并选择用于TWA的存在的测量的波形。
首先,TWA测量部118选择对于128次以上的搏动的波形、例如150次搏动的波形。
如图7所示,在对于150次搏动的波形的选择中,从由心电图仪控制部117获取的矢量幅度VM中仪选择其中HR是从105到110并且其中似乎容易产生TWA的搏动。优选地选择连续搏动。然而,在连续150次搏动的情况下,同样可能存在其中很难测量TWA的存在的嘈杂的搏动的可能性。这影响测量精度。因此,不连续地选择其中似乎容易产生TWA的150次搏动。
在如此选择的对于150次搏动的矢量幅度波形中,存在可能包括在形状方面很大程度地不同的波形。因此,接下来,从对于150次搏动的矢量幅度波形中选择其中波形之间的相关性等于或高于给定阈值的、对于128次搏动的波形。
具体地,首先,获取通过在首先指定的时间周期中将对于150次搏动的矢量幅度波形平均化而获得的平均波形。然后,确定其中各个波形与平均波形之间的相关性等于或高于给定阈值的对于128次搏动的波形。
然后,将如此选择的、对于128次搏动的矢量幅度波形的ST-T段选为分析段。
<步骤S9-2>
在对于128次搏动的分析段中包括离群值(异常值)的情况下,TWA测量部118进行校正离群值的处理。
似乎不利地影响TWA的存在的测量的离群值在被分组成奇数搏动和偶数搏动的同时被检测,并且对检测到的搏动进行校正。
具体地,对于奇数搏动和偶数搏动获得心电图波形的中间值,并且对于奇数搏动和偶数搏动获得标准偏差。计算用于确定离群值并且其中标准偏差用作为参数的阈值,确定关于该阈值的幅度关系,并确定离群值。在其中标准偏差用作为参数的校正函数中代入被确定为离群值的值,并计算校正值。已知确定这样的阈值的各种方法和校正函数。在实施例中,可使用任何已知的方法或函数。
图8示出了在校正离群值的处理的进行之前和之后的波形图。在尚末经历校正处理的原始信号的波形中,如图8的上方波形图所示,在中间部分中存在电平很大程度地下降的部分。相反,在已经历校正处理的校正信号的波形中,如图8的下方波形图所示,在中间部分中不存在这样的电平很大程度地下降的部分。可以看出,校正离群值的处理达到了显著效果。
<步骤S9-3>
在离群值校正之后,如图6所示,TWA测量部118对128次搏动的分析段进行FFT处理,并计算周期图。
图9和10是用于图示周期图的视图。图9和10的上方波形图分别示出了奇数搏动和偶数搏动的平均波形。图9和10的下方波形图示出了在计算周期图之后所获得的波形。
在计算周期图之后所获得的图9的波形示出当频率(循环/搏动)(cycle/beat)为0.5时,矢量幅度的值是小的。这意味着没有TWA。相反,在计算周期图之后获得的图10的波形示出当频率(循环/搏动)为0.5时,矢量幅度的值是大的。这意味着TWA的值是大的,并且TWA存在的可能性是大的。
如上所述,当计算周期图时,能在某种程度上从矢量幅度的值的大小预测TWA的存在。
<步骤S9-4>
接下来,TWA测量部118计算交替(alternance)。在图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中计算的区域的平均数SNB和交替Valt的值,计算以下公式5,从而计算交替率k。
k=(Valt)2/σNB ......公式5
然后,从交替Valt的值和交替率k判定TWA的存在。用于判定TWA的存在的条件是交替Valt>1.9μV和交替率k>3。当满足该判定条件时,判定TWA存在。
如上所述,当分析从12导联心电图获得的Frank矢量心电图的波形形状时,能够判定其中具有不同形状的T波在每次搏动时交替(ABABAB…)出现的TWA的存在。
在不存储患者的个人系数的情况下,接下来,实施例1的TWA测量心电图仪110使显示部116显示用于提示个人系数的获取的消息。此时,TWA测量心电图仪110的操作者获取患者的个人系数,并将该个人系数存储在TWA测量心电图仪110的个人系数数据库113A中。
将参考图11的操作流程图详细描述该情况的过程。在图11的流程图中,由TWA测量心电图仪110的操作者进行步骤S10至S12的操作,并且由心电图仪控制部117进行步骤S13至S15的操作。步骤S13至S15的操作对应于实施例1的TWA测量方法的过程。
<步骤S10>
为了获取12导联心电图和Frank矢量心电图,首先,操作者将测量电极附接到患者的特定部分。
具体地,将测量电极115附接到以便产生12导联心电图(导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6)所确定的患者的特定部分,并附接到以便产生Frank矢量心电图(导联I、导联E、导联C、导联A、导联H和导联M)所确定的患者的特定部分。
<步骤S11>
接下来,操作者通过患者信息输入部111输入患者信息。例如,输入作为患者ID的特定ID“C123”和单独ID“A123”,然后输入患者的姓名和年龄。
<步骤S12>
然后,操作者接通TWA测量心电图仪110的测量开关(未示出)。
<步骤S13>
心电图仪控制部117从附接到患者的测量电极115获取测量部分的心电信号。
<步骤S14>
心电图仪控制部117从测量电极115的心电信号产生12导联心电图,并且还从测量电极115的心电信号产生Frank矢量心电图。在以上的示例中,在以下的公式6中代入所产生的12导联心电图(导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6)和所产生的Frank矢量心电图(导联I、导联E、导联C、导联A、导联H和导联M),并计算对于患者特定的个人系数来作为变换系数α。
α=(LTL)-1LTF ......公式6
其中LT表示12导联心电图L的转置矩阵,并且(LTL)-w表示(LTL)的逆矩阵。
<步骤S15>
心电图仪控制部117将计算的个人系数存储在个人系数数据库113A中。当要存储个人系数时,患者的特定ID“C123”、患者的单独ID“A123”和患者的姓名与年龄被添加为个人系数的附加信息,并以时间序列存储。
在上文中,已描述了实施例1的TWA测量心电图仪和TWA测量方法的操作。
(TWA测量系统的构造)
接下来,将描述实施例1中的TWA测量系统的构造。实施例1中的TWA测量系统的TWA测量心电图仪110、120、130、210、220、230是已在以上的实施例中描述的TWA测量心电图仪。即,心电图仪可从诸如12导联心电图的标量心电图合成Frank矢量心电图,并且可从Frank矢量心电图测量TWA的存在。
图12是实施例1的TWA测量系统的方框图。
TWA测量系统500具有:心电图管理设备100,该心电图管理设备100设置在医院A中;心电图管理设备200,该心电图管理设备200设置在医院B中;患者ID管理部300,该患者ID管理部300管理医院A和B中的患者ID(特定ID和单独ID);以及连接线400。心电图管理设备100、200和患者ID管理部300彼此通过连接线400连接。
在实施例1中,将以心电图管理设备100、200设置在医院中的情况为例。其中设置心电图管理设备100、200的机构不局限于医院。例如,它们可设置在诸如体格检查中心、诊所、学校和进行健康检查的高级保健设施的、医院以外的机构中。
尽管患者ID管理部300设置在心电图管理设备100、200外部,但患者ID管理部300例如可设置在医院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测量心电图仪110、120、130发送与TWA测量相关的信息诸如患者信息、产生的心电图数据和个人系数,以及从该TWA测量心电图仪110、120、130接收该信息,并且管理该信息。
心电图管理系统240通过医院内部网络260连接到TWA测量心电图仪210、220、230。心电图管理系统240向TWA测量心电图仪210、220、230发送与TWA测量相关的信息诸如患者信息、产生的心电图数据和个人系数,以及从该TWA测量心电图仪110、120、130接收该信息,并且管理该信息。
心电图管理系统140、240具有下述功能:检索由TWA测量心电图仪110、120、130、210、220或230通过连接线400请求的患者的个人系数;以及朝请求该患者的个人系数的所述TWA测量心电图仪输出检索的个人系数。因此,将用于实现所述功能的程序安装在心电图管理系统140、240中。
显示装置150从心电图管理系统140接收由TWA测量心电图仪110、120、130产生并与患者的TWA的测量相关的信息,并显示由操作者指定的患者的信息。此外,显示装置显示出TWA测量心电图仪110、120、130获取并与TWA测量相关的信息。
显示装置250从心电图管理系统240接收由TWA测量心电图仪210、220、230产生并与患者的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是实施例1中的心电图管理系统的方框图。在下文中,将描述心电图管理系统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将患者I D发送到患者ID管理部300,并接收变换的患者ID。通过医院内部网络160进行患者信息和个人系数的发送和接收,并且通过连接线400进行患者ID的发送和接收。
图14是实施例1中的患者ID管理部的方框图。
患者ID管理部300包括患者ID数据库310、患者ID变换部320和通信部330。
患者ID数据库310对应地存储患者的特定ID和单独的ID。例如,在特定ID是“C123”、医院A中的单独ID是“A123”和医院B中的单独ID是“B456”的情况下,存储“C123”-“A123”和“C123”-“B456”。
患者ID变换部320在参考患者的特定ID的同时将单独ID变换成相应医院中的单独ID。
例如,在医院A要检索医院B中的个人系数的情况下,在参考“C123”的特定ID的同时访问患者ID数据库310,并且将医院A中的“A123”的单独ID变换成医院B中的“B456”的单独ID。
通信部330接收从心电图管理系统140、240发送的患者的特定ID和单独的ID,并将在患者ID变换部320中变换的单独ID发送到心电图管理系统140、240。
实施例1中的TWA测量系统具有上述构造。
(TWA测量系统的操作)
接下来,将参考图15至23的操作流程图描述实施例1中的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>
如果在个人系数数据库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通过利用获取的群组系数处理由多个测量电极115检测的心电信号,并测量TWA的存在。
<步骤S30>
心电图仪控制部117将在步骤S27或步骤S29中进行的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(no good,否)的信号,并且心电图仪控制部117使显示部116显示消息“要产生个人系数,请进行测量和登记”。
即,如果输入的患者ID和患者姓名与登记在TWA测量心电图仪110或心电图管理系统140中的患者ID和患者姓名一致,则判定所述核对为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一致的个人系数。在实施例1中,从心电图管理系统240检索与TWA测量的结果的患者ID和类型一致的个人系数。在利用患者的特定ID的同时通过利用由患者ID变换部320变换的个人ID进行在心电图管理系统140、240中的个人系数的检索。
<步骤S25-7、S25-3>
如果在心电图管理系统240中存在患者的个人系数(步骤S25-7:是),则心电图仪控制部117从心电图管理系统240的个人系数数据库获取个人系数。
<步骤S25-8>
如果在心电图管理系统240的个人系数数据库中不存在患者的个人系数,则心电图仪控制部117使显示部116显示消息“要产生个人系数,请进行测量和登记”。
图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)而确定的患者的特定部分,并附接到以便产生Frank矢量心电图(导联I、导联E、导联C、导联A、导联H和导联M)而确定的患者的特定部分。心电图仪控制部117从测量电极115测量用于12导联心电图的导联I和II、导联V1、导联V2、导联V3、导联V4、导联V5和导联V6,以及用于Frank矢量心电图的导联I、导联E、导联C、导联A、导联H和导联M。
<步骤S33>
心电图仪控制部117在上述公式2中代入利用测量电极115而测量的导联,以计算患者的个人系数α。
<步骤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)和Frank矢量心电图(导联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中。
如上所述,根据实施例1中的TWA测量心电图仪、TWA测量方法和TWA测量系统,变换系数的使用使程Frank矢量心电图能够从12导联心电图合成。因此,测量人可将测量电极附接到与在要获取12导联心电图的情况下的位置相同的位置。
由于可从合成的12导联心电图测量TWA的存在,所以能够与以在要获取12导联心电图的情况相同的方式测量TWA的存在。此外,使用对于患者特定的个人系数来测量患者的TWA的存在,因此,能获取TWA测量的优化结果,使得能提高测量精度。
已描述了其中对每个患者并以时间序列存储个人系数的实施例1。在TWA测量心电图仪的存储容量小的情况下,可能不以时间序列存储个人系数,而是仅更新地存储最新的个人系数。
实施例1构造成使得在特定患者的个人系数存在的情况下,即使当例如个人系数在两或三年前产生时,也可使用该个人系数。然而,为了能高精确地检查患者的TWA的存在,优选的是通过利用在尽可能接近当前时刻的时刻所获取的个人系数来产生TWA测量的结果。
因此,替代性地,可设定个人系数的可用周期(例如,从获取起一年)。在可用周期期满的情况下,可以以与在个人系数不存在的情况下相类似的方式提示个人系数的获取,并可进行诸如图20的流程图所示的处理。
在实施例1中,进行单独ID的变换,以便在医院之间检索和获取个人系数。根据特定ID进行单独ID的变换。已知用于ID变换的各种技术。用于ID变换的技术不局限于在实施例1中例证的技术,并且可采用各种已知技术中的任何一种技术。例如,称作OpenPIXPDQ的开放软件应用被认为是用于医院之间的ID变换的软件应用。可通过利用这样的开放软件应用进行ID的变换。
(实施例2)
接下来,将描述实施例2。在实施例1中,合成Frank矢量心电图,由合成的Frank矢量心电图计算矢量幅度,并从计算的矢量幅度测量TWA的存在。相反,在实施例2中,由从测量电极的心电信号直接获得的Frank矢量心电图计算矢量幅度,并且从计算的矢量幅度测量TWA的存在。在实施例2中,通过利用独特技术进行TWA的存在的测量。因此,能够以与相关技术的测量精度显著不同的精度测量TWA的存在。
在下文中,将描述实施例2的TWA测量心电图仪和TWA测量方法。
(TWA测量心电图仪的构造)
首先,将描述实施例2的TWA测量心电图仪的构造。图24是实施例2的TWA测量心电图仪的方框图。
TWA测量心电图仪510包括患者信息输入部511、患者信息存储部512、测量电极515、显示部516、心电图仪控制部517和TWA测量部518。
患者信息输入部511用于利用测量人的键操作输入患者信息。具体地,患者信息包含患者ID、患者的姓名、患者的年龄和患者的性别。
患者信息存储部512存储通过患者信息输入部511输入的患者信息。例如,将患者的特定ID“C123”、患者的姓名、患者的年龄和患者的性别存储为患者ID。
测量电极515是要附接到患者的体表的电极。测量电极515附接到为了产生Frank矢量心电图而预先确定的患者的特定位置。
显示部516在显示装置上显示或打印出通过患者信息输入部511输入的患者信息和由TWA测量部518测量的TWA的存在的测量结果。
心电图仪控制部517从测量电极515的心电信号产生Frank矢量心电图,并计算矢量幅度。心电图仪控制部517产生Frank矢量心电图,并总体地控制TWA测量心电图仪510的所有操作。心电图仪控制部517包括用于实现TWA测量心电图仪510的所有操作的程序。稍后将详细地描述心电图仪控制部517的操作。
TWA测量部518从由心电图仪控制部517计算的矢量幅度测量TWA的存在。在TWA的存在的测量中,TWA测量部518从由心电图仪控制部517计算的矢量幅度的多个波形中选择能有助于TWA的存在的测量的波形,并从所选择的波形中测量TWA的存在。这是因为当选择能有助于TWA的存在的测量的波形时,提高了TWA的测量的精度。稍后将详细地描述TWA测量部518的操作。
(TWA测量心电图仪的操作)
接下来,将参考图25的操作流程图描述实施例2的TWA测量心电图仪的操作。
在图25的操作流程图中,步骤S101至S103的操作由TWA测量心电图仪的操作者(测量人)执行,并且步骤S104的操作由心电图仪控制部517执行。步骤S101至S105的操作还对应于实施例2的TWA测量方法的过程。步骤S105的操作由TWA测量部518执行。
<步骤S101>
为了产生Frank矢量心电图,图24所示的TWA测量心电图仪510的操作者将测量电极515附接到患者的体表的预定部分。由于实施例2的TWA测量心电图仪510以Frank矢量心电图为目标,所以测量电极515分别附接到为获取Frank矢量心电图而确定的患者的测量部分。
<步骤S102>
接下来,操作者通过患者信息输入部511输入患者信息。例如,输入作为患者ID的特定ID“C123”和单独ID“A123”,然后输入患者的姓名、患者的年龄和患者的性别。
<步骤S103>
然后,操作者接通TWA测量心电图仪510的测量开关(未示出)。当测量开关接通时,开始TWA的存在的测量。
<步骤S104>
心电图仪控制部517从在步骤S101中附接到患者的测量电极515的心电信号产生Frank矢量心电图。
<步骤S105>
心电图仪控制部517从在步骤S104中产生的Frank矢量心电图计算矢量幅度。TWA测量部518从矢量幅度的波形测量TWA的存在。将参考图26的流程图具体地描述测量TWA的存在的处理。
如上所述,在基于于美国专利号4,802,491中公开的频谱分析的技术中,需要其中HR在当负荷施加到被检对象的周期期间等于或小于110的、对于128次连续搏动的矢量幅度波形。
因此,当波形受噪声的影响或者不满足HR=110或更低的条件时,不开始搏动计数。考虑到ST与T之间的时间周期取决于HR而变化,所以合乎需要的是尽可能保持HR恒定。然而,在实践中,HR很难保持恒定。即使在HR能保持恒定的情况下,在128次搏动中也可能包括异常值。
此外,在实施例2中,为了去除这样的基于频谱分析的技术中的问题,首先,从固定区中的搏动中选择其中HR为从105到110的128次以上的搏动的波形,例如150次搏动的波形。基本上,选择的150次搏动是连续的。在搏动不连续的情况下,进行选择,以便维持交替。例如,通过将夹在所选择的搏动之间的、未选择的搏动的数量设定成偶数或者利用例如所选择的搏动的形状相似性来进行其中维持交替的选择。
当通过这样的技术提取对于150次搏动的矢量幅度时,排除了在施加负荷期间保持HR恒定的必要性,因此能容易地测量TWA的存在。
接下来,从所选择的对于150次搏动的矢量幅度波形中选择其中波形之间的相关性等于或高于预定阈值的对于128次搏动的波形。
此外,在对于128次搏动的矢量幅度波形中包括异常值(离群值)的情况下,通过利用另一值校正该异常值。
作为这样的处理的结果,能尽可能减轻异常值(离群值)的影响,并且提高最终获得的TWA的存在的测量的精度。
图26的操作流程图示出了用于测量TWA的存在的过程。该流程图是图25中的步骤S105的子过程流程图。
<步骤S105-1>
心电图仪控制部517从装接到患者的测量电极515的心电信号计算诸如图27的上图所示的Frank矢量心电图的导联VX、导联VY和导联VZ。
<步骤S105-2>
然后,心电图仪控制部517对计算的导联VX、导联VY和导联VZ的值的进行求平方,将平方值加到一起,并计算平方值的和的平方根,从而计算诸如图27的下图所示的矢量幅度VM的值。
<步骤S105-3>
TWA测量部518从心电图仪控制部517获取如上所述获得的Frank矢量心电图的矢量幅度VM的心电图波形。为了提高TWA的存在的测量的精度,TWA测量部518从矢量幅度VM的获取的心电图波形中去除异常值等,并选择用于TWA的存在的测量的波形。
首先,TWA测量部518选择对于128次以上的搏动的波形,例如150次搏动的波形。
如图7所示,在对于150次搏动的波形的选择中,从心电图仪控制部517所获取的矢量幅度VM中仅选择其中HR是从105到110并且其中似乎容易产生TWA的搏动。优选地选择连续搏动。然而,在连续150次搏动的情况下,同样可能存在其中很难测量TWA的存在的嘈杂的搏动的可能性。这影响测量精度。因此,不连续地选择其中似乎容易产生TWA的150次搏动。
在如此选择的、对于150次搏动的矢量幅度波形中,存在可能包括在形状上很大程度地不同的波形。因此,接下来,从对于150次搏动的矢量幅度波形中选择其中波形之间的相关性等于或高于给定阈值的、对于128次搏动的波形。
具体地,首先,获取通过在首先指定的时间周期中使对于150次搏动的矢量幅度波形平均化而获得的平均波形。然后,确定其中各个波形与平均波形之间的相关性等于或高于给定阈值的、对于128次搏动的波形。
然后,将如此选择的、对于128次搏动的矢量幅度波形的ST-T段选为分析段。
<步骤S105-4>
在对于128次搏动的分析段中包括离群值(异常值)的情况下,TWA测量部518进行校正离群值的处理。
似乎不利地影响TWA的存在的测量的离群值在被分组成奇数搏动和偶数搏动的同时被检测,并且对检测到的搏动进行校正。
具体地,对于奇数搏动和偶数搏动获得心电图波形的中间值,并且对于奇数搏动和偶数搏动获得标准偏差。计算用于确定离群值的并且其中标准偏差被用作为参数的阈值,确定关于阈值的幅度关系,并确定离群值。在其中标准偏差用作为参数的校正函数中代入被确定为离群值的值,并计算校正值。已知确定这样的阈值的各种方法和校正函数。在实施例中,可使用任何已知的方法或函数。
在尚末经历校正处理的原始信号的波形中,如图8的上方波形图所示,在中间部分中存在电平很大程度地下降的部分。相反,在已经历校正处理的校正信号的波形中,如图8的下方波形图所示,在中间部分中不存在这样的电平很大程度地下降的部分。可以看出,校正离群值的处理取得了显著效果。
<步骤S105-5>
在离群值校正之后,如图6所示,TWA测量部518对128次搏动的分析段进行FFT处理,并计算周期图。
图9和10是用于图示周期图的视图。图9和10的上方波形图分别示出了奇数搏动和偶数搏动的平均波形。图9和10的下方波形图示出了在计算周期图之后获得的波形。
在计算周期图之后获得的图9的波形表明当频率(循环/搏动)为0.5时,矢量幅度的值是小的。这意味着没有TWA。相比之下,在计算周期图之后获得的图10的波形表明当频率(循环/搏动)为0.5时,矢量幅度的值是大的。这意味着TWA的值是大的,并且TWA存在的可能性是大的。
如上所述,当计算周期图时,能在某种程度上从矢量幅度的值的大小预测TWA的存在。
<步骤S105-6>
接下来,TWA测量部518计算交替。在图9和10所示并且在计算周期图之后获得的波形图中,将循环/搏动频率从0.44到0.49的区域定义为噪声频带,并获得该区域的平均数SNB和标准偏差σNB。当循环/搏动频率为0.5时所获得的值由S0.5指示,并计算以下公式4,从而计算交替Valt。
Valt=(S0.5-SNB)1/2 ......公式4
<步骤S105-7>
TWA测量部518确定TWA的存在。通过利用在步骤S105-6中计算的区域的平均数SNB和交替Valt的值,计算以下公式5,从而计算交替率k。
k=(Valt)2/σNB ......公式5
然后,从交替Valt的值和交替率k判定TWA的存在。用于判定TWA的存在的条件是交替Valt>1.9μV和交替率k>3。当满足判定条件时,判定TWA存在。
如上所述,当分析Frank矢量心电图的波形形状时,能够判定其中具有不同形状的T波在每次搏动时交替(ABABAB…)出现的TWA的存在。
如上所述,根据实施例2,可从当前使用的Frank矢量心电图精确测量TWA的存在。
替代性地,显示部516可比较地显示通过其中从Frank矢量心电图获取TWA的存在的相关技术所获得的测量结果,以及通过其中根据实施例2获取TWA的存在的技术所获得的测量结果。
(实施例3)
接下来,将描述实施例3。在实施例2中,由从测是电极的心电信号直接获得的Frank矢量心电图计算矢量幅度,并且从所计算的矢量幅度测量TWA的存在。相比之下,在实施例3中,由从测量电极的心电信号合成或直接产生的标量心电图测量TWA的存在。因此,能够从标量心电图测量TWA的存在,虽然在相关技术中并不能进行这样的测量。标量心电图的示例是12导联心电图、Holter心电图、事件心电图、运动心电图和监测心电图。在实施例3中,可从心电图中的任何一种测量TWA的存在。
在下文中,将描述实施例3的TWA测量心电图仪和TWA测量方法。
(TWA测量心电图仪的构造)
首先,将描述实施例3的TWA测量心电图仪的构造。尽管图24是实施例2的TWA测量心电图仪的方框图,实施例3的TWA测量心电图仪的直观构造和实施例2的TWA测量心电图仪的方框图所示的直观构造一样。
TWA测量心电图仪510包括患者信息输入部511、患者信息存储部512、测量电极515、显示部516、心电图仪控制部517和TWA测量部518。
患者信息输入部511用于利用测量人的键操作输入患者信息。具体地,患者信息包含患者ID、患者的姓名、患者的年龄和患者的性别。
患者信息存储部512存储通过患者信息输入部511输入的患者信息。例如,患者的特定ID“C123”、患者的姓名、患者的年龄和患者的性别被存储为患者ID。
测量电极515是要附接到患者的体表的电极。测量电极515附接到为产生标量心电图而预先确定的患者的特定位置。
显示部516在显示装置上显示或打印出通过患者信息输入部511输入的患者信息和由TWA测量部518测量的TWA的存在的测量结果。
心电图仪控制部517从测量电极515的心电信号产生标量心电图。心电图仪控制部517产生标量心电图,并总体地控制TWA测量心电图510的所有操作。心电图仪控制部517包括用于实现TWA测量心电图仪510的所有操作的程序。稍后将详细地描述心电图仪控制部517的操作。
TWA测量部518从由心电图仪控制部517产生的标量心电图测量TWA的存在。在TWA的存在的测量中,TWA测量部518从由心电图仪控制部517产生的多个标量心电图的波形中选择能有助于TWA的存在的测量的波形,并从所选择的波形中测量TWA的存在。这是因为当选择能有助于TWA的存在的测量的波形时,提高了TWA的测量的精度。稍后将详细地描述TWA测量部518的操作。
(TWA测量心电图仪的操作)
接下来,将参考图28的操作流程图描述实施例3的TWA测量心电图仪的操作。
在图28的操作流程图中,步骤S201至S203的操作由TWA测量心电图仪的操作者(测量人)执行,并且步骤S204的操作由心电图仪控制部517执行。步骤S201至S205的操作还对应于实施例2的TWA测量方法的过程。步骤S205的操作由TWA测量部518执行。
<步骤S201>
为了计算标量心电图,图24所示的TWA测量心电图仪510的操作者将测量电极515附接到患者的体表的预定部分。例如,在12导联心电图用作标量心电图的情况下,将测量电极分别附接到为测量12导联心电图而确定的患者的测量部分。在采用Holter心电图、事件心电图、运动心电图或监测心电图的情况下,将测量电极分别附接到为测量心电图而确定的患者的测量部分。
<步骤S202>
接下来,操作者通过患者信息输入部511输入患者信息。例如,输入作为患者ID的特定ID“C123”和单独ID“A123”,然后输入患者的姓名、患者的年龄和患者的性别。
<步骤S203>
然后,操作者接通TWA测量心电图仪510的测量开关(未示出)。当测量开关接通时,开始TWA的存在的测量。
<步骤S204>
心电图仪控制部517从在步骤S201中附接到患者的测量电极515的心电信号产生标量心电图。
<步骤S205>
TWA测量部518从在步骤S204中产生的标量心电图的波形测量TWA的存在。将参考图29的流程图具体地描述测量TWA的存在的处理。
如上所述,在基于于美国专利No.4,802,491中公开的频谱分析的技术中,需要其中HR在当负荷施加到被检对象的周期期间等于或小于110的、对于128次连续搏动的矢量幅度波形。
因此,当波形受噪声的影响或者不满足HR=110或更低的条件时,不开始搏动计数。考虑到ST与T之间的时间周期取决于HR而变化,所以合乎需要的是尽可能保持HR恒定。然而,在实践中,HR很难保持恒定。即使在HR能保持恒定的情况下,在128次搏动中也可能包括异常值。
而且,在实施例3中,为了去除这样的基于频谱分析的技术中的问题,首先,从固定区中的搏动选择其中HR是从105到110的、例如150次搏动的128次以上搏动的波形。基本上,选择的150次搏动是连续的。在搏动不连续的情况下,进行选择,以便维持交替。例如,通过将夹在所选择的搏动之间的、未选择的搏动的数量设定成偶数或者利用例如所选择的搏动的形状相似性来进行其中维持交替的选择。
当通过这样的技术提取对于150次搏动的标量心电图波形时,排除了在施加负荷期间保持HR恒定的必要性,因此能容易地测量TWA的存在。
接下来,从所选择的、对于150次搏动的标量心电图波形中选择其中波形之间的相关性等于或高于预定阈值的、对于128次搏动的波形。
此外,在对于128次搏动的标量心电图波形中包括异常值(离群值)的情况下,通过利用另一值校正该异常值。
由于这样的处理,能尽可能减轻异常值(离群值)的影响,并且提高最后获得的TWA的存在的测量的精度。
图29的操作流程图示出了用于测量TWA的存在的过程。该流程图是图28中的步骤S205的子过程流程图。
<步骤S205-1>
心电图仪控制部517从附接到患者的测量电极515的心电信号产生通常从测量电极515获得的标量心电图。TWA测量部518从心电图仪控制部517获取如此产生的标量心电图波形。为了提高TWA的存在的测量的精度,TWA测量部518从获取的标量心电图波形中去除异常值等,并选择用于TWA的存在的测量的波形。
首先,TWA测量部518选择对于128次以上的搏动的波形,例如对于150次搏动的波形。如图7所示,在对于150次搏动的波形的选择中,从心电图仪控制部517获取的标量心电图波形中仪选择其中HR是从105到110并且其中似乎容易产生TWA的搏动。优选地选择连续搏动。然而,在连续150次搏动的情况下,同样可能存在其中很难测量TWA的存在的嘈杂的搏动的可能性。这影响测量精度。因此,不连续地选择其中似乎容易产生TWA的150次搏动。
在如此选择的、对于150次搏动的标量心电图波形中,存在可能包括在形状上很大程度地不同的波形。因此,接下来,从对于150次搏动的矢量幅度波形中选择其中波形之间的相关性等于或高于给定阈值的、对于128次搏动的波形。
具体地,首先,获取通过在首先指定的时间周期中使对于150次搏动的波形平均化而获得的平均波形。然后,判定其中各个波形与平均波形之间的相关性等于或高于给定阈值的、对于128次搏动的波形。
然后,将如此选择的、对于128次搏动的标量心电图波形的ST-T段选为分析段。
<步骤S205-2>
在对于128次搏动的分析段中包括离群值(异常值)的情况下,TWA测量部518进行校正离群值的过程。
似乎不利地影响TWA的存在的测量的离群值在被分组成技术搏动和偶数搏动的同时被检测,并且对检测到的搏动进行校正。
具体地,对于奇数搏动和偶数搏动获得心电图波形的中间值,并且对于奇数搏动和偶数搏动获得标准偏差。计算用于确定离群值的并且其中标准偏差被用作参数的阈值,确定关于阈值的幅度关系,并确定离群值。在其中标准偏差用作参数的校正函数中代入被确定为离群值的值,并计算校正值。已知确定这样的阈值的各种方法和校正函数。在实施例中,可使用任何已知的方法或函数。
在尚末经历校正处理的原始信号的波形中,如图8的上方波形图所示,在中间部分中存在电平很大程度地下降的部分。相反,在已经历校正处理的校正信号的波形中,如图8的下方波形图所示,在中间部分中不存在这样的电平很大程度地下降的部分。可以看出,校正离群值的处理取得了显著效果。
<步骤S205-3>
在离群值校正之后,如图6所示,TWA测量部518对128次搏动的分析段进行FFT处理,并计算周期图。
图9和10是用于图示周期图的视图。图9和10的上方波形图分别示出了奇数搏动和偶数搏动的平均波形。图9和10的下方波形图示出了在计算周期图之后获得的波形。
在计算周期图之后所获得的图9的波形示出当频率(循环/搏动)为0.5时,矢量幅度的值是小的。这意味着没有TWA。相反,在计算周期图之后所获得的图10的波形示出当频率(循环/搏动)为0.5时,矢量幅度的值是大的。这意味着TWA的值是大的,并且TWA存在的可能性是大的。
如上所述,当计算周期图时,能在某种程度上从矢量幅度的值的大小预测TWA的存在。
<步骤S205-4>
接下来,TWA测量部518计算交替。在图9和10所示并且在计算周期图之后获得的波形图中,将循环/搏动频率从0.44到0.49的区域定义为噪声频带,并获得该区域的平均数SNB和标准偏差σNB。当循环/搏动频率为0.5时获得的值由S0.5表示,并计算以下公式4,从而计算交替Valt。
Valt=(S0.5-SNB)1/2 ......公式4
<步骤S205-5>
TWA测量部518判定TWA的存在。通过利用在步骤S205-4中计算的区域的平均数SNB和交替Valt的值,计算以下公式5,从而计算交替率k。
k=(Valt)2/σNB ......公式5
然后,从交替Valt的值和交替率k判定TWA的存在。用于判定TWA的存在的条件是交替Valt>1.9μV和交替率k>3。当满足所述判定条件时,判定TWA存在。
如上所述,当分析标量心电图的波形形状时,能够判定其中具有不同形状的T波在每次搏动时交替(ABABAB…)出现的TWA的存在。
如上所述,根据实施例3,可从当前使用的标量心电图精确测量TWA的存在。
如上所述,根据当前公开的主题的一方面,通过利用从相关技术中使用的Frank矢量心电图和标量心电图中的任何一种所产生的心电图,能精确地测量TWA的存在。
根据当前公开的主题的一方面,心电图仪控制部可通过利用各种测量方法中的任何一种产生心电图。要附接到被检对象的测量电极的数量和位置取决于用于获取心电图的测量方法而不同。
例如,根据Frank矢量心电图或者通常的标量心电图,即,标准的12导联心电图、合成的12导联心电图、Holter心电图、事件心电图、运动心电图或监测心电图,要附接到被检对象的测量电极的数量和位置是不同的。TWA测量部基于由心电图仪控制部产生的矢量幅度或心电图测量TWA的存在。
根据当前公开的主题的一方面,可从由各种测量方法中的任何一种获取的心电图测量TWA的存在。此外,能够测量TWA的存在,而不受条件限制,并且因此能减轻在测量人和被检对象上的负担。
Claims (18)
1.一种TWA测量装置,包括:
心电图仪控制部,该心电图仪控制部构造成从被检对象的心电信号产生心电图;以及
TWA测量部,该TWA测量部构造成从所述心电图中选择有助于TWA的测量的至少两个波形,并且该TWA测量部构造成通过利用所选择的所述至少两个波形测量TWA的存在。
2.根据权利要求1所述的TWA测量装置,还包括变换系数存储部,该变换系数存储部构造成存储变换系数,其中,
所述心电图是Frank矢量心电图,
所述心电图仪控制部通过利用所述变换系数从所述心电信号合成所述Frank矢量心电图,以及
所述TWA测量部从由所述Frank矢量心电图计算的矢量幅度的波形中选择所述至少两个波形。
3.根据权利要求2所述的TWA测量装置,其中,
所述变换系数存储部包括:
个人系数数据库,该个人系数数据库构造成将对于所述被检对象特定的并且从所述被检对象获取的个人系数存储为所述变换系数;以及
群组系数数据库,该群组系数数据库构造成将群组系数存储为所述变换系数以便合成所述被检对象的所述Frank矢量心电图,该群组系数是从统计有效的人口的非特定数量的人获取的多个变换系数的平均数。
4.根据权利要求3所述的FWA测量装置,其中,
当从所述被检对象获取的所述个人系数存在于所述个人系数数据库中时,所述心电图仪控制部通过利用该个人系数作为所述变换系数来合成所述Frank矢量心电图,以及
当所述个人系数不存在于所述个人系数数据库中时,所述心电图仪控制部通过利用在所述群组系数数据库中存在的所述群组系数作为所述变换系数的来合成所述Frank矢量心电图。
5.根据权利要求3所述的TWA测量装置,其中,
当所述被检对象的所述个人系数不存于所述个人系数数据库中时,所述心电图仪控制部使显示部显示用于提示所述个人系数的获取的消息。
6.根据权利要求3所述的TWA测量装置,其中,
所述心电图仪控制部从所述心电信号计算所述被检对象的所述个人系数,并将所述计算的个人系数存储在所述个人系数数据库中。
7.根据权利要求1所述的TWA测量装置,其中,
所述心电图是Frank矢量心电图,
所述心电图仪控制部从所述心电信号直接产生所述Frank矢量心电图,以及
所述TWA测量部从由所述Frank矢量心电图计算的矢量幅度的波形中选择所述至少两个波形。
8.根据权利要求1所述的TWA测量装置,其中,
所述心电图是标量心电图,
所述心电图仪控制部从所述心电信号合成或直接产生所述标量心电图,以及
所述TWA测量部从所述标量心电图的波形中选择所述至少两个波形。
9.根据权利要求8所述的TWA测量装置,其中,
所述标量心电图是12导联心电图、Holter心电图、事件心电图、运动心电图和监测心电图中的一种。
10.根据权利要求1所述的FWA测量装置,其中,
所述TWA测量部从多次心搏的波形中选择有助于所述TWA的测量的所述至少两个波形。
11.根据权利要求10所述的TWA测量装置,其中,
在有助于所述TWA的测量的所述至少两个波形中包括异常值的情况下,所述TWA测量部校正所述异常值。
12.一种TWA测量方法,包括:
从被检对象的心电信号产生心电图;
从所述心电图中选择有助于TWA的测量的至少两个波形;以及
通过利用所选择的所述至少两个波形测量TWA的存在。
13.根据权利要求12所述的TWA测量方法,其中,
所述心电图是Frank矢量心电图,
通过利用变换系数从所述心电信号合成所述Frank矢量心电图,以及
从由所述Frank矢量心电图计算的矢量幅度的波形中选择所述至少两个波形。
14.根据权利要求12所述的TWA测量方法,其中,
所述心电图是Frank矢量心电图,
从所述心电信号直接产生所述Frank矢量心电图,以及
从由所述Frank矢量心电图计算的矢量幅度的波形中选择所述至少两个波形。
15.根据权利要求12所述的TWA测量方法,其中,
所述心电图是标量心电图,
所述标量心电图从所述心电信号合成或直接产生,以及
从所述标量心电图的波形中选择所述所述至少两个波形。
16.根据权利要求15所述的TWA测量方法,其中,
所述标量心电图是12导联心电图、Holter心电图、事件心电图、运动心电图和监测心电图中的一种。
17.根据权利要求12所述的TWA测量方法,其中,
从多次心搏的波形中选择有助于所述TWA的测量的所述至少两个波形。
18.根据权利要求17所述的TWA测量方法,其中,
在有助于所述TWA的测量的所述至少两个波形中包括异常值的情况下,校正所述异常值。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012049525A JP5926074B2 (ja) | 2012-03-06 | 2012-03-06 | Twa測定装置及びtwa測定装置の作動方法 |
JP2012-049525 | 2012-03-06 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103300847A true CN103300847A (zh) | 2013-09-18 |
Family
ID=47827034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013100854515A Pending CN103300847A (zh) | 2012-03-06 | 2013-03-06 | Twa测量装置和twa测量方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8843193B2 (zh) |
EP (1) | EP2636363A3 (zh) |
JP (1) | JP5926074B2 (zh) |
CN (1) | CN103300847A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106073760A (zh) * | 2016-04-29 | 2016-11-09 | 北京蓬阳丰业医疗设备有限公司 | 基于威尔逊导联的18导联动态心电图分析方法及系统 |
CN107072582A (zh) * | 2014-10-22 | 2017-08-18 | 株式会社心电技术研究所 | 脑电波测量装置、脑电波测量方法以及脑电波测量系统 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6173905B2 (ja) * | 2013-12-16 | 2017-08-02 | 大名 魏 | Twa測定装置及びtwa測定装置の作動方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1520777A (zh) * | 2003-02-14 | 2004-08-18 | GEҽҩϵͳ��Ϣ�Ƽ���˾ | 用于改善t波交替的测量的方法和系统 |
CN1682654A (zh) * | 2004-04-15 | 2005-10-19 | Ge医药系统信息科技公司 | 用于确定心电图信号的交替数据的方法和装置 |
US20060116596A1 (en) * | 2004-12-01 | 2006-06-01 | Xiaohong Zhou | Method and apparatus for detection and monitoring of T-wave alternans |
US20070260151A1 (en) * | 2006-05-03 | 2007-11-08 | Clifford Gari D | Method and device for filtering, segmenting, compressing and classifying oscillatory signals |
US20080015452A1 (en) * | 2006-06-30 | 2008-01-17 | Ricci Carlos A | Method of processing electrocardiogram waveform |
US20090306526A1 (en) * | 2008-06-09 | 2009-12-10 | David Cuesta Frau | Method and apparatus for automatic analysis of t-wave alternans |
CN102160784A (zh) * | 2010-02-16 | 2011-08-24 | 通用电气公司 | 用于患者评估的方法和系统 |
CN102342830A (zh) * | 2010-07-30 | 2012-02-08 | 日本光电工业株式会社 | 推演心电图生成系统以及推演心电图生成方法 |
CN102727193A (zh) * | 2011-03-29 | 2012-10-17 | 魏大名 | Twa测量心电图仪、twa测量方法和twa测量系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4802491A (en) | 1986-07-30 | 1989-02-07 | Massachusetts Institute Of Technology | Method and apparatus for assessing myocardial electrical stability |
EP0746229B1 (en) * | 1994-01-26 | 2003-09-10 | Cambridge Heart, Inc. | Apparatus for measuring and assessing cardiac electrical stability |
US5935082A (en) | 1995-01-26 | 1999-08-10 | Cambridge Heart, Inc. | Assessing cardiac electrical stability |
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 |
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 |
-
2012
- 2012-03-06 JP JP2012049525A patent/JP5926074B2/ja active Active
-
2013
- 2013-03-05 US US13/784,894 patent/US8843193B2/en active Active
- 2013-03-06 EP EP13158014.4A patent/EP2636363A3/en not_active Withdrawn
- 2013-03-06 CN CN2013100854515A patent/CN103300847A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1520777A (zh) * | 2003-02-14 | 2004-08-18 | GEҽҩϵͳ��Ϣ�Ƽ���˾ | 用于改善t波交替的测量的方法和系统 |
CN1682654A (zh) * | 2004-04-15 | 2005-10-19 | Ge医药系统信息科技公司 | 用于确定心电图信号的交替数据的方法和装置 |
US20060116596A1 (en) * | 2004-12-01 | 2006-06-01 | Xiaohong Zhou | Method and apparatus for detection and monitoring of T-wave alternans |
US20070260151A1 (en) * | 2006-05-03 | 2007-11-08 | Clifford Gari D | Method and device for filtering, segmenting, compressing and classifying oscillatory signals |
US20080015452A1 (en) * | 2006-06-30 | 2008-01-17 | Ricci Carlos A | Method of processing electrocardiogram waveform |
US20090306526A1 (en) * | 2008-06-09 | 2009-12-10 | David Cuesta Frau | Method and apparatus for automatic analysis of t-wave alternans |
CN102160784A (zh) * | 2010-02-16 | 2011-08-24 | 通用电气公司 | 用于患者评估的方法和系统 |
CN102342830A (zh) * | 2010-07-30 | 2012-02-08 | 日本光电工业株式会社 | 推演心电图生成系统以及推演心电图生成方法 |
CN102727193A (zh) * | 2011-03-29 | 2012-10-17 | 魏大名 | Twa测量心电图仪、twa测量方法和twa测量系统 |
Non-Patent Citations (1)
Title |
---|
MAN S等: "Individually improved VCG synthesis", 《COMPUTERS IN CARDIOLOGY》, vol. 36, 31 December 2009 (2009-12-31), XP031656322 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107072582A (zh) * | 2014-10-22 | 2017-08-18 | 株式会社心电技术研究所 | 脑电波测量装置、脑电波测量方法以及脑电波测量系统 |
US10835144B2 (en) | 2014-10-22 | 2020-11-17 | Ekg Technology Lab, Inc. | Electroencephalograph, brain wave measurement method, brain wave measurement system |
CN107072582B (zh) * | 2014-10-22 | 2020-11-27 | 株式会社心电技术研究所 | 脑电波测量装置、脑电波测量方法以及脑电波测量系统 |
CN106073760A (zh) * | 2016-04-29 | 2016-11-09 | 北京蓬阳丰业医疗设备有限公司 | 基于威尔逊导联的18导联动态心电图分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
EP2636363A2 (en) | 2013-09-11 |
US20130237871A1 (en) | 2013-09-12 |
US8843193B2 (en) | 2014-09-23 |
JP2013183801A (ja) | 2013-09-19 |
EP2636363A3 (en) | 2014-05-14 |
JP5926074B2 (ja) | 2016-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9131864B2 (en) | System and method for evaluating an electrophysiological signal | |
Loring et al. | A detailed guide for quantification of myocardial scar with the Selvester QRS score in the presence of electrocardiogram confounders | |
CN101057781B (zh) | 用于分析和编辑ecg形态和时间序列的方法和设备 | |
KR101892176B1 (ko) | 심전도 디스플레이 방법 및 이를 위한 사용자 인터페이스 장치 | |
Hubley-Kozey et al. | Spatial features in body-surface potential maps can identify patients with a history of sustained ventricular tachycardia | |
CN109411041A (zh) | 心电信息处理方法和心电工作站系统 | |
CN102727193B (zh) | Twa测量心电图仪、twa测量方法和twa测量系统 | |
WO2008064682A2 (en) | A system and a method for dimensionality reduction of ecg signals for ecg analysis and presentation | |
CN103300847A (zh) | Twa测量装置和twa测量方法 | |
Sammani et al. | Life-threatening ventricular arrhythmia prediction in patients with dilated cardiomyopathy using explainable electrocardiogram-based deep neural networks | |
Huebner et al. | Electrocardiologic and related methods of non-invasive detection and risk stratification in myocardial ischemia: state of the art and perspectives | |
Shi et al. | Correlation between vectorcardiographic measures and cardiac magnetic resonance imaging of the left ventricle in an implantable cardioverter defibrillator population | |
US20190175118A1 (en) | Ecg training and skill enhancement | |
CN103356184B (zh) | Twa测量装置和twa测量方法 | |
Waseem et al. | Using evolutionary algorithms for ECG Arrhythmia detection and classification | |
JP2013183801A5 (zh) | ||
Bashir et al. | Highlighting the current issues with pride suggestions for improving the performance of real time cardiac health monitoring | |
Jaber et al. | Association between T wave morphology parameters and abnormal cardiac SPECT imaging | |
Elmansouri et al. | Developing a real time electrocardiogram system using virtual bio-instrumentation | |
Soares et al. | Studying Arrhythmic Risk with In-Silico Programmed Ventricular Stimulation and Patient-Specific Computational Models | |
Chouvarda et al. | Time-domain analysis of the electrocardiogram | |
US11568991B1 (en) | Medical diagnostic tool with neural model trained through machine learning for predicting coronary disease from ECG signals | |
Bouwmeester et al. | The relationship between initial ST-segment deviation and final QRS complex changes related to the posterolateral wall in acute inferior myocardial infarction | |
Pattacini | Custom attention and data division strategies in transformer architectures for ECG beat classification | |
Bergquist et al. | Body Surface Potential Mapping: Contemporary Applications and Future Perspectives. Hearts 2021, 2, 514–542 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20130918 |
|
RJ01 | Rejection of invention patent application after publication |