CN112464902A - 基于多通道多维特征优化的脑电眨眼伪迹检测方法 - Google Patents
基于多通道多维特征优化的脑电眨眼伪迹检测方法 Download PDFInfo
- Publication number
- CN112464902A CN112464902A CN202011485932.1A CN202011485932A CN112464902A CN 112464902 A CN112464902 A CN 112464902A CN 202011485932 A CN202011485932 A CN 202011485932A CN 112464902 A CN112464902 A CN 112464902A
- Authority
- CN
- China
- Prior art keywords
- channel
- mean
- mean value
- coefficient
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 45
- 238000005457 optimization Methods 0.000 title claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 42
- 238000001914 filtration Methods 0.000 claims abstract description 17
- 238000012706 support-vector machine Methods 0.000 claims abstract description 16
- 238000000605 extraction Methods 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000012549 training Methods 0.000 claims description 15
- 230000003595 spectral effect Effects 0.000 claims description 11
- 238000012360 testing method Methods 0.000 claims description 10
- 238000002474 experimental method Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000003776 cleavage reaction Methods 0.000 claims description 3
- 230000000193 eyeblink Effects 0.000 claims description 3
- 238000009408 flooring Methods 0.000 claims description 3
- 230000007017 scission Effects 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 description 6
- 238000013145 classification model Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012880 independent component analysis Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 210000001061 forehead Anatomy 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- 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
- 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/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/253—Fusion techniques of extracted features
-
- 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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Epidemiology (AREA)
- Signal Processing (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Databases & Information Systems (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于多通道多维特征优化的脑电眨眼伪迹检测方法，首先对多通道脑电EEG信号进行滤波处理以及类别划分；然后对处理后的多通道脑电EEG信号进行16维特征提取；针对提取到的16维特征，采用方差过滤式方法做特征选择；再通过支持向量机算法构建针对眨眼伪迹和非眨眼伪迹样本数据集的伪迹检测模型；最后通过建立的伪迹检测模型进行脑电眨眼伪迹检测。本发明既可以克服临床脑电EEG信号多通道数据损坏问题，又可以解决现有模型复杂性高实时性低的问题，还能够实现眨眼伪迹的精准自动化检测。
Description
技术领域
本发明属于脑电信号处理及智慧医疗领域，涉及一种针对实时多通道强噪声脑电信号的基于支持向量机并融合多通道多维特征表示(FE)与特征优化(FO)的眨眼伪迹检测方法，本发明是基于多通道脑电EEG数据的眨眼伪迹自动化检测方法。
背景技术
脑电图在患者健康检测与诊断等方面应用广泛，然而眨眼伪迹会严重污染采集的脑电信号，并会曲解潜在的神经信号信息而严重阻碍脑电信号后续的分析研究和应用。有效的眨眼伪迹检测算法对脑科学研究和临床应用都具有重大的理论和实际意义。传统的眨眼伪迹检测模型通常是利用多通道脑电EEG信号借助独立成分分析(ICA)方法来构造，其建模的前提往往都假设脑电数据所有通道是干净无其他外部因素干扰、多通道数据线性混合、伪迹成分是非高斯源的。而在实际临床应用中，脑电信号并不满足这些基本假设，并且由于患者的年龄和个体差异导致很多通道数据噪声极大，无法反映真实脑电情况。针对单个个体尤其是儿童来说，多通道脑电信号采集会受到各种因素干扰，多数通道无法使用。并且还会由于年龄差异等原因导致针对个体化的特征失效。传统的基于单一特征的检测算法无法解决这类问题。针对这个问题，本发明提出了一种针对实时多通道强噪声脑电信号的融合多通道多维特征表示(FE)与特征优化(FO)的眨眼伪迹检测方法；同时，针对传统眨眼伪迹检测算法复杂性高、非自动化的弊端，利用支持向量机(SVM)算法实现脑电信号中眨眼伪迹的自动化精准检测。
发明内容
针对现有技术中存在的不足，本发明提供一种基于多通道多维特征优化的脑电眨眼伪迹检测方法。
本发明针对脑电信号中包含大量眨眼伪迹并主要出现在前额通道以及多数通道噪声信号大这一特点，可以利用少数通道数据处理该通道数据损坏问题，提出了可以在少数通道下实现的多维特征表示(FE)与特征优化(FO)方法，并结合支持向量机(SVM)算法构建了眨眼伪迹和非眨眼伪迹信号的精准分类模型。
本发明的技术方案主要包括如下步骤：
步骤1、对多通道脑电EEG信号进行滤波处理以及类别划分。
步骤2、对步骤1处理后的多通道脑电EEG信号进行16维特征提取。
步骤3、针对提取到的16维特征，采用方差过滤式方法做特征选择。
步骤4、通过支持向量机算法构建针对眨眼伪迹和非眨眼伪迹样本数据集的伪迹检测模型。
步骤5、通过建立的伪迹检测模型进行脑电眨眼伪迹检测。
所述步骤1，具体步骤如下：
将21通道、1000Hz的原始脑电EEG信号进行滤波处理。每通道分别经过50Hz陷波滤波器和0.5-30Hz带通滤波器，得到滤波信号；由于眨眼伪迹电位信号持续时间为0.5s，故采用滑动窗口方法进行数据切割，窗口大小设置为1s，步长为0.5s。由此得到切割实验信号。然后将每个切割信号进行类别划分。将包含眨眼伪迹的信号记为eyeblink类，类别标签为1；将不包含眨眼伪迹的信号记为normal类，类别标签为0，并从两类样本中各取出20％作为测试样本，剩下的80％为训练样本。
所述步骤2，具体步骤如下：
2-1计算训练样本的通道1(Fp1通道)和通道2(Fp2通道)的极差均值MAD(mean ofamplitude difference)、平均功率均值MAP(mean of average power)、方差均值MV(meanof variance)、峰态系数均值MKC(mean of kurtosis coefficient)、偏态系数均值MSC(mean of skewness coefficient)、极值点数均值MEN(mean of extremum number)、分形维数均值MFD(mean of fractal dimension)、相关系数CC(correlation coeffcient)、样本熵均值MSE(mean of sample entropy)、多尺度熵均值MME(mean of multiscaleentrop)以及平滑非线性能量算子峰态系数均值MNKC(mean of smooth nonlinear energyoperator kurtosis coefficient)、偏态系数均值MNSC(mean of smooth nonlinearenergy operator skewness coefficient)、相关系数NCC(smooth nonlinear energyoperator correlation coefficient)，并计算通道1(Fp1通道)、通道2(Fp2通道)、通道5(C3通道)和通道6(C4通道)的平均功率差值APD(average power difference)，通道1(Fp1通道)、通道2(Fp2通道)、通道9(O1通道)和通道10(O2通道)的相关系数差值CCD(correlation coefficient difference)以及通道1(Fp1通道)、2(Fp2通道)、3(F3通道)、4(F4通道)四个通道的功率谱密度均值与通道7(P3通道)、8(P4通道)、9(O1通道)、10(O2通道)四个通道的功率谱密度均值的的差值均值PSDD(mean value of power spectraldensity difference)。以上共16维特征。
假设样本数据为S21*1000，Si(i＝1，2，...，21)表示第i个通道的所有数据，N表示数据长度。
2-2通过如下公式计算极差均值MAD：
式中max()代表取最大值，min()代表取最小值。
2-3通过如下公式计算方差均值MV：
式中μi表示第i个通道所有数据的均值，Si(t)表示第i个通道的第t个数据。
2-4通过如下公式计算平均功率均值MAP：
2-5通过如下公式计算峰态系数均值MKC：
其中σi表示第i个通道所有数据的标准差。
2-6通过如下公式计算偏态系数均值MSC：
2-7通过如下方法计算极值点数均值MEN：
首先计算Si(t)，i＝1，2的极大值以及极小值点，然后统计极大值点个数与极小值点个数之和，分别记为eni，i＝1，2，则
2-8通过如下公式与方法计算分形维数均值MFD：
其中FD()表示分形维数值，ε表示计盒维数计算时的格子边长，G()表示信号数据的格子数累计，通过上式求得前两通道的分形维数值并求平均得到MFD。
2-9通过如下公式相关系数CC：
其中var()表示方差计算，Cov()表示协方差计算。
2-10通过如下公式与方法计算样本熵均值MSE：
式中SE()表示信号样本熵，m表示重构信号维数，δ表示时间延迟，r表示预定义的阈值参数，n()表示匹配m维向量对的数目。通过上式求得前两通道的样本熵值并求平均得到MSE。
2-11通过如下公式与方法计算多尺度熵均值MME：
ME(S，m，τ，r)＝SE(S′，m，δ＝τ，r)
式中ME()表示多尺度熵值，
通过上式求得前两通道的多尺度熵值并计算平均得到MME。
2-12通过如下公式计算第i个通道Si的所有数据的平滑非线性能量算子数据，记为ESi：
ESi(t)＝ω(t)*(Si(t)*Si(t)-Si(t-1)*Si(t+1))
式中ω(t)表示窗函数，ESi(t)表示ESi的第t个数据。然后分别计算平滑非线性能量算子峰态系数均值MNKC、平滑非线性能量算子偏态系数均值MNSC以及平滑非线性能量算子相关系数NCC，公式如下：。
式中Eμi和Eσi分别表示ESi的均值和标准差。
2-13为了弥补多通道脑电信号的空间信息，通过如下公式与方法计算平均功率差值APD：
APD＝2*(P(S1)+P(S2))-(P(S5)+P(S6)+P(S7)+P(S8))
式中P()为计算平均功率值，其计算公式为：
2-14通过如下公式与方法计算相关系数差值CCD：
CCD＝2*CC(S1，S2)-CC(S1，S9)-CC(S2，S10)
式中CC()为利用步骤2-9计算相关系数值。
2-15通过如下公式与方法计算功率谱密度差值PSDD：
式中PSD()代表通道信号的功率谱密度,mean()代表均值计算。
2-16假设两类训练样本所有数目为T，对所有训练样本进行步骤2-2至2-15的16维特征提取，并构成T*16维的特征向量FVs。
所述步骤3，具体步骤如下：
针对特征向量FVs,采用方差过滤式方法做特征选择。
3-1通过如下公式计算特征向量FVs每一维度方差Vj(j＝1，2，...，16)。
其中L为每一维度所有数据个数，vμj为第j个维度所有数据均值，Vj(t)为第j个维度的第t个数据。
3-2将得到的方差结果进行降序排序，输出排序后的方差值对应的特征维度索引。
3-3根据3-2得到的特征维度索引结果，从前往后选取数目为M的特征维度索引对应的特征作为新的T*M维的特征向量NFVs。
所述步骤4，具体步骤如下：
4-1利用步骤3得到的特征向量NFVs结合两个类别的类别标签和支持向量机算法进行模型训练得到分类器模型。
4-2对测试样本进行特征提取，提取的特征为步骤3特征选择后的M维特征，然后得到测试样本特征向量TFVs，并利用分类器模型对TFVs进行逐个预测分类，计算检测结果。
4-3根据检测结果，通过调节分类器模型的核函数类型、核函数参数构建出最优分类器模型，即伪迹检测模型。
本发明有益效果如下
本发明针对个人尤其是儿童脑电信号多通道数据极易损坏以及个体间和同一个体不同年龄间信号差异性大的问题，可以利用少数通道数据，融合多维特征表示(FE)与特征优化(FO)，并结合支持向量机算法构建眨眼伪迹和非眨眼伪迹信号的精准分类模型，既可以克服临床脑电EEG信号多通道数据损坏问题，又可以解决现有模型复杂性高实时性低的问题，还能够实现眨眼伪迹的精准自动化检测。
附图说明
图1本发明的方法流程图；
图2本发明多维特征提取与选择示意图。
具体实施方式
下面结合附图和具体实施方式对本发明作详细说明。
如图1和2所示，基于支持向量机的眨眼伪迹检测方法，多维特征表示(FE)与特征优化(FO)方法的实现步骤在发明内容内已有详细的介绍，即本发明的技术方案主要包括如下步骤：
步骤1、对多通道脑电EEG信号进行滤波处理以及类别划分。
步骤2、对步骤1处理后的多通道脑电EEG信号进行16维特征提取。
步骤3、针对提取到的16维特征，采用方差过滤式方法做特征选择。
步骤4、通过支持向量机算法构建针对眨眼伪迹和非眨眼伪迹样本数据集的伪迹检测模型。
步骤5、通过建立的伪迹检测模型进行脑电眨眼伪迹检测。
所述步骤1，具体步骤如下：
将21通道、1000Hz的原始脑电EEG信号进行滤波处理。每通道分别经过50Hz陷波滤波器和0.5-30Hz带通滤波器，得到滤波信号；由于眨眼伪迹电位信号持续时间为0.5s，故采用滑动窗口方法进行数据切割，窗口大小设置为1s，步长为0.5s。由此得到切割实验信号。然后将每个切割信号进行类别划分。将包含眨眼伪迹的信号记为eyeblink类，类别标签为1；将不包含眨眼伪迹的信号记为normal类，类别标签为0，并从两类样本中各取出20％作为测试样本，剩下的80％为训练样本。
所述步骤2，具体步骤如下：
2-1计算训练样本的通道1(Fp1通道)和通道2(Fp2通道)的极差均值MAD(mean ofamplitude difference)、平均功率均值MAP(mean of average power)、方差均值MV(meanof variance)、峰态系数均值MKC(mean of kurtosis coefficient)、偏态系数均值MSC(mean of skewness coefficient)、极值点数均值MEN(mean of extremum number)、分形维数均值MFD(mean of fractal dimension)、相关系数CC(correlation coeffcient)、样本熵均值MSE(mean of sample entropy)、多尺度熵均值MME(mean of multiscaleentrop)以及平滑非线性能量算子峰态系数均值MNKC(mean of smooth nonlinear energyoperator kurtosis coefficient)、偏态系数均值MNSC(mean of smooth nonlinearenergy operator skewness coefficient)、相关系数NCC(smooth nonlinear energyoperator correlation coefficient)，并计算通道1(Fp1通道)、通道2(Fp2通道)、通道5(C3通道)和通道6(C4通道)的平均功率差值APD(average power difference)，通道1(Fp1通道)、通道2(Fp2通道)、通道9(O1通道)和通道10(O2通道)的相关系数差值CCD(correlation coefficient difference)以及通道1(Fp1通道)、2(Fp2通道)、3(F3通道)、4(F4通道)四个通道的功率谱密度均值与通道7(P3通道)、8(P4通道)、9(O1通道)、10(O2通道)四个通道的功率谱密度均值的的差值均值PSDD(mean value of power spectraldensity difference)。以上共16维特征。
假设样本数据为S21*1000，Si(i＝1，2，...，21)表示第i个通道的所有数据，N表示数据长度。
2-2通过如下公式计算极差均值MAD：
式中max()代表取最大值，min()代表取最小值。
2-3通过如下公式计算方差均值MV：
式中μi表示第i个通道所有数据的均值，Si(t)表示第i个通道的第t个数据。
2-4通过如下公式计算平均功率均值MAP：
2-5通过如下公式计算峰态系数均值MKC：
其中σi表示第i个通道所有数据的标准差。
2-6通过如下公式计算偏态系数均值MSC：
2-7通过如下方法计算极值点数均值MEN：
首先计算Si(t)，i＝1,2的极大值以及极小值点，然后统计极大值点个数与极小值点个数之和，分别记为eni，i＝1，2，则
2-8通过如下公式与方法计算分形维数均值MFD：
其中FD()表示分形维数值，ε表示计盒维数计算时的格子边长，G()表示信号数据的格子数累计，通过上式求得前两通道的分形维数值并求平均得到MFD。
2-9通过如下公式相关系数CC：
其中var()表示方差计算，Cov()表示协方差计算。
2-10通过如下公式与方法计算样本熵均值MSE：
式中SE()表示信号样本熵，m表示重构信号维数，δ表示时间延迟，r表示预定义的阈值参数，n()表示匹配m维向量对的数目。通过上式求得前两通道的样本熵值并求平均得到MSE。
2-11通过如下公式与方法计算多尺度熵均值MME：
ME(S，m，τ，r)＝SE(S′，m，δ＝τ，r)
式中ME()表示多尺度熵值，
通过上式求得前两通道的多尺度熵值并计算平均得到MME。
2-12通过如下公式计算第i个通道Si的所有数据的平滑非线性能量算子数据，记为ESi：
ESi(t)＝ω(t)*(Si(t)*Si(t)-Si(t-1)*Si(t+1))
式中ω(t)表示窗函数，ESi(t)表示ESi的第t个数据。然后分别计算平滑非线性能量算子峰态系数均值MNKC、平滑非线性能量算子偏态系数均值MNSC以及平滑非线性能量算子相关系数NCC，公式如下：。
式中Eμi和Eσi分别表示ESi的均值和标准差。
2-13为了弥补多通道脑电信号的空间信息，通过如下公式与方法计算平均功率差值APD：
APD＝2*(P(S1)+P(S2))-(P(S5)+P(S6)+P(S7)+P(S8))
式中P()为计算平均功率值，其计算公式为：
2-14通过如下公式与方法计算相关系数差值CCD：
CCD＝2*CC(S1，S2)-CC(S1，S9)-CC(S2，S10)
式中CC()为利用步骤2-9计算相关系数值。
2-15通过如下公式与方法计算功率谱密度差值PSDD：
式中PSD()代表通道信号的功率谱密度，mean()代表均值计算。
2-16假设两类训练样本所有数目为T，对所有训练样本进行步骤2-2至2-15的16维特征提取，并构成T*16维的特征向量FVs。
所述步骤3，具体步骤如下：
针对特征向量FVs，采用方差过滤式方法做特征选择。
3-1通过如下公式计算特征向量FVs每一维度方差Vj(j＝1，2，...，16)。
其中L为每一维度所有数据个数，vμj为第j个维度所有数据均值，Vj(t)为第j个维度的第t个数据。
3-2将得到的方差结果进行降序排序，输出排序后的方差值对应的特征维度索引。
3-3根据3-2得到的特征维度索引结果，从前往后选取数目为M的特征维度索引对应的特征作为新的T*M维的特征向量NFVs。
所述步骤4，具体步骤如下：
4-1利用步骤3得到的特征向量NFVs结合两个类别的类别标签和支持向量机算法进行模型训练得到分类器模型。
4-2对测试样本进行特征提取，提取的特征为步骤3特征选择后的M维特征，然后得到测试样本特征向量TFVs，并利用分类器模型对TFVs进行逐个预测分类，计算检测结果。
4-3根据检测结果，通过调节分类器模型的核函数类型、核函数参数构建出最优分类器模型，即伪迹检测模型。
为了达到更好的眨眼伪迹和非眨眼伪迹信号的精准检测效果，以下将从实际应用时参数的选择与设计方面展开介绍，以作为该发明用于其他应用的参考:
本方法在处理多通道脑电EEG信号时，设置每一段输入数据为1秒，数据的重叠率设为50％。针对临床上个体差异导致个体化特征失效以及儿童信号采集噪声大、多通道数据损坏等问题，本发明利用信号多维特征提取以及特征选择算法处理个体差异问题。步骤2-2至2-12提取的多维特征包含了脑电信号的时域、频域等最主要信息，而步骤2-12至2-15提取的特征则是为了保证多通道脑电信息的空间性以及完全性。
在特征选择过程中，保留方差排序在前8位的特征。最后采用的支持向量机算法设置了rbf核函数以实现非线性分类，参数sigma设置为3，参数C设置为2。
为了真实测试本发明对眨眼伪迹信号的检测效果，在浙江大学医学院附属儿童医院真实病患EEG数据上与多种目前主流检测算法进行了对比实验：
实验数据为21通道，采样频率为1000Hz，数据共分为5个不同患病个体，平均数据长度为18分钟。LeorShoker等人在2005年提出的BSS-SVM算法在数据集上的平均灵敏度(sensitivity)为61.47％，平均特异度(specificity)为70.98％；ChongYehSai等人在2018年提出的Wavelet-ICA-SVM算法在数据集上的平均灵敏度为67.66％,平均特异度为75.78％；RajdeepGhosh等人在2018年提出的SVM-Autoencoder算法在数据集上的平均灵敏度为91.34％，平均特异度为96.80％；本发明在数据集上的平均灵敏度为96.48％，平均特异度为98.62％。本发明相对于上述三种对比算法平均灵敏度分别提高了35.01％、28.82％、5.14％，平均特异度分别提高了27.64％、22.84％、1.82％。这充分证明了本发明在真实数据上的相较于目前眨眼伪迹检测模型的优异性能。
本发明的提出可用于解决临床上个体差异导致个体化特征失效以及儿童信号采集噪声大、多通道数据损坏的问题，通过信号多维特征表示(FE)与特征优化(FO)算法得到特征向量，结合支持向量机分类模型能够达到精准的眨眼伪迹与非伪迹信号的精准检测。
Claims (5)
1.基于多通道多维特征优化的脑电眨眼伪迹检测方法，其特征在于，步骤如下：
步骤1、对多通道脑电EEG信号进行滤波处理以及类别划分；
步骤2、对步骤1处理后的多通道脑电EEG信号进行16维特征提取；
步骤3、针对提取到的16维特征，采用方差过滤式方法做特征选择；
步骤4、通过支持向量机算法构建针对眨眼伪迹和非眨眼伪迹样本数据集的伪迹检测模型；
步骤5、通过建立的伪迹检测模型进行脑电眨眼伪迹检测。
2.根据权利要求1所述的基于多通道多维特征优化的脑电眨眼伪迹检测方法，其特征在于，所述步骤1，具体步骤如下：
将21通道、1000Hz的原始脑电EEG信号进行滤波处理；每通道分别经过50Hz陷波滤波器和0.5-30Hz带通滤波器，得到滤波信号；由于眨眼伪迹电位信号持续时间为0.5s，故采用滑动窗口方法进行数据切割，窗口大小设置为1s，步长为0.5s；由此得到切割实验信号；然后将每个切割信号进行类别划分；将包含眨眼伪迹的信号记为eyeblink类，类别标签为1；将不包含眨眼伪迹的信号记为normal类，类别标签为0，并从两类样本中各取出20％作为测试样本，剩下的80％为训练样本。
3.根据权利要求2所述的基于多通道多维特征优化的脑电眨眼伪迹检测方法，其特征在于，所述步骤2，具体步骤如下：
2-1计算训练样本的通道1(Fp1通道)和通道2(Fp2通道)的极差均值MAD(mean ofamplitude difference)、平均功率均值MAP(mean of average power)、方差均值MV(meanof variance)、峰态系数均值MKC(mean of kurtosis coefficient)、偏态系数均值MSC(mean of skewness coefficient)、极值点数均值MEN(mean of extremum number)、分形维数均值MFD(mean of fractal dimension)、相关系数CC(correlation coeffcient)、样本熵均值MSE(mean of sample entropy)、多尺度熵均值MME(mean of multiscaleentrop)以及平滑非线性能量算子峰态系数均值MNKC(mean of smooth nonlinear energyoperator kurtosis coefficient)、偏态系数均值MNSC(mean of smooth nonlinearenergy operator skewness coefficient)、相关系数NCC(smooth nonlinear energyoperator correlation coefficient)，并计算通道1(Fp1通道)、通道2(Fp2通道)、通道5(C3通道)和通道6(C4通道)的平均功率差值APD(average power difference)，通道1(Fp1通道)、通道2(Fp2通道)、通道9(O1通道)和通道10(O2通道)的相关系数差值CCD(correlation coefficient difference)以及通道1(Fp1通道)、2(Fp2通道)、3(F3通道)、4(F4通道)四个通道的功率谱密度均值与通道7(P3通道)、8(P4通道)、9(O1通道)、10(O2通道)四个通道的功率谱密度均值的的差值均值PSDD(mean value of power spectraldensity difference)；以上共16维特征；
假设样本数据为S21*1000，Si(i＝1，2，...，21)表示第i个通道的所有数据，N表示数据长度；
2-2通过如下公式计算极差均值MAD：
式中max()代表取最大值，min()代表取最小值；
2-3通过如下公式计算方差均值MV：
式中μi表示第i个通道所有数据的均值，Si(t)表示第i个通道的第t个数据；
2-4通过如下公式计算平均功率均值MAP：
2-5通过如下公式计算峰态系数均值MKC：
其中σi表示第i个通道所有数据的标准差；
2-6通过如下公式计算偏态系数均值MSC：
2-7通过如下方法计算极值点数均值MEN：
首先计算Si(t)，i＝1，2的极大值以及极小值点，然后统计极大值点个数与极小值点个数之和，分别记为eni，i＝1，2，则
2-8通过如下公式与方法计算分形维数均值MFD：
其中FD()表示分形维数值，ε表示计盒维数计算时的格子边长，G()表示信号数据的格子数累计，通过上式求得前两通道的分形维数值并求平均得到MFD；
2-9通过如下公式相关系数CC：
其中var()表示方差计算，Cov()表示协方差计算；
2-10通过如下公式与方法计算样本熵均值MSE：
式中SE()表示信号样本熵，m表示重构信号维数，δ表示时间延迟，r表示预定义的阈值参数，n()表示匹配m维向量对的数目；通过上式求得前两通道的样本熵值并求平均得到MSE；
2-11通过如下公式与方法计算多尺度熵均值MME：
ME(S，m，τ，r)＝SE(S′，m，δ＝τ，r)
式中ME()表示多尺度熵值，
通过上式求得前两通道的多尺度熵值并计算平均得到MME；
2-12通过如下公式计算第i个通道Si的所有数据的平滑非线性能量算子数据，记为ESi：
ESi(t)＝ω(t)*(Si(t)*Si(t)-Si(t-1)*Si(t+1))
式中ω(t)表示窗函数，ESi(t)表示ESi的第t个数据；然后分别计算平滑非线性能量算子峰态系数均值MNKC、平滑非线性能量算子偏态系数均值MNSC以及平滑非线性能量算子相关系数NCC，公式如下：；
式中Eμi和Eσi分别表示ESi的均值和标准差；
2-13为了弥补多通道脑电信号的空间信息，通过如下公式与方法计算平均功率差值APD：
APD＝2*(P(S1)+P(S2))-(P(S5)+P(S6)+P(S7)+P(S8))
式中P()为计算平均功率值，其计算公式为：
2-14通过如下公式与方法计算相关系数差值CCD：
CCD＝2*CC(S1，S2)-CC(S1，S9)-CC(S2，S10)
式中CC()为利用步骤2-9计算相关系数值；
2-15通过如下公式与方法计算功率谱密度差值PSDD：
式中PSD()代表通道信号的功率谱密度，mean()代表均值计算；
2-16假设两类训练样本所有数目为T，对所有训练样本进行步骤2-2至2-15的16维特征提取，并构成T*16维的特征向量FVs。
4.根据权利要求3所述的基于多通道多维特征优化的脑电眨眼伪迹检测方法，其特征在于，所述步骤3，具体步骤如下：
针对特征向量FVs，采用方差过滤式方法做特征选择；
3-1通过如下公式计算特征向量FVs每一维度方差Vj(j＝1，2，...，16)；
其中L为每一维度所有数据个数，vμj为第j个维度所有数据均值，Vj(t)为第j个维度的第t个数据；
3-2将得到的方差结果进行降序排序，输出排序后的方差值对应的特征维度索引；
3-3根据3-2得到的特征维度索引结果，从前往后选取数目为M的特征维度索引对应的特征作为新的T*M维的特征向量NFVs。
5.根据权利要求4所述的基于多通道多维特征优化的脑电眨眼伪迹检测方法，其特征在于，所述步骤4，具体步骤如下：
4-1利用步骤3得到的特征向量NFVs结合两个类别的类别标签和支持向量机算法进行模型训练得到分类器模型；
4-2对测试样本进行特征提取，提取的特征为步骤3特征选择后的M维特征，然后得到测试样本特征向量TFVs，并利用分类器模型对TFVs进行逐个预测分类，计算检测结果；
4-3根据检测结果，通过调节分类器模型的核函数类型、核函数参数构建出最优分类器模型，即伪迹检测模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011485932.1A CN112464902B (zh) | 2020-12-16 | 2020-12-16 | 基于多通道多维特征优化的脑电眨眼伪迹检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011485932.1A CN112464902B (zh) | 2020-12-16 | 2020-12-16 | 基于多通道多维特征优化的脑电眨眼伪迹检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112464902A true CN112464902A (zh) | 2021-03-09 |
| CN112464902B CN112464902B (zh) | 2024-02-09 |
Family
ID=74803017
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011485932.1A Active CN112464902B (zh) | 2020-12-16 | 2020-12-16 | 基于多通道多维特征优化的脑电眨眼伪迹检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112464902B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113662558A (zh) * | 2021-08-19 | 2021-11-19 | 杭州电子科技大学 | 一种区分脑电眨眼伪迹与额极痫样放电的智能分类方法 |
| CN113835068A (zh) * | 2021-09-22 | 2021-12-24 | 南京信息工程大学 | 一种基于独立成分分析的盲源分离实时抗主瓣干扰方法 |
| WO2023116263A1 (zh) * | 2021-12-24 | 2023-06-29 | 南京邮电大学 | 一种基于脑电图的眨眼咬牙判断方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104586387A (zh) * | 2015-01-19 | 2015-05-06 | 秦皇岛市惠斯安普医学系统有限公司 | 一种时、频、空域多参数脑电特征提取与融合方法 |
| CN107260166A (zh) * | 2017-05-26 | 2017-10-20 | 昆明理工大学 | 一种实用化在线脑电伪迹剔除方法 |
| CN109924990A (zh) * | 2019-03-27 | 2019-06-25 | 兰州大学 | 一种基于emd算法的脑电信号抑郁症识别系统 |
| CN111134667A (zh) * | 2020-01-19 | 2020-05-12 | 中国人民解放军战略支援部队信息工程大学 | 基于脑电信号的时间迁移情绪识别方法及系统 |
| CN111956217A (zh) * | 2020-07-15 | 2020-11-20 | 山东师范大学 | 一种面向实时脑电信号的眨眼伪迹识别方法及系统 |
-
2020
- 2020-12-16 CN CN202011485932.1A patent/CN112464902B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104586387A (zh) * | 2015-01-19 | 2015-05-06 | 秦皇岛市惠斯安普医学系统有限公司 | 一种时、频、空域多参数脑电特征提取与融合方法 |
| CN107260166A (zh) * | 2017-05-26 | 2017-10-20 | 昆明理工大学 | 一种实用化在线脑电伪迹剔除方法 |
| CN109924990A (zh) * | 2019-03-27 | 2019-06-25 | 兰州大学 | 一种基于emd算法的脑电信号抑郁症识别系统 |
| CN111134667A (zh) * | 2020-01-19 | 2020-05-12 | 中国人民解放军战略支援部队信息工程大学 | 基于脑电信号的时间迁移情绪识别方法及系统 |
| CN111956217A (zh) * | 2020-07-15 | 2020-11-20 | 山东师范大学 | 一种面向实时脑电信号的眨眼伪迹识别方法及系统 |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113662558A (zh) * | 2021-08-19 | 2021-11-19 | 杭州电子科技大学 | 一种区分脑电眨眼伪迹与额极痫样放电的智能分类方法 |
| CN113662558B (zh) * | 2021-08-19 | 2024-01-30 | 杭州电子科技大学 | 一种区分脑电眨眼伪迹与额极痫样放电的智能分类方法 |
| CN113835068A (zh) * | 2021-09-22 | 2021-12-24 | 南京信息工程大学 | 一种基于独立成分分析的盲源分离实时抗主瓣干扰方法 |
| CN113835068B (zh) * | 2021-09-22 | 2023-06-20 | 南京信息工程大学 | 一种基于独立成分分析的盲源分离实时抗主瓣干扰方法 |
| WO2023116263A1 (zh) * | 2021-12-24 | 2023-06-29 | 南京邮电大学 | 一种基于脑电图的眨眼咬牙判断方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112464902B (zh) | 2024-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112464902B (zh) | 基于多通道多维特征优化的脑电眨眼伪迹检测方法 | |
| Sangaiah et al. | An intelligent learning approach for improving ECG signal classification and arrhythmia analysis | |
| CN104586387B (zh) | 一种时、频、空域多参数脑电特征提取与融合方法 | |
| CN110269609B (zh) | 基于单通道的脑电信号中眼电伪迹分离方法 | |
| CN107024987A (zh) | 一种基于eeg的实时人脑注意力测试和训练系统 | |
| CN114886388B (zh) | 一种麻醉深度监测过程中脑电信号质量的评估方法及装置 | |
| CN107260166A (zh) | 一种实用化在线脑电伪迹剔除方法 | |
| CN106236080B (zh) | 基于多通道的脑电信号中肌电噪声的消除方法 | |
| CN106485208B (zh) | 单通道脑电信号中眼电干扰的自动去除方法 | |
| CN110059564B (zh) | 基于功率谱密度和互相关熵谱密度融合的特征提取方法 | |
| CN104305993A (zh) | 一种基于格兰杰因果性的脑电源定位方法 | |
| CN113850192A (zh) | 基于emd-csp脑电特征学习与智能融合的眨眼伪迹检测方法 | |
| CN104905786A (zh) | 一种心电伪迹在线去除算法 | |
| Ai et al. | Multi-features fusion diagnosis of tremor based on artificial neural network and D–S evidence theory | |
| Liu et al. | Identification of anisomerous motor imagery EEG signals based on complex algorithms | |
| Shoker et al. | Removal of eye blinking artifact from the electro-encephalogram, incorporating a new constrained blind source separation algorithm | |
| Navarro et al. | A comparison of time, frequency and ICA based features and five classifiers for wrist movement classification in EEG signals | |
| Ahmed et al. | Effective hybrid method for the detection and rejection of electrooculogram (EOG) and power line noise artefacts from electroencephalogram (EEG) mixtures | |
| CN105395185A (zh) | 基于多通道柔性融合的生理信号检测系统 | |
| Priyadharsini et al. | An Efficient method for the removal of ECG artifact from measured EEG Signal using PSO algorithm | |
| Li et al. | An improved MAMA-EMD for the automatic removal of EOG artifacts | |
| Yun-Mei et al. | The abnormal detection of electroencephalogram with three-dimensional deep convolutional neural networks | |
| Bhimraj et al. | Autonomous noise removal from EEG signals using independent component analysis | |
| CN111671421A (zh) | 一种基于脑电图的儿童需求感知方法 | |
| Turnip et al. | Utilization of EEG-SSVEP method and ANFIS classifier for controlling electronic wheelchair |
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 |