CN110141216A - 一种心电信号qrs特征波的识别方法、训练方法及系统 - Google Patents

一种心电信号qrs特征波的识别方法、训练方法及系统 Download PDF

Info

Publication number
CN110141216A
CN110141216A CN201910455091.0A CN201910455091A CN110141216A CN 110141216 A CN110141216 A CN 110141216A CN 201910455091 A CN201910455091 A CN 201910455091A CN 110141216 A CN110141216 A CN 110141216A
Authority
CN
China
Prior art keywords
electrocardiosignal
neural network
network model
signal
class sample
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
CN201910455091.0A
Other languages
English (en)
Other versions
CN110141216B (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.)
Dongguan Jianda Information Technology Co ltd
Shenzhen International Graduate School of Tsinghua University
Original Assignee
Dongguan Jianda Information Technology Co Ltd
Shenzhen Graduate School Tsinghua University
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 Dongguan Jianda Information Technology Co Ltd, Shenzhen Graduate School Tsinghua University filed Critical Dongguan Jianda Information Technology Co Ltd
Priority to CN201910455091.0A priority Critical patent/CN110141216B/zh
Publication of CN110141216A publication Critical patent/CN110141216A/zh
Application granted granted Critical
Publication of CN110141216B publication Critical patent/CN110141216B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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]
    • 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • 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/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • 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)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Artificial Intelligence (AREA)
  • Signal Processing (AREA)
  • Cardiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (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特征波的识别系统、神经网络模型的训练装置、计算机设备以及计算机可读存储介质，涉及生物医学技术领域。其中，心电信号QRS特征波的识别方法包括：获取待识别的心电信号；将所述心电信号进行预处理，得到信号片段；基于残差网络构造神经网络模型并对所述神经网络模型进行训练；基于所述神经网络模型对所述信号片段进行识别，得到所述心电信号中的QRS特征波。本发明由于在模型训练过程中有效避免了梯度消失的问题，因此具有执行效率较高的特点，在整个算法流程中有效减少了特征工程的步骤。

Description

一种心电信号QRS特征波的识别方法、训练方法及系统
技术领域
本发明关于生物医学技术领域，特别是关于生物医学中信号的处理技术，具体的讲是一种心电信号QRS特征波的识别方法、神经网络模型的训练方法、心电信号QRS特征波的识别系统、神经网络模型的训练装置、计算机设备以及计算机可读存储介质。
背景技术
本部分旨在为权利要求书中陈述的本发明的实施方式提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
心电图反映了心脏的肌电活动，QRS特征波是心电信号中反映心室收缩最明显的复合波。QRS特征波具有特殊的形状，是心电信号中特征最为突出的波形之一，它为心跳的自动检测提供了最重要的信息，检测QRS特征波也成为了研究已久的重要课题。因此，提高自动心跳检测的质量对于从准确的心跳定位中获益良多的功能非常重要，准确获取心率和心率变异性具有重要的临床意义。
以卷积神经网络为代表的深度学习技术发展迅速，利用深度学习技术对医疗影像中的潜在疾病进行识别和睡眠分期等领域基于神经网络的深度学习算法都取得了比较大的成功，有一些甚至超过了人类专家医生的诊疗结果。所以基于神经网络的深度学习算法对生物医疗信号的处理具有很大的潜力。残差网络是在2015年被提出，曾经在多个国际知名分类比赛任务上夺冠，因为它“简单与实用”并存，所以在图像分类、目标检测中得到了广泛的应用。残差网络有效解决了在神经网络模型训练过程中的梯度消失问题，能够有效提高模型准确率。
现阶段在检测QRS特征波方面主要是利用数字信号处理的相关知识，结合了滤波、小波变换和动态阈值等方法。根据QRS特征波的形态特征，对信号做一阶差分得到其导数信息，之后对一阶差分后的信号平方，对平方后的一阶差分信号做滑动窗积分；然后通过设置动态阈值找到QRS特征波。此过程有很多特征工程的步骤。
因此，如何提供一种新的心电信号的识别方案，其能够简化特征工程的步骤，实现端到端的心电特征波检测是本领域亟待解决的技术难题。
发明内容
有鉴于此，本发明实施例提供了一种心电信号QRS特征波的识别方法、神经网络模型的训练方法、心电信号QRS特征波的识别系统、神经网络模型的训练装置、计算机设备以及计算机可读存储介质，基于残差网络构造了神经网络模型，通过对心电信号进行样本分类以及切片处理后对神经网络模型进行训练，之后应用训练后的神经网络模型对待识别的心电信号进行识别，对神经网络模型的输出结果进行聚类即可识别出心电信号中的QRS特征波，由于在模型训练过程中有效避免了梯度消失的问题，因此本发明具有执行效率较高的特点，在整个算法流程中有效减少了特征工程的步骤。
本发明的目的之一是，提供一种心电信号QRS特征波的识别方法，包括：
获取待识别的心电信号；
将所述心电信号进行预处理，得到信号片段；
基于残差网络构造神经网络模型并对所述神经网络模型进行训练；
基于所述神经网络模型对所述信号片段进行识别，得到所述心电信号中的QRS特征波。
本发明的目的之一是，提供一种神经网络模型的训练方法，包括：
获取心电信号；
将所述心电信号进行预处理，得到训练阶段的信号片段；
基于残差网络构造神经网络模型；
根据所述信号片段对所述神经网络模型进行训练。
本发明的目的之一是，提供一种心电信号QRS特征波的识别系统，包括：
心电信号获取装置，用于获取待识别的心电信号；
心电信号预处理装置，用于将所述心电信号进行预处理，得到信号片段；
网络模型训练装置，用于基于残差网络构造神经网络模型并对所述神经网络模型进行训练；
信号片段识别装置，用于基于神经网络模型对所述信号片段进行识别，得到所述心电信号中的QRS特征波。
本发明的目的之一是，提供一种神经网络模型的训练装置，所述训练装置包括：
心电信号获取模块，用于获取心电信号；
心电信号预处理模块，用于将所述心电信号进行预处理，得到训练阶段的信号片段；
网络模型构造模块，用于基于残差网络构造神经网络模型；
网络模型训练模块，用于根据所述信号片段对所述神经网络模型进行训练。
本发明的目的之一是，提供一种计算机设备，包括：适于实现各指令的处理器以及存储设备，所述存储设备存储有多条指令，所述指令适于由处理器加载并执行一种心电信号QRS特征波的识别方法或如一种神经网络模型的训练方法。
本发明的目的之一是，提供一种计算机可读存储介质，存储有计算机程序，所述计算机程序用于执行一种心电信号QRS特征波的识别方法或一种神经网络模型的训练方法。
本发明的有益效果在于，本发明实施例提供了一种心电信号QRS特征波的识别方法、神经网络模型的训练方法、心电信号QRS特征波的识别系统、神经网络模型的训练装置、计算机设备以及计算机可读存储介质，基于残差网络构造了神经网络模型，通过对心电信号进行样本分类以及切片处理后对神经网络模型进行训练，之后应用训练后的神经网络模型对待识别的心电信号进行识别，对神经网络模型的输出结果进行聚类即可识别出心电信号中的QRS特征波，由于在模型训练过程中有效避免了梯度消失的问题，因此本发明具有执行效率较高的特点，在整个算法流程中有效减少了特征工程的步骤。
为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种心电信号QRS特征波的识别系统的结构示意图；
图2为本发明实施例提供的一种心电信号QRS特征波的识别系统中心电信号预处理装置的结构示意图；
图3为本发明实施例提供的一种心电信号QRS特征波的识别系统中滤波去燥处理装置的结构示意图；
图4为本发明实施例提供的一种心电信号QRS特征波的识别系统中归一化处理装置的结构示意图；
图5为本发明实施例提供的一种心电信号QRS特征波的识别系统中上下界确定单元的结构示意图；
图6为本发明实施例提供的一种心电信号QRS特征波的识别系统中切片处理装置的结构示意图；
图7为本发明实施例提供的一种心电信号QRS特征波的识别系统中信号片段识别装置的结构示意图；
图8为本发明实施例提供的一种神经网络模型的训练装置的结构示意图；
图9为本发明实施例提供的一种神经网络模型的训练装置中心电信号预处理模块的实施方式一的结构示意图；
图10为本发明实施例提供的一种神经网络模型的训练装置中样本分类模块的结构示意图；
图11为本发明实施例提供的一种神经网络模型的训练装置中心电信号预处理模块的实施方式二的结构示意图；
图12为本发明实施例提供的一种神经网络模型的训练装置中样本平衡处理模块的结构示意图；
图13为本发明实施例提供的一种神经网络模型的训练装置中样本切片模块的结构示意图；
图14为本发明实施例提供的一种神经网络模型的训练装置中网络模型构造模块的结构示意图；
图15为本发明实施例提供的一种心电信号QRS特征波的识别方法的流程示意图；
图16为图15中的步骤S102的具体流程示意图；
图17为图16中的步骤S201的具体流程示意图；
图18为图16中的步骤S202的具体流程示意图；
图19为图18中的步骤S401的具体流程示意图；
图20为图16中的步骤S203的具体流程示意图；
图21为图15中的步骤S104的具体流程示意图；
图22为本发明实施例提供的一种神经网络模型的训练方法的流程示意图；
图23为图22中的步骤S12的实施方式一的具体流程示意图；
图24为图23中的步骤S23的具体流程示意图；
图25为图22中的步骤S12的实施方式二的具体流程示意图；
图26为图25中的步骤S45的具体流程示意图；
图27为图25中的步骤S46的具体流程示意图；
图28为图22中的步骤S13的具体流程示意图；
图29为本发明提供的具体实施例中原始的心电信号的波形示意图；
图30为本发明提供的具体实施例中利用零相移滤波器去噪后的心电信号的波形示意图；
图31为本发明提供的具体实施例中归一化后的心电信号的波形示意图；
图32是本发明提供的具体实施例中对心电信号采样点进行样本分类的示意图。
图33是本发明提供的具体实施例中基于残差网络构建的神经网络模型示意图。
具体实施方式
下面将结合附图，对本发明实施例中的技术方案进行清楚、完整地描述参考在附图中示出并在以下描述中详述的非限制性示例实施例，更加全面地说明本发明的示例实施例和它们的多种特征及有利细节。应注意的是，图中示出的特征不是必须按照比例绘制。本发明省略了已知材料、组件和工艺技术的描述，从而不使本发明的示例实施例模糊。所给出的示例仅旨在有利于理解本发明示例实施例的实施，以及进一步使本领域技术人员能够实施示例实施例。因而，这些示例不应被理解为对本发明的实施例的范围的限制。
除非另外特别定义，本发明使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。本发明中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性，而只是用来区分不同的组成部分。此外，在本发明各个实施例中，相同或类似的参考标号表示相同或类似的构件。
在本发明中，首先基于残差网络构造识别心电QRS特征波的神经网络模型，其次构建神经网络模型的训练样本，以对该神经网络模型模型进行训练和测试。在构建神经网络模型的训练样本的过程中，根据是否属于QRS特征波对心电信号中的采样点进行正负样本分类，对归一化后的心电信号进行切片，信号片段的标签等于该信号片段中给定采样点的标签。训练好识别心电QRS特征波的神经网络模型后，首先获取待识别的心电信号，对心电信号进行滤波去噪、归一化、切片等预处理后，基于神经网络模型对心电信号进行识别，即可得到待识别的心电信号中的QRS特征波。该方案是基于残差网络实现，在模型训练过程中可以有效避免梯度消失的问题，具有执行效率较高的特点，在整个算法流程中有效减少了特征工程的步骤。
下面参考本发明的若干代表性实施方式，详细阐释本发明的原理和精神。图1为本发明实施例提供的一种心电信号QRS特征波的识别系统的结构示意图，请参见图1，所述心电信号QRS特征波的识别系统包括：
心电信号获取装置101，用于获取待识别的心电信号。原始的心电信号的波形如图29所示。
心电信号预处理装置102，用于将所述心电信号进行预处理，得到信号片段。
图2为心电信号预处理装置102的结构示意图，请参阅图2，心电信号预处理装置102包括：
滤波去燥处理装置201，用于将心电信号进行滤波去噪处理。在心电信号的采集过程中通常掺杂着各种噪声，比较常见的有由呼吸造成的基线漂移、工频干扰和肌电干扰等噪声。利用零相移滤波器去噪后的心电信号如图30所示，请参阅图29、图30，在本发明中，对心电信号进行滤波去噪处理可以去除心电采集过程中的基线漂移和工频干扰等噪声。
在本发明的一种实施方式中，滤波去燥处理装置201可利用双向零相移滤波器同时去除基线低频的基线漂移和各种高频噪声。去噪是处理心电信号过程中非常重要的一步，与普通的滤波器相比，零相移滤波器输出的信号不会造成信号采样点的偏移，这样能使得结果更准确；而且实现步骤简单，有利于提高算法的执行效率。
具体的，图3为滤波去燥处理装置201的结构示意图，请参阅图3，滤波去燥处理装置201包括：
第一滤波单元301，用于将所述心电信号进行滤波，得到第一信号；
第一逆转单元302，用于将所述第一信号进行逆转，得到第二信号；
第二滤波单元303，用于将所述第二信号进行滤波，得到第三信号；
第二逆转单元304，用于将所述第三信号进行逆转后输出，即得到滤波去噪处理之后的心电信号。
也即，首先将心电信号序列按顺序滤波，此时滤波器的输出是为Y(z)＝X(z)H(z)，然后将所得结果逆转，此时Y(z)＝X(z-1)H(z-1)，其次将反转后的信号通过滤波器，此时滤波器的输出是：Y(z)＝X(z-1)H(z-1)H(z)，最后再将滤波器的输出结果逆转输出，零相移滤波器的最终输出是Y(z)＝X(z)H(z)H(z-1)。零相移滤波在频域中最终输出可描述为Y(ejw)＝X(ejw)|H(ejw)|2，由此表达式可知输出和输入之间不存在附加相位。此滤波器有两个特点，首先是滤波后的心电信号波形更平滑，其次是滤波后的心电信号无位移。在常用的各类滤波器中，巴特沃斯滤波器具有在通带内幅频曲线幅度最平坦的特点，在实现零相移滤波器的过程中采用2阶巴特沃斯滤波器。
请参阅图2，心电信号预处理装置102还包括：
归一化处理装置202，用于将滤波去燥处理后的心电信号进行归一化处理。不同个体、设备采集的信号振幅不同，为了提高算法鲁棒性和泛化能力，需要对信号做归一化操作，同时在信号归一化的过程中需要避免极端异常值的影响。归一化后的心电信号如图31所示。
图4为本发明的一种实施方式中归一化处理装置202的结构示意图，请参阅图4，归一化处理装置202包括：
上下界确定单元401，用于确定归一化上界以及归一化下界。图5为本发明的一种实施方式中上下界确定单元401的结构示意图，请参阅图5，上下界确定单元401包括：
采样点排序单元501，用于将滤波去噪处理后的心电信号按照采样点的幅值从小到大进行排序，确定出上界索引以及下界索引；
上界确定单元502，用于根据所述上界索引确定归一化上界；
下界确定单元503，用于根据所述下界索引确定所述归一化下界。
上界索引为Upper_bound_index＝Num_samples×0.9974，则
归一化上界为Upper_bound＝ECG_sorted[Upper_bound_index]。
下界索引为Lower_bound_index＝Num_samples×0.0026，则
归一化下界值为Lower_bound＝ECG_sorted[Lower_bound_index]。
其中，Num_samples为此段心电信号中采样点的总个数，ECG_sorted为所有采样点的幅值从小到大排序后的心电信号，系数0.9974、0.0026可预先设定，符号[]代表数组下标。这种取值方式可以有效减少心电信号中极端异常值的影响。
请参阅图4，归一化处理装置202还包括：
归一化处理单元402，用于根据所述归一化上界以及归一化下界对所述滤波去噪处理后的心电信号进行归一化处理。在本发明的一种实施方式中，归一化处理可通过下述公式实现：
其中，ECG_nor和ECG_ori分别表示归一化后的心电信号和原始信号。
请参阅图2，心电信号预处理装置还包括：
切片处理装置203，用于将归一化处理后的心电信号进行切片处理，得到信号片段。信号片段的标签等于该信号片段中给定采样点的标签。
图6为本发明的一种实施方式中切片处理装置203的结构示意图，请参阅图6，切片处理装置203包括：
第一基准确定单元601，用于将归一化处理后的心电信号的采样点作为第一基准；
时间阈值获取单元602，用于获取预先设定的第一时间阈值以及第二时间阈值；
切片处理单元603，用于将所述第一基准前第一时间阈值时间和后第二时间阈值时间内的采样点作为一信号片段，所述信号片段开始的采样点和结束的采样点通过零填充的方式补充。
在具体的实施例中，对于给定的采样点，以该采样点为第一基准，取其前0.1s(即第一时间阈值)和后0.3s(即第二时间阈值)的采样点为一信号片段，开始和结束的采样点采用0填充的方式补充。
请参见图1，所述心电信号QRS特征波的识别系统还包括：
网络模型训练装置103，用于基于残差网络构造神经网络模型并对所述神经网络模型进行训练。神经网络模型的构造、训练过程详见后文，此处不再赘述。
信号片段识别装置104，用于基于神经网络模型对所述信号片段进行识别，得到所述心电信号中的QRS特征波。图7为本发明的一种实施方式中信号片段识别装置104的结构示意图，请参阅图7，信号片段识别装置104包括：
信号片段输出单元701，用于将所述信号片段输入至所述神经网络模型，得到正负标签序列；
标签序列聚类单元702，用于通过滑动窗口对所述正负标签序列进行聚类，识别出所述心电信号中的QRS特征波。
如上即为本发明实施例提供的一种心电信号QRS特征波的识别系统，基于残差网络构造了神经网络模型，并对神经网络模型进行训练，应用训练后的神经网络模型对待识别的心电信号进行识别，对神经网络模型的输出结果进行聚类即可识别出心电信号中的QRS特征波，由于在模型训练过程中有效避免了梯度消失的问题，因此本发明具有执行效率较高的特点，在整个算法流程中有效减少了特征工程的步骤。
图8为本发明实施例提供的一种神经网络模型的训练装置的结构示意图，请参见图8，所述神经网络模型的训练装置包括：
心电信号获取模块11，用于获取待心电信号。原始的心电信号的波形如图29所示。
心电信号预处理模块12，用于将所述心电信号进行预处理，得到训练阶段的信号片段。
图9为心电信号预处理模块12的实施方式一的结构示意图，请参阅图9，心电信号预处理模块12在实施方式一中包括：
滤波去燥处理模块21，用于将心电信号进行滤波去噪处理。在心电信号的采集过程中通常掺杂着各种噪声，比较常见的有由呼吸造成的基线漂移、工频干扰和肌电干扰等噪声。利用零相移滤波器去噪后的心电信号如图30所示，请参阅图29、图30，在本发明中，对心电信号进行滤波去噪处理可以去除心电采集过程中的基线漂移和工频干扰等噪声。
在本发明的一种实施方式中，滤波去燥处理模块21可利用双向零相移滤波器同时去除基线低频的基线漂移和各种高频噪声。去噪是处理心电信号过程中非常重要的一步，与普通的滤波器相比，零相移滤波器输出的信号不会造成信号采样点的偏移，这样能使得结果更准确；而且实现步骤简单，有利于提高算法的执行效率。
归一化处理模块22，用于将滤波去燥处理后的心电信号进行归一化处理。为减小不同个体信号之间的差异并增强算法的泛化能力，需要对滤波后的心电信号做归一化处理。归一化后的心电信号如图31所示。归一化处理的具体过程与归一化处理装置类似，此处不再赘述。
样本分类模块23，用于将归一化处理后的心电信号进行样本分类，得到正类样本以及负类样本。图10为本发明的一种实施方式中样本分类模块23的结构示意图，请参阅图10，样本分类模块23包括：
第三阈值获取模块31，用于获取预先设定的第三时间阈值；
第二基准确定模块32，用于将QRS复合波中R波波峰设为第二基准；
正类样本确定模块33，用于将所述第二基准前后的所述第三时间阈值时间内的采样点确定为正类样本；
负类样本确定模块34，用于将所述正类样本以外的采样点确定为负类样本。
在具体的实施例中，以QRS波中R波波峰索引为第二基准，前后各0.5s(即第三时间阈值)时间内信号的采样点为正类样本，其他的为负类样本。对心电信号采样点进行样本分类后的示意图如图32所示，图中标记黑点的是正类采样点，其余是负类采样点，“X”表示R波波峰。
请参阅图9，心电信号预处理模块12还包括：
样本切片模块24，用于将所述正类样本以及负类样本进行切片处理，得到信号片段。信号片段的标签等于该信号片段中给定采样点的标签。
图13为本发明的一种实施方式中样本切片模块24的结构示意图，请参阅图13，样本切片模块24包括：
第一基准确定模块51，用于将所述正类样本以及负类样本的采样点作为第一基准；
时间阈值获取模块52，用于获取预先设定的第一时间阈值以及第二时间阈值；
信号片段确定模块53，用于将所述第一基准前第一时间阈值时间和后第二时间阈值时间内的采样点作为一信号片段，所述信号片段开始的采样点和结束的采样点通过零填充的方式补充。
在具体的实施例中，对于给定的采样点，以该采样点为第一基准，取其前0.1s(即第一时间阈值)和后0.3s(即第二时间阈值)的采样点为一信号片段，开始和结束的采样点采用0填充的方式补充。
图11为本发明实施例中心电信号预处理模块的实施方式二的结构示意图，请参阅图11，在实施方式二中心电信号预处理模块还包括：
样本数量确定模块25，用于确定所述正类样本的数量以及负类样本的数量；
样本平衡处理模块26，用于根据所述正类样本的数量以及负类样本的数量进行样本平衡处理。
图12为本发明的一种实施方式中样本平衡处理模块26的结构示意图，请参阅图12，样本平衡处理模块26包括：
样本数量判断模块41，用于判断所述正类样本的数量是否小于所述负类样本的数量；
样本修正模块42，用于当所述判断模块判断为是时，在所述正类样本中加入低振幅噪声，得到修正后的正类样本。
在本发明的一种实施方式中，由于心电信号的信号片段中正负样本数量差距较大，正类样本数量较小，在正样本中加入低振幅噪声增加正样本数量。
请参见图8，神经网络模型的训练装置还包括：
网络模型构造模块13，用于基于残差网络构造神经网络模型。图14为本发明的一种实施方式中网络模型构造模块13的结构示意图，请参阅图14，网络模型构造模块13包括：
神经网络构建模块61，用于构建并行的第一卷积神经网络以及第二卷积神经网络；
第一卷积核设置模块62，用于设置所述第一卷积神经网络的卷积核，称为第一卷积核；
第二卷积核设置模块63，用于设置所述第二卷积神经网络的卷积核，称为第二卷积核，所述第一卷积核的尺寸大于所述第二卷机核的尺寸；
粗糙特征学习模块64，用于使用所述第一卷积神经网络学习所述心电信号的粗糙特征；
精细特征学习模块65，用于使用所述第二卷积神经网络学习所述心电信号的精细特征。
图33是本发明提供的具体实施例中基于残差网络构建的神经网络模型示意图，其中Input为卷积神经网络中的输入层，conv1d为卷积核，是一维的卷积层，BN即batchnormalization，是批归一化，PReLU为激活函数，max_pooling为最大池化层，Fullyconnected layer为全连接层，Softmax为此网络结构中的输出层。
请参见图8，神经网络模型的训练装置还包括：
网络模型训练模块14，用于根据所述信号片段对所述神经网络模型进行训练。
如上即为本发明实施例提供的一种神经网络模型的训练装置，基于残差网络构造了神经网络模型，并对神经网络模型进行训练，后续应用训练后的神经网络模型对待识别的心电信号进行识别，对神经网络模型的输出结果进行聚类即可识别出心电信号中的QRS特征波，由于在模型训练过程中有效避免了梯度消失的问题，因此本发明具有执行效率较高的特点，在整个算法流程中有效减少了特征工程的步骤。
此外，尽管在上文详细描述中提及了系统的若干单元模块，但是这种划分仅仅并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。同样，上文描述的一个单元的特征和功能也可以进一步划分为由多个单元来具体化。以上所使用的术语“模块”和“单元”，可以是实现预定功能的软件和/或硬件。尽管以下实施例所描述的模块较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
在介绍了本发明示例性实施方式的心电信号QRS特征波的识别系统、神经网络模型的训练装置之后，接下来，参考附图对本发明示例性实施方式的方法进行介绍。该方法的实施可以参见上述整体的实施，重复之处不再赘述。
图15为本发明实施例提供的一种心电信号QRS特征波的识别方法的流程意图，请参见图15，该方法包括：
S101：获取待识别的心电信号。原始的心电信号的波形如图29所示。
S102：将所述心电信号进行预处理，得到信号片段。
图16为步骤S102的流程示意图，请参阅图16，该步骤包括：
S201：将心电信号进行滤波去噪处理。在心电信号的采集过程中通常掺杂着各种噪声，比较常见的有由呼吸造成的基线漂移、工频干扰和肌电干扰等噪声。利用零相移滤波器去噪后的心电信号如图30所示，请参阅图29、图30，在本发明中，对心电信号进行滤波去噪处理可以去除心电采集过程中的基线漂移和工频干扰等噪声。
在本发明的一种实施方式中，该步骤可利用双向零相移滤波器同时去除基线低频的基线漂移和各种高频噪声。去噪是处理心电信号过程中非常重要的一步，与普通的滤波器相比，零相移滤波器输出的信号不会造成信号采样点的偏移，这样能使得结果更准确；而且实现步骤简单，有利于提高算法的执行效率。
具体的，图17为步骤S201的结构示意图，请参阅图17，该步骤包括：
S301：将所述心电信号进行滤波，得到第一信号；
S302：将所述第一信号进行逆转，得到第二信号；
S303：将所述第二信号进行滤波，得到第三信号；
S304：将所述第三信号进行逆转后输出，即得到滤波去噪处理之后的心电信号。
请参阅图15，S102还包括：
S202：将滤波去燥处理后的心电信号进行归一化处理。不同个体、设备采集的信号振幅不同，为了提高算法鲁棒性和泛化能力，需要对信号做归一化操作，同时在信号归一化的过程中需要避免极端异常值的影响。归一化后的心电信号如图31所示。
图18为S202的流程示意图，请参阅图18，该步骤包括：
S401，用于确定归一化上界以及归一化下界。图19为S401的流程示意图，请参阅图19，该步骤包括：
S501：将滤波去噪处理后的心电信号按照采样点的幅值从小到大进行排序，确定出上界索引以及下界索引；
S502：根据所述上界索引确定归一化上界；
S503：根据所述下界索引确定所述归一化下界。
上界索引为Upper_bound_index＝Num_samples×0.9974，则
归一化上界为Upper_bound＝ECG_sorted[Upper_bound_index]。
下界索引为Lower_bound_index＝Num_samples×0.0026，则
归一化下界值为Lower_bound＝ECG_sorted[Lower_bound_index]。
其中，Num_samples为此段心电信号中采样点的总个数，ECG_sorted为所有采样点的幅值从小到大排序后的心电信号，系数0.9974、0.0026可预先设定，符号[]代表数组下标。这种取值方式可以有效减少心电信号中极端异常值的影响。
请参阅图16，S202还包括：
S402：根据所述归一化上界以及归一化下界对所述滤波去噪处理后的心电信号进行归一化处理。在本发明的一种实施方式中，归一化处理可通过下述公式实现：
其中，ECG_nor和ECG_ori分别表示归一化后的心电信号和原始信号。
请参阅图15，S102还包括：
S203：将归一化处理后的心电信号进行切片处理，得到信号片段。信号片段的标签等于该信号片段中给定采样点的标签。
图20为本发明的一种实施方式中S203的流程示意图，请参阅图20，该步骤包括：
S601：将归一化处理后的心电信号的采样点作为第一基准；
S602：获取预先设定的第一时间阈值以及第二时间阈值；
S603：将所述第一基准前第一时间阈值时间和后第二时间阈值时间内的采样点作为一信号片段，所述信号片段开始的采样点和结束的采样点通过零填充的方式补充。
在具体的实施例中，对于给定的采样点，以该采样点为第一基准，取其前0.1s(即第一时间阈值)和后0.3s(即第二时间阈值)的采样点为一信号片段，开始和结束的采样点采用0填充的方式补充。
请参见图15，所述心电信号QRS特征波的识别方法还包括：
S103：基于残差网络构造神经网络模型并对所述神经网络模型进行训练。神经网络模型的构造、训练过程详见后文，此处不再赘述。
S104：基于神经网络模型对所述信号片段进行识别，得到所述心电信号中的QRS特征波。图21为本发明的一种实施方式中S104的流程示意图，请参阅图21，步骤S104包括：
S701：将所述信号片段输入至所述神经网络模型，得到正负标签序列；
S702：通过滑动窗口对所述正负标签序列进行聚类，识别出所述心电信号中的QRS特征波。
如上即为本发明实施例提供的一种心电信号QRS特征波的识别方法，基于残差网络构造了神经网络模型，并对神经网络模型进行训练，应用训练后的神经网络模型对待识别的心电信号进行识别，对神经网络模型的输出结果进行聚类即可识别出心电信号中的QRS特征波，由于在模型训练过程中有效避免了梯度消失的问题，因此本发明具有执行效率较高的特点，在整个算法流程中有效减少了特征工程的步骤。
图22为本发明实施例提供的一种神经网络模型的训练方法的流程示意图，请参见图22，所述神经网络模型的训练方法包括：
S11：获取待心电信号。原始的心电信号的波形如图29所示。
S12：将所述心电信号进行预处理，得到训练阶段的信号片段。
图23为S12的实施方式一的流程示意图，请参阅图23，S12在实施方式一中包括：
S21：将心电信号进行滤波去噪处理。在心电信号的采集过程中通常掺杂着各种噪声，比较常见的有由呼吸造成的基线漂移、工频干扰和肌电干扰等噪声。利用零相移滤波器去噪后的心电信号如图30所示，请参阅图29、图30，在本发明中，对心电信号进行滤波去噪处理可以去除心电采集过程中的基线漂移和工频干扰等噪声。
在本发明的一种实施方式中，该步骤可利用双向零相移滤波器同时去除基线低频的基线漂移和各种高频噪声。去噪是处理心电信号过程中非常重要的一步，与普通的滤波器相比，零相移滤波器输出的信号不会造成信号采样点的偏移，这样能使得结果更准确；而且实现步骤简单，有利于提高算法的执行效率。
S22：将滤波去燥处理后的心电信号进行归一化处理。为减小不同个体信号之间的差异并增强算法的泛化能力，需要对滤波后的心电信号做归一化处理。归一化后的心电信号如图31所示。归一化处理的具体过程与归一化处理装置类似，此处不再赘述。
S23：将归一化处理后的心电信号进行样本分类，得到正类样本以及负类样本。图24为本发明的一种实施方式中S23的流程示意图，请参阅图24，该步骤包括：
S31：获取预先设定的第三时间阈值；
S32：将QRS复合波中R波波峰设为第二基准；
S33：将所述第二基准前后的所述第三时间阈值时间内的采样点确定为正类样本；
S34：将所述正类样本以外的采样点确定为负类样本。
在具体的实施例中，以QRS波中R波波峰索引为第二基准，前后各0.5s(即第三时间阈值)时间内信号的采样点为正类样本，其他的为负类样本。对心电信号采样点进行样本分类后的示意图如图32所示，图中标记黑点的是正类采样点，其余是负类采样点，“X”表示R波波峰。
请参阅图23，步骤S12还包括：
S24：将所述正类样本以及负类样本进行切片处理，得到信号片段。信号片段的标签等于该信号片段中给定采样点的标签。
图25为S12的实施方式二的流程示意图，请参阅图25，在实施方式二中该步骤包括：
S41：将心电信号进行滤波去噪处理。
S42：将滤波去燥处理后的心电信号进行归一化处理。为减小不同个体信号之间的差异并增强算法的泛化能力，需要对滤波后的心电信号做归一化处理。归一化后的心电信号如图31所示。归一化处理的具体过程与归一化处理装置类似，此处不再赘述。
S43：将归一化处理后的心电信号进行样本分类，得到正类样本以及负类样本；
S44：确定所述正类样本的数量以及负类样本的数量；
S45：根据所述正类样本的数量以及负类样本的数量进行样本平衡处理。
S46：将所述正类样本以及负类样本进行切片处理，得到信号片段。信号片段的标签等于该信号片段中给定采样点的标签。
图26为图25中的步骤S45的流程示意图，请参阅图27，步骤S45包括：
S51：判断所述正类样本的数量是否小于所述负类样本的数量；
S52：当所述判断模块判断为是时，在所述正类样本中加入低振幅噪声，得到修正后的正类样本。
在本发明的一种实施方式中，由于心电信号的信号片段中正负样本数量差距较大，正类样本数量较小，在正样本中加入低振幅噪声增加正样本数量。
图27为图25中的步骤S46的流程示意图，请参阅图27，步骤S46包括：
S61：将所述正类样本以及负类样本的采样点作为第一基准；
S62：获取预先设定的第一时间阈值以及第二时间阈值；
S63：将所述第一基准前第一时间阈值时间和后第二时间阈值时间内的采样点作为一信号片段，所述信号片段开始的采样点和结束的采样点通过零填充的方式补充。
在具体的实施例中，对于给定的采样点，以该采样点为第一基准，取其前0.1s(即第一时间阈值)和后0.3s(即第二时间阈值)的采样点为一信号片段，开始和结束的采样点采用0填充的方式补充。请参见图22，神经网络模型的训练方法还包括：
S13：基于残差网络构造神经网络模型。图28为本发明的一种实施方式中S13的流程示意图，请参阅图28，步骤S13包括：
S71：构建并行的第一卷积神经网络以及第二卷积神经网络；
S72：设置所述第一卷积神经网络的卷积核，称为第一卷积核；
S73：设置所述第二卷积神经网络的卷积核，称为第二卷积核，所述第一卷积核的尺寸大于所述第二卷机核的尺寸；
S74：使用所述第一卷积神经网络学习所述心电信号的粗糙特征；
S75：使用所述第二卷积神经网络学习所述心电信号的精细特征。
图33是本发明提供的具体实施例中基于残差网络构建的神经网络模型示意图，其中Input为卷积神经网络中的输入层，conv1d为卷积核，是一维的卷积层，BN即batchnormalization，是批归一化，PReLU为激活函数，max_pooling为最大池化层，Fullyconnected layer为全连接层，Softmax为此网络结构中的输出层。
请参见图22，神经网络模型的训练方法还包括：
S14：根据所述信号片段对所述神经网络模型进行训练。
本发明还提供一种计算机设备，包括：适于实现各指令的处理器以及存储设备，所述存储设备存储有多条指令，所述指令适于由处理器加载并执行一种心电信号QRS特征波的识别方法或如一种神经网络模型的训练方法。
本发明还提供一种计算机可读存储介质，存储有计算机程序，所述计算机程序用于执行一种心电信号QRS特征波的识别方法或一种神经网络模型的训练方法。
综上所述，本发明实施例提供了一种心电信号QRS特征波的识别方法、神经网络模型的训练方法、心电信号QRS特征波的识别系统、神经网络模型的训练装置、计算机设备以及计算机可读存储介质，基于残差网络构造了神经网络模型，通过对心电信号进行样本分类以及切片处理后对神经网络模型进行训练，之后应用训练后的神经网络模型对待识别的心电信号进行识别，对神经网络模型的输出结果进行聚类即可识别出心电信号中的QRS特征波，由于在模型训练过程中有效避免了梯度消失的问题，因此本发明具有执行效率较高的特点，在整个算法流程中有效减少了特征工程的步骤。
在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现客户端和服务器以外，完全可以通过将方法步骤进行逻辑编程来使得客户端和服务器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种客户端和服务器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施方式或者实施方式的某些部分所述的方法。
本说明书中的各个实施方式均采用递进的方式描述，各个实施方式之间相同相似的部分互相参见即可，每个实施方式重点说明的都是与其他实施方式的不同之处。尤其，针对客户端和服务器的实施方式来说，均可以参照前述方法的实施方式的介绍对照解释。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施方式描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (30)

1.一种心电信号QRS特征波的识别方法，其特征在于，所述方法包括：
获取待识别的心电信号；
将所述心电信号进行预处理，得到信号片段；
基于残差网络构造神经网络模型并对所述神经网络模型进行训练；
基于所述神经网络模型对所述信号片段进行识别，得到所述心电信号中的QRS特征波。
2.根据权利要求1所述的方法，其特征在于，将所述心电信号进行预处理，得到信号片段包括：
将所述心电信号进行滤波去噪处理；
将滤波去燥处理后的心电信号进行归一化处理；
将归一化处理后的心电信号进行切片处理，得到信号片段。
3.根据权利要求2所述的方法，其特征在于，将所述心电信号进行滤波去噪处理包括：
将所述心电信号进行滤波，得到第一信号；
将所述第一信号进行逆转，得到第二信号；
将所述第二信号进行滤波，得到第三信号；
将所述第三信号进行逆转后输出，即得到滤波去噪处理之后的心电信号。
4.根据权利要求2所述的方法，其特征在于，对滤波去燥处理后的心电信号进行归一化处理包括：
确定归一化上界以及归一化下界；
根据所述归一化上界以及归一化下界对所述滤波去噪处理后的心电信号进行归一化处理。
5.根据权利要求4所述的方法，其特征在于，所述确定归一化上界以及归一化下界包括：
将滤波去噪处理后的心电信号按照采样点的幅值从小到大进行排序，确定出上界索引以及下界索引；
根据所述上界索引确定归一化上界；
根据所述下界索引确定所述归一化下界。
6.根据权利要求2所述的方法，其特征在于，对归一化处理后的心电信号进行切片处理包括：
将归一化处理后的心电信号的采样点作为第一基准；
获取预先设定的第一时间阈值以及第二时间阈值；
将所述第一基准前第一时间阈值时间和后第二时间阈值时间内的采样点作为一信号片段，所述信号片段开始的采样点和结束的采样点通过零填充的方式补充。
7.根据权利要求1所述的方法，其特征在于，基于神经网络模型对所述信号片段进行识别，得到所述心电信号对应的QRS特征波包括：
将所述信号片段输入至所述神经网络模型，得到正负标签序列；
通过滑动窗口对所述正负标签序列进行聚类，识别出所述心电信号中的QRS特征波。
8.一种神经网络模型的训练方法，其特征在于，所述方法包括：
获取心电信号；
将所述心电信号进行预处理，得到训练阶段的信号片段；
基于残差网络构造神经网络模型；
根据所述信号片段对所述神经网络模型进行训练。
9.根据权利要求8所述的方法，其特征在于，将所述心电信号进行预处理，得到训练阶段的信号片段包括：
将所述心电信号进行滤波去噪处理；
将滤波去燥处理后的心电信号进行归一化处理；
将归一化处理后的心电信号进行样本分类，得到正类样本以及负类样本；
将所述正类样本以及负类样本进行切片处理，得到训练阶段的信号片段。
10.根据权利要求9所述的方法，其特征在于，将归一化处理后的心电信号进行样本分类包括：
获取预先设定的第三时间阈值；
将QRS复合波中R波波峰设为第二基准；
将所述第二基准前后的所述第三时间阈值时间内的采样点确定为正类样本；
将所述正类样本以外的采样点确定为负类样本。
11.根据权利要求10所述的方法，其特征在于，将所述心电信号进行预处理还包括：
确定所述正类样本的数量以及负类样本的数量；
根据所述正类样本的数量以及负类样本的数量进行样本平衡处理。
12.根据权利要求11所述的方法，其特征在于，根据所述正类样本的数量以及负类样本的数量进行样本平衡处理包括：
判断所述正类样本的数量是否小于所述负类样本的数量；
当判断为是时，在所述正类样本中加入低振幅噪声，得到修正后的正类样本。
13.根据权利要求10或12所述的方法，其特征在于，将所述正类样本以及负类样本进行切片处理包括：
将所述正类样本以及负类样本的采样点作为第一基准；
获取预先设定的第一时间阈值以及第二时间阈值；
将所述第一基准前第一时间阈值时间和后第二时间阈值时间内的采样点作为一信号片段，所述信号片段开始的采样点和结束的采样点通过零填充的方式补充。
14.根据权利要求8所述的方法，其特征在于，所述基于残差网络构造神经网络模型包括：
构建并行的第一卷积神经网络以及第二卷积神经网络；
设置所述第一卷积神经网络的卷积核，称为第一卷积核；
设置所述第二卷积神经网络的卷积核，称为第二卷积核，所述第一卷积核的尺寸大于所述第二卷机核的尺寸；
使用所述第一卷积神经网络学习所述心电信号的粗糙特征；
使用所述第二卷积神经网络学习所述心电信号的精细特征。
15.一种心电信号QRS特征波的识别系统，其特征在于，所述系统包括：
心电信号获取装置，用于获取待识别的心电信号；
心电信号预处理装置，用于将所述心电信号进行预处理，得到信号片段；
网络模型训练装置，用于基于残差网络构造神经网络模型并对所述神经网络模型进行训练；
信号片段识别装置，用于基于神经网络模型对所述信号片段进行识别，得到所述心电信号中的QRS特征波。
16.根据权利要求15所述的系统，其特征在于，所述心电信号预处理装置包括：
滤波去燥处理装置，用于将心电信号进行滤波去噪处理；
归一化处理装置，用于将滤波去燥处理后的心电信号进行归一化处理；
切片处理装置，用于将归一化处理后的心电信号进行切片处理，得到信号片段。
17.根据权利要求16所述的系统，其特征在于，所述滤波去燥处理装置包括：
第一滤波单元，用于将所述心电信号进行滤波，得到第一信号；
第一逆转单元，用于将所述第一信号进行逆转，得到第二信号；
第二滤波单元，用于将所述第二信号进行滤波，得到第三信号；
第二逆转单元，用于将所述第三信号进行逆转后输出，即得到滤波去噪处理之后的心电信号。
18.根据权利要求16所述的系统，其特征在于，所述归一化处理装置包括：
上下界确定单元，用于确定归一化上界以及归一化下界；
归一化处理单元，用于根据所述归一化上界以及归一化下界对所述滤波去噪处理后的心电信号进行归一化处理。
19.根据权利要求18所述的系统，其特征在于，所述上下界确定单元包括：
采样点排序单元，用于将滤波去噪处理后的心电信号按照采样点的幅值从小到大进行排序，确定出上界索引以及下界索引；
上界确定单元，用于根据所述上界索引确定归一化上界；
下界确定单元，用于根据所述下界索引确定所述归一化下界。
20.根据权利要求16所述的系统，其特征在于，所述切片处理装置包括：
第一基准确定单元，用于将归一化处理后的心电信号的采样点作为第一基准；
时间阈值获取单元，用于获取预先设定的第一时间阈值以及第二时间阈值；
切片处理单元，用于将所述第一基准前第一时间阈值时间和后第二时间阈值时间内的采样点作为一信号片段，所述信号片段开始的采样点和结束的采样点通过零填充的方式补充。
21.根据权利要求15所述的系统，其特征在于，所述信号片段识别装置包括：
信号片段输出单元，用于将所述信号片段输入至所述神经网络模型，得到正负标签序列；
标签序列聚类单元，用于通过滑动窗口对所述正负标签序列进行聚类，识别出所述心电信号中的QRS特征波。
22.一种神经网络模型的训练装置，其特征在于，所述训练装置包括：
心电信号获取模块，用于获取心电信号；
心电信号预处理模块，用于将所述心电信号进行预处理，得到训练阶段的信号片段；
网络模型构造模块，用于基于残差网络构造神经网络模型；
网络模型训练模块，用于根据所述信号片段对所述神经网络模型进行训练。
23.根据权利要求22所述的装置，其特征在于，所述心电信号预处理模块包括：
滤波去燥处理模块，用于将心电信号进行滤波去噪处理；
归一化处理模块，用于将滤波去燥处理后的心电信号进行归一化处理；
样本分类模块，用于将归一化处理后的心电信号进行样本分类，得到正类样本以及负类样本；
样本切片模块，用于将所述正类样本以及负类样本进行切片处理，得到信号片段。
24.根据权利要求23所述的装置，其特征在于，所述样本分类模块包括：
第三阈值获取模块，用于获取预先设定的第三时间阈值；
第二基准确定模块，用于将QRS复合波中R波波峰设为第二基准；
正类样本确定模块，用于将所述第二基准前后的所述第三时间阈值时间内的采样点确定为正类样本；
负类样本确定模块，用于将所述正类样本以外的采样点确定为负类样本。
25.根据权利要求24所述的装置，其特征在于，所述心电信号预处理模块还包括：
样本数量确定模块，用于确定所述正类样本的数量以及负类样本的数量；
样本平衡处理模块，用于根据所述正类样本的数量以及负类样本的数量进行样本平衡处理。
26.根据权利要求25所述的装置，其特征在于，所述样本平衡处理模块包括：
样本数量判断模块，用于判断所述正类样本的数量是否小于所述负类样本的数量；
样本修正模块，用于当所述判断模块判断为是时，在所述正类样本中加入低振幅噪声，得到修正后的正类样本。
27.根据权利要求23所述的装置，其特征在于，所述样本切片模块包括：
第一基准确定模块，用于将所述正类样本以及负类样本的采样点作为第一基准；
时间阈值获取模块，用于获取预先设定的第一时间阈值以及第二时间阈值；
信号片段确定模块，用于将所述第一基准前第一时间阈值时间和后第二时间阈值时间内的采样点作为一信号片段，所述信号片段开始的采样点和结束的采样点通过零填充的方式补充。
28.根据权利要求22所述的装置，其特征在于，所述网络模型构造模块包括：
神经网络构建模块，用于构建并行的第一卷积神经网络以及第二卷积神经网络；
第一卷积核设置模块，用于设置所述第一卷积神经网络的卷积核，称为第一卷积核；
第二卷积核设置模块，用于设置所述第二卷积神经网络的卷积核，称为第二卷积核，所述第一卷积核的尺寸大于所述第二卷机核的尺寸；
粗糙特征学习模块，用于使用所述第一卷积神经网络学习所述心电信号的粗糙特征；
精细特征学习模块，用于使用所述第二卷积神经网络学习所述心电信号的精细特征。
29.一种计算机设备，其特征在于，包括：适于实现各指令的处理器以及存储设备，所述存储设备存储有多条指令，所述指令适于由处理器加载并执行如权利要求1至7任意一项所述的一种心电信号QRS特征波的识别方法或如权利要求8至14任意一项所述的一种神经网络模型的训练方法。
30.一种计算机可读存储介质，其特征在于，存储有计算机程序，所述计算机程序用于执行如权利要求1至7任意一项所述的一种心电信号QRS特征波的识别方法或如权利要求8至14任意一项所述的一种神经网络模型的训练方法。
CN201910455091.0A 2019-05-29 2019-05-29 一种心电信号qrs特征波的识别方法、训练方法及系统 Active CN110141216B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910455091.0A CN110141216B (zh) 2019-05-29 2019-05-29 一种心电信号qrs特征波的识别方法、训练方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910455091.0A CN110141216B (zh) 2019-05-29 2019-05-29 一种心电信号qrs特征波的识别方法、训练方法及系统

Publications (2)

Publication Number Publication Date
CN110141216A true CN110141216A (zh) 2019-08-20
CN110141216B CN110141216B (zh) 2022-09-30

Family

ID=67592042

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910455091.0A Active CN110141216B (zh) 2019-05-29 2019-05-29 一种心电信号qrs特征波的识别方法、训练方法及系统

Country Status (1)

Country Link
CN (1) CN110141216B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110403601A (zh) * 2019-08-27 2019-11-05 安徽心之声医疗科技有限公司 基于深度学习的心电信号qrs波群识别方法
CN111938640A (zh) * 2020-08-05 2020-11-17 深圳大学 一种心阻抗微分信号的c点定位方法、装置及存储介质
CN113081025A (zh) * 2021-03-10 2021-07-09 华中科技大学 一种基于肺音的新冠肺炎智能诊断系统和装置
CN113197582A (zh) * 2021-04-27 2021-08-03 浙江大学 一种高通用性的心电数据压缩感知系统、终端和存储介质
CN113647959A (zh) * 2021-07-27 2021-11-16 东软集团股份有限公司 一种心电波形信号的波形识别方法、装置及设备
CN113940686A (zh) * 2021-10-20 2022-01-18 武汉中旗生物医疗电子有限公司 一种ecg信号的工频滤除方法、装置、电子设备及存储介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101361650A (zh) * 2007-08-07 2009-02-11 深圳市理邦精密仪器有限公司 Iir滤波器的零相位实现方法及零相位iir滤波装置
CN105550653A (zh) * 2015-12-14 2016-05-04 广州视源电子科技股份有限公司 心电导联智能选择方法和系统
CN105769173A (zh) * 2016-02-29 2016-07-20 浙江铭众科技有限公司 一种具有心电信号去噪功能的心电监测系统
CN106667439A (zh) * 2016-12-30 2017-05-17 包磊 一种心电信号的降噪方法及装置
CN106725428A (zh) * 2016-12-19 2017-05-31 中国科学院深圳先进技术研究院 一种心电信号分类方法及装置
CN107657318A (zh) * 2017-11-13 2018-02-02 成都蓝景信息技术有限公司 一种基于深度学习模型的心电图分类方法
CN108229673A (zh) * 2016-12-27 2018-06-29 北京市商汤科技开发有限公司 卷积神经网络的处理方法、装置和电子设备
CN108416277A (zh) * 2018-02-11 2018-08-17 广州市碳码科技有限责任公司 一种心电监测方法、装置、终端及计算机可读存储介质
CN109480826A (zh) * 2018-12-14 2019-03-19 东软集团股份有限公司 一种心电信号处理方法、装置及设备
CN109508726A (zh) * 2017-09-15 2019-03-22 北京京东尚科信息技术有限公司 数据处理方法及其系统

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101361650A (zh) * 2007-08-07 2009-02-11 深圳市理邦精密仪器有限公司 Iir滤波器的零相位实现方法及零相位iir滤波装置
CN105550653A (zh) * 2015-12-14 2016-05-04 广州视源电子科技股份有限公司 心电导联智能选择方法和系统
CN105769173A (zh) * 2016-02-29 2016-07-20 浙江铭众科技有限公司 一种具有心电信号去噪功能的心电监测系统
CN106725428A (zh) * 2016-12-19 2017-05-31 中国科学院深圳先进技术研究院 一种心电信号分类方法及装置
CN108229673A (zh) * 2016-12-27 2018-06-29 北京市商汤科技开发有限公司 卷积神经网络的处理方法、装置和电子设备
CN106667439A (zh) * 2016-12-30 2017-05-17 包磊 一种心电信号的降噪方法及装置
CN109508726A (zh) * 2017-09-15 2019-03-22 北京京东尚科信息技术有限公司 数据处理方法及其系统
CN107657318A (zh) * 2017-11-13 2018-02-02 成都蓝景信息技术有限公司 一种基于深度学习模型的心电图分类方法
CN108416277A (zh) * 2018-02-11 2018-08-17 广州市碳码科技有限责任公司 一种心电监测方法、装置、终端及计算机可读存储介质
CN109480826A (zh) * 2018-12-14 2019-03-19 东软集团股份有限公司 一种心电信号处理方法、装置及设备

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
项延德: "基于卷积神经网络的心电信号检测和分类研究", 《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110403601A (zh) * 2019-08-27 2019-11-05 安徽心之声医疗科技有限公司 基于深度学习的心电信号qrs波群识别方法
CN111938640A (zh) * 2020-08-05 2020-11-17 深圳大学 一种心阻抗微分信号的c点定位方法、装置及存储介质
CN111938640B (zh) * 2020-08-05 2023-02-07 深圳大学 一种心阻抗微分信号的c点定位方法、装置及存储介质
CN113081025A (zh) * 2021-03-10 2021-07-09 华中科技大学 一种基于肺音的新冠肺炎智能诊断系统和装置
CN113081025B (zh) * 2021-03-10 2022-04-29 华中科技大学 一种基于肺音的新冠肺炎智能诊断系统和装置
CN113197582A (zh) * 2021-04-27 2021-08-03 浙江大学 一种高通用性的心电数据压缩感知系统、终端和存储介质
CN113647959A (zh) * 2021-07-27 2021-11-16 东软集团股份有限公司 一种心电波形信号的波形识别方法、装置及设备
CN113647959B (zh) * 2021-07-27 2024-01-12 东软集团股份有限公司 一种心电波形信号的波形识别方法、装置及设备
CN113940686A (zh) * 2021-10-20 2022-01-18 武汉中旗生物医疗电子有限公司 一种ecg信号的工频滤除方法、装置、电子设备及存储介质

Also Published As

Publication number Publication date
CN110141216B (zh) 2022-09-30

Similar Documents

Publication Publication Date Title
CN110141216A (zh) 一种心电信号qrs特征波的识别方法、训练方法及系统
CN108714026B (zh) 基于深度卷积神经网络和在线决策融合的细粒度心电信号分类方法
CN105899268B (zh) 基于gpu的并行心电信号分析方法
US20200237246A1 (en) Automatic recognition and classification method for electrocardiogram heartbeat based on artificial intelligence
CN109645989B (zh) 一种麻醉深度估计系统
CN108874149B (zh) 一种基于表面肌电信号连续估计人体关节角度的方法
CN110555468A (zh) 一种联合递归图和cnn的脑电信号识别方法及系统
CN110946566B (zh) 一种基于U-Net网络的心拍分类方法、装置、设备及存储介质
CN109567789A (zh) 心电图数据的分割处理方法、装置及可读存储介质
Wei et al. Motor imagery EEG signal classification based on deep transfer learning
CN111291727A (zh) 一种光体积变化描记图法信号质量检测方法和装置
CN110141215A (zh) 降噪自编码器的训练方法、心电信号的降噪方法及相关装置、设备
Asghar et al. Semi-skipping layered gated unit and efficient network: hybrid deep feature selection method for edge computing in EEG-based emotion classification
CN109859850B (zh) 基于变分模态分解和极限学习机的血糖预测方法
US10535001B2 (en) Reducing problem complexity when analyzing 3-D images
Gabardi et al. A multi-artifact EEG denoising by frequency-based deep learning
CN112790756A (zh) 任务态功能磁共振扫描数据去噪方法、装置、设备及介质
Golgowski et al. Anomaly detection in ECG using wavelet transformation
CN107837083A (zh) 基于最小二乘支持向量机的j波自动检测方法
Okreghe et al. A deep neural network-based spike sorting with improved channel selection and artefact removal
CN105662391B (zh) 一种胃磁信号的特征提取与分类系统及方法
CN114692682A (zh) 一种基于图嵌入表示的运动想象分类方法及系统
Sury et al. Brain-Computer Interface of Motor Imagery and Emotion Using Multiple Recurrent Neural Networks
Wang et al. An improved Generative Adversarial Network for Denoising EEG signals of brain-computer interface systems
CN109325402A (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
CP01 Change in the name or title of a patent holder

Address after: Second floor, building a, Tsinghua campus, Shenzhen University Town, Xili street, Nanshan District, Shenzhen City, Guangdong Province

Patentee after: Tsinghua Shenzhen International Graduate School

Patentee after: Dongguan JIANDA INFORMATION Technology Co.,Ltd.

Address before: Second floor, building a, Tsinghua campus, Shenzhen University Town, Xili street, Nanshan District, Shenzhen City, Guangdong Province

Patentee before: GRADUATE SCHOOL AT SHENZHEN, TSINGHUA University

Patentee before: DONGGUAN JIANDA INFORMATION TECHNOLOGY Co.,Ltd.

CP01 Change in the name or title of a patent holder
PP01 Preservation of patent right

Effective date of registration: 20240328

Granted publication date: 20220930

PP01 Preservation of patent right