CN114145752B - 一种基于小波变换的多模态脑机接口数据融合方法 - Google Patents
一种基于小波变换的多模态脑机接口数据融合方法 Download PDFInfo
- Publication number
- CN114145752B CN114145752B CN202111234957.9A CN202111234957A CN114145752B CN 114145752 B CN114145752 B CN 114145752B CN 202111234957 A CN202111234957 A CN 202111234957A CN 114145752 B CN114145752 B CN 114145752B
- Authority
- CN
- China
- Prior art keywords
- signal
- wavelet
- signals
- matrix
- fusion
- 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
- 230000009466 transformation Effects 0.000 title claims abstract description 13
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 11
- 230000004927 fusion Effects 0.000 claims abstract description 50
- 238000000034 method Methods 0.000 claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims description 52
- 238000002474 experimental method Methods 0.000 claims description 16
- 238000000354 decomposition reaction Methods 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 11
- 230000000004 hemodynamic effect Effects 0.000 claims description 10
- 210000004556 brain Anatomy 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 6
- 230000002087 whitening effect Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 5
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 claims description 4
- 108010064719 Oxyhemoglobins Proteins 0.000 claims description 4
- 108010002255 deoxyhemoglobin Proteins 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000002329 infrared spectrum Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 7
- 108010054147 Hemoglobins Proteins 0.000 description 4
- 102000001554 Hemoglobins Human genes 0.000 description 4
- 230000003340 mental effect Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000002610 neuroimaging Methods 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000004497 NIR spectroscopy Methods 0.000 description 1
- 208000012902 Nervous system disease Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001054 cortical effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000002582 magnetoencephalography Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 210000004761 scalp Anatomy 0.000 description 1
- 230000002123 temporal effect Effects 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/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14542—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
- A61B5/14551—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
- A61B5/14553—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases specially adapted for cerebral tissue
-
- 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/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
-
- 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
- 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
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/251—Fusion techniques of input or preprocessed data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Optics & Photonics (AREA)
- Psychology (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Neurology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Mathematical Physics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Fuzzy Systems (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于小波变换的多模态脑机接口数据融合方法,本发明首先对EEG信号、fNIRS信号进行小波分解,再将两组小波系数通过基于Fisher值的融合规则生成新的小波系数,然后通过小波重构得到融合信号,最后提取融合信号的共空间模式(CSP)特征,利用线性判别分析(LDA)进行分类。本发明能够有效提升BCI系统的计算效率,一定程度上提高了系统的分类精度,为脑电与功能性近红外光谱多模态数据融合提供了新思路。
Description
技术领域
本发明属于模式识别领域,涉及一种数据融合方法,特别是基于小波变换的多模态数据融合方法。
背景技术
脑机接口(Brain-computer interface,BCI)提供了大脑和外部设备之间的连接。它被广泛应用于机器人控制、精神压力检测、神经系统疾病检测等领域。BCI技术通常分为侵入式BCI和非侵入式BCI。常见的非侵入式的BCI信号采集成像方式包括脑电图(Electroencephalogram,EEG),皮层区脑电图(Electrocardiogram,ECoG),功能性近红外光谱(Functional near-infrared spectroscopy,fNIRS),功能磁共振成像(Functionalmagnetic resonance imaging,fMRI)和脑磁图(Magnetoencephalography,MEG)等。
每种脑成像模态信号都有自身的优势与局限性,比如,EEG是一种非侵入式脑成像技术,它利用头皮电极来测量神经元大量电活动引起的电压波动。这项技术具有很好的时间分辨率。但是,脑电信号极易受到运动伪影的影响,并且由于容积传导效应,导致EEG信号的空间分辨率较低。这些将会降低脑机接口的准确性。fNIRS是一种非侵入式光学成像技术,通常利用两个不同的波长(600~1000nm)来测量氧合血红蛋白(OxygenatedHemoglobin,HbO)和脱氧血红蛋白(Deoxygenated Hemoglobin,HbR)的浓度变化,并以此来测量大脑神经元的代谢活动。基于fNIRS的BCI的主要局限性在于血流动力学响应延迟较长,这导致了较差的时间分辨率,被认为是基于实时fNIRS的BCI应用的主要障碍。另一方面,fNIRS具有较好的空间分辨率,血流动力学信号不像电信号那样在通道之间传播,而且对运动伪影具有显著的鲁棒性。
因此,研究人员试图结合不同模态的优势来克服单模态的缺陷,开展多模态脑机接口系统的研究。为了将多种模态的信号进行融合,研究者们提出了许多的融合方法,比如利用联合互信息(JMI)准则对EEG特征和fNIRS特征进行组合优化的方法;通过不同的决策路径从特定方面分析大脑信息,最后通过独立决策路径融合(IDPF)的方法;通过张量映射将EEG和fNIRS时间序列的组合转换为三阶张量实现融合的方法。尽管通过这些方法融合的混合模态BCI在准确性和稳定性方面都优于单一模态,但由于不同模态之间存在数据的采样率,幅值尺度等性质不同以及多源信号的时空配准等问题,多模态的数据级融合是一个具有挑战性的问题。
发明内容
本发明针对现有技术的不足,提出了一种基于小波变换的多模态脑机接口数据融合方法。本发明通过对两种或多种模态数据进行小波分解,然后计算小波系数的Fisher值,按照Fishe值赋予小波系数不同的权重,再对小波系数加权求和得到一组新的小波系数,最后通过小波重构从而实现不同模态数据的融合。本发明包括以下步骤:
步骤(1):多通道多模态信号的采集与预处理,对两种模态中每个通道信号的预处理是进行适带通滤波;
1):频带处理,对EEG信号进行0.5-35Hz的滤波,目的是为了消除干扰和噪声;
2):将光强度测量值转换为血液动力学变化,使用六阶零相位巴特沃斯滤波器对fNIRS信号进行0.01-0.1Hz的滤波操作,同时采用实验前N秒的HbR和HbO浓度数据的均值进行基线校正;
3):插值与放缩:由于EEG和fNIRS信号的采样率不同,且两者的幅度差异较大,为了保证小波变换的兼容性和融合的效果,需要对fNIRS信号进行插值,使其与EEG的采样率保持一致,同时为了更好地体现fNIRS携带的信息,将两种信号放缩到同一尺度;
步骤(2):使用sym3小波基函数对两种模态信号分别进行l层的小波分解,分别得到两种信号的一个近似系数al和细节系数dl;
两种信号离散小波变换是通过正交镜像滤波器组实现的;低通和高通滤波器分支的输出被视为近似系数和细节系数;通过将前一级的系数作为输入,同样的过程重复多次以提高频率分辨率;每次分解后,通过降采样使时间分辨率减半,而通过滤波使频率分辨率加倍;分解的第l层,其输出的小波系数为
其中g[k]与h[k]分别为低通和高通滤波器,m,k表示数据点;由于fNIRS信号是一种低频信号,在小波分解后的高频细节系数的幅值非常小,为了更好地体现fNIRS高频细节系数所携带的信息,本发明将两种信号的细节系数放缩到了同一尺度;
步骤(3):对每一层的小波系数计算其Fisher值,根据Fisher值给小波系数赋予权重,并将两种信号对应层的小波系数加权求和,从而得到一组新的小波系数,然后对新的小波系数进行小波重构,得到融合信号;
1):选择时间窗长度与EEG频带;
2):计算小波系数的Fisher值;
对不同的样本分别进行了I1与I2次实验,通过计算样本信号的各层小波系数的类间欧式距离与类内方差的比值,定义每个小波系数表征二分类能力的值为:
其中,F代表Fisher值,n代表通道数量,Xsi代表第s个通道第i次试验的小波系数,Msi代表Xsi的平均值;
3):融合规则
基于Fisher值加权与小波重构融合规则,首先计算两种信号小波系数的Fisher值,然后通过
来计算给各层小波系数分配的权重w*,再利用融合规则得到新的小波系数
Cnew=wEEGCEEG+wfNIRSCfNIRS (6)
最后使用新的小波系数Cnew进行小波重构得到融合信号;
步骤(4):提取融合信号的CSP特征,并使用LDA分类。
作为优选,所述的将光强度测量值转换为血液动力学变化;具体为:将在实验中接收到的近红外光的强度信号,通过修正的Beer-Lamber定律转换为含氧血红蛋白HbO和脱氧血红蛋白HbR的血液动力学变化,从而获取到大脑信息。
作为优选,采用实验前5秒的HbR和HbO浓度数据的均值进行基线校正。
作为优选,步骤(3)中的时间窗长度与EEG频带的选择,具体为:时间窗长度选择7s,EEG信号选择4-13Hz频带的信号。
作为优选,所述的提取融合信号的CSP特征,并使用LDA分类;具体为:
1):提取融合信号的CSP空域特征:通过一个投影矩阵使得不同类别的信号在投影空间有很大差异,然后求得每次实验信号X的投影信号Z,并用log-variances方法提取每次实验的CSP特征向量f;
为了最大化的区分信号,就需要找出能使两种信号经过投影后信号间的方差最大化的投影矩阵;假设两类任务的信号矩阵分别为X1和X2,构造投影矩阵;首先对两类样本数据归一化,然后求两类信号矩阵协方差Rj,得到:
其中,trace(·)是矩阵的迹,然后计算出混合空间的协方差矩阵R,并对其进行特征分解:
其中,是两类任务的协方差矩阵,j=1,2,Λ是总协方差矩阵的特征值,U是特征值对应的特征向量;对方差做均化处理,定义白化矩阵P为:
用该白化矩阵分别处理和/>
经过处理后的协方差矩阵用以下式子表示:
S1=UΛ1UT,S2=UΛ2UT (11)
其中,Λ1和Λ2分别为S1和S2对应的特征值,且满足Λ1+Λ2=I,I为单位矩阵;因此当S1的特征值取最大值的时候,S2中的对应的特征值也就最小,这样就使得两类信号能最大化的投影在不同的区域,实现分类;分别从两个特征向量中提取前后t个特征向量形成矩阵B,于是得到最后的投影矩阵W:
W=BTP (12)
通过投影矩阵得到两个类别信号投影后的信号矩阵:
Z=WX (13)
最后取对数得到最后的特征向量:
其中,Zu表示Z的第u行;
2):使用线性判别分析方法进行分类;
通过求解最优化问题来找出最佳投影矩阵W,即:
其中,trace(·)表示矩阵的迹,Sb和SW分别为类间和类内散布矩阵;使用从融合信号提取出来的CSP特征训练LDA分类器。
本发明的有益效果:本发明通过基于小波变换的方法实现的多种模态的数据级融合能够同时兼顾多种模态的信息,即从多个角度获取大脑信息,在一定程度上提高了BCI系统的分类性能。另外,与其他的多模态融合方法相比,比如特征级融合与决策级融合,数据级融合在特征提取之前就将多种模态融合,这能在一定程度上提高系统的计算效率。
附图说明
图1为本发明的流程图;
图2为单次实验时间轴;
图3为四层小波分解图;
图4为信号融合可视化图;
图5为不同时间窗长度与EEG频带分类比较;
图6为整体框架图;
具体实施方案
本发明使用由30个通道的EEG和36个通道的fNIRS信号所构成的数据集进行说明,数据集包含从29位健康受试者(14位男性和15位女性,平均年龄28.5±3.7)收集的心理算数实验信号。MA数据集有心算和基线任务,且每个受试者进行60次试验。图1为一次实验的具体时间轴,包括两秒钟的视觉介绍,10s的任务时间段和14-16s的休息时间。在MA任务的过程中,要求受试者从“三位数”中减去一位数字(例如123-9),并从先前的减法结果中重复减去一位数字,直到任务完成。在MA任务中,通过休息来执行基线任务。
下面结合附图详细描述本发明基于小波变换的多模态数据融合方法。本发明的基本流程如图2所示,具体包括4个步骤:
步骤(1):多通道EEG与fNIRS信号的采集与预处理。
4):频带处理,由于脑电信号中通常包含有非常多的干扰与噪声,其中一部分可以在信号采集中尽量避免,而有一部分则只能在后续的处理中去消除,如工频干扰等。所以首先对EEG信号进行0.5-35Hz的滤波,目的是为了消除干扰和噪声。
5):将光强度测量值转换为血液动力学变化,由于近红外信号是一种光信号,我们需要将在实验中接收到的近红外光的强度信号,通过修正的Beer-Lamber定律转换为含氧血红蛋白(Oxygenated Hemoglobin,HbO)和脱氧血红蛋白(Deoxygenated Hemoglobin,HbR)的血液动力学变化,从而获取到大脑信息。并且使用六阶零相位巴特沃斯滤波器对fNIRS信号进行0.01-0.1Hz的滤波操作,同时采用实验前5s的HbR和HbO浓度数据的均值进行基线校正。
6):插值与放缩:由于EEG和fNIRS信号的采样率不同,且两者的幅度差异较大,为了保证小波变换的兼容性和融合的效果,需要对fNIRS信号进行插值,使其与EEG的采样率保持一致,同时为了更好地体现fNIRS携带的信息,需要将两信号放缩到同一尺度。
步骤(2):使用sym3小波基函数对EEG和fNIRS信号分别进行4层的小波分解。
信号X(n)离散小波变换是通过正交镜像滤波器组实现的,如图3所示。低通和高通滤波器分支的输出被视为近似(al)和细节(dl)系数,其中i表示小波网络的层数。通过将前一级的系数作为输入,同样的过程可以重复多次以提高频率分辨率。每次分解后,通过降采样使时间分辨率减半,而通过滤波使频率分辨率加倍。在分解的第m层,其输出的小波系数为
其中g[k]与h[k]分别为低通和高通滤波器。由于fNIRS信号是一种低频信号,在小波分解后的高频细节系数的幅值非常小,为了更好地体现fNIRS高频细节系数所携带的信息,本发明将两种信号的细节系数放缩到了同一尺度。
步骤(3):按照融合规则进行融合。
1):计算小波系数的Fisher值。
Fisher准则是一种将高维参数投影成一维来测量类判别属性的统计量参数,Fisher值代表了两类样本的可分程度,Fisher值越大表示可分性越强,。若两种不同的样本分别进行了I1与I2次实验,通过计算样本信号的各层小波系数的类间欧式距离与类内方差的比值,定义每个小波系数表征二分类能力的值为:
其中,F代表Fisher值,n代表通道数量,Xsi代表第s个通道第i次试验的小波系数,Msi代表Xsi的平均值。图4将信号融合过程进行了可视化。
2):融合规则。
基于Fisher值加权与小波重构(Fisher weighted and Wavelet reconstruction,FW-WR)融合规则,首先计算两种信号小波系数的Fisher值,然后通过
来计算给各层小波系数分配的权重w*,再利用融合规则得到新的小波系数
Cnew=wEEGCEEG+wfNIRSCfNIRS (6)
最后使用新的小波系数Cnew进行小波重构得到融合信号,其中F*代表Fisher值。
3):时间窗长度与EEG频带的选择。
由于近红外信号具有较高的延迟性,为了选取到数据质量最高、可分性最好的的时间段,对融合信号进行滑动时间窗处理,时间窗长度的选择十分重要,太短则包含的信息量不足,太长则可能造成冗余。图5是时间窗长度与EEG频带不同的组合的分类效果,从中可以看出在时间窗长度选择7s,EEG选择4-13Hz也就是thete和alpha频段时的融合信号的分类效果是最好的,所以本发明最终选择7s的时间窗长度,EEG信号选择4-13Hz频带的信号与fNIRS融合进行研究。
步骤(4):提取融合信号的CSP特征,并使用LDA分类。
1):CSP特征提取:本发明提取的是融合信号的CSP空域特征,CSP是一种高效的空域滤波特征提取方法,已经广泛应用在人脸识别、疲劳检测、脑电信号分析、语音识别等领域。CSP算法提取特征的核心在于通过一个投影矩阵使得不同类别的信号在投影空间有很大差异,然后求得每次实验信号X的投影信号Z,并用log-variances方法提取每次实验的CSP特征向量f。
为了最大化的区分信号,就需要找出能使两种信号经过投影后信号间的方差最大化的投影矩阵。假设两类任务的信号矩阵分别为X1和X2,构造投影矩阵;首先求两类信号矩阵协方差,并归一化,得到:
其中,trace(·)是矩阵的迹,然后计算出混合空间的协方差矩阵R,并对其进行特征分解:
其中,是两类任务的协方差矩阵,j=1,2,Λ是总协方差矩阵的特征值,
U是特征值对应的特征向量。对方差做均化处理,定义白化矩阵P为:
用该白化矩阵分别处理和/>
经过处理后的协方差矩阵也可以用以下式子表示:
S1=UΛ1UT,S2=UΛ2UT (11)
其中,Λ1和Λ2分别为S1和S2对应的特征值,且满足Λ1+Λ2=I,I为单位矩阵。因此当S1的特征值取最大值的时候,S2中的对应的特征值也就最小,这样就使得两类信号能最大化的投影在不同的区域,实现分类。分别从两个特征向量中提取前后m个特征向量形成矩阵B,于是得到最后的投影矩阵W:
W=BTP (12)
通过投影矩阵就可以得到两个类别信号投影后的信号矩阵:
Z=WX (13)
最后取对数就可以得到最后的特征向量:
其中,Zu表示Z的第u行;
2):LDA分类:本发明使用的分类方法为线性判别分析(Linear discriminantanalysis,LDA)。LDA的主要思想是将一个高维空间中的数据投影到一个较低维的空间中,且投影后使样本的类内距离最小而类间距离最大,可以通过求解最优化问题来找出这样的一个最佳投影矩阵W,即:
其中,trace(·)表示矩阵的迹,Sb和SW分别为类间和类内散布矩阵。使用从融合信号提取出来的CSP特征训练LDA分类器,本发明的整体框架如图6所示。本发明所提出的方法与单模态EEG,fNIRS,多模态的特征级融合进行了分类比较,实验精度如表1所示:
表1 10名受试者不同模态特征的分类精度表
针对上述数据集,基于小波变换的数据级融合表现出来最佳的分类性能。来自数据集的十位受试者:基于小波变换的数据级融合在S1、S2、S3、S6和S7五位受试者身上获得了最高的分类进度,而特征级融合在S5、S9和S10身上获得了最高分类精度,单模态EEG则在S4和S8身上获得最高分类精度。基于小波变换的数据级多模态融合的方法作为一种结合了多种模态的方法,结合了各种模态的优势,从不同的角度获取大脑信息斌融合,增强了获取到的信息的丰富性,除此之外,本发明使用了基于Fisher值加权的融合规则,给可分性强的小波系数赋予更高的权重,可分性低的小波系数赋予小的权重,这即结合了两种模态的信息,同时又避免了不必要信息的冗余,这些都有利于在一定程度上提高BCI系统的分类性能。
从整体上看,本发明提出了一种基于小波变换的多模态信号融合方法用于BCI系统的任务分类。利用两种信号各层小波系数的Fisher值衡量各信号成分对于心算任务可分性的作用,然后采用基于Fisher值加权的融合规则分配各小波系数的权重,可分性越强的小波系数权重越大,由这种方式所得到小波系数重构出的信号同时兼顾了两种信号的信息,因而具有良好的可分性。最后的分类精度结果表明该方法为的多模态数据级融合提供了新思路,在一定程度上能够提高BCI系统的分类性能。
Claims (5)
1.一种基于小波变换的多模态脑机接口数据融合方法,其特征在于,该方法具体包括以下步骤:
步骤(1):多通道多模态信号的采集与预处理,对两种模态中每个通道信号的预处理是进行适带通滤波;
1):频带处理,对EEG信号进行0.5-35Hz的滤波,目的是为了消除干扰和噪声;
2):将光强度测量值转换为血液动力学变化,使用六阶零相位巴特沃斯滤波器对fNIRS信号进行0.01-0.1Hz的滤波操作,同时采用实验前N秒的HbR和HbO浓度数据的均值进行基线校正;
3):插值与放缩:对fNIRS信号进行插值,使其与EEG的采样率保持一致,将两种信号放缩到同一尺度;
步骤(2):使用sym3小波基函数对两种模态信号分别进行l层的小波分解,分别得到两种信号的一个近似系数al和细节系数dl;
两种信号离散小波变换是通过正交镜像滤波器组实现的;低通和高通滤波器分支的输出被视为近似系数和细节系数;通过将前一级的系数作为输入,同样的过程重复多次以提高频率分辨率;每次分解后,通过降采样使时间分辨率减半,而通过滤波使频率分辨率加倍;分解的第l层,其输出的小波系数为
其中g[k]与h[k]分别为低通和高通滤波器,m,k表示数据点;将两种信号的细节系数放缩到了同一尺度;
步骤(3):对每一层的小波系数计算其Fisher值,根据Fisher值给小波系数赋予权重,并将两种信号对应层的小波系数加权求和,从而得到一组新的小波系数,然后对新的小波系数进行小波重构,得到融合信号;
1):选择时间窗长度与EEG频带;
2):计算小波系数的Fisher值;
对不同的样本分别进行了I1与I2次实验,通过计算样本信号的各层小波系数的类间欧式距离与类内方差的比值,定义每个小波系数表征二分类能力的值为:
其中,F代表Fisher值,n代表通道数量,Xsi代表第s个通道第i次试验的小波系数,Msi代表Xsi的平均值;
3):融合规则
基于Fisher值加权与小波重构融合规则,首先计算两种信号小波系数的Fisher值,然后通过
来计算给各层小波系数分配的权重w*,再利用融合规则得到新的小波系数
Cnew=wEEGCEEG+wfNIRSCfNIRS (6)
最后使用新的小波系数Cnew进行小波重构得到融合信号;
步骤(4):提取融合信号的CSP特征,并使用LDA分类。
2.根据权利要求1所述的一种基于小波变换的多模态脑机接口数据融合方法,其特征在于:所述的将光强度测量值转换为血液动力学变化;具体为:将在实验中接收到的近红外光的强度信号,通过修正的Beer-Lamber定律转换为含氧血红蛋白HbO和脱氧血红蛋白HbR的血液动力学变化,从而获取到大脑信息。
3.根据权利要求1所述的一种基于小波变换的多模态脑机接口数据融合方法,其特征在于:采用实验前5秒的HbR和HbO浓度数据的均值进行基线校正。
4.根据权利要求1所述的一种基于小波变换的多模态脑机接口数据融合方法,其特征在于:步骤(3)中的时间窗长度与EEG频带的选择,具体为:时间窗长度选择7s,EEG信号选择4-13Hz频带的信号。
5.根据权利要求1所述的一种基于小波变换的多模态脑机接口数据融合方法,其特征在于:所述的提取融合信号的CSP特征,并使用LDA分类;具体为:
1):提取融合信号的CSP空域特征:通过一个投影矩阵使得不同类别的信号在投影空间有很大差异,然后求得每次实验信号X的投影信号Z,并用log-variances方法提取每次实验的CSP特征向量f;
为了最大化的区分信号,就需要找出能使两种信号经过投影后信号间的方差最大化的投影矩阵;假设两类任务的信号矩阵分别为X1和X2,构造投影矩阵;首先对两类样本数据归一化,然后求两类信号矩阵协方差Rj,得到:
其中,trace(·)是矩阵的迹,然后计算出混合空间的协方差矩阵R,并对其进行特征分解:
其中,是两类任务的协方差的均值,j=1,2,Λ是总协方差矩阵的特征值,U是特征值对应的特征向量;对方差做均化处理,定义白化矩阵P为:
用该白化矩阵分别处理和/>
经过处理后的协方差矩阵用以下式子表示:
S1=UΛ1UT,S2=UΛ2UT (11)
其中,Λ1和Λ2分别为S1和S2对应的特征值,且满足Λ1+Λ2=I,I为单位矩阵;因此当S1的特征值取最大值的时候,S2中的对应的特征值也就最小,这样就使得两类信号能最大化的投影在不同的区域,实现分类;分别从两个特征向量中提取前后t个特征向量形成矩阵B,于是得到最后的投影矩阵W:
W=BTP (12)
通过投影矩阵得到两个类别信号投影后的信号矩阵:
Z=WX (13)
最后取对数得到最后的特征向量:
其中,Zu表示Z的第u行;
2):使用线性判别分析方法进行分类;
通过求解最优化问题来找出最佳投影矩阵W,即:
其中,trace(·)表示矩阵的迹,Sb和SW分别为类间和类内散布矩阵;使用从融合信号提取出来的CSP特征训练LDA分类器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111234957.9A CN114145752B (zh) | 2021-10-22 | 2021-10-22 | 一种基于小波变换的多模态脑机接口数据融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111234957.9A CN114145752B (zh) | 2021-10-22 | 2021-10-22 | 一种基于小波变换的多模态脑机接口数据融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114145752A CN114145752A (zh) | 2022-03-08 |
CN114145752B true CN114145752B (zh) | 2024-03-29 |
Family
ID=80458638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111234957.9A Active CN114145752B (zh) | 2021-10-22 | 2021-10-22 | 一种基于小波变换的多模态脑机接口数据融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114145752B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020042511A1 (zh) * | 2018-08-28 | 2020-03-05 | 天津大学 | 基于空间滤波与模版匹配的运动电位脑机接口编解码方法 |
CN111931576A (zh) * | 2020-07-07 | 2020-11-13 | 杭州电子科技大学 | 基于ewt和改进csp的癫痫信号分类方法 |
CN112560703A (zh) * | 2020-12-18 | 2021-03-26 | 杭州电子科技大学 | 一种基于pf系数的多模态bci特征提取方法 |
CN113158793A (zh) * | 2021-03-15 | 2021-07-23 | 东北电力大学 | 一种基于多特征融合的多类运动想象脑电信号识别方法 |
CN113208613A (zh) * | 2021-04-21 | 2021-08-06 | 杭州电子科技大学 | 基于fhls特征选择的多模态bci时序优化方法 |
-
2021
- 2021-10-22 CN CN202111234957.9A patent/CN114145752B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020042511A1 (zh) * | 2018-08-28 | 2020-03-05 | 天津大学 | 基于空间滤波与模版匹配的运动电位脑机接口编解码方法 |
CN111931576A (zh) * | 2020-07-07 | 2020-11-13 | 杭州电子科技大学 | 基于ewt和改进csp的癫痫信号分类方法 |
CN112560703A (zh) * | 2020-12-18 | 2021-03-26 | 杭州电子科技大学 | 一种基于pf系数的多模态bci特征提取方法 |
CN113158793A (zh) * | 2021-03-15 | 2021-07-23 | 东北电力大学 | 一种基于多特征融合的多类运动想象脑电信号识别方法 |
CN113208613A (zh) * | 2021-04-21 | 2021-08-06 | 杭州电子科技大学 | 基于fhls特征选择的多模态bci时序优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114145752A (zh) | 2022-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Barachant et al. | Riemannian geometry applied to BCI classification | |
Shahbakhti et al. | Simultaneous eye blink characterization and elimination from low-channel prefrontal EEG signals enhances driver drowsiness detection | |
KR20190128978A (ko) | 인간 감정 인식을 위한 딥 생리적 정서 네트워크를 이용한 인간 감정 추정 방법 및 그 시스템 | |
Sawangjai et al. | EEGANet: Removal of ocular artifacts from the EEG signal using generative adversarial networks | |
Asgher et al. | Classification of mental workload (MWL) using support vector machines (SVM) and convolutional neural networks (CNN) | |
CN114533086B (zh) | 一种基于空域特征时频变换的运动想象脑电解码方法 | |
Wang et al. | A 2D CNN-LSTM hybrid algorithm using time series segments of EEG data for motor imagery classification | |
CN113191395A (zh) | 一种基于双脑的多层次信息融合的目标检测方法 | |
Ergün et al. | Decoding of binary mental arithmetic based near-infrared spectroscopy signals | |
Song et al. | Common spatial generative adversarial networks based EEG data augmentation for cross-subject brain-computer interface | |
CN113208613B (zh) | 基于fhls特征选择的多模态bci时序优化方法 | |
Chaddad | Brain function diagnosis enhanced using denoised fNIRS raw signals | |
Wang et al. | Multiband decomposition and spectral discriminative analysis for motor imagery BCI via deep neural network | |
CN112560703B (zh) | 一种基于pf系数的多模态bci特征提取方法 | |
Ghonchi et al. | Spatio-temporal deep learning for EEG-fNIRS brain computer interface | |
CN114145752B (zh) | 一种基于小波变换的多模态脑机接口数据融合方法 | |
Ergün et al. | A Hybrid BCI Using Singular Value Decomposition Values of the Fast Walsh–Hadamard Transform Coefficients | |
CN111096730A (zh) | 基于自发动力学活动的波动熵的自闭症分类方法 | |
Khoa et al. | Recognizing brain activities by functional near-infrared spectroscope signal analysis | |
Silva et al. | Towards a LMI approach to feature extraction improvements and classification by Riemann distance | |
Meng et al. | Crossing time windows optimization based on mutual information for hybrid BCI | |
Liu et al. | fNIRS feature extraction and classification in grip-force tasks | |
Hanafi et al. | Brain Computer Interface (BCI) Machine Learning Process: A Review | |
Jayashekar et al. | Hybrid Feature Extraction for EEG Motor Imagery Classification Using Multi-Class SVM. | |
CN110363242A (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 |