CN103892826A - 基于心电图的心率分析方法及设备 - Google Patents

基于心电图的心率分析方法及设备 Download PDF

Info

Publication number
CN103892826A
CN103892826A CN201410146200.8A CN201410146200A CN103892826A CN 103892826 A CN103892826 A CN 103892826A CN 201410146200 A CN201410146200 A CN 201410146200A CN 103892826 A CN103892826 A CN 103892826A
Authority
CN
China
Prior art keywords
heart rate
electrocardiogram
vector
time
scatterplot
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.)
Granted
Application number
CN201410146200.8A
Other languages
English (en)
Other versions
CN103892826B (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.)
Medex Tech-Trade Corp
Original Assignee
Medex Tech-Trade 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 Medex Tech-Trade Corp filed Critical Medex Tech-Trade Corp
Priority to CN201410146200.8A priority Critical patent/CN103892826B/zh
Publication of CN103892826A publication Critical patent/CN103892826A/zh
Application granted granted Critical
Publication of CN103892826B publication Critical patent/CN103892826B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明涉及一种基于心电图的心率分析方法及设备,其中,该方法包括:根据心电图绘制心率散点图;分别连接所述心率散点图中时间上相邻的两个散点,以得到多个矢量;以及按所述矢量出现的时间顺序绘制所述矢量的栅状图,其中所述矢量出现的时间顺序与用于构成所述矢量的散点出现的时间顺序一致。根据本发明实施例的基于心电图的心率分析方法及设备,通过分别连接心率散点图中时间上相邻的两个散点以得到多个矢量,并按所述矢量出现的时间顺序绘制所述矢量的栅状图,可以在栅状图中清楚的反映心率与时间点的对应关系,这样,医生根据矢量的栅状图在进行心律失常等疾病诊断时,可以清楚、直观的看到疾病所发生的时间点。

Description

基于心电图的心率分析方法及设备
技术领域
本发明涉及心电信号分析技术,尤其涉及一种基于心电图的心率分析方法及设备。
背景技术
目前,已知有以原始心电图(下称心电图)为基础,绘制出心率散点图,以辅助医生进行如心律失常、心肌梗死、心肌缺血等疾病的诊断。
例如,可以采用心电图中的RR间期数据来绘制心率散点图。图1示出了截取的部分心电图的示意图。如图1所示,心电图的波形通常可以包括P波、QRS波群(Q波、R波、S波)、T波等波形。其中,RR间期表示心电图中两个相邻QRS波群中R波之间的时间差,如图1示出的R1R2(x1)、R2R3(x2)、R3R4(x3)。根据心电图中的RR间期数据来绘制心率散点图的具体方法主要可以包括:首先,假设心电图中有M段连续的RR间期x1,x2,x3,x4,……,xM,M为正整数,提取这M段连续的RR间期;然后,建立如图2所示的X轴(横轴)、Y轴(纵轴)均为RR间期的二维平面直角坐标系;最后,在上述二维平面直角坐标系中,按RR间期的时间先后顺序,连续绘制出横坐标为xi、纵坐标为xi+1的散点,1<i<M且i为整数,从而构成如图3所示的心率散点图。如此绘制出的心率散点图中的散点依次为(x1,x2)、(x2,x3)、(x3,x4)、……、(xM-1,xM),其中(xi-1,xi)与(xi,xi+1)可称为相邻的散点。
一般来说,正常人的RR间期是均匀的,这使得在基于心电图的RR间期绘制出的心率散点图上,散点基本上分布在通过坐标原点且斜率为1的直线附近,也即如图3所示,这些散点在X轴上的坐标与在Y轴上的坐标大致相等。
然而,如果某人存在某种心脏疾病,则心率散点图上的散点可能出现其它的分布。例如:
室上性早搏可能出现如图4所示的心率散点图。
又如,室性早搏可能出现如图5所示的心率散点图。
又如,房颤可能出现如图6所示的心率散点图:
图3、图4、图5以及图6所示的心率散点图虽然在一定程度上能给医生对疾病的诊断提供帮助,但明显缺失了时间信息,这有可能对医生的诊断造成一定的障碍。
发明内容
技术问题
有鉴于此,本发明要解决的技术问题是,如何根据心电图获得结合时间的心率信息。
解决方案
为了解决上述技术问题,根据本发明一实施例,在第一方面,提供一种基于心电图的心率分析方法,包括:
根据心电图绘制心率散点图;
分别连接所述心率散点图中时间上相邻的两个散点,以得到多个矢量;以及
按所述矢量出现的时间顺序绘制所述矢量的栅状图,其中所述矢量出现的时间顺序与用于构成所述矢量的散点出现的时间顺序一致。
结合第一方面,在一种可能的实现方式中,按所述矢量出现的时间顺序绘制所述矢量的栅状图,包括:
计算各所述矢量的大小以及角度;
按所述矢量的出现顺序绘制所述栅状图,其中,所述栅状图的横轴表示所述矢量出现的时间顺序,所述栅状图的纵轴表示所述矢量的大小以及角度中任一。
结合第一方面,在一种可能的实现方式中,根据心电图绘制心率散点图,包括:
对预定时间段内的心电图上的波形进行采样,获得多个心动周期数据;
将时间上连续的N个心动周期数据分别作为一个散点的一维数据,绘制N维的所述心率散点图,其中N为大于或等于2的整数。
结合第一方面,在一种可能的实现方式中,绘制N维的所述心率散点图,包括:
按所述散点出现的时间顺序,对所述散点的绘制颜色进行编码;
以编码后的颜色,绘制各所述散点。
结合第一方面,在一种可能的实现方式中,按所述矢量的出现顺序绘制所述矢量的栅状图,还包括:
在所述栅状图的横轴上,标识与所述矢量相关联的心动周期数据出现的时间,其中与所述矢量相关联的心动周期数据表示与用于构成所述矢量的任一散点相关联的心动周期数据中任一。
结合第一方面,在一种可能的实现方式中,所述心动周期数据包括以下数据中任一:
RR间期,表示心电图中两个相邻QRS波群中R波之间的时间差;
PP间期,表示心电图中两个相邻P波之间的时间差;
QRS波宽度,表示心电图中一个QRS波群中Q波起始到S波结束的时间差;
PR间期,表示心电图中P波起始到QRS波起始的时间差;
R波幅值,表示心电图中一次R波高度的最大值;
T波幅值,表示心电图中一次T波高度的最大值;
QT间期,表示心电图中QRS波起始到T波结束的时间差。
为了解决上述技术问题,根据本发明另一实施例,在第二方面,提供一种基于心电图的心率分析设备,包括:
散点图绘制模块,用于根据心电图绘制心率散点图;
矢量生成模块,与所述散点图绘制模块连接,用于分别连接所述心率散点图中时间上相邻的两个散点,以得到多个矢量;
栅状图绘制模块,与所述矢量生成模块连接,用于按所述矢量出现的时间顺序绘制所述矢量的栅状图,其中所述矢量出现的时间顺序与用于构成所述矢量的散点出现的时间顺序一致。
结合第二方面,在一种可能的实现方式中,所述栅状图绘制模块用于:
计算各所述矢量的大小以及角度;
按所述矢量的出现顺序绘制所述栅状图,其中,所述栅状图的横轴表示所述矢量出现的时间顺序,所述栅状图的纵轴表示所述矢量的大小以及角度中任一。
结合第二方面,在一种可能的实现方式中,所述散点图绘制模块用于:
对预定时间段内的心电图上的波形进行采样,获得多个心动周期数据;
将时间上连续的N个心动周期数据分别作为一个散点的一维数据,绘制N维的所述心率散点图,其中N为大于或等于2的整数。
结合第二方面,在一种可能的实现方式中,所述散点图绘制模块还用于:
按所述散点出现的时间顺序,对所述散点的绘制颜色进行编码;
以编码后的颜色,绘制各所述散点。
结合第二方面,在一种可能的实现方式中,所述栅状图绘制模块还用于:
在所述栅状图的横轴上,标识与所述矢量相关联的心动周期数据出现的时间,其中与所述矢量相关联的心动周期数据表示与用于构成所述矢量的任一散点相关联的心动周期数据中任一。
结合第二方面,在一种可能的实现方式中,所述心动周期数据包括以下数据中任一:
RR间期,表示心电图中两个相邻QRS波群中R波之间的时间差;
PP间期,表示心电图中两个相邻P波之间的时间差;
QRS波宽度,表示心电图中一个QRS波群中Q波起始到S波结束的时间差;
PR间期,表示心电图中P波起始到QRS波起始的时间差;
R波幅值,表示心电图中一次R波高度的最大值;
T波幅值,表示心电图中一次T波高度的最大值;
QT间期,表示心电图中QRS波起始到T波结束的时间差。
有益效果
本发明实施例的基于心电图的心率分析方法及设备,通过分别连接心率散点图中时间上相邻的两个散点以得到多个矢量,并按所述矢量出现的时间顺序绘制所述矢量的栅状图,可以在栅状图中清楚的反映心率与时间点的对应关系,这样,医生根据矢量的栅状图在进行心律失常等疾病诊断时,可以清楚、直观的看到疾病所发生的时间点。
根据下面参考附图对示例性实施例的详细说明,本发明的其它特征及方面将变得清楚。
附图说明
包含在说明书中并且构成说明书的一部分的附图与说明书一起示出了本发明的示例性实施例、特征和方面,并且用于解释本发明的原理。
图1示出截取的部分心电图的示意图;
图2示出以X轴、Y轴均为RR间期的二维平面直角坐标系绘制现有心率散点图的示意图;
图3示出绘制出的正常状态下的心率散点图;
图4示出绘制出的室上性早搏下的心率散点图;
图5示出绘制出的室性早搏下的心率散点图;
图6示出绘制出的房颤下的心率散点图;
图7示出根据本发明实施例的一种基于心电图的心率分析方法的流程图;
图8示出根据本发明实施例的另一种基于心电图的心率分析方法的流程图;
图9a示出根据本发明实施例的基于心电图的心率分析方法的矢量大小和方向的示意图;
图9b示出根据本发明实施例的基于矢量大小绘制栅状图的方法的示意图;
图9c示出根据本发明实施例的基于矢量方向绘制栅状图的方法的示意图;
图10a、10b分别示出根据本发明实施例绘制出的正常状态下的长度栅状图、角度栅状图;
图11a、11b分别示出根据本发明实施例绘制出的室上性早搏的长度栅状图、角度栅状图;
图12a、12b分别示出根据本发明实施例绘制出的室性早搏的长度栅状图、角度栅状图;
图13a、13b分别示出根据本发明实施例绘制出的房颤的长度栅状图、角度栅状图;
图14示出根据本发明一实施例的基于心电图的心率分析设备的结构框图。
具体实施方式
以下将参考附图详细说明本发明的各种示例性实施例、特征和方面。附图中相同的附图标记表示功能相同或相似的元件。尽管在附图中示出了实施例的各种方面,但是除非特别指出,不必按比例绘制附图。
在这里专用的词“示例性”意为“用作例子、实施例或说明性”。这里作为“示例性”所说明的任何实施例不必解释为优于或好于其它实施例。
另外,为了更好的说明本发明,在下文的具体实施方式中给出了众多的具体细节。本领域技术人员应当理解,没有某些具体细节,本发明同样可以实施。在一些实例中,对于本领域技术人员熟知的方法、手段、元件和电路未作详细描述,以便于凸显本发明的主旨。
如背景技术和发明内容所述,已知有以心电图为基础绘制相应的心率散点图,以辅助医生进行疾病的诊断。然而,由于心率散点图上缺失时间信息,这有可能对医生的诊断造成一定的障碍。对此,本发明人独创性地提出,可通过将心率散点图中时间上相邻的散点连接起来构成心率的矢量,并求出该矢量的大小(还可称为“长度”)以及方向(还可称为“角度”),然后分别以矢量的大小或者方向为纵轴并且以与该矢量对应的散点出现的时间(顺序)为横轴绘制出一种栅状图,以能够体现心率变化与时间的对应关系,即能够直观的体现某些疾病、如心律失常具体发生的时间,从而更有助于医生进行疾病诊断,具体方案阐述如下。
图7示出根据本发明实施例的一种基于心电图的心率分析方法的流程图。如图7所示,该方法主要可以包括以下步骤:
步骤S110、根据心电图绘制心率散点图;
步骤S120、分别连接所述心率散点图中时间上相邻的两个散点,以得到多个矢量;
步骤S130、按所述矢量出现的时间顺序绘制所述矢量的栅状图,其中所述矢量出现的时间顺序与用于构成所述矢量的散点出现的时间顺序一致。
图8示出根据本发明实施例的另一种基于心电图的心率分析方法的流程图。如图8所示,对于上述步骤S110,在一种可能的实现方式中,主要可以包括以下步骤:
步骤S1101、对预定时间段内的心电图上的波形进行采样,获得多个心动周期数据;
步骤S1102、将时间上连续的N个心动周期数据分别作为一个散点的一维数据,绘制N维的所述心率散点图,其中N为大于或等于2的整数。
具体而言,首先获取预定时间段内的心电图上的波形,然后对波形进行采样,获得多个心动周期数据。最后将时间上连续的N个心动周期数据分别作为一个散点的一维数据,绘制N维的所述心率散点图,其中N为大于或等于2的整数。
例如,以RR间期数据作为心动周期数据,获得的时间是连续的M段RR间期数据按时间顺序依次为x1,x2,x3,x4,……,xM。绘制出的心率散点图是二维(N=2)的,其散点的坐标为(xi,xi+1),其中i为整数且1<i<M。即心率散点图中的散点的坐标依次为(x1,x2)、(x2,x3)、(x3,x4)、……、(xM-1,xM)。其中,(xi-1,xi)与(xi,xi+1)称为相邻的散点。
在一种可能的实现方式中,上述步骤S1102主要可以包括:
按所述散点出现的时间顺序,对所述散点的绘制颜色进行编码以及以编码后的颜色,绘制各所述散点。
具体地,在按照所得到的散点出现的时间顺序,以颜色渐变的方式例如颜色渐深的方式对进行编码,然后在坐标系中以编码后的各散点的颜色绘制各散点。
对于上述步骤S120,在根据心电图绘制出对应的心率散点图后,可以分别连接所得到的心率散点图中时间上相邻的两个散点,以得到多个矢量。
对于上述步骤S130,在一种可能的实现方式中,主要可以包括以下步骤:
步骤S1301、计算各所述矢量的大小以及角度;
步骤S1302、按所述矢量出现的时间顺序绘制所述栅状图,其中,所述栅状图的横轴表示所述矢量出现的时间顺序,所述栅状图的纵轴表示所述矢量的大小以及角度中任一。
具体而言,在得到多个矢量后,可以分别计算出得到的多个矢量的大小以及角度,如图9a的D1、θ1。
如图9b所示,在一种可能的实现方式中,可以以所得到的矢量的大小D为纵轴,所述矢量出现的时间顺序为横轴建立平面直角坐标系,根据所得到的矢量的大小及其对应出现的时间顺序在所述平面直角坐标系绘制矢量的第一栅状图即长度栅状图。
在一种可能的实现方式中,可以预定一个阈值,然后将所得到的矢量的大小与预定的阈值相比较,从而区分出心动周期变化大的和变化小的心搏,以用于辅助对心律失常的诊断。
如图9c所示,在一种可能的实现方式中,可以以所得到的矢量的角度θ为纵轴,所述矢量出现的时间顺序为横轴建立平面直角坐标系,根据所得到的矢量的角度及其对应出现的时间顺序在所述平面直角坐标系绘制矢量的第二栅状图即角度栅状图。
例如,心率散点图中的散点的坐标为(xi,xi+1),即按出现的时间顺序依次为(x1,x2)、(x2,x3)、(x3,x4)、……、(xM-1,xM)。则连接在时间上相邻的散点,得到的矢量的坐标为(xi-xi-1,xi+1-xi),即按出现顺序依次为(x2-x1,x3-x2)、(x3-x2,x4-x3)、……、(xM-1-xM-2,xM-xM-1)。则矢量的大小Di-1和方向θi-1可以分别用式1和式2表示。
D i - 1 = ( x i - x i - 1 ) 2 + ( x i + 1 - x i ) 2     式1
&theta; i - 1 = arctan | x i - x i - 1 | ( x i - x i - 1 ) 2 + ( x i + 1 - x i ) 2     式2
其中,1<i<M,且i、M为整数。
这样,就能够以所得到的矢量的大小D为纵轴,以得到的矢量出现的时间顺序为横轴建立平面直角坐标系,根据所得到的矢量的大小Di-1及其对应出现的时间顺序在所述平面直角坐标系绘制长度栅状图。并且,能够以的矢量的角度θ为纵轴,以得到的矢量出现的时间顺序为横轴建立平面直角坐标系,根据所得到的矢量的角度θi-1及其对应出现的时间顺序在所述平面直角坐标系绘制角度栅状图。
对于上述步骤S1302,在一种可能的实现方式中,还可以在上述长度栅状图和角度栅状图的横轴上,标识与所述矢量相关联的心动周期数据出现的时间,其中与所述矢量相关联的心动周期数据表示与用于构成所述矢量的任一散点相关联的心动周期数据中任一。
图10a、10b分别示出根据本发明实施例绘制出的正常状态下的长度栅状图、角度栅状图;图11a、11b分别示出根据本发明实施例绘制出的室上性早搏的长度栅状图、角度栅状图;图12a、12b分别示出根据本发明实施例绘制出的室性早搏的长度栅状图、角度栅状图;图13a、13b分别示出根据本发明实施例绘制出的房颤的长度栅状图、角度栅状图;
需要说明的是,根据心率散点图可以单独绘制出长度栅状图,也可以单独绘制出角度栅状图,还可以同时绘制出长度栅状图和角度栅状图,以供医生参考,从而有利于相关疾病的诊断。
根据图10a、10b、图11a、11b、图12a、12b以及图13a、13b所示的栅状图,不仅能辅助医生进行心律失常等疾病诊断,而且在栅状图中能清楚的反映心率与时间点的对应关系,从而在栅状图上能够直观的体现某些疾病如心律失常具体发生的时间。
还需要说明的是,本实施例中的基于心电图的心率分析方法是基于二维的心率散点图,但本领域技术人员能够理解,本发明应不限于此。事实上,用户完全可根据三维、四维……的心率散点,利用本发明实施例中的心电图的栅状图绘制方法,绘制出心电图的栅状图。
例如,如果需要绘制出的心率散点图是三维(N=3)。则同样以RR间期数据作为心动周期数据,获得的时间是连续的M段RR间期数据按时间顺序依次为x1,x2,x3,x4,……,xM。各散点的坐标为(xi-1,xi,xi+1),其中i为整数且1<i<M-1。即心率散点图中的散点的坐标依次为(x1,x2,x3)、(x2,x3,x4)、……、(xM-2,xM-1,xM)。其中,(xi-1,xi,xi+1)与(xi,xi+1,xi+2)称为相邻的散点。连接在时间上相邻的散点,得到的矢量的坐标为(xi-xi-1,xi+1-xi,xi+2-xi+1),即按出现顺序依次为(x2-x1,x3-x2,x4-x3)、(x3-x2,x5-x4,x4-x3)、……、(xM-2-xM-3,xM-1-xM-2,xM-xM-1)。
则矢量的大小用Di-1表示和矢量的方向用θi-1
Figure BDA0000489987380000115
表示,Di-1、θi-1可以分别用式3、式4和式5表示。
D i - 1 = ( x i - x i - 1 ) 2 + ( x i + 1 - x i ) 2 + x i + 2 - x i - 1     式3
&theta; i - 1 = arctan | x i - x i - 1 | ( x i - x i - 1 ) 2 + ( x i + 1 - x i ) 2 + ( x i + 2 - x i + 1 ) 2     式4
Figure BDA0000489987380000113
    式5
其中,1<i<M-1,且i、M为整数。
这样,就能够以所得到的矢量的大小D为纵轴,以得到的矢量出现的时间顺序为横轴建立平面直角坐标系,根据所得到的矢量的大小Di-1及其对应出现的时间顺序在所述平面直角坐标系绘制第一栅状图。并且,能够分别以的矢量的角度θi-1为纵轴,以得到的矢量出现的时间顺序为横轴建立平面直角坐标系,根据所得到的矢量的角度θi-1及其对应出现的时间顺序在所述平面直角坐标系绘制第二栅状图。并且,能够分别以的矢量的角度
Figure BDA0000489987380000114
为纵轴,以得到的矢量出现的时间顺序为横轴建立平面直角坐标系,根据所得到的矢量的角度
Figure BDA0000489987380000121
及其对应出现的时间顺序在所述平面直角坐标系绘制第三栅状图。
根据第一栅状图、第二栅状图和/或第三栅状图,不仅能辅助医生进行心律失常等疾病诊断,而且在栅状图中能清楚的反映心率与时间点的对应关系,从而在第一栅状图、第二栅状图和/或第三栅状图上能够直观的体现某些疾病如心律失常具体发生的时间。
在一种可能的实现方式中,心动周期数据包括以下数据中任一:
RR间期,表示心电图中两个相邻QRS波群中R波之间的时间差;
PP间期,表示心电图中两个相邻P波之间的时间差;
QRS波宽度,表示心电图中一个QRS波群中Q波起始到S波结束的时间差;
PR间期,表示心电图中P波起始到QRS波起始的时间差;
R波幅值,表示心电图中一次R波高度的最大值;
T波幅值,表示心电图中一次T波高度的最大值;
QT间期,表示心电图中QRS波起始到T波结束的时间差。
本发明实施例的基于心电图的心率分析方法,可以在栅状图中清楚的反映心率与时间点的对应关系,这样,医生根据矢量的栅状图在进行心律失常等疾病诊断时,可以清楚、直观的看到疾病所发生的时间点。
图14示出根据本发明一实施例的基于心电图的心率分析设备的结构框图,该基于心电图的心率分析设备可以用来执行本发明实施例的上述任何一种基于心电图的心率分析方法。
如图14所示,上述心率分析设备10主要可以包括散点图绘制模块11、矢量生成模块12以及栅状图绘制模块13。
在一种可能的实现方式中,散点图绘制模块11用于根据心电图绘制心率散点图;矢量生成模块12,与所述散点图绘制模块11连接,用于分别连接所述心率散点图中时间上相邻的两个散点,以得到多个矢量;栅状图绘制模块13,与所述矢量生成模块12连接,用于按所述矢量出现的时间顺序绘制所述矢量的栅状图,其中所述矢量出现的时间顺序与用于构成所述矢量的散点出现的时间顺序一致。
在一种可能的实现方式中,所述栅状图绘制模块13用于:计算各所述矢量的大小以及角度;按所述矢量的出现顺序绘制所述栅状图,其中,所述栅状图的横轴表示所述矢量出现的时间顺序,所述栅状图的纵轴表示所述矢量的大小以及角度中任一。
在一种可能的实现方式中,所述散点图绘制模块11用于:对预定时间段内的心电图上的波形进行采样,获得多个心动周期数据;将时间上连续的N个心动周期数据分别作为一个散点的一维数据,绘制N维的所述心率散点图,其中N为大于或等于2的整数。
在一种可能的实现方式中,所述散点图绘制模块11还用于:按所述散点出现的时间顺序,对所述散点的绘制颜色进行编码;以编码后的颜色,绘制各所述散点。
在一种可能的实现方式中,所述栅状图绘制模块13还用于:在所述栅状图的横轴上,标识与所述矢量相关联的心动周期数据出现的时间。其中与所述矢量相关联的心动周期数据表示与用于构成所述矢量的任一散点相关联的心动周期数据中任一。
在一种可能的实现方式中,所述心动周期数据包括以下数据中任一:
RR间期,表示心电图中两个相邻QRS波群中R波之间的时间差;
PP间期,表示心电图中两个相邻P波之间的时间差;
QRS波宽度,表示心电图中一个QRS波群中Q波起始到S波结束的时间差;
PR间期,表示心电图中P波起始到QRS波起始的时间差;
R波幅值,表示心电图中一次R波高度的最大值;
T波幅值,表示心电图中一次T波高度的最大值;
QT间期,表示心电图中QRS波起始到T波结束的时间差。
本发明实施例的基于心电图的心率分析设备,可以在栅状图中清楚的反映心率与时间点的对应关系,这样,医生根据矢量的栅状图在进行心律失常等疾病诊断时,可以清楚、直观的看到疾病所发生的时间点。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求所述的保护范围为准。

Claims (12)

1.一种基于心电图的心率分析方法,其特征在于,包括:
根据心电图绘制心率散点图;
分别连接所述心率散点图中时间上相邻的两个散点,以得到多个矢量;以及
按所述矢量出现的时间顺序绘制所述矢量的栅状图,其中所述矢量出现的时间顺序与用于构成所述矢量的散点出现的时间顺序一致。
2.根据权利要求1所述的基于心电图的心率分析方法,其特征在于,按所述矢量出现的时间顺序绘制所述矢量的栅状图,包括:
计算各所述矢量的大小以及角度;
按所述矢量的出现顺序绘制所述栅状图,其中,所述栅状图的横轴表示所述矢量出现的时间顺序,所述栅状图的纵轴表示所述矢量的大小以及角度中任一。
3.根据权利要求2所述的基于心电图的心率分析方法,其特征在于,根据心电图绘制心率散点图,包括:
对预定时间段内的心电图上的波形进行采样,获得多个心动周期数据;
将时间上连续的N个心动周期数据分别作为一个散点的一维数据,绘制N维的所述心率散点图,其中N为大于或等于2的整数。
4.根据权利要求3所述的基于心电图的心率分析方法,其特征在于,绘制N维的所述心率散点图,包括:
按所述散点出现的时间顺序,对所述散点的绘制颜色进行编码;
以编码后的颜色,绘制各所述散点。
5.根据权利要求4所述的基于心电图的心率分析方法,其特征在于,按所述矢量的出现顺序绘制所述矢量的栅状图,还包括:
在所述栅状图的横轴上,标识与所述矢量相关联的心动周期数据出现的时间,其中与所述矢量相关联的心动周期数据表示与用于构成所述矢量的任一散点相关联的心动周期数据中任一。
6.根据权利要求1-5中任一项所述的基于心电图的心率分析方法,其特征在于,所述心动周期数据包括以下数据中任一:
RR间期,表示心电图中两个相邻QRS波群中R波之间的时间差;
PP间期,表示心电图中两个相邻P波之间的时间差;
QRS波宽度,表示心电图中一个QRS波群中Q波起始到S波结束的时间差;
PR间期,表示心电图中P波起始到QRS波起始的时间差;
R波幅值,表示心电图中一次R波高度的最大值;
T波幅值,表示心电图中一次T波高度的最大值;
QT间期,表示心电图中QRS波起始到T波结束的时间差。
7.一种基于心电图的心率分析设备,其特征在于,包括:
散点图绘制模块,用于根据心电图绘制心率散点图;
矢量生成模块,与所述散点图绘制模块连接,用于分别连接所述心率散点图中时间上相邻的两个散点,以得到多个矢量;
栅状图绘制模块,与所述矢量生成模块连接,用于按所述矢量出现的时间顺序绘制所述矢量的栅状图,其中所述矢量出现的时间顺序与用于构成所述矢量的散点出现的时间顺序一致。
8.根据权利要求7所述的基于心电图的心率分析设备,其特征在于,所述栅状图绘制模块用于:
计算各所述矢量的大小以及角度;
按所述矢量的出现顺序绘制所述栅状图,其中,所述栅状图的横轴表示所述矢量出现的时间顺序,所述栅状图的纵轴表示所述矢量的大小以及角度中任一。
9.根据权利要求8所述的基于心电图的心率分析设备,其特征在于,所述散点图绘制模块用于:
对预定时间段内的心电图上的波形进行采样,获得多个心动周期数据;
将时间上连续的N个心动周期数据分别作为一个散点的一维数据,绘制N维的所述心率散点图,其中N为大于或等于2的整数。
10.根据权利要求9所述的基于心电图的心率分析设备,其特征在于,所述散点图绘制模块还用于:
按所述散点出现的时间顺序,对所述散点的绘制颜色进行编码;
以编码后的颜色,绘制各所述散点。
11.根据权利要求10所述的基于心电图的心率分析设备,其特征在于,所述栅状图绘制模块还用于:
在所述栅状图的横轴上,标识与所述矢量相关联的心动周期数据出现的时间,其中与所述矢量相关联的心动周期数据表示与用于构成所述矢量的任一散点相关联的心动周期数据中任一。
12.根据权利要求6-11中任一项所述的基于心电图的心率分析设备,其特征在于,所述心动周期数据包括以下数据中任一:
RR间期,表示心电图中两个相邻QRS波群中R波之间的时间差;
PP间期,表示心电图中两个相邻P波之间的时间差;
QRS波宽度,表示心电图中一个QRS波群中Q波起始到S波结束的时间差;
PR间期,表示心电图中P波起始到QRS波起始的时间差;
R波幅值,表示心电图中一次R波高度的最大值;
T波幅值,表示心电图中一次T波高度的最大值;
QT间期,表示心电图中QRS波起始到T波结束的时间差。
CN201410146200.8A 2014-04-11 2014-04-11 基于心电图的心率分析方法及设备 Active CN103892826B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410146200.8A CN103892826B (zh) 2014-04-11 2014-04-11 基于心电图的心率分析方法及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410146200.8A CN103892826B (zh) 2014-04-11 2014-04-11 基于心电图的心率分析方法及设备

Publications (2)

Publication Number Publication Date
CN103892826A true CN103892826A (zh) 2014-07-02
CN103892826B CN103892826B (zh) 2016-01-13

Family

ID=50984643

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410146200.8A Active CN103892826B (zh) 2014-04-11 2014-04-11 基于心电图的心率分析方法及设备

Country Status (1)

Country Link
CN (1) CN103892826B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104586385A (zh) * 2015-02-12 2015-05-06 中国人民解放军总医院 基于心电图的心率分析方法及设备
CN104866579A (zh) * 2015-05-26 2015-08-26 北京海思敏医疗技术有限公司 动态心电图图形数据处理方法、客户端和服务器
CN105844075A (zh) * 2015-01-30 2016-08-10 三星电子株式会社 心血管周期、等值线分析和心率的低功耗、鲁棒估计方法
CN107212854A (zh) * 2016-03-21 2017-09-29 北京麦迪克斯科技有限公司 可穿戴设备及其信息显示方法和装置
KR101799194B1 (ko) 2016-04-18 2017-11-17 계명대학교 산학협력단 심전도 신호를 이용한 부정맥 진단 장치
WO2018082190A1 (zh) * 2016-11-03 2018-05-11 华为技术有限公司 一种ecg信号处理方法及装置
CN108403105A (zh) * 2017-02-09 2018-08-17 深圳市理邦精密仪器股份有限公司 一种心电散点的展示方法及展示装置
CN108784680A (zh) * 2018-03-19 2018-11-13 武汉海星通技术股份有限公司 基于散点图的心电图智能分析方法及系统
CN109567790A (zh) * 2018-12-24 2019-04-05 深圳市瑞康宏业科技开发有限公司 一种心电散点图快速局部分析方法、装置和计算机装置
CN109953754A (zh) * 2019-03-08 2019-07-02 东南大学 一种室性并行心律快速识别的方法以及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03109038A (ja) * 1989-09-23 1991-05-09 Minoru Tawara 運動負荷による心臓診断装置
US5622178A (en) * 1994-05-04 1997-04-22 Spacelabs Medical, Inc. System and method for dynamically displaying cardiac interval data using scatter-plots
CN1843293A (zh) * 2006-03-21 2006-10-11 李方洁 一种长程心电数据分析的方法
US20100106035A1 (en) * 2005-10-31 2010-04-29 Medtronic, Inc. Method of and apparatus for classifying arrhythmias using scatter plot analysis
US20100198092A1 (en) * 2009-01-30 2010-08-05 Isabel Jimenez-Acquarone Method and apparatus for physiological assessment in mammals
CN103006207A (zh) * 2013-01-11 2013-04-03 山东师范大学 基于形态学的t波交替散点图法的心电信号分析方法
CN103271737A (zh) * 2013-05-23 2013-09-04 山东师范大学 基于云模型和散点图的窦性心率震荡趋势提取方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03109038A (ja) * 1989-09-23 1991-05-09 Minoru Tawara 運動負荷による心臓診断装置
US5622178A (en) * 1994-05-04 1997-04-22 Spacelabs Medical, Inc. System and method for dynamically displaying cardiac interval data using scatter-plots
US20100106035A1 (en) * 2005-10-31 2010-04-29 Medtronic, Inc. Method of and apparatus for classifying arrhythmias using scatter plot analysis
CN1843293A (zh) * 2006-03-21 2006-10-11 李方洁 一种长程心电数据分析的方法
US20100198092A1 (en) * 2009-01-30 2010-08-05 Isabel Jimenez-Acquarone Method and apparatus for physiological assessment in mammals
CN103006207A (zh) * 2013-01-11 2013-04-03 山东师范大学 基于形态学的t波交替散点图法的心电信号分析方法
CN103271737A (zh) * 2013-05-23 2013-09-04 山东师范大学 基于云模型和散点图的窦性心率震荡趋势提取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
孙丽琴等: "心电散点图与心率变异分析", 《中国现代药物应用》 *
张芳芳等: "时间散点图诊断心房颤动伴II度房室传导阻滞的探讨", 《中国循环杂志 》 *
贾新未等: "急性心肌梗塞患者心率变异的时域指标、频域指标与非线性分析Lorenz散点图指标的相关性研究", 《河北职工医学院学报》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105844075A (zh) * 2015-01-30 2016-08-10 三星电子株式会社 心血管周期、等值线分析和心率的低功耗、鲁棒估计方法
CN104586385A (zh) * 2015-02-12 2015-05-06 中国人民解放军总医院 基于心电图的心率分析方法及设备
CN104586385B (zh) * 2015-02-12 2017-07-18 中国人民解放军总医院 基于心电图的心率分析方法及设备
CN104866579A (zh) * 2015-05-26 2015-08-26 北京海思敏医疗技术有限公司 动态心电图图形数据处理方法、客户端和服务器
CN107212854A (zh) * 2016-03-21 2017-09-29 北京麦迪克斯科技有限公司 可穿戴设备及其信息显示方法和装置
KR101799194B1 (ko) 2016-04-18 2017-11-17 계명대학교 산학협력단 심전도 신호를 이용한 부정맥 진단 장치
WO2018082190A1 (zh) * 2016-11-03 2018-05-11 华为技术有限公司 一种ecg信号处理方法及装置
US11337638B2 (en) 2016-11-03 2022-05-24 Huawei Technologies Co., Ltd. ECG signal processing method and apparatus
CN108403105A (zh) * 2017-02-09 2018-08-17 深圳市理邦精密仪器股份有限公司 一种心电散点的展示方法及展示装置
CN108784680A (zh) * 2018-03-19 2018-11-13 武汉海星通技术股份有限公司 基于散点图的心电图智能分析方法及系统
CN109567790A (zh) * 2018-12-24 2019-04-05 深圳市瑞康宏业科技开发有限公司 一种心电散点图快速局部分析方法、装置和计算机装置
CN109953754A (zh) * 2019-03-08 2019-07-02 东南大学 一种室性并行心律快速识别的方法以及装置
CN109953754B (zh) * 2019-03-08 2022-02-08 东南大学 一种室性并行心律快速识别的方法以及装置

Also Published As

Publication number Publication date
CN103892826B (zh) 2016-01-13

Similar Documents

Publication Publication Date Title
CN103892826B (zh) 基于心电图的心率分析方法及设备
CN104586385B (zh) 基于心电图的心率分析方法及设备
CN105496402B (zh) 基于散点图和符号动力学的心电特征分析方法
Krishnamoorthi et al. Simulation methods and validation criteria for modeling cardiac ventricular electrophysiology
Allen et al. Signal analysis: time, frequency, scale, and structure
US5483968A (en) Method and apparatus for analyzing the electrical activity of the heart
Aston et al. Beyond HRV: attractor reconstruction using the entire cardiovascular waveform data for novel feature extraction
CN102073451B (zh) 呈现心电图数据的方法
CN103284702A (zh) 心电图和脉搏波关系分析及融合分析的方法和装置
CN102579039B (zh) 一种心电图中的t波交替检测的方法
CN105962926B (zh) 基于心电图的心拍特征分析方法和设备
US20150025403A1 (en) Mood analysis method, system, and apparatus
US20140148719A1 (en) Multi-touch approach for ECG measurement
CN105125208A (zh) 一种心电图绘制方法及心电图绘制工具
CN105595971A (zh) 一种基于视频的脉象信息采集系统及其采集方法
KR20180037525A (ko) 심전도 신호로부터 중요 포인트를 검출하는 장치 및 방법
CN103914871B (zh) 基于点云数据的交互式选取物体表面坐标点的方法
CN105534519A (zh) 彩色心电向量图显示方法
CN104866724A (zh) 基于图像的心电分析方法及装置
CN107908456A (zh) 一种收音机显示面板背景切换方法
CN105054925A (zh) 特征点的获取及波形绘制
US20170143311A1 (en) Systems and methods for pulsed wave predictive processing
CN107374610B (zh) 心磁图生成方法及生成系统
US9901259B2 (en) System for generating cardiac and other electric waveforms
CN101836860B (zh) 一种心脏磁场最大值和最小值的点位确定方法

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
CB03 Change of inventor or designer information

Inventor after: Zou Fang

Inventor after: Gao Xiaofeng

Inventor after: Shen Yanni

Inventor before: Gao Xiaofeng

Inventor before: Shen Yanni

COR Change of bibliographic data