CN114224360B - 一种基于改进emd-ica的eeg信号处理方法、设备及存储介质 - Google Patents
一种基于改进emd-ica的eeg信号处理方法、设备及存储介质 Download PDFInfo
- Publication number
- CN114224360B CN114224360B CN202111614480.7A CN202111614480A CN114224360B CN 114224360 B CN114224360 B CN 114224360B CN 202111614480 A CN202111614480 A CN 202111614480A CN 114224360 B CN114224360 B CN 114224360B
- Authority
- CN
- China
- Prior art keywords
- max
- ica
- signal
- emd
- electroencephalogram
- 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
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 29
- 238000010586 diagram Methods 0.000 claims abstract description 24
- 238000013527 convolutional neural network Methods 0.000 claims abstract description 22
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 11
- 230000009466 transformation Effects 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 238000000926 separation method Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 24
- 230000006870 function Effects 0.000 claims description 21
- 238000011176 pooling Methods 0.000 claims description 10
- 210000004556 brain Anatomy 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000012880 independent component analysis Methods 0.000 description 34
- 238000004458 analytical method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000007177 brain activity Effects 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 238000002565 electrocardiography Methods 0.000 description 1
- 238000002570 electrooculography Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000337 motor cortex Anatomy 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000012360 testing method Methods 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- 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/7253—Details of waveform analysis characterised by using transforms
- A61B5/726—Details of waveform analysis characterised by using transforms using Wavelet transforms
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Psychiatry (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychology (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于改进EMD‑ICA的EEG信号处理方法、设备及存储介质,涉及信号处理技术领域,该方法利用改进经验模式分解法EMD对脑电信号EEG进行分解,得到本征模态分量IMFS和余量;对得到的所有符合要求的本征模态分量IMFS进行叠加重构,对重构后的脑电信号通过ICA法进行信号分离,消除噪声,得到去噪后的脑电信号;利用连续小波变换将去噪后的脑电信号变成二维时频图,并将生成的二维时频图输入卷积神经网络模型;采用卷积神经网络模型对步骤S3中的二维时频图进行特征提取,并进行分类。该方法基于改进的EMD‑ICA对脑电信号进行去噪,采用去噪后的信号再进行特征提取和分类,使得运动想象脑电信号分类准确,准确率明显高于现有的其他的分类方法。
Description
技术领域
本发明涉及EEG信号处理方法技术领域,尤其涉及一种基于改进EMD-ICA的EEG信号处理方法、设备及存储介质。
背景技术
EEG信号易受噪声影响,例如眼电(EOG),肌电(EMG),心电(ECG),及电源线干扰。为了滤除接收信号中的噪声并恢复源信号,多年来开发的技术主要包括时域分析,频域分析和时频分析方法,例如Wiener(WD)分布和小波变换。
ICA是一种无监督的统计学习方法,可以将复杂的混合信号分解为独立的分量。在BCI-P300系统中,已经使用独立成分分析(ICA)将P300与背景噪声区分开来,以实现对运动图像任务的伪影识别。但是,由于ICA算法无法获得时域噪声信号的特征,如果将其视为噪声,则会去除某些区域的大脑活动成分,也就是说,一些有价值的大脑活动信息可能会丢失。
经验模式分解法(Empirical Mode Decomposition,简称EMD)对于分析非线性、非平稳信号序列具有很大的优势,同时具有高信噪比的特点。该方法的核心是经验模式分解,即将复杂的信号分解为按频率由高到低顺序排列的有限个本征模函数(Intrinsic ModeFunction,简称IMF),分解出来的IMF分量包含了原信号的不同时间尺度的局部特征信号。EMD的关键是使复杂信号分解为有限个本征模函数(IMF),所分解出来的各IMF分量包含了原信号的不同时间尺度的局部特征信号。
小波变换是1980年代后期发展起来的应用数学的一个分支。它是一种时频分解技术,并已广泛应用于信号处理,图像处理,语音识别等领域。由于通过WT后信号和噪声的统计特性不同,因此在多尺度分析中它们表现出不同的传播行为。该特性可用于对噪声信号进行降噪。但是,WT不能有效地保留噪声的时频结构,也无法单独恢复隐藏在噪声成分中的神经活动。
卷积神经网络通常由多个卷积层和池化层组成。典型的CNN网络有以下几个层次:输入层、卷积层、池化层、完全连接层和输出层。深度学习算法提取特征的能力与传统算法相比有很大的提升,而且通常网络越复杂提取的特征越充分,得到的分类器效果越好。作为一种多层神经网络,卷积神经网络凭借其强大的特征提取能力成功地应用在计算机视觉和图像处理等领域。
总之,由于EEG信号本身的特性,仅使用上述任一传统的去噪方法及特征提取分类难以实现理想的分类效果。因此,本文结合EMD、ICA算法、连续小波变换及CNN算法,提出一种基于改进的EMD-ICA的EEG信号处理方法。
发明内容
本发明的目的在于提供一种基于改进EMD-ICA的EEG信号处理方法,从而解决现有技术中存在的前述问题。
为了实现上述目的,本发明采用的技术方案如下:
一种基于改进EMD-ICA的EEG信号处理方法,包括以下步骤:
S1,利用改进经验模式分解法EMD对脑电信号EEG进行分解,得到本征模态分量IMFS和余量;
S2,对得到的所有符合要求的本征模态分量IMFS进行叠加重构,对重构后的脑电信号通过ICA法进行信号分离,消除噪声,得到去噪后的脑电信号;
S3,利用连续小波变换将去噪后的脑电信号变成二维时频图,并将生成的二维时频图输入卷积神经网络模型;
S4,采用卷积神经网络模型对步骤S3中的二维时频图进行特征提取,并进行分类。
优选的,步骤S1中具体包括:
S11,把原始脑电信号X(t)的起始端点作为一个极大值起点记为M(0),横坐标为t(0),选取一个相邻的极大值点记为M(1),横坐标为tm(1),两点连线,计算其斜率S1;
S12,计算起始端点处延展得到的极小值N(0);
N(0)=N(1)-S1·[tn(1)-tn(0)] (2)(公式变化了)
N(1)指的是距离起始端点处最近的极小值,tn(1)指的是N(1)的横坐标,tn(0)=t(0);
S13,相应的,把原始脑电信号的终端端点作为极小值点记为N(Qmin),横坐标为tn(Qmin),选取一个相邻的极小值点记为N(Qmin-1),横坐标为tn(Qmin-1),采用公式(3)计算两点连线的斜率S2:
S14,计算终端端点延展得到的极大值M(Qmax):
M(Qmax)=M(Qmax-1)-S2·[tm(Qmax-1)-t(N)] (4)
选取一个相邻的极大值点记为M(Qmax-1),横坐标为tm(Qmax-1),Tm(Qmax)=t(N);
S15,重复步骤S14,由此找出原始脑电信号数据序列内的所有的极大值点M1(Qmax)和极小值点N1(0),得到极大值序列;
找出极大值序列中的最大值MAX和最小值MIN,若M(Qmax)<MIN,则定义该延展极大值为:M1(Qmax)=(终端端点值+MAX)/2;若N(0)>MAX,则定义该延展极大值为:N1(0)=(起始端点值+MIN)/2;
S16,采用三次样条插值函数拟合形成原数据的上包络线Umax(t)和下包络线Umin(t);
根据上、下包络线计算原始数据X(t)的局部平均值,记作公式(5):
m1(t)=(Umax(t)+Umin(t))/2 (5)
S17,计算原始脑电信号与局部平均值的差值h1(t),记为公式(6):
h1(t)=X(t)-m1(t) (6)
S18,判断h1(t)是否满足IMF的两个条件,若满足,则h1(t)为第一个IMF分量,从原始信号X(t)中减去h1(t)得到剩余信号r1(t);否则就把h1(t)当成原始脑电信号转到步骤S16,重复步骤S16-S17,继续寻找原数据序列h1(t)信号内的所有的极大值点M1(Qmax)和极小值点N1(0),直至满足IMF的两个条件;
r1(t)=X(t)-h1(t) (7)
S19,将剩余信号r1(t)看作新的脑电信号重复步骤S11-S18,可筛选得到n个IMF分量:
rn(t)=rn-1(t)-hn(t) (8)
S110,对n个IMF分量h1(t),h2(t),…,hn(t)做阈值处理,将大于阈值V的分量置为零,如下式所示:
其中,其中,F为原始脑电信号数据X(t)的长度;
Median()函数,返回一组已知数字的中值,中值是一组数的中间数。
优选的,满足IMF的两个条件具体为:1)在整个数据段内,极值点和零点数量相等或相差至多1个;2)在任意时刻,由局部极大值点形成的上包络和由局部极小值点形成的下包络的平均值为0。
优选的,步骤S110后还包括:分解后的IMF分量和余量叠加重构为去噪后的重构脑电信号
优选的,步骤S2具体包括:
S21采用下式,对重构脑电信号去均值,得到处理后的脑电观测信号X’(t);
其中:E[·]表示数学期望;
S22,采用下式,将处理后的脑电观测信号X’(t),分解为各分量间互不相关的脑电信号Z(t)=(z1(t),…,zm(t))T:
其中:代表投影因子;
Ds是以X’(t)的协方差矩阵CX=E[X’(t)*X’(t)T]特征值为对角元素的对角矩阵;
Us是以CX的单位范数特性向量为列的矩阵;
I是单位矩阵;
σ表示X’(t)的噪声方差,
S23,定义Y(t)=WZ(t),W为解混矩阵;
初始化各参数:令循环变量i=1,初始解混矩阵为:
E[·]表示数学期望;
S24,计算Y(t)=WZ(t);
ΔW=λ[I-Ktanh(Y(t))(Y(t))T-Y(t)(Y(t))T]W
其中,λ为学习步长,本文取λ=0.001,I为单位矩阵;K为对角矩阵,其对角元素kii=sgn{E[sech2yi(t)]-E[yi(t)tanhyi(t)]}
Sgn()符号函数,返回参数的正负,sech()为双曲正割函数,tanh()为双曲正切函数。
S25,W=W+ΔW
返回步骤S24继续计算,直到W收敛得到独立分量Y(t)=WZ(t);处理后的解混矩阵W的无穷范数是否小于10-6,如果是,则代表收敛。
优选的,步骤S3中绘制二维时频图的具体过程包括:
S31,从C3和C4通道的脑电信号中提取时间、频率和位置信息;
S32,利用连续小波转换将步骤S31中获得的时间、频率和位置信息结合绘制出二维时频图。
优选的,步骤S4中采用的卷积神经网络模型采用2层卷积神经网络模型,具体包括2个卷积层,2个池化层,2个全连接层及1个softmax输出层。
优选的,步骤S4具体包括:
S41,将绘制的二维时频图首先传入第一卷积层,提取局部特征
S42,将提取局部特征后的特征图传入第一池化层,经过特征选择和信息过滤后再依次传入第二卷积层和第二池化层;
S43,进入全连接层,最后进入输出层,采用softmax分类函数即可得到左右手的分类。
本发明的另一个目的在于提供一种计算机可读存储介质,所述计算机可读存储介质包括处理器,所述处理器用于实现基于改进EMD-ICA的EEG信号处理方法。
本发明的最后一个目的在于提供一种基于改进EMD-ICA的EEG信号处理设备,用于执行基于改进EMD-ICA的EEG信号处理方法。
本发明的有益效果是:
本发明公开了一种基于改进EMD-ICA的EEG信号处理方法、设备及存储介质,本方法基于改进的EMD-ICA对脑电信号进行去噪,采用去噪后的信号再进行特征提取和分类,使得脑电运动想象信号进行准确分类,准确率明显高于现有的其他的分类方法。
附图说明
图1是实施例1中提供的基于改进EMD-ICA的EEG信号处理方法流程图;
图2是实施例2中提供的将脑电数据进行EMD-ICA分解的实验显示结果;
图3是利用连续小波变换将去噪后的C3脑电信号变成二维时频图;
图4是利用连续小波变换将去噪后的C4脑电信号变成二维时频图;
图5是图3和图4中的图像组合而成的最终输入图像;
图6是实施例2中采用的卷积神经网络结构模型示意图;
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
实施例1
本实施例提供一种基于改进EMD-ICA的EEG信号处理方法,包括以下步骤:
S1,利用改进经验模式分解法EMD对脑电信号EEG进行分解,得到本征模态分量IMFS和余量;
S2,对得到的所有符合要求的本征模态分量IMFS进行叠加重构,对重构的脑电信号通过ICA法进行信号分离,消除噪声,得到去噪后的脑电信号;
S3,利用连续小波变换将去噪后的脑电信号变成二维时频图,并将生成的二维时频图输入卷积神经网络模型;
S4,采用卷积神经网络模型对步骤S3中的二维时频图进行特征提取,并进行分类。
本实施例中的步骤S1中具体包括:
S11,把原始脑电信号X(t)的起始端点作为一个极大值起点记为M(0),横坐标为t(0),选取一个相邻的极大值点记为M(1),横坐标为tm(1),两点连线,计算其斜率S1;
S12,计算起始端点处延展得到的极小值N(0);
N(0)=N(1)-S1·[tn(1)-tn(0)] (2)(公式变化了)
N(1)指的是距离起始端点处最近的极小值,tn(1)指的是N(1)的横坐标,tn(0)=t(0);
S13,相应的,把原始脑电信号的终端端点作为极小值点记为N(Qmin),横坐标为tn(Qmin),选取一个相邻的极小值点记为N(Qmin-1),横坐标为tn(Qmin-1),采用公式(3)计算两点连线的斜率S2:
S14,计算终端端点延展得到的极大值M(Qmax):
M(Qmax)=M(Qmax-1)-S2·[tm(Qmax-1)-t(N)] (4)
选取一个相邻的极大值点记为M(Qmax-1),横坐标为tm(Qmax-1),Tm(Qmax)=t(N);
S15,重复步骤S14,由此找出原始脑电信号数据序列内的所有的极大值点M1(Qmax)和极小值点N1(0),得到极大值序列;
找出极大值序列中的最大值MAX和最小值MIN,若M(Qmax)<MIN,则定义该延展极大值为:M1(Qmax)=(终端端点值+MAX)/2;若N(0)>MAX,则定义该延展极大值为:N1(0)=(起始端点值+MIN)/2;
S16,采用三次样条插值函数拟合形成原数据的上包络线Umax(t)和下包络线Umin(t);
根据上、下包络线计算原始数据X(t)的局部平均值,记作公式(5):
m1(t)=(Umax(t)+Umin(t))/2 (5)
S17,计算原始脑电信号与局部平均值的差值h1(t),记为公式(6):
h1(t)=X(t)-m1(t) (6)
S18,判断h1(t)是否满足IMF的两个条件,若满足,则h1(t)为第一个IMF分量,从原始信号X(t)中减去h1(t)得到剩余信号r1(t);否则就把h1(t)当成原始脑电信号转到步骤S16,重复步骤S16-S17,继续寻找原数据序列h1(t)信号内的所有的极大值点M1(Qmax)和极小值点N1(0),直至满足IMF的两个条件;
r1(t)=X(t)-h1(t) (7)
S19,将剩余信号r1(t)看作新的脑电信号重复步骤S11-S18,可筛选得到n个IMF分量:
rn(t)=rn-1(t)-hn(t) (8)
S110,对n个IMF分量h1(t),h2(t),…,hn(t)做阈值处理,将大于阈值V的分量置为零,如下式所示:
其中,其中,F为原始脑电信号数据X(t)的长度;
Median()函数,返回一组已知数字的中值,中值是一组数的中间数。
本实施例中满足IMF的两个条件具体为:1)在整个数据段内,极值点和零点数量相等或相差至多1个;2)在任意时刻,由局部极大值点形成的上包络和由局部极小值点形成的下包络的平均值为0。
步骤S110后还包括:分解后的IMF分量和余量rn叠加重构为去噪后的重构脑电信号
本实施例中的步骤S2具体包括:
S21采用下式,对重构脑电信号去均值,得到处理后的脑电观测信号X’(t);
其中:E[·]表示数学期望;
S22,采用下式,将处理后的脑电观测信号X’(t),分解为各分量间互不相关的脑电信号Z(t)=(z1(t),…,zm(t))T:
其中:代表投影因子;
Ds是以X’(t)的协方差矩阵CX=E[X’(t)*X’(t)T]特征值为对角元素的对角矩阵;
Us是以CX的单位范数特性向量为列的矩阵;
I是单位矩阵;
σ表示X’(t)的噪声方差,
S23,定义Y(t)=WZ(t),W为解混矩阵;
初始化各参数:令循环变量i=1,初始解混矩阵为:
E[·]表示数学期望;
S24,计算Y(t)=WZ(t);
ΔW=λ[I-Ktanh(Y(t))(Y(t))T-Y(t)(Y(t))T]W
其中,λ为学习步长,本文取λ=0.001,I为单位矩阵;K为对角矩阵,其对角元素kii=sgn{E[sech2yi(t)]-E[yi(t)tanhyi(t)]}
Sgn()符号函数,返回参数的正负,sech()为双曲正割函数,tanh()为双曲正切函数。
S25,W=W+ΔW
返回步骤S24继续计算,直到W收敛得到独立分量Y(t)=WZ(t);处理后的解混矩阵W的无穷范数是否小于10-6,如果是,则代表收敛。
本实施例中的步骤S3中绘制二维时频图的具体过程包括:
S31,从C3和C4通道的脑电信号中提取时间、频率和位置信息;
S32,利用连续小波转换将步骤S31中获得的时间、频率和位置信息结合绘制出二维时频图。
本实施例中步骤S4中采用的卷积神经网络模型采用2层卷积神经网络模型,具体包括2个卷积层,2个池化层,2个全连接层及1个softmax输出层。
步骤S4具体包括:
S41,将绘制的二维时频图首先传入第一卷积层,提取局部特征;
S42,将提取局部特征后的特征图传入第一池化层,经过特征选择和信息过滤后再依次传入第二卷积层和第二池化层;
S43,进入全连接层,最后进入输出层,采用softmax分类函数即可得到左右手的分类。
实施例2
本实施例采用实施例1中所记载的基于改进EMD-ICA的EEG信号处理方法对一组EEG信号进行处理,包括以下步骤:
(1)EMD-ICA去噪
采用第四届国际脑机接口竞赛——格茨数据集2a,选取第一位受试者的脑电数据进行EMD-ICA分解,实验结果如图2所示,原始脑电信号不平滑,混杂着各种各样的噪声,对原始信号进行EMD-ICA分解后,得到的结果是比较理想的,该方法获得的去噪信号更符合原始数据的自然走向,由此可见EMD-ICA算法不仅精确地消除伪迹成分,更很好地保留原始EEG的局部特性。
(2)利用连续小波变换将去噪后的脑电信号变成二维时频图
左手和右手运动运动想象任务分别导致运动皮层左右两侧出现ERD和ERS现象(事件相关去同步化(event-related desynchronization,ERD)和事件相关同步化(eventrelated synchronization,ERS)),影响C4和C3电极的脑电图信号。图3和图4是在执行左手运动想象任务时C3和C4通道的小波时频图像。从图3、4中可以看出,在实验开始后3秒,C4通道的能量显著下降,一段时间后恢复,即ERD现象发生。然而,C3通道的能量保持在高水平而不是降低,这被称为ERS现象。将C3、C4的图像组合为最终输入图像,如图5所示。
(3)利用卷积神经网络对图像进行特征提取并分类
采用如图6所示的卷积神经网络模型,对图5进行特征提取并分类。
(4)对比试验
为了对算法进行定量评价,对9名受试者的数据分别进行ICA、ICA-WT和本发明中提供的EMD-ICA三种不同方法去噪处理,对比分类结果如表1所示,由表1中显示的内容可以知晓,本文提出的EMD-ICA方法的平均精度结果要高于ICA和ICA-WT。可以分析出,使用本文提出的EMD-ICA算法提高了信号的信噪比,很好的保留了信号的有效成分,得到了更理想的消噪效果,可以有效地提高分类精度。
表1 ICA、ICA-WT和EMD-ICA去噪算法的分类结果
kappa值是去除随机分类准确性影响的分类性能度量,其中,p0为分类精度,pe为随机分类的结果,两类分类时为0.5。
将本文方法的Kappa结果与当前的研究现状进行了比较,如表2所示。
表2 CNN、SAE和SVM的Kappa值
从表2可以看出来自每个受试者的运动想象脑电信号的质量不同,导致预测的准确性有很大差异。
本文提出的EMD-ICA算法对相同受试者的脑电图记录在不同时段的变化具有相当的鲁棒性,对比其他方法对受试者之间的差异更具有鲁棒性。结果表明,本文提出的EMD-ICA去噪方法具有较高的准确率和可靠性。
通过采用本发明公开的上述技术方案,得到了如下有益的效果:
本发明公开了一种基于改进EMD-ICA的EEG信号处理方法、设备及存储介质,本方法基于改进的EMD-ICA对脑电信号进行去噪,采用去噪后的信号再进行特征提取和分类,使得运动想象脑电信号进行准确分类,准确率明显高于现有的其他的分类方法。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。
Claims (8)
1.一种基于改进EMD-ICA的EEG信号处理方法,其特征在于,包括以下步骤:
S1,利用改进经验模式分解法EMD对脑电信号EEG进行分解,得到本征模态分量IMFS和余量;
S2,对得到的所有符合要求的本征模态分量IMFS进行叠加重构,对重构后的脑电信号通过改进ICA法进行信号分离,消除噪声,得到去噪后的脑电信号;
S3,利用连续小波变换将去噪后的脑电信号变成二维时频图,并将生成的二维时频图输入卷积神经网络模型;
S4,采用卷积神经网络模型对步骤S3中的二维时频图进行特征提取,并进行分类;
步骤S1中具体包括:
S11,把原始脑电信号X(t)的起始端点作为一个极大值起点记为M(0),横坐标为t(0),选取一个相邻的极大值点记为M(1),横坐标为tm(1),两点连线,计算其斜率S1;
S12,计算起始端点处延展得到的极小值N(0);
N(0)=N(1)-S1·[tn(1)-tn(0)] (2)
N(1)指的是:距离起始端点处最近的极小值,tn(1)指的是:N(1)的横坐标,tn(0)=t(0);
S13,相应的,把原始脑电信号的终端端点作为极小值点记为N(Qmin),横坐标为tn(Qmin),选取一个相邻的极小值点记为N(Qmin-1),横坐标为tn(Qmin-1),采用公式(3)计算两点连线的斜率S2:
S14,计算终端端点延展得到的极大值M(Qmax):
M(Qmax)=M(Qmax-1)-S2·[tm(Qmax-1)-t(N)] (4)
选取一个相邻的极大值点记为M(Qmax-1),横坐标为tm(Qmax-1),tm(Qmax)=t(N);
S15,重复步骤S14,由此找出原始脑电信号数据序列X(t)内的所有的极大值点M1(Qmax)和极小值点N1(0),得到极大值序列;
找出极大值序列中的最大值MAX和最小值MIN,若M(Qmax)<MIN,则定义该延展极大值为:M1(Qmax)=(终端端点值+MAX)/2;若N(0)>MAX,则定义该延展极大值为:N1(0)=(起始端点值+MIN)/2;
S16,采用三次样条插值函数拟合形成原数据的上包络线Umax(t)和下包络线Umin(t);
根据上、下包络线计算原始数据X(t)的局部平均值,记作公式(5):
m1(t)=(Umax(t)+Umin(t))/2 (5)
S17,计算原始脑电信号与局部平均值的差值h1(t),记为公式(6):
h1(t)=X(t)-m1(t) (6)
S18,判断h1(t)是否满足IMF的两个条件,若满足,则h1(t)为第一个IMF分量,从原始信号X(t)中减去h1(t)得到剩余信号r1(t);否则就把h1(t)当成原始脑电信号转到步骤S16,重复步骤S16-S17,继续寻找原数据序列h1(t)信号内的所有的极大值点M1(Qmax)和极小值点N1(0),直至满足IMF的两个条件
r1(t)=X(t)-h1(t) (7)
S19,将剩余信号r1(t)看作新的脑电信号重复步骤S11-S18,可筛选得到n个IMF分量:
rn(t)=rn-1(t)-hn(t) (8)
S110,对n个IMF分量h1(t),h2(t)…hn(t)做阈值处理,将大于阈值V的分量置为零,如下式所示:
其中,其中,F为原始脑电信号数据X(t)的长度;
Median()函数,返回一组已知数字的中值,中值是一组数的中间数,i的值为1~n;
步骤S2具体包括:
S21,采用下式,对重构脑电信号去均值,得到处理后的脑电观测信号X’(t);
其中:E[·]表示数学期望;
S22,采用下式,将处理后的脑电观测信号X’(t),分解为各分量间互不相关的脑电信号Z(t)=(z1(t),L,zm(t))T:
其中:代表投影因子;
Ds是以X’(t)的协方差矩阵CX=E[X’(t)*X’(t)T]特征值为对角元素的对角矩阵;
Us是以CX的单位范数特性向量为列的矩阵;
I是单位矩阵;
σ表示X’(t)的噪声方差,
S23,定义Y(t)=WZ(t),W为解混矩阵;
初始化各参数:令循环变量i=1,初始解混矩阵为:
E[·]表示数学期望;
S24,计算Y(t)=WZ(t);
ΔW=λ[I-K tanh(Y(t))(Y(t))T-Y(t)(Y(t))T]W
其中,λ为学习步长,本文取λ=0.001,I为单位矩阵;K为对角矩阵,其对角元素kii=sgn{E[sech2 yi(t)]-E[yi(t)tanh yi(t)]}
Sgn()符号函数,返回参数的正负,sech()为双曲正割函数,tanh()为双曲正切函数;
S25,W=W+ΔW
返回步骤S24继续计算,直到W收敛得到独立分量Y(t)=WZ(t);处理后的解混矩阵W的无穷范数是否小于10-6,如果是,则代表收敛。
2.根据权利要求1所述的基于改进EMD-ICA的EEG信号处理方法,其特征在于,满足IMF的两个条件具体为:1)在整个数据段内,极值点和零点数量相等或相差至多1个;2)在任意时刻,由局部极大值点形成的上包络和由局部极小值点形成的下包络的平均值为0。
3.根据权利要求1所述的基于改进EMD-ICA的EEG信号处理方法,其特征在于,步骤S110后还包括:分解后的IMF分量和余量rn叠加重构为去噪后的重构脑电信号
4.根据权利要求1所述的基于改进EMD-ICA的EEG信号处理方法,其特征在于,步骤S3中绘制二维时频图的具体过程包括:
S31,从C3和C4通道的脑电信号中提取时间、频率和位置信息;
S32,利用连续小波转换将步骤S31中获得的时间、频率和位置信息结合绘制出二维时频图。
5.根据权利要求1所述的基于改进EMD-ICA的EEG信号处理方法,其特征在于,步骤S4中采用的卷积神经网络模型采用2层卷积神经网络模型,具体包括2个卷积层,2个池化层,2个全连接层及1个softmax输出层。
6.根据权利要求5所述的基于改进EMD-ICA的EEG信号处理方法,其特征在于,步骤S4具体包括:
S41,将绘制的二维时频图首先传入第一卷积层,提取局部特征
S42,将提取局部特征后的特征图传入第一池化层,经过特征选择和信息过滤后再依次传入第二卷积层和第二池化层;
S43,进入全连接层,最后进入输出层,采用softmax分类函数即可得到左右手的分类。
7.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质包括处理器,所述处理器用于实现权利要求1-6任一所述的基于改进EMD-ICA的EEG信号处理方法。
8.一种基于改进EMD-ICA的EEG信号处理设备,其特征在于,用于执行权利要求1-6任一所述的基于改进EMD-ICA的EEG信号处理方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111614480.7A CN114224360B (zh) | 2021-12-27 | 2021-12-27 | 一种基于改进emd-ica的eeg信号处理方法、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111614480.7A CN114224360B (zh) | 2021-12-27 | 2021-12-27 | 一种基于改进emd-ica的eeg信号处理方法、设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114224360A CN114224360A (zh) | 2022-03-25 |
CN114224360B true CN114224360B (zh) | 2023-10-10 |
Family
ID=80763485
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111614480.7A Active CN114224360B (zh) | 2021-12-27 | 2021-12-27 | 一种基于改进emd-ica的eeg信号处理方法、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114224360B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116738221B (zh) * | 2023-08-15 | 2023-10-20 | 湖南天联城市数控有限公司 | 一种带压管道气体分析方法及系统 |
CN117349603B (zh) * | 2023-12-06 | 2024-03-12 | 小舟科技有限公司 | 脑电信号的自适应降噪方法及装置、设备、存储介质 |
CN117390373B (zh) * | 2023-12-13 | 2024-05-17 | 广东企禾科技有限公司 | 一种通信传输设备调测维修管理方法及系统 |
CN117860273B (zh) * | 2024-03-13 | 2024-06-11 | 小舟科技有限公司 | 脑电信号智能降噪方法、系统、设备及介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108042132A (zh) * | 2017-12-27 | 2018-05-18 | 南京邮电大学 | 基于dwt和emd融合csp的脑电特征提取方法 |
CN108888264A (zh) * | 2018-05-03 | 2018-11-27 | 南京邮电大学 | Emd和csp融合功率谱密度脑电特征提取方法 |
CN109924990A (zh) * | 2019-03-27 | 2019-06-25 | 兰州大学 | 一种基于emd算法的脑电信号抑郁症识别系统 |
CN110163128A (zh) * | 2019-05-08 | 2019-08-23 | 南京邮电大学 | 改进的emd算法结合小波包变换及csp算法的脑电信号分类方法 |
CN111413095A (zh) * | 2020-04-13 | 2020-07-14 | 西安电子科技大学 | 基于瞬时角速度的行星轴承分布式故障诊断分析方法 |
CN112107310A (zh) * | 2020-09-30 | 2020-12-22 | 西安理工大学 | 基于iwt与aga-bp模型的ecg身份识别方法 |
WO2021253094A1 (en) * | 2020-06-19 | 2021-12-23 | The Florey Institute Of Neuroscience And Mental Health | System and device for quantifying motor control disorder |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7108192B2 (en) * | 1999-09-17 | 2006-09-19 | Silverbrook Research Pty Ltd | Rotationally symmetric tags |
US8321017B2 (en) * | 2009-07-08 | 2012-11-27 | Pacesetter, Inc. | Electromechanical delay (EMD) monitoring devices, systems and methods |
WO2011127609A1 (en) * | 2010-04-16 | 2011-10-20 | Applied Brain And Vision Sciences Inc. | Encephalography method and apparatus incorporating independent component analysis and a spectral shaping filter |
US10606353B2 (en) * | 2012-09-14 | 2020-03-31 | Interaxon Inc. | Systems and methods for collecting, analyzing, and sharing bio-signal and non-bio-signal data |
EP3424418B1 (en) * | 2017-07-05 | 2023-11-08 | Stichting IMEC Nederland | A method and a system for detecting a vital sign of a subject |
US11406316B2 (en) * | 2018-02-14 | 2022-08-09 | Cerenion Oy | Apparatus and method for electroencephalographic measurement |
WO2020112877A1 (en) * | 2018-11-28 | 2020-06-04 | Easyg Llc | Systems and methods for digitally processing biopotential signal |
FR3095546A1 (fr) * | 2019-04-26 | 2020-10-30 | Thales | Dispositif et procédé d’analyse de l’état d’un système en contexte bruité |
-
2021
- 2021-12-27 CN CN202111614480.7A patent/CN114224360B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108042132A (zh) * | 2017-12-27 | 2018-05-18 | 南京邮电大学 | 基于dwt和emd融合csp的脑电特征提取方法 |
CN108888264A (zh) * | 2018-05-03 | 2018-11-27 | 南京邮电大学 | Emd和csp融合功率谱密度脑电特征提取方法 |
CN109924990A (zh) * | 2019-03-27 | 2019-06-25 | 兰州大学 | 一种基于emd算法的脑电信号抑郁症识别系统 |
CN110163128A (zh) * | 2019-05-08 | 2019-08-23 | 南京邮电大学 | 改进的emd算法结合小波包变换及csp算法的脑电信号分类方法 |
CN111413095A (zh) * | 2020-04-13 | 2020-07-14 | 西安电子科技大学 | 基于瞬时角速度的行星轴承分布式故障诊断分析方法 |
WO2021253094A1 (en) * | 2020-06-19 | 2021-12-23 | The Florey Institute Of Neuroscience And Mental Health | System and device for quantifying motor control disorder |
CN112107310A (zh) * | 2020-09-30 | 2020-12-22 | 西安理工大学 | 基于iwt与aga-bp模型的ecg身份识别方法 |
Non-Patent Citations (3)
Title |
---|
Bono, V ; Das, S ; Jamal, W ; 等.Hybrid wavelet and EMD/ICA approach for artifact suppression in pervasive EEG.《JOURNAL OF NEUROSCIENCE METHODS》.2016,全文. * |
张碧薇.基于EEMD与平稳小波变换的脉搏波形特征分析研究.《中国优秀硕士学位论文全文数据库信息科技辑》.2014,全文. * |
耿晓中 ; 李得志.基于独立分量分析的眼电伪迹去除方法研究.《长春工程学院学报(自然科学版)》.2020,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN114224360A (zh) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114224360B (zh) | 一种基于改进emd-ica的eeg信号处理方法、设备及存储介质 | |
CN107844755B (zh) | 一种结合dae和cnn的脑电信号特征提取与分类方法 | |
CN111310570B (zh) | 一种基于vmd和wpd的脑电信号情感识别方法及系统 | |
Roy et al. | Image denoising by data adaptive and non-data adaptive transform domain denoising method using EEG signal | |
WO2018120088A1 (zh) | 情感识别模型生成方法及装置 | |
CN114533086A (zh) | 一种基于空域特征时频变换的运动想象脑电解码方法 | |
CN115211869A (zh) | 一种基于经验模态分解和卡尔曼滤波的脑电信号去噪方法、装置及系统 | |
CN113486752A (zh) | 基于心电信号的情感识别方法及系统 | |
Mishra et al. | Noise removal in EEG signals using SWT–ICA combinational approach | |
Prakash et al. | Noise reduction of ecg using chebyshev filter and classification using machine learning algorithms | |
Asghar et al. | Semi-skipping layered gated unit and efficient network: hybrid deep feature selection method for edge computing in EEG-based emotion classification | |
Geng et al. | [Retracted] A Fusion Algorithm for EEG Signal Processing Based on Motor Imagery Brain‐Computer Interface | |
Hai et al. | Wavelet-based kernel construction for heart disease classification | |
CN115017960B (zh) | 一种基于时空联合mlp网络的脑电信号分类方法及应用 | |
CN115462803A (zh) | 一种基于BG-Attention的脑电信号去噪方法、装置及存储介质 | |
CN116098637A (zh) | 一种基于ica优化校正脑电微状态的大脑功能评测装置 | |
CN114795247A (zh) | 脑电信号分析方法、装置、电子设备和存储介质 | |
CN114692682A (zh) | 一种基于图嵌入表示的运动想象分类方法及系统 | |
CN114757236A (zh) | 基于tqwt与svmd的脑电信号去噪优化方法及系统 | |
Wang et al. | Improved brain–computer interface signal recognition algorithm based on few-channel motor imagery | |
CN111008590B (zh) | 基于多维特征提取无监督eeg信号分类方法 | |
CN111493864B (zh) | Eeg信号混合噪声处理方法、设备和存储介质 | |
Dong et al. | An Approach for EEG Denoising Based on Wasserstein Generative Adversarial Network | |
Mir et al. | Adaptive data analysis methods for biomedical signal processing applications | |
Srinivasulu et al. | Basis pursuit sparse decomposition using tunable-Q wavelet transform (BPSD-TQWT) for denoising of electrocardiograms |
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 |