CN110522443B - 基于心电信号的房室传导阻滞检测方法、装置和电子设备 - Google Patents

基于心电信号的房室传导阻滞检测方法、装置和电子设备 Download PDF

Info

Publication number
CN110522443B
CN110522443B CN201910817863.0A CN201910817863A CN110522443B CN 110522443 B CN110522443 B CN 110522443B CN 201910817863 A CN201910817863 A CN 201910817863A CN 110522443 B CN110522443 B CN 110522443B
Authority
CN
China
Prior art keywords
atrioventricular conduction
conduction block
interval
power spectral
spectral density
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910817863.0A
Other languages
English (en)
Other versions
CN110522443A (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.)
Guangzhou Shiyuan Electronics Thecnology Co Ltd
Original Assignee
Guangzhou Shiyuan Electronics Thecnology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangzhou Shiyuan Electronics Thecnology Co Ltd filed Critical Guangzhou Shiyuan Electronics Thecnology Co Ltd
Priority to CN201910817863.0A priority Critical patent/CN110522443B/zh
Publication of CN110522443A publication Critical patent/CN110522443A/zh
Application granted granted Critical
Publication of CN110522443B publication Critical patent/CN110522443B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7282Event detection, e.g. detecting unique waveforms indicative of a medical condition

Landscapes

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

Abstract

本申请涉及一种基于心电信号的房室传导阻滞检测方法、装置和电子设备。包括:在心电信号中,提取信号波形;获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种;将所述房室传导阻滞特征输入至检测器;获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果;所述信号类型包括房室传导阻滞和非房室传导阻滞。采用本方案,可以更准确地识别出房室传导阻滞,避免因个体差异不同错将非房室传导阻滞识别为房室传导阻滞的情况,提升了房室传导阻滞检测的准确性。

Description

基于心电信号的房室传导阻滞检测方法、装置和电子设备
技术领域
本申请涉及医疗器械与医疗产品领域,特别是涉及一种基于心电信号的房室传导阻滞检测方法、装置、电子设备和存储介质。
背景技术
房室传导阻滞(atrioventricular block,AVB)是一类心脏电传导系统阻滞疾病,主要包括第一级阻滞、第二级阻滞和第三极阻滞,阻滞级别由心房到心室的传导严重程度确定,而这类心房心室去极化不匹配给出了心电图变化,因此可通过心电特征提取来识别房室传导阻滞。
但是,由于房室传导阻滞反映在心电图上识别较为复杂,仅靠手工提取方式提取少数特征,并结合传统机器学习方法,在可靠性和准确性上都难以达到临床诊断标准。
因此,现有技术中的多房室传导阻滞检测方法存在着检测准确度低的问题。
发明内容
基于此,有必要针对上述技术问题,提供一种基于心电信号的房室传导阻滞检测方法、装置、电子设备和存储介质。
第一方面,提供了一种基于心电信号的房室传导阻滞检测方法,包括:
在心电信号中,提取信号波形;
获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种;
将所述房室传导阻滞特征输入至检测器;
获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果;所述信号类型包括房室传导阻滞和非房室传导阻滞。
上述基于心电信号的房室传导阻滞检测方法中,由于房室传导阻滞特征可以反映出房室传导阻滞的波形特点,因此考虑采用房室传导阻滞特征进行信号类型甄别,可以更准确地识别出房室传导阻滞,避免因个体差异不同错将非房室传导阻滞识别为房室传导阻滞的情况,提升了房室传导阻滞检测的准确性。
在一个实施例中,所述检测器包括密度检测器、时间检测器、变异检测器和综合检测器,所述将所述房室传导阻滞特征输入至检测器,包括:
若所述房室传导阻滞特征包括所述功率谱密度关联值,则将所述功率谱密度关联值输入至所述密度检测器;若所述房室传导阻滞特征包括所述PR间期持续时间,则将所述PR间期持续时间输入至所述时间检测器;若所述房室传导阻滞特征包括所述PR间期变异参数,则将所述PR间期变异参数输入至所述变异检测器;若所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数,则将所述房室传导阻滞特征输入至所述综合检测器。
上述实施例中,不同的房室传导阻滞特征对应识别房室传导阻滞会存在不同的检测标准,因此为不同的房室传导阻滞特征设置不同的独立检测器,使之针对性地进行特征识别、信号检测,大大提升了房室传导阻滞检测的准确性。
在一个实施例中,所述检测结果包括密度关联检测结果、持续时间检测结果以及变异性检测结果,所述获取所述检测器输出的检测结果,得到所述心电信号的信号类型,包括:
统计相匹配结果的结果数量;所述相匹配结果包括房室传导阻滞和非房室传导阻滞;若所述房室传导阻滞的结果数量大于所述非房室传导阻滞的结果数量,则确定所述房室传导阻滞,作为所述心电信号的信号类型;若所述房室传导阻滞的结果数量小于所述非房室传导阻滞的结果数量,则确定所述非房室传导阻滞,作为所述心电信号的信号类型。
上述实施例中,各个独立检测器依据各房室传导阻滞特征进行分析,并输出得到各房室传导阻滞特征的检测结果之后,将统计相同检测结果的结果数量,以此判定检测结果存在相同数量多的即为心电信号的信号类型。相对独立检测器所检测的结果而言,该统计判定的方式进一步增强了对最终信号类型确定的可靠性。
在一个实施例中,在所述在心电信号中,提取信号波形之前,还包括:
采集原始心电信号;对所述原始心电信号进行低通滤波,得到高频噪声滤除信号,作为所述心电信号;
所述在心电信号中,提取信号波形,包括:通过小波变换技术,在所述心电信号中提取P波和QRS波,得到所述信号波形。
上述实施例中,通过对原始心电信号进行低通滤波,可避免心电信号受噪声干扰,从而提高对信号类型确定的准确性。
在一个实施例中,当所述房室传导阻滞特征包括所述功率谱密度关联值时,所述获取所述信号波形的房室传导阻滞特征,包括:
获取第一功率谱密度,以及,获取第二功率谱密度;所述第一功率谱密度为房室传导阻滞心电信号的功率谱密度;所述第二功率谱密度为待测心电信号的功率谱密度;计算功率谱密度欧范积,以及,计算交叉功率谱密度;所述交叉功率谱密度为所述第一功率谱密度与所述第二功率谱密度的交叉功率谱密度;所述功率谱密度欧范积为所述第一功率谱密度的欧几里得范数与所述第二功率谱密度的欧几里得范数之积;计算所述交叉功率谱密度与所述功率谱密度欧范积的比值,作为所述功率谱密度关联值。
上述实施例中,房室传导阻滞特征可通过功率谱密度关联值体现,因此要获取房室传导阻滞特征,可先计算功率谱密度关联值,利用功率谱密度关联值的大小变化分析房室传导阻滞的有无,使房室传导阻滞可实际追溯,实施性强。
在一个实施例中,当所述房室传导阻滞特征包括所述PR间期持续时间时,所述获取所述信号波形的房室传导阻滞特征,包括:
获取P波序列,以及,获取R波序列;计算所述R波序列与所述P波序列的差值,得到PR间期持续时间分子;计算所述PR间期持续时间分子与预设采样频率的比值,作为所述PR间期持续时间。
上述实施例中,房室传导阻滞特征可通过PR间期持续时间体现,即PR间隔的时长,因此要获取房室传导阻滞特征,可先计算PR间期持续时间,利用PR间期持续时间的大小变化分析房室传导阻滞的有无,使房室传导阻滞可实际追溯,实施性强。
在一个实施例中,当所述房室传导阻滞特征包括所述PR间期变异参数时,所述获取所述信号波形的房室传导阻滞特征,包括:
获取PR间期序列;计算所述PR间期的相位空间,得到PR间期相位空间;当所述PR间期相位空间为PR间期延迟相位时,计算所述PR间期延迟相位的欧几里得距离,得到PR间期欧几里得距离;根据所述PR间期欧几里得距离,获取所述PR间期变异参数。
上述实施例中,房室传导阻滞特征可通过PR间期变异参数体现,因此要获取房室传导阻滞特征,可先计算PR间期变异参数,利用PR间期变异参数大小变化分析房室传导阻滞的有无,使房室传导阻滞可实际追溯,实施性强。
在一个实施例中,当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,所述获取所述信号波形的房室传导阻滞特征,包括:
计算所述功率谱密度关联值、所述PR间期持续时间与所述PR间期变异参数之和,得到综合指标分母;计算所述功率谱密度关联值、所述PR间期持续时间与所述PR间期变异参数之积,得到综合指标分子;计算所述综合指标分子与所述综合指标分母的比值,得到综合指标特征;确定所述综合指标特征,作为所述房室传导阻滞特征。
上述实施例中,由于功率谱密度关联值、PR间期持续时间与PR间期变异参数均可体现出信号波形中所存在的房室传导阻滞,因此综合三者之间的关联性,设计一个综合指标特征,利用综合指标特征反映房室传导阻滞,相当于将房室传导阻滞在信号波形中的特征放大,即使是在微弱信号下,也能够检测出房室传导阻滞,使房室传导阻滞的检测灵敏度增大。
第二方面,提供了一种基于心电信号的房室传导阻滞检测装置,包括:
信号波形提取模块,用于在心电信号中,提取信号波形;
阻滞特征获取模块,用于获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种;
阻滞特征输入模块,用于将所述房室传导阻滞特征输入至检测器;
信号类型确定模块,用于获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果;所述信号类型包括房室传导阻滞和非房室传导阻滞。
上述基于心电信号的房室传导阻滞检测装置中,采用房室传导阻滞特征进行信号类型甄别,可以更准确地识别出房室传导阻滞,避免因个体差异不同错将非房室传导阻滞识别为房室传导阻滞的情况,提升了房室传导阻滞检测的准确性。
第三方面,提供了一种电子设备,其特征在于,包括:存储器、以及一个或多个处理器;
所述存储器,用于存储一个或多个程序;
在心电信号中,提取信号波形;
获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种;
将所述房室传导阻滞特征输入至检测器;
获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果;所述信号类型包括房室传导阻滞和非房室传导阻滞。
上述电子设备中,采用房室传导阻滞特征进行信号类型甄别,可以更准确地识别出房室传导阻滞,避免因个体差异不同错将非房室传导阻滞识别为房室传导阻滞的情况,提升了房室传导阻滞检测的准确性。
第四方面,提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
在心电信号中,提取信号波形;
获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种;
将所述房室传导阻滞特征输入至检测器;
获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果;所述信号类型包括房室传导阻滞和非房室传导阻滞。
上述计算机可读存储介质中,采用房室传导阻滞特征进行信号类型甄别,可以更准确地识别出房室传导阻滞,避免因个体差异不同错将非房室传导阻滞识别为房室传导阻滞的情况,提升了房室传导阻滞检测的准确性。
附图说明
图1是本申请实施例一提供的一种房室传导阻滞检测方法的流程图;
图2是一个实施例的一种包含房室传导阻滞的心电信号示意图;
图3是一种P波、QRS波、T波、U波的信号波形的示意图;
图4是本申请实施例二提供的一种房室传导阻滞检测方法的流程图;
图5是本申请实施例三提供的一种房室传导阻滞检测装置的结构示意图;
图6是本申请实施例四提供的一种电子设备的结构示意图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
实施例一
图1是本申请实施例一提供的一种基于心电信号的房室传导阻滞检测方法的流程图,该房室传导阻滞检测方法,可以具体包括以下步骤:
步骤S110,在心电信号中,提取信号波形。
具体实现中,可以利用多通道同步数据,采集人体心脏信号、背景噪声和心电信号。更具体地,首先,可以通过心电导联和传感器采集心电信号,通过模拟电路对采集的心电信号进行阻抗匹配、过滤、放大等处理。然后,由模数转换器,将人体生理参数的模拟信号转化为数字信号。然后,通过低通滤波技术获取滤波后的心电信号。最后,通过小波变换技术从滤波后的心电信号中提取出信号波形,该信号波形如图2所示。
需要说明的是,房室传导阻滞的波形特征体现在心电图中时,存在以下波形特点:(1)I度房室传导阻滞,PR间隔>200ms;(2)II度房室传导阻滞,PR间隔>200ms,且偶尔有P波无QRS波;(3)III度房室传导阻滞,P波与P波间隔固,QRS与QRS间隔固定。
基于房室传导阻滞的波形特点,在得到如图2所示包含房室传导阻滞的心电信号示意图之后,本申请提出采用Daubechies4(Db4)小波变换对步骤S110中进行预处理后的心电信号进行离散小波分解,分解层数为N。在本实施例中,优选N=7,得到小波分解系数D1、D2、D3、D4、D5、D6、D7和A7,由于小波系数D3、D4和D5包含最大信息和最大能量,因此采用小波系数D3、D4和D5进行QRS复合波检测,检测得到的QRS复合波如图3所示。
图3是一种P波、QRS波、T波、U波的信号波形的示意图。如图所示,以时间t作为X轴,电压V作为Y轴,构建一个XY坐标系,根据心电信号中在每个时间点tn对应的电压Vn,生成心电信号的信号波形。
步骤S120,获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种。
其中,功率谱密度关联值可以反映房室传导阻滞心电信号与待测心电信号的相关程度。
其中,PR间期持续时间可以反映PR间期的波形持续时间。
其中,PR间期变异参数可以反映PR间期的离散程度。
具体实现中,可以获取信号波形的功率谱密度关联值、PR间期持续时间以及PR间期变异参数,作为信号波形的房室传导阻滞特征。
获取信号波形的功率谱密度关联值的步骤可以是:首先计算心电信号的功率谱密度(PSD),表示为S,功率谱密度S的计算公式如下:
S=G(ω)=∑Pnδ(ω-nΩ0)
其中,nΩ0表示周期信号中频率为nΩ0的谐波分量,Pn=An2/2是频率为nΩ0的谐波分量的功率,An是周期信号中频率为nΩ0的谐波分量的幅值,δ为预设参数,Ω0=2π/T0表示基频,即间隔长度,G(ω)的物理意义为频率ω上的功率密度,即功率谱密度S。
其次,在计算得到功率谱密度之后,可计算参考信号的功率谱密度Sx1,即参考信号为房室传导阻滞心电信号的功率谱密度Sx1,以及,计算目标信号的功率谱密度Sx2,即目标信号为待测心电信号的功率谱密度Sx2,计算得到Sx1和Sx2之后,可计算目标信号与参考信号的交叉功率谱密度Sx1x2
最后,采用Welch方法分析目标信号和参考信号的相关性,即可得到功率谱密度关联值Cx1x2,功率谱密度关联值Cx1x2越高,表示目标信号与房室传导阻滞的心电图波形越相似,反之则表示房室传导阻滞不存在,该功率谱密度关联值Cx1x2的计算公式如下所示:
Figure BDA0002186789440000081
获取信号波形的PR间期持续时间的步骤可以是:考虑到发生房室传导阻滞时,PR间期>200ms,该现象可以用PR间期持续时间PR_interval来表示,因此令X_P、X_Q、X_R、X_S和X_T分别表示P、Q、R、S和T波的序列,fs表示为采样频率,则PR间期持续时间PR_interval的计算公式如下所示:
PR_interval=(X_R-X_P)/fs
获取信号波形的PR间期变异参数的步骤可以是:考虑到发生房室传导阻滞时,偶尔有P波而无QRS波,该现象可以用PR间期离散程度体现,即PR间期变异参数PRIV,因此令x(n),n=1,2……m为PR间期的序列,则PR间期变异参数PRIV的计算公式如下所示:
y(n)=(x(n),x(n+1),……x(n+(m-1)t)),n=1,2,……,m
Figure BDA0002186789440000082
其中,y(n)表示为PR间期的相位空间,||.||表示为欧几里得距离,h为阶梯函数,m为嵌入维数,t为延迟时间,C为组合运算,r为预设参数。
步骤S130,将所述房室传导阻滞特征输入至检测器。
其中,检测器可以是针对单个房室传导阻滞特征所设计的独立检测器,也可以是针对多个房室传导阻滞特征所设计的综合检测器。
具体实现中,将单独的房室传导阻滞特征输入至对应的独立检测器之外,还可将多个所获取到的房室传导阻滞特征输入至综合检测器,以便根据多个房室传导阻滞特征组合甄别以实现对房室传导阻滞的判断。
步骤S140,获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果;所述信号类型包括房室传导阻滞和非房室传导阻滞。
具体实现中,可分别获取不同检测器所输出的检测结果,该检测结果可以是根据功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种特征所检测的结果,也可以是根据所有房室传导阻滞特征所检测的结果。心电信号的信号类型可根据单个或多个检测结果进行确定,以确保在心电信号微弱的情况下,仍可检测出房室传导阻滞。
本申请实施例提供的技术方案,首先在心电信号中提取信号波形,然后获取信号波形的房室传导阻滞特征,将房室传导阻滞特征输入至检测器之后,得到检测器输出的检测结果,进而根据检测结果确定心电信号的信号类型。因为房室传导阻滞特征可以反映出房室传导阻滞的波形特点,采用房室传导阻滞特征进行信号类型甄别,可以更准确地识别出房室传导阻滞,避免因个体差异不同错将非房室传导阻滞识别为房室传导阻滞的情况,提升了房室传导阻滞检测的准确性。
实施例二
图4是本申请实施例二提供的一种房室传导阻滞检测方法的流程图。具体的,参考图4,本申请实施例二的房室传导阻滞检测方法,具体包括:
步骤S210,采集原始心电信号,并对所述原始心电信号进行低通滤波,得到高频噪声滤除信号,作为所述心电信号。
具体实现中,可以采用低通数字滤波器进行低通滤波,滤除高频噪声(如300Hz以上),得到滤波后的心电信号。其中,低通数字滤波器可以具体为巴特沃斯滤波器。
步骤S220,在心电信号中,提取信号波形。
具体实现中,通过对原始心电信号进行低通滤波预处理,可避免心电信号受噪声干扰,从而提高对信号类型确定的准确性。
进一步地,所述步骤S220,包括:通过小波变换技术,在所述心电信号中提取P波、QRS波和T波,得到所述信号波形。
具体实现中,采用小波变换技术,提取心电信号中的P波、QRS波和T波的波形信息。
步骤S230,获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种。
步骤S240,将所述房室传导阻滞特征输入至检测器。
步骤S250,获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果;所述信号类型包括房室传导阻滞和非房室传导阻滞。
上述实施例中,利用小波变换技术提取心电信号中的信号波形,可确保波形稳定性,提高信号波形的可靠性。
在一个实施例中,所述检测器包括密度检测器、时间检测器、变异检测器和综合检测器,所述步骤S240,包括:
若所述房室传导阻滞特征包括所述功率谱密度关联值,则将所述功率谱密度关联值输入至所述密度检测器;若所述房室传导阻滞特征包括所述PR间期持续时间,则将所述PR间期持续时间输入至所述时间检测器;若所述房室传导阻滞特征包括所述PR间期变异参数,则将所述PR间期变异参数输入至所述变异检测器;若所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数,则将所述房室传导阻滞特征输入至所述综合检测器。
具体实现中,若房室传导阻滞特征仅包括功率谱密度关联值,则仅需将该特征输入至密度监测器,若房室传导阻滞特征仅包括PR间期持续时间,则仅需将该特征输入至时间检测器,若房室传导阻滞特征仅包括PR间期变异参数,则仅需将该特征输入至变异检测器,若房室传导阻滞特征包括上述特征中的至少两个,则不仅需将多个特征单独地输入至对应的检测器,还需将多个特征一同输入至综合检测器,以确保在提取到不同房室传导阻滞特征的情况下,均可获取检测结果,并提高检测结果的可靠性。
上述实施例中,不同房室传导阻滞特征对应识别房室传导阻滞会存在不同的检测标准,因此为不同房室传导阻滞特征设置不同的独立检测器,使之针对性地进行特征识别、信号检测,大大提升了左束支传导阻滞检测的准确性。同时,为多个不同的房室传导阻滞特征设置综合检测器,利用综合检测器同时分析多个房室传导阻滞特征,放大房室传导阻滞信号的检测灵敏度,这在上述实施例仅设置独立检测器的基础上又增强了最终输出结果的稳定性,使得房室传导阻滞的检测结果可靠性更高。
在一个实施例中,所述检测结果包括密度关联检测结果、持续时间检测结果以及变异性检测结果,所述根据所述检测结果,确定所述心电信号的信号类型,包括:
统计相匹配结果的结果数量;所述相匹配结果包括房室传导阻滞和非房室传导阻滞;若所述房室传导阻滞的结果数量大于所述非房室传导阻滞的结果数量,则确定所述房室传导阻滞,作为所述心电信号的信号类型;若所述房室传导阻滞的结果数量小于所述非房室传导阻滞的结果数量,则确定所述非房室传导阻滞,作为所述心电信号的信号类型。
其中,结果数量是指相同信号类型的结果数量,例如,在三种单独检测结果中,存在两种检测结果均为“房室传导阻滞”、而另一种检测结果为“非房室传导阻滞”,则“房室传导阻滞”的结果数量为二,“非房室传导阻滞”的结果数量为一。
具体实现中,确定心电信号的信号类型,不仅需判断三种独立检测结果,还需结合综合检测结果进行判断,即将功率谱密度关联值、PR间期持续时间和PR间期变异参数等房室传导阻滞特征均输入至综合检测器之后,得到不同的独立检测结果,利用多个独立检测结果进行综合判断,最终可确定以统计结果数量最多的类型作为最终心电信号的信号类型。
上述实施例中,各个独立检测器依据各房室传导阻滞特征进行分析,并输出得到各房室传导阻滞特征的检测结果之后,将统计相同检测结果的结果数量,以此判定检测结果存在相同数量多的即为心电信号的信号类型。相对独立检测器所检测的结果而言,该统计判定的方式进一步增强了对最终信号类型确定的可靠性。
在一个实施例中,当所述房室传导阻滞特征包括所述功率谱密度关联值时,所述步骤S230,包括:
获取第一功率谱密度,以及,获取第二功率谱密度;所述第一功率谱密度为房室传导阻滞心电信号的功率谱密度;所述第二功率谱密度为待测心电信号的功率谱密度;计算功率谱密度欧范积,以及,计算交叉功率谱密度;所述交叉功率谱密度为所述第一功率谱密度与所述第二功率谱密度的交叉功率谱密度;所述功率谱密度欧范积为所述第一功率谱密度的欧几里得范数与所述第二功率谱密度的欧几里得范数之积;计算所述交叉功率谱密度与所述功率谱密度欧范积的比值,作为所述功率谱密度关联值。
其中,第一功率谱密度可以是参考信号设定为房室传导阻滞心电信号功率谱密度Sx1
其中,第二功率谱密度可以是目标信号设定为待测心电信号的功率谱密度Sx2
其中,功率谱密度欧范积可以是Sx1的欧几里得范数||Sx1||2,与Sx2的欧几里得范数||Sx2||2的点积。
具体实现中,功率谱密度关联值Cx1x2的计算公式可参考上述实施例一中所示公式,在此不再赘述。
上述实施例中,房室传导阻滞特征可通过功率谱密度关联值体现,因此要获取房室传导阻滞特征,可先计算功率谱密度关联值,利用功率谱密度关联值的大小变化分析房室传导阻滞的有无,使房室传导阻滞可实际追溯,实施性强。
在一个实施例中,所述获取第一功率谱密度,以及,获取第二功率谱密度,包括:
获取心电信号的谐波分量差值,以及,获取心电信号的谐波分量功率;所述谐波分量差值为预设频率与谐波分量的差值;所述谐波分量功率为所述谐波分量的功率;计算所述谐波分量功率与所述谐波分量差值之和,得到功率谱密度;若所述心电信号为房室传导阻滞心电信号,则确定所述功率谱密度为所述第一功率谱密度;若所述心电信号为待测心电信号,则确定所述功率谱密度为所述第二功率谱密度。
其中,谐波分量差值可以是频率ω与谐波分量nΩ0的差值“ω-nΩ0”。
具体实现中,第一功率谱密度Sx1与第二功率谱密度Sx2的计算公式可参考上述实施例一中所示公式,在此不再赘述。
上述实施例中,第一功率谱密度以及第二功率谱密度均为计算功率谱密度关联值的基础参数,因此要获取功率谱密度关联值,首先需要通过对房室传导阻滞心电信号的分析获取第一功率谱密度,继而引入作为待测心电信号的功率谱密度,两者的相似性比较可进一步通过功率谱密度关联值的大小体现,即可得到待测心电信号是否与房室传导阻滞心电信号相似,进而检测出当前心电信号中是否存在有房室传导阻滞。
在一个实施例中,当所述房室传导阻滞特征包括所述功率谱密度关联值时,所述检测器包括密度检测器,所述步骤S250,包括:
在所述密度检测器中,提取M个导联的功率谱密度关联值,得到多导联功率谱密度关联值;M≥1;计算所述多导联功率谱密度关联值的平均值,得到多导联功率谱密度关联平均值;根据所述多导联功率谱密度关联平均值,确定所述心电信号的信号类型。
具体实现中,由于心电数据可由多个导联设定获取,如十二导联I、Ⅱ、Ⅲ、aVR、aVL、aVF、V1-V6,因此本实施例中设定导联数为M,则将提取到M个功率谱密度关联值Cx1x2,进而计算其平均值Cmean,由多导联功率谱密度关联平均值Cmean确定心电信号的信号类型。
更具体地,利用多导联功率谱密度关联平均值Cmean确定信号类型的方式可以是通过下述公式计算得出:
Figure BDA0002186789440000141
其中,Cs表示第s个导联的功率谱密度关联值Cx1x2值,Cmeanq表示连续监测过程中第q个计算得到的Cmean值,设定Cmean的初始阈值优选为CmeanD0,则CmeanDq表示第q个功率谱密度关联阈值。
当实时检测得到的Cmeanq值大于CmeanDq时,可确定该信号类型为房室传导阻滞,反之,即可认为该信号类型为非房室传导阻滞。
上述实施例中,考虑在多导联设置获取得到多个功率谱密度关联值的情况下,需针对多个功率谱密度关联值进行综合分析,则计算多个功率谱密度关联值的平均值,根据多导联功率谱密度关联平均值的大小判定心电信号是否存在房室传导阻滞,检测结果可靠性强。
在一个实施例中,所述根据所述多导联功率谱密度关联平均值,确定所述心电信号的信号类型,包括:
确定所述多导联功率谱密度关联平均值的历史获取数量;根据所述历史获取数量,确定所述多导联功率谱密度关联平均值的第一预设阈值;当所述多导联功率谱密度关联平均值大于所述第一预设阈值时,确定所述信号类型为房室传导阻滞。
其中,历史获取数量可以是多导联功率谱密度关联平均值的已获取数量,例如,3、5、7等。
具体实现中,将多导联功率谱密度关联平均值与其顺序对应的阈值进行比较,首先需确定多导联功率谱密度关联平均值的当前获取次序,进而获取相同次序对应的第一预设阈值,随即比较两者大小,若多导联功率谱密度关联平均值大于第一预设阈值,则可确定检测器输出的检测结果为房室传导阻滞。第一预设阈值的更新公式如下:
CmeanDq+1=λ1CmeanDq1CmeanDq_1
其中,CmeanDq+1表示为多导联功率谱密度关联平均值Cmean的第q+1个阈值、CmeanDq表示为Cmean的第q个阈值、CmeanDq-1表示为Cmean的第q-1个阈值,λ1表示为固定参数0.75、μ1表示为固定参数0.25,λ11=1。
上述实施例中,提出根据多导联功率谱密度关联平均值的历史获取数量来设置对应需匹配的“自适应阈值检测器”,“自适应”匹配各房室传导阻滞特征对应的阈值,由此提升房室传导阻滞检测的准确性。
在一个实施例中,当所述房室传导阻滞特征包括所述PR间期持续时间时,所述步骤S230,包括:
获取P波序列,以及,获取R波序列;计算所述R波序列与所述P波序列的差值,得到PR间期持续时间分子;计算所述PR间期持续时间分子与预设采样频率的比值,作为所述PR间期持续时间。
其中,预设采样频率表示为fs,优选实验最佳参数。
具体实现中,PR间期持续时间PR_interval的计算公式可参考上述实施例一中所示公式,在此不再赘述。
上述实施例中,房室传导阻滞特征可通过PR间期持续时间体现,即PR间隔的时长,因此要获取房室传导阻滞特征,可先计算PR间期持续时间,利用PR间期持续时间的大小变化分析房室传导阻滞的有无,使房室传导阻滞可实际追溯,实施性强。
在一个实施例中,当所述房室传导阻滞特征包括所述PR间期持续时间时,所述检测器包括时间检测器,所述步骤S250,包括:
在所述时间检测器中,提取M个导联的PR间期持续时间,得到多导联PR间期持续时间;M≥1;计算所述多导联PR间期持续时间的平均值,得到多导联PR间期持续时间平均值;根据所述多导联PR间期持续时间平均值,确定所述心电信号的信号类型。
具体实现中,利用多导联PR间期持续时间平均值PR_intervalmean确定信号类型的方式可以是通过下述公式计算得出:
Figure BDA0002186789440000151
其中,PR_intervals表示第s个导联的PR间期持续时间PR_interval,PR_intervalmeanq表示连续监测过程中第q个计算得到的多导联PR间期持续时间平均值PR_intervalmean。
上述实施例中,考虑在多导联设置获取得到多个PR间期持续时间的情况下,需针对多个PR间期持续时间进行综合分析,则计算多个PR间期持续时间的平均值,根据多导联PR间期持续时间平均值的大小判定心电信号是否存在房室传导阻滞,检测结果可靠性强。
在一个实施例中,所述根据所述多导联PR间期持续时间平均值,确定所述心电信号的信号类型,包括:
确定所述多导联PR间期持续时间平均值的历史获取数量;根据所述历史获取数量,确定所述多导联PR间期持续时间平均值的第二预设阈值;当所述多导联PR间期持续时间平均值大于所述第二预设阈值时,确定所述信号类型为房室传导阻滞。
其中,历史获取数量可以是多导联PR间期持续时间平均值的已获取数量,例如,3、5、7等。
具体实现中,将多导联PR间期持续时间平均值与其顺序对应的阈值进行比较,首先需确定多导联PR间期持续时间平均值的当前获取次序,进而获取相同次序对应的第二预设阈值,随即比较两者大小,若多导联PR间期持续时间平均值大于第二预设阈值,则可确定检测器输出的检测结果为房室传导阻滞。第二预设阈值的更新公式如下:
PR_intervalmeanDq+1=PR_intervalmeanD0=120ms
其中,设定PR_intervalmean的初始阈值为PR_intervalmeanD0=120ms(该阈值设定为成人阈值120ms,在其他实施例中可根据具体情况进行调整),PR_intervalmeanDq+1表示第q+1个PR_intervalmean的阈值。
上述实施例中,提出根据多导联PR间期持续时间平均值的历史获取数量来设置对应需匹配的“自适应阈值检测器”,“自适应”匹配各房室传导阻滞特征对应的阈值,由此提升房室传导阻滞检测的准确性。
在一个实施例中,当所述房室传导阻滞特征包括所述PR间期变异参数时,所述步骤S230,包括:
获取PR间期序列;计算所述PR间期的相位空间,得到PR间期相位空间;当所述PR间期相位空间为PR间期延迟相位时,计算所述PR间期延迟相位的欧几里得距离,得到PR间期欧几里得距离;根据所述PR间期欧几里得距离,获取所述PR间期变异参数。
其中,PR间期序列表示为x(n),PR间期相位空间表示为y(n),PR间期延迟相位表示为y(t),PR间期欧几里得距离表示为||y(t)||。
具体实现中,计算PR间期变异参数PRIV的公式可参考上述实施例一中所示计算公式,在此不再赘述,但需说明的是,虽然在本实施中PR间期欧几里得距离表示为||y(t)||,但实际计算PR间期变异参数PRIV,需首先计算相位空间差值yi(t)-yj(t),进而计算该相位空间差值的欧几里得距离。
上述实施例中,房室传导阻滞特征可通过PR间期变异参数体现,因此要获取房室传导阻滞特征,可先计算PR间期变异参数,利用PR间期变异参数大小变化分析房室传导阻滞的有无,使房室传导阻滞可实际追溯,实施性强。
在一个实施例中,所述根据所述PR间期欧几里得距离,获取所述PR间期变异参数,包括:
计算所述PR间期欧几里得距离的阶梯函数,得到PR间期变异分量分子;计算所述PR间期变异分量分子与预设组合参数的比值,得到PR间期变异分量;计算所述PR间期变异分量的多项分量之和,得到PR间期变异分量总和;确定所述PR间期变异分量总和,作为所述PR间期变异参数。
其中,预设组合参数在本实施例及上述实施例中表示为
Figure BDA0002186789440000171
具体实现中,计算PR间期欧几里得距离的阶梯函数,首先需计算PR间期欧几里得距离,即在实施例一中所示公式||yi(t)-yj(t)||,进而计算其阶梯函数h(r||yi(t)-yj(t)||),作为PR间期变异分量分子。计算PR间期变异分量的多项分量之和,即首先计算PR间期变异分量
Figure BDA0002186789440000172
的和
Figure BDA0002186789440000173
在此基础上进一步计算得到PR间期变异分量总和
Figure BDA0002186789440000174
上述实施例中,PR间期变异参数由PR间期欧几里得距离确定,因此要获取PR间期变异参数,首先需计算PR间期欧几里得距离的阶梯函数,将PR间期变异参数细致拆分为多个参数组成,对PR间期变异参数的采集更加精细化,使得房室传导阻滞的检测可靠性更高。
在一个实施例中,当所述房室传导阻滞特征包括所述PR间期变异参数时,所述检测器包括变异检测器,所述步骤S250,包括:
在所述变异检测器中,提取M个导联的PR间期变异参数,得到多导联PR间期变异参数;M≥1;计算所述多导联PR间期变异参数的平均值,得到多导联PR间期变异平均值;根据所述多导联PR间期变异平均值,确定所述心电信号的信号类型。
具体实现中,利用多导联PR间期变异平均值PRIVmean确定信号类型的方式可以是通过下述公式计算得出:
Figure BDA0002186789440000181
其中,PRIVs表示第s个导联的PR间期变异参数PRIV,PRIVmeanq表示连续监测过程中第q个计算得到的多导联PR间期变异平均值PRIVmean。
上述实施例中,考虑在多导联设置获取得到多个PR间期变异参数的情况下,需针对多个PR间期变异参数进行综合分析,则计算多个PR间期变异参数的平均值,根据多导联PR间期变异平均值的大小判定心电信号是否存在房室传导阻滞,检测结果可靠性强。
在一个实施例中,所述根据所述多导联PR间期变异平均值,确定所述心电信号的信号类型,包括:
确定所述多导联PR间期变异平均值的历史获取数量;根据所述历史获取数量,确定所述多导联PR间期变异平均值的第三预设阈值;当所述多导联PR间期变异平均值大于所述第三预设阈值时,确定所述信号类型为房室传导阻滞。
其中,历史获取数量可以是多导联PR间期变异平均值的已获取数量,例如,3、5、7等。
具体实现中,将多导联PR间期变异平均值与其顺序对应的阈值进行比较,首先需确定多导联PR间期变异平均值的当前获取次序,进而获取相同次序对应的第三预设阈值,随即比较两者大小,若多导联PR间期变异平均值大于第三预设阈值,则可确定检测器输出的检测结果为房室传导阻滞。第三预设阈值的更新公式如下:
PRIVmeanDq+1=λ2PRIVmeanDq2PRIVmeanDq_1
其中,设定PRIVmean的初始阈值为PRIVmeanD0,PRIVmeanDq+1表示第q+1个PRIVmean的阈值,PRIVmeanDq-1表示第q-1个PRIVmean的阈值,PRIVmeanDq表示第q个PRIVmean的阈值,λ2表示为固定参数0.85、μ2表示为固定参数0.15,λ22=1。
上述实施例中,提出根据多导联PR间期变异平均值的历史获取数量来设置对应需匹配的“自适应阈值检测器”,“自适应”匹配各房室传导阻滞特征对应的阈值,由此提升房室传导阻滞检测的准确性。
在一个实施例中,当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,所述步骤S230,包括:
计算所述功率谱密度关联值、所述PR间期持续时间与所述PR间期变异参数之和,得到综合指标分母;计算所述功率谱密度关联值、所述PR间期持续时间与所述PR间期变异参数之积,得到综合指标分子;计算所述综合指标分子与所述综合指标分母的比值,得到综合指标特征;确定所述综合指标特征,作为所述房室传导阻滞特征。
具体实现中,除利用上述独立检测器(密度检测器、时间检测器和变异检测器)获取检测结果之外,还可利用综合检测器获取输出的检测结果,该综合检测器需输入功率谱密度关联值Cx1x2、PR间期持续时间PR_interval以及PR间期变异参数PRIV方可输出检测结果。具体地,综合检测器输出检测结果所需利用的计算公式如下:
C_PR=Cx1,x2×PR_interval×PRIV/(Cx1,x2+PR_interval+PRIV)
其中,综合指标特征C_PR可以是房室传导阻滞特征的一个综合指标特征,该特征结合了所有房室传导阻滞特征(功率谱密度关联值、PR间期持续时间以及PR间期变异参数)的特点,即使在微弱信号下,也能够检测出房室传导阻滞。
上述实施例中,由于功率谱密度关联值、PR间期持续时间与PR间期变异参数均可体现出信号波形中所存在的房室传导阻滞,因此综合三者之间的关联性,设计一个综合指标特征,利用综合指标特征反映房室传导阻滞,相当于将房室传导阻滞在信号波形中的特征放大,即使是在微弱信号下,也能够检测出房室传导阻滞,使房室传导阻滞的检测灵敏度增大。
在一个实施例中,当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,所述检测器包括综合检测器,所述步骤S250,包括:
在所述综合检测器中,提取M个导联的综合指标特征,得到多导联综合指标特征;M≥1;计算所述多导联综合指标特征的平均值,得到多导联综合指标特征平均值;根据所述多导联综合指标特征平均值,确定所述心电信号的信号类型。
具体实现中,利用多导联综合指标特征平均值C_PRmean确定信号类型的方式可以是通过下述公式计算得出:
Figure BDA0002186789440000201
其中,C_PRs表示第s个导联的综合指标特征C_PR,C_PRmeanq表示连续监测过程中第q个计算得到的多导联综合指标特征平均值C_PRmean。
上述实施例中,功率谱密度关联值、PR间期持续时间和PR间期变异参数,均可作为房室传导阻滞的判定因素,但将其所有结合作为一个综合指标参数,可通过对综合指标参数的数值大小来判定当前出现的特征是否为房室传导阻滞特征,进而有效提升房室传导阻滞检测的准确性。
在一个实施例中,所述根据所述多导联综合指标特征平均值,确定所述心电信号的信号类型,包括:
确定所述多导联综合指标特征平均值的历史获取数量;根据所述历史获取数量,确定所述多导联综合指标特征平均值的第四预设阈值;当所述多导联综合指标特征平均值大于所述第四预设阈值时,确定所述信号类型为房室传导阻滞。
其中,历史获取数量可以是多导联综合指标特征平均值的已获取数量,例如,3、5、7等。
具体实现中,将多导联综合指标特征平均值与其顺序对应的阈值进行比较,首先需确定多导联综合指标特征平均值的当前获取次序,进而获取相同次序对应的第四预设阈值,随即比较两者大小,若多导联综合指标特征平均值大于第四预设阈值,则可确定检测器输出的检测结果为房室传导阻滞。第四预设阈值的更新公式如下:
C_PRmeanDq+1=λ3C_PRmeanDq3C_PRmeanDq_1
其中,设定C_PRmean的初始阈值为C_PRmeanD0,C_PRmeanDq+1表示第q+1个C_PRmean的阈值,C_PRmeanDq-1表示第q-1个C_PRmean的阈值,C_PRmeanDq表示第q个C_PRmean的阈值。
上述实施例中,提出根据多导联综合指标特征平均值的历史获取数量来设置对应需匹配的综合指标阈值,该综合指标阈值的作用在于可应用于上述多个实施例中的“自适应阈值检测器”,“自适应”匹配各综合指标参数对应的阈值,由此提升房室传导阻滞检测的准确性。
在一个实施例中,所述检测器包括卷积神经网络检测器、循环神经网络检测器、卷积循环神经网络检测器、深度神经网络检测器中的至少一种。
上述实施例中,提供了多种检测器可应用于本申请用以获取检测结果,增强了房室传导阻滞检测装置的可实施性。
本申请实施例提供的技术方案,由于不同导联的房室传导阻滞特征可以准确反映房室传导阻滞,针对不同导联的房室传导阻滞设定阻滞阈值,通过比较特征参数与预设阈值的大小关系,即可得到检测结果,进而确定心电信号的信号类型。采用本方法,可以更准确地识别出房室传导阻滞,避免因个体差异不同错将非房室传导阻滞识别为房室传导阻滞的情况,提升了房室传导阻滞检测的准确性。
应该理解的是,虽然图1、图4流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图1、图4中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
实施例三
图5是本申请实施例三提供的一种房室传导阻滞检测装置的结构示意图。参考图5,本实施例三提供的房室传导阻滞检测装置具体包括:信号波形提取模块510、阻滞特征获取模块520、阻滞特征输入模块530和信号类型确定模块540,其中:
信号波形提取模块510,用于在心电信号中,提取信号波形;
阻滞特征获取模块520,用于获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种;
阻滞特征输入模块530,用于将所述房室传导阻滞特征输入至检测器;
信号类型确定模块540,用于获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果;所述信号类型包括房室传导阻滞和非房室传导阻滞。
本申请实施例提供的技术方案,首先在心电信号中提取信号波形,然后获取信号波形的房室传导阻滞特征,将房室传导阻滞特征输入至检测器之后,得到检测器输出的检测结果,进而根据检测结果确定心电信号的信号类型。因为房室传导阻滞特征可以反映出房室传导阻滞的波形特点,采用房室传导阻滞特征进行信号类型甄别,可以更准确地识别出房室传导阻滞,避免因个体差异不同错将非房室传导阻滞识别为房室传导阻滞的情况,提升了房室传导阻滞检测的准确性。
在一个实施例中,所述检测器包括密度检测器、时间检测器、变异检测器和综合检测器,所述阻滞特征输入模块530,还用于若所述房室传导阻滞特征包括所述功率谱密度关联值,则将所述功率谱密度关联值输入至所述密度检测器;若所述房室传导阻滞特征包括所述PR间期持续时间,则将所述PR间期持续时间输入至所述时间检测器;若所述房室传导阻滞特征包括所述PR间期变异参数,则将所述PR间期变异参数输入至所述变异检测器;若所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数,则将所述房室传导阻滞特征输入至所述综合检测器。
在一个实施例中,所述检测结果包括密度关联检测结果、持续时间检测结果以及变异性检测结果,所述信号类型确定模块540,还用于将所述密度关联检测结果、所述持续时间检测结果和所述变异性检测结果进行匹配;统计相匹配结果的结果数量;所述相匹配结果包括房室传导阻滞和非房室传导阻滞;若所述房室传导阻滞的结果数量大于所述非房室传导阻滞的结果数量,则确定所述房室传导阻滞,作为所述心电信号的信号类型;若所述房室传导阻滞的结果数量小于所述非房室传导阻滞的结果数量,则确定所述非房室传导阻滞,作为所述心电信号的信号类型。
在一个实施例中,所述房室传导阻滞检测装置,还包括:
原始信号采集模块,用于采集原始心电信号,并对所述原始心电信号进行低通滤波,得到高频噪声滤除信号,作为所述心电信号;
在一个实施例中,所述信号波形提取模块510,还用于通过小波变换技术,在所述心电信号中提取P波和QRS波,得到所述信号波形。
在一个实施例中,当所述房室传导阻滞特征包括所述功率谱密度关联值时,所述阻滞特征获取模块520,还用于获取第一功率谱密度,以及,获取第二功率谱密度;所述第一功率谱密度为房室传导阻滞心电信号的功率谱密度;所述第二功率谱密度为待测心电信号的功率谱密度;计算功率谱密度欧范积,以及,计算交叉功率谱密度;所述交叉功率谱密度为所述第一功率谱密度与所述第二功率谱密度的交叉功率谱密度;所述功率谱密度欧范积为所述第一功率谱密度的欧几里得范数与所述第二功率谱密度的欧几里得范数之积;计算所述交叉功率谱密度与所述功率谱密度欧范积的比值,作为所述功率谱密度关联值。
在一个实施例中,所述阻滞特征获取模块520,还用于获取心电信号的谐波分量差值,以及,获取心电信号的谐波分量功率;所述谐波分量差值为预设频率与谐波分量的差值;所述谐波分量功率为所述谐波分量的功率;计算所述谐波分量功率与所述谐波分量差值之和,得到功率谱密度;若所述心电信号为房室传导阻滞心电信号,则确定所述功率谱密度为所述第一功率谱密度;若所述心电信号为待测心电信号,则确定所述功率谱密度为所述第二功率谱密度。
在一个实施例中,当所述房室传导阻滞特征包括所述功率谱密度关联值时,所述检测器包括密度检测器,所述信号类型确定模块540,还用于在所述密度检测器中,提取M个导联的功率谱密度关联值,得到多导联功率谱密度关联值;M≥1;计算所述多导联功率谱密度关联值的平均值,得到多导联功率谱密度关联平均值;根据所述多导联功率谱密度关联平均值,确定所述心电信号的信号类型。
在一个实施例中,所述信号类型确定模块540,还用于确定所述多导联功率谱密度关联平均值的历史获取数量;根据所述历史获取数量,确定所述多导联功率谱密度关联平均值的第一预设阈值;当所述多导联功率谱密度关联平均值大于所述第一预设阈值时,确定所述信号类型为房室传导阻滞。
在一个实施例中,当所述房室传导阻滞特征包括所述PR间期持续时间时,所述阻滞特征获取模块520,还用于获取P波序列,以及,获取R波序列;计算所述R波序列与所述P波序列的差值,得到PR间期持续时间分子;计算所述PR间期持续时间分子与预设采样频率的比值,作为所述PR间期持续时间。
在一个实施例中,当所述房室传导阻滞特征包括所述PR间期持续时间时,所述检测器包括时间检测器,所述信号类型确定模块540,还用于在所述时间检测器中,提取M个导联的PR间期持续时间,得到多导联PR间期持续时间;M≥1;计算所述多导联PR间期持续时间的平均值,得到多导联PR间期持续时间平均值;根据所述多导联PR间期持续时间平均值,确定所述心电信号的信号类型。
在一个实施例中,所述信号类型确定模块540,还用于确定所述多导联PR间期持续时间平均值的历史获取数量;根据所述历史获取数量,确定所述多导联PR间期持续时间平均值的第二预设阈值;当所述多导联PR间期持续时间平均值大于所述第二预设阈值时,确定所述信号类型为房室传导阻滞。
在一个实施例中,当所述房室传导阻滞特征包括所述PR间期变异参数时,所述阻滞特征获取模块520,还用于获取PR间期序列;计算所述PR间期的相位空间,得到PR间期相位空间;当所述PR间期相位空间为PR间期延迟相位时,计算所述PR间期延迟相位的欧几里得距离,得到PR间期欧几里得距离;根据所述PR间期欧几里得距离,获取所述PR间期变异参数。
在一个实施例中,所述阻滞特征获取模块520,还用于计算所述PR间期欧几里得距离的阶梯函数,得到PR间期变异分量分子;计算所述PR间期变异分量分子与预设组合参数的比值,得到PR间期变异分量;计算所述PR间期变异分量的多项分量之和,得到PR间期变异分量总和;确定所述PR间期变异分量总和,作为所述PR间期变异参数。
在一个实施例中,当所述房室传导阻滞特征包括所述PR间期变异参数时,所述检测器包括变异检测器,所述信号类型确定模块540,还用于在所述变异检测器中,提取M个导联的PR间期变异参数,得到多导联PR间期变异参数;M≥1;计算所述多导联PR间期变异参数的平均值,得到多导联PR间期变异平均值;根据所述多导联PR间期变异平均值,确定所述心电信号的信号类型。
在一个实施例中,所述信号类型确定模块540,还用于确定所述多导联PR间期变异平均值的历史获取数量;根据所述历史获取数量,确定所述多导联PR间期变异平均值的第三预设阈值;当所述多导联PR间期变异平均值大于所述第三预设阈值时,确定所述信号类型为房室传导阻滞。
在一个实施例中,当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,所述阻滞特征获取模块520,还用于计算所述功率谱密度关联值、所述PR间期持续时间与所述PR间期变异参数之和,得到综合指标分母;计算所述功率谱密度关联值、所述PR间期持续时间与所述PR间期变异参数之积,得到综合指标分子;计算所述综合指标分子与所述综合指标分母的比值,得到综合指标特征;确定所述综合指标特征,作为所述房室传导阻滞特征。
在一个实施例中,当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,所述检测器包括综合检测器,所述信号类型确定模块540,还用于在所述综合检测器中,提取M个导联的综合指标特征,得到多导联综合指标特征;M≥1;计算所述多导联综合指标特征的平均值,得到多导联综合指标特征平均值;根据所述多导联综合指标特征平均值,确定所述心电信号的信号类型。
在一个实施例中,所述信号类型确定模块540,还用于确定所述多导联综合指标特征平均值的历史获取数量;根据所述历史获取数量,确定所述多导联综合指标特征平均值的第四预设阈值;当所述多导联综合指标特征平均值大于所述第四预设阈值时,确定所述信号类型为房室传导阻滞。
在一个实施例中,所述检测器包括卷积神经网络检测器、循环神经网络检测器、卷积循环神经网络检测器、深度神经网络检测器中的至少一种。
本申请实施例提供的技术方案,由于不同导联的房室传导阻滞特征可以准确反映房室传导阻滞,针对不同导联的房室传导阻滞设定阻滞阈值,通过比较特征参数与预设阈值的大小关系,即可得到检测结果,进而确定心电信号的信号类型。采用本方法,可以更准确地识别出房室传导阻滞,避免因个体差异不同错将非房室传导阻滞识别为房室传导阻滞的情况,提升了房室传导阻滞检测的准确性。
实施例四
图6是本申请实施例四提供的一种电子设备的结构示意图。如图中所示,该电子设备包括:处理器60、存储器61、具有触摸功能的显示屏62、输入装置63、输出装置64以及通信装置65。该电子设备中处理器60的数量可以是一个或者多个,图中以一个处理器60为例。该电子设备中存储器61的数量可以是一个或者多个,图中以一个存储器61为例。该电子设备的处理器60、存储器61、显示屏62、输入装置63、输出装置64以及通信装置65可以通过总线或者其他方式连接,图中以通过总线连接为例。实施例中,电子设备可以是电脑,手机,平板,投影仪或交互智能平板等。实施例中,以电子设备为交互智能平板为例,进行描述。
存储器61作为一种计算机可读存储介质,可用于存储软件程序、计算机可执行程序以及模块,如本申请任意实施例所述的房室传导阻滞检测方法对应的程序指令/模块。存储器61可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序;存储数据区可存储根据设备的使用所创建的数据等。此外,存储器61可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中,存储器61可进一步包括相对于处理器60远程设置的存储器,这些远程存储器可以通过网络连接至设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
显示屏62为具有触摸功能的显示屏62,其可以是电容屏、电磁屏或者红外屏。一般而言,显示屏62用于根据处理器60的指示显示数据,还用于接收作用于显示屏62的触摸操作,并将相应的信号发送至处理器60或其他装置。可选的,当显示屏62为红外屏时,其还包括红外触摸框,该红外触摸框设置在显示屏62的四周,其还可以用于接收红外信号,并将该红外信号发送至处理器60或者其他设备。
通信装置65,用于与其他设备建立通信连接,其可以是有线通信装置和/或无线通信装置。
输入装置63可用于接收输入的数字或者字符信息,以及产生与电子设备的用户设置以及功能控制有关的键信号输入,还可以是用于获取图像的摄像头以及获取音频数据的拾音设备。输出装置64可以包括扬声器等音频设备。需要说明的是,输入装置63和输出装置64的具体组成可以根据实际情况设定。
处理器60通过运行存储在存储器61中的软件程序、指令以及模块,从而执行设备的各种功能应用以及数据处理,即实现上述的房室传导阻滞检测方法。
具体的,实施例中,处理器60执行存储器61中存储的一个或多个程序时,具体实现如下操作:
在心电信号中,提取信号波形;
获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种;
将所述房室传导阻滞特征输入至检测器;
获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果;所述信号类型包括房室传导阻滞和非房室传导阻滞。
在上述实施例的基础上,所述检测器包括密度检测器、时间检测器、变异检测器和综合检测器,一个或多个处理器60还实现如下操作:
若所述房室传导阻滞特征包括所述功率谱密度关联值,则将所述功率谱密度关联值输入至所述密度检测器;若所述房室传导阻滞特征包括所述PR间期持续时间,则将所述PR间期持续时间输入至所述时间检测器;若所述房室传导阻滞特征包括所述PR间期变异参数,则将所述PR间期变异参数输入至所述变异检测器;若所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数,则将所述房室传导阻滞特征输入至所述综合检测器。
在上述实施例的基础上,所述检测结果包括密度关联检测结果、持续时间检测结果以及变异性检测结果,一个或多个处理器60还实现如下操作:
将所述密度关联检测结果、所述持续时间检测结果和所述变异性检测结果进行匹配;统计相匹配结果的结果数量;所述相匹配结果包括房室传导阻滞和非房室传导阻滞;若所述房室传导阻滞的结果数量大于所述非房室传导阻滞的结果数量,则确定所述房室传导阻滞,作为所述心电信号的信号类型;若所述房室传导阻滞的结果数量小于所述非房室传导阻滞的结果数量,则确定所述非房室传导阻滞,作为所述心电信号的信号类型。
在上述实施例的基础上,一个或多个处理器60还实现如下操作:
采集原始心电信号;对所述原始心电信号进行低通滤波,得到高频噪声滤除信号,作为所述心电信号;所述在心电信号中,提取信号波形,包括:通过小波变换技术,在所述心电信号中提取P波和QRS波,得到所述信号波形。
在上述实施例的基础上,当所述房室传导阻滞特征包括所述功率谱密度关联值时,一个或多个处理器60还实现如下操作:
获取第一功率谱密度,以及,获取第二功率谱密度;所述第一功率谱密度为房室传导阻滞心电信号的功率谱密度;所述第二功率谱密度为待测心电信号的功率谱密度;计算功率谱密度欧范积,以及,计算交叉功率谱密度;所述交叉功率谱密度为所述第一功率谱密度与所述第二功率谱密度的交叉功率谱密度;所述功率谱密度欧范积为所述第一功率谱密度的欧几里得范数与所述第二功率谱密度的欧几里得范数之积;计算所述交叉功率谱密度与所述功率谱密度欧范积的比值,作为所述功率谱密度关联值。
在上述实施例的基础上,一个或多个处理器60还实现如下操作:
获取心电信号的谐波分量差值,以及,获取心电信号的谐波分量功率;所述谐波分量差值为预设频率与谐波分量的差值;所述谐波分量功率为所述谐波分量的功率;计算所述谐波分量功率与所述谐波分量差值之和,得到功率谱密度;若所述心电信号为房室传导阻滞心电信号,则确定所述功率谱密度为所述第一功率谱密度;若所述心电信号为待测心电信号,则确定所述功率谱密度为所述第二功率谱密度。
在上述实施例的基础上,当所述房室传导阻滞特征包括所述功率谱密度关联值时,所述检测器包括密度检测器,一个或多个处理器60还实现如下操作:
在所述密度检测器中,提取M个导联的功率谱密度关联值,得到多导联功率谱密度关联值;M≥1;计算所述多导联功率谱密度关联值的平均值,得到多导联功率谱密度关联平均值;根据所述多导联功率谱密度关联平均值,确定所述心电信号的信号类型。
在上述实施例的基础上,一个或多个处理器60还实现如下操作:
确定所述多导联功率谱密度关联平均值的历史获取数量;根据所述历史获取数量,确定所述多导联功率谱密度关联平均值的第一预设阈值;当所述多导联功率谱密度关联平均值大于所述第一预设阈值时,确定所述信号类型为房室传导阻滞。
在上述实施例的基础上,当所述房室传导阻滞特征包括所述PR间期持续时间时,一个或多个处理器60还实现如下操作:
获取P波序列,以及,获取R波序列;计算所述R波序列与所述P波序列的差值,得到PR间期持续时间分子;计算所述PR间期持续时间分子与预设采样频率的比值,作为所述PR间期持续时间。
在上述实施例的基础上,当所述房室传导阻滞特征包括所述PR间期持续时间时,所述检测器包括时间检测器,一个或多个处理器60还实现如下操作:
在所述时间检测器中,提取M个导联的PR间期持续时间,得到多导联PR间期持续时间;M≥1;计算所述多导联PR间期持续时间的平均值,得到多导联PR间期持续时间平均值;根据所述多导联PR间期持续时间平均值,确定所述心电信号的信号类型。
在上述实施例的基础上,一个或多个处理器60还实现如下操作:
确定所述多导联PR间期持续时间平均值的历史获取数量;根据所述历史获取数量,确定所述多导联PR间期持续时间平均值的第二预设阈值;当所述多导联PR间期持续时间平均值大于所述第二预设阈值时,确定所述信号类型为房室传导阻滞。
在上述实施例的基础上,当所述房室传导阻滞特征包括所述PR间期变异参数时,一个或多个处理器60还实现如下操作:
获取PR间期序列;计算所述PR间期的相位空间,得到PR间期相位空间;当所述PR间期相位空间为PR间期延迟相位时,计算所述PR间期延迟相位的欧几里得距离,得到PR间期欧几里得距离;根据所述PR间期欧几里得距离,获取所述PR间期变异参数。
在上述实施例的基础上,一个或多个处理器60还实现如下操作:
计算所述PR间期欧几里得距离的阶梯函数,得到PR间期变异分量分子;计算所述PR间期变异分量分子与预设组合参数的比值,得到PR间期变异分量;计算所述PR间期变异分量的多项分量之和,得到PR间期变异分量总和;确定所述PR间期变异分量总和,作为所述PR间期变异参数。
在上述实施例的基础上,当所述房室传导阻滞特征包括所述PR间期变异参数时,所述检测器包括变异检测器,一个或多个处理器60还实现如下操作:
在所述变异检测器中,提取M个导联的PR间期变异参数,得到多导联PR间期变异参数;M≥1;计算所述多导联PR间期变异参数的平均值,得到多导联PR间期变异平均值;根据所述多导联PR间期变异平均值,确定所述心电信号的信号类型。
在上述实施例的基础上,一个或多个处理器60还实现如下操作:
确定所述多导联PR间期变异平均值的历史获取数量;根据所述历史获取数量,确定所述多导联PR间期变异平均值的第三预设阈值;当所述多导联PR间期变异平均值大于所述第三预设阈值时,确定所述信号类型为房室传导阻滞。
在上述实施例的基础上,当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,一个或多个处理器60还实现如下操作:
计算所述功率谱密度关联值、所述PR间期持续时间与所述PR间期变异参数之和,得到综合指标分母;计算所述功率谱密度关联值、所述PR间期持续时间与所述PR间期变异参数之积,得到综合指标分子;计算所述综合指标分子与所述综合指标分母的比值,得到综合指标特征;确定所述综合指标特征,作为所述房室传导阻滞特征。
在上述实施例的基础上,当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,所述检测器包括综合检测器,一个或多个处理器60还实现如下操作:
在所述综合检测器中,提取M个导联的综合指标特征,得到多导联综合指标特征;M≥1;计算所述多导联综合指标特征的平均值,得到多导联综合指标特征平均值;根据所述多导联综合指标特征平均值,确定所述心电信号的信号类型。
在上述实施例的基础上,一个或多个处理器60还实现如下操作:
确定所述多导联综合指标特征平均值的历史获取数量;根据所述历史获取数量,确定所述多导联综合指标特征平均值的第四预设阈值;当所述多导联综合指标特征平均值大于所述第四预设阈值时,确定所述信号类型为房室传导阻滞。
在上述实施例的基础上,所述检测器包括卷积神经网络检测器、循环神经网络检测器、卷积循环神经网络检测器、深度神经网络检测器中的至少一种。
实施例五
本申请实施例五还提供一种包含计算机可执行指令的存储介质,所述计算机可执行指令在由计算机处理器执行时用于执行一种房室传导阻滞检测方法,包括:
在心电信号中,提取信号波形;
获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种;
将所述房室传导阻滞特征输入至检测器;
获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果;所述信号类型包括房室传导阻滞和非房室传导阻滞。
当然,本申请实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上所述的房室传导阻滞检测方法操作,还可以执行本申请任意实施例所提供的房室传导阻滞检测方法中的相关操作,且具备相应的功能和有益效果。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(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 (20)

1.一种基于心电信号的房室传导阻滞检测装置,其特征在于,所述装置包括:
信号波形提取模块,用于在心电信号中,提取信号波形;
阻滞特征获取模块,用于获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种或多种;
阻滞特征输入模块,用于将所述房室传导阻滞特征输入至检测器;
输出结果获取模块,用于获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果,所述检测结果包括密度关联检测结果、持续时间检测结果以及变异性检测结果;所述信号类型包括房室传导阻滞和非房室传导阻滞;
所述输出结果获取模块包括:
结果数量统计子模块,用于在所述密度关联检测结果、所述持续时间检测结果和所述变异性检测结果中,统计相匹配结果的结果数量;所述相匹配结果包括房室传导阻滞和非房室传导阻滞;
第一类型确定子模块,用于若所述房室传导阻滞的结果数量大于所述非房室传导阻滞的结果数量,则确定所述房室传导阻滞,作为所述心电信号的信号类型;
第二类型确定子模块,用于若所述房室传导阻滞的结果数量小于所述非房室传导阻滞的结果数量,则确定所述非房室传导阻滞,作为所述心电信号的信号类型。
2.根据权利要求1所述的装置,其特征在于,所述检测器包括密度检测器、时间检测器、变异检测器和综合检测器;所述阻滞特征输入模块包括:
功率谱密度输入子模块,用于将所述功率谱密度关联值输入至所述密度检测器;
持续时间输入子模块,用于将所述PR间期持续时间输入至所述时间检测器;
变异参数输入子模块,用于将所述PR间期变异参数输入至所述变异检测器;
综合参数输入子模块,用于将所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数均输入至所述综合检测器。
3.根据权利要求1所述的装置,其特征在于,所述装置还包括:
心电信号采集模块,用于采集原始心电信号;
低通滤波模块,用于对所述原始心电信号进行低通滤波,得到高频噪声滤除信号,作为所述心电信号;
所述信号波形提取模块包括:
小波变换子模块,用于通过小波变换技术,在所述心电信号中提取P波和QRS波,得到所述信号波形。
4.根据权利要求1所述的装置,其特征在于,所述阻滞特征获取模块包括:
功率谱密度获取子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值时,获取第一功率谱密度,以及,获取第二功率谱密度;所述第一功率谱密度为房室传导阻滞心电信号的功率谱密度;所述第二功率谱密度为待测心电信号的功率谱密度;
密度参数计算子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值时,计算功率谱密度欧范积,以及,计算交叉功率谱密度;所述交叉功率谱密度为所述第一功率谱密度与所述第二功率谱密度的交叉功率谱密度;所述功率谱密度欧范积为所述第一功率谱密度的欧几里得范数与所述第二功率谱密度的欧几里得范数之积;
密度比值计算子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值时,计算所述交叉功率谱密度与所述功率谱密度欧范积的比值,作为所述功率谱密度关联值。
5.根据权利要求4所述的装置,其特征在于,所述功率谱密度获取子模块包括:
谐波分量获取单元,用于获取心电信号的谐波分量差值,以及,获取心电信号的谐波分量功率;所述谐波分量差值为预设频率与谐波分量的差值;所述谐波分量功率为所述谐波分量的功率;
功率谱密度计算单元,用于计算所述谐波分量功率与所述谐波分量差值之和,得到功率谱密度;
第一密度确定单元,用于若所述心电信号为房室传导阻滞心电信号,则确定所述功率谱密度为所述第一功率谱密度;
第二密度确定单元,用于若所述心电信号为待测心电信号,则确定所述功率谱密度为所述第二功率谱密度。
6.根据权利要求1所述的装置,其特征在于,所述检测器包括密度检测器;所述输出结果获取模块包括:
密度关联值提取子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值时,在所述密度检测器中,提取M个导联的功率谱密度关联值,得到多导联功率谱密度关联值;M≥1;
关联平均值计算子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值时,计算所述多导联功率谱密度关联值的平均值,得到多导联功率谱密度关联平均值;
信号类型确定子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值时,根据所述多导联功率谱密度关联平均值,确定所述心电信号的信号类型。
7.根据权利要求6所述的装置,其特征在于,所述信号类型确定子模块包括:
第一数量确定单元,用于确定所述多导联功率谱密度关联平均值的历史获取数量;
第一预设阈值确定单元,用于根据所述历史获取数量,确定所述多导联功率谱密度关联平均值的第一预设阈值;
第一信号类型确定单元,用于当所述多导联功率谱密度关联平均值大于所述第一预设阈值时,确定所述信号类型为房室传导阻滞。
8.根据权利要求1所述的装置,其特征在于,所述阻滞特征获取模块包括:
序列获取子模块,用于当所述房室传导阻滞特征包括所述PR间期持续时间时,获取P波序列,以及,获取R波序列;
持续时间分子计算子模块,用于当所述房室传导阻滞特征包括所述PR间期持续时间时,计算所述R波序列与所述P波序列的差值,得到PR间期持续时间分子;
时序时间计算子模块,用于当所述房室传导阻滞特征包括所述PR间期持续时间时,计算所述PR间期持续时间分子与预设采样频率的比值,作为所述PR间期持续时间。
9.根据权利要求1所述的装置,其特征在于,所述检测器包括时间检测器;所述输出结果获取模块包括:
持续时间提取子模块,用于当所述房室传导阻滞特征包括所述PR间期持续时间时,在所述时间检测器中提取M个导联的PR间期持续时间,得到多导联PR间期持续时间;M≥1;
时间平均值计算子模块,用于当所述房室传导阻滞特征包括所述PR间期持续时间时,计算所述多导联PR间期持续时间的平均值,得到多导联PR间期持续时间平均值;
心电信号类型确定子模块,用于当所述房室传导阻滞特征包括所述PR间期持续时间时,根据所述多导联PR间期持续时间平均值,确定所述心电信号的信号类型。
10.根据权利要求9所述的装置,其特征在于,所述心电信号类型确定子模块包括:
第二数量确定单元,用于确定所述多导联PR间期持续时间平均值的历史获取数量;
第二预设阈值确定单元,用于根据所述历史获取数量,确定所述多导联PR间期持续时间平均值的第二预设阈值;
第二信号类型确定单元,用于当所述多导联PR间期持续时间平均值大于所述第二预设阈值时,确定所述信号类型为房室传导阻滞。
11.根据权利要求1所述的装置,其特征在于,所述阻滞特征获取模块包括:
间期序列获取子模块,用于当所述房室传导阻滞特征包括所述PR间期变异参数时,获取PR间期序列;
相位空间计算子模块,用于当所述房室传导阻滞特征包括所述PR间期变异参数时,计算所述PR间期的相位空间,得到PR间期相位空间;
欧几里得距离计算子模块,用于当所述房室传导阻滞特征包括所述PR间期变异参数时,当所述PR间期相位空间为PR间期延迟相位时,计算所述PR间期延迟相位的欧几里得距离,得到PR间期欧几里得距离;
间期变异参数获取子模块,用于当所述房室传导阻滞特征包括所述PR间期变异参数时,根据所述PR间期欧几里得距离,获取所述PR间期变异参数。
12.根据权利要求11所述的装置,其特征在于,所述间期变异参数获取子模块包括:
分量分子计算单元,用于计算所述PR间期欧几里得距离的阶梯函数,得到PR间期变异分量分子;
变异分量计算单元,用于计算所述PR间期变异分量分子与预设组合参数的比值,得到PR间期变异分量;
分量总和计算单元,用于计算所述PR间期变异分量的多项分量之和,得到PR间期变异分量总和;
变异参数确定单元,用于确定所述PR间期变异分量总和,作为所述PR间期变异参数。
13.根据权利要求1所述的装置,其特征在于,所述检测器包括变异检测器;所述输出结果获取模块包括:
多导联变异参数提取子模块,用于当所述房室传导阻滞特征包括所述PR间期变异参数时,在所述变异检测器中提取M个导联的PR间期变异参数,得到多导联PR间期变异参数;M≥1;
多导联变异平均值计算子模块,用于当所述房室传导阻滞特征包括所述PR间期变异参数时,计算所述多导联PR间期变异参数的平均值,得到多导联PR间期变异平均值;
第一阻滞信号类型确定子模块,用于当所述房室传导阻滞特征包括所述PR间期变异参数时,根据所述多导联PR间期变异平均值,确定所述心电信号的信号类型。
14.根据权利要求13所述的装置,其特征在于,所述阻滞信号类型确定子模块包括:
第三数量确定单元,用于确定所述多导联PR间期变异平均值的历史获取数量;
第三预设阈值单元,用于根据所述历史获取数量,确定所述多导联PR间期变异平均值的第三预设阈值;
第三信号类型单元,用于当所述多导联PR间期变异平均值大于所述第三预设阈值时,确定所述信号类型为房室传导阻滞。
15.根据权利要求1所述的装置,其特征在于,所述阻滞特征获取模块包括:
综合指标分母计算子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,计算所述功率谱密度关联值、所述PR间期持续时间与所述PR间期变异参数之和,得到综合指标分母;
综合指标分子计算子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,计算所述功率谱密度关联值、所述PR间期持续时间与所述PR间期变异参数之积,得到综合指标分子;
综合指标特征计算子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,计算所述综合指标分子与所述综合指标分母的比值,得到综合指标特征;
综合指标特征确定子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,确定所述综合指标特征,作为所述房室传导阻滞特征。
16.根据权利要求1所述的装置,其特征在于,所述检测器包括综合检测器;所述输出结果获取模块包括:
多导联综合指标特征提取子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,在所述综合检测器中提取M个导联的综合指标特征,得到多导联综合指标特征;M≥1;
指标特征平均值计算子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,计算所述多导联综合指标特征的平均值,得到多导联综合指标特征平均值;
第二阻滞信号类型确定子模块,用于当所述房室传导阻滞特征包括所述功率谱密度关联值、所述PR间期持续时间和所述PR间期变异参数时,根据所述多导联综合指标特征平均值,确定所述心电信号的信号类型。
17.根据权利要求16所述的装置,其特征在于,所述第二阻滞信号类型确定子模块包括:
第四数量确定单元,用于确定所述多导联综合指标特征平均值的历史获取数量;
第四预设阈值确定单元,用于根据所述历史获取数量,确定所述多导联综合指标特征平均值的第四预设阈值;
第四信号类型确定单元,用于当所述多导联综合指标特征平均值大于所述第四预设阈值时,确定所述信号类型为房室传导阻滞。
18.根据权利要求1所述的装置,其特征在于,所述检测器包括卷积神经网络检测器、循环神经网络检测器、卷积循环神经网络检测器、深度神经网络检测器中的至少一种。
19.一种电子设备,其特征在于,包括:存储器、一个或多个处理器;
所述存储器,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器执行如下步骤:
在心电信号中,提取信号波形;
获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种或多种;
将所述房室传导阻滞特征输入至检测器;
获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果,所述检测结果包括密度关联检测结果、持续时间检测结果以及变异性检测结果;所述信号类型包括房室传导阻滞和非房室传导阻滞;
所述获取所述检测器输出的检测结果,得到所述心电信号的信号类型的步骤包括:在所述密度关联检测结果、所述持续时间检测结果和所述变异性检测结果中,统计相匹配结果的结果数量;所述相匹配结果包括房室传导阻滞和非房室传导阻滞;若所述房室传导阻滞的结果数量大于所述非房室传导阻滞的结果数量,则确定所述房室传导阻滞,作为所述心电信号的信号类型;若所述房室传导阻滞的结果数量小于所述非房室传导阻滞的结果数量,则确定所述非房室传导阻滞,作为所述心电信号的信号类型。
20.一种包含计算机可执行指令的存储介质,其特征在于,所述计算机可执行指令在由计算机处理器执行时用于执行如下步骤:
在心电信号中,提取信号波形;
获取所述信号波形的房室传导阻滞特征;所述房室传导阻滞特征包括功率谱密度关联值、PR间期持续时间以及PR间期变异参数中的任意一种或多种;
将所述房室传导阻滞特征输入至检测器;
获取所述检测器输出的检测结果,得到所述心电信号的信号类型;所述检测结果为所述检测器根据所述房室传导阻滞特征进行检测的结果,所述检测结果包括密度关联检测结果、持续时间检测结果以及变异性检测结果;所述信号类型包括房室传导阻滞和非房室传导阻滞;
所述获取所述检测器输出的检测结果,得到所述心电信号的信号类型的步骤包括:在所述密度关联检测结果、所述持续时间检测结果和所述变异性检测结果中,统计相匹配结果的结果数量;所述相匹配结果包括房室传导阻滞和非房室传导阻滞;若所述房室传导阻滞的结果数量大于所述非房室传导阻滞的结果数量,则确定所述房室传导阻滞,作为所述心电信号的信号类型;若所述房室传导阻滞的结果数量小于所述非房室传导阻滞的结果数量,则确定所述非房室传导阻滞,作为所述心电信号的信号类型。
CN201910817863.0A 2019-08-30 2019-08-30 基于心电信号的房室传导阻滞检测方法、装置和电子设备 Active CN110522443B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910817863.0A CN110522443B (zh) 2019-08-30 2019-08-30 基于心电信号的房室传导阻滞检测方法、装置和电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910817863.0A CN110522443B (zh) 2019-08-30 2019-08-30 基于心电信号的房室传导阻滞检测方法、装置和电子设备

Publications (2)

Publication Number Publication Date
CN110522443A CN110522443A (zh) 2019-12-03
CN110522443B true CN110522443B (zh) 2022-05-03

Family

ID=68665708

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910817863.0A Active CN110522443B (zh) 2019-08-30 2019-08-30 基于心电信号的房室传导阻滞检测方法、装置和电子设备

Country Status (1)

Country Link
CN (1) CN110522443B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111449646B (zh) * 2020-03-18 2023-05-02 广州视源电子科技股份有限公司 一种心电图识别方法、装置、存储介质及电子设备
CN111449647B (zh) * 2020-03-18 2023-05-02 广州视源电子科技股份有限公司 一种心电图识别方法、装置、存储介质及电子设备
CN114424944B (zh) * 2021-12-31 2024-05-03 纳龙健康科技股份有限公司 一种房室传导阻滞快速识别方法、终端设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109077720A (zh) * 2018-07-05 2018-12-25 广州视源电子科技股份有限公司 信号处理方法、装置、设备和存储介质
CN109077721A (zh) * 2018-07-20 2018-12-25 广州视源电子科技股份有限公司 房颤检测装置及存储介质
CN109770861A (zh) * 2019-03-29 2019-05-21 广州视源电子科技股份有限公司 心电节律模型的训练及检测方法、装置、设备和存储介质
CN109864739A (zh) * 2019-03-19 2019-06-11 苏州哈特智能医疗科技有限公司 一种从10秒心电图中识别二度房室传导阻滞的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120259234A1 (en) * 2009-12-22 2012-10-11 St. Jude Medical Ab Method and device for detecting incipient a-v node malfunction

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109077720A (zh) * 2018-07-05 2018-12-25 广州视源电子科技股份有限公司 信号处理方法、装置、设备和存储介质
CN109077721A (zh) * 2018-07-20 2018-12-25 广州视源电子科技股份有限公司 房颤检测装置及存储介质
CN109864739A (zh) * 2019-03-19 2019-06-11 苏州哈特智能医疗科技有限公司 一种从10秒心电图中识别二度房室传导阻滞的方法
CN109770861A (zh) * 2019-03-29 2019-05-21 广州视源电子科技股份有限公司 心电节律模型的训练及检测方法、装置、设备和存储介质

Also Published As

Publication number Publication date
CN110522443A (zh) 2019-12-03

Similar Documents

Publication Publication Date Title
CN110522442B (zh) 多导联心电异常检测装置、电子设备和存储介质
CN108478209B (zh) 心电信息动态监护方法和动态监护系统
CN110226921B (zh) 心电信号检测分类方法、装置、电子设备和存储介质
CN109770862B (zh) 心电信号分类方法、装置、电子设备和存储介质
CN110522443B (zh) 基于心电信号的房室传导阻滞检测方法、装置和电子设备
EP2688468B1 (en) Apparatus and method for measuring physiological signal quality
CN109411041B (zh) 心电信息处理方法和心电工作站系统
US7300405B2 (en) Analysis of auscultatory sounds using single value decomposition
CN110236573B (zh) 心理压力状态的检测方法及相关装置
CN109411042B (zh) 心电信息处理方法和心电工作站
Andersen et al. A novel approach for automatic detection of Atrial Fibrillation based on Inter Beat Intervals and Support Vector Machine
CN104173043A (zh) 一种适合于移动平台的心电数据分析方法
Oster et al. Open source Java-based ECG analysis software and Android app for atrial fibrillation screening
Kelwade et al. Radial basis function neural network for prediction of cardiac arrhythmias based on heart rate time series
CN103970975A (zh) 心电数据处理方法和系统
CN110693488A (zh) 心电信号处理系统、方法、电子设备及可读存储介质
CN111839494A (zh) 一种心率监测方法及系统
KR102216047B1 (ko) 자동 결합 뉴런 네트워크를 이용한 심장부정맥 분류 방법 및 장치
Liao et al. Accurate ECG R-peak detection for telemedicine
CN115486849A (zh) 心电信号质量评估方法、装置及设备
AU2018287078A1 (en) System and methods for classifying arrhythmia-related heartbeats
CN110179451B (zh) 心电信号质量检测方法、装置、计算机设备和存储介质
CN114246569A (zh) 检测早搏类型的方法、装置及介质
CN111345815B (zh) 一种心电信号中qrs波的检测方法、装置、设备和存储介质
CN110477904B (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