CN112006682B - 一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法 - Google Patents
一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法 Download PDFInfo
- Publication number
- CN112006682B CN112006682B CN202010491390.2A CN202010491390A CN112006682B CN 112006682 B CN112006682 B CN 112006682B CN 202010491390 A CN202010491390 A CN 202010491390A CN 112006682 B CN112006682 B CN 112006682B
- Authority
- CN
- China
- Prior art keywords
- frequency
- data
- channel
- electroencephalogram
- user
- 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
Images
Classifications
-
- 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- 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)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Power Engineering (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明设计了一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法,该系统的硬件设备包括脑电信号采集装置和显示屏;软件使用Anaconda3开源平台Python3.6编程环境。系统实现方法如下:受试者佩戴采集装置并跟显示屏指示进行运动想象任务,采集脑电信号;对实验数据进行滤波和独立成分分析,去除噪声;提取三组主要电极信号,分别使用Multitaper的方法计算其功率谱密度,用于时间‑频率二维特征的分类;选择二维分类准确率大于75%的频率特征,输入SVM进行调参、训练,得到的分类器用于左右手运动想象脑电信号的在线识别。本发明针对用户提取高质量的频率特征,去除冗余信息,提高识别效率,旨在应用于脑卒中等神经疾病的智能康复机器人的脑电识别中。
Description
技术领域
本发明属于运动想象脑电信号识别领域,尤其涉及一种基于多通道频率特征定制的左右手运动想象脑电信号分类系统。
背景技术
脑卒中,俗称“中风”,是最典型的脑血管病之一。调查显示,世界受中风影响的绝对人数正在上升。而我国脑卒中的患病率一直高于全球,是国民第一死亡原因。这些患者通过正常的大脑输出途径来完成与外界环境的交流十分困难。因此,为了提高脑卒中患者的生活质量,减轻疾病对家庭与社会的负担,我们期望可以通过人工智能技术辅助并提高脑卒中等手部偏瘫患者的康复治疗。
脑机接口(Brain-Computer-Interface,BCI)是目前国际上广泛研究的一种脑卒中康复方式,使用BCI技术的机器人康复设备,利用脑信号采集装置和智能算法自动检测人脑主动活动的电信号,并通过外接机电设备对人体相应的部位产生刺激或引导活动,人工建立起神经-肌肉的运动通路,刺激原神经回路的重建。这种方式通过提高对剩余脑回路的利用辅助运动功能的恢复与改善。采集头皮表面脑电信号的EEG(Electroencephalography)技术是一种比较方便、安全和低廉的BCI技术。基于运动想象(Motor Imagery,MI)的EEG-BCI系统由于其依赖于使用者想象肢体的运动动作,涉及了实际运动的控制机制和神经基质,可以促进大脑运动感觉皮层(Sensorimotor cortex,SM1)的激活,尤其适用于脑卒中患者的手部康复治疗。实现一个运动想象分类系统包括信号采集设计、脑电信号数据预处理、脑电信号特征提取以及分类算法设计等关键技术。但由于EEG信号的非稳态特性且存在个性化差异的问题,目前大多数方法忽略了用户数据的多样性,往往在实验中达到的高效率无法推广到实际应用中。
发明内容
针对现有技术存在的以上问题,本发明旨在设计一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法,通过将相关通道去噪信号的功率谱密度作为特征在时频二维平面上分区,并分别进行SVM分类,根据分类的交叉验证得分来评估特征区的分类贡献度,并选择贡献度大于阈值0.75的通道频率作为定制频率,相应的通道频率特征作为定制特征用于SVM分类器的训练和在线识别。由此改善脑电信号的个性化差异对运动想象分类在实际应用中的影响。
为实现本发明的目的,本发明所采用的技术方案如下:一种基于多通道频率特征定制的左右手运动想象脑电信号分类系统,具体包括:通过FIR滤波器和独立成分分析去除采集信号的噪声,对其进行epochs数据段截取后,提取10-20电极分布系统中FC(FC1、FC2和Cz)、C(C3、C4和Cz)和CP(CP1、CP2和Cz)三组电极通道的信号数据,并使用Multitaper与Morlet小波结合的方法计算其在分析频率段FOI的功率谱密度,作为信号特征在FOI-epochs时频平面上进行分区;分别做SVM分类得到每个特征区的分类得分用来衡量其分类贡献度,标记贡献度大于0.75的通道频率作为定制通道频率,并将定制通道频率特征输入SVM分类器进行调参和模型训练。在线识别直接提取定制通道频率特征输入训练好的SVM分类器,得到识别结果。
为了实现上述目的,本发明的技术方案如下:一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法,该方法包括以下几个步骤:
步骤1:采集用户的脑电信号;
步骤1.1:开启脑电数据采集设备,将电脑USB接口连接数据采集设备的Digital口,打开电脑实验范式程序(即实验引导GUI);程序运行过程中,电脑通过USB接口发送两种串口信号给数据采集设备,标记两种想象状态的触发时间点;
步骤1.2:用户正确佩戴上符合10-20电极分布系统的脑电帽,并在电极腔中注射导电膏,使得电极导通(电阻≤5KΩ);在用户手腕、眼眶上下分别贴上电极片,作为参考电极记录伪迹信号;
步骤1.3:采集实验中,用户须端正地坐在扶手靠椅上,身体放松,双手自然搭在扶手上方,掌心向上;实验过程中须注视距离其头部约40-60cm远的显示屏,按照屏幕指令操作;当用户准备好时,即可开始实验;一轮(Run)实验包含40次想象尝试(Trials),每次尝试包含1个手部运动想象任务,左/右手握拳运动想象任务各20个,出现次序随机;执行运动想象任务之前,屏幕上呈现3s“十”字提示用户集中注意力;之后屏幕上显示1s“准备”并伴有“滴”提示音;提示音停止随即开始4s的运动想象任务,此时屏幕的提示为“十”字上添加向左或向右的箭头,表示受试者进行左/右手运动想象任务;运动想象任务结束之后是1.5s的休息时间;在进行下一次尝试之前,还有小于1s的随机时间间隙;离线训练需要每位用户完成3轮数据采集,共120个trials;
步骤2:对采集数据进行预处理,即FIR滤波和独立成分分析(ICA);
步骤2.1:使用窗函数为Hamming窗的有限冲击响应(FIR)滤波器分别对采集到的原始数据进行1Hz高通滤波和45Hz的低通滤波,其中低通过渡带为min(max(l_freq*0.25,2),l_freq),高通过渡带为这里的l_freq、h_freq和sfreq分别表示下通带频率、上通带频率和采样频率;滤波器长度设定为最短过渡带倒数的6.6倍;这里可以去除电源噪声(50Hz或60Hz)和高频率的肌电伪迹(约60Hz);
步骤2.2:滤波后,使用独立成分分析(ICA)的FastICA算法对32个脑电(EEG)通道的信号提取10个独立成分,将10个独立成分与眼部伪迹电极信号中的眼动成分进行相关性分析,设定阈值2.5,自适应z得分高于阈值的独立成分被认为眼动伪迹并剔除,将剩下的独立成分逆映射回电极空间;至此,预处理完成,将得到相对洁净的数据;
步骤3:对步骤2中得到的去噪数据根据同步信号的标记点分割出运动想象开始前2s至想象任务结束后2s时间内的数据段epochs,并提取出位于大脑运动感觉区的FC、C和CP三组电极通道的epoch数据,每组数据都包括10-20电极分布系统中的左右两个电极(FC1和FC2、C3和C4、CP1和CP2)和中心电极Cz;对三组数据计算PSD,这里使用的Multitaper算法其实是Morlet小波与倍数DPSS(Discrete Prolate Spheroidal)Multitaper窗相结合的;令采样周期为Δ,则算法的具体步骤如下:
(1)选择F=[5Hz,37Hz]作为分析频率段(Frequencies of Interest,FOI),每个频率段Morlet小波的长度ncycles设定为7.0,就是说对于频率f∈F,小波窗(和taper向量)的长度为
(2)对每个频率f∈F求DPSS序列:对ΦR进行主成分分析(PCA),ΦR为Nf×Nf的矩阵,ΦR的每个元素
所以L取3,则得到的3个主成分即为3个相互正交的tapers向量,它们组成的序列就是DPSS,或称Slepian序列;
(3)将每个频率f∈F且长度为Nf的Morlet小波函数分别与DPSS序列中每个taper向量做点乘,得到关于每个频率f的Morlet-DPSS序列,序列维度为(3×33×Nf);
(4)对每个通道的每段epoch数据的快速傅里叶变换(FFT)与Morlet-DPSS序列中每个taper向量的快速傅里叶变换相乘并做快速傅里叶逆变换(iFFT),得到每个频率f的Multitaper系数(C表示通道序号,k表示epoch序号);
至此,便得到三组数据在FOI-epoch时间频率二维平面上的功率谱密度PSD分布;
步骤4:对步骤3中得到的三张FOI-epoch时频二维特征平面,将每组通道数据的特征分为面积为1Hz×0.125s的2112个bins,共6336个bins;每一个bin的特征分别输入SVM分类器中进行训练,并使用5×5fold cross-validation验证每个bin特征的分类AUC得分;这里,SVM核函数为径向基函数RBF,n_features表示特征数量,当采样频率sfreq=1000Hz时,n_features=125,惩罚因子C设定为15;通过上述过程,可以得到三组通道数据在FOI-epoch特征平面的SVM分类得分分布Feature Scores,其维度为(3×33×64);
步骤5:对步骤4中得到的三组通道数据FOI-epoch特征平面的SVM分类得分Feature Scores设定阈值thresholdscores=0.75,筛选出每组通道数据大于阈值的通道频率特征,组成定制通道频率特征向量Z0,并标记定制的通道特征频率fo;可以注意到,由于在线识别无法对时间特征进行选择,所以这里仅仅需要提取符合要求的频率特征;令定制通道频率特征Z为输入,以RBF为核函数的SVM的参数(C,gamma)作为调参对象,选择在5×5fold cross-validation中AUC得分最高的一组参数(Co,gammao),并将此时训练的SVM模型作为在线识别的分类器;
步骤6:离线训练模型用于在线识别,即实时采集用户脑电数据,依次进行预处理、提取定制通道频率fo的数据并计算通道频率特征Z、输入SVM模型进行识别,最终得到识别结果反馈给用户;
步骤6.1:用户佩戴好采集装置后,跟随显示屏指示进行左右手握拳的运动想象;电脑实时接收采集设备输送的脑电信号,并在接收到触发信号后,对每个采样时刻的数据进行FIR高通(1Hz)和低通(45Hz)滤波后,进行独立成分分析剔除伪迹噪声;
步骤6.2:对去噪信号进行定制特征提取和分类识别操作:
(1)根据定制通道特征频率fo,对特定电极通道的信号做带通滤波,中心频率为f∈fo,带宽为1Hz,提取全部通道特征频率的数据;
(2)对(1)中提取的数据使用Multitaper算法计算PSD特征,即定制通道频率特征Z;
(3)将Z输入SVM分类器识别想象内容;
(4)将识别结果输出给反馈装置,形成闭环。
与现有技术相比,本发明的技术方案具有以下的有益技术效果:
1、本发明提出了一种在运动想象最相关的三组通道(10-20电极分布系统中的FC(FC1、FC2和Cz)、C(C3、C4和Cz)和CP(CP1、CP2和Cz))上做PSD特征提取,能够有效消除信号的冗余,降低数据计算量,提高在线识别的速度;
2、针对运动想象脑电信号的个性化差异问题,对三组通道的信号特征在时频平面上作分区,对特征区的分类贡献打分,筛选得分大于0.75的特征频率作为定制通道频率,仅计算该频率范围内的特征,使得特征的维数降低,但质量提高,在提高在线识别的速度同时,改善在线识别的效果;
3、不同用户可以得到不同的定制通道频率,但由于阈值的限定,定制通道频率特征的质量不会受到用户差异的影响。因此,在兼顾脑电信号个性化差异的同时,也能保证系统分类的效果。
4、定制特征从通道和频率两方面出发,与离线训练的数据无关,避免了很多现有技术处理脑电信号个性化问题时由于过分关注训练数据而导致的过拟合问题,本方法具有较好且稳定的泛化能力。
附图说明
图1是本发明一种基于多通道频率特征定制的左右手运动想象脑电信号分类系统的框架示意图;
图2是本发明采用通道频率特征定制的方法流程图。
具体实施方式
下面结合附图和实施例对本发明的技术方案作进一步的说明。
实施例1:如图1所示,本发明提出了一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法,该方法的详细步骤为:
步骤1:采集用户的脑电信号;
步骤1.1:开启脑电数据采集设备,将电脑USB接口连接数据采集设备的Digital口,打开电脑实验范式程序(即实验引导GUI)。程序运行过程中,电脑通过USB接口发送两种串口信号给数据采集设备,标记两种想象状态的触发时间点。
步骤1.2:用户正确佩戴上符合10-20电极分布系统的脑电帽,并在电极腔中注射导电膏,使得电极导通(电阻≤5KΩ)。在用户手腕、眼眶上下分别贴上电极片,作为参考电极记录伪迹信号。
步骤1.3:采集实验中,用户须端正地坐在扶手靠椅上,身体放松,双手自然搭在扶手上方,掌心向上。实验过程中须注视距离其头部约40-60cm远的显示屏,按照屏幕指令操作。当用户准备好时,即可开始实验。一轮(Run)实验包含40次想象尝试(Trials),每次尝试包含1个手部运动想象任务,左/右手握拳运动想象任务各20个,出现次序随机;执行运动想象任务之前,屏幕上呈现3s“十”字提示用户集中注意力;之后屏幕上显示1s“准备”并伴有“滴”提示音;提示音停止随即开始4s的运动想象任务,此时屏幕的提示为“十”字上添加向左或向右的箭头,表示受试者进行左/右手运动想象任务;运动想象任务结束之后是1.5s的休息时间;在进行下一次尝试之前,还有小于1s的随机时间间隙。离线训练需要每位用户完成3轮数据采集,共120个trials。
步骤2:对采集数据进行预处理,即FIR滤波和独立成分分析(ICA);
步骤2.1:使用窗函数为Hamming窗的有限冲击响应(FIR)滤波器分别对采集到的原始数据进行1Hz高通滤波和45Hz的低通滤波,其中低通过渡带为min(max(l_freq*0.25,2),l_freq),高通过渡带为这里的l_freq、h_freq和sfreq分别表示下通带频率、上通带频率和采样频率;滤波器长度设定为最短过渡带倒数的6.6倍。这里可以去除电源噪声(50Hz或60Hz)和高频率的肌电伪迹(约60Hz)
步骤2.2:滤波后,使用独立成分分析(ICA)的FastICA算法对32个脑电(EEG)通道的信号提取10个独立成分,将10个独立成分与眼部伪迹电极信号中的眼动成分进行相关性分析,设定阈值2.5,自适应z得分高于阈值的独立成分被认为眼动伪迹并剔除,将剩下的独立成分逆映射回电极空间。至此,预处理完成,将得到相对洁净的数据。
提取10-20电极分布系统中的三组通道数据,即电极FC1、FC2和Cz、电极C3、C4和Cz,以及电极CP1、CP2和Cz,使用Multitaper方法计算信号功率谱密度(PSD);
步骤3:对步骤(2)中得到的去噪数据根据同步信号的标记点分割出运动想象开始前2s至想象任务结束后2s时间内的数据段epochs,并提取出位于大脑运动感觉区的FC、C和CP三组通道的epoch数据,每组数据都包括左右两个电极(FC1和FC2、C3和C4、CP1和CP2)和中心电极Cz。对三组数据计算PSD,这里使用的Multitaper算法其实是Morlet小波与倍数DPSS(Discrete Prolate Spheroidal)Multitaper窗相结合的。令采样周期为Δ,则算法的具体步骤如下:
1)选择F=[5Hz,37Hz]作为分析频率段(Frequencies of Interest,FOI),每个频率段Morlet小波的长度ncycles设定为7.0,就是说对于频率f∈F,小波窗(和taper向量)的长度为
2)对每个频率f∈F求DPSS序列:对ΦR进行主成分分析(PCA),ΦR为Nf×Nf的矩阵,ΦR的每个元素
所以L取3,则得到的3个主成分即为3个相互正交的tapers向量,它们组成的序列就是DPSS,或称Slepian序列;
3)将每个频率f∈F且长度为Nf的Morlet小波函数分别与DPSS序列中每个taper向量做点乘,得到关于每个频率f的Morlet-DPSS序列,序列维度为(3×33×Nf);
4)对每个通道的每段epoch数据的快速傅里叶变换(FFT)与Morlet-DPSS序列中每个taper向量的快速傅里叶变换相乘并做快速傅里叶逆变换(iFFT),得到每个频率f的Multitaper系数(C表示通道序号,k表示epoch序号);
至此,便得到三组数据在FOI-epoch时间频率二维平面上的功率谱密度PSD分布。
将三组通道信号的PSD在时-频二维平面上进行分区,并逐一输入SVM分类器,得到PSD特征在时-频平面上的分类得分分布;
步骤4:对步骤(3)中得到的三张FOI-epoch时频二维特征平面,将每组通道数据的特征分为面积为1Hz×0.125s的2112个bins,共6336个bins。每一个bin的特征分别输入SVM分类器中进行训练,并使用5×5fold cross-validation验证每个bin特征的分类AUC得分。这里,SVM核函数为径向基函数RBF,n_features表示特征数量,当采样频率sfreq=1000Hz时,n_features=125,惩罚因子C设定为15。
通过上述过程,可以得到三组通道数据在FOI-epoch特征平面的SVM分类得分分布Feature Scores,其维度为(3×33×64)。
设定得分阈值thresholdscores=0.75,获得分类贡献较高的一批通道频率特征Z0,并标记定制的通道特征频率fo,将Z0输入SVM分类器进行调参,训练得到最终的SVM分类器模型;
步骤5:对步骤(4)中得到的三组通道数据FOI-epoch特征平面的SVM分类得分Feature Scores设定阈值thresholdscores=0.75,筛选出每组通道数据大于阈值的通道频率特征,组成定制通道频率特征向量Z0,并标记定制的通道特征频率fo。可以注意到,由于在线识别无法对时间特征进行选择,所以这里仅仅需要提取符合要求的频率特征。具体操作过程由图2所示。
令定制通道频率特征Z为输入,以RBF为核函数的SVM的参数(C,gamma)作为调参对象,选择在5×5fold cross-validation中AUC得分最高的一组参数(Co,gammao),并将此时训练的SVM模型作为在线识别的分类器。
步骤6:离线训练模型用于在线识别,即实时采集用户脑电数据,依次进行预处理、提取定制通道频率fo的数据并计算通道频率特征Z、输入SVM模型进行识别,最终得到识别结果反馈给用户。
步骤6.1:用户佩戴好采集装置后,跟随显示屏指示进行左右手握拳的运动想象。电脑实时接收采集设备输送的脑电信号,并在接收到触发信号后,对每个采样时刻的数据进行FIR高通(1Hz)和低通(45Hz)滤波后,进行独立成分分析剔除伪迹噪声。
步骤6.2:对去噪信号进行定制特征提取和分类识别操作:
1)根据定制通道特征频率fo,对特定电极通道的信号做带通滤波,中心频率为f∈fo,带宽为1Hz,提取全部通道特征频率的数据;
2)对(1)中提取的数据使用Multitaper算法计算PSD特征,即定制通道频率特征Z;
3)将Z输入SVM分类器识别想象内容;
4)将识别结果输出给反馈装置,形成闭环。
需要说明的是上述实施例仅仅是本发明的较佳实施例,并没有用来限定本发明的保护范围,在上述技术方案的基础上做出的等同替换或者替代均属于本发明的保护范围。
Claims (5)
1.一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法,其特征在于:该方法包括以下几个步骤:
步骤1:采集用户的脑电信号;
步骤2:对采集数据进行预处理,即FIR滤波和独立成分分析(ICA);
步骤3:对步骤2中得到的去噪数据根据同步信号的标记点分割出运动想象开始前2s至想象任务结束后2s时间内的数据段epochs,并提取出位于大脑运动感觉区的FC、C和CP三组电极通道的epoch数据,每组数据都包括10-20电极分布系统中的左右两个电极(FC1和FC2、C3和C4、CP1和CP2)和中心电极Cz;使用Multitaper算法对三组数据计算PSD,这里使用的Multitaper算法其实是Morlet小波与倍数DPSS(Discrete Prolate Spheroidal)Multitaper窗相结合的;
步骤4:对步骤3中得到的三张FOI-epoch时频二维特征平面,将每组通道数据的特征分为面积为1Hz×0.125s的2112个bins,共6336个bins;每一个bin的特征分别输入SVM分类器中进行训练,并使用5×5fold cross-validation验证每个bin特征的分类AUC得分;这里,SVM核函数为径向基函数RBF,n_features表示特征数量,当采样频率sfreq=1000Hz时,n_features=125,惩罚因子C设定为15;由此得到三组通道数据在FOI-epoch特征平面的SVM分类得分分布Feature Scores,其维度为(3×33×64);
步骤5:对步骤4中得到的三组通道数据FOI-epoch特征平面的SVM分类得分FeatureScores设定阈值thresholdscores=0.75,筛选出每组通道数据大于阈值的通道频率特征,组成定制通道频率特征向量Z0,并标记定制的通道特征频率fo;令定制通道频率特征Z0为输入,以RBF为核函数的SVM的参数(C,gamma)作为调参对象,选择在5×5fold cross-validation中AUC得分最高的一组参数(Co,gammao),并将此时训练的SVM模型作为在线识别的分类器;
步骤6:离线训练模型用于在线识别,即实时采集用户脑电数据,依次进行预处理、提取定制通道频率fo的数据并计算通道频率特征Z、输入SVM模型进行识别,最终得到识别结果反馈给用户。
2.根据权利要求1所述的基于多通道频率特征定制的左右手运动想象脑电信号分类方法,其特征在于:所述步骤1采集用户的脑电信号;具体如下:
步骤1.1:开启脑电数据采集设备,将电脑USB接口连接数据采集设备的Digital口,打开电脑实验范式程序(即实验引导GUI);程序运行过程中,电脑通过USB接口发送两种串口信号给数据采集设备,标记两种想象状态的触发时间点;
步骤1.2:用户正确佩戴上符合10-20电极分布系统的脑电帽,并在电极腔中注射导电膏,使得电极导通;在用户手腕、眼眶上下分别贴上电极片,作为参考电极记录伪迹信号;
步骤1.3:采集实验中,用户须端正地坐在扶手靠椅上,身体放松,双手自然搭在扶手上方,掌心向上;实验过程中须注视距离其头部约40-60cm远的显示屏,按照屏幕指令操作;当用户准备好时,即可开始实验;一轮(Run)实验包含40次想象尝试(Trials),每次尝试包含1个手部运动想象任务,左/右手握拳运动想象任务各20个,出现次序随机;执行运动想象任务之前,屏幕上呈现3s“十”字提示用户集中注意力;之后屏幕上显示1s“准备”并伴有“滴”提示音;提示音停止随即开始4s的运动想象任务,此时屏幕的提示为“十”字上添加向左或向右的箭头,表示受试者进行左/右手运动想象任务;运动想象任务结束之后是1.5s的休息时间;在进行下一次尝试之前,还有小于1s的随机时间间隙;离线训练需要每位用户完成3轮数据采集,共120个trials。
3.根据权利要求1所述的基于多通道频率特征定制的左右手运动想象脑电信号分类方法,其特征在于:步骤2:对采集数据进行预处理,即FIR滤波和独立成分分析(ICA);具体如下:
步骤2.1:使用窗函数为Hamming窗的有限冲击响应(FIR)滤波器分别对采集到的原始数据进行1Hz高通滤波和45Hz的低通滤波,其中低通过渡带为min(max(l_freq*0.25,2),l_freq),高通过渡带为这里的l_freq、h_freq和sfreq分别表示下通带频率、上通带频率和采样频率;滤波器长度设定为最短过渡带倒数的6.6倍;这里能够去除电源噪声和高频率的肌电伪迹;
步骤2.2:滤波后,使用独立成分分析(ICA)的FastICA算法对32个脑电(EEG)通道的信号提取10个独立成分,将10个独立成分与眼部伪迹电极信号中的眼动成分进行相关性分析,设定阈值2.5,自适应z得分高于阈值的独立成分被认为眼动伪迹并剔除,将剩下的独立成分逆映射回电极空间;至此,预处理完成,将得到相对洁净的数据。
4.根据权利要求1所述的基于多通道频率特征定制的左右手运动想象脑电信号分类方法,其特征在于:步骤3中Multitaper算法其实是Morlet小波与倍数DPSS(DiscreteProlate Spheroidal)Multitaper窗相结合的;令采样周期为Δ,则算法的具体步骤如下:
(1)选择F=[5Hz,37Hz]作为分析频率段(Frequencies of Interest,FOI),每个频率段Morlet小波的长度ncycles设定为7.0,就是说对于频率f∈F,小波窗和taper向量的长度为
(2)对每个频率f∈F求DPSS序列:对ΦR进行主成分分析(PCA),ΦR为Nf×Nf的矩阵,ΦR的每个元素
所以L取3,则得到的3个主成分即为3个相互正交的tapers向量,它们组成的序列就是DPSS;
(3)将每个频率f∈F且长度为Nf的Morlet小波函数分别与DPSS序列中每个taper向量做点乘,得到关于每个频率f的Morlet-DPSS序列,序列维度为(3×33×Nf);
(4)对每个通道的每段epoch数据的快速傅里叶变换(FFT)与Morlet-DPSS序列中每个taper向量的快速傅里叶变换相乘并做快速傅里叶逆变换(iFFT),得到每个频率f的Multitaper系数C表示通道序号,k表示epoch序号;
5.根据权利要求1所述的基于多通道频率特征定制的左右手运动想象脑电信号分类方法,其特征在于:所述步骤6离线训练模型用于在线识别,具体如下:
步骤6.1:用户佩戴好采集装置后,跟随显示屏指示进行左右手握拳的运动想象;电脑实时接收采集设备输送的脑电信号,并在接收到触发信号后,对每个采样时刻的数据进行FIR高通和低通滤波后,进行独立成分分析剔除伪迹噪声;
步骤6.2:对去噪信号进行定制特征提取和分类识别操作:
(1)根据定制通道特征频率fo,对特定电极通道的信号做带通滤波,中心频率为f∈fo,带宽为1Hz,提取全部通道特征频率的数据;
(2)对(1)中提取的数据使用Multitaper算法计算PSD特征,即定制通道频率特征Z;
(3)将Z输入SVM分类器识别想象内容;
(4)将识别结果输出给反馈装置,形成闭环。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010491390.2A CN112006682B (zh) | 2020-06-02 | 2020-06-02 | 一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010491390.2A CN112006682B (zh) | 2020-06-02 | 2020-06-02 | 一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112006682A CN112006682A (zh) | 2020-12-01 |
CN112006682B true CN112006682B (zh) | 2022-09-30 |
Family
ID=73506343
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010491390.2A Active CN112006682B (zh) | 2020-06-02 | 2020-06-02 | 一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112006682B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112617851A (zh) * | 2021-01-06 | 2021-04-09 | 北京航空航天大学 | 一种基于心电信号的脑力负荷分类方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102613972A (zh) * | 2012-03-28 | 2012-08-01 | 西安电子科技大学 | 运动想象脑电信号特征的提取方法 |
CN103793058A (zh) * | 2014-02-13 | 2014-05-14 | 山西大学 | 一种主动式脑机交互系统运动想象任务分类方法及装置 |
CN108523883A (zh) * | 2018-03-19 | 2018-09-14 | 天津大学 | 一种基于实际动作建模的左右食指连续运动想象识别系统 |
CN108710895A (zh) * | 2018-04-23 | 2018-10-26 | 安徽大学 | 基于独立分量分析的运动想象脑电信号分类方法 |
CN108875799A (zh) * | 2018-05-29 | 2018-11-23 | 齐鲁工业大学 | 一种基于改进s变换的运动想象分类识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10986992B2 (en) * | 2016-08-17 | 2021-04-27 | Massachusetts Institute Of Technology | Dynamic display system and method for customizing a controller in a display system |
-
2020
- 2020-06-02 CN CN202010491390.2A patent/CN112006682B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102613972A (zh) * | 2012-03-28 | 2012-08-01 | 西安电子科技大学 | 运动想象脑电信号特征的提取方法 |
CN103793058A (zh) * | 2014-02-13 | 2014-05-14 | 山西大学 | 一种主动式脑机交互系统运动想象任务分类方法及装置 |
CN108523883A (zh) * | 2018-03-19 | 2018-09-14 | 天津大学 | 一种基于实际动作建模的左右食指连续运动想象识别系统 |
CN108710895A (zh) * | 2018-04-23 | 2018-10-26 | 安徽大学 | 基于独立分量分析的运动想象脑电信号分类方法 |
CN108875799A (zh) * | 2018-05-29 | 2018-11-23 | 齐鲁工业大学 | 一种基于改进s变换的运动想象分类识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112006682A (zh) | 2020-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113398422B (zh) | 基于运动想象-脑机接口和虚拟现实的康复训练系统及方法 | |
CN108304068B (zh) | 一种基于脑机接口的上肢康复训练机器人控制系统及方法 | |
CN109656365B (zh) | 一种基于实时闭环振动刺激增强的脑机接口方法及系统 | |
CN109992113B (zh) | 一种基于多场景诱发的mi-bci系统及其控制方法 | |
Lemm et al. | Spatio-spectral filters for improving the classification of single trial EEG | |
Citi et al. | P300-based BCI mouse with genetically-optimized analogue control | |
CN108681396A (zh) | 基于脑-肌电双模态神经信号的人机交互系统及其方法 | |
CN110442244A (zh) | 一种基于脑机接口的虚拟现实游戏交互方法及系统 | |
Ge et al. | A double-partial least-squares model for the detection of steady-state visual evoked potentials | |
CN109247917A (zh) | 一种空间听觉诱发p300脑电信号识别方法及装置 | |
CN109144238B (zh) | 一种基于眼电编码的人机交互系统及其交互方法 | |
CN106502404A (zh) | 一种基于稳态体感诱发电位的新型脑‑机接口方法及系统 | |
CN114021604A (zh) | 一种基于3d虚拟现实技术实时反馈的运动想象训练系统 | |
CN112006682B (zh) | 一种基于多通道频率特征定制的左右手运动想象脑电信号分类方法 | |
CN112987917A (zh) | 运动想象增强方法、装置、电子设备和存储介质 | |
Feng et al. | Feature extraction algorithm based on csp and wavelet packet for motor imagery eeg signals | |
Li et al. | Optimal combination of channels selection based on common spatial pattern algorithm | |
CN108509869A (zh) | 基于OpenBCI的特征集优化在线训练方法 | |
CN110688013A (zh) | 一种基于ssvep的英文键盘拼写系统及拼写方法 | |
Liu et al. | Detection of lower-limb movement intention from EEG signals | |
CN116994697A (zh) | 基于完全性脊髓损伤患者训练评估的脑机交互方法 | |
CN115474948A (zh) | 基于运动幻觉增强下肢运动想象的脑机接口方法及系统 | |
Petoku et al. | Object movement motor imagery for EEG based BCI system using convolutional neural networks | |
Molina | Direct brain-computer communication through scalp recorded EEG signals | |
CN114936574A (zh) | 一种基于bci的高灵活度机械手系统及其实现方法 |
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 |