CN113796873B - 穿戴式动态心电信号分类方法及*** - Google Patents
穿戴式动态心电信号分类方法及*** Download PDFInfo
- Publication number
- CN113796873B CN113796873B CN202111247597.6A CN202111247597A CN113796873B CN 113796873 B CN113796873 B CN 113796873B CN 202111247597 A CN202111247597 A CN 202111247597A CN 113796873 B CN113796873 B CN 113796873B
- Authority
- CN
- China
- Prior art keywords
- scattering
- wavelet
- order
- signals
- electrocardiosignals
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000006870 function Effects 0.000 claims abstract description 56
- 239000011159 matrix material Substances 0.000 claims abstract description 45
- 238000013145 classification model Methods 0.000 claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 23
- 238000012549 training Methods 0.000 claims description 42
- 239000006185 dispersion Substances 0.000 claims description 21
- 238000004590 computer program Methods 0.000 claims description 18
- 238000000605 extraction Methods 0.000 claims description 18
- 230000015654 memory Effects 0.000 claims description 18
- 238000013441 quality evaluation Methods 0.000 claims description 18
- 201000010099 disease Diseases 0.000 claims description 14
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 14
- 238000003745 diagnosis Methods 0.000 claims description 12
- 238000002372 labelling Methods 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 10
- 238000001303 quality assessment method Methods 0.000 claims description 9
- 238000001228 spectrum Methods 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 14
- 238000010801 machine learning Methods 0.000 abstract description 5
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 238000013135 deep learning Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 3
- 208000024172 Cardiovascular disease Diseases 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 3
- 230000006403 short-term memory Effects 0.000 description 3
- 238000013519 translation Methods 0.000 description 3
- 230000002457 bidirectional effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000007787 long-term memory Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000003183 myoelectrical effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000011425 standardization method Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6802—Sensor mounted on worn items
-
- 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/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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
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)
- Artificial Intelligence (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Cardiology (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
本发明提供一种穿戴式动态心电信号分类方法及***，属于心电信号分析技术领域，获取待分类的原始心电信号；利用预先构建的小波散射网络提取待分类的心电信号的小波散射系数，利用小波散射系数计算散射特征矩阵，并转化为散射特征图；利用预先训练好的分类模型对散射特征图进行处理，得到心电信号的分类结果。本发明将穿戴式动态信号分为三类，构建了信号分类数据库；提出了导联脱落和纯噪声的判断准则，降低了计算成本；利用Gabor小波函数和尺度算子搭建了小波散射网络，利用小波散射网络对三类信号提取小波散射系数，构建散射特征矩阵，利用机器学习方法LSTM提取散射矩阵特征并进行分类，实现了穿戴式动态心电信号的自动精准分类。
Description
技术领域
本发明涉及心电信号分析技术领域，具体涉及一种穿戴式动态心电信号分类方法及***。
背景技术
对于心血管病的早期预防、连续监测，动态长时间监测环境下，不同运动状态体内外噪声、运动伪迹变化多样，信号噪声严重，且噪声种类和程度不可预期，穿戴式信号质量复杂多变难以控制，易受诸如身体移动、环境噪声、电磁等干扰，从而降低信号质量，产生误报警。海量数据中的冗余信息与居高不下的误警率不仅导致医疗资源浪费，且会造成医生懈怠甚至误诊。因此，动态生理数据必须先进行质量评估，将其分为不可用于和可用于临床的信号，达到去粗留精、去伪存真的目的，从而提高效率和诊断质量。
目前，关于心电信号质量评估主要基于质量评估指标的提取，质量评估指标提取主要涉及到时域、频域，非线性域和空间形态信息等。Newcastle大学Murray课题组提出的频域测度，Langley等人使用直线、基线抬高、基线漂移、过低振幅、过高振幅、峭度等6个波形特征进行质量评估，Zaunseder等人结合功率谱提出35个频域指标。这些方法在实际应用中通常很难合理设定各种阈值参数，导致应用于不同场景和不同导联时泛化能力较弱。
针对心电信号的处理，小波变换利用时频域精准定位特性能够有效的分析非平稳的心电信号，但是无法保持尺度变换的稳定性。
发明内容
本发明的目的在于提供一种基于小波散射网络的穿戴式动态心电信号分类方法及***，以解决上述背景技术中存在的至少一项技术问题。
为了实现上述目的，本发明采取了如下技术方案：
一方面，本发明提供一种穿戴式动态心电信号分类方法，包括：
获取待分类的原始心电信号；
利用预先构建的小波散射网络提取待分类的心电信号的小波散射系数，利用小波散射系数计算散射特征矩阵，并转化为散射特征图；
利用预先训练好的分类模型对散射特征图进行处理，得到心电信号的分类结果；其中，所述预先训练好的分类模型由训练集训练得到，所述训练集包括多个心电信号以及标注心电信号的类别的标签。
优选的，基于尺度函数和小波函数构建了三层小波散射网络，生成0阶散射系数、1阶散射系数和2阶散射系数。
优选的，待分类的心电信号与尺度函数卷积，获得0阶小波散射系数；待分类的心电信号与一阶小波函数卷积，并经过非线性模运算生成一阶散射传播算子；一阶散射传播算子与尺度函数卷积获得一阶小波散射系数；一阶散射传播算子与二阶小波函数卷积，并经过非线性模运算生成二阶散射传播算子；二阶散射传播算子与尺度函数卷积获得二阶小波散射系数。
优选的，基于信号的长度与采样频率，确定尺度函数保持平移不变性的最大尺度因子，即时间支持度。
优选的，训练分类模型的基础网络为长短期记忆神经网络。
优选的，所述训练集来自于质量评估数据库，所述质量评估数据库中包括人工标注完成的多条心电信号，其中，心电信号的类别包括可用于疾病检测和诊断的干净信号、可用于心率提取的轻度污染信号、需要剔除的噪声污染严重信号。
优选的，对质量评估数据库中所有的数据首先进行预处理，包括：数据标准化处理、导联脱落判断以及纯噪声判断。
第二方面，本发明提供一种穿戴式动态心电信号分类***，包括：
获取模块，用于获取待分类的原始心电信号；
提取模块，用于利用预先构建的小波散射网络提取待分类的心电信号的小波散射系数，利用小波散射系数计算散射特征矩阵，并转化为散射特征图；
分类模块，用于利用预先训练好的分类模型对散射特征图进行处理，得到心电信号的分类结果；其中，所述预先训练好的分类模型由训练集训练得到，所述训练集包括多个心电信号以及标注心电信号的类别的标签。
第三方面，本发明提供一种非暂态计算机可读存储介质，所述非暂态计算机可读存储介质用于存储计算机指令，所述计算机指令被处理器执行时，实现如上所述的穿戴式动态心电信号分类方法。
第四方面，本发明提供一种电子设备，包括：处理器、存储器以及计算机程序；其中，处理器与存储器连接，计算机程序被存储在存储器中，当电子设备运行时，所述处理器执行所述存储器存储的计算机程序，以使电子设备执行实现如上所述的穿戴式动态心电信号分类方法的指令。
本发明有益效果：
将穿戴式动态信号分为三类：可用于疾病检测和诊断的比较干净的信号、仅可用于心率提取的轻度污染信号、噪声污染严重的需要剔除的信号；利用Gabor小波函数和尺度算子搭建了小波散射网络，利用小波散射网络提取小波散射系数，利用散射系数构建散射特征矩阵；提出了导联脱落和纯噪声的判断准则，大大降低计算成本；利用小波散射对三类信号提取散射特征矩阵，利用机器学习方法LSTM提取散射矩阵的特征并进行分类，实现了穿戴式动态心电信号的自动精准分类。
本发明附加的方面和优点将在下面的描述中部分给出，这些将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例所述的基于小波散射网络的穿戴式动态心电信号的分类方法流程图。
图2为本发明实施例所述的原始心电信号和添加噪声后的心电信号示意图。
具体实施方式
下面详细叙述本发明的实施方式，所述实施方式的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过附图描述的实施方式是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。
本技术领域技术人员可以理解，除非另外定义，这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。
还应该理解的是，诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义，并且除非像这里一样定义，不会用理想化或过于正式的含义来解释。
本技术领域技术人员可以理解，除非特意声明，这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是，本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件，但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件和/或它们的组。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
为便于理解本发明，下面结合附图以具体实施例对本发明作进一步解释说明，且具体实施例并不构成对本发明实施例的限定。
本领域技术人员应该理解，附图只是实施例的示意图，附图中的部件并不一定是实施本发明所必须的。
实施例1
利用非线性模运算和尺度算子级联小波变换构建具有深度卷积网络的小波散射网络，小波散射变换在时频域具有非常高的分辨率，并且具有计算平移不变性，对于形变能够保持较高的稳定性，可以保留高频信息。另外，小波散射网络在提取数据特征时具有深度学习架构的特性：多尺度收缩、层次对称的线性化和系数表征。小波散射网络保持了传统方法的优势，同时融入了深度学习网络的特性，在图形识别，音频分析，信号处理等领域具有较高的特性。因此，小波散射网络对于混入多种复杂噪声的心电质量评估工作具有非常大的优势。
因此，本实施例1提供一种基于小波散射网络的穿戴式动态心电信号分类***及方法。
首先，本实施例1中提供的穿戴式动态心电信号分类***，包括：
获取模块，用于获取待分类的原始心电信号；
提取模块，用于利用预先构建的小波散射网络提取待分类的心电信号的小波散射系数，利用小波散射系数计算散射特征矩阵，并转化为散射特征图；
分类模块，用于利用预先训练好的分类模型对散射特征图进行处理，得到心电信号的分类结果；其中，所述预先训练好的分类模型由训练集训练得到，所述训练集包括多个心电信号以及标注心电信号的类别的标签。
利用上述的穿戴式动态心电信号分类***，实现了穿戴式动态心电信号分类方法，包括：
获取待分类的原始心电信号；
利用预先构建的小波散射网络提取待分类的心电信号的小波散射系数，利用小波散射系数计算散射特征矩阵，并转化为散射特征图；
利用预先训练好的分类模型对散射特征图进行处理，得到心电信号的分类结果；其中，所述预先训练好的分类模型由训练集训练得到，所述训练集包括多个心电信号以及标注心电信号的类别的标签。
在本实施例1中，基于尺度函数和小波函数构建了三层小波散射网络，生成0阶散射系数、1阶散射系数和2阶散射系数。
具体的，待分类的心电信号与尺度函数卷积，获得0阶小波散射系数；待分类的心电信号与一阶小波函数卷积，并经过非线性模运算生成一阶散射传播算子；一阶散射传播算子与尺度函数卷积获得一阶小波散射系数；一阶散射传播算子与二阶小波函数卷积，并经过非线性模运算生成二阶散射传播算子；二阶散射传播算子与尺度函数卷积获得二阶小波散射系数。
在构建小波散射网络时，基于信号的长度与采样频率，确定尺度函数保持平移不变性的最大尺度因子，即时间支持度。
本实施例1中，训练分类模型所使用的基础网络为长短期记忆神经网络。
其中，所述训练集来自于质量评估数据库，所述质量评估数据库中包括人工标注完成的多条心电信号，其中，心电信号的类别包括可用于疾病检测和诊断的干净信号、可用于心率提取的轻度污染信号、需要剔除的噪声污染严重信号。对长时动态穿戴式动态心电信号利用10s滑动窗进行分割，并通过添加噪声的方式弥补数据不均衡问题。
对质量评估数据库中所有的数据首先进行预处理，包括：数据标准化处理、导联脱落判断以及纯噪声判断。
分割后的10s信号首先进行是否导联脱落的判断，如果为导联脱落信号，直接剔除，程序结束，可以节省计算成本；然后进行纯高斯噪声的判断，如果为纯高斯噪声信号，直接剔除，程序结束，可以节省计算成本。
实施例2
由于身体移动、环境噪声、电磁干扰等，穿戴式动态心电信号质量复杂多变难以控制，为医生诊断带来很大的干扰，如何对穿戴式心电信号质量评估，去粗留精、去伪存真，是提高心血管早期检测和诊断效率的亟待解决的关键问题。本发明实施例2中，提出了一种穿戴式动态心电信号的分类方法，充分利用小波散射网络的优势，提取穿戴式动态心电信号的散射特征，构建散射特征矩阵，利用长短期记忆网络(LSTM，Long Short-Term Memory)对信号进行分类，基于分类结果踢除噪声污染严重的信号，并将保留信号分为仅用于心率提取的轻度污染信号和可用于疾病分类的干净信号，为心血管疾病的实时精准检测提供更有力的数据。
本实施例2所述的穿戴式动态心电信号的分类方法，需要建立干净信号、轻度污染信号和严重污染信号数据库，并需要构建合适的小波散射网络，从而准确获取数据库中各类数据的散射特征矩阵特征，同时需要确定三类数据的质量分类效果的评估指标，测评模型的分类效果。
本实施例2中，利用Gabor小波函数和尺度算子搭建了包含41个1阶小波函数和7个2阶小波函数的小波散射网络，利用小波散射网络提取0阶、1阶和2阶的小波散射系数，利用散射系数构建81×20维的散射特征矩阵，最后利用深度学习LSTM方法提取散射特征图中的特征并分类，将穿戴式动态信号分为三类，一类是需要剔除的噪声污染严重信号，一类是可用于心率提取的轻度污染信号，一类是可用于疾病检测和诊断的比较干净的信号。
小波散射具有平移不变性，特征稳定性以及信息丰富性等特点，这些特点让它对于穿戴式动态心电信号的形变具有较高的辨别性。本实施例2中，构建三层小波散射网络主要包括以下步骤：
选择时间支持度为I(Invariance scale)的尺度函数φI和Morlet小波函数ψ构建三层小波散射网络，生成0阶、1阶和2阶散射系数，该网络可以覆盖信号的整个频域范围。网络搭建步骤如下：
(1)X(t)代表需要分析的心电信号，X(t)与尺度函数φI卷积，获得0阶小波散射系数S0：S0X(t)＝X(t)*φI；
(2)心电信号X(t)与一阶小波函数卷积，并经过非线性模运算生成一阶散射传播算子/>
(3)一阶传播算子与尺度函数卷积获得一阶小波散射系数S1：
(4)一阶散射传播算子与二阶小波函数/>卷积，并经过非线性模运算生成二阶散射传播算子/>
(5)二阶散射传播算子与尺度函数卷积获得二阶小波散射系数S2：
为了实现三类质量数据精准分类，本实施例2中，构建了数据库，并对数据库中的的数据进行了预处理：对长时动态穿戴式动态心电信号利用10s滑动窗进行分割，并通过添加噪声的方式弥补数据不均衡问题。分割后的10s信号首先进行是否导联脱落的判断，如果为导联脱落信号，直接剔除，程序结束，可以节省计算成本；然后进行纯高斯噪声的判断，如果为纯高斯噪声信号，直接剔除，程序结束，可以节省计算成本。
通过以上两种判断的信号输入构建的三层小波散射网络计算小波散射特征矩阵。首先需要确定尺度函数的时间支持度I，即保持平移不变性的最大尺度因子,基于信号的长度与采样频率Fs确定合适的I值，实验证明I一般取值为2*Fs～10*Fs。计算0阶、1阶和2阶散射系数，将三层小波网络的散射系数提取构建散射特征矩阵，并转化为散射特征图。
将训练数据小波散射网络提取的散射特征矩阵输入长短期记忆神经网络LSTM，求解器采用Adam算法，确定合适的梯度阈值、最大轮数、小批量数等参数。序列输入层的维度为81×20，经过双向LSTM层，指定隐藏节点，输出最后一个分类值，然后进入全连接层，指出输出类别，后面是softmax层，输出各个类别分类的概率，最后经过分类层输出最后的分类结果。
本实施例2中，提供的基于小波散射网络的穿戴式动态心电质量评估方法，将穿戴式动态信号分为A、B、C三类，A类是可用于疾病检测和诊断的比较干净的信号，B类是仅可用于心率提取的轻度污染信号，C类是噪声污染严重的需要剔除的信号。利用Gabor小波函数和尺度算子搭建了0阶、1阶和2阶的小波散射网络，利用小波散射网络提取3层小波散射系数，利用散射系数构建散射特征矩阵，最后利用深度学习LSTM方法提取散射矩阵的特征并分类。
目前关于穿戴式动态心电质量评估的研究大部分将数据质量等级划分为两类，本实施例2中将信号分类三类，并构建质量评估数据库，数据库共有数据31111条10s信号，干净信号11709条，轻度污染信号7860条，严重污染信号11542条。其中，提出的导联脱落和纯噪声的判断准则，大大降低计算成本。利用小波散射这一信号分析的利器，对A、B、C三类信号提取散射特征矩阵，然后利用机器学习方法LSTM提取散射矩阵的特征并进行分类，从而实现了穿戴式动态心电数据自动分类。现有数据库中的31111个数据进行训练和测试，正确率在95％左右。
实施例3
本实施例3中，提供了一种基于小波散射网络的穿戴式设备获取的动态心电信号的分类方法，如图1所示，主要步骤，共分为四个步骤。第一步构建穿戴式动态心电质量评估数据库，将根据专家标注的信息将长时的动态心电划分为长度为10s的信号，并确定每条信号的标签；第二步对数据进行预处理，首先是标准化处理，然后剔除导联脱落的信号和纯高斯噪声信号，用于降低计算成本；第三步将剩余的信号输入构建好的小波散射网络，生成散射特征矩阵；第四步利用数据库70％的数据作为训练数据训练LSTM网络，30％的数据作为测试数据，测试模型的分类效果。
本实施例3中，构建穿戴式动态心电质量评估数据库包括：
现有穿戴式动态心电信号18组，每组数据长度约为24小时，数据来自于15个测试人员(男6人，女9人)，年龄分布为21-83岁，利用穿戴式心电设备记录了测试者日常生活环境下心电信号，信号的采样频率为1000Hz。为降低计算成本，本实施例中，将信号降采样频率为250Hz，根据专家标注的质量标注信息将长时的动态信号划分为长度为10s信号，如表1所示，并确定每条信号的质量标签。信号质量等级划分为三类：A类可用于疾病检测和诊断的比较干净的信号11709组，B类为可用以心率提取的轻度污染信号7860组，C类为需要剔除的噪声污染严重信号2687组。
由于噪声数据量较小，将B类信号随机分为6组，添加六种噪声(极限偏移2组，肌电噪声2组，电极移动噪声2组)，另外从A类数据中随机选取1000条数据添加高斯噪声，添加噪声的数据与原始的2687条噪声数据共同构成噪声数据11542条。三类数据共31111条，构成信号质量评估数据库，如图2所示，N1-N6是来自B类轻度污染的6条信号，M1-6是该6条信号添加6种噪声后的变成的C类严重污染信号，N7是来自于A类的干净信号，M7是该条信号添加高斯噪声后变成的严重污染信号。
表1.数据库三类信号
对质量评估数据库中所有的数据首先进行预处理，数据预处理阶段分为三个部分。约定一组ECG信号S＝{s1,s2,...,sn-1,sn}共有n个样本数量，该信号的一阶差分diff(S)＝{s2-s1,s3-s2,...,sn-sn-1}。
数据标准化处理：
为了消除不同种类噪声之间的影响，首先需要进行数据标准化处理，利用min-max标准化方法，对原始数据线性变换映射到[0-1]之间，公式如下：
导联脱落判断：
标准化之后的10s信号进行导联脱落的判断，如果为导联脱落信号直接剔除，程序结束，可以节省计算成本。对于一条信号如果相邻两个样本幅值相等的样本数量大于总样本数量的60％，该信号认为是导联脱落信号，计算公式如下：
对于样本si：
纯噪声判断：
第三步进行纯噪声的判断，如果为纯噪声信号直接剔除，程序结束，可以节省计算成本。基于心电信号频谱范围为0-40Hz，如果信号在0-40Hz范围内的功率谱能量占总能量的比值小于30％，说明该信号主要成分不是心电信号，属于噪声信号，可直接舍弃。计算公式如下：
本实施例3中，构建小波散射网络计算小波散射特征矩阵包括：
通过以上两种判断的信号输入构建的三层小波散射网络计算小波散射特征矩阵。首先需要确定尺度函数的时间支持度I，即保持平移不变性的最大尺度因子,基于信号的长度与采样频率Fs，确定合适的I值，实验证明I一般取值为2*Fs～10*Fs。计算0阶、1阶和2阶散射系数，将三层小波网络的散射系数提取构建散射特征矩阵。网络搭建步骤如下：
(1)X(t)代表需要分析的心电信号，X(t)与尺度函数φI卷积，获得0阶小波散射系数S0：S0X(t)＝X(t)*φI；
(2)心电信号X(t)与一阶小波函数卷积，并经过非线性模运算生成一阶散射传播算子/>
(3)一阶传播算子与尺度函数卷积获得一阶小波散射系数S1：
(4)一阶散射传播算子与二阶小波函数/>卷积，并经过非线性模运算生成二阶散射传播算子/>
(5)二阶散射传播算子与尺度函数卷积获得二阶小波散射系数S2：
本实施例3中，利用机器学习长短期记忆网络LSTM建立分类模型包括：
将训练数据小波散射网络提取的散射特征矩阵输入LSTM，求解器采用Adam算法，确定合适的梯度阈值、最大轮数、小批量数等参数。序列输入层的维度为81×20，经过双向LSTM层，指定隐藏节点，输出最后一个分类值，然后进入全连接层，指出输出类别，后面是softmax层，输出各个类别分类的概率，最后经过分类层输出最后的分类结果。
本实施例3中，数据库总数量为31111，随机挑选70％的数据作为训练数据，剩余30％作为测试数据。训练过程设置参数设定如下，采用Adam作为求解器，maxEpoch为1000，minBatchsize为1000，InitialLength为0.001，执行环境为单GPU。
表2给出测试结果的混淆矩阵。测试训练集共9333组数据，其中A类数据是3513组，B类数据是2352组，C类数据是3495组。
本实施例3中，用的评估指标有敏感度Sensitiviety(Se),精确度Precision(+P),三种质量等级的精确度和敏感度平衡综合性指标F1测度：F1A,F1B,F1C,以及准确率Accuracy(ACC)。
敏感度Se：被真实预测正确的某类样本占训练集中所有某类样本总数的比例，以A类样本为例：
精确度+P：被真实预测正确的某类样本占训练集中所有被预测某类样本总数的比例。以A类样本为例：
表2
精确度和敏感度平衡综合性指标F1测度：
其中TNA、TNB、TNC为分别被准确预测为A、B、C类信号的组数，NA、NB、NC为训练集中被预测为A、B、C类信号的组数，TA、TB、TC为训练集中为A、B、C类信号的组数。
准确率代表被真实预测正确的所有类样本占训练集中所有样本总数的比例，计算公式如下：
最终，三种质量等级信号的分类结果如表3所示。
表3
综上，本实施例3中，所提出的穿戴式心电信号分类方法，构建了一个包含三类质量等级的穿戴式动态心电质量评估数据库；确定导联脱落和纯噪声信号的判断标准；构建了适应于提取穿戴式动态心电散射矩阵的三层小波散射网络；确定适了用于三类质量等级心电信号特征提取的合理散射路径；确定了适用于三类质量等级心电信号分类效果的评估指标；最终利用机器学习方法LSTM提取A、B、C三个质量等级的动态心电散射矩阵的特征，实现穿戴式动态心电数据质量的三分类，分类准确率达到95.44％。剔除污染严重的数据，筛选适用于疾病诊断的干净信号，为医生节省诊断时间。
实施例4
本发明实施例4提供一种非暂态计算机可读存储介质，所述非暂态计算机可读存储介质用于存储计算机指令，所述计算机指令被处理器执行时，实现如上所述的穿戴式动态心电信号分类方法，该方法包括：
获取待分类的原始心电信号；
利用预先构建的小波散射网络提取待分类的心电信号的小波散射系数，利用小波散射系数计算散射特征矩阵，并转化为散射特征图；
利用预先训练好的分类模型对散射特征图进行处理，得到心电信号的分类结果；其中，所述预先训练好的分类模型由训练集训练得到，所述训练集包括多个心电信号以及标注心电信号的类别的标签。
实施例5
本发明实施例5提供一种计算机程序(产品)，包括计算机程序，所述计算机程序当在一个或多个处理器上运行时，用于实现如上所述的穿戴式动态心电信号分类方法，该方法包括：
获取待分类的原始心电信号；
利用预先构建的小波散射网络提取待分类的心电信号的小波散射系数，利用小波散射系数计算散射特征矩阵，并转化为散射特征图；
利用预先训练好的分类模型对散射特征图进行处理，得到心电信号的分类结果；其中，所述预先训练好的分类模型由训练集训练得到，所述训练集包括多个心电信号以及标注心电信号的类别的标签。
实施例6
本发明实施例6提供一种电子设备，包括：处理器、存储器以及计算机程序；其中，处理器与存储器连接，计算机程序被存储在存储器中，当电子设备运行时，所述处理器执行所述存储器存储的计算机程序，以使电子设备执行实现如上所述的穿戴式动态心电信号分类方法的指令，该方法包括：
获取待分类的原始心电信号；
利用预先构建的小波散射网络提取待分类的心电信号的小波散射系数，利用小波散射系数计算散射特征矩阵，并转化为散射特征图；
利用预先训练好的分类模型对散射特征图进行处理，得到心电信号的分类结果；其中，所述预先训练好的分类模型由训练集训练得到，所述训练集包括多个心电信号以及标注心电信号的类别的标签。
综上所述，本发明实施例所述的穿戴式动态心电信号分类方法及***，充分利用小波散射网络的优势，提取了穿戴式动态心电信号的散射特征，构建散射特征矩阵，利用长短期记忆网络(LSTM，Long Short-Term Memory)对信号进行分类，基于分类结果踢除噪声污染严重的信号，并将保留信号分为仅用于心率提取的轻度污染信号和可用于疾病分类的干净信号，为心血管疾病的实时精准检测提供更有力的数据。此方法需要建立干净信号、轻度污染信号和严重污染信号数据库，并需要构建合适的小波散射网络，从而准确获取数据库中各类数据的散射特征矩阵特征，同时需要确定三类数据的质量分类效果的评估指标，测评模型的分类效果。
本领域内的技术人员应明白，本发明的实施例可提供为方法、***、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明公开的技术方案的基础上，本领域技术人员在不需要付出创造性劳动即可做出的各种修改或变形，都应涵盖在本发明的保护范围之内。
Claims (4)
1.一种穿戴式动态心电信号分类方法，其特征在于，包括：
获取待分类的原始心电信号；
利用预先构建的小波散射网络提取待分类的心电信号的小波散射系数，利用小波散射系数构建散射特征矩阵，并转化为散射特征图；其中，基于尺度函数和小波函数构建三层小波散射网络，生成0阶散射系数、1阶散射系数和2阶散射系数；
利用预先训练好的分类模型对散射特征图进行处理，得到心电信号的分类结果，基于分类结果踢除噪声污染严重的信号，并将保留信号分为仅用于心率提取的轻度污染信号和可用于疾病分类的干净信号；其中，所述预先训练好的分类模型由训练集训练得到，所述训练集包括多个心电信号以及标注心电信号的类别的标签，该模型的训练过程包括：
构建训练集，所述训练集来自于质量评估数据库，所述质量评估数据库中包括人工标注完成的多条心电信号，其中，心电信号的类别包括可用于疾病检测和诊断的干净信号、可用于心率提取的轻度污染信号、需要剔除的噪声污染严重信号；对质量评估数据库中所有的数据首先进行预处理，包括：数据标准化处理、导联脱落判断以及纯噪声判断；
具体的，获取设定采样频率下的长时动态穿戴式动态心电信号，利用滑动窗对心电信号进行分割，确定每条限号的质量标签，并通过添加噪声来弥补数据不均衡问题，构成信号质量评估数据库；对心电信号进行导联脱落判断，对于一组ECG信号，若相邻两个样本幅值相等的样本数量大于总样本数量的设定值，则认为该信号是导联脱落信号，剔除该信号；对心电信号进行纯高斯噪声判断，基于心电信号的设定频谱范围，若信号在该设定频谱范围内的功率谱能量占总能量的比值小于预设值，则认为该信号属于噪声信号，直接舍弃；
将通过预处理的信号作为三层小波散射网络的输入信号，计算0阶、1阶和2阶散射系数；其中，基于信号的长度与采样频率Fs，确定尺度函数保持平移不变性的最大尺度因子，即时间支持度I，该I取值为2*Fs～10*Fs；输入的待分类的心电信号与尺度函数卷积，获得0阶小波散射系数；待分类的心电信号与一阶小波函数卷积，并经过非线性模运算生成一阶散射传播算子；一阶散射传播算子与尺度函数卷积获得一阶小波散射系数；一阶散射传播算子与二阶小波函数卷积，并经过非线性模运算生成二阶散射传播算子；二阶散射传播算子与尺度函数卷积获得二阶小波散射系数；
将训练数据小波散射网络提取的散射特征矩阵输入基于长短期记忆网络LSTM的分类模型中，进行训练，得到训练完成的分类模型。
2.一种穿戴式动态心电信号分类***，其特征在于，包括：
获取模块，用于获取待分类的原始心电信号；
提取模块，用于利用预先构建的小波散射网络提取待分类的心电信号的小波散射系数，利用小波散射系数计算散射特征矩阵，并转化为散射特征图；其中，基于尺度函数和小波函数构建三层小波散射网络，生成0阶散射系数、1阶散射系数和2阶散射系数；
分类模块，用于利用预先训练好的分类模型对散射特征图进行处理，得到心电信号的分类结果，基于分类结果踢除噪声污染严重的信号，并将保留信号分为仅用于心率提取的轻度污染信号和可用于疾病分类的干净信号；其中，所述预先训练好的分类模型由训练集训练得到，所述训练集包括多个心电信号以及标注心电信号的类别的标签，该模型的训练过程包括：
构建训练集，所述训练集来自于质量评估数据库，所述质量评估数据库中包括人工标注完成的多条心电信号，其中，心电信号的类别包括可用于疾病检测和诊断的干净信号、可用于心率提取的轻度污染信号、需要剔除的噪声污染严重信号；对质量评估数据库中所有的数据首先进行预处理，包括：数据标准化处理、导联脱落判断以及纯噪声判断；
具体的，获取设定采样频率下的长时动态穿戴式动态心电信号，利用滑动窗对心电信号进行分割，确定每条限号的质量标签，并通过添加噪声来弥补数据不均衡问题，构成信号质量评估数据库；对心电信号进行导联脱落判断，对于一组ECG信号，若相邻两个样本幅值相等的样本数量大于总样本数量的设定值，则认为该信号是导联脱落信号，剔除该信号；对心电信号进行纯高斯噪声判断，基于心电信号的设定频谱范围，若信号在该设定频谱范围内的功率谱能量占总能量的比值小于预设值，则认为该信号属于噪声信号，直接舍弃；
将通过预处理的信号作为三层小波散射网络的输入信号，计算0阶、1阶和2阶散射系数；其中，基于信号的长度与采样频率Fs，确定尺度函数保持平移不变性的最大尺度因子，即时间支持度I，该I取值为2*Fs～10*Fs；输入的待分类的心电信号与尺度函数卷积，获得0阶小波散射系数；待分类的心电信号与一阶小波函数卷积，并经过非线性模运算生成一阶散射传播算子；一阶散射传播算子与尺度函数卷积获得一阶小波散射系数；一阶散射传播算子与二阶小波函数卷积，并经过非线性模运算生成二阶散射传播算子；二阶散射传播算子与尺度函数卷积获得二阶小波散射系数；
将训练数据小波散射网络提取的散射特征矩阵输入基于长短期记忆网络LSTM的分类模型中，进行训练，得到训练完成的分类模型。
3.一种非暂态计算机可读存储介质，其特征在于，所述非暂态计算机可读存储介质用于存储计算机指令，所述计算机指令被处理器执行时，实现如权利要求1所述的穿戴式动态心电信号分类方法。
4.一种电子设备，其特征在于，包括：处理器、存储器以及计算机程序；其中，处理器与存储器连接，计算机程序被存储在存储器中，当电子设备运行时，所述处理器执行所述存储器存储的计算机程序，以使电子设备执行实现如权利要求1所述的穿戴式动态心电信号分类方法的指令。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111247597.6A CN113796873B (zh) | 2021-10-26 | 2021-10-26 | 穿戴式动态心电信号分类方法及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111247597.6A CN113796873B (zh) | 2021-10-26 | 2021-10-26 | 穿戴式动态心电信号分类方法及*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113796873A CN113796873A (zh) | 2021-12-17 |
| CN113796873B true CN113796873B (zh) | 2024-05-03 |
Family
ID=78898296
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111247597.6A Active CN113796873B (zh) | 2021-10-26 | 2021-10-26 | 穿戴式动态心电信号分类方法及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113796873B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115517685A (zh) * | 2022-09-15 | 2022-12-27 | 山东大学 | 基于动态学习和多尺度分析的心电信号检测方法及*** |
| CN116584902B (zh) * | 2023-07-18 | 2024-05-07 | 山东大学 | 基于特征选优及可视化的心音分类装置 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1427333A2 (en) * | 2001-09-21 | 2004-06-16 | Cardiomag Imaging | Nonlinear noise reduction for magnetocardiograms using wavelet transforms |
| WO2013036718A1 (en) * | 2011-09-08 | 2013-03-14 | Isis Innovation Ltd. | Determining acceptability of physiological signals |
| CN109117769A (zh) * | 2018-07-31 | 2019-01-01 | 东南大学 | 一种针对穿戴式心电采集的实时质量评估反馈方法 |
| WO2019075035A1 (en) * | 2017-10-12 | 2019-04-18 | Zoll Medical Corporation | SYSTEM AND METHOD FOR DISTINGUISHING A CARDIAC EVENT OF NOISE IN AN ELECTROCARDIOGRAM (ECG) SIGNAL |
| CN109745033A (zh) * | 2018-12-25 | 2019-05-14 | 东南大学 | 基于时频二维图像和机器学习的动态心电质量评估方法 |
| WO2019100565A1 (zh) * | 2017-11-27 | 2019-05-31 | 乐普（北京）医疗器械股份有限公司 | 基于人工智能自学习的动态心电图分析方法和装置 |
| CN111160139A (zh) * | 2019-12-13 | 2020-05-15 | 中国科学院深圳先进技术研究院 | 心电信号的处理方法、装置及终端设备 |
| CN111956221A (zh) * | 2020-09-07 | 2020-11-20 | 南京医科大学 | 一种基于小波散射因子与lstm神经网络模型的颞叶癫痫分类方法 |
| CN113343805A (zh) * | 2021-05-26 | 2021-09-03 | 南京医科大学 | 一种基于rr间期心电数据和集成学习的房颤节律识别方法 |
| CN113361385A (zh) * | 2021-06-03 | 2021-09-07 | 山东大学 | 心音分类方法及***、可读存储介质、电子设备 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060079773A1 (en) * | 2000-11-28 | 2006-04-13 | Allez Physionix Limited | Systems and methods for making non-invasive physiological assessments by detecting induced acoustic emissions |
| US7022077B2 (en) * | 2000-11-28 | 2006-04-04 | Allez Physionix Ltd. | Systems and methods for making noninvasive assessments of cardiac tissue and parameters |
| CN113727638A (zh) * | 2018-12-26 | 2021-11-30 | 生命解析公司 | 用于自动量化信号质量的方法和*** |
-
2021
- 2021-10-26 CN CN202111247597.6A patent/CN113796873B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1427333A2 (en) * | 2001-09-21 | 2004-06-16 | Cardiomag Imaging | Nonlinear noise reduction for magnetocardiograms using wavelet transforms |
| WO2013036718A1 (en) * | 2011-09-08 | 2013-03-14 | Isis Innovation Ltd. | Determining acceptability of physiological signals |
| WO2019075035A1 (en) * | 2017-10-12 | 2019-04-18 | Zoll Medical Corporation | SYSTEM AND METHOD FOR DISTINGUISHING A CARDIAC EVENT OF NOISE IN AN ELECTROCARDIOGRAM (ECG) SIGNAL |
| WO2019100565A1 (zh) * | 2017-11-27 | 2019-05-31 | 乐普（北京）医疗器械股份有限公司 | 基于人工智能自学习的动态心电图分析方法和装置 |
| CN109117769A (zh) * | 2018-07-31 | 2019-01-01 | 东南大学 | 一种针对穿戴式心电采集的实时质量评估反馈方法 |
| CN109745033A (zh) * | 2018-12-25 | 2019-05-14 | 东南大学 | 基于时频二维图像和机器学习的动态心电质量评估方法 |
| CN111160139A (zh) * | 2019-12-13 | 2020-05-15 | 中国科学院深圳先进技术研究院 | 心电信号的处理方法、装置及终端设备 |
| CN111956221A (zh) * | 2020-09-07 | 2020-11-20 | 南京医科大学 | 一种基于小波散射因子与lstm神经网络模型的颞叶癫痫分类方法 |
| CN113343805A (zh) * | 2021-05-26 | 2021-09-03 | 南京医科大学 | 一种基于rr间期心电数据和集成学习的房颤节律识别方法 |
| CN113361385A (zh) * | 2021-06-03 | 2021-09-07 | 山东大学 | 心音分类方法及***、可读存储介质、电子设备 |
Non-Patent Citations (3)
| Title |
|---|
| 《Wavelet Scattering Transform for ECG Beat Classification》;Liu Z,Yao G,Zhang J,Zeng X;《Comput Math Methods Med》;第7卷;正文1-10页 * |
| Liu Z,Yao G,Zhang J,Zeng X.《Wavelet Scattering Transform for ECG Beat Classification》.《Comput Math Methods Med》.2020,第7卷正文1-10页. * |
| 基于脑电（EEG）的帕金森疾病早期诊断的研究;骆春桃;中国优秀硕士学位论文全文数据库;20200615(第06期);E070-30 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113796873A (zh) | 2021-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Übeyli | Combined neural network model employing wavelet coefficients for EEG signals classification | |
| Sridhar et al. | Brain tumor classification using discrete cosine transform and probabilistic neural network | |
| CN110432898A (zh) | 一种基于非线性动力学特征的癫痫发作脑电信号分类*** | |
| CN113796873B (zh) | 穿戴式动态心电信号分类方法及*** | |
| CN103800011B (zh) | 一种基于功能磁共振成像的脑区效应连接分析*** | |
| CN105956623A (zh) | 一种基于模糊熵的癫痫脑电信号分类方法 | |
| CN113069117A (zh) | 一种基于时间卷积神经网络的脑电情绪识别方法及*** | |
| CN112641451A (zh) | 基于单通道脑电信号多尺度残差网络睡眠分期方法及*** | |
| Haider et al. | Evaluation of p300 based lie detection algorithm | |
| CN113076878A (zh) | 基于注意力机制卷积网络结构的体质辨识方法 | |
| Mouleeshuwarapprabu et al. | Nonlinear vector decomposed neural network based EEG signal feature extraction and detection of seizure | |
| Qin et al. | Deep multi-scale feature fusion convolutional neural network for automatic epilepsy detection using EEG signals | |
| Ravi et al. | Epileptic seizure detection using convolutional neural networks and recurrence plots of EEG signals | |
| Singh et al. | An intelligent method for epilepsy seizure detection based on hybrid nonlinear EEG data features using adaptive signal decomposition methods | |
| Sadiq et al. | Auto-correlation based feature extraction approach for EEG alcoholism identification | |
| Lin et al. | Modified multiscale sample entropy and cross-sample entropy based on horizontal visibility graph | |
| Hasanah et al. | Classification of brain tumor on magnetic resonance imaging using support vector machine | |
| CN108597615A (zh) | 一种轻度认知障碍患者痴呆转化的筛选参考方法 | |
| CN115862897B (zh) | 一种基于临床数据的症候群监测方法及*** | |
| CN114464319B (zh) | 一种基于慢特征分析和深度神经网络的ams易感性评估*** | |
| AlSharabi et al. | EEG-based clinical decision support system for Alzheimer's disorders diagnosis using EMD and deep learning techniques | |
| Satapathy et al. | A deep neural model CNN-LSTM network for automated sleep staging based on a single-channel EEG signal | |
| Islam et al. | Genetic algorithm based optimal feature selection extracted by time-frequency analysis for enhanced sleep disorder diagnosis using eeg signal | |
| CN113558643A (zh) | 基于vmd和nltwsvm的多特征癫痫信号分类方法 | |
| Wasatkar et al. | Case Study Medical Images Analysis and Classification with Data-Centric Approach |
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 |