CN110543832A - 基于随机森林和卷积神经网络的脑电数据分类方法 - Google Patents
基于随机森林和卷积神经网络的脑电数据分类方法 Download PDFInfo
- Publication number
- CN110543832A CN110543832A CN201910745397.XA CN201910745397A CN110543832A CN 110543832 A CN110543832 A CN 110543832A CN 201910745397 A CN201910745397 A CN 201910745397A CN 110543832 A CN110543832 A CN 110543832A
- Authority
- CN
- China
- Prior art keywords
- neural network
- electroencephalogram data
- convolutional neural
- random forest
- data
- 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.)
- Pending
Links
- 238000007637 random forest analysis Methods 0.000 title claims abstract description 83
- 238000013527 convolutional neural network Methods 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000012549 training Methods 0.000 claims abstract description 27
- 238000012216 screening Methods 0.000 claims abstract description 19
- 238000012360 testing method Methods 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000003066 decision tree Methods 0.000 claims description 38
- 210000002569 neuron Anatomy 0.000 claims description 6
- 238000011176 pooling Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 6
- 101100119135 Mus musculus Esrrb gene Proteins 0.000 claims description 4
- 238000013528 artificial neural network Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 210000004556 brain Anatomy 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 5
- 238000013135 deep learning Methods 0.000 description 4
- 238000010801 machine learning Methods 0.000 description 4
- 238000002790 cross-validation Methods 0.000 description 3
- 238000012800 visualization Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000000007 visual effect Effects 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/369—Electroencephalography [EEG]
-
- 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
-
- 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/7271—Specific aspects of physiological measurement analysis
- A61B5/7282—Event detection, e.g. detecting unique waveforms indicative of a medical condition
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- 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
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Psychiatry (AREA)
- Pathology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Mathematical Physics (AREA)
- Signal Processing (AREA)
- Computing Systems (AREA)
- Physiology (AREA)
- Computational Linguistics (AREA)
- Software Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Fuzzy Systems (AREA)
- Psychology (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本发明涉及一种基于随机森林和卷积神经网络的脑电数据分类方法,包括:S1、采集原始的时序脑电数据,并将时序脑电数据转换为频域脑电数据;S2、基于随机森林方法,对频域脑电数据进行导联筛选,以对频域脑电数据进行降维处理;S3、构建卷积神经网络,将降维之后的频域脑电数据输入给卷积神经网络,对卷积神经网络进行训练和测试,得到训练好的卷积神经网络;S4、结合随机森林方法以及训练好的卷积神经网络,对经过傅里叶转换的实际脑电数据依次进行导联筛选和分类。与现有技术相比,本发明利用随机森林筛选重要的特征,降低了脑电数据维度,减少了计算量;结合卷积神经网络对降维的脑电数据进行拟合,保证了脑电数据分类准确率。
Description
技术领域
本发明涉及脑电信号的数据分析技术领域,尤其是涉及一种基于随机森林和卷积神经网络的脑电数据分类方法。
背景技术
利用机器学习和深度学习方法对脑电信号进行分类是近年来脑电信号研究领域的一个热门方向,但由于脑电数据的维度高、计算量大,导致影响分类结果的数据特征过多,因此,采用机器学习方法对脑电数据进行分析时,通常无法找到一个合适的模型对脑电数据进行分类;采用深度学习的方法虽然能很好的拟合模型,可是随着数据规模的增大,必然会使计算量呈指数级增加。
综上所述,在面对维度高、计算量大的脑电数据时,现有技术没有从脑电数据中提取出对分类结果产生重要作用的特征,机器学习和深度学习方法均不能在保证分类准确率的前提下,实现快速输出分类结果的目的。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于随机森林和卷积神经网络的脑电数据分类方法。
本发明的目的可以通过以下技术方案来实现:一种基于随机森林和卷积神经网络的脑电数据分类方法,包括以下步骤:
S1、采集原始的时序脑电数据,对原始时序脑电数据进行预处理,之后再通过傅里叶转换,将时序脑电数据转换为频域脑电数据;
S2、基于随机森林方法,对频域脑电数据进行导联筛选,得到脑电数据中各导联上所有特征的重要程度,以对频域脑电数据进行降维处理;
S3、构建卷积神经网络,将降维之后的频域脑电数据输入给卷积神经网络,对卷积神经网络进行训练和测试,得到训练好的卷积神经网络;
S4、结合随机森林方法以及训练好的卷积神经网络,对经过傅里叶转换的实际脑电数据依次进行导联筛选和分类。
优选地,所述步骤S1中对原始时序脑电数据进行预处理包括对脑电数据进行滤波和去除伪迹的操作。
优选地,所述步骤S2具体包括以下步骤:
S21、建立随机森林;
S22、获取频域脑电数据的灰度图,并将每个频率下的灰度图转换成对应的一维向量,以作为随机森林的输入;
S23、随机森林对所有的输入进行导联筛选,计算特征的重要程度;
S24、根据特征的重要程度,对频域脑电数据进行降维处理。
优选地,所述步骤S21具体包括以下步骤:
S211、定义样本数量为T,从所有样本中放回的随机抽样T个样本,将随机抽样的T个样本作为决策树根节点处的样本,以训练决策树;
S212、定义样本集中每个样本的特征数量为F,则每次决策树的节点需要进行分裂时,首先从这F个特征中选出f个特征,然后采用信息增益、信息增益比或基尼指数这三种策略,从f个特征中选出一个特征作为该节点的分裂属性,其中,f<<F;
S213、重复步骤S212的过程,以进行决策树节点的分裂,直到无法再继续分裂为止,完成一颗决策树的训练;
S214、完成T棵决策树的训练后,由这T棵训练好的决策树构成随机森林。
优选地,所述步骤S22中随机森林的输入具体为W维的向量:
W=A×B
其中,A表示导联的个数,B表示一个导联上所有频率的数量。
优选的,所述步骤S23和步骤S24中特征具体为导联上对应频率的数据。
优选地,所述步骤S23中计算特征的重要程度的具体过程为:
S231、使用袋外数据计算随机森林中每一棵决策树的误差err1;
S232、在所有袋外数据的相应特征位置加入噪声扰动,并计算此时决策树的误差err2;
S233、计算特征的重要程度为:
其中,wf表示特征f的重要程度,T为随机森林中决策树的总棵数,i表示随机森林中的第i棵决策树,i_err2表示在袋外数据相应特征位置加入噪声扰动后计算的第i棵决策树的误差,i_err1表示使用袋外数据计算的第i棵决策树的误差。
优选地,所述步骤S3具体包括以下步骤:
S31、构建卷积神经网络,所述卷积神经网络包括依次连接的卷积层、最大池化层和全连接层;
S32、将降维之后的频域脑电数据分为训练集和测试集;
S33、将训练集输入卷积神经网络,以预设的训练轮数进行训练;
S34、每完成一次预设的训练轮数,则将测试集输入卷积神经网络,记录一次卷积神经网络输出分类结果的准确率;
S35、重复步骤S34,并判断卷积神经网络输出分类结果的准确率是否收敛,若收敛,则该卷积神经网络已经训练好,否则返回步骤S33。
优选地,所述步骤S31中卷积层由三个子卷积层构成,三个子卷积层的卷积核大小与步长均一致,全连接层由第一全连接层和第二全连接层组成,所述第一全连接层的神经元个数与随机森林的输入的向量维数一致,所述第二全连接层的神经元个数为脑电数据分类结果的类别数。
与现有技术相比,本发明使用随机森林与卷积神经网络的组合模型对脑电数据进行分类,首先采用随机森林计算所有特征对于分类结果的重要程度,并根据计算的特征的重要程度进行导联筛选,筛选出重要的特征,降低了数据维度,大大减少了计算量;然后基于卷积神经网络优良的数据拟合性,利用反向传播算法去学习适合脑电数据的参数,以实现脑电数据分类的目的,通过结合随机森林与卷积神经网络,能够在保证分类准确率的前提下,降低脑电数据的维度与计算量,从而提高脑电数据分类的速度。
附图说明
图1为本发明的方法流程图;
图2为随机森林的工作流程示意图;
图3为实施例中八个随机森林对应的分类准确率示意图;
图4为实施例中包含100棵决策树的随机森林特征重要程度的可视化示意图;
图5为实施例中八个随机森林的特征重要程度平均值的可视化示意图;
图6为实施例中求取特征重要程度平均值后的导联重要程度示意图;
图7为实施例中导联筛选之后八个随机森林对应的分类准确率示意图;
图8为实施例中卷积神经网络的结构图;
图9为实施例中导联筛选之后卷积神经网络的输出分类准确率示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
如图1所示,一种基于随机森林和卷积神经网络的脑电数据分类方法,包括以下步骤:
S1、采集原始的时序脑电数据,对原始时序脑电数据进行预处理,之后再通过傅里叶转换,将时序脑电数据转换为频域脑电数据;
S2、基于随机森林方法,对频域脑电数据进行导联筛选,得到脑电数据中各导联上所有特征的重要程度,以对频域脑电数据进行降维处理;
S3、构建卷积神经网络,将降维之后的频域脑电数据输入给卷积神经网络,对卷积神经网络进行训练和测试,得到训练好的卷积神经网络;
S4、结合随机森林方法以及训练好的卷积神经网络,对经过傅里叶转换的实际脑电数据依次进行导联筛选和分类。
本实施例中,采集的原始时序脑电数据包括56个正样本和54个负样本,在数据采集条件相同的条件下,均为闭眼时的静息态脑电数据,其中,正样本为受试者在正常情况下的脑电数据,负样本为受试者在有外界噪声干扰情况下的脑电数据。每个受试者均采集180秒的数据,采集完原始时序脑电数据之后,首先对数据进行滤波、去除伪迹等预处理操作,之后使用快速傅里叶变换提取脑电数据中的频率信息,针对频率提取之后的脑电数据,本实施例还进行了min-max归一化处理,即:
其中,min表示脑电数据集中特征x的最小值,max表示脑电数据集中特征x的最大值,使用min-max的方法将数据进行归一化,使得预处理的数据被限定在一定的范围内(比如[0,1]或者[-1,1]),从而消除奇异样本数据导致的不良影响。
步骤S2中随机森林的工作流程如图2所示,随机森林由很多棵决策树所构成,森林的建立方式完全随机,因此被称作随机森林。森林构建完成之后,每次有新样本输入时,森林中的每棵决策树都会对这个样本进行分类判断,并以一种投票的方式来进行决策,将票数最多的类别作为最后的分类结果;
随机森林的基分类器是决策树,决策树是一种树型结构,决策树的叶节点代表了相应的分类类别。使用决策树进行分类的过程就是从树的根节点出发,通过特征中的不同属性选出决策树的不同分支,最终到达叶节点,结束分类的过程。
本实施例使用随机森林方法进行导联筛选的具体过程为:
一、随机森林输入的每个样本都是一个由多个特征组成的向量,因此需要将傅里叶变换后的灰度图拉成一维向量,作为随机森林的输入。直观上来看,随机森林的输入为1024维的向量(32导联x 32个频率值),每个特征代表了每个导联上对应频率的数据;
本实施例中,建立了八个不同决策树棵数的随机森林,森林中决策树的棵树从40棵开始,每次增加20棵,一直增加到180棵,图3中所示八个随机森林对应的分类准确率,横坐标为随机森林中决策树的棵树,纵坐标为分类的准确率,图3中这八个随机森林模型的分类准确率均在98%以上,当设置决策树棵树为80时,其对应的分类准确率最高,达到99%,这是由于随机森林参数少,从而能在小型数据集上产生优良的表现;
使用随机森林对所有数据进行分类之后,可以得到训练好的模型对每个特征的敏感程度,将上述八个随机森林对每个特征计算的重要程度进行可视化,获取到每个特征的重要程度之后,由于输入数据自身的特点,每个特征都有其特定的空间和频域上的意义,将所有的特征重要程度还原成32x 32的图片之后可将其进行可视化,图4为100棵决策树的随机森林特征重要程度的可视化示意图,由图4可知,第2、5、7、20、22这5个导联在特征的分类中起到了主导作用,图中颜色较暗的位置对应的值较小,说明输入的样本中对应位置的特征对分类结果的贡献不大,即特征重要程度不高,而颜色较亮的位置是随机森林中对分类结果影响较大的特征;
二、随机森林由于其初始化构建的过程中存在随机性,这可能导致多次重复运行得到的特征重要程度不同,因此较好的方式是使用不同配置的随机森林求取特征重要程度的平均值,因此在得到不同特征对分类结果的重要程度之后,可以利用此结果对输入的数据进行降维,将降维之后的数据再次使用随机森林进行分类,通过将降维前后分类结果进行对比就可以评估出选取的特征对分类结果的重要程度,即由于随机森林中决策树自身存在的随机性,每个随机森林模型给出的特征重要程度可能不同,为了降低随机误差对特征重要程度计算带来的影响,本实施例将上述八个随机森林给出的特征重要程度进行了平均,并从中找出了对最终分类结果影响较大的5个导联,图5所示为八个随机森林平均特征重要程度的可视化示意图,横坐标表示不同的频率,纵坐标表示不同的导联,同样沿着第一个维度求和,算出每个导联上的所有特征对分类结果的重要程度,图6所示为求取平均值之后各个导联上的特征对分类结果的贡献,横坐标表示不同的导联,纵坐标表示导联对分类结果的重要程度,从图6可以明显看出,第2、5、7、20、22这5个导联对分类结果的影响较大,这其中第20个导联最为重要,从图5中也能看到,第20个导联的亮度最大;另一方面,观察图5中这5行相应位置的亮度,我们可以看到这5行中频率数值小的像素点亮度更大,这说明导联上的低频成份对分类的影响更大;
三、选出这5个导联之后,对脑电数据进行降维处理,从1024维的特征中筛选出了160维(5个导联x 32维特征),再次使用随机森林对筛选后的数据进行分类,图7所示为导联筛选之后八个随机森林对应的分类准确率,图7中横坐标为随机森林中决策树的棵树,纵坐标为分类的准确率,从图7中可以看出,八个随机森林的分类准确率均在98%以上,对比图3未经过导联筛选的八个随机森林的分类准确率,仅降低了0.4%左右,这说明使用随机森林进行特征筛选方法是合理的,通过这种方式,特征数量由1024维降低至160维,同时保证了分类准确率。
步骤S3中构建的卷积神经网络结构如图8所示,网络中包含三个卷积层和两个全连接层,relu被用作网络中的激活函数,每个卷积层之前都使用batchnorm(批标准化层)进行标准化,网络中所有卷积层中卷积核大小都为3,步长设置为1,降采样操作使用最大池化来完成,池化操作的大小为2x 2,第一个全连接层中神经元的个数为1024,第二个全连接层中神经元的个数为2(即设置的分类结果的类别数)。
本实施例在网络的末端并没有使用全局平均池化,而是在最后一个卷积层之后连接了两个全连接层,主要是出于两方面考虑:首先,采集的原始脑电数据包含了54个正样本和56个负样本,因此虽然全连接层会增加一些参数量,但是训练数据的大量增加使得这些参数并不会影响模型训练;此外,由于分类结果的类别数被设置为2,如果使用全局平均池化的方式进行二分类,那么最后一个卷积层中卷积核的数量就要设置为2,通过图8可以看到,经过S6之后输出的特征图大小为4x 4x 128,若直接使用卷积层将特征图减少到两张,信息量损失较大,因此在最后的特征图之后使用了两个全连接层。
步骤S3中对卷积神经网络进行训练和测试的具体过程为:将降维之后的正样本和负样本数据分别以3:1的比例划分到训练集和测试集中,卷积神经网络使用Adam优化方法进行参数更新,学习率为0.05,L2权重衰减系数设置为0.5,同时使用指数衰减对学习率进行衰减,每200次迭代对学习率进行一次衰减,衰减系数为0.95,每训练3个epoch后,在测试集上进行一次分类准确率的测试,模型训练过程中测试结果如图9所示,具体为为卷积神经网络迭代次数的分类准确率,图9的横坐标代表分类准确率记录的次数,每三轮记录一次结果,纵坐标代表分类准确率。从图9可以看到,训练过程中分类准确率在95%±0.4%上下波动,经过5次记录之后,此时卷积神经网络收敛,在测试集上的准确率可以达到97%。
为了进一步验证本发明方法的有效性,实施例还使用k折交叉验证的方法划分数据集,并统计最终对应的分类准确率:即将数据集等分成k份,将其中的k-1份用来训练,剩余的1份用来测试,然后取k次实验结果的均值作为最终的分类结果。本实施例中k的取值范围为2到7,表1所示为经过k折交叉验证的实施例的结果数据:
表1
k折 | 2 | 3 | 4 | 5 | 6 | 7 |
分类准确率 | 94.2% | 94.7% | 94.6% | 94.7% | 94.6% | 94.5% |
从表1可以看出,交叉验证的结果较为稳定,准确率均在94%以上,从实施例的结果表明,相对于其他机器学习和深度学习的脑电分类算法,本发明提出的组合模型,首先利用随机森林筛选出重要的特征,降低了数据维度,减少了计算量,并且准确度只是些许下降;又利用了神经网络自身的拟合性,利用反向传播算法去学习适合脑电数据的参数,克服了根据数据选择模型对分类结果影响的未知性。
Claims (9)
1.一种基于随机森林和卷积神经网络的脑电数据分类方法,其特征在于,包括以下步骤:
S1、采集原始的时序脑电数据,对原始时序脑电数据进行预处理,之后再通过傅里叶转换,将时序脑电数据转换为频域脑电数据;
S2、基于随机森林方法,对频域脑电数据进行导联筛选,得到脑电数据中各导联上所有特征的重要程度,以对频域脑电数据进行降维处理;
S3、构建卷积神经网络,将降维之后的频域脑电数据输入给卷积神经网络,对卷积神经网络进行训练和测试,得到训练好的卷积神经网络;
S4、结合随机森林方法以及训练好的卷积神经网络,对经过傅里叶转换的实际脑电数据依次进行导联筛选和分类。
2.根据权利要求1所述的一种基于随机森林和卷积神经网络的脑电数据分类方法,其特征在于,所述步骤S1中对原始时序脑电数据进行预处理包括对脑电数据进行滤波和去除伪迹的操作。
3.根据权利要求1所述的一种基于随机森林和卷积神经网络的脑电数据分类方法,其特征在于,所述步骤S2具体包括以下步骤:
S21、建立随机森林;
S22、获取频域脑电数据的灰度图,并将每个频率下的灰度图转换成对应的一维向量,以作为随机森林的输入;
S23、随机森林对所有的输入进行导联筛选,计算特征的重要程度;
S24、根据特征的重要程度,对频域脑电数据进行降维处理。
4.根据权利要求3所述的一种基于随机森林和卷积神经网络的脑电数据分类方法,其特征在于,所述步骤S21具体包括以下步骤:
S211、定义样本数量为T,从所有样本中放回的随机抽样T个样本,将随机抽样的T个样本作为决策树根节点处的样本,以训练决策树;
S212、定义样本集中每个样本的特征数量为F,则每次决策树的节点需要进行分裂时,首先从这F个特征中选出f个特征,然后采用信息增益、信息增益比或基尼指数这三种策略,从f个特征中选出一个特征作为该节点的分裂属性,其中,f<<F;
S213、重复步骤S212的过程,以进行决策树节点的分裂,直到无法再继续分裂为止,完成一颗决策树的训练;
S214、完成T棵决策树的训练后,由这T棵训练好的决策树构成随机森林。
5.根据权利要求4所述的一种基于随机森林和卷积神经网络的脑电数据分类方法,其特征在于,所述步骤S22中随机森林的输入具体为W维的向量:
W=A×B
其中,A表示导联的个数,B表示一个导联上所有频率的数量。
6.根据权利要求5所述的一种基于随机森林和卷积神经网络的脑电数据分类方法,其特征在于,所述步骤S23和步骤S24中特征具体为导联上对应频率的数据。
7.根据权利要求6所述的一种基于随机森林和卷积神经网络的脑电数据分类方法,其特征在于,所述步骤S23中计算特征的重要程度的具体过程为:
S231、使用袋外数据计算随机森林中每一棵决策树的误差err1;
S232、在所有袋外数据的相应特征位置加入噪声扰动,并计算此时决策树的误差err2;
S233、计算特征的重要程度为:
其中,wf表示特征f的重要程度,T为随机森林中决策树的总棵数,i表示随机森林中的第i棵决策树,i_err2表示在袋外数据相应特征位置加入噪声扰动后计算的第i棵决策树的误差,i_err1表示使用袋外数据计算的第i棵决策树的误差。
8.根据权利要求1所述的一种基于随机森林和卷积神经网络的脑电数据分类方法,其特征在于,所述步骤S3具体包括以下步骤:
S31、构建卷积神经网络,所述卷积神经网络包括依次连接的卷积层、最大池化层和全连接层;
S32、将降维之后的频域脑电数据分为训练集和测试集;
S33、将训练集输入卷积神经网络,以预设的训练轮数进行训练;
S34、每完成一次预设的训练轮数,则将测试集输入卷积神经网络,记录一次卷积神经网络输出分类结果的准确率;
S35、重复步骤S34,并判断卷积神经网络输出分类结果的准确率是否收敛,若收敛,则该卷积神经网络已经训练好,否则返回步骤S33。
9.根据权利要求8所述的一种基于随机森林和卷积神经网络的脑电数据分类方法,其特征在于,所述步骤S31中卷积层由三个子卷积层构成,三个子卷积层的卷积核大小与步长均一致,全连接层由第一全连接层和第二全连接层组成,所述第一全连接层的神经元个数与随机森林的输入的向量维数一致,所述第二全连接层的神经元个数为脑电数据分类结果的类别数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910745397.XA CN110543832A (zh) | 2019-08-13 | 2019-08-13 | 基于随机森林和卷积神经网络的脑电数据分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910745397.XA CN110543832A (zh) | 2019-08-13 | 2019-08-13 | 基于随机森林和卷积神经网络的脑电数据分类方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110543832A true CN110543832A (zh) | 2019-12-06 |
Family
ID=68711491
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910745397.XA Pending CN110543832A (zh) | 2019-08-13 | 2019-08-13 | 基于随机森林和卷积神经网络的脑电数据分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110543832A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111671423A (zh) * | 2020-06-18 | 2020-09-18 | 四川大学 | 一种eeg信号的表示方法、分类方法、可视化方法及介质 |
CN111768389A (zh) * | 2020-07-01 | 2020-10-13 | 武汉楚精灵医疗科技有限公司 | 基于卷积神经网络和随机森林的消化道操作自动计时方法 |
CN113539414A (zh) * | 2021-07-30 | 2021-10-22 | 中电药明数据科技(成都)有限公司 | 一种抗生素用药合理性预测方法及系统 |
CN113762104A (zh) * | 2021-08-23 | 2021-12-07 | 北京航空航天大学 | 基于脑电信号的情景意识评估方法、装置及电子设备 |
CN117589951A (zh) * | 2023-12-08 | 2024-02-23 | 山东工商学院 | 一种生鲜食品新鲜度检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102058413A (zh) * | 2010-12-03 | 2011-05-18 | 上海交通大学 | 基于连续小波变换的脑电信号警觉度检测方法 |
CN109493886A (zh) * | 2018-12-13 | 2019-03-19 | 西安电子科技大学 | 基于特征选择和优化的语音情感识别方法 |
CN109934089A (zh) * | 2018-10-31 | 2019-06-25 | 北京航空航天大学 | 基于监督梯度提升器的多级癫痫脑电信号自动识别方法 |
CN110070133A (zh) * | 2019-04-24 | 2019-07-30 | 北京工业大学 | 一种基于深度森林的脑功能网络分类方法 |
-
2019
- 2019-08-13 CN CN201910745397.XA patent/CN110543832A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102058413A (zh) * | 2010-12-03 | 2011-05-18 | 上海交通大学 | 基于连续小波变换的脑电信号警觉度检测方法 |
CN109934089A (zh) * | 2018-10-31 | 2019-06-25 | 北京航空航天大学 | 基于监督梯度提升器的多级癫痫脑电信号自动识别方法 |
CN109493886A (zh) * | 2018-12-13 | 2019-03-19 | 西安电子科技大学 | 基于特征选择和优化的语音情感识别方法 |
CN110070133A (zh) * | 2019-04-24 | 2019-07-30 | 北京工业大学 | 一种基于深度森林的脑功能网络分类方法 |
Non-Patent Citations (1)
Title |
---|
HONGYU LIU、LIANGHUA HE: "Convolutional Neural Network for Resting State EEG Signal Analysis", 《2019 IEEE 16TH INTERNATIONAL CONFERENCE ON NETWORKING》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111671423A (zh) * | 2020-06-18 | 2020-09-18 | 四川大学 | 一种eeg信号的表示方法、分类方法、可视化方法及介质 |
CN111671423B (zh) * | 2020-06-18 | 2022-02-18 | 四川大学 | 一种eeg信号的表示方法、分类方法、可视化方法及介质 |
CN111768389A (zh) * | 2020-07-01 | 2020-10-13 | 武汉楚精灵医疗科技有限公司 | 基于卷积神经网络和随机森林的消化道操作自动计时方法 |
CN113539414A (zh) * | 2021-07-30 | 2021-10-22 | 中电药明数据科技(成都)有限公司 | 一种抗生素用药合理性预测方法及系统 |
CN113762104A (zh) * | 2021-08-23 | 2021-12-07 | 北京航空航天大学 | 基于脑电信号的情景意识评估方法、装置及电子设备 |
CN117589951A (zh) * | 2023-12-08 | 2024-02-23 | 山东工商学院 | 一种生鲜食品新鲜度检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110543832A (zh) | 基于随机森林和卷积神经网络的脑电数据分类方法 | |
CN110163813B (zh) | 一种图像去雨方法、装置、可读存储介质及终端设备 | |
CN111460956B (zh) | 一种基于数据增强与损失加权的不平衡心电样本分类方法 | |
CN112232476A (zh) | 更新测试样本集的方法及装置 | |
CN112508110A (zh) | 一种基于深度学习的心电信号图的分类方法 | |
CN110929581A (zh) | 一种基于时空特征加权卷积神经网络的脑电信号识别方法 | |
CN113488060B (zh) | 一种基于变分信息瓶颈的声纹识别方法及系统 | |
CN111127387B (zh) | 无参考图像的质量的评价方法 | |
CN113723171B (zh) | 基于残差生成对抗网络的脑电信号去噪方法 | |
CN110070116A (zh) | 基于深度树状训练策略的分段式选择集成图像分类方法 | |
CN113295702B (zh) | 电气设备故障诊断模型训练方法和电气设备故障诊断方法 | |
CN111681182A (zh) | 图片修复方法、装置、终端设备以及存储介质 | |
CN113987236B (zh) | 基于图卷积网络的视觉检索模型的无监督训练方法和装置 | |
CN114037001A (zh) | 基于wgan-gp-c和度量学习的机械泵小样本故障诊断方法 | |
CN111368648B (zh) | 雷达辐射源个体识别方法、装置、电子设备及其存储介质 | |
CN117473414A (zh) | 一种基于少噪声时频图像的轴承故障位置识别方法 | |
CN113378620B (zh) | 监控视频噪声环境下跨摄像头行人重识别方法 | |
CN113283519A (zh) | 一种基于离散系数的深度神经网络近似模型分析方法 | |
CN114186593B (zh) | 噪声条件下轴承故障识别方法、装置、设备和存储介质 | |
CN116350239A (zh) | 一种脑电信号专注度分类方法及系统 | |
CN114818969A (zh) | 变点位置及类别的检测方法及装置 | |
CN114611560A (zh) | 基于卷积神经网络的ssvep脑电信号分类方法 | |
CN115017960A (zh) | 一种基于时空联合mlp网络的脑电信号分类方法及应用 | |
CN112446345A (zh) | 一种低质量三维人脸识别方法、系统、设备和存储介质 | |
CN112488238A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191206 |