技术领域
[0001] 该
发明涉及一种脑电通道选择方法，特别是涉及一种功能磁共振成像辅助的脑电通道选择方法。
背景技术
[0002] 近二十年来，随着脑影像技术的进步，脑科学的研究进入了一个高速发展时期。作为无创的进行脑功能研究的两种主要神经成像工具，脑电(Electroencephalography，EEG)和功能磁共振成像(functional magnetic resonance imaging,fMRI)分别以其高时间
分辨率和高空间分辨率的优势在临床诊断和学术研究中得到了高度关注和广泛应用。
[0003] 在基于EEG的特定脑功能研究中，利用少量通道实现特定脑功能数据
信号可靠有效的获取是EEG研究向便携化、实用化发展的关键。利用特征选择的方法虽然可以挑选出对大脑相关功能信息获取较为有效的通道，但需要提前进行一定量的
数据采集和分析，且基于特征选择的通道选择方法理论依据并不充分。作为神经影像领域的一项重要技术，fMRI高空间分辨率的特性，可以在一定程度上弥补EEG所获取空间信息的不足。如何利用由fMRI得到的大脑激活空间信息来辅助EEG通道的选择，是本
专利要尝试解决的问题。
[0004] 目前，EEG与fMRI的融合研究已成为了前沿热点。这方面的研究主要分为三个方向：利用fMRI信息对EEG信号的分析进行空间约束；利用EEG信息对fMRI信号的分析进行时间预测；建立共同的生成模型对两者数据进行对称融合。
[0005] 在基于fMRI信息的EEG分析中，大多数的研究针对的是EEG溯源问题(EEG反问题)，即以EEG正演模型为
基础，采用fMRI的空间信息来重建EEG的源分布。其基本思想由如下公式描述：Y＝Lθ+ε，其中Y∈Rn×s代表n个
电极s个
采样点的数据，L∈Rn×d代表由EEG正演模型得到的传递矩阵，θ∈Rd×s表示d个EEG信号源的动态过程。EEG溯源问题，则是已知Y和L，求θ。这是一个典型的不适定问题。利用fMRI所提供的空间信息，可以为θ提供一定的先验信息，从而使求得的解更加的符合真实情况。
[0006] 根据利用fMRI空间信息方式的不同，基于fMRI信息的EEG分析方法主要分为三种：一是直接利用fMRI的激活脑区
位置直接约束EEG的源位置，或者对偶极子的
种子点进行初始化；二是引入贝叶斯
框架来放宽激活脑区对EEG源位置的约束；三是利用fMRI的功能网络作为EEG源
定位的先验。
发明内容
[0007] 本发明克服了
现有技术中，单纯依靠脑电数据进行脑电通道选择所面临的空间分辨率低的限制，提供一种效果优越的功能磁共振成像辅助的脑电通道选择方法。含有以下四个步骤：
[0008] (1)由fMRI实验数据获取相关功能脑区的激活情况；
[0009] (2)由大脑标准结构像构建EEG正演模型；
[0010] (3)通过EEG正演模型计算各个通道和特定脑功能的相关程度；
[0011] (4)根据得到的脑功能相关程度地形图对EEG通道进行选择。
[0012] 所述步骤(1)使用任务态的fMRI
数据处理方法，得到被试在做特定脑功能任务时相关脑区的激活情况。
[0013] 所述步骤(2)通过大脑标准结构像建立一个标准的头模型，每个个体的数据图像则通过SPM的标准化步骤配准到标准脑上；即通过SPM
软件得到标准大脑的头皮、颅骨和皮层三层球模型，并对皮层进行网格化，取8196个网格，每个格点代表一个偶极子，电极配准到头表；每个网格对应了Talairach标准解剖空间中的不同位置；利用边界元法，将得到的n×d三层头模型输入到Fieldtrip软件中，就得到从皮层网格到头表的传递矩阵L∈R ，n表示电极个数，d表示EEG信号源个数。
[0014] 所述步骤(3)由步骤(1)fMRI任务数据的分析，可得到与特定脑功能相关脑区的激活程度，使用t-检验的t值表示，由步骤(2)构建的EEG正演模型，可得到从皮层到头表电极处的传递矩阵；在EEG正演模型的建立过程中，对皮层进行了离散网格化处理，对于每个离散网格，使用其邻域3×3
体素t值绝对值的均值代表网格的激活情况，可得到8196维代表皮层网格激活程度的向量a，通过传递矩阵，就可得到各个通道和特定脑功能的相关程度值E＝|L·a|。
[0015] 所述步骤(4)根据步骤(3)得到各个通道和特定脑功能的相关程度值，对每个通道和特定脑功能的相关程度进行排序，依据相关程度值从大到小的顺序进行选择。
[0016] 与现有技术相比，本发明功能磁共振成像辅助的脑电通道选择方法具有以下优点：利用了fMRI技术高空间分辨率的优势，在一定程度上突破了脑电通道选择中EEG空间分辨率低的限制；相比于传统的依靠经验或数据分析进行通道选择，更有理论依据；可针对不同人，制定出不同的通道选择方法。此外，本发明方法可适用于不同的脑功能研究，具有一定的通用性。
附图说明
[0017] 图1是本发明功能磁共振成像辅助的脑电通道选择方法中基于EEG正演模型的脑电特定脑功能通道选择
流程图；
[0018] 图2是本发明功能磁共振成像辅助的脑电通道选择方法中的大脑三层头模型结构示意图；
[0019] 图3是本发明功能磁共振成像辅助的脑电通道选择方法中的皮层网格对应到Talairach标准解剖空间位置图；
[0020] 图4是本发明功能磁共振成像辅助的脑电通道选择方法中的情绪识别实验范式设计示意图；
[0021] 图5是本发明功能磁共振成像辅助的脑电通道选择方法中被试1强相关性通道筛选前后及对比实验的分类正确率示意图；
[0022] 图6是本发明功能磁共振成像辅助的脑电通道选择方法中被试2强相关性通道筛选前后及对比实验的分类正确率示意图；
[0023] 图7是本发明功能磁共振成像辅助的脑电通道选择方法中被试3强相关性通道筛选前后及对比实验的分类正确率示意图。
具体实施方式
[0024] 下面结合附图和具体实施方式对本发明功能磁共振成像辅助的脑电通道选择方法作进一步说明，本发明的技术方案为：EEG正演模型反映的是皮层处信号源电位变化与头表处电位信号变化之间的关系，若假设由fMRI得到的相关功能激活脑区便是主要对头表处电位变化有贡献的脑区，那么由脑电的正演模型，便可以将fMRI得到的相关功能激活脑区映射到头表处，以反映不同电极位置处与特定脑功能相关程度的大小，从而来辅助进行EEG通道的选择。
[0025] 方法思想可以由如下公式描述：E＝L·a，其中L∈Rn×d为由皮层到头表的传递矩阵，n为电极数目，d为皮层上的网格数，a为大脑皮层网格处的激活程度向量，E为最终得到的各个电极与特定脑功能的相关性程度。首先利用fMRI实验数据得到的特定脑功能相关脑区的激活情况，再通过由大脑标准结构像构建的EEG正演模型得到传递矩阵，从而得到反映与特定脑功能相关程度的脑电地形图，为基于脑电脑特定功能研究的通道选择提供指导。
[0026] 方法具体流程如附图1所示，主要分为四个步骤：
[0027] (1)由fMRI实验数据获取相关功能脑区的激活情况
[0028] 设计反映大脑相关功能的fMRI任务实验，使用任务态的fMRI数据处理方法，得到被试在做特定脑功能任务时相关脑区的激活情况。
[0029] (2)由大脑标准结构像构建EEG正演模型
[0030] 对于每个个体而言，通过对其MRI结构像中进行
图像分割来精确地建立3层边界元头模型并不容易。本发明的做法是通过大脑标准结构像建立一个标准的头模型，而每个个体的数据图像则通过SPM的标准化步骤配准到标准脑上。首先，通过SPM软件得到标准大脑的头皮、颅骨和皮层三层球模型，并对皮层进行网格化(本文取了8196个网格，每个格点代表一个偶极子)，电极配准到头表(如附图2)。每个网格对应了Talairach标准解剖空间中的不同位置(图3中红色体素位置)。利用边界元法(Boundary Element Method，BEM)，将得到的三层头模型输入到Fieldtrip软件中，就可得到从皮层网格到头表的传递矩阵L∈Rn×d。
[0031] (3)通过EEG正演模型计算各个通道和特定脑功能的相关程度
[0032] 由第(1)步fMRI任务数据的分析，可得到与特定脑功能相关脑区的激活程度(使用t-检验的t值表示)，由第(2)步构建的EEG正演模型，可得到从皮层到头表电极处的传递矩阵。在EEG正演模型的建立过程中，为了便于模型的求解，对皮层进行了离散网格化处理，所以这里需要将特定脑功能相关脑区的激活情况映射到离散的网格上。对于每个离散网格，使用其邻域3×3体素t值绝对值的均值代表网格的激活情况，可得到8196维代表皮层网格激活程度的向量a，通过传递矩阵，就可得到各个通道和特定脑功能的相关程度值E＝L·a，这里取绝对值表示只关注相关程度的大小。
[0033] (4)根据得到的脑功能相关程度地形图对EEG通道进行选择
[0034] 步骤(4)根据步骤(3)得到各个通道和特定脑功能的相关程度值，对每个通道和特定脑功能的相关程度进行排序，依据相关程度值从大到小的顺序进行选择。代表n个电极sn×d d×s个采样点的数据，L∈R 代表由EEG正演模型得到的传递矩阵，θ∈R 表示d个EEG信号源的动态过程。
[0035] 将本发明方法应用在具体情绪识别实验中，如下：在中国情感图片系统(Chinese Affective Picture System,CAPS)中依据愉悦度选取强正，弱正，中性，弱负，强负5种类型图片作为情绪刺激材料，使得五类图片在愉悦度上有统计意义上的显著差异(两两间双样本F检验，p<0.05)。
[0036] 整个情绪识别实验包括150个trial，每张图片呈现一次。一个trial的具体设计如附图4所示，首先是4s的静息，屏幕中出现“静息”提示，接着是2s的准备阶段，期间要求被试盯着屏幕中间的十字并保持注意
力。随后是4s的情绪诱发图片呈现阶段，要求被试在此期间充分感受图片中的内容，以便在接下来4s的评估阶段给出图片的愉悦度评分，强负、弱负、中性、弱正、强正5种情绪体验分别对应1～5数字键。被试按键反应后，屏幕中“评估”两字会消失，表明按键成功。整个实验刺激序列通过E-prime软件编程完成。
[0037] 选取3名被试分别在进行EEG和fMRI实验。EEG数据采集使用奥地利gtec脑电公司的16通道g-USBamp系统，fMRI数据采集使用
磁场强度为3.0T的GE Discovery MRI 750磁共振
扫描仪获取。
[0038] 在情绪识别实验数据的分析中，选择在情绪识别中最为广泛应用的
功率谱特征。使用相关性特征选择方法和基于RBF核的SVM分类器进行5折交叉检验。
[0039] 运用本发明方法进行EEG通道选择。比较了三种情况下的情绪情绪状态分类正确率：1、未进行通道筛选前；2、运用本发明方法筛选7个和情绪最相关通道；3、在剩余通道中随机选择7个通道(进行3次，求均值)。结果如附图5、6、7所示，表明由本发明方法筛选通道后的分类正确率虽然比为筛选前有所下降，但是明显高于随机选择其他通道下的结果，即表明由本发明方法筛选出来的通道确实比其他通道表现出来了较高的情绪相关性。
