CN114081502B - 一种基于机器学习的非侵入性心脏病诊断方法和装置 - Google Patents
一种基于机器学习的非侵入性心脏病诊断方法和装置 Download PDFInfo
- Publication number
- CN114081502B CN114081502B CN202111331695.8A CN202111331695A CN114081502B CN 114081502 B CN114081502 B CN 114081502B CN 202111331695 A CN202111331695 A CN 202111331695A CN 114081502 B CN114081502 B CN 114081502B
- Authority
- CN
- China
- Prior art keywords
- feature vector
- sample space
- qrs complex
- signal
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000003745 diagnosis Methods 0.000 title claims abstract description 26
- 238000010801 machine learning Methods 0.000 title claims abstract description 15
- 208000019622 heart disease Diseases 0.000 title abstract description 19
- 239000013598 vector Substances 0.000 claims abstract description 36
- 238000000718 qrs complex Methods 0.000 claims abstract description 32
- 239000011159 matrix material Substances 0.000 claims abstract description 7
- 238000013186 photoplethysmography Methods 0.000 claims description 24
- 238000004422 calculation algorithm Methods 0.000 claims description 16
- 238000000605 extraction Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 14
- 230000015654 memory Effects 0.000 claims description 13
- 230000006870 function Effects 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 11
- 230000004927 fusion Effects 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000012706 support-vector machine Methods 0.000 claims description 6
- 238000000513 principal component analysis Methods 0.000 claims description 5
- 238000012216 screening Methods 0.000 claims description 5
- 238000005520 cutting process Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000003064 k means clustering Methods 0.000 claims description 2
- 230000002708 enhancing effect Effects 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 18
- 230000008569 process Effects 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 6
- 210000004027 cell Anatomy 0.000 description 6
- 238000012549 training Methods 0.000 description 6
- 238000013528 artificial neural network Methods 0.000 description 5
- 230000000747 cardiac effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000036541 health Effects 0.000 description 4
- 230000002861 ventricular Effects 0.000 description 4
- 208000020446 Cardiac disease Diseases 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000028161 membrane depolarization Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 206010003658 Atrial Fibrillation Diseases 0.000 description 2
- 208000024172 Cardiovascular disease Diseases 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008602 contraction Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 241000544061 Cuculus canorus Species 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000001746 atrial effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000004195 computer-aided diagnosis Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 210000005240 left ventricle Anatomy 0.000 description 1
- 230000007787 long-term memory Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 230000034225 regulation of ventricular cardiomyocyte membrane depolarization Effects 0.000 description 1
- 230000002336 repolarization Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 210000005241 right ventricle Anatomy 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 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/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
-
- 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/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/366—Detecting abnormal QRS complex, e.g. widening
-
- 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
-
- 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
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/726—Details of waveform analysis characterised by using transforms using Wavelet transforms
-
- 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
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- 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
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification 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)
- Cardiology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于机器学习的非侵入性心脏病诊断方法和装置。该方法包括:采集目标的心电图信号并检测对应的QRS波群;采集该目标的脉搏波信号;基于所述QRS波群提取心电信号特征,构建为第一特征向量的样本空间;提取所述脉搏波信号的特征,将得到的特征向量矩阵构建为第二特征向量的样本空间;将所述第一特征向量的样本空间和所述第二特征向量的样本空间进行融合,利用融合特征输入到预训练的分类器,获得心脏诊断结果。本发明提高了心脏识别的准确率和效率。
Description
技术领域
本发明涉及医学图像分析技术领域,更具体地,涉及一种基于机器学习的非侵入性心脏病诊断方法和装置。
背景技术
随着心脏病患者的不断增加,有限的医疗资源无法满足日益增长的心血管疾病患者的就医需求。目前,心脏病已经成为危害人类健康的重大疾病。近年来,人工智能在医疗行业中逐渐受到人们关注,大量机器学习算法在医疗保健领域广泛应用,涉及智能风险评估、辅助诊断、医学影像、药物研发等多个领域,由此可见“互联网+医疗健康”模式将在未来成为人们日常生活的一部分。心电图(ECG)是医生诊断心血管疾病的重要依据,能真实地反映心脏的健康状况。同时,随着传感器技术的进步,基于脉搏波信号(PPG)的可穿戴设备越来越多,人们可以通过这些可穿戴设备收集人体生理体征信号,随时监控自身健康状况。
在传统的心脏诊断中,医生通过查看心电图(Electrocardiography,ECG)判断心脏病,但这种诊断方式带有主观性,且耗时耗力。Kora等基于ECG使用Cuckoo Search算法实现自动识别心脏病;王凯基于ECG的自动心律失常检测模型能够有效监测心脏异常信号的ECG信号特征。心脏病监测的困难主要在于阵发性房颤,此类心脏病的特点是突发突止、持续时间短,24h动态心电图常作为首选的监测方法。李芳等通过分析总结心脏病的心电散点图特征,分析24h动态心电图来诊断心脏病。武求花等运用Lorenz-RR散点图(Lorenz散点图)与时间-RR间期散点图(时间散点图)及逆向技术相结合的方法,快速识别不同发作时间的PAF,为临床提供更快更准确的心电诊断技术。目前国内外对于心脏识别的研究大多基于ECG,但ECG设备普遍有操作复杂、携带困难、价格较高、多通道等局限性。光电容积脉搏波描记(PPG)是使用激光对待测流体或固体的体积进行检测与记录的方法,可以用来监测毛细血管中血液的体积变化。有研究表明。PPG方法采集信息容易、方便携带,相较于ECG方法,PPG更加擅长长期监测。但是PPG可能没有ECG准确。
综上,在现有技术中,通常仅采用ECG或PPG进行心脏诊断分析,其准确性和效率还有待改进。
发明内容
本发明的目的是克服上述现有技术的缺陷,提供一种基于机器学习的非侵入性心脏病诊断方法和装置。
根据本发明的第一方面,提供一种基于机器学习的非侵入性心脏病诊断方法。该方法包括以下步骤:
采集目标的心电图信号并检测对应的QRS波群;
采集该目标的脉搏波信号;
基于所述QRS波群提取心电图信号特征,构建为第一特征向量的样本空间;
提取所述脉搏波信号的特征,将得到的特征向量矩阵构建为第二特征向量的样本空间;
将所述第一特征向量的样本空间和所述第二特征向量的样本空间进行融合,利用融合特征输入到预训练的分类器,获得心脏诊断结果。
根据本发明的第二方面,提供一种基于机器学习的非侵入性心脏病诊断装置。该装置包括:
第一数据采集单元:用于采集目标的心电图信号并检测对应的QRS波群;
第二数据采集单元:用于采集该目标的脉搏波信号;
第一特征提取单元:用于基于所述QRS波群提取心电图信号特征,构建为第一特征向量的样本空间;
第二特征提取单元:用于提取所述脉搏波信号的特征,将得到的特征向量矩阵构建为第二特征向量的样本空间;
分类识别单元:用于将所述第一特征向量的样本空间和所述第二特征向量的样本空间进行融合,利用融合特征输入到预训练的分类器,获得心脏诊断结果。
与现有技术相比,本发明的优点在于,提出了基于ECG和PPG相结合的新型心脏识别方法,通过将PPG信号的特征和ECG信号的特征融合,放入分类器进行识别,提高了利用机器学习的方法对心脏识别的准确率。
通过以下参照附图对本发明的示例性实施例的详细描述,本发明的其它特征及其优点将会变得清楚。
附图说明
被结合在说明书中并构成说明书的一部分的附图示出了本发明的实施例,并且连同其说明一起用于解释本发明的原理。
图1是根据本发明一个实施例的基于机器学习的非侵入性心脏病诊断方法的流程图;
图2是根据本发明一个实施例的心电周期示意图;
图3是根据本发明一个实施例的循环神经网络与长短时记忆网络的结构对比图;
图4是根据本发明一个实施例的长短时记忆网络的结构示意图;
图5是根据本发明一个实施例的基于多集判别相关分析的ECG和PPG特征层融合识别流程图。
具体实施方式
现在将参照附图来详细描述本发明的各种示例性实施例。应注意到:除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本发明的范围。
以下对至少一个示例性实施例的描述实际上仅仅是说明性的,决不作为对本发明及其应用或使用的任何限制。
对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论,但在适当情况下,所述技术、方法和设备应当被视为说明书的一部分。
在这里示出和讨论的所有例子中,任何具体值应被解释为仅仅是示例性的,而不是作为限制。因此,示例性实施例的其它例子可以具有不同的值。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步讨论。
本发明利用计算机辅助诊断心脏病,通过机器学习相关的知识分别对于ECG和PPG的分类进行了相关的研究。例如,首先,利用LSTM提取PPG信号的特征;然后,根据特征提取的结果进行了样本的截裁;利用Pan-Tompkins算法,使用收集的心电图数据集对于心电信号QRS波群进行特征提取,进行样本的聚类分析,以剔除了部分特殊的截裁小样;在数据筛选完成后,进行样本特征的选择,利用主成分分析法(PCA),进行了降低维度的处理;最后,利用样本特征训练分类器,进而利用经训练的分类器实现对目标患者的心脏诊断结果。
具体地,结合图1和图5所示,所提供的基于机器学习的非侵入性心脏病诊断方法包括以下步骤。
步骤S110,提取心电图信号中QRS波群信息。
心电图信号记录一定时间内心脏不同部位去极化产生的电压。QRS波群反映了左右心室的快速去极化的过程。心电信号具有周期性,相对于其它电信号例如脑电波,肌肉电波,属于比较容易判断分辨特征的一类电波。分析心电信号,必要环节是提取其QRS波群,基于QRS波群进行心电信号分析,有利于提高鲁棒性和稳定性。
QRS波群是心室去极化的电位变化,如图2所示,在心电图信号中,在P波结束后,进入PR区间,之后Q波会呈现出一个向下的趋势,R波则会快速向上,S波紧接着向下,这三个波所对应的连接电波,本文统称为QRS波群。
由图2可以明显看出,R波是所有波形中最窄,幅度最大的波。T波对应的是人体心室肌肉的复极过程。心脏去极化的过程对应P波,因此它是钝圆形的图像。PR对应的是传导过程,在理想状态下实际上是没有电位变化的,即心电信号为平。QT期间对应人体心脏心室收缩的一个过程,因此通过QT区间的变化可以有效的判断人的心率,QT期间的大小与人心率的快慢呈正比关系。
在正常的心电信号的周期中,P波和T波没有明显的特征,为了明显分辨心电信号,基本上还是需要通过QRS波群的位置来对其进行判断。选择QRS波群还有一个非常重要的原因是因为QRS波群是心电信号中最明显的信号,正确选择了QRS波群才能进一步完成对其它信号波的选择。
对于QRS波群的提取和检测方法包括小波变换、神经网络和Pan-Tompkins算法及其改良算法、增加了陷波器的小波改良算法等。
小波变换方法准确性高,但是缺点在于计算步骤非常繁琐,因此计算效率较低,并且在小波变换过程中会导致数据的冗余度增加。
神经网络一定程度上可以准确检测QRS波群,但训练神经网络的过程非常漫长的,并且需要采集大量的样本,此外不同的训练样本产生的训练结果也并不相同,而由于心脏病患者的心电图信号差异较大,经训练的神经网络模型并不适用于实际的应用。
为了平衡心电信号分析的准确性和计算效率,在一个实施例中,采用Pan-Tompkins算法检测心电图中QRS波群信息。相对于其他的QRS波群提取算法,Pan-Tompkins算法计算效率高并且不需要模型训练。简言之,利用Pan-Tompkins法检测R波峰值包括:利用滤波器进行滤波;对滤波后的信号求一阶导数;对求导之后的信号进行平方运算;将信号通过滑动窗口进行积分;使用阈值法检测经过处理之后的R波峰值等过程。
具体地,基于Pan-Tompkins算法提取心电图信号的特征包括:
步骤S111,利用滤波器进行降噪处理。
因为信号采集设备和环境的不同,在心电信号采集过程中会引入不同的噪声,所以对于心电信号数据中的噪声处理是数据预处理的核心内容。
在一个实施例中,采用高通滤波器和低通滤波器来滤除信号中的噪声,从而提升心电信号的信噪比。
心电信号中主要的噪声为工频干扰和基线漂移。电力***等带入的干扰信号称为工频干扰,频率一般集中在50Hz。而皮肤的接触,人的呼吸等都会产生所谓的基线漂移,频率一般集中在0.05Hz至2Hz。而QRS波群的频率主要集中在5Hz至11Hz。
以T为采样周期,低通滤波器的截止频率为11Hz,差分方程可表示为:
y(nT)=2y(nT-T)-y(nT-2T)
-2x(nT-6T)+x(nT-12T)+x(nT) (1)
高通滤波器的截至频率为2Hz,差分方程可表示为:
y(nT)=y(nT-T)-x(nT)
+32x(nT-16T)+x(nT-32T) (2)
步骤S112,进行放大处理以检测出R波。
QRS波群中R波是位于一个心电信号周期中最陡的一个波。为了方便将QRS波从其他的心电信号波中区分出来,需要利用相应的微分方法对于QRS波群进行陡坡特性的放大,具体的差分方程如下:
y(nT)=x(nT)-x(xT-2T) (3)
高频特性可以通过非线性平方函数来增强,其方程可以表示为:
y(nT)=[x(nT)]2 (4)
步骤S113,确定R波外的其他特征
上述的方法基本都是围绕R波来展开进行的,但单纯的对于R波的处理对于心电信号的异常判断无疑是不够完整的。因此对于波群的其它特征同样也需要进行确定。
在一个实施例中,采取了一个移动窗积分来确定其它相关的特征。目的是为了提取波群的其他特征。设积分窗内样本数为N,差分方程表示为:
步骤S114,阈值动态调整和搜索
阈值动态调整和搜索属于该算法的核心部分,通过其搜索心电信号峰值,通过一定的计算式更新阈值,并结合回溯检测,双向判别等内容来检测。该过程可采用现有技术实现,在此不再赘述。
步骤S120,利用长短时记忆网络提取脉搏波特征。
与传统RNN不同,长短时记忆网络(LSTM)不存在消失梯度问题,可记忆之前数千个离散时间步长的输入,与传统RNN相比,LSTM可在更长的序列中有更好表现。LSTM将神经元的激活功能替换为一个具有“门”控制结构的单元,每个单元称为一个细胞(Cell),本发明采用LSTM使用具有遗忘门的存储单元。
如图3所示,相比于RNN只有一个传递状态h(t),LSTM则有两个传输状态c(t)(cellstate,细胞状态)和h(t)(hidden state,隐藏状态),RNN中的h(t)相当于LSTM的c(t)。LSTM对于传递下去的c(t)改变很慢,通常输出的c(t)是上一个状态传过来的c(t-1)加上一些数值,而h(t)在不同节点下往往有很大区别,正是这些波动较大的h(t)产生忘记和选择记忆的效果。
首先使用当前输入x(t)和上一个状态传递下来的h(t-1)拼接训练得到4中状态(如图4所示):z、zi、zf、zo,其中zf、zi、zo分别作为遗忘门(forget gate),输入门(inputgate)和输出门(output gate)的控制信号,z是神经网络输入层的输入。4种状态计算公式分别为:
zi(t)=σ(Wxi⊙x(t)+Whi⊙h(t-1)) (6)
zo(t)=σ(Wxo⊙x(t)+Who⊙h(t-1)) (8)
zf(t)=σ(Wxf⊙x(t)+Whf⊙h(t-1)) (9)
其中,W表示对应项的权重。
在该步骤中,对于光电容积脉搏波数据(PPG)序列通过忘记、选择记忆、输出阶段等得出符合生理规律的特征参数。PPG信号中包含大量例如与心率相关的特征,这些特征有利于LSTM网络提高预测精度。需说明的是,也可采用GRU(门控循环单元)替代LSTM。
步骤S130,融合基于QRS波群提取的心电信号特征和脉搏波特征。
在获得ECG信号特征和PPG信号特征后,可以在截裁样本后使用K-Means等聚类进行数据筛选,以进一步进行特征融合。
在一个实施例中,特征融合过程包括:
步骤S131,对ECG信号特征提取后得到的特征向量组成第一特征的样本空间X∈Rp×n。
例如,ECG信号特征反映波群形态特征和间期特征等。心电信号特征包括心电信号的时域特征、心电信号的变换域特征等。
步骤S132,对PPG信号特征提取后得到的特征向量矩阵组成第二特征的样本空间Y∈Rq×n。
步骤S133,采用最大最小值法,对特征向量进行归一化操作,使得两组特征向量的数量级接近。
优选地,通过对特征向量进行归一化处理,以消除奇异样本数据对后续预测结果的不良影响。
步骤S140,利用融合特征获得心脏病分类结果。
在该步骤中,利用融合特征训练分类器,进而利用经训练的分类器进行实时心脏病诊断分析。
在一个实施例中,选用支持向量机(SVM)作为分类器。这是因为SVM是定义在特征空间上的间隔最大的分类器,并更适用于小样本学习,并且SVM最终的决策结果取决于少数支持向量,从而可以“剔除”大量冗余样本,计算效率较高并具有较好的“鲁棒”性。由于大量心脏病患者的心电图样本和PPG信号样本可能不易获得,并且在实际应用中,患者的相关数据可能差异较大,因此选用SVM可获得更准确的诊断结果,如正常心跳、左心室动脉受阻心率、右心室动脉受阻心率、房颤、心室收缩和心房收缩等。
此外,为了进一步降维特征数据并保证后续诊断结果的准确性,在此步骤中,还可以利用主成分分析选择关键特征作为分类器的输入。
相应地,本发明还提供一种基于机器学习的非侵入性心脏病诊断装置,用于实现上述方法的一个方面或多个方面。例如,该装置包括:第一数据采集单元,其用于采集目标的心电图信号并检测对应的QRS波群;第二数据采集单元,其用于采集该目标的脉搏波信号;第一特征提取单元,其用于基于所述QRS波群提取心电图信号特征,构建为第一特征向量的样本空间;第二特征提取单元,其用于提取所述脉搏波信号的特征,将得到的特征向量矩阵构建为第二特征向量的样本空间;分类识别单元,其用于将所述第一特征向量的样本空间和所述第二特征向量的样本空间进行融合,利用融合特征输入到预训练的分类器,获得心脏诊断结果。该装置所涉及的各单元可采用处理器、FPGA或专用硬件实现。
综上,本发明使用医院患者的数据,分别收集患者的脉搏波信号和心电图,首先利用LSTM提取PPG信号的特征。其次利用Pan-Tompkins算法,使用收集的心电图数据集对于心电信号QRS波群进行特征提取,在截裁样本后使用K-Means聚类进行数据筛选。再其次将提取的PPG信号特征和ECG信号特征进行融合,最后利用分类器对融合的特征进行分类。本发明提出了一种基于ECG和PPG相结合的新型心脏识别方法,首次将ECG和PPG相结合,为心脏识别提供了补充的解决方案,并提高了利用机器学习对心脏识别的准确率,获得的心脏诊断分类结果可用于临床指示。
本发明可以是***、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质,其上载有用于使处理器实现本发明的各个方面的计算机可读程序指令。
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是但不限于电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身,诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如,通过光纤电缆的光脉冲)、或者通过电线传输的电信号。
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备,或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令,并转发该计算机可读程序指令,以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本发明操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码,所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++、Python等,以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中,通过利用计算机可读程序指令的状态信息来个性化定制电子电路,例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA),该电子电路可以执行计算机可读程序指令,从而实现本发明的各个方面。
这里参照根据本发明实施例的方法、装置(***)和计算机程序产品的流程图和/或框图描述了本发明的各个方面。应当理解,流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合,都可以由计算机可读程序指令实现。
这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器,从而生产出一种机器,使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时,产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中,这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作,从而,存储有指令的计算机可读介质则包括一个制造品,其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。
也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上,使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤,以产生计算机实现的过程,从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。
附图中的流程图和框图显示了根据本发明的多个实施例的***、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分,所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的***来实现,或者可以用专用硬件与计算机指令的组合来实现。对于本领域技术人员来说公知的是,通过硬件方式实现、通过软件方式实现以及通过软件和硬件结合的方式实现都是等价的。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。本发明的范围由所附权利要求来限定。
Claims (5)
1.一种基于机器学习的非侵入性心脏病诊断装置,包括:
第一数据采集单元:用于采集目标的心电图信号并检测对应的QRS波群;
第二数据采集单元:用于采集该目标的脉搏波信号;
第一特征提取单元:用于基于所述QRS波群提取心电图信号特征,构建为第一特征向量的样本空间;
第二特征提取单元,用于利用长短时记忆网络提取所述脉搏波信号的特征,将得到的特征向量矩阵构建为第二特征向量的样本空间;
分类识别单元:用于将所述第一特征向量的样本空间和所述第二特征向量的样本空间归一化处理后进行融合,并利用主成分分析选择关键的融合特征,利用选出的融合特征输入到预训练的分类器,获得心脏诊断结果;
其中,所述采集目标的心电图信号并检测对应的QRS波群包括:
对于采集的心电图信号,利用Pan-Tompkins算法进行R波检测;
基于检测到的R波对信号进行切割,获得RR间期特征,并在截裁样本使用聚类算法进行数据筛选;
其中,检测对应的QRS波群包括:
分别采用低通滤波器和高通滤波器对采集的心电图信号进行滤波,其中低通滤波器的截止频率设置为11Hz ,高通滤波器的截至频率设置为2Hz ;
利用微分方法对于QRS波群进行陡坡特性的放大并通过非线性平方函数增强心电图信号中的高频特性;
利用移动窗积分确定波群中R波之外的其他特征;
利用阈值动态调整和搜索确定心电图信号的峰值。
2.根据权利要求1所述的装置,其特征在于,所述分类器是支持向量机。
3.根据权利要求1所述的装置,其特征在于,所述脉搏波信号利用可穿戴设备基于光电容积脉搏波描记法获得。
4.根据权利要求1所述的装置,其特征在于,所述聚类算法是K-Means聚类。
5.一种计算机可读存储介质,其上存储有计算机程序,其中,该程序被处理器执行时实现以下步骤:
采集目标的心电图信号并检测对应的QRS波群;
采集该目标的脉搏波信号;
基于所述QRS波群提取心电图信号特征,构建为第一特征向量的样本空间;
利用长短时记忆网络提取所述脉搏波信号的特征,将得到的特征向量矩阵构建为第二特征向量的样本空间;
将所述第一特征向量的样本空间和所述第二特征向量的样本空间归一化处理后进行融合,并利用主成分分析选择关键的融合特征,利用选出的融合特征输入到预训练的分类器,获得心脏诊断结果;
其中,所述采集目标的心电图信号并检测对应的QRS波群包括:
对于采集的心电图信号,利用Pan-Tompkins算法进行R波检测;
基于检测到的R波对信号进行切割,获得RR间期特征,并在截裁样本使用聚类算法进行数据筛选;
其中,检测对应的QRS波群包括:
分别采用低通滤波器和高通滤波器对采集的心电图信号进行滤波,其中低通滤波器的截止频率设置为11Hz ,高通滤波器的截至频率设置为2Hz ;
利用微分方法对于QRS波群进行陡坡特性的放大并通过非线性平方函数增强心电图信号中的高频特性;
利用移动窗积分确定波群中R波之外的其他特征;
利用阈值动态调整和搜索确定心电图信号的峰值。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111331695.8A CN114081502B (zh) | 2021-11-11 | 2021-11-11 | 一种基于机器学习的非侵入性心脏病诊断方法和装置 |
PCT/CN2021/137920 WO2023082402A1 (zh) | 2021-11-11 | 2021-12-14 | 一种基于机器学习的非侵入性心脏病诊断方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111331695.8A CN114081502B (zh) | 2021-11-11 | 2021-11-11 | 一种基于机器学习的非侵入性心脏病诊断方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114081502A CN114081502A (zh) | 2022-02-25 |
CN114081502B true CN114081502B (zh) | 2022-09-27 |
Family
ID=80299785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111331695.8A Active CN114081502B (zh) | 2021-11-11 | 2021-11-11 | 一种基于机器学习的非侵入性心脏病诊断方法和装置 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN114081502B (zh) |
WO (1) | WO2023082402A1 (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104706349A (zh) * | 2015-04-13 | 2015-06-17 | 大连理工大学 | 一种基于脉搏波信号的心电信号构建方法 |
CN108113663A (zh) * | 2016-11-30 | 2018-06-05 | 安徽华米信息科技有限公司 | 心脏状况检测***及方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103284702A (zh) * | 2013-06-17 | 2013-09-11 | 中国科学院苏州纳米技术与纳米仿生研究所 | 心电图和脉搏波关系分析及融合分析的方法和装置 |
EP3381356A1 (en) * | 2017-03-28 | 2018-10-03 | IMEC vzw | System and method for cuffless blood pressure estimation |
EP3449820B1 (en) * | 2017-08-30 | 2023-05-24 | Qompium | Computer-implemented method and system for direct photoplethysmography (ppg) |
EP3856016A4 (en) * | 2018-09-24 | 2022-06-15 | Sotera Wireless, Inc. | METHOD AND SYSTEM FOR MONITORING A PATIENT FOR ATRIAL FIBRILLATION AND/OR ASYSTOL |
CN109497992A (zh) * | 2019-01-04 | 2019-03-22 | 济南汇医融工科技有限公司 | 基于机器学习方法的冠心病智能筛查装置 |
CN111358450B (zh) * | 2020-03-17 | 2023-04-07 | 乐普(北京)医疗器械股份有限公司 | 一种血压分类方法和装置 |
CN111839505A (zh) * | 2020-07-02 | 2020-10-30 | 复旦大学 | 一种基于心血管***电-机械活动信息的房颤分级方法 |
CN111657925A (zh) * | 2020-07-08 | 2020-09-15 | 中国科学院苏州生物医学工程技术研究所 | 基于机器学习的心电信号分类方法、***、终端以及存储介质 |
CN112869717B (zh) * | 2021-02-25 | 2023-02-24 | 佛山科学技术学院 | 一种基于bl-cnn的脉搏特征识别分类***及方法 |
CN113349753A (zh) * | 2021-07-19 | 2021-09-07 | 成都芯跳医疗科技有限责任公司 | 一种基于便携式动态心电监护仪的心律失常检测方法 |
-
2021
- 2021-11-11 CN CN202111331695.8A patent/CN114081502B/zh active Active
- 2021-12-14 WO PCT/CN2021/137920 patent/WO2023082402A1/zh unknown
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104706349A (zh) * | 2015-04-13 | 2015-06-17 | 大连理工大学 | 一种基于脉搏波信号的心电信号构建方法 |
CN108113663A (zh) * | 2016-11-30 | 2018-06-05 | 安徽华米信息科技有限公司 | 心脏状况检测***及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114081502A (zh) | 2022-02-25 |
WO2023082402A1 (zh) | 2023-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mousavi et al. | ECGNET: Learning where to attend for detection of atrial fibrillation with deep visual attention | |
Lin | Frequency-domain features for ECG beat discrimination using grey relational analysis-based classifier | |
Oweis et al. | QRS detection and heart rate variability analysis: A survey | |
Gupta et al. | Efficient R-peak detection in electrocardiogram signal based on features extracted using Hilbert transform and Burg method | |
US20170112401A1 (en) | Automatic method to delineate or categorize an electrocardiogram | |
Haleem et al. | Time adaptive ECG driven cardiovascular disease detector | |
Zhang et al. | Localization of myocardial infarction with multi-lead bidirectional gated recurrent unit neural network | |
Sumathi et al. | A wavelet transform based feature extraction and classification of cardiac disorder | |
Gupta | Wavelet transform and vector machines as emerging tools for computational medicine | |
Singh et al. | Short and noisy electrocardiogram classification based on deep learning | |
Moghadam et al. | Automatic diagnosis and localization of myocardial infarction using morphological features of ECG signal | |
Patel et al. | Machine learning for biomedical signal processing | |
Reddy et al. | Classification of arrhythmia disease through electrocardiogram signals using sampling vector random forest classifier | |
Paul et al. | Automated detection of cardinal points of ECG signal for feature extraction using a single median filter | |
KR20140097039A (ko) | 자동 결합 뉴런 네트워크를 이용한 심장부정맥 분류 방법 및 장치 | |
Simjanoska et al. | ECG-derived Blood Pressure Classification using Complexity Analysis-based Machine Learning. | |
CN114081502B (zh) | 一种基于机器学习的非侵入性心脏病诊断方法和装置 | |
Huo et al. | ECG segmentation algorithm based on bidirectional hidden semi-Markov model | |
Hasan et al. | Cardiac arrhythmia detection in an ECG beat signal using 1D convolution neural network | |
Agrawal et al. | A machine learning approach for classifying healthy and infarcted patients using heart rate variabilities derived vector magnitude | |
Maghawry et al. | An efficient approach for paroxysmal atrial fibrillation events prediction using extreme learning machine | |
Wang et al. | A high precision real-time premature ventricular contraction assessment method based on the complex feature set | |
Kammath et al. | Detection of bundle branch blocks using machine learning techniques | |
Wadhwani et al. | IOT based biomedical wireless sensor networks and machine learning algorithms for detection of diseased conditions | |
Sanamdikar et al. | Using the GAN method, analysis several characteristics of the ECG signal in order to detect cardiac arrhythmia |
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 |