CN111920399A - 心率变异性的分析方法及装置 - Google Patents
心率变异性的分析方法及装置 Download PDFInfo
- Publication number
- CN111920399A CN111920399A CN202010924321.6A CN202010924321A CN111920399A CN 111920399 A CN111920399 A CN 111920399A CN 202010924321 A CN202010924321 A CN 202010924321A CN 111920399 A CN111920399 A CN 111920399A
- Authority
- CN
- China
- Prior art keywords
- electrocardiogram
- clustering
- dimensional
- heart rate
- rate variability
- 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
- 238000004458 analytical method Methods 0.000 title claims description 19
- 238000010586 diagram Methods 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 40
- 230000000747 cardiac effect Effects 0.000 claims abstract description 24
- 230000001746 atrial effect Effects 0.000 claims description 9
- 230000002028 premature Effects 0.000 claims description 9
- 208000000418 Premature Cardiac Complexes Diseases 0.000 claims description 8
- 230000002861 ventricular Effects 0.000 claims description 8
- 206010015856 Extrasystoles Diseases 0.000 claims description 6
- 239000008188 pellet Substances 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000003086 colorant Substances 0.000 claims description 4
- 238000009877 rendering Methods 0.000 claims description 4
- 238000002372 labelling Methods 0.000 claims 1
- 206010003119 arrhythmia Diseases 0.000 abstract description 10
- 230000006793 arrhythmia Effects 0.000 abstract description 10
- 238000004445 quantitative analysis Methods 0.000 abstract description 8
- 238000011160 research Methods 0.000 description 4
- 238000000926 separation method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 208000004301 Sinus Arrhythmia Diseases 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000005183 dynamical system Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 206010003130 Arrhythmia supraventricular Diseases 0.000 description 1
- 206010047281 Ventricular arrhythmia Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002567 autonomic effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000016507 interphase Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 230000002889 sympathetic effect Effects 0.000 description 1
- 210000001186 vagus nerve Anatomy 0.000 description 1
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/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02405—Determining heart rate variability
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Engineering & Computer Science (AREA)
- Physiology (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Artificial Intelligence (AREA)
- Pulmonology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本申请公开了一种心率变异性的分析方法及装置。其中，该方法包括：采集心电图信号；对采集的心电图信号进行分析，得到心电图信号的RR间期，依据心电图信号的RR间期绘制三维心电散点图，其中，RR间期为心电图上两个相邻QRS波的R波时长；将三维心电散点图中的各聚类团块划分为不同属性的聚类团块，其中，各聚类团块为心脏起搏点构成的点集；分别确定不同属性的聚类团块的心率变异性指标。本申请解决了采用传统心电散点图对HRV进行分析时，利用相邻的两个RR间期在平面进行绘制，在心律失常多样化的情况下易出现心搏点集重叠情况发生，导致研究者难以对心电散点图进行彻底有效心率变异性的定量分析的技术问题。
Description
技术领域
本申请涉及心律失常分析及心电信号处理领域，具体而言，涉及一种心率变异性的分析方法及装置。
背景技术
心率变异性(Heart rate variability，HRV)是指窦性心律一定时间内逐次心搏间期时间长度的差异变化。HRV定量分析是评估自主神经和心脏状态的方法。通过对HRV信号的分析处理，可以得到心脏、交感神经、迷走神经等的状态以及相互的制衡情况。近年来众多方法应用到HRV研究中，但是目前广泛应用的主要是传统的时域、频域分析方法。由于人体是一个复杂的非线性系统，心脏本质上也是一个非线性动力系统，因此非线性分析方法更有助于揭示心脏动力系统的本质。
心电散点图作为HRV非线性分析的一种重要手段，早已开始应用到对心率变异性的辅助分析中。而且心电散点图的HRV分析日益受到重视，心电散点图作为一种非线性的分析方法，是描述心率变异性的重要指标。常用的心电散点图使用连续相邻的2个RR间期即(RRn，R Rn+1)作为散点横纵坐标在平面直角坐标系中做散点图。它能反映相邻2个RR间期的变化，在显示HRV整体特征的同时，又能显示出逐次心搏之间的变化，揭示出心率变异性的非线性特征。但是常用的心电散点图是利用相邻的两个RR间期在平面进行绘制，在心律失常多样化的情况下易出现心搏点集重叠情况发生，导致研究者难以对心电散点图进行彻底有效心率变异性的定量分析。此外，目前的心率变异性分析指对窦性心搏间期的变化进行研究，而对室性、房性等心搏并没有开展研究工作。
针对采用传统心电散点图对HRV进行分析时，利用相邻的两个RR间期在平面进行绘制，在心律失常多样化的情况下易出现心搏点集重叠情况发生，导致研究者难以对心电散点图进行彻底有效心率变异性的定量分析的问题，目前尚未提出有效的解决方案。
发明内容
本申请实施例提供了一种心率变异性的分析方法及装置，以至少解决采用传统心电散点图对HRV进行分析时，利用相邻的两个RR间期在平面进行绘制，在心律失常多样化的情况下易出现心搏点集重叠情况发生，导致研究者难以对心电散点图进行彻底有效心率变异性的定量分析的技术问题。
根据本申请实施例的一个方面，提供了一种心率变异性的分析方法，包括：采集心电图信号；对采集的心电图信号进行分析，得到心电图信号的RR间期，依据心电图信号的RR间期绘制三维心电散点图，其中，RR间期为心电图上两个相邻QRS波的R波时长；将三维心电散点图中的各聚类团块划分为不同属性的聚类团块，其中，各聚类团块为心脏起搏点构成的点集；分别确定不同属性的聚类团块的心率变异性指标。
可选地，对采集的心电图信号进行分析，得到心电图信号的RR间期，包括：确定心电图信号的R波的位置，计算相邻两个R波的时长，得到心电图信号的相邻的RR间期；依据心电图信号的RR间期绘制三维心电散点图，包括：依据相邻的三个RR间期迭代绘制三维心电散点图。
可选地，依据相邻的三个RR间期迭代绘制三维心电散点图之前，上述方法还包括：对采集的心电图信号进行分析，将心电图信号至少分为如下类型：正常心搏的心电图信号、室性早搏的心电图信号以及房性早搏的心电图信号。
可选地，依据相邻的三个RR间期迭代绘制三维心电散点图，包括：以当前间期R Rn为坐标X，以间期R Rn+1为坐标Y，以间期R Rn+2为坐标Z，在三维空间坐标系中绘制三维心电散点图；在三维心电散点图中分别以不同的颜色标注不同类型的心电图信号。
可选地，将三维心电散点图中的各聚类团块划分为不同属性的聚类团块，包括：将三维心电散点图中的各聚类团块划分为以下属性的聚类团块：室早团块、房早团块以及正常团块；分别从不同属性的聚类团块中分离出不相容的小聚类团块。
可选地，分别确定不同属性的聚类团块的心率变异性指标，包括：分别计算不同属性的聚类团块及小聚类团块的心率变异性指标；将心率变异性指标对应的计算结果输入至目标表格。
可选地，分别计算不同属性的聚类团块及小聚类团块的心率变异性指标，包括：确定时间范围；旋转三维心电散点图，从旋转之后的三维心电散点图中确定目标聚类团块，目标聚类团块包括从聚类团块和小聚类团块中选定的团块；计算目标聚类团块在时间范围内的各项心率变异性指标。
根据本申请实施例的另一方面，还提供了一种心率变异性的分析装置，包括：采集模块，用于采集心电图信号；绘制模块，用于对采集的心电图信号进行分析，得到心电图信号的RR间期，依据心电图信号的RR间期绘制三维心电散点图，其中，RR间期为心电图上两个相邻QRS波的R波时长；划分模块，用于将三维心电散点图中的各聚类团块划分为不同属性的聚类团块，其中，各聚类团块为心脏起搏点构成的点集；确定模块，用于分别确定不同属性的聚类团块的心率变异性指标。
根据本申请实施例的再一方面，还提供了一种计算机可读存储介质，存储介质包括存储的程序，其中，在程序运行时控制存储介质所在设备执行以上的心率变异性的分析方法。
根据本申请实施例的再一方面，还提供了一种处理器，处理器用于运行存储在存储器中的程序，其中，程序运行时执行以上的心率变异性的分析方法。
在本申请实施例中，采用采集心电图信号；对采集的心电图信号进行分析，得到心电图信号的RR间期，依据心电图信号的RR间期绘制三维心电散点图，其中，RR间期为心电图上两个相邻QRS波的R波时长；将三维心电散点图中的各聚类团块划分为不同属性的聚类团块，其中，各聚类团块为心脏起搏点构成的点集；分别确定不同属性的聚类团块的心率变异性指标的方式，通过对于同一份动态心电数据采用连续相邻的3个RR间期即(R Rn，R Rn+1，RRn+2)在空间坐标系绘制三维心电散点图，在三维心电散点图中计算各团块的心率变异性指标，从而实现了提高心率变异性分析的准确性的技术效果，进而解决了采用传统心电散点图对HRV进行分析时，利用相邻的两个RR间期在平面进行绘制，在心律失常多样化的情况下易出现心搏点集重叠情况发生，导致研究者难以对心电散点图进行彻底有效心率变异性的定量分析技术问题。
附图说明
此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
图1是根据本申请实施例的一种二维心电散点图；
图2是根据本申请实施例的一种三维心电散点图；
图3是根据本申请实施例的一种心率变异性的分析方法的流程图；
图4是根据本申请实施例的一种三维心电散点图；
图5是根据本申请实施例的一种三维心电散点图空间心率变异性分析结果示意图；
图6是根据本申请实施例的一种指定时间的三维心电散点图；
图7是根据本申请实施例的一种指定时间段的三维心电散点图空间心率变异性分析结果示意图；
图8是根据本申请实施例的一种心率变异性的分析装置的结构图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本申请的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
根据本申请实施例，提供了一种心率变异性的分析方法的实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
常用的心电散点图是利用相邻的两个RR间期在平面直角坐标系下绘制的二维心电散点图，在心律失常多样化的情况下易出现不同属性心搏点集的重叠，导致研究者难以对二维心电散点图进行彻底有效心率变异性定量分析，如图1，在二维心电散点图上，室性团块只能观察到一个。为了解决这个问题，对于同一份动态心电数据采用连续相邻的3个RR间期即(R Rn，R Rn+1，R Rn+2)在空间坐标系中绘制三维心电散点图，如图2所示，在三维心电散点图上，通过旋转，可以观察到三个不同的室性聚类团块。三维心电散点图在空间中具备团块分离的优势更易于使用模糊识别算法分离各聚类空间团块，并对分离后的各团块计算心率变异性指标，其计算的心率变异性指标具备空间向量的特性，稳定性更强。利用该方法可进行相关疾病的检测和分类，结论可用于临床监护和远程医疗中。同时本发明也针对室性、房性聚类团块计算各心率变异性指标，用以评估各心律失常事件的心率变异性特征，方便临床对其进行研究。
鉴于HRV信号所具有的非线性特征，以及传统心电散点图分析方法会出现点集重叠的情况发生，本申请提供了一种采用三维心电散点图对空间各聚类团块进行模糊分离并计算各聚类团块心率变异性指标的新方法：
图3是根据本申请实施例的一种心率变异性的分析方法的流程图，如图3所示，该方法包括如下步骤：
步骤S302，采集心电图信号。
根据本申请的一个可选的实施例，在本步骤中对患者进行动态心电数据采集，注意进行皮肤清洁处理，保证采集来的数据干扰较小。
步骤S304，对采集的心电图信号进行分析，得到心电图信号的RR间期，依据心电图信号的RR间期绘制三维心电散点图，其中，RR间期为心电图上两个相邻QRS波的R波时长。
步骤S306，将三维心电散点图中的各聚类团块划分为不同属性的聚类团块，其中，各聚类团块为心脏起搏点构成的点集。
步骤S308，分别确定不同属性的聚类团块的心率变异性指标。
通过上述步骤，通过对于同一份动态心电数据采用连续相邻的3个RR间期即(RRn，R Rn+1，R Rn+2)在空间坐标系绘制三维心电散点图，在三维心电散点图中计算各团块的心率变异性指标，从而实现了提高心率变异性分析的准确性的技术效果。
根据本申请的一个可选的实施例，执行步骤S102时，通过以下方法对采集的心电图信号进行分析，得到心电图信号的RR间期：确定心电图信号的R波的位置，计算相邻两个R波的时长，得到心电图信号的相邻的RR间期。然后依据相邻的三个RR间期迭代绘制三维心电散点图。
在本步骤中，自动分析提取R波位置并计算相邻的RR间期，根据相邻的三个RR间期迭代绘制三维心电散点图。
根据本申请的一个可选的实施例，依据相邻的三个RR间期迭代绘制三维心电散点图之前，还需要对采集的心电图信号进行分析，将心电图信号至少分为如下类型：正常心搏的心电图信号、室性早搏的心电图信号以及房性早搏的心电图信号。
对采集来的心电数据进行自动分析或人工分析，正确标识各心搏属性如：正常心搏、室性早搏、房性早搏等。
在本申请的另一个可选的实施例中，依据相邻的三个RR间期迭代绘制三维心电散点图，包括以下步骤：以当前间期R Rn为坐标X，以间期R Rn+1为坐标Y，以间期R Rn+2为坐标Z，在三维空间坐标系中绘制三维心电散点图；在三维心电散点图中分别以不同的颜色标注不同类型的心电图信号。
图4是根据本申请实施例的一种三维心电散点图，如图4所示，以当前心搏间期RRn为横坐标X，以心搏间期R Rn+1为纵坐Y，以心搏间期R Rn+2为竖坐标Z，在三维空间中绘制三维心电散点图，并用不同的颜色表示不同的心搏属性，如用绿色代表室早心搏，粉色代表房早心搏、黑色代表正常心搏等。
根据本申请的一个可选的实施例，步骤S306通过以下方式实现：将三维心电散点图中的各聚类团块划分为以下属性的聚类团块：室早团块、房早团块以及正常团块；分别从不同属性的聚类团块中分离出不相容的小聚类团块。
在本步骤中，对空间中的各聚类团块根据属性划分为室早团块、房早团块、正常团块，然后对各属性团块使用模糊识别算法分离其中的各聚类小团块。三维心电散点图由连续3心搏RR间期迭代作图生成，其在三维空间中由于吸引子作用，形成了多个既融合又不完全融合的聚类空间团块，本发明就是将这些不同的聚类空间团块通过模糊识别算法自动识别，并分别计算其三维空间中的空间心率变异性指标，这些空间心率变异性指标将用于心脏突发事件及其他相关心律失常的预测。
如图4所示，首先，根据自动识别的心搏属性，将所有室早属性的心搏点集归为室早团块、将所有房早属性的心搏点集归为房早团块、将所有正常属性的心搏点集归为正常团块。然后，以各属性团块点集分别为识别目标，根据团块在三维空间中的分布，先检定团块中的多个高密度分布区域，然后对这些高密度分布区域使用模糊算法检测相容性，最终分离出若干不相容的小聚类团块。
在本申请的另一个可选的实施例中，执行步骤S308时，分别计算不同属性的聚类团块及小聚类团块的心率变异性指标；将心率变异性指标对应的计算结果输入至目标表格。
对各属性团块及各聚类小团块分别进行多项心率变异指标计算，包括但不限于以下指标：三维空间向量值R3d、SDNN3d、rMSSD3d、pNN503d、三角指数3d、VLI3d、VAI3d、TotalPower3d、VLF3d、LF3d、HF3d等指标，并将结果输出为表格。
对各属性团块进行多项心率变异指标进行计算，如
其中，Ra(i)为三维散点图中的空间半径，
对各属性团块分离出来的小聚类团块也分别进行各项心率变异指标的计算，并将结果输出为表格，如图5所示，输出结果每次以心搏修改校准后为基础。
根据本申请的一个可选的实施例，分别计算不同属性的聚类团块及小聚类团块的心率变异性指标，包括：确定时间范围；旋转三维心电散点图，从旋转之后的三维心电散点图中确定目标聚类团块，目标聚类团块包括从聚类团块和小聚类团块中选定的团块；计算目标聚类团块在时间范围内的各项心率变异性指标。
在本步骤中，选择时间范围及指定团块，计算各项心率变异性指标并将结果输出为表格。具体实施时，时间散点图上，选定要分析的时间范围，旋转三维心电散点图，用鼠标圈选任意三维团块，计算所选团块的各项心率变异性指标，并输出为表格。
图6为在时间散点图上随机选择进行绘制的三维心电散点图(不局限在时间散点图这一种选择方式)，通过该选定的时间范围的三维心电散点图去计算团块的各项心率变异性指标，并输出为表格如图7所示。
在三维心电散点图上，通过旋转，可以观察到三个不同的室性聚类团块。三维心电散点图在空间中具备团块分离的优势更易于使用模糊识别算法分离各聚类空间团块，并对分离后的各团块计算心率变异性指标，其计算的心率变异性指标具备空间向量的特性，稳定性更强。
对比现有的基于二维心电散点图的非线性定量分析方法，本申请采用三维心电散点图作为研究对象，利用三维心电散点图具有的团块分离特性，可避免由于二维心电散点图上点集重叠导致的研究目标不准确的缺点。此外在三维心电散点图的计算上，采用分属性、模糊识别各属性小的聚类团块，分别输出各自的心率变异性指标，可以对各团块的空间心率变异指标做出更准确的研究。这些聚类团块包含着窦性、房性、室性心律失常的各种属性，对于当下只研究窦性心率变异的现状提供了更广的研究范围，在过去的研究中已经证实，散点图中某些区域的散点可能与心律失常发作呈相关性，而另一些区域散点集可能具有保护作用。在HRV分析中如果采用空间散点计算各团块的变异性，可能会给某些类型心律失常提早预测从而得到提前干预的效果。
图8是根据本申请实施例的一种心率变异性的分析装置的结构图，如图8所示，该装置包括：
采集模块80，用于采集心电图信号；
绘制模块82，用于对采集的心电图信号进行分析，得到心电图信号的RR间期，依据心电图信号的RR间期绘制三维心电散点图，其中，RR间期为心电图上两个相邻QRS波的R波时长；
划分模块84，用于将三维心电散点图中的各聚类团块划分为不同属性的聚类团块，其中，各聚类团块为心脏起搏点构成的点集；
确定模块86，用于分别确定不同属性的聚类团块的心率变异性指标。
需要说明的是，图8所示实施例的优选实施方式可以参见图3所示实施例的相关描述，此处不再赘述。
本申请实施例还提供了一种计算机可读存储介质，该计算机可读存储介质包括存储的程序，其中，在程序运行时控制计算机可读存储介质所在设备执行以上的心率变异性的分析方法。
计算机可读存储介质用于存储执行以下功能的程序：采集心电图信号；对采集的心电图信号进行分析，得到心电图信号的RR间期，依据心电图信号的RR间期绘制三维心电散点图，其中，RR间期为心电图上两个相邻QRS波的R波时长；将三维心电散点图中的各聚类团块划分为不同属性的聚类团块，其中，各聚类团块为心脏起搏点构成的点集；分别确定不同属性的聚类团块的心率变异性指标。
本申请实施例还提供了一种处理器，处理器用于运行程序，其中，在程序运行时执行以上的心率变异性的分析方法。
处理器用于运行执行以下功能的程序采集心电图信号；对采集的心电图信号进行分析，得到心电图信号的RR间期，依据心电图信号的RR间期绘制三维心电散点图，其中，RR间期为心电图上两个相邻QRS波的R波时长；将三维心电散点图中的各聚类团块划分为不同属性的聚类团块，其中，各聚类团块为心脏起搏点构成的点集；分别确定不同属性的聚类团块的心率变异性指标。
上述本申请实施例序号仅仅为了描述，不代表实施例的优劣。
在本申请的上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
在本申请所提供的几个实施例中，应该理解到，所揭露的技术内容，可通过其它的方式实现。其中，以上所描述的装置实施例仅仅是示意性的，例如所述单元的划分，可以为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，单元或模块的间接耦合或通信连接，可以是电性或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、只读存储器(ROM，ReDMRJd-Only Memory)、随机存取存储器(RDMRJM，RDMRJndom DMRJccess Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅是本申请的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本申请原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本申请的保护范围。
Claims (10)
1.一种心率变异性的分析方法，其特征在于，包括：
采集心电图信号；
对采集的所述心电图信号进行分析，得到所述心电图信号的RR间期，依据所述心电图信号的RR间期绘制三维心电散点图，其中，所述RR间期为心电图上两个相邻QRS波的R波时长；
将所述三维心电散点图中的各聚类团块划分为不同属性的聚类团块，其中，所述各聚类团块为心脏起搏点构成的点集；
分别确定所述不同属性的聚类团块的心率变异性指标。
2.根据权利要求1所述的方法，其特征在于，
对采集的所述心电图信号进行分析，得到所述心电图信号的RR间期，包括：确定所述心电图信号的R波的位置，计算相邻两个R波的时长，得到所述心电图信号的相邻的RR间期；
依据所述心电图信号的RR间期绘制三维心电散点图，包括：依据相邻的三个RR间期迭代绘制所述三维心电散点图。
3.根据权利要求2所述的方法，其特征在于，依据相邻的三个RR间期迭代绘制所述三维心电散点图之前，所述方法还包括：
对采集的所述心电图信号进行分析，将所述心电图信号至少分为如下类型：正常心搏的心电图信号、室性早搏的心电图信号以及房性早搏的心电图信号。
4.根据权利要求3所述的方法，其特征在于，依据相邻的三个RR间期迭代绘制所述三维心电散点图，包括：
以当前间期R Rn为坐标X，以间期R Rn+1为坐标Y，以间期R Rn+2为坐标Z，在三维空间坐标系中绘制所述三维心电散点图；
在所述三维心电散点图中分别以不同的颜色标注不同所述类型的心电图信号。
5.根据权利要求1所述的方法，其特征在于，将所述三维心电散点图中的各聚类团块划分为不同属性的聚类团块，包括：
将所述三维心电散点图中的各聚类团块划分为以下属性的聚类团块：室早团块、房早团块以及正常团块；
分别从不同属性的所述聚类团块中分离出不相容的小聚类团块。
6.根据权利要求5所述的方法，其特征在于，分别确定所述不同属性的聚类团块的心率变异性指标，包括：
分别计算不同属性的所述聚类团块及所述小聚类团块的心率变异性指标；
将所述心率变异性指标对应的计算结果输入至目标表格。
7.根据权利要求6所述的方法，其特征在于，分别计算不同属性的所述聚类团块及所述小聚类团块的心率变异性指标，包括：
确定时间范围；
旋转所述三维心电散点图，从旋转之后的所述三维心电散点图中确定目标聚类团块，所述目标聚类团块包括从所述聚类团块和所述小聚类团块中选定的团块；
计算所述目标聚类团块在所述时间范围内的各项心率变异性指标。
8.一种心率变异性的分析装置，其特征在于，包括：
采集模块，用于采集心电图信号；
绘制模块，用于对采集的所述心电图信号进行分析，得到所述心电图信号的RR间期，依据所述心电图信号的RR间期绘制三维心电散点图，其中，所述RR间期为心电图上两个相邻QRS波的R波时长；
划分模块，用于将所述三维心电散点图中的各聚类团块划分为不同属性的聚类团块，其中，所述各聚类团块为心脏起搏点构成的点集；
确定模块，用于分别确定所述不同属性的聚类团块的心率变异性指标。
9.一种计算机可读存储介质，其特征在于，所述存储介质包括存储的程序，其中，在所述程序运行时控制所述存储介质所在设备执行权利要求1至7中任意一项所述的心率变异性的分析方法。
10.一种处理器，其特征在于，所述处理器用于运行存储在存储器中的程序，其中，所述程序运行时执行权利要求1至7中任意一项所述的心率变异性的分析方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010924321.6A CN111920399A (zh) | 2020-09-04 | 2020-09-04 | 心率变异性的分析方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010924321.6A CN111920399A (zh) | 2020-09-04 | 2020-09-04 | 心率变异性的分析方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111920399A true CN111920399A (zh) | 2020-11-13 |
Family
ID=73308676
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010924321.6A Pending CN111920399A (zh) | 2020-09-04 | 2020-09-04 | 心率变异性的分析方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111920399A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112842350A (zh) * | 2021-02-20 | 2021-05-28 | 无锡市中健科仪有限公司 | 基于正反传播算法获得心拍的主导心率的方法 |
| CN113273986A (zh) * | 2021-05-19 | 2021-08-20 | 于杨 | 一种心率变异性分析方法、装置、智能终端及存储介质 |
| CN113647960A (zh) * | 2021-07-22 | 2021-11-16 | 厦门纳龙科技有限公司 | 一种心率失常的快速识别方法、终端设备及存储介质 |
| CN117958835A (zh) * | 2024-03-28 | 2024-05-03 | 北华大学 | 用于患者术后康复的心电图分析方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10233149A1 (de) * | 2002-07-22 | 2004-02-12 | Ganimed Gmbh | Digitales, elektronisches Diagnose-Verfahren und Meßeinrichtung zu seiner Durchführung |
| CN102314684A (zh) * | 2011-07-18 | 2012-01-11 | 复旦大学 | 评价Lorenz散点图离散程度的量化方法 |
| CN206792400U (zh) * | 2017-01-16 | 2017-12-26 | 吉林东华原医疗设备有限责任公司 | 心率变异性检测装置 |
| CN108836317A (zh) * | 2018-06-27 | 2018-11-20 | 广州迪茂信息科技有限公司 | 动态心电时间散点图逆向分析方法 |
| CN110403598A (zh) * | 2019-08-05 | 2019-11-05 | 迪姆软件(北京)有限公司 | 三维心电散点图的处理方法及装置 |
-
2020
- 2020-09-04 CN CN202010924321.6A patent/CN111920399A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10233149A1 (de) * | 2002-07-22 | 2004-02-12 | Ganimed Gmbh | Digitales, elektronisches Diagnose-Verfahren und Meßeinrichtung zu seiner Durchführung |
| CN102314684A (zh) * | 2011-07-18 | 2012-01-11 | 复旦大学 | 评价Lorenz散点图离散程度的量化方法 |
| CN206792400U (zh) * | 2017-01-16 | 2017-12-26 | 吉林东华原医疗设备有限责任公司 | 心率变异性检测装置 |
| CN108836317A (zh) * | 2018-06-27 | 2018-11-20 | 广州迪茂信息科技有限公司 | 动态心电时间散点图逆向分析方法 |
| CN110403598A (zh) * | 2019-08-05 | 2019-11-05 | 迪姆软件(北京)有限公司 | 三维心电散点图的处理方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| 刘克钦: "心率变异性与Lorenz图分析研究", 中国优秀硕士学位论文全文数据库, no. 1, pages 33 - 53 * |
| 胡敏 等: "心电三维RR间期散点图的构建及识别", 中国心脏起搏与心电生理杂志, vol. 30, no. 5, pages 455 - 457 * |
| 胡敏，江成璠，王素霞: "三维心电散点图吸引子的析取和识别", 生物医学工程学杂志, vol. 35, no. 1, pages 25 - 30 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112842350A (zh) * | 2021-02-20 | 2021-05-28 | 无锡市中健科仪有限公司 | 基于正反传播算法获得心拍的主导心率的方法 |
| CN113273986A (zh) * | 2021-05-19 | 2021-08-20 | 于杨 | 一种心率变异性分析方法、装置、智能终端及存储介质 |
| CN113273986B (zh) * | 2021-05-19 | 2021-12-21 | 于杨 | 一种心率变异性分析方法、装置、智能终端及存储介质 |
| CN113647960A (zh) * | 2021-07-22 | 2021-11-16 | 厦门纳龙科技有限公司 | 一种心率失常的快速识别方法、终端设备及存储介质 |
| CN113647960B (zh) * | 2021-07-22 | 2024-05-14 | 纳龙健康科技股份有限公司 | 一种心律失常的快速识别方法、终端设备及存储介质 |
| CN117958835A (zh) * | 2024-03-28 | 2024-05-03 | 北华大学 | 用于患者术后康复的心电图分析方法 |
| CN117958835B (zh) * | 2024-03-28 | 2024-06-11 | 北华大学 | 用于患者术后康复的心电图分析方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11147500B2 (en) | Electrocardiogram processing system for delineation and classification | |
| CN111920399A (zh) | 心率变异性的分析方法及装置 | |
| JP6396253B2 (ja) | 電気生理学データの視覚化 | |
| EP3612091A1 (en) | Electroanatomical mapping tools facilitated by activation waveforms | |
| CN110464333B (zh) | 一种心电图数据的存储方法及装置 | |
| CN113876335B (zh) | 心电图心搏分类方法、装置、电子设备和介质 | |
| CN108403107B (zh) | 一种心律失常判别方法及系统 | |
| TW202100095A (zh) | 心臟圖形顯示系統 | |
| CN105962926B (zh) | 基于心电图的心拍特征分析方法和设备 | |
| JP2023067876A (ja) | 心臓グラフィックディスプレイシステム | |
| CN109044348A (zh) | 房颤检测装置及存储介质 | |
| Ganesan et al. | Computer-aided clinical decision support systems for atrial fibrillation | |
| US11672464B2 (en) | Electrocardiogram processing system for delineation and classification | |
| Gillette et al. | Medalcare-xl: 16,900 healthy and pathological 12 lead ecgs obtained through electrophysiological simulations | |
| CN109431488B (zh) | 一种心电图动态实时分析数据中典型数据的识别方法 | |
| CN109044338A (zh) | 房颤检测装置及存储介质 | |
| WO2011094105A2 (en) | System and associated methods for quantitative assessment of cardiac electrical events | |
| Bergquist et al. | Body Surface Potential Mapping: Contemporary Applications and Future Perspectives. Hearts 2021, 2, 514–542 | |
| US12011273B2 (en) | Method and device for the technical support of the analysis of signals acquired by measurement, the signals having a time- and space-dependent signal characteristic | |
| Sun et al. | Multi-scale Joint Recurrence Quantification Analysis Integrating ECG Spatiotemporal and Dynamic Information for Cardiopathy Detection | |
| Winslow et al. | OPEN ACCESS EDITED AND REVIEWED BY | |
| Doste Beltrán et al. | Training machine learning models with synthetic data improves the prediction of ventricular origin in outflow tract ventricular arrhythmias | |
| Rao et al. | Atrial Fibrillation and Atrial Flutter detection using Electrocardiomatrix technique from ECG signal | |
| CN115990024A (zh) | 使用致心律失常激活之间的相似性对电生理(ep)信号进行集群 | |
| Timofeeva et al. | System for Automatic Representation of Correlation Rhythmogram Dynamics of Long-Term ECG Signal Recordings |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |