CN111436926B - 基于统计特征和卷积循环神经网络的房颤信号检测方法 - Google Patents

基于统计特征和卷积循环神经网络的房颤信号检测方法 Download PDF

Info

Publication number
CN111436926B
CN111436926B CN202010257349.9A CN202010257349A CN111436926B CN 111436926 B CN111436926 B CN 111436926B CN 202010257349 A CN202010257349 A CN 202010257349A CN 111436926 B CN111436926 B CN 111436926B
Authority
CN
China
Prior art keywords
atrial fibrillation
neural network
signals
convolution
calculating
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
CN202010257349.9A
Other languages
English (en)
Other versions
CN111436926A (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.)
Shandong Computer Science Center National Super Computing Center in Jinan
Shandong Institute of Artificial Intelligence
Original Assignee
Shandong Computer Science Center National Super Computing Center in Jinan
Shandong Institute of Artificial Intelligence
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 Shandong Computer Science Center National Super Computing Center in Jinan, Shandong Institute of Artificial Intelligence filed Critical Shandong Computer Science Center National Super Computing Center in Jinan
Priority to CN202010257349.9A priority Critical patent/CN111436926B/zh
Priority to NL2025957A priority patent/NL2025957B9/en
Publication of CN111436926A publication Critical patent/CN111436926A/zh
Application granted granted Critical
Publication of CN111436926B publication Critical patent/CN111436926B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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
    • 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/347—Detecting the frequency distribution of signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/353—Detecting P-waves
    • 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/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
    • A61B5/7253—Details of waveform analysis characterised by using transforms
    • A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02405—Determining heart rate variability
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • 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/355—Detecting T-waves
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/366—Detecting abnormal QRS complex, e.g. widening

Landscapes

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

Abstract

一种基于统计特征和卷积循环神经网络的房颤信号检测方法，通过对两种不同类型不同维度的特征进行融合，得到特征集，使用粒子群优化算法训练支持向量机，并使用带权重的支持向量机对心电信号进行分类，将统计特征和卷积循环神经网络结合起来，有效解决了目前房颤信号检测存在的问题，更加全面的概括了房颤信号的特征，提高了房颤限号检测的精确度。

Description

基于统计特征和卷积循环神经网络的房颤信号检测方法
技术领域
本发明涉及心电信号处理技术领域，具体涉及一种基于统计特征和卷积循环神经网络的房颤信号检测方法。
背景技术
心房颤动(Atrial Fibrillation)是最常见的持续性心律失常，房颤信号的早发现、早治疗具有十分重要的意义。房颤的主要表现为两点，即P波缺失和RR间期不规则，基于这两点已有研究者使用统计特征或者深度学习从心电信号中检测房颤，但是现有的基于深度学习的房颤信号检测方法可解释性差并不能挖掘房颤信号检测的潜力，且由于它们的泛化能力很差，需要通过精心挑选的数据集，除此之外由于P波存在定位困难，易受噪声影响等问题，而仅基于RR间期进行房颤信号检测，又会出现其他疾病也存在RR间期不规则的问题，因此，对非精心挑选的房颤数据库进行分类仍是一个较大的挑战。
发明内容
本发明为了克服以上技术的不足，提供了一种分别提取一维统计特征和二维卷积循环神经网络特征并将两者融合以提高房颤信号检测精度的方法。
本发明克服其技术问题所采用的技术方案是：
一种基于统计特征和卷积循环神经网络的房颤信号检测方法，包括如下步骤：
a)对心电信号进行切片，得到相同长度的心电信号s(τ)；
b)使用公式
Figure GDA0002960941210000011
通过短时傅里叶变换的方式计算获得具有二维结构的心电信号的频谱图，式中γ()为窗函数，j为复变量，τ与t均为时间变量，ω为频率；
c)使用Pan-Tompkins算法识别切片后的各心电信号s(τ)中的R波位点，根据心电信号的特异性阈值进行QRS波群检测，基于QRS波和P波、Q波、R波、S波及T波计算统计特征，构建一个特征向量S；
d)使用2维卷积循环神经网络对步骤b)得到频谱图提取特征信息，得到特征向量M；
e)将特性向量S与特征向量M加权后进行特征融合；
f)将融合后的每个心电信号对应的特征向量输入到经过训练的带权重的支持向量机中判断心电信号的类型；
g)使用带有权重的支持向量机对心电信号进行分类。
进一步的，步骤a)具体操作步骤如下：
a-1)按房颤数据库分类标准将心电信号分为正常、房颤、噪声、其他四类信号；
a-2)将正常类的信号从中间截取一段长度为9s的信号；
a-3)将持续时间少于20s的其他类信号从中间截取一段长度为9s的信号，对于持续时间大于20s的其他类信号，将其从中间部分随机截取2个长度为9s的信号；
a-4)通过重叠长度为6s的方式对房颤类的信号截取长度为9s的信号；
a-5)通过重叠长度为8s的方式对噪声类的信号截取长度为9s的信号。
步骤b)中短时傅里叶变换的窗口采用汉宁窗，汉宁窗大小为64个采样点，汉宁窗重叠区域大小为32个采样点。
进一步的，步骤c)构建特征向量S包括如下步骤：
c-1)计算心电信号RR间期标准差与方差，计算心电信号最大RR间期、最小RR间期以及平均RR间期，计算心电信号的pNN50、RMSSD、SDSD；
c-2)计算心电信号样本熵、样本熵系数及变异系数；
c-3)在心电信号R峰确定后得到P波所在位置，将P波平均分为六个片段，计算每个片段的平均值，计算心电信号P波的方差、P波偏度以及P波斜度。
进一步的，步骤d)中2维卷积循环神经网络使用6组且每组含4个卷积层的卷积块，在每组卷积块的所有卷积层在一个常数3×3接受域上应用若干个过滤器，之后使用批处理化和Relu激活操作，每组卷积块的前三个卷积层的过滤器数量设置为32，前三个卷积层的步长为1，每组卷积块的最后一层过滤器数量增加32，最后一层步长为2，将经过卷积块处理输出的特征图通过flatten层一维化处理，一维化处理后的特征图使用两个神经元为200的双向LSTM进行时间上的聚合，聚合后的输出结果为特征向量M。
进一步的，步骤e)的操作步骤如下：
e-1)S∈Rp，M∈Rq，R为向量空间，p与q为特征维数，寻求两组投影方向α∈Rp和β∈Rq，使投影后的矩阵aTS与βTM之间的相关性最大，T为矩阵转置，用遍历加权方式确定投影向量α的权重ws以及投影向量β的权重wm，ws+wm＝1，从0.1至0.9遍历权重ws与wm，通过公式
Figure GDA0002960941210000031
计算正确率A，选择A最大时的权重分配，式中TP为房颤正确分为房颤，FP为余类错误分为房颤；
e-2)通过公式
Figure GDA0002960941210000032
计算wsα和wmβ，式中λ1与λ2为拉格朗日乘子，Sxx为wsα和wmβ所形成的矩阵X的协方差，Syy为wsα和wmβ所形成的矩阵Y的协方差，
Figure GDA0002960941210000033
式中η为可调参数，Sw为类内相关矩阵，Sb为类间相关矩阵；
e-3)将步骤e-2)中求wsα和wmβ的问题即化为求特征值问题，再将求得的的非零特征值以非递增的方式排列，选择非零特征值对应的特征向量作为所求典型相关投影矢量集，之后使用投影矢量集变换后的特征向量进行求和，得到融合后的特征向量集Z。
进一步的，步骤f)中带权重的支持向量机使用二进制粒子群优化算法，粒子的适合度为分类的正确结果，选择FM值作为粒子的评价标准。
本发明的有益效果是：通过对两种不同类型不同维度的特征进行融合，得到特征集，使用粒子群优化算法训练支持向量机，并使用带权重的支持向量机对心电信号进行分类，将统计特征和卷积循环神经网络结合起来，有效解决了目前房颤信号检测存在的问题，更加全面的概括了房颤信号的特征，提高了房颤限号检测的精确度。
附图说明
图1为本发明的方法流程图；
图2维本发明的对频谱图使用卷积循环神经网络进行训练的模型结构图。
具体实施方式
下面结合附图1、附图2对本发明做进一步说明。
一种基于统计特征和卷积循环神经网络的房颤信号检测方法，包括如下步骤：
a)对心电信号进行切片，得到相同长度的心电信号s(τ)；
b)使用公式
Figure GDA0002960941210000041
通过短时傅里叶变换的方式计算获得具有二维结构的心电信号的频谱图，式中γ()为窗函数，j为复变量，τ与t均为时间变量，ω为频率；
c)使用Pan-Tompkins算法识别切片后的各心电信号s(τ)中的R波位点，根据心电信号的特异性阈值进行QRS波群检测，基于QRS波和P波、Q波、R波、S波及T波计算统计特征，构建一个特征向量S；
d)使用2维卷积循环神经网络对步骤b)得到频谱图提取特征信息，得到特征向量M；
e)将特性向量S与特征向量M加权后进行特征融合；
f)将融合后的每个心电信号对应的特征向量输入到经过训练的带权重的支持向量机中判断心电信号的类型；
g)使用带有权重的支持向量机对心电信号进行分类。
通过对两种不同类型不同维度的特征进行融合，得到特征集，使用粒子群优化算法训练支持向量机，并使用带权重的支持向量机对心电信号进行分类，将统计特征和卷积循环神经网络结合起来，有效解决了目前房颤信号检测存在的问题，更加全面的概括了房颤信号的特征，提高了房颤限号检测的精确度。
使用不同受试者的单导联大数据集心电信号，增强了算法的泛化能力，使得从心电信号中检测房颤更加具有现实意义。使用统计特征和卷积循环神经网络这两种方法提取特征，分别考虑了房颤的统计特征的意义和卷积循环神经网络模型的优势，将两者结合起来，更加全面的概括了房颤信号的特征。使用前期特征融合方法，将不同维度不同大小的特征使用加权后的判别典型相关分析特征融合，在减少了数据冗余的同时最大化了不同类别之间的差异，对房颤信号检测精度的提升有重大的意义。
进一步的，步骤a)具体操作步骤如下：
a-1)按房颤数据库分类标准将心电信号分为正常、房颤、噪声、其他四类信号；
a-2)将正常类的信号从中间截取一段长度为9s的信号；
a-3)将持续时间少于20s的其他类信号从中间截取一段长度为9s的信号，对于持续时间大于20s的其他类信号，将其从中间部分随机截取2个长度为9s的信号；
a-4)通过重叠长度为6s的方式对房颤类的信号截取长度为9s的信号；
a-5)通过重叠长度为8s的方式对噪声类的信号截取长度为9s的信号。
进一步的，步骤b)中短时傅里叶变换的窗口采用汉宁窗，汉宁窗大小为64个采样点，汉宁窗重叠区域大小为32个采样点，在转换心电频谱图使用频谱对数转化。
进一步的，步骤c)构建特征向量S包括如下步骤：
c-1)计算心电信号RR间期标准差与方差，计算心电信号最大RR间期、最小RR间期以及平均RR间期，计算心电信号的pNN50、RMSSD、SDSD；
c-2)计算心电信号样本熵、样本熵系数及变异系数；
c-3)在心电信号R峰确定后得到P波所在位置，将P波平均分为六个片段，计算每个片段的平均值，计算心电信号P波的方差、P波偏度以及P波斜度。
最终得到如表1所示的具体特征
Figure GDA0002960941210000051
Figure GDA0002960941210000061
表1
进一步的，如附图2所示，步骤d)中2维卷积循环神经网络使用6组且每组含4个卷积层的卷积块，在每组卷积块的所有卷积层在一个常数3×3接受域上应用若干个过滤器，之后使用批处理化和Relu激活操作，每组卷积块的前三个卷积层的过滤器数量设置为32，前三个卷积层的步长为1，每组卷积块的最后一层过滤器数量增加32，最后一层步长为2，将经过卷积块处理输出的特征图通过flatten层一维化处理，一维化处理后的特征图使用两个神经元为200的双向LSTM进行时间上的聚合，聚合后的输出结果为特征向量M。
进一步的，步骤e)的操作步骤如下：
e-1)S∈Rp，M∈Rq，R为向量空间，p与q为特征维数，寻求两组投影方向α∈Rp和β∈Rq，使投影后的矩阵aTS与βTM之间的相关性最大，T为矩阵转置，用遍历加权方式确定投影向量α的权重ws以及投影向量β的权重wm，ws+wm＝1，从0.1至0.9遍历权重ws与wm，通过公式
Figure GDA0002960941210000062
计算正确率A，选择A最大时的权重分配，式中TP为房颤正确分为房颤，FP为余类错误分为房颤；
e-2)获得权重分配后的投影特征向量使用判别典型相关分析，即使aTS与βTM的相关性最大，考虑同类相关关系最大与不同类相关关系最小，使用类内相关矩阵和类间相关矩阵，使类之间的区分度最大，求解相关系数最大化使用拉格朗日乘子法，其具体公式如下所示：
Figure GDA0002960941210000063
通过上述公式计算wsα和wmβ，式中λ1与λ2为拉格朗日乘子，Sxx为wsα和wmβ所形成的矩阵X的协方差，Syy为wsα和wmβ所形成的矩阵Y的协方差，
Figure GDA0002960941210000071
式中η为可调参数，Sw为类内相关矩阵，Sb为类间相关矩阵；
e-3)将步骤e-2)中求wsα和wmβ的问题即化为求特征值问题，再将求得的的非零特征值以非递增的方式排列，选择非零特征值对应的特征向量作为所求典型相关投影矢量集，之后使用投影矢量集变换后的特征向量进行求和，得到融合后的特征向量集Z。
在SVM分类器部分，指定每个分类的权重，以避免训练决策过于偏向样本容量大的分类的现象。设定权值的原则是使样本权值对小样本的响应(噪声)越高越好，使大样本对应(正常)的样本权值越低越好。根据每个类别的大小，我们为每个类的心电图样本设定一个粗略的权重范围，因此进一步的，步骤f)中带权重的支持向量机使用二进制粒子群优化算法(BPSO)，在BPSO算法中，粒子的适合度为分类的正确结果，由于心电数据集不是均衡分布的，一次你准确性或错误率并不代表模型在不均衡数据下的性能，因此选择FM值作为粒子的评价标准，它同时考虑了精度和召回率。

Claims (6)

1.一种基于统计特征和卷积循环神经网络的房颤信号检测方法，其特征在于，包括如下步骤：
a)对心电信号进行切片，得到相同长度的心电信号s(τ)；
b)使用公式
Figure FDA0002960941200000011
通过短时傅里叶变换的方式计算获得具有二维结构的心电信号的频谱图，式中γ()为窗函数，j为复变量，τ与t均为时间变量，ω为频率；
c)使用Pan-Tompkins算法识别切片后的各心电信号s(τ)中的R波位点，根据心电信号的特异性阈值进行QRS波群检测，基于QRS波和P波、Q波、R波、S波及T波计算统计特征，构建一个特征向量S；
d)使用2维卷积循环神经网络对步骤b)得到频谱图提取特征信息，得到特征向量M；
e)将特性向量S与特征向量M加权后进行特征融合；
f)将融合后的每个心电信号对应的特征向量输入到经过训练的带权重的支持向量机中判断心电信号的类型；
g)使用带有权重的支持向量机对心电信号进行分类；
步骤e)包括如下步骤：
e-1)S∈Rp，M∈Rq，R为向量空间，p与q为特征维数，寻求两组投影方向α∈Rp和β∈Rq，使投影后的矩阵aTS与βTM之间的相关性最大，T为矩阵转置，用遍历加权方式确定投影向量α的权重ws以及投影向量β的权重wm，ws+wm＝1，从0.1至0.9遍历权重ws与wm，通过公式
Figure FDA0002960941200000012
计算正确率A，选择A最大时的权重分配，式中TP为房颤正确分为房颤，FP为余类错误分为房颤；
e-2)通过公式
Figure FDA0002960941200000021
计算wsα和wmβ，式中λ1与λ2为拉格朗日乘子，Sxx为wsα和wmβ所形成的矩阵X的协方差，Syy为wsα和wmβ所形成的矩阵Y的协方差，
Figure FDA0002960941200000022
式中η为可调参数，Sw为类内相关矩阵，Sb为类间相关矩阵；e-3)将步骤e-2)中求wsα和wmβ的问题即化为求特征值问题，再将求得的的非零特征值以非递增的方式排列，选择非零特征值对应的特征向量作为所求典型相关投影矢量集，之后使用投影矢量集变换后的特征向量进行求和，得到融合后的特征向量集Z。
2.根据权利要求1所述的基于统计特征和卷积循环神经网络的房颤信号检测方法，其特征在于，步骤a)具体操作步骤如下：
a-1)使用PhysioNet/Cinc Challenge 2017数据库将心电信号分为正常、房颤、噪声、其他四类信号；
a-2)将正常类的信号从中间截取一段长度为9s的信号；
a-3)将持续时间少于20s的其他类信号从中间截取一段长度为9s的信号，对于持续时间大于20s的其他类信号，将其从中间部分随机截取2个长度为9s的信号；
a-4)通过重叠长度为6s的方式对房颤类的信号截取长度为9s的信号；
a-5)通过重叠长度为8s的方式对噪声类的信号截取长度为9s的信号。
3.根据权利要求1所述的基于统计特征和卷积循环神经网络的房颤信号检测方法，其特征在于，步骤b)中短时傅里叶变换的窗口采用汉宁窗，汉宁窗大小为64个采样点，汉宁窗重叠区域大小为32个采样点。
4.根据权利要求1所述的基于统计特征和卷积循环神经网络的房颤信号检测方法，其特征在于，步骤c)构建特征向量S包括如下步骤：
c-1)计算心电信号RR间期标准差与方差，计算心电信号最大RR间期、最小RR间期以及平均RR间期，计算心电信号的pNN50、RMSSD、SDSD；
c-2)计算心电信号样本熵、样本熵系数及变异系数；
c-3)在心电信号R峰确定后得到P波所在位置，将P波平均分为六个片段，计算每个片段的平均值，计算心电信号P波的方差、P波偏度以及P波斜度。
5.根据权利要求1所述的基于统计特征和卷积循环神经网络的房颤信号检测方法，其特征在于步骤d)中2维卷积循环神经网络使用6组且每组含4个卷积层的卷积块，在每组卷积块的所有卷积层在一个常数3×3接受域上应用若干个过滤器，之后使用批处理化和Relu激活操作，每组卷积块的前三个卷积层的过滤器数量设置为32，前三个卷积层的步长为1，每组卷积块的最后一层过滤器数量增加32，最后一层步长为2，将经过卷积块处理输出的特征图通过flatten层一维化处理，一维化处理后的特征图使用两个神经元为200的双向LSTM进行时间上的聚合，聚合后的输出结果为特征向量M。
6.根据权利要求1所述的基于统计特征和卷积循环神经网络的房颤信号检测方法，其特征在于:步骤f)中带权重的支持向量机使用二进制粒子群优化算法，粒子的适合度为分类的正确结果，选择FM值作为粒子的评价标准。
CN202010257349.9A 2020-04-03 2020-04-03 基于统计特征和卷积循环神经网络的房颤信号检测方法 Active CN111436926B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010257349.9A CN111436926B (zh) 2020-04-03 2020-04-03 基于统计特征和卷积循环神经网络的房颤信号检测方法
NL2025957A NL2025957B9 (en) 2020-04-03 2020-06-30 Atrial fibrillation signal detection method based on statistical features and a convolutional recurrent neural network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010257349.9A CN111436926B (zh) 2020-04-03 2020-04-03 基于统计特征和卷积循环神经网络的房颤信号检测方法

Publications (2)

Publication Number Publication Date
CN111436926A CN111436926A (zh) 2020-07-24
CN111436926B true CN111436926B (zh) 2021-04-20

Family

ID=71649854

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010257349.9A Active CN111436926B (zh) 2020-04-03 2020-04-03 基于统计特征和卷积循环神经网络的房颤信号检测方法

Country Status (2)

Country Link
CN (1) CN111436926B (zh)
NL (1) NL2025957B9 (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112022149B (zh) * 2020-09-04 2022-10-04 无锡博智芯科技有限公司 一种基于心电信号的房颤检测方法
CN111920405A (zh) * 2020-09-15 2020-11-13 齐鲁工业大学 一种房颤信号识别系统及方法
CN112244861B (zh) * 2020-10-09 2021-08-10 广东工业大学 一种单导联心电信号f波提取方法
CN112971800B (zh) * 2021-01-27 2023-01-03 东南大学 基于深度学习的心电早搏信号处理方法、模块和装置
CN112842348B (zh) * 2021-02-07 2021-09-14 山东省人工智能研究院 一种基于特征提取与深度学习的心电信号自动分类方法
CN113261973B (zh) * 2021-04-13 2022-09-02 山东师范大学 加权多尺度有限穿越可视图辨识房颤心电信号方法及系统
CN113349753A (zh) * 2021-07-19 2021-09-07 成都芯跳医疗科技有限责任公司 一种基于便携式动态心电监护仪的心律失常检测方法
CN114764581B (zh) * 2022-06-16 2022-09-09 合肥心之声健康科技有限公司 一种基于rr间期空间特征的房颤分类方法、装置及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105686801A (zh) * 2015-12-25 2016-06-22 四川大学 基于固有时间尺度分解的房颤自发终止预测方法
CN106066933A (zh) * 2016-05-27 2016-11-02 成都信汇聚源科技有限公司 获取房颤识别人工神经网络权重值矩阵的方法
CN107203692A (zh) * 2017-05-09 2017-09-26 哈尔滨工业大学（威海） 基于深度卷积神经网络的房颤检测的实现方法
CN108968951A (zh) * 2018-08-15 2018-12-11 武汉中旗生物医疗电子有限公司 心电图检测方法、装置及系统
CN109124620A (zh) * 2018-06-07 2019-01-04 深圳市太空科技南方研究院 一种房颤检测方法、装置及设备
CN109412900A (zh) * 2018-12-04 2019-03-01 腾讯科技（深圳）有限公司 一种网络状态识别的方法、模型训练的方法及装置
CN110840402A (zh) * 2019-11-19 2020-02-28 山东大学 一种基于机器学习的房颤信号识别方法及系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101908138B (zh) * 2010-06-30 2012-05-23 北京航空航天大学 基于噪声独立成分分析的合成孔径雷达图像目标识别方法
CN105320969A (zh) * 2015-11-20 2016-02-10 北京理工大学 基于多尺度Renyi熵的心率变异性特征分类方法
CN108852308B (zh) * 2018-06-22 2021-03-30 苏州高新区人民医院 一种基于多参数影像技术优化检测体系的超时间窗ais临床再灌注方法
CN109359578A (zh) * 2018-10-09 2019-02-19 四川师范大学 加权融合三通道特征步态表征方法
CN110289081B (zh) * 2019-05-14 2021-11-02 杭州电子科技大学 深度网络堆栈模型自适应加权特征融合的癫痫检测方法
CN110664395B (zh) * 2019-09-29 2022-04-12 京东方科技集团股份有限公司 图像处理方法、图像处理装置以及存储介质

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105686801A (zh) * 2015-12-25 2016-06-22 四川大学 基于固有时间尺度分解的房颤自发终止预测方法
CN106066933A (zh) * 2016-05-27 2016-11-02 成都信汇聚源科技有限公司 获取房颤识别人工神经网络权重值矩阵的方法
CN107203692A (zh) * 2017-05-09 2017-09-26 哈尔滨工业大学（威海） 基于深度卷积神经网络的房颤检测的实现方法
CN109124620A (zh) * 2018-06-07 2019-01-04 深圳市太空科技南方研究院 一种房颤检测方法、装置及设备
CN108968951A (zh) * 2018-08-15 2018-12-11 武汉中旗生物医疗电子有限公司 心电图检测方法、装置及系统
CN109412900A (zh) * 2018-12-04 2019-03-01 腾讯科技（深圳）有限公司 一种网络状态识别的方法、模型训练的方法及装置
CN110840402A (zh) * 2019-11-19 2020-02-28 山东大学 一种基于机器学习的房颤信号识别方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Classification of Atrial Fibrillation with Pre-Trained Convolutional Neural Network Models;Abdul Qayyum et al;《2018 IEEE-EMBS Conference on Biomedical Engineering and Sciences》;20181231;全文 *
基于多特征融合与卷积神经网络的房颤检测;魏晓玲等;《激光医学》;20170517;第38卷(第5期);全文 *

Also Published As

Publication number Publication date
NL2025957B9 (en) 2021-05-12
NL2025957B1 (en) 2021-04-20
CN111436926A (zh) 2020-07-24

Similar Documents

Publication Publication Date Title
CN111436926B (zh) 基于统计特征和卷积循环神经网络的房颤信号检测方法
CN105411565A (zh) 基于广义尺度小波熵的心率变异性特征分类方法
CN105760889A (zh) 一种高效的不均衡数据集分类方法
CN105320969A (zh) 基于多尺度Renyi熵的心率变异性特征分类方法
US8010296B2 (en) Apparatus and method for removing non-discriminatory indices of an indexed dataset
CN108846346B (zh) 基于密度峰值聚类算法的神经元峰电位在线分类方法
CN109829412B (zh) 基于动态模式分解分形特征的局部放电模式识别方法
CN112819059A (zh) 一种基于流行保持迁移学习的滚动轴承故障诊断方法
CN109875548B (zh) 一种基于多导联综合分析的心电图特征波形聚类方法
CN112070008A (zh) 高光谱图像特征识别方法、装置、设备及存储介质
CN110292377A (zh) 基于瞬时频率和功率谱熵融合特征的脑电信号分析方法
CN114237046A (zh) 一种基于sift数据特征提取算法与bp神经网络模型的局部放电模式识别方法
CN108985462B (zh) 基于互信息和分形维数的无监督特征选择方法
CN112883895B (zh) 基于自适应加权pca的非法电磁信号检测方法及其实现系统
CN110458071A (zh) 一种基于dwt-dfpa-gbdt的光纤振动信号特征提取与分类方法
CN116584902B (zh) 基于特征选优及可视化的心音分类装置
CN111601358B (zh) 一种多阶段分层分簇空间相关性温度感知数据去冗余方法
CN112287036A (zh) 一种基于谱聚类的离群点检测方法
CN111666999A (zh) 一种遥感图像的分类方法
CN116776245A (zh) 一种基于机器学习的三相逆变器设备故障诊断方法
CN113298138B (zh) 一种雷达辐射源个体识别方法及系统
De Chazal et al. Selection of optimal parameters for ECG diagnostic classification
Köhler et al. Unsupervised feature selection for pattern search in seismic time series
CN115363598B (zh) 一种心电信号处理方法及系统
Hao et al. Omics feature learning for cross individual ALS disease identification with EMG signal

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