CN114237383B - 基于前额单导联脑电信号的多状态识别方法 - Google Patents
基于前额单导联脑电信号的多状态识别方法 Download PDFInfo
- Publication number
- CN114237383B CN114237383B CN202111322156.8A CN202111322156A CN114237383B CN 114237383 B CN114237383 B CN 114237383B CN 202111322156 A CN202111322156 A CN 202111322156A CN 114237383 B CN114237383 B CN 114237383B
- Authority
- CN
- China
- Prior art keywords
- signals
- signal
- forehead
- electroencephalogram
- minimum value
- 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
- 210000001061 forehead Anatomy 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000004886 head movement Effects 0.000 claims abstract description 25
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 19
- 210000004556 brain Anatomy 0.000 claims abstract description 13
- 238000013145 classification model Methods 0.000 claims abstract description 10
- 230000008569 process Effects 0.000 claims abstract description 9
- 230000000875 corresponding effect Effects 0.000 claims description 38
- 238000012549 training Methods 0.000 claims description 23
- 230000004397 blinking Effects 0.000 claims description 14
- 230000009471 action Effects 0.000 claims description 13
- 239000000203 mixture Substances 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 230000001055 chewing effect Effects 0.000 claims description 7
- 238000003491 array Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 238000011897 real-time detection Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 230000002040 relaxant effect Effects 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 5
- 230000003993 interaction Effects 0.000 description 7
- 238000002570 electrooculography Methods 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 230000005611 electricity Effects 0.000 description 2
- 230000004424 eye movement Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 230000003183 myoelectrical effect Effects 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- 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
-
- 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/25—Bioelectric electrodes therefor
- A61B5/251—Means for maintaining electrode contact with the body
- A61B5/256—Wearable electrodes, e.g. having straps or bands
-
- 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/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/291—Bioelectric electrodes therefor specially adapted for particular uses for 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/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/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/384—Recording apparatus or displays specially adapted therefor
-
- 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/398—Electrooculography [EOG], e.g. detecting nystagmus; Electroretinography [ERG]
-
- 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
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- 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/013—Eye tracking input arrangements
-
- 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
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Human Computer Interaction (AREA)
- Psychology (AREA)
- General Physics & Mathematics (AREA)
- Neurosurgery (AREA)
- Neurology (AREA)
- Ophthalmology & Optometry (AREA)
- Dermatology (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种基于前额单导联脑电信号的多状态识别方法。适用于脑电领域。本发明所采用的技术方案是:一种基于前额单导联脑电信号的多状态识别方法,其特征在于:获取经滤除漂移伪迹及高频伪迹的前额叶脑电信号;求取脑电信号的极值点,当依次检测到一个极小值‑极大值‑极小值对时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个特征点;采用经训练的分类模型对特征点进行粗分类;当粗分类检测到特征点是眨眼或头动伪迹后,采用模板匹配算法精确的区分头动信号和眼电信号。本发明通过前额单导联电极同时采集脑电信号,并能从脑电信号中精确识别出眼电信号,提高了检测过程中设备佩戴的舒适性。
Description
技术领域
本发明涉及一种基于前额单导联脑电信号的多状态识别方法。适用于脑电领域。
背景技术
脑电信号具有较高的时间分辨率和精度,能够准确地实时反映人的生理状态,已经在疲劳检测、状态识别、人机交互等方面得到了广泛推广。脑电信号是由大脑神经元放电产生的,极为微弱,在采集过程中容易受到各种噪声的干扰,如低频噪声、高频噪声、工频噪声、眼动伪迹、肌电伪迹等。在大多数的脑电研究中,此类伪迹相对于目标脑电信号而言是无用的,并且为此开发了不少的处理算法而予以降低伪迹,也收获比较好的效果。然而,伪迹本身也是人行为的一种表征,比如检测眨眼伪迹可以研究人的眨眼模式,通过眨眼频率等参数就可以建立起眨眼和疲劳等相关关系;检测低频高幅伪迹可以刻画人的动作行为等。
已有研究中,如[Ma,J.,et al.,A novel EOG/EEG hybrid human–machineinterface adopting eye movements and ERPs:Application to robot control.]中,作者结合眼电与脑电开发了一套控制系统获得了比较好的效果,但是由于采用了垂直眼电,眼电电极需要放置在眼睛的正下方,这影响了一定的美观度和佩戴舒适度。
在文献[姚娟娟,路堃,马兴宇,等.基于脑电和眼电的疲劳检测方法的研究[J].电子设计工程,2020,v.28;No.428(06):121-126]中,作者采用前额单通道脑电信号以及脑电信号中的眼电信号特征来检测疲劳状态,取得93.10%的疲劳识别准确率和3.79%的误报率。这说明合理利用相关伪迹信号可以更多元地获取人的生理信号,从而评估人的行为,更好得解决人机交互问题。
现实与生理相关信号的处理中,随着人的状态变化及传感器的氧化等问题,不同时间段采集到的信号大多会有差异。尤其是对于非平稳的脑电等信号,受人的状态和所处的环境等影响较大。
发明内容
本发明要解决的技术问题是:针对上述存在的问题,提供一种基于前额单导联脑电信号的多状态识别方法。
本发明所采用的技术方案是:一种基于前额单导联脑电信号的多状态识别方法,其特征在于:
获取经滤除漂移伪迹及高频伪迹的前额叶脑电信号;
求取脑电信号的极值点,当依次检测到一个极小值-极大值-极小值对时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个特征点;
采用经训练的分类模型对特征点进行粗分类;
当粗分类检测到特征点是眨眼或头动伪迹后,采用模板匹配算法精确的区分头动信号和眼电信号;
所述分类模型的训练,包括:
获取用户分别完成放松、咀嚼、眨眼和头动动作任务时的前额叶脑电信号,并对脑电信号进行滤除漂移伪迹及高频伪迹的处理;
求取信号极值点,当依次检测到一个极小值-极大值-极小值对时且该极大值在相应动作任务进行过程中时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个与该相应动作任务对应的特征点;
基于各动作任务对应的特征点训练所述分类模型。
所述分类模型采用混合高斯模型,基于各动作任务对应的特征点使用EM算法估计所述混合高斯模型的参数。
所述采用模板匹配算法精确的区分头动信号和眼电信号,包括:
选取放松任务对应特征点中的幅值组成数组A,对A中幅值从小到大进行排序,选取其中数组A长度95%的位置点的数值为阈值Thr1;
对属于眨眼任务的特征点进行筛选,去掉其中幅值小于Thr1的特征点;
基于筛选后的眨眼任务特征点生成眼电模板,如下:
1)选取特征点对应脑电信号段的前0.15秒至特征点后0.15秒为可能的眨眼信号段;
2)将第一个特征点形成的信号段为类似眼电模板,并记此模型的计数值为1;
当第二个特征点形成的信号段与类似眼电模板的相似度大于预先设定的相似度阈值Thr_sim,则认为这两段信号为同一类信号,此时对模板进行更新,同时计数值加1;
当相似度不大于Thr_sim时,则认为这两段信号分属不同的模板,此时把这段信号当做一个新的眼电模板;
3)对所有特征点对应脑电信号段按步骤2)对特征点依次进行处理后再对得到的模板进行更新,直到任意两个模板之间的相似度小于Thr_sim则停止模板的更新;
4)将步骤3)得到的计数值最大的眼电模板作为最终的眼电模板。
所述相似度的计算包括:
假设存在两同长度的一维数组,A和B;数据的下标从1开始,特征点的下标为i;
a.上下平移使得将两数组的特征点位置处重叠,假设对A数组进行平移,则平移后的数组A’有:
A’=A-(A(i)–B(i))
b.计算Sshare,Sshare为A’与B两曲线围成的面积;
c.计算Stotal,Stotal为A’与B两曲线围成的面积以及两曲线处于同侧时与横坐标围成的公共区域面积之和;
d.则相似度的计算为sim=(Stotal-Sshare)/Stotal。
所述模板进行更新,包括:
假设源模板数组A,模板的计数为idx,目标模板数组B;
For ii=1:length(A)
A(ii)=[A(ii)*idx+B(ii)]/(idx+1)
Idx=idx+1。
所述混合高斯模型的训练发生在该识别方法执行开始时、用户主动选择训练时、按期定时训练时。
一种基于前额单导联脑电信号的多状态识别装置,其特征在于:包括模型训练模块和实时检测模块;
其中,实时检测模块包括:
信号获取模块,用于获取经滤除漂移伪迹及高频伪迹的前额叶脑电信号;
特征提取模块,用于求取脑电信号的极值点,当依次检测到一个极小值-极大值-极小值对时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个特征点;
粗分类模块,用于采用经训练的混合高斯模型对特征点进行粗分类;
精分类模块,用于当粗分类检测到特征点是眨眼或头动伪迹后,采用模板匹配算法精确的区分头动信号和眼电信号;
所述模型训练模块,用于获取用户分别完成放松、咀嚼、眨眼和头动动作任务时的前额叶脑电信号,并对脑电信号进行滤除漂移伪迹及高频伪迹的处理;求取信号极值点,当依次检测到一个极小值-极大值-极小值对时且该极大值在相应动作任务进行过程中时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个与该相应动作任务对应的特征点;基于各动作任务对应的特征点使用EM算法估计所述混合高斯模型的参数。
一种存储介质,其上存储有能被处理器执行的计算机程序,其特征在于:所述计算机程序被执行时实现所述基于前额单导联脑电信号的多状态识别方法的步骤。
一种基于前额单导联脑电信号的多状态识别系统,其特征在于,包括:
信号采集模块,用于采集前额叶的脑电信号,采集到信号首先经过4阶的50Hz的陷波器,之后经过10阶的0.1~30Hz的带通滤波器滤除漂移伪迹以及高频的伪迹;
信号处理模块,具有存储器和处理器,存储器上存储有能被处理器执行的计算机程序,其特征在于:所述计算机程序被执行时实现权利要求1~5任意一项所述基于前额单导联脑电信号的多状态识别方法的步骤;
人机交互模块,用于脑电信号和信号处理模块从脑电信号中识别出的眼电信号进行交互控制。
本发明的有益效果是:本发明通过前额单导联电极同时采集脑电信号,并能从脑电信号中精确识别出眼电信号,提高了检测过程中设备佩戴的舒适性。本发明将模板匹配算法结合到模型训练中,提高系统识别运算的鲁棒性和准确性。本发明提供的识别系统能够在每次启动时进行模型自学习,实时更新模型,减少环境或人为因素对脑电信号采集的影响,提高了特征识别的准确性。本发明包括多种状态的识别,提高系统的适应性和可扩展性,在人机交互领域具有广泛的应用潜力。
附图说明
图1为实施例的流程图。
具体实施方式
本实施例为一种基于前额单导联脑电信号的多状态识别方法,具体包括以下步骤:
S1、获取信号采集模块采集的前额叶的脑电信号,信号采集模块采集到信号首先经过4阶的50Hz的陷波器,之后经过10阶的0.1~30Hz的带通滤波器滤除漂移伪迹以及高频的伪迹;
S2、求取脑电信号的极值点,当依次检测到一个极小值-极大值-极小值对时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个特征点;
S3、采用经训练的混合高斯模型对特征点进行粗分类;
S4、当步骤S3粗分类检测到特征点是眨眼或头动伪迹后,采用模板匹配算法精确的区分头动信号和眼电信号。
本实施例中混合高斯模型的训练发生在该识别方法执行开始时、用户主动选择训练时和按期定时训练时,混合高斯模型的训练,包括:
A、获取用户分别完成放松、咀嚼、眨眼和头动动作任务时的前额叶脑电信号,并对脑电信号进行滤除漂移伪迹及高频伪迹的处理。
本实施例中首先给出15秒的准备时间,提示用户按照需要进入模型训练中,然后依次呈现给用户放松、咀嚼、眨眼、头动等四种指令,每种指令持续15秒,每俩指令间给出8秒的放松时间,在每种指令出现时,要求用户按照指令不断重复指令的动作。当用户完成4种指令的动作任务后,自动进行四种动作任务特征的混合高斯模型训练。
B、求取信号极值点,当依次检测到一个极小值-极大值-极小值对时且该极大值在相应动作任务进行过程中时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个与该相应动作任务对应的特征点;
C、基于各动作任务对应的特征点使用EM算法估计所述混合高斯模型的参数,分别计算出四种特征的期望值及方差值。
本实施例中采用模板匹配算法精确的区分头动信号和眼电信号,包括:
选取放松任务对应特征点中的幅值组成数组A,对A中幅值从小到大进行排序,选取其中数组A长度95%的位置点的数值为阈值Thr1;
在实际训练中,当系统要求训练者进行放松任务时,训练者可能由于无意的出现眨眼或者其他大幅动作,这些实际是会影响后面分类效果。这里选择的95%的位置点是在放松线索下排除幅度最大的5%特征值,从而减少因大幅干扰造成的影响,该数值是根据实际操作中一个比较理想的比例值;
对属于眨眼任务的特征点进行筛选,去掉其中幅值小于Thr1的特征点;
基于筛选后的眨眼任务特征点生成眼电模板,如下:
1)选取特征点对应脑电信号段的前0.15秒至特征点后0.15秒为可能的眨眼信号段(选取该特征点前0.15秒至特征点后0.15秒的数据段作为疑似眨眼信号段,之所以选择0.15秒是因为普通人眨一次眼的时间约0.3秒);
2)将第一个特征点形成的信号段为类似眼电模板,并记此模型的计数值为1;
当第二个特征点形成的信号段与类似眼电模板的相似度大于设定的相似度阈值Thr_sim,则认为这两段信号为同一类信号,此时对模板进行更新,同时计数值加1;
当相似度不大于Thr_sim时,则认为这两段信号分属不同的模板,此时把这段信号当做一个新的眼电模板;
3)对所有特征点对应脑电信号段按步骤2)对特征点依次进行处理后再对得到的模板进行更新,直到任意两个模板之间的相似度小于Thr_sim则停止模板的更新;
4)将步骤3)得到的计数值最大的眼电模板作为最终的眼电模板。
本实施例中相似度的计算包括:
假设存在两同长度的一维数组,A和B;数据的下标从1开始,特征点的下标为i;
a、上下平移使得将两数组的特征点位置处重叠,假设对A数组进行平移,则平移后的数组A’有:
A’=A-(A(i)–B(i))
b、计算Sshare,Sshare为A’与B两曲线围成的面积;
c、计算Stotal,Stotal为A’与B两曲线围成的面积以及两曲线处于同侧时与横坐标围成的公共区域面积之和;
d、则相似度的计算为sim=(Stotal-Sshare)/Stotal。
本实施例中模板进行更新的方法,包括:
假设源模板数组A,模板的计数为idx,目标模板数组B;
For ii=1:length(A)
A(ii)=[A(ii)*idx+B(ii)]/(idx+1)
Idx=idx+1。
本实施例还提供一种基于前额单导联脑电信号的多状态识别装置,主要包括模型训练模块和实时检测模块,其中实时检测模块包括信号获取模块、特征提取模块、粗分类模块和精分类模块。
本例中信号获取模块用于获取经滤除漂移伪迹及高频伪迹的前额叶脑电信号;特征提取模块用于求取脑电信号的极值点,当依次检测到一个极小值-极大值-极小值对时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个特征点;粗分类模块用于采用经训练的混合高斯模型对特征点进行粗分类;精分类模块,用于当粗分类检测到特征点是眨眼或头动伪迹后,采用模板匹配算法精确的区分头动信号和眼电信号。
本例中模型训练模块用于获取用户分别完成放松、咀嚼、眨眼和头动动作任务时的前额叶脑电信号,并对脑电信号进行滤除漂移伪迹及高频伪迹的处理;求取信号极值点,当依次检测到一个极小值-极大值-极小值对时且该极大值在相应动作任务进行过程中时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个与该相应动作任务对应的特征点;基于各动作任务对应的特征点使用EM算法估计混合高斯模型的参数。
本实施例还提供一种存储介质,其上存储有能被处理器执行的计算机程序,该计算机程序被执行时实现本例中基于前额单导联脑电信号的多状态识别方法的步骤。
本实施例还提供一种基于前额单导联脑电信号的多状态识别系统,包括信号采集模块、信号处理模块和人机交互模块。
信号采集模块主要为一个脑电采集设备(头带或者脑电帽等),主要采集前额叶的脑电信号[FP1/FP2]。采样率为250,采集到信号首先经过4阶的50Hz的陷波器,之后经过10阶的0.1~30Hz的带通滤波器滤除漂移伪迹以及高频的伪迹。
信号处理模块具有存储器和处理器,存储器上存储有能被处理器执行的计算机程序,该计算机程序被执行时实现本例中基于前额单导联脑电信号的多状态识别方法的步骤。
人机交互模块在每次启动时根据使用场景进行模型的自学习,实时更新模型,进行四类特征的识别,本例应用这四类特征进行简单unity3d的模型识别。
Claims (8)
1.一种基于前额单导联脑电信号的多状态识别方法,其特征在于:
获取经滤除漂移伪迹及高频伪迹的前额叶脑电信号;
求取脑电信号的极值点,当依次检测到一个极小值-极大值-极小值对时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个特征点;
采用经训练的分类模型对特征点进行粗分类;
当粗分类检测到特征点是眨眼或头动伪迹后,采用模板匹配算法精确的区分头动信号和眼电信号;
所述分类模型的训练,包括:
获取用户分别完成放松、咀嚼、眨眼和头动动作任务时的前额叶脑电信号,并对脑电信号进行滤除漂移伪迹及高频伪迹的处理;
求取信号极值点,当依次检测到一个极小值-极大值-极小值对时且该极大值在相应动作任务进行过程中时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个与该相应动作任务对应的特征点;
基于各动作任务对应的特征点训练所述分类模型;
所述采用模板匹配算法精确的区分头动信号和眼电信号,包括:
选取放松任务对应特征点中的幅值组成数组A,对A中幅值从小到大进行排序,选取其中数组A长度95%的位置点的数值为阈值Thr1;
对属于眨眼任务的特征点进行筛选,去掉其中幅值小于Thr1的特征点;
基于筛选后的眨眼任务特征点生成眼电模板,如下:
1)选取特征点对应脑电信号段的前0.15秒至特征点后0.15秒为可能的眨眼信号段;
2)将第一个特征点形成的信号段为类似眼电模板,并记此模型的计数值为1;
当第二个特征点形成的信号段与类似眼电模板的相似度大于预先设定的相似度阈值Thr_sim,则认为这两段信号为同一类信号,此时对模板进行更新,同时计数值加1;
当相似度不大于Thr_sim时,则认为这两段信号分属不同的模板,此时把这段信号当做一个新的眼电模板;
3)对所有特征点对应脑电信号段按步骤2)对特征点依次进行处理后再对得到的模板进行更新,直到任意两个模板之间的相似度小于Thr_sim则停止模板的更新;
4)将步骤3)得到的计数值最大的眼电模板作为最终的眼电模板。
2.根据权利要求1所述的基于前额单导联脑电信号的多状态识别方法,其特征在于:所述分类模型采用混合高斯模型,基于各动作任务对应的特征点使用EM算法估计所述混合高斯模型的参数。
3.根据权利要求1所述的基于前额单导联脑电信号的多状态识别方法,其特征在于,所述相似度的计算包括:
假设存在两同长度的一维数组,A和B;数据的下标从1开始,特征点的下标为i;
a.上下平移使得将两数组的特征点位置处重叠,假设对A数组进行平移,则平移后的数组A’有:
A’=A-(A(i)–B(i))
b.计算Sshare,Sshare为A’与B两曲线围成的面积;
c.计算Stotal,Stotal为A’与B两曲线围成的面积以及两曲线处于同侧时与横坐标围成的公共区域面积之和;
d.则相似度的计算为sim=(Stotal-Sshare)/Stotal。
4.根据权利要求1所述的基于前额单导联脑电信号的多状态识别方法,其特征在于,所述模板进行更新,包括:
假设源模板数组A,模板的计数为idx,目标模板数组B;
For ii=1:length(A)
A(ii)=[A(ii)*idx+B(ii)]/(idx+1)
Idx=idx+1。
5.根据权利要求1所述的基于前额单导联脑电信号的多状态识别方法,其特征在于,所述分类模型的训练发生在该识别方法执行开始时、用户主动选择训练时、按期定时训练时。
6.一种基于前额单导联脑电信号的多状态识别装置,其特征在于:包括模型训练模块和实时检测模块;
其中,实时检测模块包括:
信号获取模块,用于获取经滤除漂移伪迹及高频伪迹的前额叶脑电信号;
特征提取模块,用于求取脑电信号的极值点,当依次检测到一个极小值-极大值-极小值对时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个特征点;
粗分类模块,用于采用经训练的混合高斯模型对特征点进行粗分类;
精分类模块,用于当粗分类检测到特征点是眨眼或头动伪迹后,采用模板匹配算法精确的区分头动信号和眼电信号;
所述模型训练模块,用于获取用户分别完成放松、咀嚼、眨眼和头动动作任务时的前额叶脑电信号,并对脑电信号进行滤除漂移伪迹及高频伪迹的处理;求取信号极值点,当依次检测到一个极小值-极大值-极小值对时且该极大值在相应动作任务进行过程中时记录该极大值的位置、幅值以及前后极小值之间的宽度,形成一个与该相应动作任务对应的特征点;基于各动作任务对应的特征点使用EM算法估计所述混合高斯模型的参数;
所述采用模板匹配算法精确的区分头动信号和眼电信号,包括:
选取放松任务对应特征点中的幅值组成数组A,对A中幅值从小到大进行排序,选取其中数组A长度95%的位置点的数值为阈值Thr1;
对属于眨眼任务的特征点进行筛选,去掉其中幅值小于Thr1的特征点;
基于筛选后的眨眼任务特征点生成眼电模板,如下:
1)选取特征点对应脑电信号段的前0.15秒至特征点后0.15秒为可能的眨眼信号段;
2)将第一个特征点形成的信号段为类似眼电模板,并记此模型的计数值为1;
当第二个特征点形成的信号段与类似眼电模板的相似度大于预先设定的相似度阈值Thr_sim,则认为这两段信号为同一类信号,此时对模板进行更新,同时计数值加1;
当相似度不大于Thr_sim时,则认为这两段信号分属不同的模板,此时把这段信号当做一个新的眼电模板;
5)对所有特征点对应脑电信号段按步骤2)对特征点依次进行处理后再对得到的模板进行更新,直到任意两个模板之间的相似度小于Thr_sim则停止模板的更新;
6)将步骤3)得到的计数值最大的眼电模板作为最终的眼电模板。
7.一种存储介质,其上存储有能被处理器执行的计算机程序,其特征在于:所述计算机程序被执行时实现权利要求1~5任意一项所述基于前额单导联脑电信号的多状态识别方法的步骤。
8.一种基于前额单导联脑电信号的多状态识别系统,其特征在于,包括:
信号采集模块,用于采集前额叶的脑电信号,采集到信号首先经过4阶的50Hz的陷波器,之后经过10阶的0.1~30Hz的带通滤波器滤除漂移伪迹以及高频的伪迹;
信号处理模块,具有存储器和处理器,存储器上存储有能被处理器执行的计算机程序,其特征在于:所述计算机程序被执行时实现权利要求1~5任意一项所述基于前额单导联脑电信号的多状态识别方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111322156.8A CN114237383B (zh) | 2021-11-09 | 2021-11-09 | 基于前额单导联脑电信号的多状态识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111322156.8A CN114237383B (zh) | 2021-11-09 | 2021-11-09 | 基于前额单导联脑电信号的多状态识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114237383A CN114237383A (zh) | 2022-03-25 |
CN114237383B true CN114237383B (zh) | 2024-03-12 |
Family
ID=80748943
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111322156.8A Active CN114237383B (zh) | 2021-11-09 | 2021-11-09 | 基于前额单导联脑电信号的多状态识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114237383B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115969398B (zh) * | 2022-12-02 | 2024-02-02 | 深湾创新技术(深圳)有限公司 | 眨眼检测方法及装置 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101599127A (zh) * | 2009-06-26 | 2009-12-09 | 安徽大学 | 眼电信号的特征提取与识别方法 |
CN102274032A (zh) * | 2011-05-10 | 2011-12-14 | 北京师范大学 | 一种基于脑电信号的驾驶员疲劳检测系统 |
WO2012151498A2 (en) * | 2011-05-04 | 2012-11-08 | Cardioinsight Technologies, Inc. | Signal averaging |
CN103340637A (zh) * | 2013-06-06 | 2013-10-09 | 同济大学 | 基于眼动与脑电融合的驾驶员警觉度智能监控系统及方法 |
WO2016201731A1 (zh) * | 2015-06-17 | 2016-12-22 | 宇龙计算机通信科技(深圳)有限公司 | 一种指纹识别方法、装置及电子设备 |
CN107510453A (zh) * | 2017-10-12 | 2017-12-26 | 北京翼石科技有限公司 | 一种前额区脑电分析方法 |
CN108420430A (zh) * | 2018-04-02 | 2018-08-21 | 东北电力大学 | 一种基于嗅觉脑电波和pso-svm的感官物质分类方法 |
CN110269609A (zh) * | 2019-07-08 | 2019-09-24 | 江苏理工学院 | 基于单通道的脑电信号中眼电伪迹分离方法 |
CN110737331A (zh) * | 2019-09-11 | 2020-01-31 | 浙江迈联医疗科技有限公司 | 基于多任务脑机接口的个性化认知训练及康复方法、装置及设备 |
CN111046779A (zh) * | 2019-12-06 | 2020-04-21 | 江苏理工学院 | 基于单通道脑电信号的有意眨眼识别方法及装置 |
CN111956217A (zh) * | 2020-07-15 | 2020-11-20 | 山东师范大学 | 一种面向实时脑电信号的眨眼伪迹识别方法及系统 |
-
2021
- 2021-11-09 CN CN202111322156.8A patent/CN114237383B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101599127A (zh) * | 2009-06-26 | 2009-12-09 | 安徽大学 | 眼电信号的特征提取与识别方法 |
WO2012151498A2 (en) * | 2011-05-04 | 2012-11-08 | Cardioinsight Technologies, Inc. | Signal averaging |
CN102274032A (zh) * | 2011-05-10 | 2011-12-14 | 北京师范大学 | 一种基于脑电信号的驾驶员疲劳检测系统 |
CN103340637A (zh) * | 2013-06-06 | 2013-10-09 | 同济大学 | 基于眼动与脑电融合的驾驶员警觉度智能监控系统及方法 |
WO2016201731A1 (zh) * | 2015-06-17 | 2016-12-22 | 宇龙计算机通信科技(深圳)有限公司 | 一种指纹识别方法、装置及电子设备 |
CN107510453A (zh) * | 2017-10-12 | 2017-12-26 | 北京翼石科技有限公司 | 一种前额区脑电分析方法 |
CN108420430A (zh) * | 2018-04-02 | 2018-08-21 | 东北电力大学 | 一种基于嗅觉脑电波和pso-svm的感官物质分类方法 |
CN110269609A (zh) * | 2019-07-08 | 2019-09-24 | 江苏理工学院 | 基于单通道的脑电信号中眼电伪迹分离方法 |
CN110737331A (zh) * | 2019-09-11 | 2020-01-31 | 浙江迈联医疗科技有限公司 | 基于多任务脑机接口的个性化认知训练及康复方法、装置及设备 |
CN111046779A (zh) * | 2019-12-06 | 2020-04-21 | 江苏理工学院 | 基于单通道脑电信号的有意眨眼识别方法及装置 |
CN111956217A (zh) * | 2020-07-15 | 2020-11-20 | 山东师范大学 | 一种面向实时脑电信号的眨眼伪迹识别方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN114237383A (zh) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106108894A (zh) | 一种提高情绪识别模型时间鲁棒性的情绪脑电识别方法 | |
Ortiz-Rosario et al. | Wavelet methodology to improve single unit isolation in primary motor cortex cells | |
CN102722727A (zh) | 基于脑功能网络邻接矩阵分解的脑电特征提取方法 | |
EP1934677A1 (en) | Method and system for detecting and classifying facial muscle movements | |
Pun et al. | Brain-computer interaction research at the Computer Vision and Multimedia Laboratory, University of Geneva | |
CN113274037B (zh) | 一种动态脑功能网络的生成方法、系统及设备 | |
Ramos-Aguilar et al. | Analysis of EEG signal processing techniques based on spectrograms | |
Fatima et al. | Towards a low cost Brain-computer Interface for real time control of a 2 DOF robotic arm | |
CN114237383B (zh) | 基于前额单导联脑电信号的多状态识别方法 | |
CN117918863A (zh) | 一种脑电信号实时伪迹处理与特征提取的方法和系统 | |
CN111144450B (zh) | 一种构建基于不同长度名字刺激的erp范式的方法 | |
CN109009098A (zh) | 一种运动想象状态下的脑电信号特征识别方法 | |
CN108491792B (zh) | 基于眼电信号的办公场景人机交互行为识别方法 | |
CN111736690A (zh) | 基于贝叶斯网络结构辨识的运动想象脑机接口 | |
CN110705656A (zh) | 一种基于eeg传感器的面部动作识别方法 | |
Noce et al. | Spike detection: the first step towards an ENG-based neuroprosheses | |
Vahdani-Manaf | Classification of eye movement signals using electrooculography in order to device controlling | |
Zacharaki et al. | One-class classification of temporal EEG patterns for K-complex extraction | |
Heger et al. | Filling a glass of water: continuously decoding the speed of 3D hand movements from EEG signals | |
KR20140009635A (ko) | Pso를 이용한 erd/ers특징 추출 | |
Wei et al. | Automatic recognition of epileptic discharges based on shape similarity in time-domain | |
Vařeka et al. | Off-line analysis of the P300 event-related potential using discrete wavelet transform | |
Kim et al. | One-Dimensional Shallow Neural Network Using Non-Fiducial Based Segmented Electrocardiogram for User Identification System | |
Krupinski et al. | Real–time low–latency estimation of the blinking and EOG signals | |
CN113017647A (zh) | 监测用户清醒状态的方法、装置、电子设备及介质 |
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 |