CN112464902B - 基于多通道多维特征优化的脑电眨眼伪迹检测方法 - Google Patents

基于多通道多维特征优化的脑电眨眼伪迹检测方法 Download PDF

Info

Publication number
CN112464902B
CN112464902B CN202011485932.1A CN202011485932A CN112464902B CN 112464902 B CN112464902 B CN 112464902B CN 202011485932 A CN202011485932 A CN 202011485932A CN 112464902 B CN112464902 B CN 112464902B
Authority
CN
China
Prior art keywords
mean
channel
following formula
data
blink
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
Application number
CN202011485932.1A
Other languages
English (en)
Other versions
CN112464902A (zh
Inventor
曹九稳
王建辉
胡丁寒
蒋铁甲
高峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN202011485932.1A priority Critical patent/CN112464902B/zh
Publication of CN112464902A publication Critical patent/CN112464902A/zh
Application granted granted Critical
Publication of CN112464902B publication Critical patent/CN112464902B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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 (3)

1.基于多通道多维特征优化的脑电眨眼伪迹检测方法，其特征在于，步骤如下：
步骤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和通道2的极差均值MAD、平均功率均值MAP、方差均值MV、峰态系数均值MKC、偏态系数均值MSC、极值点数均值MEN、分形维数均值MFD、相关系数CC、样本熵均值MSE、多尺度熵均值MME以及平滑非线性能量算子峰态系数均值MNKC、偏态系数均值MNSC、相关系数NCC，并计算通道1、通道2、通道5和通道6的平均功率差值APD，通道1、通道2、通道9和通道10的相关系数差值CCD以及通道1、2、3、4四个通道的功率谱密度均值与通道7、8、9、10四个通道的功率谱密度均值的的差值均值PSDD；共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()表示多尺度熵值，1≤j≤N-τ+1，通过上式求得前两通道的多尺度熵值并计算平均得到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。
2.根据权利要求1所述的基于多通道多维特征优化的脑电眨眼伪迹检测方法，其特征在于，所述步骤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。
3.根据权利要求2所述的基于多通道多维特征优化的脑电眨眼伪迹检测方法，其特征在于，所述步骤4，具体步骤如下：
4-1利用步骤3得到的特征向量NFVs结合两个类别的类别标签和支持向量机算法进行模型训练得到分类器模型；
4-2对测试样本进行特征提取，提取的特征为步骤3特征选择后的M维特征，然后得到测试样本特征向量TFVs，并利用分类器模型对TFVs进行逐个预测分类，计算检测结果；
4-3根据检测结果，通过调节分类器模型的核函数类型、核函数参数构建出最优分类器模型，即伪迹检测模型。
CN202011485932.1A 2020-12-16 2020-12-16 基于多通道多维特征优化的脑电眨眼伪迹检测方法 Active CN112464902B (zh)

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 CN112464902A (zh) 2021-03-09
CN112464902B true 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)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113662558B (zh) * 2021-08-19 2024-01-30 杭州电子科技大学 一种区分脑电眨眼伪迹与额极痫样放电的智能分类方法
CN113835068B (zh) * 2021-09-22 2023-06-20 南京信息工程大学 一种基于独立成分分析的盲源分离实时抗主瓣干扰方法
CN114431879B (zh) * 2021-12-24 2024-04-16 南京邮电大学 一种基于脑电图的眨眼咬牙判断方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
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 山东师范大学 一种面向实时脑电信号的眨眼伪迹识别方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
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 山东师范大学 一种面向实时脑电信号的眨眼伪迹识别方法及系统

Also Published As

Publication number Publication date
CN112464902A (zh) 2021-03-09

Similar Documents

Publication Publication Date Title
CN112464902B (zh) 基于多通道多维特征优化的脑电眨眼伪迹检测方法
CN104586387B (zh) 一种时、频、空域多参数脑电特征提取与融合方法
Sangaiah et al. An intelligent learning approach for improving ECG signal classification and arrhythmia analysis
US11529103B2 (en) Artificial intelligence self-learning-based automatic electrocardiography analysis method and apparatus
CN107951485B (zh) 基于人工智能自学习的动态心电图分析方法和装置
CN108577835B (zh) 一种基于微状态的脑功能网络构建方法
CN103971124B (zh) 一种基于相位同步的多类别运动想象脑电信号分类方法
CN108511055B (zh) 基于分类器融合及诊断规则的室性早搏识别系统及方法
CN110269609B (zh) 基于单通道的脑电信号中眼电伪迹分离方法
CN107909117A (zh) 一种基于脑功能网络特征对早晚期轻度认知障碍的分类方法及装置
CN110353673B (zh) 一种基于标准互信息的脑电通道选择方法
CN103092971B (zh) 一种用于脑机接口中的分类方法
CN114533086B (zh) 一种基于空域特征时频变换的运动想象脑电解码方法
Chen et al. A Kalman filtering and nonlinear penalty regression approach for noninvasive anemia detection with palpebral conjunctiva images
CN114469124B (zh) 一种运动过程中异常心电信号的识别方法
Diery et al. Automated ECG diagnostic P-wave analysis using wavelets
CN106236080A (zh) 基于多通道的脑电信号中肌电噪声的消除方法
WO2022166401A1 (zh) 基于eemd-pca去除eeg信号中运动伪迹的方法及装置
CN113191395A (zh) 一种基于双脑的多层次信息融合的目标检测方法
Gavas et al. Multivariate variational mode decomposition based approach for blink removal from EEG signal
Vijayasankar et al. Correction of blink artifacts from single channel EEG by EMD-IMF thresholding
CN114886388A (zh) 一种麻醉深度监测过程中脑电信号质量的评估方法及装置
CN103720471A (zh) 一种基于因子分析的眼电伪迹去除方法
CN113128384B (zh) 一种基于深度学习的脑卒中康复系统脑机接口软件关键技术方法
CN112438741B (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
GR01 Patent grant
GR01 Patent grant