注:本教程使用了Matlab中的Brainstorm工具箱,依据Brainstorn官网教程进行,更多具体细节可参考Brainstorm官网。教程中使用个人采集的同步EEG-fMRI数据,fMRI数据包括T1结构像数据和EPI功能像数据。其中涉及到的工具包:CAT12、OpenMEEG均可于Brainstorm中自动下载。
一、使用CAT12分割MRI
1、在Matlab中运行Brainstorm,在解剖端创建一个新的Protocol:点击File->New protocol。由于采集了被试的结构像数据,选择个人模板。
2、点击Create,此时界面增加了一些默认的结构数据,如下图所示:
3、右键点击新建的Protocol,选择New subject
4、新建Subject,设置被试名(这里以8liping为例),后面的选择与建立Protocol时一致。
5、点击Save,右击8liping,选择Import MRI。导入其结构像数据nii格式文件。
6、计算MNI归一化。即点击下图绿框处。
7、点击绿框处,得到如下图,选择maff8,稍等,然后Save。
8、右击s-8liping-3D->MRI segmentation->CAT12。后面默认OK。
9、这步时间较长,生成的结果将自动导入,结果如下:
二、使用基于OpenMEEG的BEM法构建真实头模型
该正演模型使用对称边界元法(对称BEM),为了计算一个真实的边界元正向模型,我们需要受试者的解剖结构包括代表模型中使用的组织之间边界的表面:头皮(头部-空气界面)、外颅骨(头皮-颅骨界面)和内颅骨(颅骨-大脑+脑脊液界面)。
1、计算前向分析模型所需的文件:MRI、头部表面、皮层表面。在协议的解剖视图中,右键单击MRI segmentation->Generate BEM surfaces,选择Brainstorm。
Brainstorm可以根据受试者的皮层和头部表面以及ICBM152的内颅面和外颅面生成粗略的近似值。生成的曲面是非相交的。
2、切换到视图“功能数据(按主题排序)”,数据库浏览上方的第二个按钮。右键单击通道文件->计算头模型,在列表中选择OpenMEEG BEM。
本例遇到如下问题:
上示问题表示偶极子外露,即通过BEM生成的头皮层与外颅骨相交了,如下图所示。
解决方法:偶极子外露解决方案
重新生成具有更多顶点的BEM曲面来更好地进行包裹。本例中,将原来1922的Scalp增加至2562.。
计算头部模型结束后,将会生成新头部模型,如下图所示。