CN113662561B - 子带级联共空间模式的脑电特征提取方法及装置 - Google Patents
子带级联共空间模式的脑电特征提取方法及装置 Download PDFInfo
- Publication number
- CN113662561B CN113662561B CN202110953882.3A CN202110953882A CN113662561B CN 113662561 B CN113662561 B CN 113662561B CN 202110953882 A CN202110953882 A CN 202110953882A CN 113662561 B CN113662561 B CN 113662561B
- Authority
- CN
- China
- Prior art keywords
- csp
- feature
- sub
- filtering
- 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
- 238000000605 extraction Methods 0.000 title claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims abstract description 102
- 238000001914 filtration Methods 0.000 claims abstract description 81
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 43
- 239000013598 vector Substances 0.000 claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 16
- 238000012549 training Methods 0.000 claims description 33
- 239000000284 extract Substances 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 description 27
- 210000004556 brain Anatomy 0.000 description 16
- 238000005070 sampling Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000007781 pre-processing Methods 0.000 description 6
- 238000000354 decomposition reaction Methods 0.000 description 5
- 210000003141 lower extremity Anatomy 0.000 description 4
- 230000002087 whitening effect Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- VZCCETWTMQHEPK-QNEBEIHSSA-N gamma-linolenic acid Chemical compound CCCCC\C=C/C\C=C/C\C=C/CCCCC(O)=O VZCCETWTMQHEPK-QNEBEIHSSA-N 0.000 description 3
- 210000003128 head Anatomy 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000007670 refining Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000004064 dysfunction Effects 0.000 description 2
- 210000000624 ear auricle Anatomy 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000005021 gait Effects 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 206010017577 Gait disturbance Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013145 classification model Methods 0.000 description 1
- 239000000084 colloidal system Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000006258 conductive agent Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000004424 eye movement Effects 0.000 description 1
- 230000004630 mental health Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 210000004761 scalp Anatomy 0.000 description 1
- 239000007787 solid Substances 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]
- A61B5/372—Analysis of electroencephalograms
-
- 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/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management
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)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Signal Processing (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
本申请涉及子带级联共空间模式的脑电特征提取方法及装置，涉及脑机接口技术领域，其中，方法包括：获取预处理后的脑电信号；将预处理后的脑电信号通过多个巴特沃夫滤波器，划分为多个频率区间不同的子带；分别将每个子带对应的脑电信号进行CSP滤波，得到每个子带对应的子特征矩阵；将各子特征矩阵进行竖向拼接，得到拼接后的特征矩阵；对拼接后的特征矩阵进行CSP滤波，得到满足预设要求的特征矩阵，并基于满足预设要求的特征矩阵，获取用于特征分类的特征向量。本申请的特征提取方法能够深度精炼特征质量和去除特征冗余性，提高基于脑电信号中的运动想象信号的特征分类的准确性，解决现有CSP算法提取特征过于冗余等问题。
Description
技术领域
本申请涉及脑机接口技术领域，尤其涉及子带级联共空间模式的脑电特征提取方法及装置。
背景技术
脑机接口技术作为一种新兴人体与外界交流的技术，通过监测和识别大脑区域的生物信号并将其转换为指令来实现对各种设备的驱动控制，已经在辅助机器人和康复工程应用中得到了初步广泛的应用。
当今全球存在大量的下肢功能障碍患者，虽然轮椅可以提供另一种行动方式，但其低效的控制方法无疑增加了使用者的负担。此外，由于其不能完全替代自然步行，自然步行能力的丧失也会对患者的心理健康产生严重影响。
而辅助行走外骨骼机器人能够有效地提供一个与外界沟通的非肌肉通道，提高了体弱老年人和步态障碍患者的生活质量。该机器人将通过捕捉受试者脑部关于下肢的运动想象信号，即想象左右腿及其肌群发生运动，但没有实际的运动输出，来驱动外骨骼机器人辅助下肢进行步行和站立。
目前，最主流且有效的脑电特征提取算法之一是共空间模式(Common SpatialPattern，CSP)算法及其变种算法，而它们在脑机接口领域的成功很大程度上取决于对脑电信号特定范围频带和特定主题特征进行正确的选择。但同时，CSP算法及其变种算法也存在着提取特征数目过多、特征过于冗余、所需计算量过大、模型收敛速度过慢等问题。
因此，有必要对现有技术予以改良以克服现有技术中的所述缺陷。
发明内容
本申请的目的在于提供一种子带级联共空间模式的脑电特征提取方法及装置，以解决现有技术中采用CSP算法对脑电信号进行特征提取时，提取特征数目过多、特征过于冗余、所需计算量过大等问题。
本申请的目的通过以下技术方案实现：
第一方面，提供一种子带级联共空间模式的脑电特征提取方法，包括：
获取运动想象任务的预处理后的脑电信号，所述预处理后的脑电信号为滤除噪声的脑电信号；
将所述预处理后的脑电信号通过多个巴特沃夫滤波器划分为多个不同频率区间的子带；
分别将每个所述子带对应的脑电信号进行CSP滤波，得到每个所述子带对应的子特征矩阵；
将各所述子特征矩阵进行竖向拼接，得到拼接后的特征矩阵；
对所述拼接后的特征矩阵进行CSP滤波，得到满足预设要求的特征矩阵，并基于所述满足预设要求的特征矩阵，获取用于特征分类的特征向量。
进一步地，在本申请第一方面一种可行的实施例中，所述对所述预处理后的脑电信号划分为多个子带，包括：
将所述预处理后的脑电信号，按照预设的频率间隔，利用巴特沃夫滤波器在频域上划分为包含预设重叠范围的多个重叠的子带
进一步地，在本申请第一方面一种可行的实施例中，所述CSP滤波，包括：
基于预先构建的子带级联CSP算法模型进行级联CSP滤波，得到满足预设要求的特征矩阵，以获取用于特征分类的特征向量。
进一步地，在本申请第一方面一种可行的实施例中，所述子带级联CSP算法模型包括至少两层CSP模型，每一层所述CSP模型的输出为下一层所述CSP模型的输入，且每一层所述CSP模型均基于对应的输入构建相应的CSP空间滤波器，所述CSP模型用于通过CSP算法对输入的脑电信号样本数据进行空间滤波。
进一步地，在本申请第一方面一种可行的实施例中，所述级联CSP滤波，包括：
将所述脑电信号样本数据，输入第一层CSP模型进行CSP滤波，得到特征矩阵；
将所述特征矩阵作为脑电信号，输入下一层所述CSP模型进行CSP滤波，得到下一层所述CSP模型输出的特征矩阵；
执行所述将所述特征矩阵作为脑电信号，输入下一层所述CSP模型进行CSP滤波，得到下一层所述CSP模型输出的特征矩阵的步骤；
直到最后一层CSP模型对输入的脑电信号进行CSP滤波，得到满足预设要求的特征矩阵。
进一步地，在本申请第一方面一种可行的实施例中，对于多分类的CSP滤波，还包括构建多分类CSP空间滤波器的步骤：
获取每一类运动想象任务对应的脑电信号的训练集样本数据，得到多类训练集样本数据；
遍历所述多类训练集样本数据，将当前训练集样本数据作为第一类数据，将其余类训练集样本数据合并为第二类数据，得到对应的多组训练集样本数据，其中，每一组训练集样本数据均由第一类数据和第二类数据组成；
对每一组训练集样本数据，分别构建对应的CSP空间滤波器，并选取所述CSP空间滤波器的前m行数据和后m行数据，组成筛选后的CSP空间滤波器；
分别选取每一组训练集样本数据对应的筛选后的CSP空间滤波器的2m行数据，组成多分类CSP空间滤波器。
第二方面，提供一种子带级联共空间模式的脑电特征提取装置，包括：
信号获取模块，用于获取运动想象任务的预处理后的脑电信号，所述预处理后的脑电信号为滤除噪声的脑电信号；
子带划分模块，用于将所述预处理后的脑电信号划分为多个子带；
CSP滤波模块，用于分别将每个所述子带对应的脑电信号进行CSP滤波，得到每个所述子带对应的子特征矩阵；
拼接模块，用于将各所述子特征矩阵进行竖向拼接，得到拼接后的特征矩阵；
级联CSP滤波模块，用于对所述拼接后的特征矩阵进行CSP滤波，得到满足预设要求的特征矩阵，并基于所述满足预设要求的特征矩阵，获取用于特征分类的特征向量。
在本申请第二方面一种可行的实施例中，所述级联滤波模块，还用于：
将所述脑电信号样本数据输入预先构建的子带级联CSP算法模型进行级联CSP滤波，得到满足预设要求的特征矩阵，以获取用于特征分类的特征向量。
第三方面，提供一种计算机设备，包括处理器和存储器，所述存储器中存储有程序，所述程序被处理器执行时，用以实现第一方面所述方法的步骤。
第四方面，提供一种计算机存储介质，所述计算机存储介质上存储有程序，所述程序在被处理器执行时，用以实现第一方面所述方法的步骤。
与现有技术相比，本申请具有如下有益效果：本申请实施例的特征提取方法，将运动想象任务的脑电信号利用巴特沃夫滤波器在频域上划分为多个子带信号，并对每个子带信号分别进行CSP滤波，得到对应的特征矩阵，然后将每个特征矩阵进行竖向拼接，得到拼接后的特征矩阵，并将拼接后的特征矩阵作为脑电信号的样本数据，进行CSP滤波，提取用于特征分类的特征向量。本申请的特征提取方法能够深度精炼特征质量和去除特征冗余性，提高基于脑电信号中的运动想象信号的特征分类的准确性，可以解决现有CSP算法提取特征数目过多、特征过于冗余、所需计算量过大等问题。
附图说明
图1是本申请一个实施例提供的实现本申请特征提取方法和装置的系统构架图；
图2是本申请一个实施例提供的特征提取方法的流程图；
图3是本申请一个实施例提供的CSP滤波算法的流程图；
图4是本申请一个实施例提供的特征提取装置的结构框图；
图5是本申请一个实施例提供的电子设备的结构框图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
针对现有技术中，通过CSP算法对脑电信号进行特征提取时，提取特征数目过多、特征过于冗余、所需计算量过大、模型收敛速度过慢等问题，本申请提供了一种基于共空间模式的脑电信号特征提取方法及装置。
图1给出了本申请一个实施例提供的实现本申请特征提取方法和装置的系统构架图，如图1所示，该系统构架包括：脑电信号采集设备101和电子设备102。其中：
脑电信号采集设备101，用于采集人脑的运动想象任务对应的原始脑电信号，并对该原始脑电信号进行预处理，该原始脑电信号为多通道的脑电信号。
示例性地，对于下肢功能障碍患者，本实施例的运动想象任务，例如可以是腿部的步行运动想象任务，例如，左腿行走、右腿行走。
脑电信号是人体重要的生理信号，是有大脑发出的电信号。本实施例的脑电信号采集设备101可以采用EMOTIV公司生产的EPOC Flex胶电极版。通过胶体导电剂作为中间介质，穿透头发，将电极和头皮连接起来，具体可以根据不同的运动想象任务，将电极放置在头部的预设位置，例如，FPz(额极中线)、C4(右中央)、FP1(左额极)和FP2(右额极)等，以采集多通道的脑电信号。
由于原始脑电信号含有不同成分伪迹干扰，例如，心电伪迹、眼动伪迹、肌电伪迹等，因此，对采集的原始脑电信号需要进行预处理，以滤除脑电信号中的伪迹信号。
电子设备102，用于获取预处理后的脑电信号，能够实现对预处理后的脑电信号进行特征提取和特征分类的功能。
本实施例的电子设备102可以为智能移动设备、计算机等，本实施例不对电子设备的设备类型作限定。
电子设备102将提取的脑电信号特征，输入训练好的分类模型进行特征分类，以确定运动类别，例如，右腿行走。然后输出对应的控制命令，驱动辅助行走设备103按照确定的运动类别辅助人体的行走动作。
本实施例中电子设备102对预处理后的脑电信号采用了深度CSP滤波的方法，实现脑电信号的特征提取。电子设备102对预处理后的脑电信号进行特征提取的具体实施方式将在下文详述。
图2给出了本申请一个实施例提供的基于共空间模式的脑电信号特征提取方法的流程图，下面以图1所示系统构架中的电子设备103为执行主体，对本申请的特征提取方法进行说明，如图2所示，所述方法包括：
S201：获取运动想象任务的预处理后的脑电信号。
具体地，通过预处理，可以滤除脑电信号中的噪声，提升脑电信号的信噪比。
本实施例可以通过EEGLAB对脑电信号进行预处理，主要包括：
FIR滤波，使用FIR滤波器过滤多余部分信号，只保留运动想象信号表征明显的4-35Hz。
基线校准，根据每一次不同运动想象任务，选取运动想象前一秒到运动想象后五秒的信号作为每一次测试的数据。并根据运动想象前一秒到运动想象开始时刻的脑电信号，基线较准该段数据。
ICA分析，通过ICA分析将脑电信号进行解构，剔除非脑电成分的伪迹信号。
关于EEGLAB对脑电信号进行预处理属于本领域比较熟知的技术，在此不再赘述。
本实施例通过电子设备102从脑电信号采集设备101中获取预处理后的脑电信号。
S202：将所述预处理后的脑电信号通过多个巴特沃夫滤波器划分为多个不同频率区间的子带。
具体地，本实施例的子带信号，指的是预处理后的脑电信号整个频率域范围内的其中一段。
可选地，本实施例将预处理后的脑电信号，按照预设的频率间隔，划分为包含预设重叠范围的多个重叠的子带信号。
其中，预设的频率间隔，例如，可以是4Hz。预设重叠范围，例如，可以是2Hz。
按照本实施例的预设频率间隔以及预设重叠范围，示例性地，将预处理后的脑电信号划分为多个子带信号，例如，可以是4-8Hz、6-10Hz、8-12Hz……等。
子带信号的划分还可以采用其他的划分方式，例如，可以将预设频率间隔设置为6Hz，将预设重叠范围设置为1Hz，具体可以是2-8Hz、9-15Hz....等，当然，也可以是没有重叠范围的，具体取决于实际的特征提取需求，本实施例对子带信号的具体划分方式不做限定。
本实施例划分的子带信号均设有重叠范围，可以在通过CSP滤波提取特征时，能够从预处理后的脑电信号中提取到足够充足的脑电信号特征。
S203：分别将每个所述子带对应的脑电信号进行CSP滤波，得到每个所述子带对应的子特征矩阵。
具体地，对每个子带信号分别输入CSP进行CSP滤波，得到每个子带信号对应的子特征矩阵，例如，4-8Hz子带信号对应的子特征矩阵，6-10Hz子带信号对应的子特征矩阵...等，从而得到与每个子带信号对应的多组子特征矩阵。
可选地，本实施例中，针对每个子带的脑电信号，均通过以下方式进行CSP滤波：
将当前子带对应的脑电信号输入预先构建的CSP空间滤波器，得到特征矩阵，具体如下：
对于二分类的运动想象任务，采用预先构建的二分类CSP空间滤波器，进行CSP滤波，得到对应的特征矩阵：
Ztest＝WsortXtest
对于多分类的运动想象任务，采用预先构建的多分类空间滤波器，进行CSP滤波，得到对应的特征矩阵：
Ztest′＝WnewXtest
其中，Wsort表示构建的二分类空间滤波器，Wnew表示构建的多分类空间滤波器，Xtest表示待提取特征的脑电信号数据的矩阵，Ztest表示二分类的运动想象任务下，待提取特征的脑电信号数据的特征矩阵，Ztest′表示多分类的运动想象任务下，待提取特征的脑电信号数据的特征矩阵。
本实施例预先构建的CSP空间滤波器，是通过对采集的脑电信号的训练集样本数据采用CSP滤波算法得到的。
可选地，图3给出了本申请一个实施例提供的CSP滤波算法的流程图，如图3所示，本实施例CSP滤波算法的步骤包括：
S301：分别计算每一类运动想象任务的脑电信号训练样本数据对应的空间协方差矩阵，以及对应的混合空间协方差矩阵。
其中，Ri表示第i类运动想象任务的脑电信号训练样本数据对应的空间协方差矩阵，Xi表示采集的脑电信号的训练集样本数据的矩阵，表示Xi的转置矩阵，/>表示Xi、/>对角线上的元素和。
从而得到混合空间协方差矩阵R为：
R＝R1+R2+...+Ri
对于二分类，i＝1,2。对于多分类，i＝1,2,3....
S302:使用奇异值分解定理对混合空间协方差矩阵R进行特征值分解，并按降序排列。
R＝UλUT
其中，U是R的特征向量矩阵，λ是对应的特征值构成的对角阵。
S303：计算白化矩阵。
通过正交白化得到白化矩阵P：
S304:构建CSP空间滤波器。
对于二分类的CSP滤波算法，i＝1,2，利用白化矩阵P分别作用于R1、R2，得到对应的两个白化后的矩阵S1、S2：
S1＝PR1PT＝Bλ1BT
S2＝PR2PT＝Bλ2BT
因为S1、S2具有公共特征向量矩阵B，而且存在两个对角矩阵λ1、λ2，通过主分量分解可得λ1+λ2等于单位矩阵I。
因此，构造空间滤波器W＝BTP，该空间滤波器W满足：当S1有最大的特征值时，S2有最小的特征值。
因为特征信息主要集中在特征矩阵的头部和尾部，而中间的特征信息不明显可以忽略，所以，选取空间滤波器W的前m行和后m行数据作为最终的二分类空间滤波器Wsort。
对于多分类的CSP滤波算法，i＝1,2,3...，针对每一类脑电信号，将该类脑电信号对应的空间协方差矩阵作为其中一类数据，并将其余类脑电信号对应的空间协方差矩阵归为另一类数据。
以R1为例,将R1归为一类数据，其余归为另一类数据，记为R′1，且R′1＝R2+R3+...。
利用白化矩阵P作用于R1、R′1，得到矩阵S1、S′1：
S1＝PR1PT＝B′λ1B′T
S′1＝PR′1PT＝B′λ1′B′T
因为S1、S′1具有公共特征向量矩阵B′，而且存在两个对角矩阵λ1、λ′1，通过主分量分解可得λ1+λ′1等于单位矩阵I。
因此，构造空间滤波器为：W′＝B′TP，该滤波器满足：当S1有最大的特征值时，S′1有最小的特征值。
因为特征信息主要集中在特征矩阵的头部和尾部，而中间的特征信息不明显可以忽略，所以选取空间滤波器W′的前m行和后m行数据得到筛选后的空间滤波器W′sort。
对于多分类的CSP滤波算法，针对每一类运动想象任务，均根据步骤S34确定筛选后的空间滤波器，并分别选取筛选后的空间滤波器W′sort的2m行数据，从而组成多分类CSP空间滤波器Wnew。
S305：CSP滤波。
构建好空间滤波器之后，对于二分类的运动想象任务，采用二分类CSP空间滤波器，对当前子带对应的训练集的脑电信号数据进行CSP滤波，得到对应的特征矩阵：
Zi＝WsortXi，i＝1,2
对于多分类的运动想象任务，采用多分类CSP空间滤波器，对当前子带对应的训练集的脑电信号数据进行CSP滤波，得到对应的特征矩阵：
Zi′＝WnewXi，i＝1,2,3....
其中，Xi表示子带对应的训练集的第i类脑电信号数据的矩阵，Zi表示二分类的运动想象任务下，子带对应的待提取特征的第i类脑电信号数据的特征矩阵，Zi′表示多分类的运动想象任务下，子带对应的第i类脑电信号数据的特征矩阵。
本实施例进一步对得到的训练集的脑电特征，进行归一化处理，得到特征向量fi：
得到的该特征向量用于训练分类器。
S204：将各所述子特征矩阵进行竖向拼接，得到拼接后的特征矩阵。
S205：对所述拼接后的特征矩阵进行CSP滤波，得到满足预设要求的特征矩阵，并基于所述满足预设要求的特征矩阵，获取用于特征分类的特征向量。
具体地，本实施例将各子带信号分别按照步骤S203进行CSP滤波后，得到对应的多组脑电信号的特征矩阵，然后将各特征矩阵全部竖向拼接在一起，对拼接后的特征矩阵再次按照步骤S203进行CSP滤波，以获得满足预设要求的特征矩阵Mtest(对应二分类运动想象任务)或M′test(对应多分类运动想象任务)，即：
Mtest＝WsortYtest
Mtest′＝WnewYtest
其中，Ytest表示拼接后的特征矩阵。
对特征矩阵Mtest或M′test，进行归一化处理，得到用于特征分类的特征向量，即：
或
其中，ftest表示二分类运动想象任务对应的特征分量，f′test表示多分类运动想象任务对应的特征分量。
本实施例能够实现更深度的精炼特征，同时去除重复冗余特征，提高了运动想象任务的脑电信号的特征分类的精度。
可选地，本实施例对拼接后的特征矩阵进行CSP滤波的过程，可以重复多次，实现级联CSP滤波，具体为：
将所述脑电信号样本数据输入预先构建的子带级联CSP算法模型，进行级联CSP滤波。
本实施例构建的子带级联CSP算法模型包括至少两层CSP模型，每一层CSP模型的输出均为下一层CSP模型的输入，本实施例基于子带级联CSP算法模型执行级联CSP滤波的步骤包括：
将所述脑电信号样本数据，输入第一层CSP模型进行CSP滤波，得到特征矩阵。将所述特征矩阵作为下一层CSP模型的输入数据，进行CSP滤波，得到下一层所述CSP模型输出的特征矩阵。执行所述将所述特征矩阵作为下一层CSP模型的输入数据，进行CSP滤波，得到下一层CSP模型输出的特征矩阵的步骤；直到最后一层CSP模型进行CSP滤波后，得到满足预设要求的特征矩阵。
本实施例的子带级联CSP算法模型中，每一层CSP模型均按照步骤S31-S35，对输入的脑电信号进行CSP滤波。
本实施例通过对拼接后的特征矩阵采用构建的子带级联CSP算法模型进行级联CSP滤波，提取特征向量后，能更进一步地增加脑电信号特征识别的准确性。
下面以一个用于对脑电信号进行二分类的特征提取的具体实例，对本申请的脑电信号特征提取方法进行说明：
步骤1，采集脑电信号样本数据。
本实施例采集的运动想象任务分别为左腿步态任务M1和右腿步态任务M2。
脑电信号样本数据是针对两名视力正常或矫正的健全受试者采集的。两名受试者在脑电信号样本数据采集开始前的几分钟内，分别练习左腿行走和右腿行走五次的真实动作，以确保运动想象任务的清晰性。
脑电采集设备是EMOTIV公司生产的EPOC Flex胶电极版。将脑电信号采集设备的6个电极遵循国际上的10-20系统，放置在FPz、FP1、FP2、Cz、C1、C2，参考电极放置在双侧耳垂A1和A2。
对两名受试者的想象的左腿行走和右腿行走进行脑电信号进行采集，以组成脑电信号样本数据。
本实施例中，取“采样频率(128Hz)*单次采样时间(1s)”，即128个采样点，作为每个脑电通道的脑电信号样本数据。然后将采集的所有脑电通道和所有试验的脑电信号汇总为三维脑电数据(即：采样点数*通道数*试验数)，从而获得两类任务标签(即左腿步行运动、右腿步行运动)的运动想象任务的脑电信号样本数据。。
步骤2，对采集的脑电信号进行预处理。
将采集的脑电信号样本数据分为训练集和测试集，并分别通过MATLAB软件的EEGLAB工具进行预处理。
将两名受试者预处理后的脑电信号样本数据，根据任务标签混在一起，得到混合数据集。
将混合数据集按7：3的比例划分训练集和测试集。例如，提取到552组数据，则取其中的400组数据作为训练集，其余的152组数据作为测试集。
步骤3，划分子带。
分别将训练集和测试集划分为多个子带，本实施例是在4-35Hz的频率范围里，以每4Hz为一个子带间隔(除最后一个为3Hz)，并加以2Hz的重叠范围，例如，子带4-8Hz、6-10Hz、8-12Hz……32-35Hz。并设置了4-35Hz的宽子带。
步骤4，CSP滤波
分别对每个子带按照步骤S31-S35进行滤波，得到每个子带对应的特征矩阵，将得到的所有特征矩阵按竖向拼接，得到拼接后的特征矩阵。
对拼接后的特征矩阵，再次按照步骤S31-S35重新进行CSP滤波，并可以执行至少两次，以提取满足预设要求的特征，实现对脑电信号特征更深度的精炼，同时去除重复冗余部分。
表1给出了本申请一个实施例的特征提取方法与CSP、MBCSP以及FBCSP等方法进行特征提取的性能比较，其中，数据里的第一项为最优值，第二项为平均值。本实施例通过分类精度和运行时间来衡量脑电信号特征分类的性能。
表1 SBCCSP和标准CSP及其它改良CSP算法的性能比较汇总
通过表1可以看出，本申请实施例的子带级联CSP算法相比于其它算法，针对运动想象任务的脑电信号的特征分类的性能，有着显著的提高，尤其是在重叠子带和宽子带的加持下，子带级联CSP算法能够从脑电信号中提取到更多有价值的脑电特征，并筛选掉冗余部分，提升了分类的能力。
另外，由于子带级联CSP算法能在提取的脑电信号特征的基础上，再次提取脑电信号特征，因此该算法在收敛速度上，也比其他算法更具有优势。
下面以一个用于多分类的特征提取的具体实例，对本申请的特征提取方法进行说明：
步骤1，采集脑电信号样本数据。
本实施例采集脑电信号样本数据的脑电采集设备同样是EMOTIV公司生产的EPOCFlex胶电极版。
以三分类为例，将脑电信号采集设备的18个电极遵循国际上的10-20系统，放置在FPz、FP1、FP2、FCz、FC1、FC2、FC3、FC4、Cz、C1、C2、C3、C4、CPz、CP1、CP2、CP3、CP4。参考电极放置在双侧耳垂A1和A2。
本实施例中，将“采样频率(128Hz)*单次采样时间(0.5s)”，即64个采样点作为每个脑电通道的脑电数据。由此获得两类任务标签(即左腿步行、右腿步行)的运动想象任务的脑电信号样本数据。
取来自提示受试者进入准备状态3s中的后2s，以每64个点作为第三个任务标签，即闲置状态(非进行运动想象任务状态)的脑电信号样本数据。
步骤2，脑电信号的预处理。
同样将采集的脑电信号样本数据分为训练集和测试集，并分别通过MATLAB软件的EEGLAB工具进行预处理。
将受试者预处理后的脑电信号样本数据，根据任务标签混在一起，得到混合数据集。
将混合数据集按7：3的比例划分训练集和测试集。例如，本实施例提取到960组数据，则取其中的672组数据作为训练集，其余的288组数据作为测试集。
步骤3，划分子带。
划分子带的方式请参见二分类实施例中的方式，在此不再赘述。
步骤4，CSP滤波
分别对每个子带按照步骤S31-S35进行滤波，得到每个子带对应的特征矩阵，将得到的所有特征矩阵按竖向拼接，得到拼接后的特征矩阵。
对拼接后的特征矩阵，再次按照步骤S31-S35重新进行CSP滤波，并可以执行至少两次，以提取满足预设要求的特征，实现对脑电信号特征更深度的精炼，同时去除重复冗余部分。
表2给出了本申请一个实施例的特征提取方法与CSP、MBCSP以及FBCSP等方法进行特征提取性能的比较结果，其中，数据里的第一项为最优值，第二项为平均值。本实施例通过分类精度和运行时间来衡量脑电信号特征分类的性能。
表2多分类SBCCSP和标准CSP及其它CSP算法的性能比较
由表2可以看出，本实施例的子带级联CSP算法相比于其它算法，针对基于运动想象任务的脑电信号的特征分类性能，有着显著的提高，尤其是在重叠子带和宽子带的加持下，子带级联CSP能够从脑电信号中提取到更多有价值的特征，并筛选掉冗余部分，提升了分类能力。
此外，由于子带级联CSP能在提取的脑电信号特征的基础上，再次提取脑电信号特征，因此该算法在收敛速度上，也比其他算法更具有优势。
本申请还提供一种双主动全桥变换器的负载稳态控制装置，图4给出了本申请一个实施例提供的双主动全桥变换器的负载稳态控制装置的框图，如图4所示，该装置包括：
信号获取模块，用于获取运动想象任务的预处理后的脑电信号，所述预处理后的脑电信号为滤除噪声的脑电信号；
子带划分模块，用于将所述预处理后的脑电信号划分为多个子带；
第一级联CSP滤波模块，用于分别将每个所述子带对应的脑电信号进行CSP滤波，得到每个所述子带对应的子特征矩阵；
拼接模块，用于将各所述子特征矩阵进行竖向拼接，得到拼接后的特征矩阵；
第二级联CSP滤波模块，用于对所述拼接后的特征矩阵进行CSP滤波，得到满足预设要求的特征矩阵，并基于所述满足预设要求的特征矩阵，获取用于特征分类的特征向量。
进一步地，级联滤波模块，还用于：
将所述脑电信号样本数据输入预先构建的子带级联CSP算法模型进行级联CSP滤波，得到满足预设要求的特征矩阵，以获取用于特征分类的特征向量。
上述实施例提供的控制装置与对应的控制方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
需要说明的是：上述实施例中提供的控制装置，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要将上述功能分配由不同的功能模块完成，即将控制装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。
本申请实施例还提供一种电子设备，图5是本申请一个实施例提供的电子设备的结构框图，如图5所示，该电子设备包括处理器和存储器，其中：
处理器，可以包括一个或多个处理核心，比如：4核心处理器、6核心处理器等。处理器可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field－Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。
存储器，可以包括高速随机存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(Secure Digital,SD)卡、闪存卡(Flash Card)、至少一个磁盘存储器件、内存器件、或其他易失性固态存储器件。
本实施例的存储器中存储有计算机程序，所述计算机程序可在所述处理器上运行，所述处理器执行所述计算机程序时，可以实现本申请特征提取方法或上述特征提取装置相关实施例中的所有或部分实施步骤，和/或文本中描述的其他内容。
本领域技术人员可以理解，图5仅仅是本申请实施例电子设备的一种可能的实现方式，其他实施方式中，还可以包括更多或更少的部件，或者组合某些部件，或者不同部件，本实施例对此不作限定。
本申请还提供一种计算机存储介质，所述计算机存储介质上存储有程序，所述程序被处理器执行时，用以实现上述特征提取方法实施例的步骤。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对申请专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请的保护范围应以所附权利要求为准。
Claims (6)
1.一种子带级联共空间模式的脑电特征提取方法，其特征在于，包括：
获取运动想象任务的预处理后的脑电信号，所述预处理后的脑电信号为滤除噪声的脑电信号；
将所述预处理后的脑电信号通过多个巴特沃夫滤波器划分为多个不同频率区间的子带；
分别将每个所述子带对应的脑电信号进行CSP滤波，得到每个所述子带对应的子特征矩阵；
将各所述子特征矩阵进行竖向拼接，得到拼接后的特征矩阵；
对所述拼接后的特征矩阵进行CSP滤波，得到满足预设要求的特征矩阵，并基于所述满足预设要求的特征矩阵，获取用于特征分类的特征向量；
所述将所述预处理后的脑电信号通过多个巴特沃夫滤波器划分为多个不同频率区间的子带，包括：
将所述预处理后的脑电信号，按照预设的频率间隔，利用巴特沃夫滤波器在频域上划分为包含预设重叠范围的多个重叠的子带；
所述对所述拼接后的特征矩阵进行CSP滤波，包括：
基于预先构建的子带级联CSP算法模型进行级联CSP滤波，得到满足预设要求的特征矩阵，以获取用于特征分类的特征向量；
所述子带级联CSP算法模型包括至少两层CSP模型，每一层所述CSP模型的输出为下一层所述CSP模型的输入，且每一层所述CSP模型均基于对应的输入构建相应的CSP空间滤波器，所述CSP模型用于通过CSP算法对输入的脑电信号样本数据进行空间滤波。
2.根据权利要求1所述的方法，其特征在于，所述级联CSP滤波，包括：
将脑电信号样本数据，输入第一层CSP模型进行CSP滤波，得到特征矩阵；
将所述特征矩阵作为脑电信号，输入下一层所述CSP模型进行CSP滤波，得到下一层所述CSP模型输出的特征矩阵；
执行所述将所述特征矩阵作为脑电信号，输入下一层所述CSP模型进行CSP滤波，得到下一层所述CSP模型输出的特征矩阵的步骤；
直到最后一层CSP模型对输入的脑电信号进行CSP滤波，得到满足预设要求的特征矩阵。
3.根据权利要求1-2任一项所述的方法，其特征在于，对于多分类的CSP滤波，还包括构建多分类CSP空间滤波器的步骤：
获取每一类运动想象任务对应的脑电信号的训练集样本数据，得到多类训练集样本数据；
遍历所述多类训练集样本数据，将当前训练集样本数据作为第一类数据，将其余类训练集样本数据合并为第二类数据，得到对应的多组训练集样本数据，其中，每一组训练集样本数据均由第一类数据和第二类数据组成；
对每一组训练集样本数据，分别构建对应的CSP空间滤波器，并选取所述CSP空间滤波器的前m行数据和后m行数据，组成筛选后的CSP空间滤波器；
分别选取每一组训练集样本数据对应的筛选后的CSP空间滤波器的2m行数据，组成多分类CSP空间滤波器。
4.一种子带级联共空间模式的脑电特征提取装置，其特征在于，包括：
信号获取模块，用于获取运动想象任务的预处理后的脑电信号，所述预处理后的脑电信号为滤除噪声的脑电信号；
子带划分模块，用于将所述预处理后的脑电信号通过多个巴特沃夫滤波器划分为多个不同频率区间的子带；所述将所述预处理后的脑电信号通过多个巴特沃夫滤波器划分为多个不同频率区间的子带，包括：
将所述预处理后的脑电信号，按照预设的频率间隔，利用巴特沃夫滤波器在频域上划分为包含预设重叠范围的多个重叠的子带；
第一级联CSP滤波模块，用于分别将每个所述子带对应的脑电信号进行CSP滤波，得到每个所述子带对应的子特征矩阵；
拼接模块，用于将各所述子特征矩阵进行竖向拼接，得到拼接后的特征矩阵；
第二级联CSP滤波模块，用于对所述拼接后的特征矩阵进行CSP滤波，得到满足预设要求的特征矩阵，并基于所述满足预设要求的特征矩阵，获取用于特征分类的特征向量；
所述对所述拼接后的特征矩阵进行CSP滤波，包括：
基于预先构建的子带级联CSP算法模型进行级联CSP滤波，得到满足预设要求的特征矩阵，以获取用于特征分类的特征向量；
所述子带级联CSP算法模型包括至少两层CSP模型，每一层所述CSP模型的输出为下一层所述CSP模型的输入，且每一层所述CSP模型均基于对应的输入构建相应的CSP空间滤波器，所述CSP模型用于通过CSP算法对输入的脑电信号样本数据进行空间滤波。
5.一种电子设备，其特征在于，包括处理器和存储器，所述存储器中存储有程序，所述程序被处理器执行时，用以实现权利要求1-3任一项所述方法的步骤。
6.一种计算机存储介质，其特征在于，所述计算机存储介质上存储有程序，所述程序在被处理器执行时，用以实现权利要求1-3任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110953882.3A CN113662561B (zh) | 2021-08-19 | 2021-08-19 | 子带级联共空间模式的脑电特征提取方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110953882.3A CN113662561B (zh) | 2021-08-19 | 2021-08-19 | 子带级联共空间模式的脑电特征提取方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113662561A CN113662561A (zh) | 2021-11-19 |
| CN113662561B true CN113662561B (zh) | 2024-03-01 |
Family
ID=78543857
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110953882.3A Active CN113662561B (zh) | 2021-08-19 | 2021-08-19 | 子带级联共空间模式的脑电特征提取方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113662561B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115208755A (zh) * | 2022-07-13 | 2022-10-18 | 北京交通大学 | 面向物联网设备资源友好的特征提取器部署方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102613972A (zh) * | 2012-03-28 | 2012-08-01 | 西安电子科技大学 | 运动想象脑电信号特征的提取方法 |
| JP2017202183A (ja) * | 2016-05-12 | 2017-11-16 | 株式会社国際電気通信基礎技術研究所 | 脳波パターン分類装置、脳波パターン分類方法、脳波パターン分類プログラムおよびニューロフィードバックシステム |
| CN109691996A (zh) * | 2019-01-02 | 2019-04-30 | 中南大学 | 一种基于混合二进制编码的脑电信号特征优选及分类器优选方法 |
| CN109858537A (zh) * | 2019-01-22 | 2019-06-07 | 南京邮电大学 | 改进的eemd与csp结合的脑电信号特征提取方法 |
| CN113255447A (zh) * | 2021-04-23 | 2021-08-13 | 桂林电子科技大学 | 脑电信号解码方法、装置、设备及计算机可读存储介质 |
-
2021
- 2021-08-19 CN CN202110953882.3A patent/CN113662561B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102613972A (zh) * | 2012-03-28 | 2012-08-01 | 西安电子科技大学 | 运动想象脑电信号特征的提取方法 |
| JP2017202183A (ja) * | 2016-05-12 | 2017-11-16 | 株式会社国際電気通信基礎技術研究所 | 脳波パターン分類装置、脳波パターン分類方法、脳波パターン分類プログラムおよびニューロフィードバックシステム |
| CN109691996A (zh) * | 2019-01-02 | 2019-04-30 | 中南大学 | 一种基于混合二进制编码的脑电信号特征优选及分类器优选方法 |
| CN109858537A (zh) * | 2019-01-22 | 2019-06-07 | 南京邮电大学 | 改进的eemd与csp结合的脑电信号特征提取方法 |
| CN113255447A (zh) * | 2021-04-23 | 2021-08-13 | 桂林电子科技大学 | 脑电信号解码方法、装置、设备及计算机可读存储介质 |
Non-Patent Citations (1)
| Title |
|---|
| 重叠特征策略与参数优化的运动想象脑电模式识别;罗天健 等;模式识别与人工智能;33(08);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113662561A (zh) | 2021-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Palaniappan | Method of identifying individuals using VEP signals and neural network | |
| CN104091172B (zh) | 一种运动想象脑电信号的特征提取方法 | |
| Kachenoura et al. | ICA: a potential tool for BCI systems | |
| CN112488002B (zh) | 基于n170的情绪识别方法及识别系统 | |
| CN109480834A (zh) | 一种基于快速多维经验模态分解的脑电信号分类方法 | |
| Zhao et al. | Deep CNN model based on serial-parallel structure optimization for four-class motor imagery EEG classification | |
| CN112270255B (zh) | 一种脑电信号识别方法、装置、电子设备及存储介质 | |
| CN109657646B (zh) | 生理时间序列的特征表示与提取方法、装置及存储介质 | |
| CN110074779A (zh) | 一种脑电信号识别方法及装置 | |
| Aljalal et al. | Feature extraction of EEG based motor imagery using CSP based on logarithmic band power, entropy and energy | |
| Lu et al. | Combined CNN and LSTM for motor imagery classification | |
| CN109858537A (zh) | 改进的eemd与csp结合的脑电信号特征提取方法 | |
| CN113662561B (zh) | 子带级联共空间模式的脑电特征提取方法及装置 | |
| Vijayendra et al. | A performance study of 14-channel and 5-channel EEG systems for real-time control of unmanned aerial vehicles (UAVs) | |
| CN117520891A (zh) | 一种运动想象脑电信号分类方法及系统 | |
| CN107016355A (zh) | 一种低误触发率运动想象双层分类识别方法 | |
| CN116919422A (zh) | 基于图卷积的多特征情感脑电识别模型建立方法及装置 | |
| WO2024040797A1 (zh) | 基于脑电数据的自闭症评估装置、方法、终端设备和介质 | |
| CN111931578A (zh) | 一种基于最小生成树和区域双层网络的脑电识别方法 | |
| CN109144277B (zh) | 一种基于机器学习实现脑控智能小车的构建方法 | |
| Velásquez-Martínez et al. | Motor imagery classification for BCI using common spatial patterns and feature relevance analysis | |
| Belwafi et al. | Online adaptive filters to classify left and right hand motor imagery | |
| CN116662782A (zh) | 一种基于MSFF-SENet的运动想象脑电图解码方法 | |
| CN115105093A (zh) | 基于功率谱密度预确定频率带的eeg信号分类识别方法 | |
| Mu et al. | A channel selection method for motor imagery EEG based on Fisher score of OVR-CSP |
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 |