CN110226921B - 心电信号检测分类方法、装置、电子设备和存储介质 - Google Patents
心电信号检测分类方法、装置、电子设备和存储介质 Download PDFInfo
- Publication number
- CN110226921B CN110226921B CN201910567569.9A CN201910567569A CN110226921B CN 110226921 B CN110226921 B CN 110226921B CN 201910567569 A CN201910567569 A CN 201910567569A CN 110226921 B CN110226921 B CN 110226921B
- Authority
- CN
- China
- Prior art keywords
- wave
- sequence
- interval
- band
- ratio
- 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
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
-
- 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
-
- 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/361—Detecting fibrillation
-
- 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)
- 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)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本申请涉及一种心电信号检测分类方法、装置、电子设备和存储介质。包括:在心电信号中,提取信号波形;获取信号波形的形态特征和深层特征;形态特征包括TR波幅度差特征、PR波数量比值特征、ST波段特征和P波变化特征中的任意一种;深层特征包括深度特征和层次特征;将形态特征和深层特征输入至分类器;获取分类器输出的分类结果,得到心电信号的信号类型;分类结果为分类器根据形态特征和深层特征进行分类的结果;信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。采用上述方案可以准确识别出各种类型的异常节律,避免了将心动过速,心动过缓,心律失常等非房颤类型的异常节律错误分类为房颤类型的情况,提升了心电信号分类的准确性。
Description
技术领域
本申请涉及医疗器械与医疗产品领域,特别是涉及一种心电信号检测分类方法、装置、电子设备和存储介质。
背景技术
心房颤动(Atrial Fibrillation,AF)简称房颤,是临床最常见的一种心律失常疾病,其特点是紊乱的心房活动和随之而来的脑卒、心肌梗死等并发症,导致较高的致残率和死亡率,严重危害着人类的健康和生命。研究心电信号是否存在房颤的算法,能及早发现并进行治疗,从而可更多的抓住治疗的最佳时机,减少房颤的发病率和死亡率,因此具有重要的临床和社会意义。
由于房颤发作时的两个重要临床表现为RR间期绝对不规则和P波消失,因此,目前较为常规的心电信号分类方法,主要是基于心房活动(P波消失)或心室活动(RR间期不规则)进行分析。
然而,心电信号中的P波或f波属于微弱信号,其波形特征较难检测。而且,许多非房颤类型的异常节律(如心动过速,心动过缓,心律失常等)表现出类似房颤的特征,如果基于常规的心电信号分类方法,可能会出现错误的识别结果。
因此,现有技术中的心电信号分类方法存在着房颤分类错误的问题。
发明内容
基于此,有必要针对上述技术问题,提供一种心电信号检测分类方法、装置、电子设备和存储介质。
第一方面,提供了一种心电信号检测分类方法,所述方法包括:
在心电信号中,提取信号波形;
获取所述信号波形的形态特征和深层特征;所述形态特征包括TR波幅度差特征、PR波数量比值特征、ST波段特征、P波变化特征中的任意一种;所述深层特征包括深度特征和层次特征;
将所述形态特征和所述深层特征输入至分类器;
获取所述分类器输出的分类结果,得到所述心电信号的信号类型;所述分类结果为所述分类器根据所述形态特征和所述深层特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
在一个实施例中,当所述形态特征为TR波幅度差特征时,所述获取所述信号波形的形态特征,包括:
获取R波信号序列,以及,获取T波信号序列;
计算所述T波信号序列与所述R波信号序列的比值,得到TR波序列比值;
计算所述TR波序列比值的中位数,得到TR波序列比值中位数;
确定所述TR波序列比值中位数和所述TR波序列比值,作为所述TR波幅度差特征。
在一个实施例中,所述将所述形态特征输入至分类器,包括:
在所述信号波形中,提取RR间期差序列;
根据所述RR间期差序列,获取RR间期差第一先验特征;所述RR间期差第一先验特征包括所述RR间期差序列的原始计数、近似熵、样本熵中的至少一种;
将所述RR间期差第一先验特征和所述TR波幅度差特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第一先验特征和所述TR波幅度差特征进行分类的结果。
在一个实施例中,当所述RR间期差第一先验特征为所述RR间期差序列的原始计数时,所述根据所述RR间期差序列,获取RR间期差第一先验特征,包括:
计算所述RR间期差序列的绝对值,得到RR间期差序列绝对值;
当所述RR间期差序列绝对值小于等于预设的常量定值时,计算所述RR间期差序列之和,得到RR间期差序列和值;
当所述RR间期差序列和值等于二时,确定所述RR间期差序列和值,作为所述RR间期差序列的原始计数。
在一个实施例中,当所述形态特征为PR波数量比值特征时,所述获取所述信号波形的形态特征,包括:
获取R波序列数量,以及,获取P波序列数量;所述P波序列数量为P波序列的总和;所述R波序列数量为R波序列的总和;
计算所述P波序列数量与所述R波序列数量的比值,得到PR波数量比值;
根据所述PR波数量比值,确定所述PR波数量比值特征;所述PR波数量比值特征包括所述PR波数量比值、PR波数量比值中位数、PR波数量比值标准差、PR波数量比值均值以及PR波数量比值方差中的至少一种。
在一个实施例中,当所述PR波数量比值特征为所述PR波数量比值方差时,所述根据所述PR波数量比值,确定所述PR波数量比值特征,包括:
计算所述PR波数量比值的均值,得到所述PR波数量比值均值,以及,计算所述PR波数量比值的标准差,得到所述PR波数量比值标准差;
计算所述PR波数量比值标准差与所述PR波数量比值均值的比值,得到所述PR波数量比值方差。
在一个实施例中,所述将所述形态特征输入至分类器,包括:
在所述信号波形中,提取RR间期序列和/或RR间期差序列;
根据所述RR间期序列和/或所述RR间期差序列,获取RR间期差第二先验特征;
将所述RR间期差第二先验特征和所述PR波数量比值特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第二先验特征和所述PR波数量比值特征进行分类的结果。
在一个实施例中,所述根据所述RR间期序列,获取RR间期差第二先验特征,包括:
计算所述RR间期序列的核密度,得到RR间期序列核密度;
计算所述RR间期序列核密度的峰度值,得到RR间期密度峰度,以及,计算所述RR间期序列核密度的偏度值,得到RR间期密度偏度;
确定所述RR间期密度偏度、所述RR间期密度峰度以及所述RR间期序列核密度中的至少一种,作为所述RR间期差第二先验特征。
在一个实施例中,所述根据所述RR间期序列和/或所述RR间期差序列,获取RR间期差第二先验特征,包括:
计算所述RR间期差序列的核密度,得到RR间期差序列核密度;
计算所述RR间期差序列和密度的峰度值,得到RR间期差密度峰度,以及,计算所述RR间期差序列核密度的偏度值,得到RR间期差密度偏度;
确定所述RR间期差密度偏度、所述RR间期差密度峰度以及所述RR间期差序列核密度中的至少一种,作为所述RR间期差第二先验特征。
在一个实施例中,当所述形态特征为ST波段特征时,所述获取所述信号波形的形态特征,包括:
获取T波序列,以及,获取S波序列;
根据所述S波序列和所述T波序列,获取ST波段斜率特征和/或ST波段基线特征;
确定所述ST波段斜率特征和/或所述ST波段基线特征,作为所述ST波段特征。
在一个实施例中,所述根据所述S波序列和所述T波序列,获取ST波段斜率特征,包括:
获取所述S波序列的S波信号序列,以及,获取所述T波序列的T波信号序列;
计算所述T波信号序列与所述S波信号序列的差值,得到ST波段斜率分子,以及,计算所述T波序列与所述S波序列的差值,得到ST波段斜率分母;
计算所述ST波段斜率分子与所述ST波段斜率分母的比值,作为所述ST波段斜率;
根据所述ST波段斜率,获取所述ST波段斜率特征;所述ST波段斜率特征包括所述ST波段斜率、ST波段斜率中位数、ST波段斜率方差以及ST波段斜率长度比值中的至少一种。
在一个实施例中,当所述ST波段斜率特征为所述ST波段斜率方差时,所述根据所述ST波段斜率,获取所述ST波段斜率特征,包括:
计算所述ST波段斜率的均值,得到ST波段斜率均值,以及,计算所述ST波段斜率的标准差,得到所述ST波段斜率标准差;
计算所述ST波段斜率标准差与所述ST波段斜率均值的比值,得到所述ST波段斜率方差。
在一个实施例中,当所述ST波段斜率特征为所述ST波段斜率长度比值时,所述根据所述ST波段斜率,获取所述ST波段斜率特征,包括:
获取QRS波段序列;所述QRS波段序列为所述S波序列与所述Q波序列的差值;
计算所述QRS波段序列的长度,得到ST波段斜率长度分母;
当所述ST波段斜率小于零时,计算所述ST波段斜率的长度,得到ST波段斜率长度分子;
计算所述ST波段斜率长度分子与所述ST波段斜率长度分母的比值,得到所述ST波段斜率长度比值。
在一个实施例中,所述根据所述S波序列和所述T波序列,获取ST波段基线特征,包括:
获取所述S波序列的S波信号序列;
当所述S波信号序列小于Q波信号序列时,确定SQ波向量;
提取所述SQ波向量中非零元素的位置索引,得到SQ波深度;
根据所述SQ波深度,获取所述ST波段基线特征;所述ST波段基线特征包括所述SQ波深度、S波序列深度、S波序列除余值、T波序列除余值以及Q波序列除余值中的至少一种。
在一个实施例中,当所述ST波段基线特征为所述S波序列深度时,所述根据所述SQ波深度,获取所述ST波段基线特征,包括:
计算所述S波序列的长度,得到S波序列深度分母,以及,计算所述SQ波深度的长度,得到S波序列深度分子;
计算所述S波序列深度分子与所述S波序列深度分母的比值,得到所述S波序列深度。
在一个实施例中,当所述ST波段基线特征为所述S波序列除余值时,所述根据所述SQ波深度,获取所述ST波段基线特征,包括:
计算所述S波序列与所述SQ波深度的比值,得到S波序列除数;
计算所述S波序列除数的余数,得到所述S波序列除余值。
在一个实施例中,当所述ST波段基线特征为所述T波序列除余值时,所述根据所述SQ波深度,获取所述ST波段基线特征,包括:
计算所述T波序列与所述SQ波深度的比值,得到T波序列除数;
计算所述T波序列除数的余数,得到所述T波序列除余值。
在一个实施例中,所述将所述形态特征输入至分类器,包括:
在所述信号波形中,提取RR间期序列;
根据所述RR间期序列,获取RR间期统计特征;所述RR间期统计特征包括所述RR间期序列的数量、均值、最大值、最小值、中位数、方差、偏度、峰度以及宽度中的至少一种;
将所述RR间期统计特征、所述ST波段斜率特征和所述ST波段基线特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期统计特征、所述ST波段斜率特征和所述ST波段基线特征进行分类的结果。
在一个实施例中,当所述RR间期统计特征为所述RR间期序列的宽度时,所述根据所述RR间期序列,获取RR间期统计特征,包括:
计算所述RR间期序列的最小值,得到RR间期序列最小值,以及,计算所述RR间期序列的最大值,得到RR间期序列最大值;
计算所述RR间期序列最大值与所述RR间期序列最小值的差值,作为所述RR间期序列的宽度。
在一个实施例中,当所述形态特征为P波变化特征时,所述获取所述信号波形的形态特征,包括:
获取P波序列;
计算所述P波序列的最小值,得到P波序列最小值,以及,计算所述P波序列的最大值,得到P波序列最大值;
计算所述P波序列最大值与所述P波序列最小值的差值,得到P波序列差值;
计算所述P波序列差值的最大值,得到P波变化特征分子,以及,计算所述P波序列的最大值,得到P波变化特征分母;
计算所述P波变化特征分子与所述P波变化特征分母的比值,作为所述P波变化特征。
在一个实施例中,所述将所述形态特征输入至分类器,包括:
在所述信号波形中,提取QRS间期信号序列;
根据所述QRS间期信号序列,获取QRS间期统计特征;所述QRS间期统计特征包括所述QRS间期信号序列的均值、标准差、方差、核密度、峰度以及偏度中的至少一种;
将所述QRS间期统计特征和所述P波变化特征输入至所述分类器;所述分类结果为所述分类器根据所述QRS间期统计特征和所述P波变化特征进行分类的结果。
在一个实施例中,所述获取所述信号波形的深层特征,包括:
采用所述信号波形训练基于深度学习的节律特征检测器;
获取训练后的节律特征检测器的最优网络结构;
获取所述最优网络结构的全连接层数据,作为所述深层特征。
在一个实施例中,所述节律特征检测器包括卷积神经网络检测器、循环神经网络检测器、卷积循环神经网络检测器、深度神经网络检测器中的至少一种。
在一个实施例中,所述分类器包括K-近邻(KNN)分类器、adboost分类器、xgboost分类器、逻辑回归分类器、支持向量机分类器中的至少一种。
在一个实施例中,所述获取所述分类器输出的分类结果,得到所述心电信号的信号类型,包括:
获取所述K-近邻(KNN)分类器、所述adboost分类器、所述xgboost分类器、所述逻辑回归分类器或所述支持向量机分类器的分类结果;
将所述分类结果输入至多层感知器;
获取所述多层感知器的多模型分类结果,得到所述心电信号的信号类型;所述多模型分类结果为所述多层感知器根据所述形态特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
第二方面,提供了一种心电信号检测分类装置,包括:
波形提取模块,用于在心电信号中,提取信号波形;
特征获取模块,用于获取所述信号波形的形态特征和深层特征;所述形态特征包括TR波幅度差特征、PR波数量比值特征、ST波段特征、P波变化特征中的任意一种;所述深层特征包括深度特征和层次特征;
特征输入模块,用于将所述形态特征和所述深层特征输入至分类器
分类模块,用于获取所述分类器输出的分类结果,得到所述心电信号的信号类型;所述分类结果为所述分类器根据所述形态特征和所述深层特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
第三方面,提供了一种电子设备,其特征在于,包括:存储器、以及一个或多个处理器;
所述存储器,用于存储一个或多个程序;
在心电信号中,提取信号波形;
获取所述信号波形的形态特征和深层特征;所述形态特征包括TR波幅度差特征、PR波数量比值特征、ST波段特征、P波变化特征中的任意一种;所述深层特征包括深度特征和层次特征;
将所述形态特征和所述深层特征输入至分类器;
获取所述分类器输出的分类结果,得到所述心电信号的信号类型;所述分类结果为所述分类器根据所述形态特征和所述深层特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
第四方面,提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
在心电信号中,提取信号波形;
获取所述信号波形的形态特征和深层特征;所述形态特征包括TR波幅度差特征、PR波数量比值特征、ST波段特征、P波变化特征中的任意一种;所述深层特征包括深度特征和层次特征;
将所述形态特征和所述深层特征输入至分类器;
获取所述分类器输出的分类结果,得到所述心电信号的信号类型;所述分类结果为所述分类器根据所述形态特征和所述深层特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
上述心电信号检测分类方法、装置、电子设备和存储介质,首先在心电信号中提取信号波形,然后获取信号波形的形态特征和深层特征,将形态特征和深层特征输入分类器,得到分类器根据形态特征和深层特征进行分类的分类结果作为心电信号的信号类型,因为形态特征和深层特征可以反映出各种类型的异常节律,采用形态特征和深层特征进行信号类型分类,可以更准确地识别出各种类型的异常节律,避免了将心动过速,心动过缓,心律失常等非房颤类型的异常节律错误分类为房颤类型的情况,提升了心电信号分类的准确性。
附图说明
图1是本申请实施例一提供的一种心电信号分类方法的流程图;
图2是一个实施例的一种原始心电信号的示意图;
图3是一种P波、QRS波、T波、U波的信号波形的示意图;
图4是一个实施例的一种卷积神经网络模型的架构示意图;
图5是一个实施例的一种节律特征检测器结构的结构示意图;
图6是一个实施例的一种卷积神经网络分类器的结构示意图;
图7是一个实施例的一种节律特征分类融合网络结构的结构示意图;
图8是本申请实施例二提供的一种心电信号检测分类装置的结构示意图;
图9是本申请实施例三提供的一种电子设备的结构示意图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
实施例一
图1是本申请实施例一提供的一种心电信号检测分类方法的流程图,该心电信号检测分类方法,可以具体包括以下步骤:
步骤S110,在心电信号中,提取信号波形。
具体实现中,可以利用多通道同步数据,采集人体心脏信号、背景噪声和心电信号。更具体地,首先,可以通过心电导联和传感器采集心电信号,通过模拟电路对采集的心电信号进行阻抗匹配、过滤、放大等处理。然后,由模数转换器,将人体生理参数的模拟信号转化为数字信号。然后,通过低通滤波技术获取滤波后的心电信号。最后,通过小波变换技术从滤波后的心电信号中提取出信号波形。
图2是一个实施例的一种原始心电信号的示意图。如图所示,采集得到的心电信号存在各种噪声,而且信号的波形粗糙和不光滑,导致难以从QRS波中提取出能反映房颤的特征。
图3是一种P波、QRS波、T波、U波的信号波形的示意图。如图所示,以时间t作为X轴,电压V作为Y轴,构建一个XY坐标系,根据心电信号中在每个时间点tn对应的电压Vn,生成心电信号的信号波形。
步骤S120,获取所述信号波形的形态特征和深层特征;所述形态特征包括TR波幅度差特征、PR波数量比值特征、ST波段特征、P波变化特征中的任意一种;所述深层特征包括深度特征和层次特征。
其中,TR波幅度差特征可以为反映TR波的幅度差相关的特征。例如,T波信号序列和R波信号序列的信号比值。
其中,PR波数量比值特征可以为反映P波数量和R波数量的比值相关特征,例如,P波数量与R波数量比值斜率相关的特征。
其中,ST波段特征可以为反映ST波段相关的特征,例如,ST波段斜率相关的特征、ST波段基线相关的特征。
其中,P波变化特征可以为反映P波的变化相关的特征。
其中,深度特征和层次特征可以反映经过深度学习网络所提取、针对房颤节律、非房颤异常节律、正常窦性心律和噪声相关的特征。
具体实现中,可以获取信号波形的TR波幅度差特征、PR波数量比值特征、ST波段特征和P波变化特征,作为信号波形的形态特征。
首先,获取信号波形的TR波幅度差特征的方式可以包括:令X_P、X_Q、X_R、X_S和X_T分别表示P、Q、R、S和T波的序列,ecg表示心电信号序列,则TR波的幅度差计算公式如下所示:
tr_amp=ecg(X_T)./ecg_f(X_R)
tr_amp_med=median(tr_amp)
其次,获取信号波形的PR波数量比值特征的方式可以包括:
(1)令X_P、X_Q、X_R、X_S和X_T分别表示P、Q、R、S和T波的序列,ecg表示心电信号序列,则P波数量与R波数量及其比值计算公式如下所示:
Num_R=Num(X_R)=sum(X_R)
Num_P=Num(X_P)=sum(X_P)
PRc=Num_P/Num_R
(2)计算P波数量与R波数量比值斜率相关的特征,即分别标记为med_prc和var_prc,计算公式如下所示:
med_prc=median(PRc)
var_prc=std(PRc)/mean(PRc)
其中,median、std、mean和var分别为计算序列中位数、标准差、均值和方差的算子。
然后,获取信号波形的ST波段特征的方式可以包括:
(1)令X_P、X_Q、X_R、X_S和X_T分别表示P、Q、R、S和T波的序列,ecg表示心电信号序列,则ST波段的斜率计算公式如下所示:
slope_ST(i)=(ecg(X_T(i))-ecg(X_S(i)))/(X_T(i)-X_S(i))
(2)计算ST波段斜率相关的特征,即分别标记为med_st、var_st和num_neg_st,计算公式如下所示:
med_st=median(slope_ST)
var_st=std(slope_ST)/mean(slope_ST)
num_neg_st=length(slope_ST<0)/length(X_QRS)
X_QRS=X_S-X_Q
其中,median、std、mean、length分别为计算序列中位数、标准差、均值和长度的算子。
(3)计算ST波段基线相关的特征,即分别标记为deep_s1、deep_s、mod_S、mod_T和mod_Q,计算公式如下所示:
deep_s1=find(ecg(X_S)<ecg(X_Q))
deep_s=length(deep_s1)/length(X_S)
mod_S=X_S(deep_s1)
mod_T=X_T(deep_s1)
mod_Q=X_Q(deep_s1)
其中,ecg表示心电信号序列,find()函数的基本功能是返回向量或者矩阵中不为0的元素的位置索引,length为计算序列长度的算子。
最后,获取信号波形的P波变化特征的方式可以包括:设P(i,j)为P波序列,计算P波序列中最大值与最小值的差PD,计算公式如下所示:
其中,i为P波序列的样本数,j表示P波序列的第j个样本。
另需说明的是,对深度特征和层次特征的获取,可以是基于不同形态特征的获取条件,进而采用不同检测器,如卷积神经网络检测器、循环神经网络检测器、卷积循环神经网络检测器、深度神经网络检测器等检测提取所得。
步骤S130,将所述形态特征和所述深层特征输入至分类器。
具体实现中,在得到形态特征和所述深层特征之后,可以将形态特征和所述深层特征输入至分类器。
需要说明的是,考虑到RR间期不规则是房颤检测的重要指标,因此提取RR间期相关的AF先验特征、心电信号中的统计特征,并将其与TR波幅度差特征、PR波数量比值特征、ST波段特征和P波变化特征等形态特征中的一个或多个同时输入至分类器,以获取分类结果显得尤为重要。其中,AF先验特征可以包括RR间期差的原始计数、改进近似熵、改进样本熵、RR间期概率密度相关指标、RR间期差概率密度相关指标等;统计特征包括时频域上形态的上升和下降状态特征、R波的数量、平均值、最大值、最小值、中位数和方差等。
除此之外,输入至分类器的特征还包括深度特征和层次特征,而深度特征和层次特征的提取首先需基于一个预设的学习模型,在学习模型中输入单个的心拍,即一个一个的RR间期信号,输出结果即为始终节律信号:房颤节律、非房颤异常节律、正常窦性心律和噪声。其中,预设的学习模型可以是卷积神经网络模型、循环神经网络模型、卷积循环神经网络模型或深度神经网络模型等。
在一个具体实施例中,若使用卷积神经网络模型(CNN)提取深度特征和层次特征,则提取步骤包括:CNN架构如图4所示,包括输入层(Input vector)、卷积层(convolutionlayer)、全连接层(fully connected NN)和输出层(Output layer);(1)卷积层具有k个滤波器或大小为q×n的内核,其中q小于长度输入向量,每个滤波器与输入进行卷积,然后进行非线性激活(sigmoid)产生大小为(m-q+1)×1的k个特征图,然后将多个滤波器的输出堆叠在一起形成单个特征向量。(2)全连接层采用有softmax激活函数,采用基于随机梯度下降法的交叉熵损失函数进行卷积层和全连接层的权重优化,训练得到最佳网络结构后,将网络保存下来,并将全连接层的数据存储下来,作为心电信号的深度特征和层次特征,并与上述实施例中所提特征一起,作为心电四种节律识别分类器的输入特征。
步骤S140,获取所述分类器输出的分类结果,得到所述心电信号的信号类型;所述分类结果为所述分类器根据所述形态特征和所述深层特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
其中,分类器可以用于对心电信号进行分类。分类器可以为K-近邻(KNN,k-NearestNeighbor)分类器、adboost分类器、xgboost(Extreme Gradient Boosting)分类器、逻辑回归分类器、支持向量机分类器中的至少一种。
具体实现中,可以将形态特征和所述深层特征输入至经过训练的分类器,分类器可以根据形态特征和所述深层特征,统计心电信号分别归类为房颤节律、非房颤异常节律、正常窦性心律和噪声的概率分布,并根据概率分布输出分类结果,将该分类结果,作为心电信号的信号类型。
本申请实施例提供的技术方案,首先在心电信号中提取信号波形,然后获取信号波形的形态特征和深层特征,将形态特征和深层特征输入分类器,得到分类器根据形态特征和深层特征进行分类的分类结果作为心电信号的信号类型,因为形态特征和深层特征可以反映出各种类型的异常节律,采用形态特征和深层特征进行信号类型分类,可以更准确地识别出各种类型的异常节律,避免了将心动过速,心动过缓,心律失常等非房颤类型的异常节律错误分类为房颤类型的情况,提升了心电信号分类的准确性。
在一个实施例中,当所述形态特征为TR波幅度差特征时,所述步骤S120,包括:
获取R波信号序列,以及,获取T波信号序列;计算所述T波信号序列与所述R波信号序列的比值,得到TR波序列比值;计算所述TR波序列比值的中位数,得到TR波序列比值中位数;确定所述TR波序列比值中位数和所述TR波序列比值,作为所述TR波幅度差特征。
具体实现中,TR波幅度差特征作为形态特征中的一种,其获取方式包括:令X_P、X_Q、X_R、X_S和X_T分别表示P、Q、R、S和T波的序列,ecg表示心电信号序列,则TR波幅度差特征的计算公式如下所示:
tr_amp=ecg(X_T)./ecg_f(X_R)
tr_amp_med=median(tr_amp)
其中,ecg(X_T)表示为T波信号序列,ecg_f(X_R)表示为R波信号序列,则tr_amp表示为TR波序列比值、tr_amp_med表示为TR波序列比值中位数。
在一个实施例中,所述将所述形态特征输入至分类器,包括:
在所述信号波形中,提取RR间期差序列;根据所述RR间期差序列,获取RR间期差第一先验特征;所述RR间期差第一先验特征包括所述RR间期差序列的原始计数、近似熵、样本熵中的至少一种;将所述RR间期差第一先验特征和所述TR波幅度差特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第一先验特征和所述TR波幅度差特征进行分类的结果。
具体实现中,若当前获取到的形态特征为TR波幅度差特征,则还需获取与该TR波幅度差特征相关联的第一先验特征,第一先验特征的获取步骤包括:令X_P、X_Q、X_R、X_S和X_T分别表示P、Q、R、S和T波的序列,ecg表示心电信号序列,第一先验特征包括RR间期差的原始计数,改进近似熵、改进样本熵。其中,(1)计算RR间期差的改进近似熵方式可以是:
步骤(1),选取一定长度N的RR间期,并别除一定个数Q的最大和最小RR间期,以减少异位心搏的干扰。
步骤(2),对剩下的RR间期进行直方图统计,具体地,求解最大、最小RR间期以获得RR间期范围,在RR间期范围内将所有RR间期分按序号顺序组成一组m维向量X(i),X(i)=[xi,xi+1,...,xi+m-1];i=1,2,...,N-m+1。
步骤(3),定义X(i)和X(j)间的距离d[X(i),X(j)]为两者对应元素中差值最大的一个,如下述公式所示:
其中,d[X(i),X(j)]表示每一个i值对应的X(i)与其余向量X(j)(j≠i,j=1,2,...,N-m+1)间的距离。
其中,r为定值,通常r=(0.1~0.25)SD(SD为序列{x1,x2,...,xN}的标准差)。
步骤(5),计算公式如下所示:
步骤(7),计算RR间期差的近似熵AE,表示为AE△RR,如下述公式所示:
步骤(8),为了避免由固定值r引起的熵估计的不可靠性,提出了一种改进的基于密度的熵估计方法,将RR间期差的近似熵AE与固定值r做比,表示为rAE△RR,如下述公式所示:
(2)计算RR间期差的改进样本熵方式可以是:
步骤(1),选取一定长度P的RR间期,并别除一定个数Q的最大和最小RR间期,以减少异位心搏的干扰。
步骤(2),对剩下的RR间期进行直方图统计,具体地,求解最大、最小RR间期以获得RR间期范围,在RR间期范围内将所有RR间期分按序号顺序组成一组m维向量X(i),X(i)=[xi,xi+1,...,xi+m-1];i=1,2,...,N-m+1。
步骤(3),定义X(i)和X(j)间的距离d[X(i),X(j)]为两者对应元素中差值最大的一个,如下述公式所示:
其中,d[X(i),X(j)]表示每一个i值对应的X(i)与其余向量X(j)(j≠i,j=1,2,...,N-m+1)间的距离。
其中,r为定值,通常r=(0.1~0.25)SD(SD为序列{x1,x2,...,xN}的标准差)。
步骤(5),计算公式如下所示:
步骤(7),计算RR间期差的样本熵SAE,表示为SAE△RR,如下述公式所示:
SAE(m,r,N)=-ln[Bm+1(r)/Bm(r)]
步骤(8),为了避免由固定值r引起的熵估计的不可靠性,提出了一种改进的基于密度的熵估计方法,将RR间期差的样本熵SAE与固定值r做比,表示为rSAE△RR,如下述公式所示:
在一个实施例中,当所述RR间期差第一先验特征为所述RR间期差序列的原始计数时,所述根据所述RR间期差序列,获取RR间期差第一先验特征,包括:
计算所述RR间期差序列的绝对值,得到RR间期差序列绝对值;当所述RR间期差序列绝对值小于等于预设的常量定值时,计算所述RR间期差序列之和,得到RR间期差序列和值;当所述RR间期差序列和值等于二时,确定所述RR间期差序列和值,作为所述RR间期差序列的原始计数。
具体实现中,计算RR间期差的原始计数公式如下:
os=sum(|dRR|<=OCmask)
OriginCount=sum(os==2)
其中,dRR表示为RR间期差序列,OCmask是一个常量定值,在本实施例中OCmask=0.02,sum为计算序列总和的算子,OriginCount表示为RR间期差的原始计数。
在一个实施例中,当所述形态特征为PR波数量比值特征时,所述步骤S120,包括:
获取R波序列数量,以及,获取P波序列数量;所述P波序列数量为P波序列的总和;所述R波序列数量为R波序列的总和;计算所述P波序列数量与所述R波序列数量的比值,得到PR波数量比值;根据所述PR波数量比值,确定所述PR波数量比值特征;所述PR波数量比值特征包括所述PR波数量比值、PR波数量比值中位数、PR波数量比值标准差、PR波数量比值均值以及PR波数量比值方差中的至少一种。
具体实现中,PR波数量比值特征作为形态特征中的一种,其获取方式包括:令X_P、X_Q、X_R、X_S和X_T分别表示P、Q、R、S和T波的序列,ecg表示心电信号序列,则PR波数量比值特征的计算公式如下所示:
Num_R=Num(X_R)=sum(X_R)
Num_P=Num(X_P)=sum(X_P)
PRc=Num_P/Num_R
其中,sum(X_R)表示为R波序列数量,sum(X_P)表示为P波序列数量,则PRc表示为PR波数量比值。
在一个实施例中,当所述PR波数量比值特征为所述PR波数量比值方差时,所述根据所述PR波数量比值,确定所述PR波数量比值特征,包括:
计算所述PR波数量比值的均值,得到所述PR波数量比值均值,以及,计算所述PR波数量比值的标准差,得到所述PR波数量比值标准差;计算所述PR波数量比值标准差与所述PR波数量比值均值的比值,得到所述PR波数量比值方差。
具体实现中,由于PR波数量比值特征包括PR波数量比值、PR波数量比值中位数、PR波数量比值标准差、PR波数量比值均值以及PR波数量比值方差中的至少一种,若当前获取到的形态特征为PR波数量比值特征,则计算PR波数量比值方差的公式如上述实施例一中的计算公式。
在一个实施例中,所述将所述形态特征输入至分类器,包括:
在所述信号波形中,提取RR间期序列和/或RR间期差序列;根据所述RR间期序列和/或所述RR间期差序列,获取RR间期差第二先验特征;将所述RR间期差第二先验特征和所述PR波数量比值特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第二先验特征和所述PR波数量比值特征进行分类的结果。
具体实现中,若当前获取到的形态特征为PR波数量比值特征,则还需获取与该PR波数量比值特征相关联的第二先验特征,第二先验特征包括RR间期密度偏度、RR间期密度峰度、RR间期序列核密度、RR间期差密度偏度、RR间期差密度峰度、RR间期差序列核密度。
在一个实施例中,所述根据所述RR间期序列,获取RR间期差第二先验特征,包括:
计算所述RR间期序列的核密度,得到RR间期序列核密度;计算所述RR间期序列核密度的峰度值,得到RR间期密度峰度,以及,计算所述RR间期序列核密度的偏度值,得到RR间期密度偏度;确定所述RR间期密度偏度、所述RR间期密度峰度以及所述RR间期序列核密度中的至少一种,作为所述RR间期差第二先验特征。
具体实现中,RR间期密度偏度、RR间期密度峰度以及RR间期序列核密度的计算公式如下:
[f,xi]=ksdensity(RR)
kurt_f=kurtosis(f)
skew_f=skewness(f)
其中,ksdensity为计算序列核密度估计的算子,用于求样本观测向量x的核密度估计,xi是在x取值范围内等间隔选取的100个点构成的向量,f是与xi相对应的核密度估计值向量,所用的核函数是Gaussian核函数,此外,length、skewness、kurtosis分别为计算序列长度、偏度和峰度的算子。
在一个实施例中,所述根据所述RR间期序列和/或所述RR间期差序列,获取RR间期差第二先验特征,包括:
计算所述RR间期差序列的核密度,得到RR间期差序列核密度;计算所述RR间期差序列和密度的峰度值,得到RR间期差密度峰度,以及,计算所述RR间期差序列核密度的偏度值,得到RR间期差密度偏度;确定所述RR间期差密度偏度、所述RR间期差密度峰度以及所述RR间期差序列核密度中的至少一种,作为所述RR间期差第二先验特征。
具体实现中,RR间期差密度偏度、RR间期差密度峰度以及RR间期差序列核密度的计算公式如下:
[f1,xi1]=ksdensity(dRR)
kurt_f1=kurtosis(f1)
skew_f1=skewness(f1)
其中,ksdensity为计算序列核密度估计的算子,用于求样本观测向量x的核密度估计,xi是在x取值范围内等间隔选取的100个点构成的向量,f是与xi相对应的核密度估计值向量,所用的核函数是Gaussian核函数,此外,length、skewness、kurtosis分别为计算序列长度、偏度和峰度的算子。
在一个实施例中,当所述形态特征为ST波段特征时,所述步骤S120,包括:
获取T波序列,以及,获取S波序列;根据所述S波序列和所述T波序列,获取ST波段斜率特征和/或ST波段基线特征;确定所述ST波段斜率特征和/或所述ST波段基线特征,作为所述ST波段特征。
具体实现中,ST波段特征作为形态特征中的一种,其计算公式包括上述实施例一中所述ST波段斜率slope_ST、ST波段斜率相关特征、ST波段基线相关特征的计算公式。
在一个实施例中,所述根据所述S波序列和所述T波序列,获取ST波段斜率特征,包括:
获取所述S波序列的S波信号序列,以及,获取所述T波序列的T波信号序列;计算所述T波信号序列与所述S波信号序列的差值,得到ST波段斜率分子,以及,计算所述T波序列与所述S波序列的差值,得到ST波段斜率分母;计算所述ST波段斜率分子与所述ST波段斜率分母的比值,作为所述ST波段斜率;根据所述ST波段斜率,获取所述ST波段斜率特征;所述ST波段斜率特征包括所述ST波段斜率、ST波段斜率中位数、ST波段斜率方差以及ST波段斜率长度比值中的至少一种。
具体实现中,ST波段斜率slope_ST、ST波段斜率中位数med_st、ST波段斜率方差var_st以及ST波段斜率长度比值num_neg_st的计算公式如上述实施例一中计算ST波段特征的公式。
在一个实施例中,当所述ST波段斜率特征为所述ST波段斜率方差时,所述根据所述ST波段斜率,获取所述ST波段斜率特征,包括:
计算所述ST波段斜率的均值,得到ST波段斜率均值,以及,计算所述ST波段斜率的标准差,得到所述ST波段斜率标准差;计算所述ST波段斜率标准差与所述ST波段斜率均值的比值,得到所述ST波段斜率方差。
具体实现中,ST波段斜率方差var_st的计算公式如上述实施例一中计算ST波段斜率相关特征的公式。
在一个实施例中,当所述ST波段斜率特征为所述ST波段斜率长度比值时,所述根据所述ST波段斜率,获取所述ST波段斜率特征,包括:
获取QRS波段序列;所述QRS波段序列为所述S波序列与所述Q波序列的差值;计算所述QRS波段序列的长度,得到ST波段斜率长度分母;当所述ST波段斜率小于零时,计算所述ST波段斜率的长度,得到ST波段斜率长度分子;计算所述ST波段斜率长度分子与所述ST波段斜率长度分母的比值,得到所述ST波段斜率长度比值。
具体实现中,ST波段斜率长度比值num_neg_st的计算公式如上述实施例一中计算ST波段斜率相关特征的公式。
在一个实施例中,所述根据所述S波序列和所述T波序列,获取ST波段基线特征,包括:
获取所述S波序列的S波信号序列;当所述S波信号序列小于Q波信号序列时,确定SQ波向量;提取所述SQ波向量中非零元素的位置索引,得到SQ波深度;根据所述SQ波深度,获取所述ST波段基线特征;所述ST波段基线特征包括所述SQ波深度、S波序列深度、S波序列除余值、T波序列除余值以及Q波序列除余值中的至少一种。
具体实现中,SQ波深度deep_s1、S波序列深度deep_s、S波序列除余值mod_S、T波序列除余值mod_T以及Q波序列除余值mod_Q的计算公式如上述实施例一中计算ST波段基线相关特征的公式。
在一个实施例中,当所述ST波段基线特征为所述S波序列深度时,所述根据所述SQ波深度,获取所述ST波段基线特征,包括:
计算所述S波序列的长度,得到S波序列深度分母,以及,计算所述SQ波深度的长度,得到S波序列深度分子;计算所述S波序列深度分子与所述S波序列深度分母的比值,得到所述S波序列深度。
具体实现中,S波序列深度的计算公式如上述实施例一中所述计算公式。
在一个实施例中,当所述ST波段基线特征为所述S波序列除余值时,所述根据所述SQ波深度,获取所述ST波段基线特征,包括:
计算所述S波序列与所述SQ波深度的比值,得到S波序列除数;计算所述S波序列除数的余数,得到所述S波序列除余值。
具体实现中,S波序列除余值的计算公式如上述实施例一中所述计算公式。
在一个实施例中,当所述ST波段基线特征为所述T波序列除余值时,所述根据所述SQ波深度,获取所述ST波段基线特征,包括:
计算所述T波序列与所述SQ波深度的比值,得到T波序列除数;计算所述T波序列除数的余数,得到所述T波序列除余值。
具体实现中,T波序列除余值的计算公式如上述实施例一中所述计算公式。
在一个实施例中,所述将所述形态特征输入至分类器,包括:
在所述信号波形中,提取RR间期序列;根据所述RR间期序列,获取RR间期统计特征;所述RR间期统计特征包括所述RR间期序列的数量、均值、最大值、最小值、中位数、方差、偏度、峰度以及宽度中的至少一种;将所述RR间期统计特征、所述ST波段斜率特征和所述ST波段基线特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期统计特征、所述ST波段斜率特征和所述ST波段基线特征进行分类的结果。
具体实现中,若当前获取到的形态特征为ST波段特征,则还需获取与该ST波段特征相关联的RR间期统计特征,RR间期统计特征中RR间期序列数量Num_RR、RR间期序列均值Mean_RR、RR间期序列最大值Max_RR、RR间期序列最小值Min_RR、RR间期序列中位数Median_RR、RR间期序列方差Var_RR、RR间期序列偏度skewness_RR、RR间期序列峰度kurtosis_RR以及RR间期序列宽度range_RR的计算公式分别如下:
Num_RR=Num(RR)=length(RR)
Mean_RR=Mean(RR)
Max_RR=Max(RR)
Min_RR=Min(RR)
Median_RR=Median(RR)
Var_RR=Var(X_RR)
skewness_RR=skewness(RR)
kurtosis_RR=kurtosis(RR)
range_RR=Max(RR)-Min(RR)
在一个实施例中,当所述RR间期统计特征为所述RR间期序列的宽度时,所述根据所述RR间期序列,获取RR间期统计特征,包括:
计算所述RR间期序列的最小值,得到RR间期序列最小值,以及,计算所述RR间期序列的最大值,得到RR间期序列最大值;计算所述RR间期序列最大值与所述RR间期序列最小值的差值,作为所述RR间期序列的宽度。
具体实现中,RR间期序列宽度range_RR的计算公式如上述实施例所示。
在一个实施例中,当所述形态特征为P波变化特征时,所述步骤S120,包括:
获取P波序列;计算所述P波序列的最小值,得到P波序列最小值,以及,计算所述P波序列的最大值,得到P波序列最大值;计算所述P波序列最大值与所述P波序列最小值的差值,得到P波序列差值;计算所述P波序列差值的最大值,得到P波变化特征分子,以及,计算所述P波序列的最大值,得到P波变化特征分母;计算所述P波变化特征分子与所述P波变化特征分母的比值,作为所述P波变化特征。
具体实现中,P波变化特征作为形态特征中的一种,P波变化特征的相关计算公式如实施例一种所述公式。
在一个实施例中,所述将所述形态特征输入至分类器,包括:
在所述信号波形中,提取QRS间期信号序列;根据所述QRS间期信号序列,获取QRS间期统计特征;所述QRS间期统计特征包括所述QRS间期信号序列的均值、标准差、方差、核密度、峰度以及偏度中的至少一种;将所述QRS间期统计特征和所述P波变化特征输入至所述分类器;所述分类结果为所述分类器根据所述QRS间期统计特征和所述P波变化特征进行分类的结果。
具体实现中,若当前获取到的形态特征为P波变化特征,则还需获取与该P波变化特征相关联的QRS间期统计特征,QRS间期统计特征的获取步骤包括:令X_P、X_Q、X_R、X_S和X_T分别表示P、Q、R、S和T波的序列,ecg表示心电信号序列,QRS表示QRS间期序列,QRS=X_Q-X_S,则QRS间期信号序列均值、QRS间期信号序列标准差、QRS间期信号序列方差、QRS间期信号序列核密度、QRS间期信号序列峰度以及QRS间期信号序列偏度的计算公式如下所示:
amplitude=ecg(QRS)
var_amp=std(amplitude)/mean(amplitude)
[f3,xi3]=skewness(amplitude)
kurt_f3=kurtosis(f3)
skew_f3=skewness(f3)
其中,length、mean、median、var、skewness、kurtosis分别为计算序列长度、均值、中位数、方差、偏度和峰度的算子,ksdensity为计算序列核密度估计的算子,用于求样本观测向量x的核密度估计,xi是在x取值范围内等间隔选取的100个点构成的向量,f是与xi相对应的核密度估计值向量,所用的核函数是Gaussian核函数。
在一个实施例中,所述获取所述信号波形的深层特征,包括:
采用所述信号波形训练基于深度学习的节律特征检测器;获取训练后的节律特征检测器的最优网络结构;获取所述最优网络结构的全连接层数据,作为所述深层特征。
具体实现中,可以采用卷积神经网络(CNNs)、卷积循环神经网络(CRNN)等提取深度特征和层次特征。例如,将单个心拍的信号输入至节律特征检测器,即一个一个的RR间期信号,节律特征检测器输出房颤节律、非房颤异常节律、正常窦性心律和噪声的分类结果,基于分类结果进行机器训练。
图5所示为一种节律特征检测器结构的结构示意图,基于图所示结构与下表1所示内容,即可分类输出始终节律信号。
表1深度特征提取器的网络结构
其中,CRNN模型有两个输入:一个是心电信号,另一个是心率,两个输入均经过3次一维卷积之后,输入循环神经网络RNN。首先,输入被发送到三个1D卷积层,在每个卷积层之间,通过最大池化(max-pooling)减少冗余和降低计算成本减少,除了max-pooling层,还添加dropout层来减少过度拟合。接下来,使用RNN提取层次特征,即LSTM(Long Short-TermMemory)层,与传统RNN不同的是,本实施例采用2层的LSTM的优势在于:(1)LSTM可以处理长时依赖并解决消失的梯度问题。(2)2个连续的LSTM层增加时间依赖性的长度。
又例如,设计一个16层卷积神经网络(CNN)进行心电节律分类,其输入为单个的心拍,即一个一个的RR间期信号,网络结构如图6所示,输出为四种节律信号。本实施例提出设计了一个16层CNN进行心电节律识别,16个卷积层顺序连接,采用batch normalization进行数据标准化处理,避免梯度爆炸问题,采用ReLU激活函数进行激活,采用Dropout层来降低CNN对训练数据的过度拟合。卷积层是提取深度特征的核心层,通过为CNN提供标记数据,该模型能够学习表征不同的重要特征类。
节律特征检测器经过训练后,可以从节律特征检测器中得到最优的网络结构,将该最优网络结构的全连接层数据保存,作为深度特征和层次特征。将深度特征、层次特征、形态统计特征一起输入至分类器的softmax层,供分类器根据深度特征、层次特征、形态统计特征进行分类。更具体地,可通过K-近邻(KNN)分类器对不同心电节律进行分类。
在一个实施例中,所述节律特征检测器包括卷积神经网络检测器、循环神经网络检测器、卷积循环神经网络检测器、深度神经网络检测器中的至少一种。
在一个实施例中,所述分类器包括K-近邻(KNN)分类器、adboost分类器、xgboost分类器、逻辑回归分类器、支持向量机分类器中的至少一种。
在一个实施例中,所述获取所述分类器输出的分类结果,得到所述心电信号的信号类型,包括:
获取所述K-近邻(KNN)分类器、所述adboost分类器、所述xgboost分类器、所述逻辑回归分类器或所述支持向量机分类器的分类结果;将所述分类结果输入至多层感知器;获取所述多层感知器的多模型分类结果,得到所述心电信号的信号类型;所述多模型分类结果为所述多层感知器根据所述形态特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
具体实现中,本实施例提出在分别获取到多个分类器求解结果之后,进行多分类器结果融合,即从心电信号中提取不同的特征集,分别输入到不同的分类器模型中,最后将各个模型的分类结果输入到具有softmax输出层的多层感知器(MLP),进行最终的分类结果判断。具体的,可参阅图7,图7为本实施例中四种心电节律分类融合的网络结构示意图。
本申请实施例提供的技术方案,由于深度特征和层次特征可以准确反映各种类型的异常节律,根据深度特征、层次特征和形态特征共同进行信号类型分类,可以更准确地识别出各种类型的异常节律,避免了将心动过速,心动过缓,心律失常等的非房颤类型的异常节律错误地分类为房颤类型的异常节律,进一步提升了心电信号分类的准确性。
应该理解的是,虽然图1流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图1中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
实施例二
图8是本申请实施例二提供的一种心电信号检测分类装置的结构示意图。参考图8,本实施例提供的心电信号检测分类装置具体包括:波形提取模块810、特征获取模块820、特征输入模块830和分类模块840,其中:
波形提取模块810,用于在心电信号中,提取信号波形;
特征获取模块820,用于获取所述信号波形的形态特征和深层特征;所述形态特征包括TR波幅度差特征、PR波数量比值特征、ST波段特征、P波变化特征中的任意一种;所述深层特征包括深度特征和层次特征;
特征输入模块830,用于将所述形态特征和所述深层特征输入至分类器;
分类模块840,用于获取所述分类器输出的分类结果,得到所述心电信号的信号类型;所述分类结果为所述分类器根据所述形态特征和所述深层特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
本申请实施例提供的技术方案,首先在心电信号中提取信号波形,然后获取信号波形的形态特征和深层特征,将形态特征和深层特征输入分类器,得到分类器根据形态特征和深层特征进行分类的分类结果作为心电信号的信号类型,因为形态特征和深层特征可以反映出各种类型的异常节律,采用形态特征和深层特征进行信号类型分类,可以更准确地识别出各种类型的异常节律,避免了将心动过速,心动过缓,心律失常等非房颤类型的异常节律错误分类为房颤类型的情况,提升了心电信号分类的准确性。
在一个实施例中,当所述形态特征为TR波幅度差特征时,所述特征获取模块820,进一步用于获取R波信号序列,以及,获取T波信号序列;计算所述T波信号序列与所述R波信号序列的比值,得到TR波序列比值;计算所述TR波序列比值的中位数,得到TR波序列比值中位数;确定所述TR波序列比值中位数和所述TR波序列比值,作为所述TR波幅度差特征。
在一个实施例中,所述特征输入模块830,进一步用于在所述信号波形中,提取RR间期差序列;根据所述RR间期差序列,获取RR间期差第一先验特征;所述RR间期差第一先验特征包括所述RR间期差序列的原始计数、近似熵、样本熵中的至少一种;将所述RR间期差第一先验特征和所述TR波幅度差特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第一先验特征和所述TR波幅度差特征进行分类的结果。
在一个实施例中,当所述RR间期差第一先验特征为所述RR间期差序列的原始计数时,所述特征输入模块830,进一步用于计算所述RR间期差序列的绝对值,得到RR间期差序列绝对值;当所述RR间期差序列绝对值小于等于预设的常量定值时,计算所述RR间期差序列之和,得到RR间期差序列和值;当所述RR间期差序列和值等于二时,确定所述RR间期差序列和值,作为所述RR间期差序列的原始计数。
在一个实施例中,当所述形态特征为PR波数量比值特征时,所述特征获取模块820,进一步用于获取R波序列数量,以及,获取P波序列数量;所述P波序列数量为P波序列的总和;所述R波序列数量为R波序列的总和;计算所述P波序列数量与所述R波序列数量的比值,得到PR波数量比值;根据所述PR波数量比值,确定所述PR波数量比值特征;所述PR波数量比值特征包括所述PR波数量比值、PR波数量比值中位数、PR波数量比值标准差、PR波数量比值均值以及PR波数量比值方差中的至少一种。
在一个实施例中,当所述PR波数量比值特征为所述PR波数量比值方差时,所述特征获取模块820,进一步用于计算所述PR波数量比值的均值,得到所述PR波数量比值均,以及,计算所述PR波数量比值的标准差,得到所述PR波数量比值标准差;计算所述PR波数量比值标准差与所述PR波数量比值均值的比值,得到所述PR波数量比值方差。
在一个实施例中,所述特征输入模块830,进一步用于在所述信号波形中,提取RR间期序列和/或RR间期差序列;根据所述RR间期序列和/或所述RR间期差序列,获取RR间期差第二先验特征;将所述RR间期差第二先验特征和所述PR波数量比值特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第二先验特征和所述PR波数量比值特征进行分类的结果。
在一个实施例中,所述特征输入模块830,进一步用于计算所述RR间期序列的核密度,得到RR间期序列核密度;计算所述RR间期序列核密度的峰度值,得到RR间期密度峰度,以及,计算所述RR间期序列核密度的偏度值,得到RR间期密度偏度;确定所述RR间期密度偏度、所述RR间期密度峰度以及所述RR间期序列核密度中的至少一种,作为所述RR间期差第二先验特征。
在一个实施例中,所述特征输入模块830,进一步用于计算所述RR间期差序列的核密度,得到RR间期差序列核密度;计算所述RR间期差序列和密度的峰度值,得到RR间期差密度峰度,以及,计算所述RR间期差序列核密度的偏度值,得到RR间期差密度偏度;确定所述RR间期差密度偏度、所述RR间期差密度峰度以及所述RR间期差序列核密度中的至少一种,作为所述RR间期差第二先验特征。
在一个实施例中,当所述形态特征为ST波段特征时,所述特征获取模块820,进一步用于获取T波序列,以及,获取S波序列;根据所述S波序列和所述T波序列,获取ST波段斜率特征和/或ST波段基线特征;确定所述ST波段斜率特征和/或所述ST波段基线特征,作为所述ST波段特征。
在一个实施例中,所述特征获取模块820,进一步用于获取所述S波序列的S波信号序列,以及,获取所述T波序列的T波信号序列;计算所述T波信号序列与所述S波信号序列的差值,得到ST波段斜率分子,以及,计算所述T波序列与所述S波序列的差值,得到ST波段斜率分母;计算所述ST波段斜率分子与所述ST波段斜率分母的比值,作为所述ST波段斜率;根据所述ST波段斜率,获取所述ST波段斜率特征;所述ST波段斜率特征包括所述ST波段斜率、ST波段斜率中位数、ST波段斜率方差以及ST波段斜率长度比值中的至少一种。
在一个实施例中,当所述ST波段斜率特征为所述ST波段斜率方差时,所述特征获取模块820,进一步用于计算所述ST波段斜率的均值,得到ST波段斜率均值,以及,计算所述ST波段斜率的标准差,得到所述ST波段斜率标准差;计算所述ST波段斜率标准差与所述ST波段斜率均值的比值,得到所述ST波段斜率方差。
在一个实施例中,当所述ST波段斜率特征为所述ST波段斜率长度比值时,所述特征获取模块820,进一步用于获取QRS波段序列;所述QRS波段序列为所述S波序列与所述Q波序列的差值;计算所述QRS波段序列的长度,得到ST波段斜率长度分母;当所述ST波段斜率小于零时,计算所述ST波段斜率的长度,得到ST波段斜率长度分子;计算所述ST波段斜率长度分子与所述ST波段斜率长度分母的比值,得到所述ST波段斜率长度比值。
在一个实施例中,所述特征获取模块820,进一步用于获取所述S波序列的S波信号序列;当所述S波信号序列小于Q波信号序列时,确定SQ波向量;提取所述SQ波向量中非零元素的位置索引,得到SQ波深度;根据所述SQ波深度,获取所述ST波段基线特征;所述ST波段基线特征包括所述SQ波深度、S波序列深度、S波序列除余值、T波序列除余值以及Q波序列除余值中的至少一种。
在一个实施例中,当所述ST波段基线特征为所述S波序列深度时,所述特征获取模块820,进一步用于计算所述S波序列的长度,得到S波序列深度分母,以及,计算所述SQ波深度的长度,得到S波序列深度分子;计算所述S波序列深度分子与所述S波序列深度分母的比值,得到所述S波序列深度。
在一个实施例中,当所述ST波段基线特征为所述S波序列除余值时,所述特征获取模块820,进一步用于计算所述S波序列与所述SQ波深度的比值,得到S波序列除数;计算所述S波序列除数的余数,得到所述S波序列除余值。
在一个实施例中,当所述ST波段基线特征为所述T波序列除余值时,所述特征获取模块820,进一步用于计算所述T波序列与所述SQ波深度的比值,得到T波序列除数;计算所述T波序列除数的余数,得到所述T波序列除余值。
在一个实施例中,所述特征输入模块830,进一步用于在所述信号波形中,提取RR间期序列;根据所述RR间期序列,获取RR间期统计特征;所述RR间期统计特征包括所述RR间期序列的数量、均值、最大值、最小值、中位数、方差、偏度、峰度以及宽度中的至少一种;将所述RR间期统计特征、所述ST波段斜率特征和所述ST波段基线特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期统计特征、所述ST波段斜率特征和所述ST波段基线特征进行分类的结果。
在一个实施例中,当所述RR间期统计特征为所述RR间期序列的宽度时,所述特征输入模块830,进一步用于计算所述RR间期序列的最小值,得到RR间期序列最小值,以及,计算所述RR间期序列的最大值,得到RR间期序列最大值;计算所述RR间期序列最大值与所述RR间期序列最小值的差值,作为所述RR间期序列的宽度。
在一个实施例中,当所述形态特征为P波变化特征时,所述特征获取模块820,进一步用于获取P波序列;计算所述P波序列的最小值,得到P波序列最小值,以及,计算所述P波序列的最大值,得到P波序列最大值;计算所述P波序列最大值与所述P波序列最小值的差值,得到P波序列差值;计算所述P波序列差值的最大值,得到P波变化特征分子,以及,计算所述P波序列的最大值,得到P波变化特征分母;计算所述P波变化特征分子与所述P波变化特征分母的比值,作为所述P波变化特征。
在一个实施例中,所述获取模块820,进一步用于在所述信号波形中,提取QRS间期信号序列;根据所述QRS间期信号序列,获取QRS间期统计特征;所述QRS间期统计特征包括所述QRS间期信号序列的均值、标准差、方差、核密度、峰度以及偏度中的至少一种;将所述QRS间期统计特征和所述P波变化特征输入至所述分类器;所述分类结果为所述分类器根据所述QRS间期统计特征和所述P波变化特征进行分类的结果。
在一个实施例中,所述特征获取模块820,进一步用于采用所述信号波形训练基于深度学习的节律特征检测器;获取训练后的节律特征检测器的最优网络结构;获取所述最优网络结构的全连接层数据,作为所述深层特征。
在一个实施例中,所述节律特征检测器包括卷积神经网络检测器、循环神经网络检测器、卷积循环神经网络检测器、深度神经网络检测器中的至少一种。
在一个实施例中,所述分类器包括K-近邻(KNN)分类器、adboost分类器、xgboost分类器、逻辑回归分类器、支持向量机分类器中的至少一种。
在一个实施例中,所述分类模块840,进一步用于获取所述K-近邻(KNN)分类器、所述adboost分类器、所述xgboost分类器、所述逻辑回归分类器或所述支持向量机分类器的分类结果;将所述分类结果输入至多层感知器;获取所述多层感知器的多模型分类结果,得到所述心电信号的信号类型;所述多模型分类结果为所述多层感知器根据所述形态特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
上述提供的心电信号检测分类装置可用于执行上述任意实施例提供的心电信号检测分类方法,具备相应的功能和有益效果。
关于心电信号检测分类装置的具体限定可以参见上文中对于心电信号检测分类方法的限定,在此不再赘述。上述心电信号检测分类装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于电子设备中的处理器中,也可以以软件形式存储于电子设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
实施例三
图9是本申请实施例三提供的一种电子设备的结构示意图。如图中所示,该电子设备包括:处理器90、存储器91、具有触摸功能的显示屏92、输入装置93、输出装置94以及通信装置95。该电子设备中处理器90的数量可以是一个或者多个,图中以一个处理器90为例。该电子设备中存储器91的数量可以是一个或者多个,图中以一个存储器91为例。该电子设备的处理器90、存储器91、显示屏92、输入装置93、输出装置94以及通信装置95可以通过总线或者其他方式连接,图中以通过总线连接为例。实施例中,电子设备可以是电脑,手机,平板,投影仪或交互智能平板等。实施例中,以电子设备为交互智能平板为例,进行描述。
存储器91作为一种计算机可读存储介质,可用于存储软件程序、计算机可执行程序以及模块,如本申请任意实施例所述的心电信号检测分类方法对应的程序指令/模块。存储器91可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序;存储数据区可存储根据设备的使用所创建的数据等。此外,存储器91可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中,存储器91可进一步包括相对于处理器90远程设置的存储器,这些远程存储器可以通过网络连接至设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
显示屏92为具有触摸功能的显示屏92,其可以是电容屏、电磁屏或者红外屏。一般而言,显示屏92用于根据处理器90的指示显示数据,还用于接收作用于显示屏92的触摸操作,并将相应的信号发送至处理器90或其他装置。可选的,当显示屏92为红外屏时,其还包括红外触摸框,该红外触摸框设置在显示屏92的四周,其还可以用于接收红外信号,并将该红外信号发送至处理器90或者其他设备。
通信装置95,用于与其他设备建立通信连接,其可以是有线通信装置和/或无线通信装置。
输入装置93可用于接收输入的数字或者字符信息,以及产生与电子设备的用户设置以及功能控制有关的键信号输入,还可以是用于获取图像的摄像头以及获取音频数据的拾音设备。输出装置94可以包括扬声器等音频设备。需要说明的是,输入装置93和输出装置94的具体组成可以根据实际情况设定。
处理器90通过运行存储在存储器91中的软件程序、指令以及模块,从而执行设备的各种功能应用以及数据处理,即实现上述的心电信号检测分类方法。
具体的,实施例中,处理器90执行存储器91中存储的一个或多个程序时,具体实现如下操作:
在心电信号中,提取信号波形;
获取所述信号波形的形态特征和深层特征;所述形态特征包括TR波幅度差特征、PR波数量比值特征、ST波段特征、P波变化特征中的任意一种;所述深层特征包括深度特征和层次特征;
将所述形态特征和所述深层特征输入至分类器;
获取所述分类器输出的分类结果,得到所述心电信号的信号类型;所述分类结果为所述分类器根据所述形态特征和所述深层特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
在上述实施例的基础上,当所述形态特征为TR波幅度差特征时,一个或多个处理器90还实现如下操作:
获取R波信号序列,以及,获取T波信号序列;
计算所述T波信号序列与所述R波信号序列的比值,得到TR波序列比值;
计算所述TR波序列比值的中位数,得到TR波序列比值中位数;
确定所述TR波序列比值中位数和所述TR波序列比值,作为所述TR波幅度差特征。
在上述实施例的基础上,一个或多个处理器90还实现如下操作:
在所述信号波形中,提取RR间期差序列;
根据所述RR间期差序列,获取RR间期差第一先验特征;所述RR间期差第一先验特征包括所述RR间期差序列的原始计数、近似熵、样本熵中的至少一种;
将所述RR间期差第一先验特征和所述TR波幅度差特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第一先验特征和所述TR波幅度差特征进行分类的结果。
在上述实施例的基础上,当所述RR间期差第一先验特征为所述RR间期差序列的原始计数时,一个或多个处理器90还实现如下操作:
计算所述RR间期差序列的绝对值,得到RR间期差序列绝对值;
当所述RR间期差序列绝对值小于等于预设的常量定值时,计算所述RR间期差序列之和,得到RR间期差序列和值;
当所述RR间期差序列和值等于二时,确定所述RR间期差序列和值,作为所述RR间期差序列的原始计数。
在上述实施例的基础上,当所述形态特征为PR波数量比值特征时,一个或多个处理器90还实现如下操作:
获取R波序列数量,以及,获取P波序列数量;所述P波序列数量为P波序列的总和;所述R波序列数量为R波序列的总和;
计算所述P波序列数量与所述R波序列数量的比值,得到PR波数量比值;
根据所述PR波数量比值,确定所述PR波数量比值特征;所述PR波数量比值特征包括所述PR波数量比值、PR波数量比值中位数、PR波数量比值标准差、PR波数量比值均值以及PR波数量比值方差中的至少一种。
在上述实施例的基础上,当所述PR波数量比值特征为所述PR波数量比值方差时,一个或多个处理器90还实现如下操作:
计算所述PR波数量比值的均值,得到所述PR波数量比值均值,以及,计算所述PR波数量比值的标准差,得到所述PR波数量比值标准差;
计算所述PR波数量比值标准差与所述PR波数量比值均值的比值,得到所述PR波数量比值方差。
在上述实施例的基础上,一个或多个处理器90还实现如下操作:
在所述信号波形中,提取RR间期序列和/或RR间期差序列;
根据所述RR间期序列和/或所述RR间期差序列,获取RR间期差第二先验特征;
将所述RR间期差第二先验特征和所述PR波数量比值特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第二先验特征和所述PR波数量比值特征进行分类的结果。
在上述实施例的基础上,一个或多个处理器90还实现如下操作:
计算所述RR间期序列的核密度,得到RR间期序列核密度;
计算所述RR间期序列核密度的峰度值,得到RR间期密度峰度,以及,计算所述RR间期序列核密度的偏度值,得到RR间期密度偏度;
确定所述RR间期密度偏度、所述RR间期密度峰度以及所述RR间期序列核密度中的至少一种,作为所述RR间期差第二先验特征。
在上述实施例的基础上,一个或多个处理器90还实现如下操作:
计算所述RR间期差序列的核密度,得到RR间期差序列核密度;
计算所述RR间期差序列和密度的峰度值,得到RR间期差密度峰度,以及,计算所述RR间期差序列核密度的偏度值,得到RR间期差密度偏度;
确定所述RR间期差密度偏度、所述RR间期差密度峰度以及所述RR间期差序列核密度中的至少一种,作为所述RR间期差第二先验特征。
在上述实施例的基础上,当所述形态特征为ST波段特征时,一个或多个处理器90还实现如下操作:
获取T波序列,以及,获取S波序列;
根据所述S波序列和所述T波序列,获取ST波段斜率特征和/或ST波段基线特征;
确定所述ST波段斜率特征和/或所述ST波段基线特征,作为所述ST波段特征。
在上述实施例的基础上,一个或多个处理器90还实现如下操作:
获取所述S波序列的S波信号序列,以及,获取所述T波序列的T波信号序列;
计算所述T波信号序列与所述S波信号序列的差值,得到ST波段斜率分子,以及,计算所述T波序列与所述S波序列的差值,得到ST波段斜率分母;
计算所述ST波段斜率分子与所述ST波段斜率分母的比值,作为所述ST波段斜率;
根据所述ST波段斜率,获取所述ST波段斜率特征;所述ST波段斜率特征包括所述ST波段斜率、ST波段斜率中位数、ST波段斜率方差以及ST波段斜率长度比值中的至少一种。
在上述实施例的基础上,当所述ST波段斜率特征为所述ST波段斜率方差时,一个或多个处理器90还实现如下操作:
计算所述ST波段斜率的均值,得到ST波段斜率均值,以及,计算所述ST波段斜率的标准差,得到所述ST波段斜率标准差;
计算所述ST波段斜率标准差与所述ST波段斜率均值的比值,得到所述ST波段斜率方差。
在上述实施例的基础上,当所述ST波段斜率特征为所述ST波段斜率长度比值时,一个或多个处理器90还实现如下操作:
获取QRS波段序列;所述QRS波段序列为所述S波序列与所述Q波序列的差值;
计算所述QRS波段序列的长度,得到ST波段斜率长度分母;
当所述ST波段斜率小于零时,计算所述ST波段斜率的长度,得到ST波段斜率长度分子;
计算所述ST波段斜率长度分子与所述ST波段斜率长度分母的比值,得到所述ST波段斜率长度比值。
在上述实施例的基础上,一个或多个处理器90还实现如下操作:
获取所述S波序列的S波信号序列;
当所述S波信号序列小于Q波信号序列时,确定SQ波向量;
提取所述SQ波向量中非零元素的位置索引,得到SQ波深度;
根据所述SQ波深度,获取所述ST波段基线特征;所述ST波段基线特征包括所述SQ波深度、S波序列深度、S波序列除余值、T波序列除余值以及Q波序列除余值中的至少一种。
在上述实施例的基础上,当所述ST波段基线特征为所述S波序列深度时,一个或多个处理器90还实现如下操作:
计算所述S波序列的长度,得到S波序列深度分母,以及,计算所述SQ波深度的长度,得到S波序列深度分子;
计算所述S波序列深度分子与所述S波序列深度分母的比值,得到所述S波序列深度。
在上述实施例的基础上,当所述ST波段基线特征为所述S波序列除余值时,一个或多个处理器90还实现如下操作:
计算所述S波序列与所述SQ波深度的比值,得到S波序列除数;
计算所述S波序列除数的余数,得到所述S波序列除余值。
在上述实施例的基础上,当所述ST波段基线特征为所述T波序列除余值时,一个或多个处理器90还实现如下操作:
计算所述T波序列与所述SQ波深度的比值,得到T波序列除数;
计算所述T波序列除数的余数,得到所述T波序列除余值。
在上述实施例的基础上,所述形态特征包括ST波段特征,一个或多个处理器90还实现如下操作:
在所述信号波形中,提取RR间期序列;
根据所述RR间期序列,获取RR间期统计特征;所述RR间期统计特征包括所述RR间期序列的数量、均值、最大值、最小值、中位数、方差、偏度、峰度以及宽度中的至少一种;
将所述RR间期统计特征、所述ST波段斜率特征和所述ST波段基线特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期统计特征、所述ST波段斜率特征和所述ST波段基线特征进行分类的结果。
在上述实施例的基础上,当所述RR间期统计特征为所述RR间期序列的宽度时,一个或多个处理器90还实现如下操作:
计算所述RR间期序列的最小值,得到RR间期序列最小值,以及,计算所述RR间期序列的最大值,得到RR间期序列最大值;
计算所述RR间期序列最大值与所述RR间期序列最小值的差值,作为所述RR间期序列的宽度。
在上述实施例的基础上,当所述形态特征为P波变化特征时,一个或多个处理器90还实现如下操作:
获取P波序列;
计算所述P波序列的最小值,得到P波序列最小值,以及,计算所述P波序列的最大值,得到P波序列最大值;
计算所述P波序列最大值与所述P波序列最小值的差值,得到P波序列差值;
计算所述P波序列差值的最大值,得到P波变化特征分子,以及,计算所述P波序列的最大值,得到P波变化特征分母;
计算所述P波变化特征分子与所述P波变化特征分母的比值,作为所述P波变化特征。
在上述实施例的基础上,一个或多个处理器90还实现如下操作:
在所述信号波形中,提取QRS间期信号序列;
根据所述QRS间期信号序列,获取QRS间期统计特征;所述QRS间期统计特征包括所述QRS间期信号序列的均值、标准差、方差、核密度、峰度以及偏度中的至少一种;
将所述QRS间期统计特征和所述P波变化特征输入至所述分类器;所述分类结果为所述分类器根据所述QRS间期统计特征和所述P波变化特征进行分类的结果。
在上述实施例的基础上,一个或多个处理器90还实现如下操作:
采用所述信号波形训练基于深度学习的节律特征检测器;
获取训练后的节律特征检测器的最优网络结构;
获取所述最优网络结构的全连接层数据,作为所述深层特征。
在上述实施例的基础上,所述节律特征检测器包括卷积神经网络检测器、循环神经网络检测器、卷积循环神经网络检测器、深度神经网络检测器中的至少一种。
在上述实施例的基础上,所述分类器包括K-近邻(KNN)分类器、adboost分类器、xgboost分类器、逻辑回归分类器、支持向量机分类器中的至少一种。
在上述实施例的基础上,一个或多个处理器90还实现如下操作:
获取所述K-近邻(KNN)分类器、所述adboost分类器、所述xgboost分类器、所述逻辑回归分类器或所述支持向量机分类器的分类结果;
将所述分类结果输入至多层感知器;
获取所述多层感知器的多模型分类结果,得到所述心电信号的信号类型;所述多模型分类结果为所述多层感知器根据所述形态特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
实施例四
本申请实施例四还提供一种包含计算机可执行指令的存储介质,所述计算机可执行指令在由计算机处理器执行时用于执行一种心电信号检测分类方法,包括:
在心电信号中,提取信号波形;
获取所述信号波形的形态特征和深层特征;所述形态特征包括TR波幅度差特征、PR波数量比值特征、ST波段特征、P波变化特征中的任意一种;所述深层特征包括深度特征和层次特征;
将所述形态特征和所述深层特征输入至分类器;
获取所述分类器输出的分类结果,得到所述心电信号的信号类型;所述分类结果为所述分类器根据所述形态特征和所述深层特征进行分类的结果;所述信号类型包括房颤节律、非房颤异常节律、正常窦性心律和噪声。
当然,本申请实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上所述的心电信号检测分类方法操作,还可以执行本申请任意实施例所提供的心电信号检测分类方法中的相关操作,且具备相应的功能和有益效果。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (22)
1.一种心电信号检测分类方法,其特征在于,包括:
在心电信号中,提取信号波形;
获取所述信号波形的形态特征和深层特征;所述形态特征包括TR波幅度差特征;所述深层特征包括深度特征和层次特征;
将所述形态特征和所述深层特征输入至分类器;
获取所述分类器输出的分类结果,得到所述心电信号的信号类型;所述分类结果为所述分类器根据所述形态特征和所述深层特征进行分类的结果;
所述获取所述信号波形的形态特征,包括:
获取R波信号序列,以及,获取T波信号序列;
计算所述T波信号序列与所述R波信号序列的比值,得到TR波序列比值;
计算所述TR波序列比值的中位数,得到TR波序列比值中位数;
确定所述TR波序列比值中位数和所述TR波序列比值,作为所述TR波幅度差特征;
当所述形态特征包括TR波幅度差特征时,将所述形态特征输入至分类器,包括:
在所述信号波形中,提取RR间期差序列;
根据所述RR间期差序列,获取所述TR波幅度差特征相关联的RR间期差第一先验特征;所述RR间期差第一先验特征包括近似熵、样本熵中的至少一种;
将所述RR间期差第一先验特征和所述TR波幅度差特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第一先验特征和所述TR波幅度差特征进行分类的结果。
2.根据权利要求1所述的方法,其特征在于,当所述形态特征还包括PR波数量比值特征,所述获取所述信号波形的形态特征,还包括:
获取R波序列数量,以及,获取P波序列数量;所述P波序列数量为P波序列的总和;所述R波序列数量为R波序列的总和;
计算所述P波序列数量与所述R波序列数量的比值,得到PR波数量比值;
根据所述PR波数量比值,确定所述PR波数量比值特征;所述PR波数量比值特征包括所述PR波数量比值、PR波数量比值中位数、PR波数量比值标准差、PR波数量比值均值以及PR波数量比值方差中的至少一种。
3.根据权利要求2所述的方法,其特征在于,当所述PR波数量比值特征为所述PR波数量比值方差时,所述根据所述PR波数量比值,确定所述PR波数量比值特征,包括:
计算所述PR波数量比值的均值,得到所述PR波数量比值均值,以及,计算所述PR波数量比值的标准差,得到所述PR波数量比值标准差;
计算所述PR波数量比值标准差与所述PR波数量比值均值的比值,得到所述PR波数量比值方差。
4.根据权利要求2所述的方法,其特征在于,当所述形态特征还包括PR波数量比值特征时,将所述形态特征输入至分类器,还包括:
在所述信号波形中,提取RR间期序列和/或RR间期差序列;
根据所述RR间期序列和/或所述RR间期差序列,获取PR波数量比值特征相关联的RR间期差第二先验特征;
将所述RR间期差第二先验特征和所述PR波数量比值特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第二先验特征和所述PR波数量比值特征进行分类的结果。
5.根据权利要求4所述的方法,其特征在于,所述根据所述RR间期序列,获取RR间期差第二先验特征,包括:
计算所述RR间期序列的核密度,得到RR间期序列核密度;
计算所述RR间期序列核密度的峰度值,得到RR间期密度峰度,以及,计算所述RR间期序列核密度的偏度值,得到RR间期密度偏度;
确定所述RR间期密度偏度、所述RR间期密度峰度以及所述RR间期序列核密度中的至少一种,作为所述RR间期差第二先验特征。
6.根据权利要求4所述的方法,其特征在于,所述根据所述RR间期序列和/或所述RR间期差序列,获取RR间期差第二先验特征,包括:
计算所述RR间期差序列的核密度,得到RR间期差序列核密度;
计算所述RR间期差序列和密度的峰度值,得到RR间期差密度峰度,以及,计算所述RR间期差序列核密度的偏度值,得到RR间期差密度偏度;
确定所述RR间期差密度偏度、所述RR间期差密度峰度以及所述RR间期差序列核密度中的至少一种,作为所述RR间期差第二先验特征。
7.根据权利要求1所述的方法,其特征在于,当所述形态特征还包括ST波段特征时,所述获取所述信号波形的形态特征,还包括:
获取T波序列,以及,获取S波序列;
根据所述S波序列和所述T波序列,获取ST波段斜率特征和/或ST波段基线特征;
确定所述ST波段斜率特征和/或所述ST波段基线特征,作为所述ST波段特征。
8.根据权利要求7所述的方法,其特征在于,所述根据所述S波序列和所述T波序列,获取ST波段斜率特征,包括:
获取所述S波序列的S波信号序列,以及,获取所述T波序列的T波信号序列;
计算所述T波信号序列与所述S波信号序列的差值,得到ST波段斜率分子,以及,计算所述T波序列与所述S波序列的差值,得到ST波段斜率分母;
计算所述ST波段斜率分子与所述ST波段斜率分母的比值,作为所述ST波段斜率;
根据所述ST波段斜率,获取所述ST波段斜率特征;所述ST波段斜率特征包括所述ST波段斜率、ST波段斜率中位数、ST波段斜率方差以及ST波段斜率长度比值中的至少一种。
9.根据权利要求8所述的方法,其特征在于,当所述ST波段斜率特征为所述ST波段斜率方差时,所述根据所述ST波段斜率,获取所述ST波段斜率特征,包括:
计算所述ST波段斜率的均值,得到ST波段斜率均值,以及,计算所述ST波段斜率的标准差,得到所述ST波段斜率标准差;
计算所述ST波段斜率标准差与所述ST波段斜率均值的比值,得到所述ST波段斜率方差。
10.根据权利要求8所述的方法,其特征在于,当所述ST波段斜率特征为所述ST波段斜率长度比值时,所述根据所述ST波段斜率,获取所述ST波段斜率特征,包括:
获取QRS波段序列;所述QRS波段序列为所述S波序列与所述Q波序列的差值;
计算所述QRS波段序列的长度,得到ST波段斜率长度分母;
当所述ST波段斜率小于零时,计算所述ST波段斜率的长度,得到ST波段斜率长度分子;
计算所述ST波段斜率长度分子与所述ST波段斜率长度分母的比值,得到所述ST波段斜率长度比值。
11.根据权利要求7所述的方法,其特征在于,当所述ST波段基线特征为S波序列深度时,所述确定所述ST波段基线特征,包括:
计算所述S波序列的长度,得到S波序列深度分母,以及,计算所述SQ波深度的长度,得到S波序列深度分子;
计算所述S波序列深度分子与所述S波序列深度分母的比值,得到所述S波序列深度。
12.根据权利要求7所述的方法,其特征在于,当所述形态特征还包括ST波段特征时,将所述形态特征输入至分类器,还包括:
在所述信号波形中,提取RR间期序列;
根据所述RR间期序列,获取所述ST波段特征相关联的RR间期统计特征;所述RR间期统计特征包括所述RR间期序列的数量、均值、最大值、最小值、中位数、方差、偏度、峰度以及宽度中的至少一种;
将所述RR间期统计特征、所述ST波段斜率特征和所述ST波段基线特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期统计特征、所述ST波段斜率特征和所述ST波段基线特征进行分类的结果。
13.根据权利要求12所述的方法,其特征在于,当所述RR间期统计特征为所述RR间期序列的宽度时,所述根据所述RR间期序列,获取RR间期统计特征,包括:
计算所述RR间期序列的最小值,得到RR间期序列最小值,以及,计算所述RR间期序列的最大值,得到RR间期序列最大值;
计算所述RR间期序列最大值与所述RR间期序列最小值的差值,作为所述RR间期序列的宽度。
14.根据权利要求1所述的方法,其特征在于,当所述形态特征还包括P波变化特征时,所述获取所述信号波形的形态特征,还包括:
获取P波序列;
计算所述P波序列的最小值,得到P波序列最小值,以及,计算所述P波序列的最大值,得到P波序列最大值;
计算所述P波序列最大值与所述P波序列最小值的差值,得到P波序列差值;
计算所述P波序列差值的最大值,得到P波变化特征分子,以及,计算所述P波序列的最大值,得到P波变化特征分母;
计算所述P波变化特征分子与所述P波变化特征分母的比值,作为所述P波变化特征。
15.根据权利要求14所述的方法,其特征在于,所述将所述形态特征输入至分类器,包括:
在所述信号波形中,提取QRS间期信号序列;
根据所述QRS间期信号序列,获取QRS间期统计特征;所述QRS间期统计特征包括所述QRS间期信号序列的均值、标准差、方差、核密度、峰度以及偏度中的至少一种;
将所述QRS间期统计特征和所述P波变化特征输入至所述分类器;所述分类结果为所述分类器根据所述QRS间期统计特征和所述P波变化特征进行分类的结果。
16.根据权利要求1所述的方法,其特征在于,所述获取所述信号波形的深层特征,包括:
采用所述信号波形训练基于深度学习的节律特征检测器;
获取训练后的节律特征检测器的最优网络结构;
获取所述最优网络结构的全连接层数据,作为所述深层特征。
17.根据权利要求16所述的方法,其特征在于,所述节律特征检测器包括卷积神经网络检测器、循环神经网络检测器、卷积循环神经网络检测器、深度神经网络检测器中的至少一种。
18.根据权利要求1所述的方法,其特征在于,所述分类器包括K-近邻(KNN)分类器、adboost分类器、xgboost分类器、逻辑回归分类器、支持向量机分类器中的至少一种。
19.根据权利要求18所述的方法,其特征在于,所述获取所述分类器输出的分类结果,得到所述心电信号的信号类型,包括:
获取所述K-近邻(KNN)分类器、所述adboost分类器、所述xgboost分类器、所述逻辑回归分类器或所述支持向量机分类器的分类结果;
将所述分类结果输入至多层感知器;
获取所述多层感知器的多模型分类结果,得到所述心电信号的信号类型;所述多模型分类结果为所述多层感知器根据所述形态特征进行分类的结果。
20.一种心电信号检测分类装置,其特征在于,所述装置包括:
波形提取模块,用于在心电信号中,提取信号波形;
特征获取模块,用于获取所述信号波形的形态特征和深层特征;所述形态特征包括TR波幅度差特征;所述深层特征包括深度特征和层次特征;
特征输入模块,用于将所述形态特征和所述深层特征输入至分类器;
分类模块,用于获取所述分类器输出的分类结果,得到所述心电信号的信号类型;所述分类结果为所述分类器根据所述形态特征和所述深层特征进行分类的结果;
所述特征获取模块,还用于获取R波信号序列,以及,获取T波信号序列;
计算所述T波信号序列与所述R波信号序列的比值,得到TR波序列比值;
计算所述TR波序列比值的中位数,得到TR波序列比值中位数;
确定所述TR波序列比值中位数和所述TR波序列比值,作为所述TR波幅度差特征;
所述特征输入模块,还用于当所述形态特征包括TR波幅度差特征时,在所述信号波形中,提取RR间期差序列;
根据所述RR间期差序列,获取所述TR波幅度差特征相关联的RR间期差第一先验特征;所述RR间期差第一先验特征包括近似熵、样本熵中的至少一种;
将所述RR间期差第一先验特征和所述TR波幅度差特征输入至所述分类器;所述分类结果为所述分类器根据所述RR间期差第一先验特征和所述TR波幅度差特征进行分类的结果。
21.一种电子设备,其特征在于,包括:存储器、具有一个或多个处理器;
所述存储器,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器执行如权利要求1-19中任一所述的心电信号检测分类方法。
22.一种包含计算机可执行指令的存储介质,其特征在于,所述计算机可执行指令在由计算机处理器执行时用于执行如权利要求1-19中任一所述的心电信号检测分类方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910567569.9A CN110226921B (zh) | 2019-06-27 | 2019-06-27 | 心电信号检测分类方法、装置、电子设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910567569.9A CN110226921B (zh) | 2019-06-27 | 2019-06-27 | 心电信号检测分类方法、装置、电子设备和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110226921A CN110226921A (zh) | 2019-09-13 |
CN110226921B true CN110226921B (zh) | 2022-07-29 |
Family
ID=67857413
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910567569.9A Active CN110226921B (zh) | 2019-06-27 | 2019-06-27 | 心电信号检测分类方法、装置、电子设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110226921B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110664395B (zh) * | 2019-09-29 | 2022-04-12 | 京东方科技集团股份有限公司 | 图像处理方法、图像处理装置以及存储介质 |
CN110751131B (zh) * | 2019-11-16 | 2022-09-02 | 李汭傧 | 一种心律失常检测装置 |
CN110840402B (zh) * | 2019-11-19 | 2021-02-26 | 山东大学 | 一种基于机器学习的房颤信号识别方法及系统 |
CN110960207A (zh) * | 2019-12-16 | 2020-04-07 | 成都天奥电子股份有限公司 | 一种基于树模型的房颤检测方法、装置、设备及存储介质 |
CN111000553B (zh) * | 2019-12-30 | 2022-09-27 | 山东省计算中心(国家超级计算济南中心) | 一种基于投票集成学习的心电数据智能分类方法 |
CN111772628B (zh) * | 2020-07-16 | 2021-10-08 | 华中科技大学 | 一种基于深度学习的心电信号房颤自动检测系统 |
CN111832537B (zh) * | 2020-07-27 | 2023-04-25 | 深圳竹信科技有限公司 | 异常心电信号识别方法及异常心电信号识别装置 |
CN111973176B (zh) * | 2020-08-03 | 2024-03-12 | 广州视源电子科技股份有限公司 | 心电信号的处理方法、装置及设备 |
CN112232253B (zh) * | 2020-10-26 | 2021-11-09 | 南通大学 | 一种粗糙合理粒度驱动的胎儿心电图异常检测方法 |
CN112971800B (zh) * | 2021-01-27 | 2023-01-03 | 东南大学 | 基于深度学习的心电早搏信号处理方法、模块和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106214123A (zh) * | 2016-07-20 | 2016-12-14 | 杨平 | 一种基于深度学习算法的心电图综合分类方法 |
CN109009073A (zh) * | 2018-07-20 | 2018-12-18 | 广州视源电子科技股份有限公司 | 房颤检测装置及存储介质 |
CN109522916A (zh) * | 2017-09-19 | 2019-03-26 | 塔塔咨询服务有限公司 | 识别单导心电图(ecg)信号中的节律的级联二元分类器 |
CN109770862A (zh) * | 2019-03-29 | 2019-05-21 | 广州视源电子科技股份有限公司 | 心电信号分类方法、装置、电子设备和存储介质 |
-
2019
- 2019-06-27 CN CN201910567569.9A patent/CN110226921B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106214123A (zh) * | 2016-07-20 | 2016-12-14 | 杨平 | 一种基于深度学习算法的心电图综合分类方法 |
CN109522916A (zh) * | 2017-09-19 | 2019-03-26 | 塔塔咨询服务有限公司 | 识别单导心电图(ecg)信号中的节律的级联二元分类器 |
CN109009073A (zh) * | 2018-07-20 | 2018-12-18 | 广州视源电子科技股份有限公司 | 房颤检测装置及存储介质 |
CN109770862A (zh) * | 2019-03-29 | 2019-05-21 | 广州视源电子科技股份有限公司 | 心电信号分类方法、装置、电子设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN110226921A (zh) | 2019-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110226921B (zh) | 心电信号检测分类方法、装置、电子设备和存储介质 | |
Li et al. | Patient-specific ECG classification by deeper CNN from generic to dedicated | |
Wang et al. | Arrhythmia classification algorithm based on multi-head self-attention mechanism | |
Yang et al. | Automatic recognition of arrhythmia based on principal component analysis network and linear support vector machine | |
Sellami et al. | A robust deep convolutional neural network with batch-weighted loss for heartbeat classification | |
Oh et al. | Automated diagnosis of arrhythmia using combination of CNN and LSTM techniques with variable length heart beats | |
Pandey et al. | Patient specific machine learning models for ECG signal classification | |
Limam et al. | Atrial fibrillation detection and ECG classification based on convolutional recurrent neural network | |
Lu et al. | Feature fusion for imbalanced ECG data analysis | |
Srivastva et al. | PlexNet: A fast and robust ECG biometric system for human recognition | |
Kaya et al. | Classification of premature ventricular contraction in ECG | |
Huanhuan et al. | Classification of electrocardiogram signals with deep belief networks | |
Bashar et al. | Novel density Poincaré plot based machine learning method to detect atrial fibrillation from premature atrial/ventricular contractions | |
CN109770862B (zh) | 心电信号分类方法、装置、电子设备和存储介质 | |
Saber et al. | Arrhythmia modern classification techniques: A review | |
Ye et al. | Combining general multi-class and specific two-class classifiers for improved customized ECG heartbeat classification | |
CN110522442B (zh) | 多导联心电异常检测装置、电子设备和存储介质 | |
Sarfraz et al. | Using independent component analysis to obtain feature space for reliable ECG Arrhythmia classification | |
Dhull et al. | ECG beat classifiers: a journey from ANN to DNN | |
CN109077720B (zh) | 信号处理方法、装置、设备和存储介质 | |
Gawande et al. | Heart diseases classification using convolutional neural network | |
Feyisa et al. | Lightweight multireceptive field CNN for 12-lead ECG signal classification | |
Xie et al. | Intelligent analysis of premature ventricular contraction based on features and random forest | |
CN110522443B (zh) | 基于心电信号的房室传导阻滞检测方法、装置和电子设备 | |
CN108537123A (zh) | 基于多特征提取的心电识别方法 |
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 |