1.基于深度自编码的脑电特征优化与癫痫发作检测方法,其特征在于包括如下步骤:
步骤1:将原始脑电信号进行预处理,然后按顺序排列脑电信号的通道,以每4s的脑电图为一个数据样本,进行离散傅里叶变换,并获得样本的幅值谱,将样本的幅值谱整理成合适的特征矩阵;
步骤2:设置卷积自编码器的网络结构;
步骤3:通过步骤1中的数据训练卷积自编码器,再将数据输入训练好的卷积自编码器的编码器对特征矩阵进行降维,得到低维特征;
步骤4:使用低维特征训练基于最近邻分类器的纠错输出编码模型;最后使用整个模型对新的脑电图进行癫痫发作检测,测试模型性能。
2.根据权利要求1所述的基于深度自编码的脑电特征优化与癫痫发作检测方法,其特征在于所述步骤1的具体实现步骤如下:
1-1.按顺序排列脑电信号的通道,将脑电信号ref导联的21个通道按照从左到右、从前到后的顺序排列,并将长时间的脑电图分割成若干时长为4秒的数据样本,每两个相邻的样本中间有50%的重叠部分;
1-2.定义癫痫发作前的十分钟为发作前期,癫痫发作时记为发作期,癫痫发作前后四小时以外的时间记为癫痫发作间期;
1-3.使用采样率为1000hz的脑电图,每一个4秒的样本为21*4000的实数矩阵,将样本中的每个通道进行离散傅里叶变换(dft):
其中,n表示每帧的长度,k=0,...,n-1;并由下式获取各通道的幅值谱p(k):
p(k)=|xk|
其中,xk为离散傅里叶变换的结果;将得到的幅值谱去除直流分量和工频干扰,并取87hz以下的信息,即可得到21*336维的原始幅值谱特征;
1-4.整理原始幅值谱特征,将21*336的特征矩阵按列平均分成4部分,即4个21*84的矩阵,再拼接成84*84的正方形特征矩阵。
3.根据权利要求1或2所述的基于深度自编码的脑电特征优化与癫痫发作检测方法,其特征在于步骤2所述的设置卷积自编码器的网络结构,具体实现如下:
2-1.卷积自编码器包含编码器和解码器;其中编码器主要由卷积层和池化层组成,解码器由反卷积层和反池化层组成;编码器将输入信号x变换成编码信号y,而解码器将编码信号y转换成输出信号z;
y=f(x)
z=g(y)=g(f(x))
2-2.编码器中包含一个输入层、3个卷积层和3个最大池化层,解码器包含4反卷积层和3个上采样层;卷积层中的计算公式为:
其中
在每个卷积层后,使用下采样降低特征维数或使用上采样还原特征维数;这里使用最大池化层完成下采样,即滤波器的输出为输入数据的最大值,最大池化层的计算公式如下:
xi=round(xo(wi/wo))
yi=round(yo(hi/ho))
其中,(xi,yi)、(xo,yo)为输入和输出矩阵元素坐标,wi、wo、hi和ho是输入和输出矩阵的行数和列数;
2-3.设置网络参数,其中矩阵边缘需要补0填充,用标记”same”表示填充后该层的输出矩阵的尺寸与输入矩阵相同。
4.根据权利要求3所述的基于深度自编码的脑电特征优化与癫痫发作检测方法,其特征在于步骤3中使用步骤1中所述的数据训练卷积自编码器,再将数据输入训练好的卷积自编码器的编码器对特征矩阵进行降维,得到低维特征,具体实现方法如下:
3-1.卷积自编码器是从训练数据中自动学习的,训练卷积自编码器时只需训练数据而不需要数据的标签;卷积自编码器训练时使用随机梯度下降优化编码/解码函数的参数以最小化特征重建损失;
l(x,z)=l(x,g(f(x)))
l表示损失函数,这里以交叉熵误差作为损失函数;使用adam作为优化器,训练卷积自编码器;卷积自编码器训练完成后,其中的编码器就是特征降维模型,解码器中第三个最大池化层的输出即特征降维模型的输出;
3-2.将步骤1中整理好的特征矩阵输入模型,进行特征重建,得到低维特征。
5.根据权利要求4所述的基于深度自编码的脑电特征优化与癫痫发作检测方法,其特征在于步骤4中使用低维特征训练基于最近邻分类器的纠错输出编码模型;最后使用整个模型对新的脑电图进行癫痫发作检测,测试模型性能,具体实现方法如下:
4-1.使用低维特征和数据的标签训练基于最近邻分类器的纠错输出编码模型;该模型中采用一对一编码,使用10个最近邻分类器进行分类,得到编码矩阵m,根据以下公式预测样本所属类别:
其中,mkl为编码矩阵m中的元素,对应第k类第l个分类器的编码;sl为第l个分类器的分类结果;函数g功能为计算mkl和sl间的海明距离;
4-2.将新的脑电图数据按照同样的方式计算幅值谱并整理成84*84的特征矩阵;输入训练好的编码器中,得到低维特征;
4-3.将降维后的低维特征输入基于最近邻分类器的纠错输出编码模型,得到测试样本的类别。