CN109934089B - 基于监督梯度提升器的多级癫痫脑电信号自动识别方法 - Google Patents
基于监督梯度提升器的多级癫痫脑电信号自动识别方法 Download PDFInfo
- Publication number
- CN109934089B CN109934089B CN201910034509.0A CN201910034509A CN109934089B CN 109934089 B CN109934089 B CN 109934089B CN 201910034509 A CN201910034509 A CN 201910034509A CN 109934089 B CN109934089 B CN 109934089B
- Authority
- CN
- China
- Prior art keywords
- classifier
- electroencephalogram
- frequency domain
- training
- electroencephalogram signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于监督梯度提升器的多级癫痫脑电信号自动识别方法,属于神经科学、智能信息处理、应用数学、人因工程技术领域。本发明首先对采集的原始脑电信号进行预处理;然后对预处理后的无噪声的脑电信号进行频域特征提取;构建分类癫痫脑电信号的分类器;通过多种评估指标来评估三个分类器的分类性能,确定最适合进行分类多种类别数据的分类器;最后将所确定的分类器的识别结果为后续使用提供决策支持。本发明不仅通过小波变换提取脑电信号的频域特征,而且通过统计技术提取统计特征;同时利用的梯度提升树分类器具有高度并行化,可提高运行效率。
Description
技术领域
本发明属于神经科学、智能信息处理、应用数学、人因工程技术领域,具体涉及一种基于监督梯度提升器的多级癫痫脑电信号自动识别方法。
背景技术
癫痫是最常见的神经系统疾病,全世界每百人中有一人患有癫痫。癫痫发作是大脑中阵发性异常超同步电活动的临床表现,其特征在于迭代、突然和暂时性。然而,脑紊乱的时间-频率是不确定的,并且它们的临床表现不容易被发现。在癫痫患者测试的早期阶段,研究人员尝试使用传感器设备收集患者表面生物学数据,包括心电图(ECG),肌电图(EMG),运动数据监测癫痫患者的数据。这些数据被收集到可穿戴传感器系统中。这些可穿戴传感器系统可以长时间非侵入地监测癫痫患者的生物数据信号。然而,这些生物数据信号具有空间分辨率不足的缺点。
因此,研究人员的目标是直接从大脑获取癫痫信息,如正电子发射断层扫描(PET),单光子发射计算机断层扫描(SPECT),磁共振成像(MRI),功能磁共振成像(fMRI)。如今,一些研究人员应用了视频脑电信号(EEG),它不仅显示时间信息,还为患者提供空间信息,作为研究癫痫的黄金标准和工具。脑电图的反映情节信息不能被其他生理脑功能方法实时替代。由于EEG信号产生的生理过程被认为是非线性、非静止和个性化信号。因此,这些特点给脑电信号分析带来了巨大挑战。
在大数据时代,机器学习技术作为脑电信号分析中一项十分重要的手段是受到高度重视的,医学和工程结合是科学发展的必然趋势,本发明的目标利用机器学习技术找到一种强适用性的、鲁棒的、分类准确率高的方法。
发明内容
本发明针对癫痫病患的脑电信号监控和检测问题,提供一种基于监督梯度提升器的多级癫痫脑电信号自动识别方法,具体涉及对癫痫患者脑电信号自动检测过程中所用到的关键特征提取和分类方法,解决人工处理分析和脑电信号的复杂弊端,有助于辅助医生进行医疗诊断,提升病患和护理者的生活质量,具体是基于计算机的癫痫患者脑电信号处理方案,本发明的目的通过以下技术方案来实现:
基于监督梯度提升器的多级癫痫脑电信号自动识别方法,具体步骤如下:
步骤一:对采集的原始脑电信号进行预处理。
采集的原始脑电信号包括三种类别,分布为连续癫痫发作脑电信号{S}、间歇癫痫发作脑电信号{F/N}和健康被试脑电信号{Z/O}。所述的预处理包括对原始脑电信号进行降采样、滤波、重参考电极、去除眼电噪声、基线校正和ICA分析的处理,最终得到无噪声的脑电信号。
对所述的无噪声的脑电信号进行存储,存储为文本格式的数据。
步骤二:对预处理后的无噪声的脑电信号进行频域特征提取,具体步骤如下;
步骤201、运用小波变换的方法进行频域分析,通过对步骤一得到无噪声的脑电信号做Symlet小波分解,提取脑电信号的五个频域特征,提取的频域特征包括alpha波、beta波、delta波、theta波和gamma波五个波段的频率子带。
步骤202、通过主成分分析法(PCA算法)进行降低维度的处理。
步骤203、对降维后的低维特征进行批量归一化处理,统一维度和格式。
步骤204、分别计算统一维度后的五个频率子带的统计特征,所述统计特征包括能量均值及能量标准方差,为导入分类器做准备。
步骤三、构建分类癫痫脑电信号的分类器,所述分类器包括梯度提升树分类器、支持向量机分类器和随机森林分类器。
步骤301、将原始脑电信号分为训练集和测试集,并以随机的方式将训练集分为十个子集,取其中一个子集作为验证集,其他九个子集作为真正用于训练的训练集。
步骤302、向每个分类器(包括梯度提升树分类器、支持向量机分类器、随机森林分类器)导入脑电数据,准备进行分类器训练。
步骤303、在分类器训练过程中采用十折交叉验证的方式进行分类器训练。
步骤304、采用每次的训练误差的平均值作为最终的交叉验证误差。
步骤305、采用网格搜索优化器对梯度提升树分类器(GBM)、支持向量机分类器(SVM)和随机森林分类器(RF)三种分类器在训练过程中产生的分类模型的参数进行寻优,运用变步长的策略对网格搜索优化器进行改造,然后将网格搜索优化器应用于分类器做优化处理。
步骤306、利用验证集在分类器上进行验证,确定最终的三个分类器。
在步骤一和步骤三之间对所处理的数据进行打标签处理。
步骤四、通过多种分类器评估指标来评估三个分类器的分类性能,确定最适合进行分类多种类别数据的分类器。
步骤401、利用测试集对步骤三中所确定的三个分类器进行测试。
步骤402、定义多个评估指标来评估分类器的性能,评估指标包括:混淆矩阵、受试者工作特征曲线、区域下面积和精确率-召回率曲线。
步骤403、对每个分类器进行性能指标对比分析,最终确定取得所需的识别结果的分类器。
步骤五、将所确定的分类器的识别结果为后续使用提供决策支持。识别结果包括三种类别:安全、低级预警和高级预警。
与现有技术相比,本发明的有益效果在于:
1、不仅通过小波变换提取脑电信号的频域特征,而且通过统计技术提取统计特征。频域特征和统计特征的实现被视为融合特征,其反映了特征的多样性并且在癫痫发作检测中表现出良好的性能。同时,主成分分析法被用于自动减少特征维度以降低硬件复杂性,同时保持高分类精度。
2、梯度提升树分类器具有高度并行化,可提高运行效率。这对于大量脑电数据的训练是有利的。然而,分类器在训练过程中产生大量超参数,依靠手动调整参数难以确定最佳参数。本文提出了网格搜索优化器来优化参数并通过以可变步长的方式重复过滤参数来确定最终的最佳分类器。为了防止梯度提升树分类器训练过程中的过度拟合,本发明使用十折交叉验证的方法,这种设计方案可以保证优化后的分类器更具有鲁棒性。
3、对比三种不同的分类器对癫痫脑电进行分类测试,得出设计的梯度提升树分类器具有最好的分类效果。它可以预测不同程度癫痫患者的病情(连续癫痫发作状态、间歇癫痫发作状态或健康状态)。同时,设计的检测方案由多个指标各方面进行评估和验证。检测策略的研究不仅限于识别的准确性,而且还包括更彻底,更清晰的分析和错误率问题。这种策略在医学筛查中具有重要意义。
附图说明
图1是本发明所提供的自动识别方法流程图;
图2是本发明所提供的Symlets小波分解的过程示意图;
图3A是本发明基于脑电信号中{S}数据做Symlet小波分解图;
图3B是本发明基于脑电信号中{F/N}数据做Symlet小波分解图;
图3C是本发明基于脑电信号中{Z/O}数据做Symlet小波分解图;
图4是本发明基于验证不同分类器产生的混淆矩阵对比图;
图5是本发明基于评估分类器产生的ROC和AUC对比图;
图6是本发明基于评估分类器产生的PR对比图。
具体实施方式
下面结合附图对本发明的具体实施方法进行详细说明。
本发明设计一种基于监督梯度提升器的多级癫痫脑电信号自动识别方法,具体步骤流程如图1所示,包括以下几个步骤:
步骤一:对采集的原始脑电信号进行预处理(脑电预处理),以德国波恩大学的癫痫脑电数据集作为原始脑电信号为实例。
波恩大学的原始脑电信号包括三种类别,具体为连续癫痫发作脑电信号{S}、间歇癫痫发作脑电信号{F/N}和健康被试脑电信号{Z/O}。
具体步骤如下:
通过MATLAB的开源工具箱EEGlab对采集来的原始脑电信号进行预处理。所述的预处理包括对原始脑电信号进行降采样、滤波、重参考电极、去除眼电噪声、基线校正和ICA(独立成分分析)分析等的处理,最终尽可能得到纯净的无噪声的脑电信号。
对所述的无噪声的脑电信号进行存储,存储为txt格式的数据。
步骤二:对预处理后的无噪声的脑电信号进行频域特征提取,该技术为脑电处理的第一项核心技术(脑电特征提取)。
特征类型主要一般分为四个类别,包括统计特征,分形维数特征,熵特征和时频域特征。一些研究利用频域特征提取解决了癫痫发作时非静止脑电信号的自动分类问题。用的最多的是采取小波变换的方式提取脑电信号的频域特征。然而在对脑电信号进行小波变换的过程中采用哪种小波进行分解至关重要。目前,在这个领域中使用最广泛的是采用Daubechies小波进行分解。但本方法调研和对比了Daubechies、Symlets、Haar、Morlet、Mexican Hat、Meyer小波特性,这些特性主要包括:正交性、紧支撑性、支持长度和对称性,如下表1所示:
表1:各小波特性示意表
对比发现Symlet小波在所有小波中具有最优的性质,首先,Symlet小波具备正交性、紧支撑性的同时还具有对称性。其次,Symlet小波是Daubechies小波的改进,弥补了Daubechies小波近似不对称的缺点,Symlet小波的支持范围和消失矩分别为2N-1和N,Symlet小波具有比Daubechies更好的规律性,可以在一定程度上减少信号分析和重构的相位失真,选择Symlet小波,它更适合于脑电信号的时频分析。
步骤201、运用小波变换的方法进行频域分析,通过对步骤一得到无噪声的脑电信号做Symlet小波分解,提取脑电信号的五个频域特征,提取的频域特征包括alpha波、beta波、delta波、theta波和gamma波五个波段的频率子带;具体为通过四级Symlet小波将{S},{F/N},{Z/O}中的数据分解为包括alpha波α(6~12Hz)、beta波β(12~25Hz)、delta波δ(0~3Hz)、theta波θ(3~6Hz)和gamma波γ(25~50Hz)的五个频率子带,本发明Symlet小波分解的过程如图2所示,其中,d1,d2,d3,d4和a1,a2,a3,a4分别为每层分解的细节系数和近似系数。{S}数据、{F/N}数据和{Z/O}数据的四级Symlet小波分解结果如图3A、图3B和图3C所示。
步骤202、由于提取的五个频率子带的特征维度过高,为减小计算机的负载,提高计算机的运行速度,希望减小频率子带特征维度降低运算的复杂性。因此,本发明对步骤201提取的五个高维的频率子带的特征通过主成分分析法(PCA算法)进行降低维度的处理,从而降低方法的开销,为输入自动分类器做准备。
主成分分析法(PCA算法)的主要思想是,通过少数不相关的几个主要的成分表示多个相关的变量的内部结构。分别对脑电信号的五个频率子带做PCA处理,将频率子带的原始脑电信号所在的坐标空间转换为新的坐标所在空间,并且该新坐标空间的选择由五个频率子带本身确定,因为频率子带的最大方差提供了脑电频率特征的重要信息,转换坐标空间时,首先计算每个子带的方差,方差最大的方向用作新坐标空间主方向。第一个新坐标空间的坐标轴选择频率子带数据中方差最大的方向,第二个新坐标空间的坐标轴选择第一个新坐标系的坐标轴和第二个最大方差方向的正交方向作为新坐标轴该过程重复多次,将4097维的频率子带特征降到64维。主成分分析法(PCA算法)的伪代码如下表2所示:
表2:主成分分析法伪代码
其中,rij为s(ti)与s(tj)的协方差,rji为s(tj)与s(ti)的协方差,rii为s(ti)与s(ti)的协方差,s(ti)是第i个脑电信号的频域特征,s(tj)为第j个脑电信号的频域特征,i=1,2,…,p,j=1,2,…,p,i*=1,2,…,m,Vars(ti)为s(ti)的方差,Vars(tj)为s(tj)的方差。
步骤203、对降维后的低维特征进行批量归一化处理,统一维度和格式。
步骤204、分别计算统一维度后的五个频率子带alpha波α(6~12Hz)、beta波β(12~25Hz)、delta波δ(0~3Hz)、theta波θ(3~6Hz)和gamma波γ(25~50Hz)的统计特征,所述统计特征包括能量均值ω*及能量标准方差δ2,为导入分类器做准备;为每个频段制造2个特征向量,一共10的特征向量,均值ω*和标准偏差δ2的定义如下,其中P(s(ti))为频域特征s(ti)的功率;
步骤三、构建分类癫痫脑电信号的分类器(脑电特征分类),所述分类器包括梯度提升树分类器、支持向量机分类器和随机森林分类器。
目前,机器学习的分类器有很多,但是如何选择一个适合数据处理的分类器是一个难以解决的问题,经过运用上述三种不同分类器做实验进行分类结果对比。
具体步骤如下:
步骤301、将原始脑电信号分为训练集和测试集,并以随机的方式将训练集分为十个子集,取其中一个子集作为验证集,其他九个子集作为真正用于训练的训练集,这样可以防止生成的分类模型过拟合。
步骤302、向每个分类器(梯度提升树分类器、支持向量机分类器、随机森林分类器)导入脑电数据。
在分类器中,支持向量机分类器SVM是研究脑电信号最常用的分类器,因为它具有良好的分类性能和出色的泛化能力。该分类器的本质在于在两种类型的数据特征向量之间找到一个分离的超平面。但SVM分类器不太适合三个或多个分类问题。主要特征是要搜索的超平面是两类数据特征向量之间的最远距离,也就是说,不仅要找到一个超平面来分离两种类型的数据,而且还要最大化边界距离。支持向量机分类器SVM构建过程的伪代码如下表3所示。
表3:支持向量机分类器伪代码
与支持向量机分类器类似,随机森林分类器也属于集成学习方法,为了解决单个模型或者某一组参数的模型所固有的缺陷,从而将多棵决策树整合成森林即随机森林分类器,取长补短,避免局限性,并合起来用来预测最终结果。
随机森林分类器的实施过程:首先,采用重采样bootstrap方法(有放回的随机采样方式)生成O个训练集,然后,对于每个训练集,构造一棵决策树,在节点找特征进行分裂的时候,并不是对所有特征找到能使得指标(如信息增益)最大的,而是在特征中随机抽取一部分特征,在抽到的特征中间找到最优解,应用于节点,进行分裂。随机森林的方法由于有了bagging,也就是集成的思想在,实际上相当于对于样本和特征都进行了采样,避免过拟合。随机森林分类器(RF)构建过程的伪代码如下表4所示。
表4:随机森林分类器伪代码
梯度提升树是逐步增强或逐渐改善错误的方法。由斯坦福大学的杰罗姆·H·弗里德曼设计;梯度提升树分类器构建过程的伪代码如下表5所示。
表5:梯度提升树分类器伪代码
综上,梯度提升树分类器、支持向量机分类器、随机森林分类器构建完成。
步骤303、由于癫痫数据集的数量限制,为防止分类器的过拟合,在分类器训练过程中采用十折交叉验证的方式进行分类器训练。
本发明采用十折交叉验证手段,来训练更可靠更稳定的分类器,对训练集以随机抽取的方式分为十个子集,选定其中的一个子集为验证集,其他九个子集为真正用于训练的训练集,进行试验。重复上面过程,再将训练集随机划分为十个子集,选择其中一个子集为验证集,其他九个子集为真正的训练集。保证训练集里的每一个数据都能够得到训练,这样训练出来的分类模型鲁棒性高,同时可以防止分类模型过拟合。
步骤304、进行10次十折交叉验证操作之后,每次验证都会得到一个训练误差。采用十次十折交叉验证平均值CVe作为最终的交叉验证误差(CV Error),定义为: eq为第q次十折交叉验证的误差,其中h是交叉验证的折数,是预测值,Yh是实际值。
步骤305、采用网格搜索优化器对梯度提升树分类器GBM、支持向量机分类器SVM和随机森林分类器RF三种分类器在训练过程中产生的分类模型的参数进行寻优(分类模型优化),运用变步长的策略对网格搜索优化器进行改造,然后将网格搜索优化器应用于分类器做优化处理。
步骤306、利用验证集在分类器上进行验证,确定最终的三个分类器。
在步骤一和步骤三之间对所处理的数据进行打标签处理。对所处理的数据中三种类型的脑电信号数据子集(连续癫痫发作脑电信号、间歇癫痫发作的脑电信号、健康被试的脑电信号)进行打标签处理,包括:连续癫痫发作脑电信号{S}打正标签“1”、间歇癫痫发作脑电信号{F/N}打标签“0”、健康被试脑电信号{Z/O}打负标签“-1”。
步骤四、通过多种分类器评估指标来评估三个分类器的分类性能,确定最适合进行分类多种类别数据的分类器(多指标评估分类器)。
虽然大多数提出的发明都侧重于追求分类器的准确性,但在医学检测过程中,其实更加注重错分率也就是对错误分类因子的分析,所以本发明对能够深入挖掘分类错误因素的指标更为感兴趣。因此,本发明联合使用混淆矩阵、精确回忆曲线、接收器工作特性曲线和曲线下面积等多项分类器评估指标来验证所提出的癫痫检测分类器的性能,而不仅仅局限于通过分类准确率来评判一个分类器的好坏,最终本发明得到梯度提升树分类器取得了最为优秀的识别效果。
步骤401、利用测试集对步骤三中所确定的梯度提升树分类器GBM、支持向量机分类器SVM和随机森林分类器RF进行测试。
步骤402、定义多个评估指标来评估分类器的性能,评估指标包括:混淆矩阵(表6虚线框里)、受试者工作特征曲线、区域下面积和精确率-召回率曲线来验证所提出的分类器的性能。
表6中定义了癫痫脑电数据进行三分类的性能指标,如准确性、敏感性、特异性、准确性和混淆矩阵等。
表6应用于三分类的分类器的评估指标定义
步骤403、对每个分类器进行性能指标对比分析,最终本发明确定梯度提升树分类器取得了所需的最为优秀的识别结果。
本发明通过比较梯度提升树分类器GBM、随机森林RF和支持向量机SVM三种分类器处理三分类癫痫脑电数据集{S},{F/N},{Z/O}的混淆矩阵,这三种分类的混淆矩阵的结果如图4所示,左侧为梯度提升树-混淆矩阵分类结果,中间为随机森林-混淆矩阵分类结果,右侧为支持向量机-混淆矩阵分类结果,以左侧梯度提升树-混淆矩阵分类结果为例,每个方格不同深度的灰度代表不同类别的脑电信号,每次从五百个数据选一百个数据进行训练,每个格子里面的数值表示不同类别脑电信号的数据个数,所有格子里面数值加和为一百。通过实验验证梯度提升树分类器GBM比随机森林RF分类器和支持向量机SVM分类器取得有更高的分类准确率(混淆矩阵对角线上的预测值越大越好,非对角线上的预测值越小越好),它不仅保证了对真阳性和真阴性的高预测,而且还避免了假阳性和真阴性率的误差。
近年来,随着机器学习技术的发展,为构建一个具有高精度或高召回率的分类器,本发明采用接收器工作特性(ROC)和评估指标的曲线下面积(AUC)来评估分类器性能。如图5所示绘制ROC曲线表示受试者工作特性曲线与区域下面积之间的关系,横坐标为假阳性,纵坐标为真阳性,分类器必须提供对每个样本判断为正或负的置信值。AUC为ROC曲线的下面积进行分类器的整体性能评估的自然度量。通常,AUC值在0.5-1之间。GBM,RF和SVM方法与使用子集{F/N},{Z/O},{S}的GSO之间的AUC值分别为0.9695,0.956和0.9538。精确回忆曲线(PRC)在分类和检索领域具有广泛的应用,以表示分类和检索的性能。图6中PR曲线表示精确率和召回率之间的关系。垂直轴的精确率意味着正确地预测阳性样品与所有阳性样品的比率;水平轴的召回率意味着正确地预测正样本与真实样本的比率。当精确率和召回率很高时,可以确定分类器是好的。可以看出,梯度提升器和网格搜索优化器(GBM+GSO)在如图6所示的三种分类器中具有最大优势和最佳性能。
步骤五、将所确定的分类器的识别结果为后续使用提供决策支持。识别结果包括三种类别:安全、低级预警和高级预警。
如果分类器输出的结果为负标签“-1”,则病人处于健康的无发作正常状态,识别结果显示为安全。
如果分类器输出的结果为标签“0”,则表明病人处于间歇癫痫发作状态,需要及时给喂药进行控制病情恶化,识别结果显示为低级预警。
如果分类器输出的结果为正标签“1”,则表明病人处于连续癫痫发作状态,需要对癫痫病人进行紧急控制和治疗,识别结果显示为高级预警。
通过设计有效地癫痫检测脑电信号方案,开发一款手机的应用或许融入一些智能睡眠分析软件中,通过无线传感器网络传递脑电信号到手机终端。特别用于在晚上在家中检测和报警,防止夜间患者或未知的婴儿癫痫的突然死亡,改善患者和护理人员的生活质量。
Claims (8)
1.基于监督梯度提升器的多级癫痫脑电信号自动识别方法,其特征在于,所述步骤如下:
步骤一:对采集的原始脑电信号进行预处理;
采集的原始脑电信号包括三种类别,具体为连续癫痫发作脑电信号{S}、间歇癫痫发作的脑电信号{F/N}和健康被试的脑电信号{Z/O};
步骤二:对预处理后的脑电信号进行频域特征提取,具体步骤如下;
步骤201、运用小波变换的方法进行频域分析,通过对步骤一得到的脑电信号做Symlet小波分解,提取脑电信号的五个频域特征;
步骤202、通过主成分分析法进行降低维度的处理;
步骤203、对降维后的低维特征进行批量归一化处理,统一维度和格式;
步骤204、分别计算统一维度后的五个频率子带的统计特征,所述统计特征包括能量均值及能量标准方差,为导入分类器做准备;
步骤三、构建分类癫痫脑电信号的分类器,所述分类器包括梯度提升树分类器、支持向量机分类器和随机森林分类器;
步骤301、将原始脑电信号分为训练集和测试集,并以随机的方式将训练集分为十个子集,取其中一个子集作为验证集,其他九个子集作为真正用于训练的训练集;
步骤302、向每个分类器导入脑电数据,准备进行分类器训练;
步骤303、在分类器训练过程中采用十折交叉验证的方式进行分类器训练;
步骤304、采用每次的训练误差的平均值作为最终的交叉验证误差;
步骤305、采用网格搜索优化器对三种分类器在训练过程中产生的分类模型的参数进行寻优,运用变步长的策略对网格搜索优化器进行改造,然后将网格搜索优化器应用于分类器做优化处理;
步骤306、利用验证集在分类器上进行验证,确定最终的三个分类器;
在步骤一和步骤三之间对所处理的数据进行打标签处理;
步骤四、通过多种分类器评估指标来评估三个分类器的分类性能,确定最适合进行分类多种类别数据的分类器;
步骤401、利用测试集对步骤三中所确定的三个分类器进行测试;
步骤402、定义多个评估指标来评估分类器的性能,评估指标包括:混淆矩阵、受试者工作特征曲线、区域下面积和精确率-召回率曲线;
步骤402、对每个分类器进行性能指标对比分析,最终确定取得所需的识别结果的分类器;
步骤五、将所确定的分类器的识别结果为后续使用提供决策支持。
2.如权利要求1所述的基于监督梯度提升器的多级癫痫脑电信号自动识别方法,其特征在于,在步骤一和步骤三之间对所处理的数据进行打标签处理;包括:连续癫痫发作脑电信号{S}打正标签“1”、间歇癫痫发作脑电信号{F/N}打标签“0”、健康被试脑电信号{Z/O}打负标签“-1”。
3.如权利要求1所述的基于监督梯度提升器的多级癫痫脑电信号自动识别方法,其特征在于,步骤一通过MATLAB的开源工具箱EEGlab对采集来的原始脑电信号进行预处理。
4.如权利要求1所述的基于监督梯度提升器的多级癫痫脑电信号自动识别方法,其特征在于,步骤二中步骤201具体为通过四级Symlet小波将{S},{F/N},{Z/O}中的数据分解为包括alpha波、beta波、delta波、theta波和gamma波的五个频率子带。
5.如权利要求1所述的基于监督梯度提升器的多级癫痫脑电信号自动识别方法,其特征在于,步骤二中步骤202主成分分析法的伪代码如下所示:
算法:主成分分析法;
输入:脑电信号的频域特征s(t1),s(t2),…,s(tp),p是频域特征的个数;
过程:
a:集中所有频域特征:
b:计算频域特征的协方差矩阵:Cov(s(ti),s(tj))
c:求解相关系数矩阵R:
R=(rij)p×p=(r(s(ti),s(tj)))p×p,(rij=rji,rii=1)
d:求解相关系数矩阵的特征值:
λ1≥λ2≥…λp≥0
e:确定主成分的数量:m
f:计算相应的特征向量:
g:计算主要成分:
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811283792 | 2018-10-31 | ||
CN2018112837922 | 2018-10-31 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109934089A CN109934089A (zh) | 2019-06-25 |
CN109934089B true CN109934089B (zh) | 2020-10-23 |
Family
ID=66984915
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910034509.0A Active CN109934089B (zh) | 2018-10-31 | 2019-01-15 | 基于监督梯度提升器的多级癫痫脑电信号自动识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109934089B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113436728A (zh) * | 2021-07-05 | 2021-09-24 | 复旦大学附属儿科医院 | 新生儿临床视频脑电图自动分析的方法及设备 |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110289097A (zh) * | 2019-07-02 | 2019-09-27 | 重庆大学 | 一种基于Xgboost神经网络堆叠模型的模式识别诊断系统 |
CN110543832A (zh) * | 2019-08-13 | 2019-12-06 | 同济大学 | 基于随机森林和卷积神经网络的脑电数据分类方法 |
CN110448273B (zh) * | 2019-08-29 | 2021-03-30 | 江南大学 | 一种基于支持向量机的低功耗癫痫预测电路 |
CN111273131A (zh) * | 2020-03-17 | 2020-06-12 | 辽宁工程技术大学 | 基于能量特性与随机森林的光伏并网发电孤岛检测方法 |
CN111553482B (zh) * | 2020-04-09 | 2023-08-08 | 哈尔滨工业大学 | 机器学习模型超参数的调优方法 |
CN111956221B (zh) * | 2020-09-07 | 2022-06-07 | 南京医科大学 | 一种基于小波散射因子与lstm神经网络模型的颞叶癫痫分类方法 |
CN112842342B (zh) * | 2021-01-25 | 2022-03-29 | 北京航空航天大学 | 一种结合希尔伯特曲线和集成学习的心电磁信号分类方法 |
CN112890834B (zh) * | 2021-03-01 | 2022-05-13 | 福建脉行人生医疗科技有限公司 | 面向注意力识别的基于机器学习的眼电信号分类器 |
CN113255591A (zh) * | 2021-06-25 | 2021-08-13 | 四川九通智路科技有限公司 | 一种基于随机森林和融合特征的轴承故障诊断方法 |
CN113723485B (zh) * | 2021-08-23 | 2023-06-06 | 天津大学 | 一种轻微肝性脑病脑影像超图处理方法 |
CN113729641A (zh) * | 2021-10-12 | 2021-12-03 | 南京润楠医疗电子研究院有限公司 | 一种基于条件对抗网络的非接触式睡眠分期系统 |
CN114093501B (zh) * | 2021-10-19 | 2024-07-02 | 杭州电子科技大学 | 基于同步视频与脑电的儿童运动性癫痫智能辅助分析方法 |
CN114176519A (zh) * | 2021-12-09 | 2022-03-15 | 复旦大学 | 一种非接触式心电信号质量分类方法 |
CN114176609A (zh) * | 2021-12-27 | 2022-03-15 | 杭州电子科技大学 | 基于脑电信号的卒中严重程度分类预测模型建立方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000126148A (ja) * | 1998-10-28 | 2000-05-09 | Nec Corp | 脳波データ処理装置及び記録媒体 |
CN104523268A (zh) * | 2015-01-15 | 2015-04-22 | 江南大学 | 一种具备迁移学习能力的脑电信号识别模糊系统方法 |
CN104887222A (zh) * | 2015-05-11 | 2015-09-09 | 重庆大学 | 可逆化脑电信号分析方法 |
CN104887224A (zh) * | 2015-05-29 | 2015-09-09 | 北京航空航天大学 | 面向癫痫脑电信号的特征提取与自动识别方法 |
CN106137185A (zh) * | 2016-06-21 | 2016-11-23 | 华南理工大学 | 一种基于模板构造小波的癫痫特征波检测方法 |
CN106874694A (zh) * | 2017-03-20 | 2017-06-20 | 江南大学 | 一种面向癫痫脑电信号识别的智能诊断系统 |
CN107153825A (zh) * | 2017-05-23 | 2017-09-12 | 哈尔滨理工大学 | 一种基于支持向量机的癫痫脑电分类方法 |
CN107530012A (zh) * | 2015-02-16 | 2018-01-02 | 内森·英特拉托 | 用于脑活动解析的系统与方法 |
CN108320800A (zh) * | 2018-04-16 | 2018-07-24 | 吉林大学 | 癫痫病发作检测和发作前预测的脑电数据分析系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10349851B2 (en) * | 2013-07-30 | 2019-07-16 | Yrt Limited | Method, non-transitory computer readable medium and apparatus for arousal intensity scoring |
-
2019
- 2019-01-15 CN CN201910034509.0A patent/CN109934089B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000126148A (ja) * | 1998-10-28 | 2000-05-09 | Nec Corp | 脳波データ処理装置及び記録媒体 |
CN104523268A (zh) * | 2015-01-15 | 2015-04-22 | 江南大学 | 一种具备迁移学习能力的脑电信号识别模糊系统方法 |
CN107530012A (zh) * | 2015-02-16 | 2018-01-02 | 内森·英特拉托 | 用于脑活动解析的系统与方法 |
CN104887222A (zh) * | 2015-05-11 | 2015-09-09 | 重庆大学 | 可逆化脑电信号分析方法 |
CN104887224A (zh) * | 2015-05-29 | 2015-09-09 | 北京航空航天大学 | 面向癫痫脑电信号的特征提取与自动识别方法 |
CN106137185A (zh) * | 2016-06-21 | 2016-11-23 | 华南理工大学 | 一种基于模板构造小波的癫痫特征波检测方法 |
CN106874694A (zh) * | 2017-03-20 | 2017-06-20 | 江南大学 | 一种面向癫痫脑电信号识别的智能诊断系统 |
CN107153825A (zh) * | 2017-05-23 | 2017-09-12 | 哈尔滨理工大学 | 一种基于支持向量机的癫痫脑电分类方法 |
CN108320800A (zh) * | 2018-04-16 | 2018-07-24 | 吉林大学 | 癫痫病发作检测和发作前预测的脑电数据分析系统 |
Non-Patent Citations (6)
Title |
---|
Review Wavelet-based EEG processing for computer-aided seizure detection and epilepsy diagnosis;Oliver Faust 等;《Seizure》;20150131;56-64 * |
基于PCA和小波变换的脑电信号的特征提取;孙宇舸 等;《2007中国控制与决策学术年会论文集》;20070701;669-672,676 * |
癫痫失神发作脑电信号小波分析与脑电信号区域特征提取;张睿;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20140215(第02期);E070-11 * |
癫痫脑电的小波方差分析;马薇 等;《科技传播》;20111031(第10期);100-101 * |
视觉诱发脑电信号的处理研究;王夏爽 等;《系统仿真学报》;20171231;第29卷;146-154 * |
高速列车牵引系统逆变器复合故障诊断;金林强;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20180315(第03期);C033-200 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113436728A (zh) * | 2021-07-05 | 2021-09-24 | 复旦大学附属儿科医院 | 新生儿临床视频脑电图自动分析的方法及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN109934089A (zh) | 2019-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109934089B (zh) | 基于监督梯度提升器的多级癫痫脑电信号自动识别方法 | |
Khalighi et al. | Automatic sleep staging: A computer assisted approach for optimal combination of features and polysomnographic channels | |
Hussein et al. | Optimized deep neural network architecture for robust detection of epileptic seizures using EEG signals | |
Hussein et al. | Epileptic seizure detection: A deep learning approach | |
CN108256629B (zh) | 基于卷积网络和自编码的eeg信号无监督特征学习方法 | |
CN111444747B (zh) | 一种基于迁移学习和空洞卷积的癫痫状态识别方法 | |
Sharmila et al. | Wavelet-based feature extraction for classification of epileptic seizure EEG signal | |
Kavitha et al. | On the use of wavelet domain and machine learning for the analysis of epileptic seizure detection from EEG signals | |
CN114366124B (zh) | 一种基于半监督深度卷积通道注意力单分类网络的癫痫脑电识别方法 | |
De Lucia et al. | A novel method for automated classification of epileptiform activity in the human electroencephalogram-based on independent component analysis | |
He et al. | A progressive deep wavelet cascade classification model for epilepsy detection | |
Yang et al. | A novel method for identifying electrocardiograms using an independent component analysis and principal component analysis network | |
Jadhav et al. | Automated sleep stage scoring using time-frequency spectra convolution neural network | |
Ibrahim et al. | Deep‐learning‐based seizure detection and prediction from electroencephalography signals | |
Liu et al. | Automatic sleep staging algorithm based on random forest and hidden Markov model | |
Raiesdana | Automated sleep staging of OSAs based on ICA preprocessing and consolidation of temporal correlations | |
Movahed et al. | A major depressive disorder diagnosis approach based on EEG signals using dictionary learning and functional connectivity features | |
Tripathi et al. | Automatic seizure detection and classification using super-resolution superlet transform and deep neural network-A preprocessing-less method | |
Abou-Abbas et al. | Patient-independent epileptic seizure detection by stable feature selection | |
Raja et al. | Existing Methodologies, Evaluation Metrics, Research Gaps, and Future Research Trends: A Sleep Stage Classification Framework | |
Jia et al. | Decision level fusion for pulse signal classification using multiple features | |
Ma et al. | Parallel Dual-Branch Fusion Network for Epileptic Seizure Prediction | |
Kinney-Lang et al. | Elucidating age-specific patterns from background electroencephalogram pediatric datasets via PARAFAC | |
Feudjio et al. | A novel use of discrete wavelet transform features in the prediction of epileptic seizures from EEG Data | |
Luo et al. | Subject-adaptive real-time sleep stage classification based on conditional random field |
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 |