CN106991409B - 一种运动想象脑电信号特征提取与分类***及方法 - Google Patents
一种运动想象脑电信号特征提取与分类***及方法 Download PDFInfo
- Publication number
- CN106991409B CN106991409B CN201710244075.8A CN201710244075A CN106991409B CN 106991409 B CN106991409 B CN 106991409B CN 201710244075 A CN201710244075 A CN 201710244075A CN 106991409 B CN106991409 B CN 106991409B
- Authority
- CN
- China
- Prior art keywords
- motor imagery
- training
- feature
- matrix
- electroencephalogram signal
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000000605 extraction Methods 0.000 title claims abstract description 46
- 238000012549 training Methods 0.000 claims abstract description 83
- 238000012360 testing method Methods 0.000 claims abstract description 80
- 239000013598 vector Substances 0.000 claims abstract description 65
- 238000004364 calculation method Methods 0.000 claims abstract description 35
- 230000009467 reduction Effects 0.000 claims abstract description 31
- 238000006243 chemical reaction Methods 0.000 claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims description 59
- 238000004070 electrodeposition Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 238000012935 Averaging Methods 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 230000000284 resting effect Effects 0.000 claims description 4
- 210000004556 brain Anatomy 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 210000003710 cerebral cortex Anatomy 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 208000012661 Dyskinesia Diseases 0.000 description 1
- 230000007177 brain activity Effects 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000003920 cognitive function Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000007659 motor function Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 210000001428 peripheral nervous system Anatomy 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 230000037152 sensory function Effects 0.000 description 1
Images
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/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/01—Input arrangements or combined input and output arrangements for interaction between user and computer
- G06F3/011—Arrangements for interaction with the human body, e.g. for user immersion in virtual reality
- G06F3/015—Input arrangements based on nervous system activity detection, e.g. brain waves [EEG] detection, electromyograms [EMG] detection, electrodermal response detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Neurology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Neurosurgery (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Dermatology (AREA)
- Human Computer Interaction (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biomedical Technology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Signal Processing (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种运动想象脑电信号特征提取与分类***及方法,该方法基于一种运动想象脑电信号特征提取与分类***,该方法的具体步骤包括:对运动想象脑电信号的训练:将脑机接口***采集的运动想象脑电信号作为训练数据进行格式转换、特征点提取、特征向量计算、降维、字典学习和类别直方图计算,得到训练字典及训练数据的类别直方图;对运动想象脑电信号的测试:将脑机接口***采集的运动想象脑电信号作为测试数据进行格式转换、特征点提取、特征向量计算、降维,利用训练模块中的训练字典得到测试数据的稀疏表达,得到测试数据的类别直方图,并根据训练数据的类别直方图与测试数据的类别直方图的对比结果进行测试数据分类。
Description
技术领域
本发明属于脑电信号处理的技术领域,尤其涉及一种运动想象脑电信号特征提取与分类***及方法。
背景技术
目前,脑机接口技术作为一种全新的通信和控制技术,可以为思维正常但有严重运动障碍的患者提供语言交流和环境控制手段。其中,脑机接口被定义为能够使人不依赖周围神经***和肌肉、而与外界进行通信或控制的设备。脑机接口技术不仅应用于未患者提供语言交流和环境控制,在自动控制、军事认知等科学领域也有潜在的应用价值。鉴于其巨大的应用前景,脑机接口已经引起国际科学界的高度重视,称为脑科学、康复工程、生物医学工程及人机控制领域的一个研究热点。
在所有能够被监测到的反映大脑活动的信号中,由于脑电信号(Electroencephalogram,EEG)具有较好的时间分辨率,监测仪器较简单,非入侵等优点,被大部分脑机接口***所采纳。
在所有基于EEG信号的脑机接口***中,能够被用作控制信号的大脑活动信号有:视觉诱发电位、P300事件相关电位、运动想象、皮层慢电位等几种。其中,针对运动想象信号,基于运动想象的脑机接口的理论基础是:当人进行某中意识活动时,大脑皮层的对应区域处于活动状态,这些区域产生的μ节律和β波的低频部分将出现幅度衰减,即事件相关去同步(event-related desynchronization,ERD);在某一时刻,大脑皮层的局部区域没有受到意识模式的激励,该区域内的脑电信号EEG局部成分将出现幅度的增强,即事件相关同步(event-related synchronization,ERS)。基于ERD/ERS的脑机接口***主要是辨别运动想象思维作业的EEG信号,比如想象左手、右手、脚、舌的运动,从而产生不同的控制命令。
目前,基于运动想象脑机接口***的研究是目前脑机接口的研究热点之一。主要原因是ERS/ERD现象的生理依据和数学模型已被许多科研机构研究和证明,并成为目前探讨正常和病理状态下识别感觉、运动和认知功能最为常用的方法。另外,该现象中的μ节律和β波变化不需要外界刺激诱发,便于受试者训练和控制,进而成为脑机接口研究最受欢迎的实现途径。
但是,目前基于运动想象的脑机接口在特征提取与分类方面存在如下问题:
(1)分类种类少:基于不同运动想象思维作业的特征提取与分类十分困难,目前能够区分的最多的运动想象思维作业是六种。分类种类少限制了脑机接口的应用范围;
(2)分类精度不高:随着运动想象种类的增加,分类精度随之下降;
(3)无法随时区分工作与休息状态:脑机接口若要用户长期佩戴,就需要持续工作,则相应脑机接口***需要能够区分用户的运动想象状态和休息状态,用户不可能长时间处于运动想象的控制状态。目前的脑机接口***还不具备该功能。
综上所述,对于基于运动想象的脑机接口在特征提取与分类方面分类种类少、分类精度低和无法随时区分工作与休息状态的问题,尚缺乏有效的解决方案。
发明内容
本发明为了解决上述问题,克服现有技术中基于运动想象的脑机接口在特征提取与分类方面存在的分类种类少、分类精度低和无法随时区分工作与休息状态的问题,提供一种运动想象脑电信号特征提取与分类***及方法。
为了实现上述目的,本发明采用如下技术方案:
一种运动想象脑电信号特征提取与分类***,该***包括训练模块和测试模块;
所述训练模块被配置为将脑机接口***采集的运动想象脑电信号作为训练数据进行格式转换、特征点提取、特征向量计算、降维、字典学习和类别直方图计算,得到训练字典及训练数据的类别直方图的模块;
所述测试模块被配置为将脑机接口***采集的运动想象脑电信号作为测试数据进行格式转换、特征点提取、特征向量计算、降维,利用训练模块中的训练字典得到测试数据的稀疏表达,得到测试数据的类别直方图,并根据训练数据的类别直方图与测试数据的类别直方图的对比结果进行测试数据分类的模块。
进一步的,所述训练模块包括第一脑电信号存储模块,所述第一脑电信号存储模块将脑机接口***采集的运动想象脑电信号作为训练数据存储;
所述测试模块包括第二脑电信号存储模块,所述第二脑电信号存储模块将脑机接口***采集的运动想象脑电信号作为测试数据存储。
本发明为了解决上述问题,克服现有技术中基于运动想象的脑机接口在特征提取与分类方面存在的分类种类少、分类精度低和无法随时区分工作与休息状态的问题,提供一种运动想象脑电信号特征提取与分类***及方法。
为了实现上述目的,本发明采用如下技术方案:
一种运动想象脑电信号特征提取与分类方法,该方法基于一种运动想象脑电信号特征提取与分类***,该方法的具体步骤包括:
(1)对运动想象脑电信号的训练:将脑机接口***采集的运动想象脑电信号作为训练数据进行格式转换、特征点提取、特征向量计算、降维、字典学习和类别直方图计算,得到训练字典及训练数据的类别直方图;
(2)对运动想象脑电信号的测试:将脑机接口***采集的运动想象脑电信号作为测试数据进行格式转换、特征点提取、特征向量计算、降维,利用训练模块中的训练字典得到测试数据的稀疏表达,得到测试数据的类别直方图,并根据训练数据的类别直方图与测试数据的类别直方图的对比结果进行测试数据分类。
进一步的,所述步骤(1)中对训练数据进行格式转换的具体步骤为:
将采集的运动想象脑电信号由二维信号转换为三维信号:将第一脑电信号存储模块中存储的脑机接口***采集的运动想象脑电信号X∈IRn×m×T构建二维空间矩阵。
在对运动想象脑电信号的分析训练中,建立能够描述脑电信号空间信息的二维矩阵,将采集的运动想象脑电信号由二维信号转换为三维信号,使得脑电信号描述过程中,不仅包含时间信息,也能够包含空间信息。
进一步的,所述步骤(1)中对特征点提取和特征向量计算的具体步骤为:
(1-1)关键点检测:第一脑电信号存储模块中存储的脑机接口***采集的运动想象脑电信号X∈IRn×m×T表示为X=[X1,X2,...,XT];
对于X中的每个值利用Harris关键点检测器检测关键点,针对检测到的所有关键点提取大小为(η×η×T)的模块;
(1-2)特征值提取:针对步骤(1-1)中的模块所有值去均值,针对模块里面的每一个像素沿着时间方向计算其方差M2,偏态M3和峰度M4,组成与模块对应的时间矩阵Mr,r={2,3,4},其中Mr=[mij],i,j=1,2,...,η
其中,vijt是位置在{i,j}的第t个时间片;
(1-3)特征向量计算:将步骤(1-2)中的每一个时间矩阵转化为一个长度为η2的向量,时间矩阵Mr,r={2,3,4}中r的三个值组成特征向量m∈IRd,其中d=3η2为关键点模块的特征向量长度,
进一步的,所述步骤(1)中降维采用随机投影(random projection)方法进行降维处理,具体步骤为:
记描述该运动想象特征向量m的矩阵为D∈IRd×p,其中,d为关键点模块的特征向量长度,p为关键点模块的数量;
将维数为d的向量投影到一个维数为n的子空间中,其中n<<d;
通过矩阵D乘以一个随机矩阵R来实现特征向量矩阵的降维处理,将特征向量矩阵维数降低为n×p:
Y=RD
其中,Y∈IRn×p为降维后的特征向量矩阵,R为随机投影矩阵,R∈IRn×d,均值为0,方差为1。
进一步的,所述步骤(1)中的字典学习根据降维后的特征向量矩阵Y,通过压缩传感得到Y的稀疏表达,具体步骤为:
其中,||·||F是Frobenius范数,||·||0是l0半范数,计算向量里包含的非零元素。
进一步的,所述步骤(2)中测试数据的稀疏表达计算的具体步骤为:
将步骤(2)中测试数据分别经过格式转换、特征点提取、特征向量计算、降维得到的降维后的特征向量矩阵Q,结合步骤(1)中得到的字典Φ,
根据:
得到测试数据的稀疏表达XQ。
进一步的,所述步骤(1)中的训练数据的类别直方图和所述步骤(2)中的测试数据的类别直方图分别根据:
得到,所述步骤(1)中的训练数据的类别直方图为hi,所述步骤(2)中的训练数据的类别直方图为hQ;
所述步骤(2)中测试数据分类的具体步骤为:
根据:
确定所述测试数据的所述类别。
本发明的有益效果:
1.分类种类得到了提高:本发明的一种运动想象脑电信号特征提取与分类***及方法,在步骤(1)的对运动想象脑电信号的训练中,不同的思维作业可以建立一个共同的字典Φ,在步骤(2)的对运动想象脑电信号的测试中,只要在训练部分出现的思维作业,在测试过程中就可以检测出该类别,因此,此种算法的分类种类得到了提高;
2.运动想象与休息状态均可进行检测:本发明的一种运动想象脑电信号特征提取与分类***及方法,将休息状态当做一种思维作业,在建立字典Φ时考虑进去该状态,则在后续的测试中,当出现休息状态时,也可以检测出用户目前处于休息状态;
3.分类速度较快:本发明的一种运动想象脑电信号特征提取与分类***及方法,不同的思维作业建立一个共同的字典,在测试中不同的思维作业可以一次性得到分类结果。
附图说明
图1是本发明的方法流程图;
图2是本发明的国际10-20电极位置示意图。
具体实施方式:
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面结合附图与实施例对本发明作进一步说明。
实施例1:
正如背景技术所介绍的,现有技术中基于运动想象的脑机接口在特征提取与分类方面存在的分类种类少、分类精度低和无法随时区分工作与休息状态的问题,提供一种运动想象脑电信号特征提取与分类***及方法。
本申请的一种运动想象脑电信号特征提取与分类***的实施方式中,
为了实现上述目的,本发明采用如下技术方案:
一种运动想象脑电信号特征提取与分类***,该***包括训练模块和测试模块;
所述训练模块被配置为将脑机接口***采集的运动想象脑电信号作为训练数据进行格式转换、特征点提取、特征向量计算、降维、字典学习和类别直方图计算,得到训练字典及训练数据的类别直方图的模块;
所述测试模块被配置为将脑机接口***采集的运动想象脑电信号作为测试数据进行格式转换、特征点提取、特征向量计算、降维,利用训练模块中的训练字典得到测试数据的稀疏表达,得到测试数据的类别直方图,并根据训练数据的类别直方图与测试数据的类别直方图的对比结果进行测试数据分类的模块。
所述训练模块包括第一脑电信号存储模块,所述第一脑电信号存储模块将脑机接口***采集的运动想象脑电信号作为训练数据存储;
所述测试模块包括第二脑电信号存储模块,所述第二脑电信号存储模块将脑机接口***采集的运动想象脑电信号作为测试数据存储。
为了实现上述目的,本发明采用如下技术方案:
一种运动想象脑电信号特征提取与分类方法,该方法基于一种运动想象脑电信号特征提取与分类***,如图1所示,该方法的具体步骤包括:
(1)对运动想象脑电信号的训练:将脑机接口***采集的运动想象脑电信号作为训练数据进行格式转换、特征点提取、特征向量计算、降维、字典学习和类别直方图计算,得到训练字典及训练数据的类别直方图;
所述步骤(1)中脑机接口***采集的运动想象脑电信号作为训练数据的具体步骤为:
用户佩戴EEG电极帽进行不同心理作业的运动想象,例如想象左手运动、右手运动,采集相应的脑电信号为X=[XL:XR]并存储于第一脑电信号存储模块中。
所述步骤(1)中对训练数据进行格式转换的具体步骤为:
将采集的运动想象脑电信号由二维信号转换为三维信号:将第一脑电信号存储模块中存储的脑机接口***采集的运动想象脑电信号X∈IRn×m×T构建二维空间矩阵。
为了体现EEG信号的空间信息,采集到的EEG信号存储为X∈IRn×m×T,其中M=n×m是为了体现EEG空间信息而构建的二维空间矩阵,其依据是国际10-20电极位置图,如图2所示。若是标准的10-20电极位置图,则构建的二维空间矩阵M为一个5*5的二维矩阵,相应位置上的值为该电极上在不同时间点测量到的脑电信号电压值,为0的位置则表示该点没有测量信号,该点的值始终为0。T是采集的时间点,若采集时间为1秒,采集频率为128,则T=128。
经过此种数据转换,测量到的脑电信号则为的三维信号,也可以看作一个的图像随着时间的变化而变化的录像。
在对运动想象脑电信号的分析训练中,建立能够描述脑电信号空间信息的二维矩阵,将采集的运动想象脑电信号由二维信号转换为三维信号,使得脑电信号描述过程中,不仅包含时间信息,也能够包含空间信息。
当大脑进行不同的运动想象时,不同电极位置上的脑电信号会随着时间的推移进行变化。同一种运动想象任务导致的这种变化是相似的,特征提取的任务就是针对某一种运动想象任务找到能够描述这种相似性的东西。
所述步骤(1)中对特征点提取和特征向量计算的具体步骤为:
(1-1)关键点检测:第一脑电信号存储模块中存储的脑机接口***采集的运动想象脑电信号X∈IRn×m×T表示为X=[X1,X2,...,XT];
针对X1,利用Harris关键点检测器检测到X1的ρ个关键点。特征提取的任务就是观察这ρ个关键点在后面的(T-1)个时间点上是如何变化的。针对这ρ个关键点,提取大小为(η×η×T)的模块。
(1-2)特征值提取:针对步骤(1-1)中的模块所有值去均值,针对模块里面的每一个像素沿着时间方向计算其方差M2,偏态M3和峰度M4,组成与模块对应的时间矩阵Mr,r={2,3,4},其中Mr=[mij],i,j=1,2,...,η
其中,vijt是位置在{i,j}的第t个时间片;
(1-3)特征向量计算:将步骤(1-2)中的每一个时间矩阵转化为一个长度为η2的向量,时间矩阵Mr,r={2,3,4}中r的三个值组成特征向量m∈IRd,其中d=3η2为关键点模块的特征向量长度,
在本实施例中,设其中一个模块为v∈IRη×η×T,首先对该模块所有值去均值,然后针对模块里面的每一个像素沿着时间方向计算其方差M2,偏态M3和峰度M4。
假设与模块v相关的时间矩阵为Mr,r={2,3,4},其中Mr=[mij],i,j=1,2,...,η
其中vijt是位置在{i,j}的第t个时间片。每一个时间矩阵都可以转化为一个长度为η2的向量,,因此r的三个值可以组成向量m∈IRd,其中d=3η2
上述方法得到的特征向量一般都是高维数的,对后续的计算十分不利。
在本实施例中,在不损失信号特征的情况下,所述步骤(1)中降维采用随机投影(random projection)方法进行降维处理,具体步骤为:
训练数据脑电信号得到了p个关键点模块,每一个关键点模块的特征向量长度为d,
记描述该运动想象特征向量m的矩阵为D∈IRd×p,其中,d为关键点模块的特征向量长度,p为关键点模块的数量;
将维数为d的向量投影到一个维数为n的子空间中,其中n<<d;
通过矩阵D乘以一个随机矩阵R来实现特征向量矩阵的降维处理,将特征向量矩阵维数降低为n×p:
Y=RD
其中,Y∈IRn×p为降维后的特征向量矩阵,R为随机投影矩阵,R∈IRn×d,均值为0,方差为1。
所述步骤(1)中的字典学习根据降维后的特征向量矩阵Y,通过压缩传感得到Y的稀疏表达,具体步骤为:
其中,||·||F是Frobenius范数,||·||0是l0半范数,计算向量里包含的非零元素。
K-SVD是一种经典的字典训练算法,依据误差最小原则,对误差项进行SVD分解,选择使误差最小的分解项作为更新的字典原子和对应的原子系数,经过不断的迭代从而得到优化的解。
在本实施例中对压缩感知技术进行应用,将降维后的特征向量矩阵Y压缩之后,得到稀疏表达的字典,后续的脑电信号利用该字典建立自己的稀疏表达,并利用:若后续的脑电信号属于训练字典时存在的某一种思维作业,则在该思维作业部分的稀疏表达系数应该包含最少的非零值,并建立稀疏表达直方图进行思维作业的分类。不同种类的脑电信号建立一个共同的字典,一旦字典建立,后续步骤(2)的分类过程将节省大量时间。
若在字典训练过程中,认为用户休息的状态为一种思维作业,则在后续分类过程中,也能随时把用户休息的状态分辨出来,因此解决了目前大多数脑机接口***无法解决的问题,实现了用户长期使用脑机接口,脑机接口***可以识别用户的工作状态和休息状态。
所述步骤(1)中的训练数据类别直方图hi=[h1,h2,...,hK]根据下式计算得到:
在本实施例中,步骤(1)训练部分结束,得到两个结果:1)学习字典Φ;2)类别直方图hi=[h1,h2,...,hK]。
(2)对运动想象脑电信号的测试:将脑机接口***采集的运动想象脑电信号作为测试数据进行格式转换、特征点提取、特征向量计算、降维,步骤(2)中测试数据的格式转换、特征点提取、特征向量计算与步骤(1)采用的方法公式完全一致,这里就不赘述了。
利用训练模块中的训练字典得到测试数据的稀疏表达,得到测试数据的类别直方图。
所述步骤(2)中测试数据的稀疏表达计算的具体步骤为:
将步骤(2)中测试数据分别经过格式转换、特征点提取、特征向量计算、降维得到的降维后的特征向量矩阵Q,结合步骤(1)中得到的字典Φ,
根据:
得到测试数据的稀疏表达XQ。
所述步骤(2)中的测试数据的类别直方图分别根据:
得到,所述步骤(2)中的测试数据的类别直方图为hQ;
根据训练数据的类别直方图与测试数据的类别直方图的对比结果进行测试数据分类。
所述步骤(2)中测试数据分类的具体步骤为:
根据:
确定所述测试数据的所述类别。
本发明的有益效果:
1.分类种类得到了提高:本发明的一种运动想象脑电信号特征提取与分类***及方法,在步骤(1)的对运动想象脑电信号的训练中,不同的思维作业可以建立一个共同的字典Φ,在步骤(2)的对运动想象脑电信号的测试中,只要在训练部分出现的思维作业,在测试过程中就可以检测出该类别,因此,此种算法的分类种类得到了提高;
2.运动想象与休息状态均可进行检测:本发明的一种运动想象脑电信号特征提取与分类***及方法,将休息状态当做一种思维作业,在建立字典Φ时考虑进去该状态,则在后续的测试中,当出现休息状态时,也可以检测出用户目前处于休息状态;
3.分类速度较快:本发明的一种运动想象脑电信号特征提取与分类***及方法,不同的思维作业建立一个共同的字典,在测试中不同的思维作业可以一次性得到分类结果。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (8)
1.一种运动想象脑电信号特征提取与分类***,该***包括训练模块和测试模块;其特征是:
所述训练模块被配置为将脑机接口***采集的运动想象脑电信号作为训练数据进行格式转换、特征点提取、特征向量计算、降维、字典学习和类别直方图计算,得到训练字典及训练数据的类别直方图的模块,所述训练模块中不同种类的脑电信号建立一个共同的字典,并且认为用户休息的状态为一种思维作业;
所述对训练数据进行格式转换的具体步骤为:
将采集的运动想象脑电信号由二维信号转换为三维信号:将第一脑电信号存储模块中存储的脑机接口***采集的运动想象脑电信号E∈IRn×m×T构建二维空间矩阵;
为了体现EEG信号的空间信息,采集到的EEG信号存储为E∈IRn×m×T,其中M=n×m是为了体现EEG空间信息而构建的二维空间矩阵,其依据是国际10-20电极位置图;若是标准的10-20电极位置图,则构建的二维空间矩阵M为一个5*5的二维矩阵,相应位置上的值为该电极上在不同时间点测量到的脑电信号电压值,为0的位置则表示该点没有测量信号,该点的值始终为0;T是采集的时间点,若采集时间为1秒,采集频率为128,则T=128;
所述训练模块中对特征点提取和特征向量计算依赖于如下方式:
关键点检测:第一脑电信号存储模块中存储的脑机接口***采集的运动想象脑电信号E∈IRn×m×T表示为E=[E1,E2,...,ET];针对E1,利用Harris关键点检测到E1的ρ个关键点,针对检测到的所有关键点提取大小为(η×η×T)的模块,设为N∈IRη×η×T;
特征值矩阵提取:针对所提取到的模块所有值去均值,针对去除均值后的V∈IRη×η×T的每一个像素沿着时间方向计算其特征矩阵,选择三个特征,M2-方差,M3-斜度,M4-峰度,三个特征矩阵的计算公式为:
其中,M2∈IRη×η,M3∈IRη×η,M4∈IRη×η;
特征向量计算:将三个特征矩阵分别按列转化为长度为η2的向量并构建成一个长度为d=3η2的列向量M,其中
所述测试模块被配置为将脑机接口***采集的运动想象脑电信号作为测试数据进行格式转换、特征点提取、特征向量计算、降维,利用训练模块中的训练字典得到测试数据的稀疏表达,得到测试数据的类别直方图,并根据训练数据的类别直方图与测试数据的类别直方图的对比结果进行测试数据分类的模块。
2.如权利要求1所述的一种运动想象脑电信号特征提取与分类***,其特征是:所述训练模块包括第一脑电信号存储模块,所述第一脑电信号存储模块将脑机接口***采集的运动想象脑电信号作为训练数据存储;
所述测试模块包括第二脑电信号存储模块,所述第二脑电信号存储模块将脑机接口***采集的运动想象脑电信号作为测试数据存储。
3.一种运动想象脑电信号特征提取与分类方法,该方法基于一种运动想象脑电信号特征提取与分类***,其特征是:该方法的具体步骤包括:
(1)对运动想象脑电信号的训练:将脑机接口***采集的运动想象脑电信号作为训练数据进行格式转换、特征点提取、特征向量计算、降维、字典学习和类别直方图计算,得到训练字典及训练数据的类别直方图,训练过程中不同种类的脑电信号建立一个共同的字典,并且认为用户休息的状态为一种思维作业;
(2)对运动想象脑电信号的测试:将脑机接口***采集的运动想象脑电信号作为测试数据进行格式转换、特征点提取、特征向量计算、降维,利用训练模块中的训练字典得到测试数据的稀疏表达,得到测试数据的类别直方图,并根据训练数据的类别直方图与测试数据的类别直方图的对比结果进行测试数据分类;
所述对训练数据进行格式转换的具体步骤为:
将采集的运动想象脑电信号由二维信号转换为三维信号:将第一脑电信号存储模块中存储的脑机接口***采集的运动想象脑电信号E∈IRn×m×T构建二维空间矩阵;
为了体现EEG信号的空间信息,采集到的EEG信号存储为E∈IRn×m×T,其中M=n×m是为了体现EEG空间信息而构建的二维空间矩阵,其依据是国际10-20电极位置图;若是标准的10-20电极位置图,则构建的二维空间矩阵M为一个5*5的二维矩阵,相应位置上的值为该电极上在不同时间点测量到的脑电信号电压值,为0的位置则表示该点没有测量信号,该点的值始终为0;T是采集的时间点,若采集时间为1秒,采集频率为128,则T=128;
所述步骤(1)中对特征点提取和特征向量计算的具体步骤为:
(1-1)关键点检测:第一脑电信号存储模块中存储的脑机接口***采集的运动想象脑电信号E∈IRn×m×T表示为E=[E1,E2,...,ET];对于E中的每个值利用Harris关键点检测器检测关键点,针对检测到的所有关键点提取大小为(η×η×T)的模块,设为N∈IRη×η×T;
(1-2)特征值矩阵提取:针对所提取到的模块所有值去均值,针对去除均值后的V∈IRη×η×T的每一个像素沿着时间方向计算其特征矩阵,选择三个特征,M2-方差,M3-斜度,M4-峰度,三个特征矩阵的计算公式为:
其中,M2∈IRη×η,M3∈IRη×η,M4∈IRη×η;
(1-3)特征向量计算:将三个特征矩阵分别按列转化为长度为η2的向量并构建成一个长度为d=3η2的列向量M,其中
4.如权利要求3所述的一种运动想象脑电信号特征提取与分类方法,其特征是:所述步骤(1)中降维采用随机投影(random projection)方法进行降维处理,具体步骤为:
记描述该运动想象特征向量m的矩阵为D∈IRd×p,其中,d为关键点模块的特征向量长度,p为关键点模块的数量;
将维数为d的向量投影到一个维数为z的子空间中,其中z<<d;
通过矩阵D乘以一个随机矩阵R来实现特征向量矩阵的降维处理,将特征向量矩阵维数降低为z×p:
Y=RD
其中,Y∈IRz×p为降维后的特征向量矩阵,R为随机投影矩阵,R∈IRz×d,均值为0,方差为1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710244075.8A CN106991409B (zh) | 2017-04-14 | 2017-04-14 | 一种运动想象脑电信号特征提取与分类***及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710244075.8A CN106991409B (zh) | 2017-04-14 | 2017-04-14 | 一种运动想象脑电信号特征提取与分类***及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106991409A CN106991409A (zh) | 2017-07-28 |
CN106991409B true CN106991409B (zh) | 2020-03-27 |
Family
ID=59416328
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710244075.8A Expired - Fee Related CN106991409B (zh) | 2017-04-14 | 2017-04-14 | 一种运动想象脑电信号特征提取与分类***及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106991409B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109784211A (zh) * | 2018-12-26 | 2019-05-21 | 西安交通大学 | 一种基于深度学习的运动想象脑电信号分类方法 |
CN109558863B (zh) * | 2019-01-07 | 2020-11-06 | 哈尔滨工业大学(深圳) | 脑电信号特征表示与提取方法、装置及存储介质 |
CN111695500A (zh) * | 2020-06-11 | 2020-09-22 | 齐鲁工业大学 | 基于迁移学习的脑卒中患者运动想象任务识别方法及*** |
CN112075932B (zh) * | 2020-10-15 | 2023-12-05 | 中国医学科学院生物医学工程研究所 | 一种诱发电位信号的高分辨率时频分析方法 |
CN114167983B (zh) * | 2021-11-24 | 2024-07-23 | 哈尔滨工业大学(深圳) | 一种基于共幅相测量的脑机接口*** |
CN115079012B (zh) * | 2022-07-22 | 2022-11-29 | 羿动新能源科技有限公司 | 大数据预测电池soc的特征数据提取***及提取方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101558997A (zh) * | 2009-05-27 | 2009-10-21 | 江西蓝天学院 | 基于二阶盲辨识的脑电信号识别方法 |
CN105468138A (zh) * | 2015-07-15 | 2016-04-06 | 武汉理工大学 | 基于脑机接口技术与激光雷达的智能车辆避障导航方法 |
CN105956624A (zh) * | 2016-05-06 | 2016-09-21 | 东南大学 | 基于空时频优化特征稀疏表示的运动想象脑电分类方法 |
CN106388814A (zh) * | 2016-10-11 | 2017-02-15 | 天津大学 | 基于最优核时频分布可视图的癫痫脑电信号识别方法 |
CN106446762A (zh) * | 2016-07-07 | 2017-02-22 | 哈尔滨工程大学 | 一种基于稀疏表征分类算法的运动想象脑电分类处理方法 |
-
2017
- 2017-04-14 CN CN201710244075.8A patent/CN106991409B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101558997A (zh) * | 2009-05-27 | 2009-10-21 | 江西蓝天学院 | 基于二阶盲辨识的脑电信号识别方法 |
CN105468138A (zh) * | 2015-07-15 | 2016-04-06 | 武汉理工大学 | 基于脑机接口技术与激光雷达的智能车辆避障导航方法 |
CN105956624A (zh) * | 2016-05-06 | 2016-09-21 | 东南大学 | 基于空时频优化特征稀疏表示的运动想象脑电分类方法 |
CN106446762A (zh) * | 2016-07-07 | 2017-02-22 | 哈尔滨工程大学 | 一种基于稀疏表征分类算法的运动想象脑电分类处理方法 |
CN106388814A (zh) * | 2016-10-11 | 2017-02-15 | 天津大学 | 基于最优核时频分布可视图的癫痫脑电信号识别方法 |
Non-Patent Citations (5)
Title |
---|
"一种联合阴影和目标区域图像的SAR目标识别方法";丁军等;《电子与信息学报》;20150331;第37卷(第3期);参见第4.2节 * |
"在线脑机接口中脑电信号的特征提取与分类方法";徐宝国等;《电子学报》;20110531(第5期);参见第2-4节 * |
"基于ICM的局部不变特征提取方法";江友谊等;《红外技术》;20120331;第34卷(第3期);参见摘要 * |
"基于内容的双字典学习及稀疏表示的图像重构";王小玉等;《仪器仪表学报》;20130831;第34卷(第8期);参见第2-3节 * |
"基于排列组合熵的脑电意识任务识别方法的研究";范影乐等;《传感技术学报》;20080131;第21卷(第1期);第74-78页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106991409A (zh) | 2017-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106991409B (zh) | 一种运动想象脑电信号特征提取与分类***及方法 | |
Zhang et al. | Sparse Bayesian classification of EEG for brain–computer interface | |
Lakshmi et al. | Survey on EEG signal processing methods | |
Pei et al. | A tensor-based frequency features combination method for brain–computer interfaces | |
Zheng et al. | Transfer components between subjects for EEG-based emotion recognition | |
Ilyas et al. | A survey of analysis and classification of EEG signals for brain-computer interfaces | |
Liyanage et al. | Dynamically weighted ensemble classification for non-stationary EEG processing | |
KR101293446B1 (ko) | 움직임을 상상하는 뇌파 분류 방법 및 그 장치 | |
CN107169462B (zh) | 一种基于层次分析的脑电信号特征分类的二分类方法 | |
Ameri et al. | Projective dictionary pair learning for EEG signal classification in brain computer interface applications | |
Dodia et al. | Lie detection using extreme learning machine: A concealed information test based on short‐time Fourier transform and binary bat optimization using a novel fitness function | |
Mousa et al. | A novel brain computer interface based on principle component analysis | |
Pfister et al. | Robustifying independent component analysis by adjusting for group-wise stationary noise | |
Jinliang et al. | EEG emotion recognition based on granger causality and capsnet neural network | |
Higashi et al. | Multilinear discriminant analysis with subspace constraints for single-trial classification of event-related potentials | |
Korats et al. | Applying ICA in EEG: Choice of the Window Length and of the Decorrelation Method | |
CN111820876A (zh) | 一种脑电空间滤波器的动态构建方法 | |
Wang et al. | A gradient-based automatic optimization CNN framework for EEG state recognition | |
Miah et al. | Movie oriented positive negative emotion classification from eeg signal using wavelet transformation and machine learning approaches | |
Gao et al. | EEG classification based on sparse representation and deep learning | |
Wang et al. | Comparison of different classification methods for EEG-based brain computer interfaces: a case study | |
Anil Kumar et al. | Performance analysis of blind source separation using canonical correlation | |
CN107423668B (zh) | 基于小波变换和稀疏表达的脑电信号分类***与方法 | |
Su et al. | Electroencephalogram classification in motor-imagery brain–computer interface applications based on double-constraint nonnegative matrix factorization | |
Cong | Blind source separation |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200616 Address after: 710000 25 / F, block D, Tsinghua Science Park, Keji 2nd Road, Zhangba Street office, hi tech Zone, Xi'an City, Shaanxi Province Patentee after: Xi'an Huinao Intelligent Technology Co.,Ltd. Address before: Feng Ming Licheng District Lingang Development Zone, 250101 Ji'nan Road, Shandong province No. 1000 Patentee before: SHANDONG JIANZHU University |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200327 |
|
CF01 | Termination of patent right due to non-payment of annual fee |