CN102835955A - 一种无需设定阈值的脑电信号中眼电伪迹自动去除的方法 - Google Patents
一种无需设定阈值的脑电信号中眼电伪迹自动去除的方法 Download PDFInfo
- Publication number
- CN102835955A CN102835955A CN201210331193XA CN201210331193A CN102835955A CN 102835955 A CN102835955 A CN 102835955A CN 201210331193X A CN201210331193X A CN 201210331193XA CN 201210331193 A CN201210331193 A CN 201210331193A CN 102835955 A CN102835955 A CN 102835955A
- Authority
- CN
- China
- Prior art keywords
- isolated component
- formula
- value
- eye
- center
- 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
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明提出一种无需人工设定阈值的脑电信号中眼电伪迹自动去除的方法,属于生物信息技术领域,主要应用于脑电信号的预处理过程中。具体包括:对采集的含有眼电伪迹的脑电信号进行独立分量分解;并且求取每个独立分量的峭度、序列renyi熵和样本熵作为特征向量,进而使用k均值聚类分析的方法自动识别出含有眼电伪迹的独立分量,并将其置零,其余分量不变,对信号进行重构,得到纯净的脑电信号。本发明解决了传统的眼电伪迹去除过程中需要人工对伪迹进行识别、费时费力、工作量大的问题,并且本方法无需人工设定阈值就可以实现自动识别并去除眼电伪迹的目的,弥补了以往方法中设定阈值时需要研究人员具备一定的先验知识、主观性强的不足。
Description
技术领域:
本发明涉及生物信息技术领域,特别涉及脑电信号(Electroenc-
ephalogram,EEG)的预处理技术。具体是涉及一种无需人工设定阈值的脑电信号中眼电伪迹(Ocular Artifact,OA)的自动去除技术。
技术背景:
脑电信号是一种微弱的电生理信号,其中包含了大量的生理与疾病信息,在临床医学、认知心理学等研究领域中发挥着重要作用。头皮表面记录脑电的方式是目前脑电信号研究中所采用的最主要的信号采集方式,该方式记录技术简单、易操作,但信号易受噪声的干扰。脑电干扰给脑电信号的解释与分析造成了很大的困难。因此,有效识别并去除脑电信号中的伪迹,得到纯净的脑电信号,对进一步研究大脑的真实活动具有重要意义。
影响脑电信号的伪迹主要包括眼电伪迹、心电伪迹、肌电伪迹、50Hz工频干扰等。相对于其它伪迹成分而言,眼电伪迹是使EEG信号失真的一种最主要的干扰成分。眼电伪迹的幅值较大,最大可达到100mV,并且眼电信号源距离测量电极较近,对脑电信号造成了很大的干扰。眼电的产生是由于跨膜间存在静息电位,角膜侧为正,巩膜侧为负,这可以看作是有一个电偶极子存在。眼球转动过程中,该偶极子改变了眼睛周围的电场,从而影响头皮电场,眼球运动方向一侧的电极记录到眼电伪迹为正,眼球运动相反方向一侧的电极记录得到的眼电伪迹为负,由于大脑是电的良导体,眼电信号从前额处向后传播,遍历整个头部,放置在头部用于采集脑电信号的电极可以采集到这些信号,造成对脑电信号的干扰。而且,眼电伪迹的频带与脑电信号的频带相混叠,很难使用传统的滤波方法进行去除,因此,脑电信号预处理阶段最重要的工作就是如何有效地去除眼电伪迹。
从脑电信号中有效地去除眼电伪迹一直是生物医学信号处理领域的研究热点,其中一些主要的研究成果有:
1实验控制法:这种方法要求研究人员给予受试者一定的实验指导,例如要求受试者在实验过程中始终保持闭眼状态或尽量避免眨眼和眼球转动以减少眼电伪迹的干扰,这样做虽然可以在一定程度上抑制干扰,但是前者会引入α波的干扰,后者在长时间的实验过程中极易引起疲劳,并且不易做到。
2伪迹剔除法:它的主要思想是将受到眼电伪迹影响脑电信号剔除以获得较纯净的脑电信号,但是这样做的同时剔除了大量的脑电信号,造成了信息的遗漏,并且人工剔除干扰信号的劳动量过大,主观性强,而且不适应在线脑机接口***的需求。
3回归方法:研究表明眼电伪迹在测量电极处的衰减程度与电极到眼部的距离有关,此方法利用这一原理通过计算得到眼电信号在每个头皮电极处的衰减因子,然后将采集到的每导脑电信号减去乘以衰减因子的眼电信号,从而得到伪迹去除后的数据。它的模型为
式中,EOG表示测量得到的眼电信号;EEGs(ε)表示原始脑电信号;EEG(ε)表示伪迹去除后的脑电信号;ε表示相应的测量通道;η表示比例系数,即眼电信号在第ε导处的衰减因子。η值常用回归分析的方法进行估计。这种伪迹去除方法的优点是直观易懂,物理意义明确,缺点是没有考虑到脑电信号的双向传播性,即实际测得的眼电信号中也包含一些脑电成分。因此,回归方法去除伪迹的同时不可避免的去除某些脑电成分。
4小波变换法:小波变换(Wavelet Transform, WT)是数学领域的一个重要分支,近年来越来越受到各个学科领域的极大重视。由于其在时频两域都具有表征信号局部特征的能力和多分辨率分析的特点,因此被誉为“数学显微镜”。小波变换的基本思想是将原始信号通过伸缩和平移后,分解为一系列具有不同时间分辨率、不同频率特性的子带信号,这些子带信号具有良好的时域、频域等局部特征,进而可以实现对信号在时域和频域的局部化分析,从而克服了傅里叶分析在处理非平稳信号和复杂图像时所存在的局限性。
但是使用小波变换法去除伪迹要求信号和噪声的频带不能混叠。而脑电信号与眼电伪迹的频带相混叠,于是研究者开始将传统的小波变换法与其它去除伪迹的方法相结合,取长补短,期望获得更好的效果。
5主成分分析法(Principal Components Analysis, PCA):主成分分析是一种多通道的信号处理方法,其基本思想是利用正交原理将一组相关变量转化成为另一组相互独立的变量(即主成分),分解出来的主成分按大小排序,通常前面几个主成分就包括了信号80%的信息,因此PCA也被用来数据降维,去除冗余。PCA可以把多导联的脑电信号分解为相互正交的成分,去掉噪声和伪迹成分,再对其余成分进行重构,就可以得到较为纯净的脑电信号。
PCA的主要缺点在于:1信号和噪音波形相似时,难以从EEG中去除伪迹;2它只对正交的信号源比较有效,在去除伪迹时只能去除与脑电信号相正交的干扰成分;3它只适合于分解低阶相关信号,难以分解高阶的相关信号。
6独立分量分析(Independent Component Analysis, ICA): 独立分量分析理论的发展可追溯到上个世纪80年代初期,法国学者J.Hernah和C.Juttne等人首次提出了ICA分析的基本概念,最初被用来解决“鸡尾酒会”问题,随着近年来在ICA方面研究兴趣的增加,该理论已被广泛应用在图像处理、天线阵列、控制科学、生物医学信号处理等领域,其研究热潮方兴未艾。
ICA是实现盲源分离 (Blind Source Separation, BSS)最主要的方法之一,是一种基于高阶统计量的分析方法,处理对象为非高斯信号。理论上可以认为脑电信号、眼电伪迹、心电伪迹、肌电伪迹以及其他干扰源所产生的干扰信号都是由相互独立的源产生的,而且可以将观测信号视为是伪迹与脑电信号的线性组合,并且伪迹源的数目通常比头皮上测量的脑电通道数要少,因而独立成分分析方法适用于分离脑电信号中的各种伪迹。在实际应用中,ICA方法可以很好的分离频带相重叠的信号,弥补了小波分析去除伪迹的不足。但是ICA方法也有缺陷:在分离出信号和伪迹后,需要人工识别并筛选出伪迹成分再进行去除,费时费力,且具有很大的主观性。
发明内容
本发明针对传统脑电信号中眼电伪迹去除的过程中需要人工对眼电伪迹进行识别、费时费力、工作量大的缺陷,提出了一种眼电伪迹自动去除的方法。
其主要特征包括:首先利用独立分量分析对含有眼电伪迹的脑电信号进行分解,得到若干个独立分量,然后求取每个独立分量的峭度、序列renyi熵和样本熵作为区分脑电信号和眼电伪迹成分的特征,进而采用k均值聚类分析——一种无监督的分类算法进行特征分类,通过将含有眼电伪迹的独立分量置零重构出纯净脑电信号。
具体步骤如下:
1独立分量分解
本发明采用1999年芬兰学者Aapo Hyvarinen等人提出的基于负熵的固定点ICA算法,简称fastica算法,进行独立分量分解。设n通道观测信号x(t)=[x1(t),x2(t),…,xn(t)],经过ICA分解得到m个独立分量y(t)=[y1(t),y2(t),…,ym(t)]。
2特征提取
经过ICA分解得到m个独立分量y(t)=[y1(t),yi(t),…,ym(t)],分别求取它们的峭度、序列renyi熵和样本熵作为分类特征。
2.1 峭度:峭度是时间序列的四阶累计量,可以表征信号的波动情况。眼电伪迹的峭度值要显著大于脑电信号,因此,利用该参数作为识别眼电伪迹的特征。
2.2 renyi熵:renyi熵是一种反映信号随机性的参数,眼电伪迹的renyi熵值明显小于脑电成分,因而使用其作为区分眼电伪迹和脑电成分的特征。
2.3 样本熵:样本熵是一种测量时间序列复杂度的方法,而且眼电伪迹的样本熵值比脑电成分低,因此使用样本熵作为特征区分眼电伪迹和脑电成分。
3特征分类:本发明使用k均值聚类分析的方法进行分类,这是一种动态的、无监督的分类方法。区分眼电伪迹和脑电成分显然是一个两分类问题,因此先随机取两个样本作为聚类中心,计算每个样本到这两个中心的距离,把样本归入最近的聚类中心,然后用本类所有样本的均值代替原有的聚类中心的值,这样就得到了新的聚类中心,再重新计算各样本到新的聚类中心的距离,重新分类、修改中心点,直到新的聚类中心与上一次的中心点一致时,算法停止,然后再按照最小距离原则将独立分量分为两类,算法流程如图2所示。
4信号的重构
将含有眼电伪迹的独立分量置零,其它独立分量不变,对信号进行重构,得到纯净的脑电成分。
具体如下:
一种无需人工设定阈值的脑电信号中眼电伪迹自动去除的方法,其特征在于:首先,对采集的脑电信号进行独立分量分解,得到m个独立分量;其次,对于每一个独立分量,求取它们的峭度、序列renyi熵和样本熵作为特征向量;然后,使用k均值聚类分析——一种无监督的分类算法自动识别出含有眼电伪迹的独立分量;最后,将含有眼电伪迹的独立分量置零,其余独立分量不变,对信号进行重构,得到纯净的脑电信号;包括以下步骤:
1)独立分量分解
使用fastICA算法对信号进行分解,它的基本模型表达式为:
式中,x(t)=[x1(t),x2(t),…,xn(t)]T∈Rn×N是n通道的观测向量,为测量得到,1≤t≤N,N表示采样点总数,s(t)∈Rm×N为源信号,混合矩阵Rn×m,独立分量分析的目标是求解分离矩阵w,通过w从观测信号x(t)中恢复出未知的源信号s(t),源信号s(t)使用下式估计得到:
式中,y(t)是源信号s(t)的估计,y(t)=[y1(t),yi(t),…,ym(t)]T∈Rm×N是求得的m个独立分量,其中yi(t)是第i个独立分量,1≤i≤m,i为整数,t表示第t个采样点,分离矩阵w∈Rm×n;
2)特征提取
2.1 峭度:对于每一个独立分量yi(t),求取它的峭度;
式中,k为独立分量yi(t)的峭度值,1≤i≤m,i为整数;mβ为yi(t)的β阶中心矩,mβ=E{(yi(t)-m1)β},β={2,4},E(·)表示对括号里的表达式求均值,m1为yi(t)的均值;
2.2 renyi熵:对于每一个独立分量yi(t),其计算公式为
2.3 样本熵:对于每一个独立分量yi(t),其样本熵的计算过程如下:
第一步:设独立分量yi(t)有N个采样点,t表示第t个采样点,1≤t≤N,t为整数;
第二步:从独立分量yi(t)中提取b维矢量,b=2,即:
式中,yi(t)和yi(t+1)是独立分量yi(t)在采样点t和t+1的采样值,1≤t≤N-1,t为整数;
第三步:使用第二步的方法得到两个b维矢量为Yb(t)和Yb(t1),定义它们之间的距离为:
式中,1≤t,t1≤N-1,t≠t1,t、t1为整数,Yb(t)和Yb(t1)为公式(5)中Yb(t)在采样点t和t1的值;yi(t+c)和yi(t1+c)为独立分量yi(t)在采样点t+c和t1+c的采样值,c∈{0,1},求取差值的最大值,即为两个b维矢量的距离;
第四步:设定一个阈值r,r=(0.1~0.25)std(yi(t)),std(yi(t))为yi(t)的标准差;对于每一个t值,统计d[Yb(t),Yb(t1)]小于r的数目,记为T,1≤t,t1≤N-1,t≠t1,t、t1为整数,t1可以取N-2个值,因此,总的距离数为N-2,计算T与总的距离数N-2的比值,记为,即
第五步:计算的均值,即
第六步:再从同一个独立分量yi(t)中提取b1维矢量,b1=3,即:
式中,yi(t)、yi(t+1)和yi(t+2)是独立分量yi(t)在采样点t、t+1和t+2的采样值,1≤t≤N-2,t为整数;
式中,1≤t,t1≤N-2,t≠t1,t、t1为整数,和为公式(9)中在采样点t、t1的采样值;yi(t+e)和yi(t1+e)为独立分量yi(t)在采样点t+e和t1+e的采样值,e∈{0,1,2},求取差值的最大值,即为两个b1维矢量的距离;
第八步:在第四步中已计算得到阈值r,对于每一个t值,统计小于r的数目,记为T1,1≤t,t1≤N-2,t≠t1,t、t1为整数,t1可以取N-3个值,因此总的距离数为N-3,计算T1与总的距离数N-3的比值,记为,即
第十步:
样本熵的计算公式为
3)特征分类
经过上面的计算,对于每个独立分量yi(t) {1≤i≤m},可计算得三个特征量,将其依序构成如下特征向量,即:
式中,f1i、f2i、f3i分别表示yi(t)的峭度、renyi熵和样本熵;分别针对每一个特征进行一次k均值聚类分析,共进行三次分类,以f1i为例,给出计算步骤:
3.1给定样本种类,选定初始聚类中心;显然,这是一个两分类问题,初始聚类中心选为f11,f12两个分量,即center1(0)=f11,center2(0)=f12; centerq(0)表示第q类初始聚类中心,q=1或2,第p次迭代计算的聚类中心用centerq(p)表示;
3.2将样本按最小距离原则分为两个类,即在第p次迭代时,若||f1i-center1(p)||<||f1i-center2(p)||,则f1i∈c1(p),否则f1i∈c2(p);||·||表示对其中的两个向量求取它们之间的欧式距离;cq(p)表示在第p次迭代计算中第q类的聚类域,f1i∈cq(p)表示将样本f1i划分到第q类,q=1或2;
3.4若两次迭代计算的聚类中心相同,即centerq(p+1)=centerq(p),则计算完毕,否则令p=p+1,返回3.2;
4信号的重构
眼电伪迹的峭度值明显高于脑电信号,而renyi熵和样本熵的数值低于脑电信号,因此,综合三次分类的结果,将所有独立分量分为c1和c2两类,设c1类表示含有眼电伪迹的独立分量,c2类表示脑电独立分量;将c1类中所有成分置零,c2不变,得到去除眼电伪迹后的独立分量,利用式(15)进行重构,即可得到纯净的脑电成分;
式中,w-1为第1部分求得的分离矩阵w的伪逆矩阵。
本方法无需人工设定阈值就可以达到自动识别并去除眼电伪迹的目的,改善了以往方法中设定阈值时需要研究人员具备一定的先验知识、主观性强的不足附图说明
图1 脑电信号中眼电伪迹自动去除的算法流程图
图2 K均值聚类分析算法流程图
图3 国际10-20***标准电极放置图
图4 C3导联眼电伪迹去除效果图
图5 Cz导联眼电伪迹去除效果图
图6 C4导联眼电伪迹去除效果图
具体实施方式
本发明提出一种无需人工设定阈值的脑电信号中眼电伪迹的自动去除方法,完整的算法流程包括下面四(1.2.3.4)个部分。其中,第1部分是现有方法,本发明申请的特征包括三(2.3.4)个部分:
完整的算法流程如图1所示。
1独立分量分解
本发明采用fastICA算法,该算法有如下优点:①收敛速度比批处理和自适应处理都快;②负熵作为高斯性度量的效果优于累计量;③采用牛顿迭代法,收敛有保证;
它的基本模型表达式为:
式中,x(t)=[x1(t),x2(t),…,xn(t)]T∈Rn×N是n通道的观测向量,为测量得到,1≤t≤N,N表示采样点总数,s(t)∈Rm×N为源信号,混合矩阵A∈Rn×m,独立分量分析的目标是求解分离矩阵w,通过w从观测信号x(t)中恢复出未知的源信号s(t),源信号s(t)使用下式估计得到:
式中,y(t)是源信号s(t)的估计,y(t)=[y1(t),yi(t),…,ym(t)]T∈Rm×N是求得的m个独立分量,其中yi(t)是第i个独立分量,1≤i≤m,i为整数,t表示第t个采样点,分离矩阵w∈Rm×n。
2特征提取
2.1 峭度:对于每一个独立分量yi(t),求取它的峭度。
式中,k为独立分量yi(t)的峭度值,1≤i≤m,i为整数;mβ为yi(t)的β阶中心矩,mβ=E{(yi(t)-m1)β},β={2,4},E(·)表示对括号里的表达式求均值,m1为yi(t)的均值。
2.2 renyi熵:对于每一个独立分量yi(t),其计算公式为
2.3 样本熵:对于每一个独立分量yi(t),其样本熵的计算过程如下:
第一步:设独立分量yi(t)有N个采样点,t表示第t个采样点,1≤t≤N, t为整数。
第二步:从独立分量yi(t)中的提取b维矢量,b=2,即:
式中,yi(t)和yi(t+1)是独立分量yi(t)在采样点t和t+1的采样值,1≤t≤N-1,t为整数。
第三步:使用第二步的方法得到两个b维矢量为Yb(t)和Yb(t1),定义它们之间的距离为:
式中,1≤t,t1≤N-1,t≠t1,t、t1为整数,Yb(t)和Yb(t1)为公式(6)中Yb(t)在采样点t和t1的值。yi(t+c)和yi(t1+c)为独立分量yi(t)在采样点t+c和t1+c的采样值,c∈{0,1},求取差值的最大值,即为两个b维矢量的距离。
第四步:设定一个阈值r,r=(0.1~0.25)std(yi(t)),std(yi(t))为yi(t)的标准差。对于每一个t值,统计d[Yb(t),Yb(t1)]小于r的数目,记为T,1≤t,t1≤N-1,t≠t1,t、t1为整数,t1可以取N-2个值,因此,总的距离数为N-2,计算T与总的距离数N-2的比值,记为,即
第六步:再从同一个独立分量yi(t)中的提取b1维矢量,b1=3,即:
式中,yi(t)、yi(t+1)和yi(t+2)是独立分量yi(t)在采样点t、t+1和t+2的采样值,1≤t≤N-2,t为整数。
第七步:使用第六步的方法得到两个b1维矢量为和,定义它们之间的距离为:
式中,1≤t,t1≤N-2,t≠t1,t、t1为整数,和为公式(10)中在采样点t、t1的采样值。yi(t+e)和yi(t1+e)为独立分量yi(t)在采样点t+e和t1+e的采样值,e∈{0,1,2},求取差值的最大值,即为两个b1维矢量的距离。
第八步:在第四步中已计算得到阈值r,对于每一个t值,统计小于r的数目,记为T1,1≤t,t1≤N-2,t≠t1,t、t1为整数,t1可以取N-3个值,因此总的距离数为N-3,计算T1与总的距离数N-3的比值,记为,即
第十步:
样本熵的计算公式为
3特征分类
经过上面的计算,对于每个独立分量yi(t) {1≤i≤m},可计算得到三个特征量,将其依序构成如下特征向量,即:
式中,f1i、f2i、f3i分别表示yi(t)的峭度、renyi熵和样本熵。分别针对每一个特征进行k均值聚类分析,共进行三次分类,以f1i为例,给出计算步骤:
3.1给定样本种类,选定初始聚类中心;显然,这是一个两分类问题,初始聚类中心选为f11,f12两个分量。即center1(0)=f11,center2(0)=f12;centerq(0); centerq(0)表示第q类初始聚类中心,q=1或2,第p次迭代计算的聚类中心用centerq(p)表示。
3.2将样本按最小距离原则分为两个类,即在第p次迭代时,若||f1i-center1(p)||<||f1i-center2(p)||,则f1i∈c1(p),否则f1i∈c2(p)。||·||表示对其中的两个向量求取它们之间的欧式距离。cq(p)表示在第p次迭代计算中第q类的聚类域,f1i∈cq(p)表示将样本f1i划分到第q类,q=1或2。
3.4若两次迭代计算的聚类中心相同,即centerq(p+1)=centerq(p),则计算完毕,否则令p=p+1,返回3.2。
4信号的重构
眼电伪迹的峭度值要明显高于脑电信号,而renyi熵和样本熵的数值要低于脑电信号,因此,综合三次分类的结果,将所有独立分量分为c1和c2两类,设c1类表示含有眼电伪迹的独立分量,c2类表示脑电独立分量。将c1类中所有成分置零,c2不变,得到去除眼电伪迹后的独立分量,使用下式进行信号重构
为了给出详细具体的实施方式和步骤,将使用本方法对国际BCI(Brain Computer Interface, BCI)竞赛数据进行处理,给出实验结果并验证本方法去眼电伪迹的有效性。
数据介绍
脑电数据来源于第四届国际BCI竞赛数据集中的Data sets 2b数据包。采集的脑电数据经过0.5Hz~100Hz的带通滤波,50Hz的工频陷波,采样频率250Hz,幅值范围为±100μV。电极按照国际标准导联10-20***安放在C3、Cz、C4三个导联,如图3所示。另外,还同步采集了3导眼电信号。
实验步骤:
1.对数据进行ICA分解:输入数据如表1所示。
表1 原始数据表
得到的6个独立分量,数据如表2所示。
表2 独立分量数据表
2、以第一个独立分量y1(t)为例,分别计算它的峭度、renyi熵和样本熵。
根据公式(4),m1=E(y1(t))=-0.0176,m2=E{(y1(t)-m1)2}=1.000,
将已知参数带入公式(5),得到y1(t)的renyi熵的计算公式为
式中,y1(j)与y1(h)是独立分量y1(t)在采样时刻j和h的值。
根据第2部分中样本熵的计算步骤可以得到独立分量y1(t)的样本熵为
分别求取每个独立分量的峭度、renyi熵和样本熵作为特征向量,数据如表3所示。
表3 各独立分量的特征提取结果
独立分量 | 峭度 | Renyi熵 | 样本熵 |
1 | 6.00678780713594 | 1.41174521697880 | 1.70030048041082 |
2 | -0.588623024270334 | 1.44455029109777 | 0.52937952495178 |
3 | 0.31548101801454 | 1.57832435990357 | 1.54817199374019 |
4 | 0.44292110594564 | 1.57436808397370 | 1.75928182468544 |
5 | 0.39335875263220 | 1.56190706120367 | 0.96911146049469 |
6 | 0.24823401414372 | 1.58037659833919 | 1.50219022964535 |
3、眼电伪迹的识别:利用k均值聚类分析进行特征分类。结果如表4所示。
表4 各独立分量的分类结果
4、眼电伪迹的峭度值要明显高于脑电信号,而renyi熵和样本熵的数值要低于脑电信号。因此将1、2、5三个独立分量判别为眼电伪迹分量并将其置零,利用式(16)进行重构。眼电伪迹去除结果如图4-6所示。
本发明结合峭度、renyi熵和样本熵三个分别表征序列峰值、随机性和复杂度的参数作为特征向量,为说明这样提取特征向量的优势,下面分别使用单一特征和三种特征对眼电伪迹进行去除,得到4段去除眼电伪迹后的脑电信号,然后分别计算这4段信号与一段该受试者在实验过程中无眼电伪迹干扰的纯净脑电信号之间的均方误差,结果如表5所示。从表5可以得出,三种特征相结合去除眼电伪迹后得到的脑电信号与纯净脑电信号之间的均方误差最小,波形最为接近,这说明三种特征相结合去除眼电伪迹最为有效,并且对脑电信号破坏较小。
表5 不同特征去除伪迹后的均方误差对比表
峭度 | Renyi熵 | 样本熵 | 三种特征 | |
C3 | 60.0824114780632 | 15.9543462159107 | 16.4989572483297 | 13.6996491011524 |
Cz | 156.878955654143 | 84.3315031246680 | 53.1042782036112 | 50.5866649065173 |
C4 | 60.1610176410670 | 22.2264340506635 | 19.6937866400084 | 18.0135610303600 |
将本发明与希尔伯特—黄变换(Hilbert—Huang Transform,HHT)去除眼电伪迹的方法做一个对比,进一步说明本发明的有效性和优势。分别使用本方法和HHT法去除同一段脑电信号中的眼电伪迹,得到两段处理后的脑电信号,然后再分别计算这两段信号与一段该受试者在实验过程中无眼电伪迹干扰的纯净脑电信号之间的均方误差,结果如表6所示。从表6中可以看出使用本方法得到的均方误差要明显低于HHT法,证明使用本方法去除眼电伪迹后的脑电信号波形与纯净的脑电信号波形比较接近,即本方法可以更有效地去除眼电伪迹成分,并且对脑电成分破坏较小。
表6 两种方法均方误差对比
C3 | Cz | C4 | |
原始信号 | 62.8809 | 159.3931 | 61.8404 |
HHT法 | 19.6196 | 64.3851 | 18.7603 |
本方法 | 13.6996 | 50.5867 | 18.0136 |
Claims (1)
1.一种无需人工设定阈值的脑电信号中眼电伪迹自动去除的方法,其特征在于:首先,对采集的脑电信号进行独立分量分解,得到m个独立分量;其次,对于每一个独立分量,求取它们的峭度、序列renyi熵和样本熵作为特征向量;然后,使用k均值聚类分析——一种无监督的分类算法自动识别出含有眼电伪迹的独立分量;最后,将含有眼电伪迹的独立分量置零,其余独立分量不变,对信号进行重构,得到纯净的脑电信号;包括以下步骤:
1)独立分量分解
使用fastICA算法对信号进行分解,它的基本模型表达式为:
式中,x(t)=[x1(t),x2(t),…,xn(t)]T∈Rn×N是n通道的观测向量,为测量得到,1≤t≤N,N表示采样点总数,s(t)∈Rm×N为源信号,混合矩阵A∈Rn×m,独立分量分析的目标是求解分离矩阵w,通过w从观测信号x(t)中恢复出未知的源信号s(t),源信号s(t)使用下式估计得到:
式中,y(t)是源信号s(t)的估计,y(t)=[y1(t),yi(t),…,ym(t)]T∈Rm×N是求得的m个独立分量,其中yi(t)是第i个独立分量,1≤i≤m,i为整数,t表示第t个采样点,分离矩阵w∈Rm×n;
2)特征提取
2.1 峭度:对于每一个独立分量yi(t),求取它的峭度;
式中,k为独立分量yi(t)的峭度值,1≤i≤m,i为整数;mβ为yi(t)的β阶中心矩,mβ=E{(yi(t)-m1)β},β={2,4},E(·)表示对括号里的表达式求均值,m1为yi(t)的均值;
2.2 renyi熵:对于每一个独立分量yi(t),其计算公式为
2.3 样本熵:对于每一个独立分量yi(t),其样本熵的计算过程如下:
第一步:设独立分量yi(t)有N个采样点,t表示第t个采样点,1≤t≤N,t为整数;
第二步:从独立分量yi(t)中提取b维矢量,b=2,即:
式中,yi(t)和yi(t+1)是独立分量yi(t)在采样点t和t+1的采样值,1≤t≤N-1,t为整数;
第三步:使用第二步的方法得到两个b维矢量为Yb(t)和Yb(t1),定义它们之间的距离为:
式中,1≤t,t1≤N-1,t≠t1,t、t1为整数,Yb(t)和Yb(t1)为公式(5)中Yb(t)在采样点t和t1的值;yi(t+c)和yi(t1+c)为独立分量yi(t)在采样点t+c和t1+c的采样值,c∈{0,1},求取差值的最大值,即为两个b维矢量的距离;
第四步:设定一个阈值r,r=(0.1~0.25)std(yi(t)),std(yi(t))为yi(t)的标准差;对于每一个t值,统计d[Yb(t),Yb(t1)]小于r的数目,记为T,1≤t,t1≤N-1,t≠t1,t、t1为整数,t1可以取N-2个值,因此,总的距离数为N-2,计算T与总的距离数N-2的比值,记为,即
第六步:再从同一个独立分量yi(t)中提取b1维矢量,b1=3,即:
式中,yi(t)、yi(t+1)和yi(t+2)是独立分量yi(t)在采样点t、t+1和t+2的采样值,1≤t≤N-2,t为整数;
式中,1≤t,t1≤N-2,t≠t1,t、t1为整数,和为公式(9)中在采样点t、t1的采样值;yi(t+e)和yi(t1+e)为独立分量yi(t)在采样点t+e和t1+e的采样值,e∈{0,1,2},求取差值的最大值,即为两个b1维矢量的距离;
第八步:在第四步中已计算得到阈值r,对于每一个t值,统计小于r的数目,记为T1,1≤t,t1≤N-2,t≠t1,t、t1为整数,t1可以取N-3个值,因此总的距离数为N-3,计算T1与总的距离数N-3的比值,记为,即
第十步:
样本熵的计算公式为
3)特征分类
经过上面的计算,对于每个独立分量yi(t) {1≤i≤m},可计算得三个特征量,将其依序构成如下特征向量,即:
式中,f1i、f2i、f3i分别表示yi(t)的峭度、renyi熵和样本熵;分别针对每一个特征进行一次k均值聚类分析,共进行三次分类,以f1i为例,给出计算步骤:
3.1给定样本种类,选定初始聚类中心;显然,这是一个两分类问题,初始聚类中心选为f11,f12两个分量,即center1(0)=f11,center2(0)=f12;centerq(0)表示第q类初始聚类中心,q=1或2,第p次迭代计算的聚类中心用centerq(p)表示;
3.2将样本按最小距离原则分为两个类,即在第p次迭代时,若||f1i-center1(p)||<||f1i-center2(p)||,则f1i∈c1(p),否则f1i∈c2(p);||·||表示对其中的两个向量求取它们之间的欧式距离;cq(p)表示在第p次迭代计算中第q类的聚类域,f1i∈cq(p)表示将样本f1i划分到第q类,q=1或2;
3.3重新计算聚类中心;以本类所有样本的均值代替原有的聚类中心值;即,n1为第q类样本总数;
3.4若两次迭代计算的聚类中心相同,即centerq(p+1)=centerq(p),则计算完毕,否则令p=p+1,返回3.2;
4信号的重构
眼电伪迹的峭度值明显高于脑电信号,而renyi熵和样本熵的数值低于脑电信号,因此,综合三次分类的结果,将所有独立分量分为c1和c2两类,设c1类表示含有眼电伪迹的独立分量,c2类表示脑电独立分量;将c1类中所有成分置零,c2不变,得到去除眼电伪迹后的独立分量,利用式(15)进行重构,即可得到纯净的脑电成分;
式中,w-1为第1部分求得的分离矩阵w的伪逆矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210331193.XA CN102835955B (zh) | 2012-09-08 | 2012-09-08 | 一种无需设定阈值的脑电信号中眼电伪迹自动去除的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210331193.XA CN102835955B (zh) | 2012-09-08 | 2012-09-08 | 一种无需设定阈值的脑电信号中眼电伪迹自动去除的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102835955A true CN102835955A (zh) | 2012-12-26 |
CN102835955B CN102835955B (zh) | 2014-02-26 |
Family
ID=47363881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210331193.XA Expired - Fee Related CN102835955B (zh) | 2012-09-08 | 2012-09-08 | 一种无需设定阈值的脑电信号中眼电伪迹自动去除的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102835955B (zh) |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103690163A (zh) * | 2013-12-21 | 2014-04-02 | 哈尔滨工业大学 | 基于ica和hht融合的自动眼电干扰去除方法 |
CN103876734A (zh) * | 2014-03-24 | 2014-06-25 | 北京工业大学 | 一种基于决策树的脑电信号特征选择方法 |
CN103939325A (zh) * | 2014-05-05 | 2014-07-23 | 重庆大学 | 一种在低速运转下的消防水泵的故障诊断方法 |
CN104188649A (zh) * | 2014-09-15 | 2014-12-10 | 南京大学 | 多点生理电监测中保障信号线性合成实时性的一种方法 |
CN104688220A (zh) * | 2015-01-28 | 2015-06-10 | 西安交通大学 | 一种去除脑电信号中眼电伪迹的方法 |
CN104714171A (zh) * | 2015-04-06 | 2015-06-17 | 长沙学院 | 基于小波变换和ica特征提取的开关电路故障分类方法 |
CN105249962A (zh) * | 2015-11-03 | 2016-01-20 | 北京联合大学 | 头皮脑电信号回顾性癫痫发作点检测方法及*** |
CN105320969A (zh) * | 2015-11-20 | 2016-02-10 | 北京理工大学 | 基于多尺度Renyi熵的心率变异性特征分类方法 |
CN105342604A (zh) * | 2015-11-10 | 2016-02-24 | 中国航天员科研训练中心 | 基于脑电幅频特性的ica伪迹识别与去除方法及装置 |
CN106859641A (zh) * | 2017-02-20 | 2017-06-20 | 华南理工大学 | 一种基于自动ica去除eeg信号中核磁伪迹的方法 |
CN107080534A (zh) * | 2017-02-24 | 2017-08-22 | 浙江大学 | 基于eeg源成像的癫痫灶点粗定位方法 |
CN107260166A (zh) * | 2017-05-26 | 2017-10-20 | 昆明理工大学 | 一种实用化在线脑电伪迹剔除方法 |
CN107374620A (zh) * | 2017-08-21 | 2017-11-24 | 南京理工大学 | 一种基于独立成分分析算法的脑电信号预处理方法 |
CN108366730A (zh) * | 2015-10-05 | 2018-08-03 | 塔塔咨询服务公司 | 对eeg信号进行预处理用于认知负荷测量的方法和*** |
CN108523873A (zh) * | 2018-01-31 | 2018-09-14 | 北京理工大学 | 基于分数阶傅里叶变换和信息熵的心电信号特征提取方法 |
CN108836321A (zh) * | 2018-05-03 | 2018-11-20 | 江苏师范大学 | 一种基于自适应噪声消除***的脑电信号预处理方法 |
CN109086686A (zh) * | 2018-07-12 | 2018-12-25 | 西安电子科技大学 | 基于自适应动量因子的时变信道下的盲源分离方法 |
CN109222967A (zh) * | 2018-10-26 | 2019-01-18 | 哈尔滨工业大学 | 基于离群点三次样条插值的脑电信号中眼电伪迹去除方法 |
CN109512424A (zh) * | 2018-11-16 | 2019-03-26 | 福州大学 | 一种高密度或多通道肌电信号的肌肉激活起点检测方法 |
CN109645988A (zh) * | 2018-11-02 | 2019-04-19 | 杭州妞诺科技有限公司 | 便携式的脑电信号监测方法及*** |
CN109993132A (zh) * | 2019-04-04 | 2019-07-09 | 北京理工大学 | 一种基于脑电信号的图形识别生成方法及*** |
CN110292376A (zh) * | 2018-03-22 | 2019-10-01 | 深圳先进技术研究院 | 去除脑电信号中眼电伪迹的方法、装置、设备及存储介质 |
CN111956217A (zh) * | 2020-07-15 | 2020-11-20 | 山东师范大学 | 一种面向实时脑电信号的眨眼伪迹识别方法及*** |
CN112580436A (zh) * | 2020-11-25 | 2021-03-30 | 重庆邮电大学 | 一种基于黎曼流形坐标对齐的脑电信号域适应方法 |
CN113425297A (zh) * | 2021-07-19 | 2021-09-24 | 山东女子学院 | 基于脑电信号的儿童注意力专注度训练方法及*** |
CN114082169A (zh) * | 2021-11-22 | 2022-02-25 | 江苏科技大学 | 基于脑电信号的伤残手软体康复机器人运动想象识别方法 |
WO2023116263A1 (zh) * | 2021-12-24 | 2023-06-29 | 南京邮电大学 | 一种基于脑电图的眨眼咬牙判断方法及*** |
CN116612310A (zh) * | 2023-07-17 | 2023-08-18 | 长春医学高等专科学校(长春职工医科大学长春市医学情报所) | 基于多媒体舞蹈动作图像分解处理方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1883384A (zh) * | 2006-06-22 | 2006-12-27 | 复旦大学 | 一种脑电信号伪迹成分的自动识别和去除方法 |
WO2010045724A1 (en) * | 2008-10-20 | 2010-04-29 | Neurochip Corporation | Method and rhythm extractor for detecting and isolating rhythmic signal features from an input signal using the wavelet packet transform |
CN101869477A (zh) * | 2010-05-14 | 2010-10-27 | 北京工业大学 | 一种自适应脑电信号中眼电伪迹的自动去除方法 |
-
2012
- 2012-09-08 CN CN201210331193.XA patent/CN102835955B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1883384A (zh) * | 2006-06-22 | 2006-12-27 | 复旦大学 | 一种脑电信号伪迹成分的自动识别和去除方法 |
WO2010045724A1 (en) * | 2008-10-20 | 2010-04-29 | Neurochip Corporation | Method and rhythm extractor for detecting and isolating rhythmic signal features from an input signal using the wavelet packet transform |
CN101869477A (zh) * | 2010-05-14 | 2010-10-27 | 北京工业大学 | 一种自适应脑电信号中眼电伪迹的自动去除方法 |
Cited By (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103690163A (zh) * | 2013-12-21 | 2014-04-02 | 哈尔滨工业大学 | 基于ica和hht融合的自动眼电干扰去除方法 |
CN103876734A (zh) * | 2014-03-24 | 2014-06-25 | 北京工业大学 | 一种基于决策树的脑电信号特征选择方法 |
CN103876734B (zh) * | 2014-03-24 | 2015-09-02 | 北京工业大学 | 一种基于决策树的脑电信号特征选择方法 |
CN103939325A (zh) * | 2014-05-05 | 2014-07-23 | 重庆大学 | 一种在低速运转下的消防水泵的故障诊断方法 |
CN104188649A (zh) * | 2014-09-15 | 2014-12-10 | 南京大学 | 多点生理电监测中保障信号线性合成实时性的一种方法 |
CN104688220A (zh) * | 2015-01-28 | 2015-06-10 | 西安交通大学 | 一种去除脑电信号中眼电伪迹的方法 |
CN104714171A (zh) * | 2015-04-06 | 2015-06-17 | 长沙学院 | 基于小波变换和ica特征提取的开关电路故障分类方法 |
CN108366730A (zh) * | 2015-10-05 | 2018-08-03 | 塔塔咨询服务公司 | 对eeg信号进行预处理用于认知负荷测量的方法和*** |
CN108366730B (zh) * | 2015-10-05 | 2022-07-05 | 塔塔咨询服务公司 | 对eeg信号进行预处理用于认知负荷测量的方法和*** |
CN105249962A (zh) * | 2015-11-03 | 2016-01-20 | 北京联合大学 | 头皮脑电信号回顾性癫痫发作点检测方法及*** |
CN105249962B (zh) * | 2015-11-03 | 2019-04-30 | 北京联合大学 | 头皮脑电信号回顾性癫痫发作点检测方法及*** |
CN105342604B (zh) * | 2015-11-10 | 2018-08-07 | 中国航天员科研训练中心 | 基于脑电幅频特性的ica伪迹识别与去除方法及装置 |
CN105342604A (zh) * | 2015-11-10 | 2016-02-24 | 中国航天员科研训练中心 | 基于脑电幅频特性的ica伪迹识别与去除方法及装置 |
CN105320969A (zh) * | 2015-11-20 | 2016-02-10 | 北京理工大学 | 基于多尺度Renyi熵的心率变异性特征分类方法 |
CN106859641B (zh) * | 2017-02-20 | 2019-08-20 | 华南理工大学 | 一种基于自动ica去除eeg信号中核磁伪迹的方法 |
CN106859641A (zh) * | 2017-02-20 | 2017-06-20 | 华南理工大学 | 一种基于自动ica去除eeg信号中核磁伪迹的方法 |
CN107080534A (zh) * | 2017-02-24 | 2017-08-22 | 浙江大学 | 基于eeg源成像的癫痫灶点粗定位方法 |
CN107260166A (zh) * | 2017-05-26 | 2017-10-20 | 昆明理工大学 | 一种实用化在线脑电伪迹剔除方法 |
CN107374620A (zh) * | 2017-08-21 | 2017-11-24 | 南京理工大学 | 一种基于独立成分分析算法的脑电信号预处理方法 |
CN108523873A (zh) * | 2018-01-31 | 2018-09-14 | 北京理工大学 | 基于分数阶傅里叶变换和信息熵的心电信号特征提取方法 |
CN108523873B (zh) * | 2018-01-31 | 2021-11-16 | 北京理工大学 | 基于分数阶傅里叶变换和信息熵的心电信号特征提取方法 |
CN110292376A (zh) * | 2018-03-22 | 2019-10-01 | 深圳先进技术研究院 | 去除脑电信号中眼电伪迹的方法、装置、设备及存储介质 |
CN108836321A (zh) * | 2018-05-03 | 2018-11-20 | 江苏师范大学 | 一种基于自适应噪声消除***的脑电信号预处理方法 |
CN109086686A (zh) * | 2018-07-12 | 2018-12-25 | 西安电子科技大学 | 基于自适应动量因子的时变信道下的盲源分离方法 |
CN109086686B (zh) * | 2018-07-12 | 2022-09-30 | 西安电子科技大学 | 基于自适应动量因子的时变信道下的盲源分离方法 |
CN109222967A (zh) * | 2018-10-26 | 2019-01-18 | 哈尔滨工业大学 | 基于离群点三次样条插值的脑电信号中眼电伪迹去除方法 |
CN109645988A (zh) * | 2018-11-02 | 2019-04-19 | 杭州妞诺科技有限公司 | 便携式的脑电信号监测方法及*** |
CN109512424A (zh) * | 2018-11-16 | 2019-03-26 | 福州大学 | 一种高密度或多通道肌电信号的肌肉激活起点检测方法 |
CN109512424B (zh) * | 2018-11-16 | 2021-07-13 | 福州大学 | 一种高密度或多通道肌电信号的肌肉激活起点检测方法 |
CN109993132A (zh) * | 2019-04-04 | 2019-07-09 | 北京理工大学 | 一种基于脑电信号的图形识别生成方法及*** |
CN111956217B (zh) * | 2020-07-15 | 2022-06-24 | 山东师范大学 | 一种面向实时脑电信号的眨眼伪迹识别方法及*** |
CN111956217A (zh) * | 2020-07-15 | 2020-11-20 | 山东师范大学 | 一种面向实时脑电信号的眨眼伪迹识别方法及*** |
CN112580436B (zh) * | 2020-11-25 | 2022-05-03 | 重庆邮电大学 | 一种基于黎曼流形坐标对齐的脑电信号域适应方法 |
CN112580436A (zh) * | 2020-11-25 | 2021-03-30 | 重庆邮电大学 | 一种基于黎曼流形坐标对齐的脑电信号域适应方法 |
CN113425297A (zh) * | 2021-07-19 | 2021-09-24 | 山东女子学院 | 基于脑电信号的儿童注意力专注度训练方法及*** |
CN114082169A (zh) * | 2021-11-22 | 2022-02-25 | 江苏科技大学 | 基于脑电信号的伤残手软体康复机器人运动想象识别方法 |
WO2023116263A1 (zh) * | 2021-12-24 | 2023-06-29 | 南京邮电大学 | 一种基于脑电图的眨眼咬牙判断方法及*** |
CN116612310A (zh) * | 2023-07-17 | 2023-08-18 | 长春医学高等专科学校(长春职工医科大学长春市医学情报所) | 基于多媒体舞蹈动作图像分解处理方法 |
CN116612310B (zh) * | 2023-07-17 | 2023-09-26 | 长春医学高等专科学校(长春职工医科大学长春市医学情报所) | 基于多媒体舞蹈动作图像分解处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102835955B (zh) | 2014-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102835955B (zh) | 一种无需设定阈值的脑电信号中眼电伪迹自动去除的方法 | |
CN102697493B (zh) | 一种快速的脑电信号中眼电伪迹自动识别和去除的方法 | |
CN104688220B (zh) | 一种去除脑电信号中眼电伪迹的方法 | |
CN107157477A (zh) | 脑电信号特征识别***及方法 | |
Jadhav et al. | Automated detection and correction of eye blink and muscular artefacts in EEG signal for analysis of Autism Spectrum Disorder | |
Ahirwal et al. | Power spectrum analysis of EEG signals for estimating visual attention | |
CN105342605B (zh) | 一种去除脑电信号中肌电伪迹的方法 | |
CN111067514B (zh) | 一种基于多尺度多变量传递熵的多通道脑电耦合分析方法 | |
CN107260166A (zh) | 一种实用化在线脑电伪迹剔除方法 | |
CN106236080B (zh) | 基于多通道的脑电信号中肌电噪声的消除方法 | |
CN103110418A (zh) | 一种脑电信号特征提取方法 | |
CN1883384A (zh) | 一种脑电信号伪迹成分的自动识别和去除方法 | |
Li et al. | An EEG-based method for detecting drowsy driving state | |
CN113723557B (zh) | 一种基于多频带时空卷积网络的抑郁脑电分类*** | |
CN109480832A (zh) | 一种单通道的脑电信号中肌电伪迹的消除方法 | |
CN106264521A (zh) | 多通道脑电信号中下颌干扰的自动去除方法 | |
Haider et al. | Performance enhancement in P300 ERP single trial by machine learning adaptive denoising mechanism | |
CN108836301A (zh) | 一种基于奇异谱分析和稀疏表示的单次诱发电位提取方法 | |
Sheoran et al. | Methods of denoising of electroencephalogram signal: A review | |
CN117462145A (zh) | 基于FastICA的眼电伪迹自识别清除方法、装置、设备及介质 | |
Zhu et al. | An ICA-based method for automatic eye blink artifact correction in multi-channel EEG | |
CN111990991A (zh) | 一种基于复杂网络的脑电信号分析方法及应用 | |
Jaffino et al. | Expectation-maximization extreme machine learning classifier for epileptic seizure detection | |
CN111671421A (zh) | 一种基于脑电图的儿童需求感知方法 | |
Assi et al. | Kmeans-ICA based automatic method for ocular artifacts removal in a motorimagery classification |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140226 Termination date: 20190908 |
|
CF01 | Termination of patent right due to non-payment of annual fee |