一种功能磁共振成像辅助的脑电通道选择方法
【技术领域】
[0001]该发明涉及一种脑电通道选择方法，特别是涉及一种功能磁共振成像辅助的脑电通道选择方法。
【背景技术】
[0002]近二十年来，随着脑影像技术的进步，脑科学的研究进入了一个高速发展时期。作为无创的进行脑功能研究的两种主要神经成像工具，脑电(Electroencephalography，EEG)和功能磁共振成像(funct1nal 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 Θ + ε，其中γ e RnXs代表η个电极s个采样点的数据，L e RnXd代表由EEG正演模型得到的传递矩阵，Θ e RdXs表示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标准解剖空间中的不同位置；利用边界元法，将得到的三层头模型输入到Fieldtrip软件中，就得到从皮层网格到头表的传递矩阵L e RnXd，η表示电极个数，d表示EEG信号源个数。
[0014]所述步骤(3)由步骤(1) fMRI任务数据的分析，可得到与特定脑功能相关脑区的激活程度，使用t-检验的t值表示，由步骤(2)构建的EEG正演模型，可得到从皮层到头表电极处的传递矩阵；在EEG正演模型的建立过程中，对皮层进行了离散网格化处理，对于每个离散网格，使用其邻域3 X 3体素t值绝对值的均值代表网格的激活情况，可得到8196维代表皮层网格激活程度的向量a，通过传递矩阵，就可得到各个通道和特定脑功能的相关程度值 E = I L.a I。
[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 e RnXd为由皮层到头表的传递矩阵，η为电极数目，d为皮层上的网格数，a为大脑皮层网格处的激活程度向量，E为最终得到的各个电极与特定脑功能的相关性程度。首先利用fMRI实验数据得到的特定脑功能相关脑区的激活情况，再通过由大脑标准结构像构建的EEG正演模型得到传递矩阵，从而得到反映与特定脑功能相关程度的脑电地形图，为基于脑电脑特定功能研究的通道选择提供指导。
[0026]方法具体流程如附图1所示，主要分为四个步骤:
[0027](1)由fMRI实验数据获取相关功能脑区的激活情况
[0028]设计反映大脑相关功能的fMRI任务实验，使用任务态的f