CN109480825B - 心电数据的处理方法及装置 - Google Patents

心电数据的处理方法及装置 Download PDF

Info

Publication number
CN109480825B
CN109480825B CN201811530312.8A CN201811530312A CN109480825B CN 109480825 B CN109480825 B CN 109480825B CN 201811530312 A CN201811530312 A CN 201811530312A CN 109480825 B CN109480825 B CN 109480825B
Authority
CN
China
Prior art keywords
data
heart beat
sample set
heartbeat
electrocardiogram
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
CN201811530312.8A
Other languages
English (en)
Other versions
CN109480825A (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.)
Wuhan Zoncare Bio Medical Electronics Co ltd
Original Assignee
Wuhan Zoncare Bio Medical Electronics Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan Zoncare Bio Medical Electronics Co ltd filed Critical Wuhan Zoncare Bio Medical Electronics Co ltd
Priority to CN201811530312.8A priority Critical patent/CN109480825B/zh
Publication of CN109480825A publication Critical patent/CN109480825A/zh
Application granted granted Critical
Publication of CN109480825B publication Critical patent/CN109480825B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques

Landscapes

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

Abstract

本发明提供了一种心电数据的处理方法及装置,涉及数据处理的技术领域,能够获取已知类别的心电数据库,其中,心电数据库包括多种类型的心电数据样本,每种类型的心电数据样本包括多个心电数据,每个心电数据携带有诊断标记;提取每种类型中每个心电数据的心拍数据,按照预先设置的数据长度对每个心拍数据进行插值和采样操作,生成每个类型对应的心拍样本集,按照预先设定的心电形态类型对心拍样本集中包括的多个心拍数据进行聚类生成聚类结果,对聚类后的心拍样本集进行重编码,生成编码心拍样本集,利用编码心拍样本集训练神经网络模型,实现对待测心电数据进行准确分类,有效缓解了当前依赖机器给出诊断结果容易造成心电数据误诊的技术问题。

Description

心电数据的处理方法及装置
技术领域
本发明涉及数据处理技术领域,尤其是涉及一种心电数据的处理方法及装置。
背景技术
心电图是由窦房结自主起搏后,一系列的电位变化通过传导系统,传导至心脏各部分形成心肌的电生理活动。按照心脏激动的时间顺序,将体表电位的变化记录下来,形成的连续曲线即为心电图。典型的心电图包括P波、QRS波、T波。P波反映心房除极过程的电位变化;P-R间期代表激动从窦房结通过房室交界区到心室肌开始除极的时限;QRS波群反映心室除极过程的电位变化;T波代表心室肌复极过程中的电位变化。
心电图具有很强的复杂性,不同种族、性别、年龄的人在各种病理情况下的差异性很大,即使同一人在不同时刻的心电图表现的类型也不一样。目前,由于医生的知识专业性和经验积累不足,常常会过于依赖机器给出的自动诊断结果,很容易造成心电数据的误诊。
发明内容
有鉴于此,本发明的目的在于提供一种心电数据的处理方法及装置,以缓解上述的技术问题。
第一方面,本发明实施例提供了一种心电数据的处理方法,其中,该方法包括:获取已知类别的心电数据库,其中,心电数据库包括多种类型的心电数据样本,每种类型的心电数据样本包括多个心电数据,每个心电数据携带有诊断标记;提取每种类型中每个心电数据的心拍数据,按照预先设置的数据长度对每个心拍数据进行插值和采样操作,生成每个类型对应的心拍样本集,其中,心拍样本集中包含多个长度一致的心拍数据;按照预先设定的心电形态类型对心拍样本集中包括的多个心拍数据进行聚类,生成聚类结果,其中,聚类结果包含心电形态类型,以及,每个心电形态类型包含的心拍数据;对聚类后的心拍样本集进行重编码,生成编码心拍样本集,以完成心电数据的处理过程。
结合第一方面,本发明实施例提供了第一方面的第一种可能的实施方式,其中,上述方法还包括:利用编码心拍样本集对神经网络模型进行训练,以生成用于对待测心电数据进行分类的心电信号分类模型。
结合第一方面的第一种可能的实施方式,本发明实施例提供了第一方面的第二种可能的实施方式,其中,利用编码心拍样本集对神经网络模型进行训练的步骤包括:设置神经网络模型的输入层、隐藏层和输出层;其中,输出层的的长度至少包括以下类型的输出层:正常心电、心律失常、室内传导阻滞、心室肥大、ST段异常、心肌梗死、心电轴偏移。
结合第一方面,本发明实施例提供了第一方面的第三种可能的实施方式,其中,提取每种类型中每个心电数据的心拍数据的步骤包括:对获取的心电数据样本中包括的多个心电数据进行滤波处理;查找滤波处理后的心电数据中的R波峰值点,以获取心电数据的R峰位置序列;提取R峰位置序列中相邻两个R波峰值点,将相邻两个R波峰值点分别作为单个心拍的首尾点,相邻两个R波峰值点之间的心电数据作为一个周期的心拍,生成每个心电数据的心拍数据。
结合第一方面,本发明实施例提供了第一方面的第四种可能的实施方式,其中,按照预先设置的数据长度对每个心拍数据进行插值和采样操作的步骤包括:对每个心拍数据进行插值处理;利用预先设置的数据长度和插值处理后心拍数据计算采样间隔;根据采样间隔对插值处理后的心拍数据进行采样处理,以生成长度一致的心拍数据。
结合第一方面,本发明实施例提供了第一方面的第五种可能的实施方式,其中,按照预先设定的心电形态类型对心拍样本集中包括的多个心拍数据进行聚类之后,该方法还包括:将心拍样本集中的心拍数量与预先设定的心拍数量进行比较;如果心拍样本集中的心拍数量等于预先设定的心拍数量,将心拍样本集中的每个心拍按照所属心电形态类型进行数值标识;如果心拍样本集中的心拍数量大于预先设定的心拍数量,按照心电数据的R峰位置序列顺序去掉心拍样本集中比预先设定的心拍数量多的心拍,使心拍样本集中的心拍数量与预先设定的心拍数量相同,并将心拍样本集中的每个心拍按照所属心电形态类型进行数值标识;如果心拍样本集中的心拍数量小于预先设定的心拍数量,按照心电数据的R峰位置序列顺序扩大心拍样本集中比预先设定的心拍数量少的心拍,使心拍样本集中的心拍数量与预先设定的心拍数量相同,并将心拍样本集中的每个心拍按照所属心电形态类型进行数值标识。
结合第一方面,本发明实施例提供了第一方面的第六种可能的实施方式,其中,对聚类后的心拍样本集进行重编码,生成编码心拍样本集的步骤包括:将每个心电形态类型包含的心拍数据进行求和平均计算,以获取每个心电形态类型包含的新的心拍数据;将每个心电形态类型包含的新的心拍数据以及心拍样本集中的每个心拍的标识数值进行融合,生成编码心拍样本集。
第二方面,本发明实施例还提供一种心电数据的处理装置,其中,该装置包括:获取模块,用于获取已知类别的心电数据库,其中,心电数据库包括多种类型的心电数据样本,每种类型的心电数据样本包括多个心电数据,每个心电数据携带有诊断标记;提取模块,用于提取每种类型中每个心电数据的心拍数据,按照预先设置的数据长度对每个心拍数据进行插值和采样操作,生成每个类型对应的心拍样本集,其中,心拍样本集中包含多个长度一致的心拍数据;聚类模块,用于按照预先设定的心电形态类型对心拍样本集中包括的多个心拍数据进行聚类,生成聚类结果,其中,聚类结果包含心电形态类型,以及,每个心电形态类型包含的心拍数据;编码模块,用于对聚类后的心拍样本集进行重编码,生成编码心拍样本集,以完成心电数据的处理过程。
结合第二方面,本发明实施例提供了第二方面的第一种可能的实施方式,其中,上述装置还包括:训练模块,用于利用编码心拍样本集对神经网络模型进行训练,以生成用于对待测心电数据进行分类的心电信号分类模型。
结合第二方面的第一种可能的实施方式,本发明实施例提供了第二方面的第二种可能的实施方式,其中,训练模型用于:设置神经网络模型的输入层、隐藏层和输出层;其中,输出层的的长度至少包括以下类型的输出层:正常心电、心律失常、室内传导阻滞、心室肥大、ST段异常、心肌梗死、心电轴偏移。
本发明实施例带来了以下有益效果:
本发明实施例提供的一种心电数据的处理方法及装置,能够获取已知类别的心电数据库,其中,心电数据库包括多种类型的心电数据样本,每种类型的心电数据样本包括多个心电数据,每个心电数据携带有诊断标记;提取每种类型中每个心电数据的心拍数据,按照预先设置的数据长度对每个心拍数据进行插值和采样操作,生成每个类型对应的心拍样本集,其中,心拍样本集中包含多个长度一致的心拍数据;按照预先设定的心电形态类型对心拍样本集中包括的多个心拍数据进行聚类,生成聚类结果,其中,聚类结果包含心电形态类型,以及,每个心电形态类型包含的心拍数据;对聚类后的心拍样本集进行重编码,生成编码心拍样本集,以完成心电数据的处理过程,利用编码心拍样本集训练神经网络模型,实现对待测心电数据进行分类,有效缓解了当前依赖机器给出诊断结果容易造成心电数据误诊的技术问题。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种心电数据的处理方法的流程图;
图2为本发明实施例提供的另一种心电数据的处理方法的流程图;
图3为本发明实施例提供的一种心电数据的处理装置的结构示意图;
图4为本发明实施例提供的另一种心电数据的处理装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前,由于对心电数据的分析过分依赖于机器给出诊断结果很容易造成误诊的情况,因此,许多学者开始采用不同的方法对心电数据进行分析,对心电数据分析的本质是对心电数据进行分类,通常是寻找不同的心电数据在时域、频域里所包含的不同特征信息,特征提取的种类越多,分类的准确率越高。采用时域、频率的特征点作为网络的输入端,根据标签样本,进行有监督的学习,如模板匹配、最邻近算法、决策树、支持向量机、BP(BackPropagation,反向传播)神经网络的方法;或无监督的训练,聚类分析、隐马尔科夫模型、高斯混合模型。常见的方法均基于MIT的心律失常数据库中标记的心拍类型识别、正常异常信号识别或针对于某一种特殊的病理心电信号进行识别分类,受限于心电数据样本的个数,方法的普适性有很大的局限性,迟迟未有较大的进展,基于此,本发明实施例提供的一种心电数据的处理方法及装置,可以实现对待测心电数据进行准确分类,有效缓解了当前依赖机器给出诊断结果容易造成心电数据误诊的技术问题。
为便于对本实施例进行理解,首先对本发明实施例所公开的一种心电数据的处理方法进行详细介绍。
实施例一:
本发明实施例提供了一种心电数据的处理方法,如图1所示的一种心电数据的处理方法的流程图,该方法包括以下步骤:
步骤S102,获取已知类别的心电数据库,其中,心电数据库包括多种类型的心电数据样本,每种类型的心电数据样本包括多个心电数据,每个心电数据携带有诊断标记。
具体实现时,选取正常心电、心律失常、室内传导阻滞、心室肥大、ST段异常、心肌梗死、心电轴偏移七大类型的心电数据样本,每种类型样本50万,并对每个心电数据样本所属诊断类型进行标记。本发明实施例,不对心电数据样本所属诊断类型以及各个诊断类型选取的样本数量进行限定。
步骤S104,提取每种类型中每个心电数据的心拍数据,按照预先设置的数据长度对每个心拍数据进行插值和采样操作,生成每个类型对应的心拍样本集,其中,心拍样本集中包含多个长度一致的心拍数据。
具体地,以12导联心电图为例进行说明,心脏是一个立体的结构,为了反应心脏不同面的电活动,在人体不同部位放置电极,以记录和反应心脏的电活动。在常规心电图检查时,通常只安放4个肢体导联电极和V1~V66个胸前导联电极,两两电极之间或电极与中央电势端之间组成一个个不同的导联,通过导联线与心电图机电流计的正负极相连,记录心脏的电活动。两个电极之间组成了双极导联,一个导联为正极,一个导联为负极。双极肢体导联包括Ⅰ导联,Ⅱ导联和Ⅲ导联;电极和中央电势端之间构成了单极导联,此时探测电极为正极,中央电势端为负极。avR、avL、avF、V1、V2、V3、V4、V5、和V6导联均为单极导联。Ⅰ导联、Ⅱ导联、Ⅲ导联、avR导联、avL导联、avF导联、V1导联、V2导联、V3导联、V4导联、V5导联和V6导联构成常规的12导联心电图,因此,每个心电数据data如下所示:
Figure BDA0001904871550000071
其中,sig1,sig2,...,sig12分别表示Ⅰ导联、Ⅱ导联、Ⅲ导联、avR导联、avL导联、avF导联、V1导联、V2导联、V3导联、V4导联、V5导联和V6导联的心电数据。在一定时间内从每个导联的心电数据中提取心拍构成单个导联的心拍样本集,12个导联的心拍样本集构成一个心电数据样本的整体心拍样本集,并对提取的每个心拍进行插值和采样处理,使得所有的心拍长度都为固定值N。在本发明实施例中,不对心电图的导联数量以及心拍长度进行限定。
步骤S106,按照预先设定的心电形态类型对心拍样本集中包括的多个心拍数据进行聚类,生成聚类结果,其中,聚类结果包含心电形态类型,以及,每个心电形态类型包含的心拍数据。
在实际应用中,为了便于对提取的心拍进行处理,首先对每个导联提取出的心拍进行聚类,现阶段聚类算法的种类有很多,如:基于划分的k-means聚类、根据层次分解的层次聚类、基于神经网络的SOM(Self-organizing Maps,自组织映射)聚类、模糊聚类中的FCM(Fuzzy C-Means,模糊聚类)算法,而在本发明实施例中采用运行效率高,原理清晰,实现方便的k-means算法对每个导联的心拍根据预先设定的聚类数量进行聚类。单个导联中每周期的心电形态类型数不大于3个,在本发明实施例中将聚类类别数统一定为3类,分别是正常心拍类型、异常心拍1类型以及异常心拍2类型,并运用k-means算法分别对每个导联中的心拍样本集进行聚类运算,得到3类经过聚类后的心拍数据以及每个心拍对应的所属聚类类型。本发明实施例,不对聚类算法进行限定。
步骤S108,对聚类后的心拍样本集进行重编码,生成编码心拍样本集,以完成心电数据的处理过程。
具体实现时,将每个导联聚类出的三种类型心拍数据与每个心拍对应的所属聚类类型对应的标记数值进行融合,构成新的心电数据,去除了12导联数据中的冗余数据,突出了心电数据的重要关注点,为后面的模型训练做了充分的数据准备,不仅大大减少了模型训练时的数据量,而且极大的降低了训练时间,提高心电数据分类的准确度。
本发明实施例提供的一种心电数据的处理方法,能够获取已知类别的心电数据库,其中,心电数据库包括多种类型的心电数据样本,每种类型的心电数据样本包括多个心电数据,每个心电数据携带有诊断标记;提取每种类型中每个心电数据的心拍数据,按照预先设置的数据长度对每个心拍数据进行插值和采样操作,生成每个类型对应的心拍样本集,其中,心拍样本集中包含多个长度一致的心拍数据;按照预先设定的心电形态类型对心拍样本集中包括的多个心拍数据进行聚类,生成聚类结果,其中,聚类结果包含心电形态类型,以及,每个心电形态类型包含的心拍数据;对聚类后的心拍样本集进行重编码,生成编码心拍样本集,以完成心电数据的处理过程,利用编码心拍样本集训练神经网络模型,实现对待测心电数据进行分类,有效缓解了当前依赖机器给出诊断结果容易造成心电数据误诊的技术问题。
进一步,基于上述心电数据的处理方法的过程,图2示出了另一种心电数据的处理方法的流程图;其中,步骤S202到步骤S208的过程可以参考上述图1中步骤S102到步骤S108对应的过程,在此不再赘述。如图2所示,该方法还包括以下步骤:
步骤S210,利用编码心拍样本集对神经网络模型进行训练,以生成用于对待测心电数据进行分类的心电信号分类模型。
通常,在机器学习中,计算的过程将输入信号进行一系列的计算操作,完成从输入空间到输出空间的映射和变换,具体为升维、降维、放大、缩小、旋转、平移、弯曲这些元素的相似、组合、非线性的计算;传统的诊断也是对输入的心电信号的抽象组合。机器学习按照计算方法的不同,将数据进过投影、变换、映射,将输出空间的计算结果与输出标签进行不断的迭代、反向传播计算、修正,最终得到一系列表征不同输出标签的数据分布组合。在本发明实施例中,利用编码心拍样本集通过对CNN(Convolution Neural Network,卷积神经网络)模型进行训练得到心电信号分类模型,以对待测心电数据进行分类。
进一步,利用编码心拍样本集对神经网络模型进行训练的步骤包括:设置神经网络模型的输入层、隐藏层和输出层;其中,输出层的的长度至少包括以下类型的输出层:正常心电、心律失常、室内传导阻滞、心室肥大、ST段异常、心肌梗死、心电轴偏移。
具体地,在CNN整个网络中,可以分为输入层、隐藏层、输出层三大的模块,其中输入层中,其卷积为64个12*32的卷积核,经过卷积后的数据进行max_pooling池化;隐藏层由两层卷积层与一层神经元个数为1024的全连接层组成,由于心电数据类型有七大类型分别为:正常心电、心律失常、室内传导阻滞、心室肥大、ST段异常、心肌梗死以及心电轴偏移,因此,设置网络的输出层长度为7。
在实际使用时,提取每种类型中每个心电数据的心拍数据的步骤包括:对获取的心电数据样本中包括的多个心电数据进行滤波处理;查找滤波处理后的心电数据中的R波峰值点,以获取心电数据的R峰位置序列;提取R峰位置序列中相邻两个R波峰值点,将相邻两个R波峰值点分别作为单个心拍的首尾点,相邻两个R波峰值点之间的心电数据作为一个周期的心拍,生成每个心电数据的心拍数据。
具体地,先对心电数据进行滤波,去除冗余数据,找出每导联的R波峰值点,进而得到每导联R峰位置序列,其中导联i的R峰位置序列为Rseqi:Rseqi=[Ri,1,...,Ri,M],其中,[Ri,1,...,Ri,M]分别表示的是导联i中各个R波峰处对应的位置;M表示导联i中R波峰数量。对每导联取R位置序列中两两相邻R波峰值点分别作为单个心拍的首尾点,两两相邻R波峰值点间数据作为一个周期的心拍。如Rseqi所示,则导联i中具有M-1个周期的心拍。则导联i中的各个心拍为:
[sigi(Rseq1:Rseq2),sigi(Rseq2:Rseq3),...,sigi(Rseqm-1:Rseqm)],其中,sigi(Rseq1:Rseq2)
表示的是导联i的信号开始于Rseq1,结束于Rseq2,此为导联i中第一个心拍的数据。
具体实现时,按照预先设置的数据长度对每个心拍数据进行插值和采样操作的步骤包括:对每个心拍数据进行插值处理;利用预先设置的数据长度和插值处理后心拍数据计算采样间隔;根据采样间隔对插值处理后的心拍数据进行采样处理,以生成长度一致的心拍数据。
在实际使用时,由于采集到的心电数据均为临床采集所得,由于人体心脏信号的固有特性,不可能每个RR间期的长度完全一致,在此对每个心拍信号的进行插值后下采样,使得到的每个心拍信号数据长度一致,单心拍数据长度不妨设为N。具体的,导联i第1个心拍数据长度为Len_Beati1=Rseq2-Rseq1,首先对心拍数据sigi(Rseq1:Rseq2)进行插值,插值后数据长度为两者的最小公倍数tmp_num:tmp_num=min_gongpei(Rseq2-Rseq1,N),其中,min_gongpei为求这两个数值的最小公倍数函数。然后对插值后数据进行下采样,下采用间隔为tmp_num/N,得到下采样后导联i第1个心拍长度为N;其他心拍重复以上操作,使得各个单心拍长度统一为N。
具体地,按照预先设定的心电形态类型对心拍样本集中包括的多个心拍数据进行聚类之后,该方法还包括:将心拍样本集中的心拍数量与预先设定的心拍数量进行比较;如果心拍样本集中的心拍数量等于预先设定的心拍数量,将心拍样本集中的每个心拍按照所属心电形态类型进行数值标识;如果心拍样本集中的心拍数量大于预先设定的心拍数量,按照心电数据的R峰位置序列顺序去掉心拍样本集中比预先设定的心拍数量多的心拍,使心拍样本集中的心拍数量与预先设定的心拍数量相同,并将心拍样本集中的每个心拍按照所属心电形态类型进行数值标识;如果心拍样本集中的心拍数量小于预先设定的心拍数量,按照心电数据的R峰位置序列顺序扩大心拍样本集中比预先设定的心拍数量少的心拍,使心拍样本集中的心拍数量与预先设定的心拍数量相同,并将心拍样本集中的每个心拍按照所属心电形态类型进行数值标识。
具体实现时,由于心电数据长短不一,有4个心拍的,有10个心拍的,为了使输入数据的大小一致,需要将每个导联的心拍设置为预先设定的心拍数量进行聚类类型标记,例如,预先设定是心拍数量为10个,首先将每个导联的心拍样本集中的心拍数量和这10个预先设定的心拍数量进行比较,如果数量相等,对每个心拍进行所属聚类类型的数值标识。比如,正常心拍对应是数值标识为1,异常心拍1对应是数值标识为2,异常心拍2对应是数值标识为3,例如,导联V1的心拍样本集中的第1个心拍、第3个心拍、第4个心拍、第6个心拍以及第9个心拍聚类为正常心拍类型,心拍样本集中的第2个心拍和第7个心拍聚类为异常心拍1类型,心拍样本集中的第5个心拍、第8个心拍以及第10个心拍聚类为异常心拍2类型,按照心电数据的R峰位置序列顺序对每个心拍进行聚类类型数值标识,标识结果为:kind_label(V1)=[1,2,1,1,3,1,2,3,1,3],在本发明实施例中需要识别出病症类别为:正常心电、心律失常、室内传导阻滞、心室肥大、ST段异常、心肌梗死、心电轴偏移总共7个类,对上述病症依次编码为:1,2,3,4,5,6,7;由于聚类只能聚出3个类,使得聚类后的编码标签与病症标签无法对应,在具体实施过程中,对数据进行聚类后,通过相关领域专家对聚类后单个心拍进行标注,使聚类后的类别与实际病症类别相对应,例如,经专家标注后可知:类1对应病症1,类2对应病症5,类3对应病症7,因此,最终聚类后标签应为:[1,5,1,1,7,1,5,7,1,7]。其它导联的心拍标识重复以上操作得到kind_label(i),其中,kind_label(i)表示导联i中每个心拍所属聚类类型的标识数值,在此不再赘述。如果心拍样本集中的心拍数量大于预先设定的心拍数量,比如,心拍样本集中的心拍数量为12个多余预先设定的心拍数量2个,这时,按照心电数据的R峰位置序列顺序去掉最后两个心拍,在对剩余其它心拍根据上述标识方法对心拍样本集中的心拍进行数值标识。如果心拍样本集中的心拍数量小于预先设定的心拍数量,比如,心拍样本集中的心拍数量为6个小余预先设定的心拍数量4个,这时,按照心电数据的R峰位置序列顺序将心拍样本集中前4个心拍不到后4个,使得心拍样本集中的心拍数量也为10个,在对这10个心拍进行数值标识。
通常,对聚类后的心拍样本集进行重编码,生成编码心拍样本集的步骤包括:将每个心电形态类型包含的心拍数据进行求和平均计算,以获取每个心电形态类型包含的新的心拍数据;将每个心电形态类型包含的新的心拍数据以及心拍样本集中的每个心拍的标识数值进行融合,生成编码心拍样本集。
具体实现时,正常心拍类型中包括V1导联中的第1个心拍、第3个心拍、第4个心拍、第6个心拍以及第9个心拍这5个心拍,设置这5个心拍包含的数据长度N为5个,也就是每个心拍包含5个数据,组成一个5*5的数据矩阵,其中,每一行向代表每个心拍包含5个数据,每一纵向代表这5个心拍对应位置的数据,将每一纵向的数据进行求和平均,得到一个1*5的数据矩阵,同理,分别可以得到异常心拍1类型和异常心拍2类型的一个1*5的数据矩阵,具体地,导联i中各个心拍经过聚类后的结果为:kind_sig(i)=[sigi1(1:N),sigi2(1:N),sigi3(1:N)]其中:kind_sig(i)为导联i聚类后3种心拍类型数据,数据长度为N,也就是上面举例中的3个1*5的数据矩阵数据。将导联i聚类后3种心拍类型数据与导联i心拍样本集中的每个心拍的标识数值进行融合,导联i在重新编码成新数据的形式为decode_sig(i):decode_sig(i)=[kind_sig(i),kind_label(i)],最后将所有导联重新编码后生成编码心拍样本集的形式为decode_data:
Figure BDA0001904871550000131
利用上述方法对心电数据进行处理带来的有益效果为:去除了12导联数据中的冗余数据,突出了心电数据的重要关注点,为后面的模型训练做了充分是数据准备,不仅大大减少了模型训练时的数据量,而且极大的降低了训练时间,提高心电数据分类的准确度。
实施例二:
在上述实施例的基础上,本发明实施例还提供了一种心电数据的处理装置,如图3所示的一种心电数据的处理装置的结构示意图,该装置包括:
获取模块302,用于获取已知类别的心电数据库,其中,心电数据库包括多种类型的心电数据样本,每种类型的心电数据样本包括多个心电数据,每个心电数据携带有诊断标记;
提取模块304,用于提取每种类型中每个心电数据的心拍数据,按照预先设置的数据长度对每个心拍数据进行插值和采样操作,生成每个类型对应的心拍样本集,其中,心拍样本集中包含多个长度一致的心拍数据;
聚类模块306,用于按照预先设定的心电形态类型对心拍样本集中包括的多个心拍数据进行聚类,生成聚类结果,其中,聚类结果包含心电形态类型,以及,每个心电形态类型包含的心拍数据;
编码模块308,用于对聚类后的心拍样本集进行重编码,生成编码心拍样本集,以完成心电数据的处理过程。
在图3的基础上,图4示出了另一种心电数据的处理装置的结构示意图,该装置还包括:
训练模块402,用于利用编码心拍样本集对神经网络模型进行训练,以生成用于对待测心电数据进行分类的心电信号分类模型。
进一步,训练模型402用于:设置神经网络模型的输入层、隐藏层和输出层;其中,输出层的的长度至少包括以下类型的输出层:正常心电、心律失常、室内传导阻滞、心室肥大、ST段异常、心肌梗死、心电轴偏移。
本发明实施例提供的心电数据的处理装置,与上述实施例提供的心电数据的处理方法具有相同的技术特征,所以也能解决相同的技术问题,达到相同的技术效果。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的装置的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
另外,在本发明实施例的描述中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
最后应说明的是:以上实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。

Claims (9)

1.一种心电数据的处理方法,其特征在于,所述方法包括:
获取已知类别的心电数据库,其中,所述心电数据库包括多种类型的心电数据样本,每种所述类型的心电数据样本包括多个心电数据,每个所述心电数据携带有诊断标记;
提取每种所述类型中每个所述心电数据的心拍数据,按照预先设置的数据长度对每个所述心拍数据进行插值和采样操作,生成每个所述类型对应的心拍样本集,其中,所述心拍样本集中包含多个长度一致的心拍数据;
按照预先设定的心电形态类型对所述心拍样本集中包括的多个心拍数据进行聚类,生成聚类结果,其中,所述聚类结果包含所述心电形态类型,以及,每个所述心电形态类型包含的心拍数据;
对聚类后的所述心拍样本集进行重编码,生成编码心拍样本集,以完成所述心电数据的处理过程;
其中,所述对聚类后的所述心拍样本集进行重编码,生成编码心拍样本集的步骤包括:
将每个所述心电形态类型包含的心拍数据进行求和平均计算,以获取每个所述心电形态类型包含的新的心拍数据;
将每个所述心电形态类型包含的新的心拍数据以及所述心拍样本集中的每个心拍的标识数值进行融合,生成所述编码心拍样本集。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
利用所述编码心拍样本集对神经网络模型进行训练,以生成用于对待测心电数据进行分类的心电信号分类模型。
3.根据权利要求2所述的方法,其特征在于,所述利用所述编码心拍样本集对神经网络模型进行训练的步骤包括:
设置所述神经网络模型的输入层、隐藏层和输出层;
其中,所述输出层的的长度至少包括以下类型的输出层:正常心电、心律失常、心室肥大、ST段异常、心肌梗死、心电轴偏移;所述心律失常包括室内传导阻滞。
4.根据权利要求1所述的方法,其特征在于,所述提取每种所述类型中每个所述心电数据的心拍数据的步骤包括:
对获取的所述心电数据样本中包括的多个所述心电数据进行滤波处理;
查找滤波处理后的所述心电数据中的R波峰值点,以获取所述心电数据的R峰位置序列;
提取所述R峰位置序列中相邻两个所述R波峰值点,将相邻两个所述R波峰值点分别作为单个心拍的首尾点,相邻两个所述R波峰值点之间的心电数据作为一个周期的心拍,生成每个所述心电数据的心拍数据。
5.根据权利要求1所述的方法,其特征在于,所述按照预先设置的数据长度对每个所述心拍数据进行插值和采样操作的步骤包括:
对每个所述心拍数据进行插值处理;
利用预先设置的数据长度和插值处理后所述心拍数据计算采样间隔;
根据所述采样间隔对插值处理后的所述心拍数据进行采样处理,以生成长度一致的心拍数据。
6.根据权利要求1所述的方法,其特征在于,所述按照预先设定的心电形态类型对所述心拍样本集中包括的多个心拍数据进行聚类之后,所述方法还包括:
将所述心拍样本集中的心拍数量与预先设定的心拍数量进行比较;
如果所述心拍样本集中的心拍数量等于所述预先设定的心拍数量,将所述心拍样本集中的每个心拍按照所属心电形态类型进行数值标识;
如果所述心拍样本集中的心拍数量大于所述预先设定的心拍数量,按照所述心电数据的R峰位置序列顺序去掉所述心拍样本集中比所述预先设定的心拍数量多的心拍,使所述心拍样本集中的心拍数量与所述预先设定的心拍数量相同,并将所述心拍样本集中的每个心拍按照所属心电形态类型进行数值标识;
如果所述心拍样本集中的心拍数量小于所述预先设定的心拍数量,按照所述心电数据的R峰位置序列顺序扩大所述心拍样本集中比所述预先设定的心拍数量少的心拍,使所述心拍样本集中的心拍数量与所述预先设定的心拍数量相同,并将所述心拍样本集中的每个心拍按照所属心电形态类型进行数值标识。
7.一种心电数据的处理装置,其特征在于,所述装置包括:
获取模块,用于获取已知类别的心电数据库,其中,所述心电数据库包括多种类型的心电数据样本,每种所述类型的心电数据样本包括多个心电数据,每个所述心电数据携带有诊断标记;
提取模块,用于提取每种所述类型中每个所述心电数据的心拍数据,按照预先设置的数据长度对每个所述心拍数据进行插值和采样操作,生成每个所述类型对应的心拍样本集,其中,所述心拍样本集中包含多个长度一致的心拍数据;
聚类模块,用于按照预先设定的心电形态类型对所述心拍样本集中包括的多个心拍数据进行聚类,生成聚类结果,其中,所述聚类结果包含所述心电形态类型,以及,每个所述心电形态类型包含的心拍数据;
编码模块,用于对聚类后的所述心拍样本集进行重编码,生成编码心拍样本集,以完成所述心电数据的处理过程;
其中,编码模块,还用于,将每个所述心电形态类型包含的心拍数据进行求和平均计算,以获取每个所述心电形态类型包含的新的心拍数据;
将每个所述心电形态类型包含的新的心拍数据以及所述心拍样本集中的每个心拍的标识数值进行融合,生成所述编码心拍样本集。
8.根据权利要求7所述的装置,其特征在于,所述装置还包括:
训练模块,用于利用所述编码心拍样本集对神经网络模型进行训练,以生成用于对待测心电数据进行分类的心电信号分类模型。
9.根据权利要求8所述的装置,其特征在于,所述训练模型用于:
设置所述神经网络模型的输入层、隐藏层和输出层;
其中,所述输出层的的长度至少包括以下类型的输出层:正常心电、心律失常、心室肥大、ST段异常、心肌梗死、心电轴偏移;所述心律失常包括室内传导阻滞。
CN201811530312.8A 2018-12-13 2018-12-13 心电数据的处理方法及装置 Active CN109480825B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811530312.8A CN109480825B (zh) 2018-12-13 2018-12-13 心电数据的处理方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811530312.8A CN109480825B (zh) 2018-12-13 2018-12-13 心电数据的处理方法及装置

Publications (2)

Publication Number Publication Date
CN109480825A CN109480825A (zh) 2019-03-19
CN109480825B true CN109480825B (zh) 2021-08-06

Family

ID=65710188

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811530312.8A Active CN109480825B (zh) 2018-12-13 2018-12-13 心电数据的处理方法及装置

Country Status (1)

Country Link
CN (1) CN109480825B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111723622B (zh) * 2019-03-22 2024-04-26 安徽华米信息科技有限公司 心拍分类方法、装置、可穿戴设备及存储介质
CN109875548B (zh) * 2019-03-24 2022-04-19 浙江好络维医疗技术有限公司 一种基于多导联综合分析的心电图特征波形聚类方法
CN110288028B (zh) * 2019-06-27 2021-11-02 北京邮电大学 心电检测方法、系统、设备及计算机可读存储介质
CN110720903B (zh) * 2019-10-22 2022-08-23 深圳旭宏医疗科技有限公司 心电信号处理方法及装置
CN111488793B (zh) * 2020-02-21 2024-04-30 广州视源电子科技股份有限公司 心电信号分类方法、装置、电子设备和存储介质
CN111738348B (zh) * 2020-06-29 2023-07-11 南方电网科学研究院有限责任公司 一种电力数据异常检测方法和装置
CN111956201B (zh) * 2020-07-22 2022-09-06 上海数创医疗科技有限公司 基于卷积神经网络的心拍类型识别方法和装置
CN112528783B (zh) * 2020-11-30 2024-04-16 深圳邦健生物医疗设备股份有限公司 心电图心搏数据聚类方法、装置、电子设备和介质
CN112716504B (zh) * 2020-12-22 2023-12-15 沈阳东软智能医疗科技研究院有限公司 心电数据的处理方法、装置、存储介质和电子设备
CN113647959B (zh) * 2021-07-27 2024-01-12 东软集团股份有限公司 一种心电波形信号的波形识别方法、装置及设备
CN114795137A (zh) * 2022-03-08 2022-07-29 哈尔滨理工大学 一种基于神经网络的患者脉象诊断系统

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102038497A (zh) * 2010-12-02 2011-05-04 广东宝莱特医用科技股份有限公司 一种心电信号噪声分析方法
CN103815897A (zh) * 2014-02-28 2014-05-28 吉林大学 一种心电图特征提取方法
CN104523266A (zh) * 2015-01-07 2015-04-22 河北大学 一种心电信号自动分类方法
CN105748063A (zh) * 2016-04-25 2016-07-13 山东大学齐鲁医院 基于多导联和卷积神经网络的心律失常智能诊断方法
CN105962926A (zh) * 2016-04-18 2016-09-28 中国人民解放军总医院 基于心电图的心拍特征分析方法和设备
CN106214145A (zh) * 2016-07-20 2016-12-14 杨平 一种基于深度学习算法的心电图分类方法
CN106805965A (zh) * 2016-12-19 2017-06-09 深圳先进技术研究院 一种心电信号分类方法及装置
CN108464827A (zh) * 2018-03-08 2018-08-31 四川大学 一种弱监督下的心电图像识别方法
CN108523869A (zh) * 2017-02-16 2018-09-14 塔塔咨询服务有限公司 用于使用融合方法来检测人中的冠状动脉疾病的方法和系统
CN108647565A (zh) * 2018-03-28 2018-10-12 浙江工业大学 一种基于深度学习模型对心电信号分类的数据预处理方法
CN108647614A (zh) * 2018-04-28 2018-10-12 吉林大学 心电图心拍分类识别方法及系统
CN108875575A (zh) * 2018-05-12 2018-11-23 鲁东大学 基于多通道信号融合深度神经网络自动心律失常分析方法
CN108937915A (zh) * 2018-07-24 2018-12-07 杭州质子科技有限公司 一种动态心电图中早搏的识别方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102688032A (zh) * 2012-05-11 2012-09-26 东华大学 基于精简公共字典库的心电信号稀疏分解压缩系统
US20160106378A1 (en) * 2014-10-21 2016-04-21 Xerox Corporation System and method for detecting an arrhythmic cardiac event from a cardiac signal
US20170039358A1 (en) * 2015-08-07 2017-02-09 Fitbit, Inc. Transaction prevention using fitness data
US10595788B2 (en) * 2016-05-31 2020-03-24 Stmicroelectronics S.R.L. Method for the detecting electrocardiogram anomalies and corresponding system
CN107456227B (zh) * 2017-08-16 2021-03-30 北京蓬阳丰业医疗设备有限公司 全导联心电图聚类模板系统及方法
CN107625519B (zh) * 2017-09-20 2020-09-08 武汉中旗生物医疗电子有限公司 心电图处理方法及装置
CN108968951B (zh) * 2018-08-15 2021-06-22 武汉中旗生物医疗电子有限公司 心电图检测方法、装置及系统

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102038497A (zh) * 2010-12-02 2011-05-04 广东宝莱特医用科技股份有限公司 一种心电信号噪声分析方法
CN103815897A (zh) * 2014-02-28 2014-05-28 吉林大学 一种心电图特征提取方法
CN104523266A (zh) * 2015-01-07 2015-04-22 河北大学 一种心电信号自动分类方法
CN105962926A (zh) * 2016-04-18 2016-09-28 中国人民解放军总医院 基于心电图的心拍特征分析方法和设备
CN105748063A (zh) * 2016-04-25 2016-07-13 山东大学齐鲁医院 基于多导联和卷积神经网络的心律失常智能诊断方法
CN106214145A (zh) * 2016-07-20 2016-12-14 杨平 一种基于深度学习算法的心电图分类方法
CN106805965A (zh) * 2016-12-19 2017-06-09 深圳先进技术研究院 一种心电信号分类方法及装置
CN108523869A (zh) * 2017-02-16 2018-09-14 塔塔咨询服务有限公司 用于使用融合方法来检测人中的冠状动脉疾病的方法和系统
CN108464827A (zh) * 2018-03-08 2018-08-31 四川大学 一种弱监督下的心电图像识别方法
CN108647565A (zh) * 2018-03-28 2018-10-12 浙江工业大学 一种基于深度学习模型对心电信号分类的数据预处理方法
CN108647614A (zh) * 2018-04-28 2018-10-12 吉林大学 心电图心拍分类识别方法及系统
CN108875575A (zh) * 2018-05-12 2018-11-23 鲁东大学 基于多通道信号融合深度神经网络自动心律失常分析方法
CN108937915A (zh) * 2018-07-24 2018-12-07 杭州质子科技有限公司 一种动态心电图中早搏的识别方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Özbay Y, 等.A fuzzy clustering neural network architecture for classification of ECG arrhythmias.《Computers in Biology and Medicine》.2006,全文. *
Patient-specific ECG classification based on recurrent neural networks and clustering technique;Zhang C,等;《2017 13th IASTED International Conference on Biomedical Engineering》;20170228;全文 *
基于卷积神经网络的ECG信号个性化分类方法;王海龙,等;《中国自动化学会、济南市人民政府:中国自动化学会》;20171020;全文 *
基于改进K均值聚类生成匹配模板的心搏分类方法;陈永波,等;《传感器与微系统》;20180430;第37卷(第4期);全文 *
心电图的诊断算法研究;张亚梅;《电子测试》;20180430;全文 *

Also Published As

Publication number Publication date
CN109480825A (zh) 2019-03-19

Similar Documents

Publication Publication Date Title
CN109480825B (zh) 心电数据的处理方法及装置
CN106725428B (zh) 一种心电信号分类方法及装置
CN110890155B (zh) 一种基于导联注意力机制的多类心律失常检测方法
US11344243B2 (en) Artificial intelligence self-learning-based static electrocardiography analysis method and apparatus
CN107822622B (zh) 基于深度卷积神经网络的心电图诊断方法和系统
CN108836302B (zh) 基于深度神经网络的心电图智能分析方法及系统
Park et al. Hierarchical support vector machine based heartbeat classification using higher order statistics and hermite basis function
CN109303559B (zh) 一种基于梯度提升决策树的动态心电图心拍分类方法
JP2020536693A (ja) 人工知能自己学習に基づく心電図自動解析方法及び装置
CN106815570B (zh) 一种基于动态模式识别的心电信号st-t段识别方法
CN110522444B (zh) 一种基于Kernel-CNN的心电信号识别分类方法
Emanet ECG beat classification by using discrete wavelet transform and Random Forest algorithm
CN115486855B (zh) 一种基于qrs波群不定次循环叠加的心电图心搏分类方法
CN111631705A (zh) 心电异常检测方法、模型训练方法、装置、设备及介质
CN109480827B (zh) 向量心电图分类方法及装置
CN111743531A (zh) 心电节律和形态异常识别及装置
JP7487965B2 (ja) グラフ畳み込みに基づく心電図心拍数マルチタイプの予測方法
CN111557660A (zh) 一种亚群体深度学习框架下的心律失常识别方法
Venton et al. Signal to image to classification: Transfer learning for ecg
Besrour et al. ECG beat classifier using support vector machine
CN110786847B (zh) 心电信号的建库方法和分析方法
CN113812962B (zh) 心电节律分类神经网络构建方法、系统及存储介质
CN114224351B (zh) 一种基于融合多个深度学习模型的房颤识别方法
CN111803062B (zh) 一种基于深度学习的房颤事件检测方法
Bashir et al. Highlighting the current issues with pride suggestions for improving the performance of real time cardiac health monitoring

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
GR01 Patent grant
GR01 Patent grant