CN106214123A - 一种基于深度学习算法的心电图综合分类方法 - Google Patents

一种基于深度学习算法的心电图综合分类方法 Download PDF

Info

Publication number
CN106214123A
CN106214123A CN201610572208.XA CN201610572208A CN106214123A CN 106214123 A CN106214123 A CN 106214123A CN 201610572208 A CN201610572208 A CN 201610572208A CN 106214123 A CN106214123 A CN 106214123A
Authority
CN
China
Prior art keywords
electrocardiogram
data
pqrst
wave
degree
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
CN201610572208.XA
Other languages
English (en)
Other versions
CN106214123B (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201610572208.XA priority Critical patent/CN106214123B/zh
Publication of CN106214123A publication Critical patent/CN106214123A/zh
Application granted granted Critical
Publication of CN106214123B publication Critical patent/CN106214123B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Landscapes

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

Abstract

本发明公开了一种基于深度学习算法的心电图综合分类方法,包括以下步骤:取得原始心电图波形数据、心电图附加信息,并获取心电图节律信息、代表PQRST波形数据;将相关信息通过训练完的第一深度学习算法进行波形分类,得到分类结果一;将相关信息通过训练完的第二深度学习算法后得到P波、QRS波、T波类型数据并计算出代表PQRST波形特征数据后输入传统心电图计算机自动分类算法,得到分类结果二;加入权重来调整分类结果评分并取评分值最大的分类结果作为最终分类结果。本发明合理地结合心电图分类的特点,并通过以上步骤对深度学习方法进行训练并用深度学习方法进行波形分类,提高了心电图解释得到的分类结果的正确率。

Description

一种基于深度学习算法的心电图综合分类方法
技术领域
本发明涉及心电图分类方法,尤其涉及一种基于深度学习算法的心电图综合分类方法。
背景技术
心电图波形数据采集及心电图分类结果是医生诊断心脏病疾患的重要辅助手段和参考信息,通常心电图波形数据采集和分类是在医院或体检中心进行,存在检测不方便、检测频率低等缺点,而且不能及时地把心电图分类信息提供给医生来做实时诊断,很难有效地预防和及时治疗心脏病病变。近年来,随着网络、移动智能手机的普及,使得携带式心电监测仪、家庭个人用心电波监测仪的推出成为可能。目前市场上推出的这类监测仪,是基于传统的测量分类法对采集到的心电图波形数据进行检测,而携带式心电监测仪、家庭个人用心电波监测仪的信号容易受外界噪音干扰,导致传统的测量分类法在对一些切分点不明显的波形进行分类时,往往会出现误分类,从而导致计算参数不正确,并最终导致误分类,其临床可靠性和准确性较低,无法满足实际的给医生提供辅助诊断信息的需要。
发明内容
本发明针对现有技术中存在的传统的心电图测量分类法在单一波形测量特征不太明显时容易出现误分类,其临床可靠性和准确性较低,无法满足实际的分类需要等缺陷,提供了一种新的基于深度学习算法的心电图综合分类方法。
为了解决上述技术问题,本发明通过以下技术方案实现:
一种基于深度学习算法的心电图综合分类方法,包括以下步骤:
a.取得测量时间在8秒以上的原始心电图波形数据、心电图附加信息,并根据原始心电图波形数据进行心电图节律信息的提取、代表PQRST波形的提取,获取心电图节律信息、代表PQRST波形数据;
b11.对第一深度学习算法的神经网络进行训练;
b12.将步骤a得到的心电图节律信息、代表PQRST波形数据、心电图附加信息排列成一维数据后通过训练完的第一深度学习算法进行波形分类,得到分类结果一;
b21.对第二深度学习算法的神经网络进行训练;
b22.将步骤a得到的代表PQRST波形数据从训练完的第二深度学习算法的输入端输入,经第二深度学习算法进行PQRST波群分类后,从第二深度学习算法的输出端得到P波类型数据、QRS波类型数据、T波类型数据;
b23.根据步骤b22得到的P波类型数据、QRS波类型数据、T波类型数据分析步骤a得到的代表PQRST波形数据,计算出代表PQRST波形特征数据;
b24.将步骤a得到的心电图附加信息、心电图节律信息以及步骤b23得到的代表PQRST波形特征数据输入传统心电图计算机自动分类算法,得到分类结果二;
c.对步骤b12得到的分类结果一加入权重W1来调整分类结果评分,对步骤b24得到的分类结果二加入权重W2来调整分类结果评分,所述的权重W1、权重W2满足条件W1+W2=1;
d.从步骤c的分类结果评分中取评分值最大的分类结果作为最终分类结果。
深度学习算法是一种人工智能领域的机器学习方法,它含有多隐层的多层感知器,是通过组合低层特征形成更加抽象的高层表示属性类别或特征,以发现数据的分布式特征表示,深度学习方法目前已经在图像识别,声音识别等应用领域里被证明了它的有效性,能够大幅度提高了传统方法的识别精度。本发明中把深度学习方法导入心电图分类领域,合理地结合心电图分类的特点,并通过以上步骤对深度学习方法进行训练并用深度学习方法进行波形分类,能够大幅度提高心电图解释得到的分类结果的正确率。
其中8秒以上的原始心电图波形数据具有数量充足的波形,使得提取得到的心电图节律信息、代表PQRST波形数据更加准确。其中代表PQRST波形的提取能够有效减小非分类要素如人体移动、电极不稳定带来的波形变化影响,同时由于代表PQRST波形的数据量相对原始心电图波形数据中的波形数据少很多,其波形数据更加稳定,可以大幅度减少后期深度学习算法的训练量,提高深度学习算法的计算效率,并提高分类准确度。心电图节律信息的提取能够用来提高相关的心电图分类信息的精确度。心电图附加信息会影响原始心电图波形数据中心电波振幅、噪音构成等特性,对心电图附加信息的使用能够进一步提升最终得到的心电图分类信息的精确度。
步骤b11至步骤b12能够将心电图节律信息、代表PQRST波形数据、心电图附加信息排列成一维数据后通过深度学习算法进行波形分类,能够让深度学习算法分析这些信息之间的关联性,从而能够得到分类正确率较高的分类结果一。步骤b21至步骤b24能够克服传统波形测量分类方法里存在的最大弱点,即因为P波、QRS波、T波类型的错误分类造成的二次参数的误算导致的最终分类错误。P波、QRS波、T波类型的错误分类是由于在传统方法里对于每个导联的代表波形进行分类时,对每个导联的代表波形间的相关关系,每个导联代表波形里的P波、QRS波、T波间的相关关系,以及每个导联的时间要素和电压要素间的相关关系,以及心电图附加信息、心电图节律信息和上述关系间的关系的分析能力不够造成的。而步骤b21至步骤b24通过深度学习算法对大量心电图波形数据进行学习后,能够自动总结出这些相关关系,从而能够提升P波、QRS波、T波的分类正确度,最终大幅度提高分类结果正确率,同时步骤b21至步骤b24保留了已被医生广为接受的代表PQRST波形特征数据和传统心电图计算机自动分类算法的各个分析过程,使得在提高最终分类正确率的同时又使得深度学习算法的引入更容易为医生们接受。
而步骤c至步骤d通过加入权重W1、权重W2来分别调整分类结果一、分类结果二,使得最终分类结果不是单一地以分类结果一或分类结果二为准,而是综合考量了分类结果一、分类结果二后获得的,从而在分类结果一、分类结果二正确率的基础上能够进一步提高最终分类结果的正确率。
作为优选,上述所述的一种基于深度学习算法的心电图综合分类方法,所述的步骤c中的权重W1、权重W2通过以下步骤得到:
c11.准备一组由n人次已知真实分类结果的心电图波形数据构成的数据库;
c12.对步骤c11中数据库里的某个心电图波形数据通过步骤b12中训练完的第一深度学习算法进行波形分类,得到分类结果评价值S1,对该心电图波形数据通过步骤b22中训练完的第二深度学习算法进行PQRST波群分类后并经步骤b23、步骤b24,得到分类结果评价值S2;
c13.在满足W1+W2=1的条件下通过公式S1×W1+S2×W2=S3调节权重W1和权重W2的值,直到综合结果评价值S3符合已知真实分类结果取值范围,得到权重W1和权重W2的值;
c14.对步骤c11中数据库里的所有心电图波形数据重复步骤c12和步骤c13;
c15.重复步骤c14直至权重W1和权重W2的值趋于稳定,得到最终的权重W1和权重W2的值。
通过以上步骤能够在确切的真实分类结果的基础上通过深度学习算法自我学习自我更新,从而得到合理的最终的权重W1和权重W2的值,从而进一步提高最终分类结果的正确率。
作为优选,上述所述的一种基于深度学习算法的心电图综合分类方法,所述的步骤a中,所述的心电图节律信息包括平均心室心率、平均RR间隔、最长RR间隔和最短RR间隔之差、RR间隔的标准方差、一致性P波信息、窦房律下各心拍的PR间隔和平均值、R波中预激波检测结果、窦房律下各心拍的QT间期和QTc间期和平均值、窦房律QRS平均波宽、窦房律P波宽和平均波宽、期外收缩信息、期外收缩类型、期外收缩形态、心房扑动的F波和心房颤动的f波的检测结果、非同期P波的检测结果。
以上信息对最终的分类结果具有较大影响,能够进一步提升最终的分类结果的正确率。
作为优选,上述所述的一种基于深度学习算法的心电图综合分类方法,所述的步骤a中,代表PQRST波形的提取包括以下步骤:
a11.通过一阶微分法和阈值法,对原始心电图波形数据进行检测,获取P波、QRS波、T波的特征点;
a12.对原始心电图波形数据中的所有PQRST波进行聚类分析,根据分析结果将拥有PQRST波数目最多的类型作为代表PQRST波形,若数目最多的类型为2个以上,选取R波平均振幅最大的类型作为代表PQRST波形,最后使用叠加平均法计算各个心拍的PQRST波的平均波形作为代表PQRST波形。
通过以上步骤,能够有效提取原始心电图波形数据中的P波、QRS波、T波的特征点,而将原始心电图波形数据中的所有PQRST波进行聚类分析,能够有效去除原始心电图波形数据中受到噪音伪差干扰的PQRST波形以及与节律相关的QRST波形,保证得到的代表PQRST波形能够传递更准确的有效信息进行心电图分类。
作为优选,上述所述的一种基于深度学习算法的心电图综合分类方法,所述的步骤a中,所述的心电图附加信息包括性别、身高、胸围、体重、脂肪率、人种。
以上信息与心电图判读基准有关,对最终的分类结果有较大影响,考虑以上信息能提升分类结果的正确率。
作为优选,上述所述的一种基于深度学习算法的心电图综合分类方法,所述的步骤b23中,所述的代表PQRST波形特征数据包括代表PQRST波形数据中每个波形的切分点、波峰电压值以及心电图节律信息中所有心拍的节律信息。具体的,代表PQRST波形特征数据包括:P波宽、P’宽、P波高、P波类型、PQ间期、QRS波宽、Q波宽、Q波振幅、R波高、R波宽、R’波高、S波振幅、S’波宽、QRS波类型、ST起点高、ST中段高、ST终点高、QT间期、T波高、修正T波高、delta波有无、U波高、P/QRS/T波电轴以及心电图节律信息中所有心拍的节律信息等信息。
以上信息是对P波、QRS波、T波分类的基础数据,能够较大程度上影响P波、QRS波、T波分类结果,能够进一步提升分类的正确率,以及最终的分类结果正确率。
作为优选,上述所述的一种基于深度学习算法的心电图综合分类方法,所述的步骤a中的原始心电波数据为单导联数据。
单导联数据一般适用便携式的心脏病检测仪器,使得本发明的适用范围更广。
作为优选,上述所述的一种基于深度学习算法的心电图综合分类方法,所述的步骤a中的原始心电图波形数据为多导联数据,所述的心电图节律信息由各个导联的心电图节律信息串联成一维数据形成,所述的代表PQRST波形数据由各个导联的代表PQRST波形数据串联成一维数据形成。
多导联的原始心电图波形数据具有更加充分的信息,能提升相关的心电图分类辅助信息的精确度,而由各个导联的心电图节律信息串联成的心电图节律信息以及由各个导联的代表PQRST波形数据串联成的代表PQRST波形数据在通过深度学习算法进行波形分类时,各个导联间的相关性在经过充分的训练后可以得到有效总结,能够进一步提升最终的分类结果的正确率。
作为优选,上述所述的一种基于深度学习算法的心电图综合分类方法,所述的第一深度学习算法、第二深度学习算法为卷积神经网络、迭代神经网络、深度神经网络中的任一种。
以上三种神经网络具有更高的准确率,能够保证最终的分类结果的正确率。
本发明的有益效果如下:
1、本发明将人工智能领域的深度学习算法和传统心电图分类算法进行了合理的结合,能够大幅度提高最终的心电图分类的正确率。本发明利用了在传统心电图分类方法里已经被证明的有效信息数据,同时又利用深度学习算法的超强学习能力、特征自动提取、特征分布关系自动分析等优势能力来弥补传统心电图分类方法里存在的特征提取不够精确、特征间相关性分析不充分的缺点,本发明通过深度学习算法能够自动学习大量的心电图波形数据,并总结出各心电图分类的敏感的特征及其分布,从而大幅度提高最终的心电图分类正确率。
2、本发明能够更有效地提供早期治疗时医生所需的心电图分类的辅助信息。传统心电图分类算法是一种静态算法,不具有自我学习能力,而本发明通过将传统心电图分类算法与深度学习算法相结合,一方面能提高心电图分类的正确率和鲁棒性,另一方面能提升对各类心电图分类的理解,提供医生解释各类心脏疾病的机理的辅助信息。
3、本发明引入深度学习算法进行P波、QRS波、T波分类,极大地提高了分类的正确度,进而极大地提高了代表PQRST波形特征数据的数据精度,使得代表PQRST波形特征数据更能反映出真实状况,从而极大地提高了代表PQRST波形特征数据通过传统心电图计算机自动分类算法分析得到的最终分类结果的准确度。
4、本发明在分类结果一、分类结果二的基础上,通过权重W1、权重W2的调整,使得最终分类结果能够进一步修正分类结果一、分类结果二之间的差异性,进一步提高最终分类结果的正确率。
附图说明
图1为本发明一种基于深度学习算法的心电图综合分类方法的流程图。
具体实施方式
下面结合附图1和具体实施方式对本发明作进一步详细描述,但它们不是对本发明的限制:
实施例1
一种基于深度学习算法的心电图综合分类方法,其流程图如图1所示,具体包括以下步骤:
a.(1)获取单导联心电图波形数据以及心电图附加信息,并根据该单导联心电图波形数据截取其中长度为10秒的数据作为原始心电图波形数据,其中单导联心电图波形数据以及心电图附加信息可以通过体检结果获得,也可以通过已有的数据库如欧共体心电图波形数据库(CSE)获得,或者通过其它途径获得,心电图附加信息中包括性别、身高、胸围、体重、脂肪率、人种。
(2)根据需要,可对步骤(1)得到的原始心电图波形数据进行去噪处理,去噪处理包括以下步骤:
i1.采用高通滤波器去除基线漂移噪音;
i2.基于PQ段信号的标准方差和阈值法确认噪音是否过高;
i3.噪音过高时使用低通巴特沃斯滤波器去除噪音干扰。
(3)根据原始心电图波形数据计算PQRST波形的区分点,从而根据PQRST波形的区分点进行心电图节律信息的提取,获取心电图节律信息,其中心电图节律信息包括平均心室心率、平均RR间隔、最长RR间隔和最短RR间隔之差、RR间隔的标准方差、一致性P波信息、窦房律下各心拍的PR间隔和平均值、R波中预激波检测结果、窦房律下各心拍的QT间期和QTc间期和平均值、窦房律QRS平均波宽、窦房律P波宽和平均波宽、期外收缩信息、期外收缩类型、期外收缩形态、心房扑动的F波和心房颤动的f波的检测结果、非同期P波的检测结果,而代表PQRST波形通过以下步骤提取:
a11.通过一阶微分法和阈值法,对原始心电图波形数据进行检测,获取P波、QRS波、T波的特征点;
a12.对原始心电图波形数据中的所有PQRST波进行聚类分析,根据分析结果将拥有PQRST波数目最多的类型作为代表PQRST波形,若数目最多的类型为2个以上,选取R波平均振幅最大的类型作为代表PQRST波形类,最后使用叠加平均法计算各个心拍的PQRST波的平均波形作为代表PQRST波形。
(4)为了对第一深度学习算法以及第二深度学习算法进行训练,还需要获取训练用数据,该训练用数据可从其它相应体检结果中获取,也可以选择从现有的数据库例如从欧共体心电图波形数据库(CSE)中获取,其中训练用数据包括其它相应单导联心电图波形数据以及心电图附加信息,以欧共体心电图波形数据库(CSE)为例,训练用数据的获取步骤为:根据步骤(1)-步骤(3)对欧共体心电图波形数据库(CSE)中各个信息进行处理,获取欧共体心电图波形数据库(CSE)中各个心电图节律信息、代表PQRST波形数据、心电图附加信息。
b11.(1)第一深度学习算法选择卷积神经网络,并设定卷积神经网络的输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据中的心电图节律信息、代表PQRST波形数据、心电图附加信息排列成一维数据后从卷积神经网络的输入端输入,并把对应的真实分类结果输入到卷积神经网络的结果端来训练卷积神经网络。
b12.在完成对卷积神经网络的训练后,将步骤a的步骤(1)中的心电图附加信息、步骤a的步骤(3)中代表PQRST波形数据以及心电图节律信息排列成一维数据后输入到卷积神经网络的输入端,即可通过卷积神经网络进行波形分类后得到分类结果一,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的排列方式可根据实际情况进行选择。
b21.(1)第二深度学习算法选择迭代神经网络,设定迭代神经网络输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据从迭代神经网络的输入端输入,经迭代神经网络进行PQRST波群分类后,从迭代神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
(3)将迭代神经网络输出端得到的P波类型数据、QRS波类型数据、T波类型数据与真实波形类型进行比较,基于输出端和真实波形类型的检测偏差,利用后向传递算法修改迭代神经网络内各节点的权值。
(4)重复步骤(2)-(3)直到迭代神经网络判读结果符合真实波形类型,即迭代神经网络结构参数收敛时,获得训练好的迭代神经网络。
b22.将步骤a中的步骤(3)得到的代表PQRST波形数据从步骤(4)训练好的迭代神经网络的输入端输入,经迭代神经网络进行PQRST波群分类后,从迭代神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
b23.使用步骤b22得到的P波类型数据、QRS波类型数据、T波类型数据分析步骤a的步骤(3)得到的代表PQRST波形数据、PQRST波形的各区分点信息,可以计算出代表PQRST波形特征数据。
b24.(1)将步骤a的步骤(1)得到的心电图附加信息、步骤a的步骤(3)得到的心电图节律信息以及步骤b23得到的代表PQRST波形特征数据组成输入数据,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的组合方式可根据实际情况进行选择。
(2)将步骤(1)得到的输入数据输入传统心电图计算机自动分类算法,得到自动分类结果二。
c.对步骤b12得到的分类结果一加入权重W1来调整分类结果评分,对步骤b24得到的分类结果二加入权重W2来调整分类结果评分,所述的权重W1、权重W2满足条件W1+W2=1。
d.从步骤c的分类结果评分中取评分值最大的分类结果作为最终分类结果。
实施例2
一种基于深度学习算法的心电图综合分类方法,其流程图如图1所示,具体包括以下步骤:
a.(1)获取单导联心电图波形数据以及心电图附加信息,并根据该单导联心电图波形数据截取其中长度为8秒的数据作为原始心电图波形数据,其中单导联心电图波形数据以及心电图附加信息可以通过体检结果获得,也可以通过已有的数据库如欧共体心电图波形数据库(CSE)获得,或者通过其它途径获得,心电图附加信息中包括性别、身高、胸围、体重、脂肪率、人种。
(2)根据需要,可对步骤(1)得到的原始心电图波形数据进行去噪处理,去噪处理包括以下步骤:
i1.采用高通滤波器去除基线漂移噪音;
i2.基于PQ段信号的标准方差和阈值法确认噪音是否过高;
i3.噪音过高时使用低通巴特沃斯滤波器去除噪音干扰。
(3)根据原始心电图波形数据计算PQRST波形的区分点,从而根据PQRST波形的区分点进行心电图节律信息的提取,获取心电图节律信息,其中心电图节律信息包括平均心室心率、平均RR间隔、最长RR间隔和最短RR间隔之差、RR间隔的标准方差、一致性P波信息、窦房律下各心拍的PR间隔和平均值、R波中预激波检测结果、窦房律下各心拍的QT间期和QTc间期和平均值、窦房律QRS平均波宽、窦房律P波宽和平均波宽、期外收缩信息、期外收缩类型、期外收缩形态、心房扑动的F波和心房颤动的f波的检测结果、非同期P波的检测结果,而代表PQRST波形通过以下步骤提取:
a11.通过一阶微分法和阈值法,对原始心电图波形数据进行检测,获取P波、QRS波、T波的特征点;
a12.对原始心电图波形数据中的所有PQRST波进行聚类分析,根据分析结果将拥有PQRST波数目最多的类型作为代表PQRST波形,若数目最多的类型为2个以上,选取R波平均振幅最大的类型作为代表PQRST波形类,最后使用叠加平均法计算各个心拍的PQRST波的平均波形作为代表PQRST波形。
(4)为了对第一深度学习算法以及第二深度学习算法进行训练,还需要获取训练用数据,该训练用数据可从其它相应体检结果中获取,也可以选择从现有的数据库例如从欧共体心电图波形数据库(CSE)中获取,其中训练用数据包括其它相应单导联心电图波形数据以及心电图附加信息,以欧共体心电图波形数据库(CSE)为例,训练用数据的获取步骤为:根据步骤(1)-步骤(3)对欧共体心电图波形数据库(CSE)中各个信息进行处理,获取欧共体心电图波形数据库(CSE)中各个心电图节律信息、代表PQRST波形数据、心电图附加信息。
b11.(1)第一深度学习算法选择迭代神经网络,并设定迭代神经网络的输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据中的心电图节律信息、代表PQRST波形数据、心电图附加信息排列成一维数据后从迭代神经网络的输入端输入,并把对应的真实分类结果输入到迭代神经网络的结果端来训练迭代神经网络。
b12.在完成对迭代神经网络的训练后,将步骤a的步骤(1)中的心电图附加信息、步骤a的步骤(3)中代表PQRST波形数据以及心电图节律信息排列成一维数据后输入到迭代神经网络的输入端,即可通过迭代神经网络进行波形分类后得到分类结果一,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的排列方式可根据实际情况进行选择。
b21.(1)第二深度学习算法选择卷积神经网络,设定卷积神经网络输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据从卷积神经网络的输入端输入,经卷积神经网络进行PQRST波群分类后,从卷积神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
(3)将卷积神经网络输出端得到的P波类型数据、QRS波类型数据、T波类型数据与真实波形类型进行比较,基于输出端和真实波形类型的检测偏差,利用后向传递算法修改卷积神经网络内各节点的权值。
(4)重复步骤(2)-(3)直到卷积神经网络判读结果符合真实波形类型,即卷积神经网络结构参数收敛时,获得训练好的卷积神经网络。
b22.将步骤a中的步骤(3)得到的代表PQRST波形数据从步骤(4)训练好的卷积神经网络的输入端输入,经卷积神经网络进行PQRST波群分类后,从卷积神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
b23.使用步骤b22得到的P波类型数据、QRS波类型数据、T波类型数据分析步骤a的步骤(3)得到的代表PQRST波形数据、PQRST波形的各区分点信息,可以计算出代表PQRST波形特征数据。
b24.(1)将步骤a的步骤(1)得到的心电图附加信息、步骤a的步骤(3)得到的心电图节律信息以及步骤b23得到的代表PQRST波形特征数据组成输入数据,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的组合方式可根据实际情况进行选择。
(2)将步骤(1)得到的输入数据输入传统心电图计算机自动分类算法,得到自动分类结果二。
c.对步骤b12得到的分类结果一加入权重W1来调整分类结果评分,对步骤b24得到的分类结果二加入权重W2来调整分类结果评分,所述的权重W1、权重W2满足条件W1+W2=1。
d.从步骤c的分类结果评分中取评分值最大的分类结果作为最终分类结果。
实施例3
一种基于深度学习算法的心电图综合分类方法,其流程图如图1所示,具体包括以下步骤:
a.(1)获取单导联心电图波形数据以及心电图附加信息,并根据该单导联心电图波形数据截取其中长度为16秒的数据作为原始心电图波形数据,其中单导联心电图波形数据以及心电图附加信息可以通过体检结果获得,也可以通过已有的数据库如欧共体心电图波形数据库(CSE)获得,或者通过其它途径获得,心电图附加信息中包括性别、身高、胸围、体重、脂肪率、人种。
(2)根据需要,可对步骤(1)得到的原始心电图波形数据进行去噪处理,去噪处理包括以下步骤:
i1.采用高通滤波器去除基线漂移噪音;
i2.基于PQ段信号的标准方差和阈值法确认噪音是否过高;
i3.噪音过高时使用低通巴特沃斯滤波器去除噪音干扰。
(3)根据原始心电图波形数据计算PQRST波形的区分点,从而根据PQRST波形的区分点进行心电图节律信息的提取,获取心电图节律信息,其中心电图节律信息包括平均心室心率、平均RR间隔、最长RR间隔和最短RR间隔之差、RR间隔的标准方差、一致性P波信息、窦房律下各心拍的PR间隔和平均值、R波中预激波检测结果、窦房律下各心拍的QT间期和QTc间期和平均值、窦房律QRS平均波宽、窦房律P波宽和平均波宽、期外收缩信息、期外收缩类型、期外收缩形态、心房扑动的F波和心房颤动的f波的检测结果、非同期P波的检测结果,而代表PQRST波形通过以下步骤提取:
a11.通过一阶微分法和阈值法,对原始心电图波形数据进行检测,获取P波、QRS波、T波的特征点;
a12.对原始心电图波形数据中的所有PQRST波进行聚类分析,根据分析结果将拥有PQRST波数目最多的类型作为代表PQRST波形,若数目最多的类型为2个以上,选取R波平均振幅最大的类型作为代表PQRST波形类,最后使用叠加平均法计算各个心拍的PQRST波的平均波形作为代表PQRST波形。
(4)为了对第一深度学习算法以及第二深度学习算法进行训练,还需要获取训练用数据,该训练用数据可从其它相应体检结果中获取,也可以选择从现有的数据库例如从欧共体心电图波形数据库(CSE)中获取,其中训练用数据包括其它相应单导联心电图波形数据以及心电图附加信息,以欧共体心电图波形数据库(CSE)为例,训练用数据的获取步骤为:根据步骤(1)-步骤(3)对欧共体心电图波形数据库(CSE)中各个信息进行处理,获取欧共体心电图波形数据库(CSE)中各个心电图节律信息、代表PQRST波形数据、心电图附加信息。
b11.(1)第一深度学习算法选择深度神经网络,并设定深度神经网络的输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据中的心电图节律信息、代表PQRST波形数据、心电图附加信息排列成一维数据后从深度神经网络的输入端输入,并把对应的真实分类结果输入到深度神经网络的结果端来训练深度神经网络。
b12.在完成对深度神经网络的训练后,将步骤a的步骤(1)中的心电图附加信息、步骤a的步骤(3)中代表PQRST波形数据以及心电图节律信息排列成一维数据后输入到深度神经网络的输入端,即可通过深度神经网络进行波形分类后得到分类结果一,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的排列方式可根据实际情况进行选择。
b21.(1)第二深度学习算法选择深度神经网络,设定深度神经网络输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据从深度神经网络的输入端输入,经深度神经网络进行PQRST波群分类后,从深度神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
(3)将深度神经网络输出端得到的P波类型数据、QRS波类型数据、T波类型数据与真实波形类型进行比较,基于输出端和真实波形类型的检测偏差,利用后向传递算法修改深度神经网络内各节点的权值。
(4)重复步骤(2)-(3)直到深度神经网络判读结果符合真实波形类型,即深度神经网络结构参数收敛时,获得训练好的深度神经网络。
b22.将步骤a中的步骤(3)得到的代表PQRST波形数据从步骤(4)训练好的深度神经网络的输入端输入,经深度神经网络进行PQRST波群分类后,从深度神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
b23.使用步骤b22得到的P波类型数据、QRS波类型数据、T波类型数据分析步骤a的步骤(3)得到的代表PQRST波形数据、PQRST波形的各区分点信息,可以计算出代表PQRST波形特征数据。
b24.(1)将步骤a的步骤(1)得到的心电图附加信息、步骤a的步骤(3)得到的心电图节律信息以及步骤b23得到的代表PQRST波形特征数据组成输入数据,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的组合方式可根据实际情况进行选择。
(2)将步骤(1)得到的输入数据输入传统心电图计算机自动分类算法,得到自动分类结果二。
c.对步骤b12得到的分类结果一加入权重W1来调整分类结果评分,对步骤b24得到的分类结果二加入权重W2来调整分类结果评分,所述的权重W1、权重W2满足条件W1+W2=1。
d.从步骤c的分类结果评分中取评分值最大的分类结果作为最终分类结果。
实施例4
一种基于深度学习算法的心电图综合分类方法,其流程图如图1所示,具体包括以下步骤:
a.(1)获取多导联心电图波形数据以及心电图附加信息,并根据该多导联心电图波形数据截取其中长度为10秒的数据作为原始心电图波形数据,其中多导联心电图波形数据以及心电图附加信息可以通过体检结果获得,也可以通过已有的数据库如欧共体心电图波形数据库(CSE)获得,或者通过其它途径获得,该多导联心电图波形数据可以为十二导联心电图波形数据、三导联心电图波形数据、六导联心电图波形数据、十八导联心电图波形数据等多导联心电图波形数据。心电图附加信息中包括性别、身高、胸围、体重、脂肪率、人种。
(2)根据需要,可对步骤(1)得到的原始心电图波形数据进行去噪处理,去噪处理包括以下步骤:
i1.采用高通滤波器去除基线漂移噪音;
i2.基于PQ段信号的标准方差和阈值法确认噪音是否过高;
i3.噪音过高时使用低通巴特沃斯滤波器去除噪音干扰。
(3)根据原始心电图波形数据计算PQRST波形的区分点,从而根据PQRST波形的区分点进行心电图节律信息的提取,获取心电图节律信息,其中心电图节律信息包括平均心室心率、平均RR间隔、最长RR间隔和最短RR间隔之差、RR间隔的标准方差、一致性P波信息、窦房律下各心拍的PR间隔和平均值、R波中预激波检测结果、窦房律下各心拍的QT间期和QTc间期和平均值、窦房律QRS平均波宽、窦房律P波宽和平均波宽、期外收缩信息、期外收缩类型、期外收缩形态、心房扑动的F波和心房颤动的f波的检测结果、非同期P波的检测结果,而代表PQRST波形通过以下步骤提取:
a11.通过一阶微分法和阈值法,对原始心电图波形数据进行检测,获取P波、QRS波、T波的特征点;
a12.对原始心电图波形数据中的所有PQRST波进行聚类分析,根据分析结果将拥有PQRST波数目最多的类型作为代表PQRST波形,若数目最多的类型为2个以上,选取R波平均振幅最大的类型作为代表PQRST波形类,最后使用叠加平均法计算各个心拍的PQRST波的平均波形作为代表PQRST波形。
(4)为了对第一深度学习算法以及第二深度学习算法进行训练,还需要获取训练用数据,该训练用数据可从其它相应体检结果中获取,也可以选择从现有的数据库例如从欧共体心电图波形数据库(CSE)中获取,其中训练用数据包括其它相应多导联心电图波形数据以及心电图附加信息,以欧共体心电图波形数据库(CSE)为例,训练用数据的获取步骤为:根据步骤a中的步骤(1)-步骤(3)对欧共体心电图波形数据库(CSE)中各个信息进行处理,获取欧共体心电图波形数据库(CSE)中各个心电图节律信息、代表PQRST波形数据、心电图附加信息并链接成一维输入数据,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的排列方式可根据实际情况进行选择。
b11.(1)第一深度学习算法选择深度神经网络,并设定深度神经网络的输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据中的心电图节律信息、代表PQRST波形数据、心电图附加信息排列成一维数据后从深度神经网络的输入端输入,并把对应的真实分类结果输入到深度神经网络的结果端来训练深度神经网络。
b12.在完成对深度神经网络的训练后,将步骤a的步骤(1)中的心电图附加信息、步骤a的步骤(3)中代表PQRST波形数据以及心电图节律信息排列成一维数据后输入到深度神经网络的输入端,即可通过深度神经网络进行波形分类后得到分类结果一,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的排列方式可根据实际情况进行选择。
b21.(1)第二深度学习算法选择迭代神经网络,设定迭代神经网络输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据从迭代神经网络的输入端输入,经迭代神经网络进行PQRST波群分类后,从迭代神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
(3)将迭代神经网络输出端得到的P波类型数据、QRS波类型数据、T波类型数据与真实波形类型进行比较,基于输出端和真实波形类型的检测偏差,利用后向传递算法修改迭代神经网络内各节点的权值。
(4)重复步骤(2)-(3)直到迭代神经网络判读结果符合真实波形类型,即迭代神经网络结构参数收敛时,获得训练好的迭代神经网络。
b22.将步骤a中的步骤(3)得到的代表PQRST波形数据从步骤(4)训练好的迭代神经网络的输入端输入,经迭代神经网络进行PQRST波群分类后,从迭代神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
b23.使用步骤b22得到的P波类型数据、QRS波类型数据、T波类型数据分析步骤a的步骤(3)得到的代表PQRST波形数据、PQRST波形的各区分点信息,可以计算出代表PQRST波形特征数据。
b24.(1)将步骤a的步骤(1)得到的心电图附加信息、步骤a的步骤(3)得到的心电图节律信息以及步骤b23得到的代表PQRST波形特征数据组成输入数据,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的组合方式可根据实际情况进行选择。
(2)将步骤(1)得到的输入数据输入传统心电图计算机自动分类算法,得到自动分类结果二。
c.对步骤b12得到的分类结果一加入权重W1来调整分类结果评分,对步骤b24得到的分类结果二加入权重W2来调整分类结果评分,所述的权重W1、权重W2满足条件W1+W2=1。
d.从步骤c的分类结果评分中取评分值最大的分类结果作为最终分类结果。
实施例5
一种基于深度学习算法的心电图综合分类方法,其流程图如图1所示,具体包括以下步骤:
a.(1)获取多导联心电图波形数据以及心电图附加信息,并根据该多导联心电图波形数据截取其中长度为8秒的数据作为原始心电图波形数据,其中多导联心电图波形数据以及心电图附加信息可以通过体检结果获得,也可以通过已有的数据库如欧共体心电图波形数据库(CSE)获得,或者通过其它途径获得,该多导联心电图波形数据可以为十二导联心电图波形数据、三导联心电图波形数据、六导联心电图波形数据、十八导联心电图波形数据等多导联心电图波形数据。心电图附加信息中包括性别、身高、胸围、体重、脂肪率、人种。
(2)根据需要,可对步骤(1)得到的原始心电图波形数据进行去噪处理,去噪处理包括以下步骤:
i1.采用高通滤波器去除基线漂移噪音;
i2.基于PQ段信号的标准方差和阈值法确认噪音是否过高;
i3.噪音过高时使用低通巴特沃斯滤波器去除噪音干扰。
(3)根据原始心电图波形数据计算PQRST波形的区分点,从而根据PQRST波形的区分点进行心电图节律信息的提取,获取心电图节律信息,其中心电图节律信息包括平均心室心率、平均RR间隔、最长RR间隔和最短RR间隔之差、RR间隔的标准方差、一致性P波信息、窦房律下各心拍的PR间隔和平均值、R波中预激波检测结果、窦房律下各心拍的QT间期和QTc间期和平均值、窦房律QRS平均波宽、窦房律P波宽和平均波宽、期外收缩信息、期外收缩类型、期外收缩形态、心房扑动的F波和心房颤动的f波的检测结果、非同期P波的检测结果,而代表PQRST波形通过以下步骤提取:
a11.通过一阶微分法和阈值法,对原始心电图波形数据进行检测,获取P波、QRS波、T波的特征点;
a12.对原始心电图波形数据中的所有PQRST波进行聚类分析,根据分析结果将拥有PQRST波数目最多的类型作为代表PQRST波形,若数目最多的类型为2个以上,选取R波平均振幅最大的类型作为代表PQRST波形类,最后使用叠加平均法计算各个心拍的PQRST波的平均波形作为代表PQRST波形。
(4)为了对第一深度学习算法以及第二深度学习算法进行训练,还需要获取训练用数据,该训练用数据可从其它相应体检结果中获取,也可以选择从现有的数据库例如从欧共体心电图波形数据库(CSE)中获取,其中训练用数据包括其它相应多导联心电图波形数据以及心电图附加信息,以欧共体心电图波形数据库(CSE)为例,训练用数据的获取步骤为:根据步骤a中的步骤(1)-步骤(3)对欧共体心电图波形数据库(CSE)中各个信息进行处理,获取欧共体心电图波形数据库(CSE)中各个心电图节律信息、代表PQRST波形数据、心电图附加信息并链接成一维输入数据,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的排列方式可根据实际情况进行选择。
b11.(1)第一深度学习算法选择迭代神经网络,并设定迭代神经网络的输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据中的心电图节律信息、代表PQRST波形数据、心电图附加信息排列成一维数据后从迭代神经网络的输入端输入,并把对应的真实分类结果输入到迭代神经网络的结果端来训练迭代神经网络。
b12.在完成对迭代神经网络的训练后,将步骤a的步骤(1)中的心电图附加信息、步骤a的步骤(3)中代表PQRST波形数据以及心电图节律信息排列成一维数据后输入到迭代神经网络的输入端,即可通过迭代神经网络进行波形分类后得到分类结果一,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的排列方式可根据实际情况进行选择。
b21.(1)第二深度学习算法选择深度神经网络,设定深度神经网络输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据从深度神经网络的输入端输入,经深度神经网络进行PQRST波群分类后,从深度神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
(3)将深度神经网络输出端得到的P波类型数据、QRS波类型数据、T波类型数据与真实波形类型进行比较,基于输出端和真实波形类型的检测偏差,利用后向传递算法修改深度神经网络内各节点的权值。
(4)重复步骤(2)-(3)直到深度神经网络判读结果符合真实波形类型,即深度神经网络结构参数收敛时,获得训练好的深度神经网络。
b22.将步骤a中的步骤(3)得到的代表PQRST波形数据从步骤(4)训练好的深度神经网络的输入端输入,经深度神经网络进行PQRST波群分类后,从深度神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
b23.使用步骤b22得到的P波类型数据、QRS波类型数据、T波类型数据分析步骤a的步骤(3)得到的代表PQRST波形数据、PQRST波形的各区分点信息,可以计算出代表PQRST波形特征数据。
b24.(1)将步骤a的步骤(1)得到的心电图附加信息、步骤a的步骤(3)得到的心电图节律信息以及步骤b23得到的代表PQRST波形特征数据组成输入数据,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的组合方式可根据实际情况进行选择。
(2)将步骤(1)得到的输入数据输入传统心电图计算机自动分类算法,得到自动分类结果二。
c.对步骤b12得到的分类结果一加入权重W1来调整分类结果评分,对步骤b24得到的分类结果二加入权重W2来调整分类结果评分,所述的权重W1、权重W2满足条件W1+W2=1。
d.从步骤c的分类结果评分中取评分值最大的分类结果作为最终分类结果。
实施例6
一种基于深度学习算法的心电图综合分类方法,其流程图如图1所示,具体包括以下步骤:
a.(1)获取多导联心电图波形数据以及心电图附加信息,并根据该多导联心电图波形数据截取其中长度为16秒的数据作为原始心电图波形数据,其中多导联心电图波形数据以及心电图附加信息可以通过体检结果获得,也可以通过已有的数据库如欧共体心电图波形数据库(CSE)获得,或者通过其它途径获得,该多导联心电图波形数据可以为十二导联心电图波形数据、三导联心电图波形数据、六导联心电图波形数据、十八导联心电图波形数据等多导联心电图波形数据。心电图附加信息中包括性别、身高、胸围、体重、脂肪率、人种。
(2)根据需要,可对步骤(1)得到的原始心电图波形数据进行去噪处理,去噪处理包括以下步骤:
i1.采用高通滤波器去除基线漂移噪音;
i2.基于PQ段信号的标准方差和阈值法确认噪音是否过高;
i3.噪音过高时使用低通巴特沃斯滤波器去除噪音干扰。
(3)根据原始心电图波形数据计算PQRST波形的区分点,从而根据PQRST波形的区分点进行心电图节律信息的提取,获取心电图节律信息,其中心电图节律信息包括平均心室心率、平均RR间隔、最长RR间隔和最短RR间隔之差、RR间隔的标准方差、一致性P波信息、窦房律下各心拍的PR间隔和平均值、R波中预激波检测结果、窦房律下各心拍的QT间期和QTc间期和平均值、窦房律QRS平均波宽、窦房律P波宽和平均波宽、期外收缩信息、期外收缩类型、期外收缩形态、心房扑动的F波和心房颤动的f波的检测结果、非同期P波的检测结果,而代表PQRST波形通过以下步骤提取:
a11.通过一阶微分法和阈值法,对原始心电图波形数据进行检测,获取P波、QRS波、T波的特征点;
a12.对原始心电图波形数据中的所有PQRST波进行聚类分析,根据分析结果将拥有PQRST波数目最多的类型作为代表PQRST波形,若数目最多的类型为2个以上,选取R波平均振幅最大的类型作为代表PQRST波形类,最后使用叠加平均法计算各个心拍的PQRST波的平均波形作为代表PQRST波形。
(4)为了对第一深度学习算法以及第二深度学习算法进行训练,还需要获取训练用数据,该训练用数据可从其它相应体检结果中获取,也可以选择从现有的数据库例如从欧共体心电图波形数据库(CSE)中获取,其中训练用数据包括其它相应多导联心电图波形数据以及心电图附加信息,以欧共体心电图波形数据库(CSE)为例,训练用数据的获取步骤为:根据步骤a中的步骤(1)-步骤(3)对欧共体心电图波形数据库(CSE)中各个信息进行处理,获取欧共体心电图波形数据库(CSE)中各个心电图节律信息、代表PQRST波形数据、心电图附加信息并链接成一维输入数据,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的排列方式可根据实际情况进行选择。
b11.(1)第一深度学习算法选择卷积神经网络,并设定卷积神经网络的输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据中的心电图节律信息、代表PQRST波形数据、心电图附加信息排列成一维数据后从卷积神经网络的输入端输入,并把对应的真实分类结果输入到卷积神经网络的结果端来训练卷积神经网络。
b12.在完成对卷积神经网络的训练后,将步骤a的步骤(1)中的心电图附加信息、步骤a的步骤(3)中代表PQRST波形数据以及心电图节律信息排列成一维数据后输入到卷积神经网络的输入端,即可通过卷积神经网络进行波形分类后得到分类结果一,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的排列方式可根据实际情况进行选择。
b21.(1)第二深度学习算法选择卷积神经网络,设定卷积神经网络输入层、隐含层、输出层的节点个数,并随机设定相邻层各节点之间的权重。
(2)将步骤a中的步骤(4)得到的训练用数据从卷积神经网络的输入端输入,经卷积神经网络进行PQRST波群分类后,从卷积神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
(3)将卷积神经网络输出端得到的P波类型数据、QRS波类型数据、T波类型数据与真实波形类型进行比较,基于输出端和真实波形类型的检测偏差,利用后向传递算法修改卷积神经网络内各节点的权值。
(4)重复步骤(2)-(3)直到卷积神经网络判读结果符合真实波形类型,即卷积神经网络结构参数收敛时,获得训练好的卷积神经网络。
b22.将步骤a中的步骤(3)得到的代表PQRST波形数据从步骤(4)训练好的卷积神经网络的输入端输入,经卷积神经网络进行PQRST波群分类后,从卷积神经网络的输出端得到P波类型数据、QRS波类型数据、T波类型数据。
b23.使用步骤b22得到的P波类型数据、QRS波类型数据、T波类型数据分析步骤a的步骤(3)得到的代表PQRST波形数据、PQRST波形的各区分点信息,可以计算出代表PQRST波形特征数据。
b24.(1)将步骤a的步骤(1)得到的心电图附加信息、步骤a的步骤(3)得到的心电图节律信息以及步骤b23得到的代表PQRST波形特征数据组成输入数据,其中心电图附加信息、代表PQRST波形数据、心电图节律信息的组合方式可根据实际情况进行选择。
(2)将步骤(1)得到的输入数据输入传统心电图计算机自动分类算法,得到自动分类结果二。
c.对步骤b12得到的分类结果一加入权重W1来调整分类结果评分,对步骤b24得到的分类结果二加入权重W2来调整分类结果评分,所述的权重W1、权重W2满足条件W1+W2=1。
d.从步骤c的分类结果评分中取评分值最大的分类结果作为最终分类结果。
实施例7
所述的步骤c中,权重W1、权重W2通过以下步骤得到:
c11.准备一组由n人次已知真实分类结果的心电图波形数据构成的数据库,设置初始的权重W1以及初始的权重W2的值,其中初始的权重W1以及初始的权重W2的值可以都设置为0.5,或者各自设置成其它值。
c12.对步骤c11中数据库里的某个心电图波形数据通过步骤b12中训练完的第一深度学习算法进行波形分类,得到分类结果评价值S1,对该心电图波形数据通过步骤b22中训练完的第二深度学习算法进行PQRST波群分类后并经步骤b23、步骤b24,得到分类结果评价值S2;
c13.在满足W1+W2=1的条件下通过公式S1×W1+S2×W2=S3调节权重W1和权重W2的值,直到综合结果评价值S3符合已知真实分类结果取值范围,得到权重W1和权重W2的值;
c14.对步骤c11中数据库里的所有心电图波形数据重复步骤c12和步骤c13;
c15.重复步骤c14直至权重W1和权重W2的值趋于稳定,得到最终的权重W1和权重W2的值。
其余实施步骤同实施例1至6任一实施例。
实施例8
本实施例将实施例1至实施例6的分类结果与传统测量法的分类结果从敏感度、特异度这两方面进行了比较,比较结果如下表所示:
从上表可知,通过本发明得到的心电图分类结果,其敏感度、特异度相比传统测量法得到的心电图分类结果提高了约10%,基本保持在98%左右,能够很好地满足实际的给医生提供辅助诊断所需要的心电图分类信息的需要。
总之,以上所述仅为本发明的较佳实施例,凡依本发明申请专利的范围所作的均等变化与修饰,皆应属本发明的涵盖范围。

Claims (9)

1.一种基于深度学习算法的心电图综合分类方法,其特征在于包括以下步骤:
a.取得测量时间在8秒以上的原始心电图波形数据、心电图附加信息,并根据原始心电图波形数据进行心电图节律信息的提取、代表PQRST波形的提取,获取心电图节律信息、代表PQRST波形数据;
b11.对第一深度学习算法的神经网络进行训练;
b12.将步骤a得到的心电图节律信息、代表PQRST波形数据、心电图附加信息排列成一维数据后通过训练完的第一深度学习算法进行波形分类,得到分类结果一;
b21.对第二深度学习算法的神经网络进行训练;
b22.将步骤a得到的代表PQRST波形数据从训练完的第二深度学习算法的输入端输入,经第二深度学习算法进行PQRST波群分类后,从第二深度学习算法的输出端得到P波类型数据、QRS波类型数据、T波类型数据;
b23.根据步骤b22得到的P波类型数据、QRS波类型数据、T波类型数据分析步骤a得到的代表PQRST波形数据,计算出代表PQRST波形特征数据;
b24.将步骤a得到的心电图附加信息、心电图节律信息以及步骤b23得到的代表PQRST波形特征数据输入传统心电图计算机自动分类算法,得到分类结果二;
c.对步骤b12得到的分类结果一加入权重W1来调整分类结果评分,对步骤b24得到的分类结果二加入权重W2来调整分类结果评分,所述的权重W1、权重W2满足条件W1+W2=1;
d.从步骤c的分类结果评分中取评分值最大的分类结果作为最终分类结果。
2.根据权利要求1所述的一种基于深度学习算法的心电图综合分类方法,其特征在于:所述的步骤c中的权重W1、权重W2通过以下步骤得到:
c11.准备一组由n人次已知真实分类结果的心电图波形数据构成的数据库;
c12.对步骤c11中数据库里的某个心电图波形数据通过步骤b12中训练完的第一深度学习算法进行波形分类,得到分类结果评价值S1,对该心电图波形数据通过步骤b22中训练完的第二深度学习算法进行PQRST波群分类后并经步骤b23、步骤b24,得到分类结果评价值S2;
c13.在满足W1+W2=1的条件下通过公式S1×W1+S2×W2=S3调节权重W1和权重W2的值,直到综合结果评价值S3符合已知真实分类结果取值范围,得到权重W1和权重W2的值;
c14.对步骤c11中数据库里的所有心电图波形数据重复步骤c12和步骤c13;
c15.重复步骤c14直至权重W1和权重W2的值趋于稳定,得到最终的权重W1和权重W2的值。
3.根据权利要求1所述的一种基于深度学习算法的心电图综合分类方法,其特征在于:所述的步骤a中,所述的心电图节律信息包括平均心室心率、平均RR间隔、最长RR间隔和最短RR间隔之差、RR间隔的标准方差、一致性P波信息、窦房律下各心拍的PR间隔和平均值、R波中预激波检测结果、窦房律下各心拍的QT间期和QTc间期和平均值、窦房律QRS平均波宽、窦房律P波宽和平均波宽、期外收缩信息、期外收缩类型、期外收缩形态、心房扑动的F波和心房颤动的f波的检测结果、非同期P波的检测结果。
4.根据权利要求1所述的一种基于深度学习算法的心电图综合分类方法,其特征在于:所述的步骤a中,代表PQRST波形的提取包括以下步骤:
a11.通过一阶微分法和阈值法,对原始心电图波形数据进行检测,获取P波、QRS波、T波的特征点;
a12.对原始心电图波形数据中的所有PQRST波进行聚类分析,根据分析结果将拥有PQRST波数目最多的类型作为代表PQRST波形,若数目最多的类型为2个以上,选取R波平均振幅最大的类型作为代表PQRST波形,最后使用叠加平均法计算各个心拍的PQRST波的平均波形作为代表PQRST波形。
5.根据权利要求1所述的一种基于深度学习算法的心电图综合分类方法,其特征在于:所述的步骤a中,所述的心电图附加信息包括性别、身高、胸围、体重、脂肪率、人种。
6.根据权利要求1所述的一种基于深度学习算法的心电图综合分类方法,其特征在于:所述的步骤b23中,所述的代表PQRST波形特征数据包括代表PQRST波形数据中每个波形的切分点、波峰电压值以及心电图节律信息中所有心拍的节律信息。
7.根据权利要求1所述的一种基于深度学习算法的心电图综合分类方法,其特征在于:所述的步骤a中的原始心电波数据为单导联数据。
8.根据权利要求1所述的一种基于深度学习算法的心电图综合分类方法,其特征在于:所述的步骤a中的原始心电图波形数据为多导联数据,所述的心电图节律信息由各个导联的心电图节律信息串联成一维数据形成,所述的代表PQRST波形数据由各个导联的代表PQRST波形数据串联成一维数据形成。
9.根据权利要求1至8中任一权利要求所述的一种基于深度学习算法的心电图综合分类方法,其特征在于:所述的第一深度学习算法、第二深度学习算法为卷积神经网络、迭代神经网络、深度神经网络中的任一种。
CN201610572208.XA 2016-07-20 2016-07-20 一种基于深度学习算法的心电图综合分类方法 Active CN106214123B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610572208.XA CN106214123B (zh) 2016-07-20 2016-07-20 一种基于深度学习算法的心电图综合分类方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610572208.XA CN106214123B (zh) 2016-07-20 2016-07-20 一种基于深度学习算法的心电图综合分类方法

Publications (2)

Publication Number Publication Date
CN106214123A true CN106214123A (zh) 2016-12-14
CN106214123B CN106214123B (zh) 2019-01-25

Family

ID=57530924

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610572208.XA Active CN106214123B (zh) 2016-07-20 2016-07-20 一种基于深度学习算法的心电图综合分类方法

Country Status (1)

Country Link
CN (1) CN106214123B (zh)

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106778685A (zh) * 2017-01-12 2017-05-31 司马大大(北京)智能系统有限公司 心电图图像识别方法、装置及服务终端
CN107451417A (zh) * 2017-09-08 2017-12-08 北京蓬阳丰业医疗设备有限公司 动态心电分析智能诊断系统及方法
CN107622259A (zh) * 2017-10-20 2018-01-23 山东大学 一种t波检测方法、心电数据分析方法及装置
CN107714023A (zh) * 2017-11-27 2018-02-23 乐普(北京)医疗器械股份有限公司 基于人工智能自学习的静态心电图分析方法和装置
CN107811649A (zh) * 2017-12-13 2018-03-20 四川大学 一种基于深度卷积神经网络的心音多分类方法
CN107832737A (zh) * 2017-11-27 2018-03-23 乐普(北京)医疗器械股份有限公司 基于人工智能的心电图干扰识别方法
CN107837082A (zh) * 2017-11-27 2018-03-27 乐普(北京)医疗器械股份有限公司 基于人工智能自学习的心电图自动分析方法和装置
CN107890348A (zh) * 2017-11-21 2018-04-10 郑州大学 一种基于深度学习法心电节拍特征自动化提取及分类方法
CN108399369A (zh) * 2018-02-02 2018-08-14 东南大学 基于分布式计算和深度学习的心电节拍分类方法
CN108509823A (zh) * 2017-02-24 2018-09-07 深圳市理邦精密仪器股份有限公司 Qrs波群的检测方法及装置
CN109009073A (zh) * 2018-07-20 2018-12-18 广州视源电子科技股份有限公司 房颤检测装置及存储介质
CN109077721A (zh) * 2018-07-20 2018-12-25 广州视源电子科技股份有限公司 房颤检测装置及存储介质
CN109259756A (zh) * 2018-09-04 2019-01-25 周军 基于非平衡训练的二级神经网络的ecg信号处理方法
CN109276242A (zh) * 2018-08-02 2019-01-29 深圳市理邦精密仪器股份有限公司 心电信号类型识别的方法和装置
CN109350032A (zh) * 2018-10-16 2019-02-19 武汉中旗生物医疗电子有限公司 一种分类方法、系统、电子设备及存储介质
CN109770861A (zh) * 2019-03-29 2019-05-21 广州视源电子科技股份有限公司 心电节律模型的训练及检测方法、装置、设备和存储介质
CN109875548A (zh) * 2019-03-24 2019-06-14 浙江好络维医疗技术有限公司 一种基于多导联综合分析的心电图特征波形聚类方法
CN109948396A (zh) * 2017-12-20 2019-06-28 深圳市理邦精密仪器股份有限公司 一种心拍分类方法、心拍分类装置及电子设备
CN110037682A (zh) * 2019-04-01 2019-07-23 上海数创医疗科技有限公司 基于改进卷积神经网络的识别心律类型的方法
CN110226917A (zh) * 2019-06-26 2019-09-13 广州视源电子科技股份有限公司 心电信号类型检测方法、装置、计算机设备和存储介质
CN110226921A (zh) * 2019-06-27 2019-09-13 广州视源电子科技股份有限公司 心电信号检测分类方法、装置、电子设备和存储介质
CN110545713A (zh) * 2017-04-19 2019-12-06 日本光电工业株式会社 生理信息波形处理方法、程序、计算机可读存储介质和生理信息波形处理装置
CN110584649A (zh) * 2019-09-09 2019-12-20 武汉大学 一种动态单导心电图中异常心电信号识别方法及系统
WO2020061987A1 (zh) * 2018-09-28 2020-04-02 深圳大学 多导联联合心电图分类方法
WO2020078200A1 (zh) * 2018-10-19 2020-04-23 中兴通讯股份有限公司 一种数据处理方法及装置、计算机可读存储介质
CN111657921A (zh) * 2020-07-01 2020-09-15 杭州脉流科技有限公司 心电异常实时监测方法、装置、计算机设备和存储介质
CN112120691A (zh) * 2020-09-22 2020-12-25 浙江智柔科技有限公司 基于深度学习的信号识别方法、装置和计算机设备
CN112244861A (zh) * 2020-10-09 2021-01-22 广东工业大学 一种单导联心电信号f波提取方法
CN113229826A (zh) * 2020-11-25 2021-08-10 联想(北京)有限公司 一种qrs波检测方法、装置及电子设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101099669A (zh) * 2007-07-13 2008-01-09 天津大学 基于最优时频空间结构码的心电数据压缩方法和解码方法
CN101969842A (zh) * 2008-01-14 2011-02-09 皇家飞利浦电子股份有限公司 房颤监测
US20110040200A1 (en) * 2009-08-12 2011-02-17 Antony Louis Piriyakumar Douglas System and method for analyzing an electrocardiogram signal
CN103038772A (zh) * 2010-03-15 2013-04-10 新加坡保健服务集团有限公司 预测患者的存活性的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101099669A (zh) * 2007-07-13 2008-01-09 天津大学 基于最优时频空间结构码的心电数据压缩方法和解码方法
CN101969842A (zh) * 2008-01-14 2011-02-09 皇家飞利浦电子股份有限公司 房颤监测
US20110040200A1 (en) * 2009-08-12 2011-02-17 Antony Louis Piriyakumar Douglas System and method for analyzing an electrocardiogram signal
CN103038772A (zh) * 2010-03-15 2013-04-10 新加坡保健服务集团有限公司 预测患者的存活性的方法

Cited By (46)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106778685A (zh) * 2017-01-12 2017-05-31 司马大大(北京)智能系统有限公司 心电图图像识别方法、装置及服务终端
CN108509823A (zh) * 2017-02-24 2018-09-07 深圳市理邦精密仪器股份有限公司 Qrs波群的检测方法及装置
CN110545713A (zh) * 2017-04-19 2019-12-06 日本光电工业株式会社 生理信息波形处理方法、程序、计算机可读存储介质和生理信息波形处理装置
CN107451417A (zh) * 2017-09-08 2017-12-08 北京蓬阳丰业医疗设备有限公司 动态心电分析智能诊断系统及方法
CN107622259A (zh) * 2017-10-20 2018-01-23 山东大学 一种t波检测方法、心电数据分析方法及装置
CN107622259B (zh) * 2017-10-20 2019-12-24 山东大学 一种t波检测方法、心电数据分析方法及装置
CN107890348A (zh) * 2017-11-21 2018-04-10 郑州大学 一种基于深度学习法心电节拍特征自动化提取及分类方法
US11324455B2 (en) 2017-11-27 2022-05-10 Shanghai Lepu CloudMed Co., LTD Artificial intelligence-based interference recognition method for electrocardiogram
US11529103B2 (en) 2017-11-27 2022-12-20 Shanghai Lepu CloudMed Co., LTD Artificial intelligence self-learning-based automatic electrocardiography analysis method and apparatus
CN107837082A (zh) * 2017-11-27 2018-03-27 乐普(北京)医疗器械股份有限公司 基于人工智能自学习的心电图自动分析方法和装置
US11344243B2 (en) 2017-11-27 2022-05-31 Shanghai Lepu CloudMed Co., LTD Artificial intelligence self-learning-based static electrocardiography analysis method and apparatus
CN107714023B (zh) * 2017-11-27 2020-09-01 上海优加利健康管理有限公司 基于人工智能自学习的静态心电图分析方法和装置
CN107832737B (zh) * 2017-11-27 2021-02-05 上海优加利健康管理有限公司 基于人工智能的心电图干扰识别方法
CN107714023A (zh) * 2017-11-27 2018-02-23 乐普(北京)医疗器械股份有限公司 基于人工智能自学习的静态心电图分析方法和装置
CN107832737A (zh) * 2017-11-27 2018-03-23 乐普(北京)医疗器械股份有限公司 基于人工智能的心电图干扰识别方法
WO2019100560A1 (zh) * 2017-11-27 2019-05-31 乐普(北京)医疗器械股份有限公司 基于人工智能自学习的心电图自动分析方法和装置
WO2019100566A1 (zh) * 2017-11-27 2019-05-31 乐普(北京)医疗器械股份有限公司 基于人工智能自学习的静态心电图分析方法和装置
CN107811649A (zh) * 2017-12-13 2018-03-20 四川大学 一种基于深度卷积神经网络的心音多分类方法
CN107811649B (zh) * 2017-12-13 2020-12-22 四川大学 一种基于深度卷积神经网络的心音多分类方法
CN109948396A (zh) * 2017-12-20 2019-06-28 深圳市理邦精密仪器股份有限公司 一种心拍分类方法、心拍分类装置及电子设备
CN108399369A (zh) * 2018-02-02 2018-08-14 东南大学 基于分布式计算和深度学习的心电节拍分类方法
CN108399369B (zh) * 2018-02-02 2021-10-19 东南大学 基于分布式计算和深度学习的心电节拍分类方法
CN109077721B (zh) * 2018-07-20 2021-03-23 广州视源电子科技股份有限公司 房颤检测装置及存储介质
CN109077721A (zh) * 2018-07-20 2018-12-25 广州视源电子科技股份有限公司 房颤检测装置及存储介质
CN109009073A (zh) * 2018-07-20 2018-12-18 广州视源电子科技股份有限公司 房颤检测装置及存储介质
CN109276242A (zh) * 2018-08-02 2019-01-29 深圳市理邦精密仪器股份有限公司 心电信号类型识别的方法和装置
CN109259756A (zh) * 2018-09-04 2019-01-25 周军 基于非平衡训练的二级神经网络的ecg信号处理方法
CN109259756B (zh) * 2018-09-04 2020-12-18 周军 基于非平衡训练的二级神经网络的ecg信号处理方法
WO2020061987A1 (zh) * 2018-09-28 2020-04-02 深圳大学 多导联联合心电图分类方法
CN109350032B (zh) * 2018-10-16 2022-03-18 武汉中旗生物医疗电子有限公司 一种分类方法、系统、电子设备及存储介质
CN109350032A (zh) * 2018-10-16 2019-02-19 武汉中旗生物医疗电子有限公司 一种分类方法、系统、电子设备及存储介质
WO2020078200A1 (zh) * 2018-10-19 2020-04-23 中兴通讯股份有限公司 一种数据处理方法及装置、计算机可读存储介质
CN109875548B (zh) * 2019-03-24 2022-04-19 浙江好络维医疗技术有限公司 一种基于多导联综合分析的心电图特征波形聚类方法
CN109875548A (zh) * 2019-03-24 2019-06-14 浙江好络维医疗技术有限公司 一种基于多导联综合分析的心电图特征波形聚类方法
CN109770861B (zh) * 2019-03-29 2022-04-12 广州视源电子科技股份有限公司 心电节律模型的训练及检测方法、装置、设备和存储介质
CN109770861A (zh) * 2019-03-29 2019-05-21 广州视源电子科技股份有限公司 心电节律模型的训练及检测方法、装置、设备和存储介质
CN110037682A (zh) * 2019-04-01 2019-07-23 上海数创医疗科技有限公司 基于改进卷积神经网络的识别心律类型的方法
CN110226917A (zh) * 2019-06-26 2019-09-13 广州视源电子科技股份有限公司 心电信号类型检测方法、装置、计算机设备和存储介质
CN110226921A (zh) * 2019-06-27 2019-09-13 广州视源电子科技股份有限公司 心电信号检测分类方法、装置、电子设备和存储介质
CN110226921B (zh) * 2019-06-27 2022-07-29 广州视源电子科技股份有限公司 心电信号检测分类方法、装置、电子设备和存储介质
CN110584649B (zh) * 2019-09-09 2021-01-15 武汉大学 一种动态单导心电图中异常心电信号识别方法及系统
CN110584649A (zh) * 2019-09-09 2019-12-20 武汉大学 一种动态单导心电图中异常心电信号识别方法及系统
CN111657921A (zh) * 2020-07-01 2020-09-15 杭州脉流科技有限公司 心电异常实时监测方法、装置、计算机设备和存储介质
CN112120691A (zh) * 2020-09-22 2020-12-25 浙江智柔科技有限公司 基于深度学习的信号识别方法、装置和计算机设备
CN112244861A (zh) * 2020-10-09 2021-01-22 广东工业大学 一种单导联心电信号f波提取方法
CN113229826A (zh) * 2020-11-25 2021-08-10 联想(北京)有限公司 一种qrs波检测方法、装置及电子设备

Also Published As

Publication number Publication date
CN106214123B (zh) 2019-01-25

Similar Documents

Publication Publication Date Title
CN106214123A (zh) 一种基于深度学习算法的心电图综合分类方法
CN106108889A (zh) 基于深度学习算法的心电图分类方法
CN106214145A (zh) 一种基于深度学习算法的心电图分类方法
CN112353407B (zh) 一种基于神经功能康复主动训练的评估系统及方法
Belgacem et al. ECG based human authentication using wavelets and random forests
CN101766484B (zh) 用于心电图识别与分类的方法和设备
CN107714023A (zh) 基于人工智能自学习的静态心电图分析方法和装置
CN106725428A (zh) 一种心电信号分类方法及装置
CN107951485A (zh) 基于人工智能自学习的动态心电图分析方法和装置
CN106725420A (zh) 室性早搏识别方法及室性早搏识别系统
CN109276242A (zh) 心电信号类型识别的方法和装置
CN110946556A (zh) 基于可穿戴式体感网的帕金森静息态震颤评估方法
CN111887858B (zh) 基于跨模态映射的心冲击图信号心率估计方法
CN112826513B (zh) 一种基于深度学习和特异性矫正在fecg上胎儿心率检测系统
CN109288515B (zh) 基于穿戴式心电信号中早搏信号的周期性监测方法及装置
CN107595276A (zh) 一种基于单导联心电信号时频特征的房颤检测方法
CN102670190A (zh) 基于心率变异性非线性特性的充盈性心衰自动诊断方法
Malek et al. Automated detection of premature ventricular contraction in ECG signals using enhanced template matching algorithm
CN109864739A (zh) 一种从10秒心电图中识别二度房室传导阻滞的方法
CN113057648A (zh) 一种基于复合lstm结构的ecg信号分类方法
CN106510687A (zh) 一种识别心电图数据异常的方法及系统
CN109893118A (zh) 一种基于深度学习的心电信号分类诊断方法
CN110175510A (zh) 基于脑功能网络特征的多模式运动想象识别方法
Costa et al. The association between reconstructed phase space and Artificial Neural Networks for vectorcardiographic recognition of myocardial infarction
Qu et al. ECG signal classification based on BPNN

Legal Events

Date Code Title Description
C06 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