本文是我进入实验室后,针对原本的配准工具fieldtrip在MNI空间映射上不准确的问题,改用的新工具Brainstorm的教程,但因为该工具基本没有可编程部分,因此依然保留了使用Fieldtrip进行前期acpc定位和批处理导出电极触点所在分区的功能,希望对有需要的同学有所帮助。实验室对于Fieldtrip的工作流程很熟悉,因此在这里没有详述,如果不了解的话请阅读

Integrated analysis of anatomical and electrophysiological human intracranial data | Nature Protocols

出于对病人隐私的考虑,原始MRI、CT图像无法在这里给出,请在其他.nii数据或这些软件给出的公共数据(不推荐,国外的IEEG电极与我们的不相同)进行尝试,在开始之前使用freesurfer进行分割重建

目录

一、前言

二、教程正文

1. fieldtrip前期配准:

2. 创建工程:

3. 解剖结构的导入及准备:

4. 导入植入后CT图像:

5. 标记植入电极:

6. 确定电极触点分区:

三、示例代码

四、 资源下载


一、前言

工具版本:Fieldtrip-20210614

Brainstorm-210711

本教程为IEEG电极植入后定位的标准化流程,为后续分析提供参考,包含MRI-CT图像的配准、IEEG电极定位、IEEG触点在不同坐标系(acpc、ctf、mni)的坐标定位及IEEG触点在不同模板(AAL、Desikan-Killiany、Destrieux等)的分区定位,在之前使用fieldtrip工作流程的基础上,针对于原方法MNI映射不准确、图像可视化差等问题,引入brainstorm工具予以解决。本教程以brainstorm为主,fieldtrip仅作为辅助作用,研究人员可在此基础上,继续使用brainstorm对SEEG信号及癫痫发病的相关问题予以进一步分析,该工具也为此提供了丰富的方法支持。本教程针对的研究人员对于之前使用Fieldtrip的工作流程有大致了解,因此仅给出示例代码予以参考,对于fieldtrip配准、MRI图像基准点的定位、freesurfer分割等具体操作不再做说明。

brainstorm下载及安装教程请参考:

EEG/MEG数据处理工具包Brainstorm的下载与安装教程_悦影科技-CSDN博客_brainstorm安装

brainstorm详细教程请参考:

Tutorials - Brainstorm

二、教程正文

1. fieldtrip前期配准:

使用fieldtrip将原始MRI图像配准到acpc坐标系

2. 创建工程:

将matlab工作目录切换至brainstorm3安装文件夹,命令行中输入“brainstorm”启动,首次启动需要进行安装

Plugins-spm12-install,安装spm工具箱。brainstorm提供了包括spm12、fieldtrip在内的多种工具的集成、调用,若您的电脑中已经安装了spm12,直接点击Load进行加载、

File-New protocol,为患者创建新的协议,在这里我们以编号-患者姓名的格式为其命名,Brainstorm数据库是由一组协议组成的,一次只能打开一个协议,保留所有默认设置创建

File-New subject创建一个新主题,一个协议中可以包含多个主题,依然保留所有默认设置创建

3. 解剖结构的导入及准备:

右键Subject-Import anatomy folder,导入freesurfer分割文件夹,文件类型选择FreeSurfer+Volume atlases(除头部表面外也导入了体积分区),系统提示输入最终皮质表面所需的顶点数,默认情况下,我们将整个大脑的这个值设置为 15000(这意味着每个半球为 7500)。在导入过程中,实际只需要freesurfer文件夹中的label、surf、mri这三个文件夹

在MRI查看器中,定位NAS(鼻根),左耳 (LPA),右耳 (RPA),前连合 (AC),后联合(PC),半球间点 (IH),我们默认您已经具备了这些定位的相关知识,在此不过多赘述,brainstorm内部已经舍弃了acpc相关坐标,仅作为ctf坐标系的辅助。Brainstorm导入的MRI图像为freesurfer文件夹中的T1.mgz,由于步骤1中,我们已将该图像配准到acpc坐标系下,因此world坐标系原点即为ac点,尽管acpc定位不会对brainstorm产生过多影响,我们仍建议您这么做

完成基准的选择后,点击Save,brainstorm将进行后续的导入工作,并生成皮层的3D视图

右键MRI-MNI normalization,将坐标注册到MNI空间中。在弹出的窗口中选择maff8,此方法嵌入在 Brainstorm 中,不需要安装 SPM12,它基于与 SPM 软件中的 MNI ICBM152 模板的仿射联合配准

这种标准化不会修改解剖结构,它只是保存了一个转换,可以在头脑风暴坐标和 MNI 坐标之间进行转换。计算此转换后,我们可以在MRI 查看器中看到某点的MNI坐标

右键皮层cortex_15002V-MRI registration-Check MRI/surface registration,我们可以检查皮层表面与MRI图像是否吻合

4. 导入植入后CT图像:

右键Subject-Import MRI,选择CT图像,弹出窗口均选择默认选项,在Registration: How to register the new volume with the previous one?窗口中,我们通常选择第一个选项SPM,用 SPM 工具箱来估计从CT图像到MRI的变换,若您已经在其他地方(如fieldtrip)已将CT与MRI共同注册配准,这里也可选择Ignore

MRI查看器会自动打开,CT图像会作为彩色层覆盖到MRI图像,我们可以在此查看两个图像的配准是否准确

5. 标记植入电极

右键CT图像-SEEG/ECOG implantation

Subject中将会创建新的implantation文件夹,并开始编辑一个空的channel文件。在brainstrom窗口的右侧,iEEG选项卡下,我们可以添加电极(+符号)并设置其模型,请根据患者电极植入的相关信息,添加相应数量的电极,并设置电极长度、触点间距等信息。在本示例中,共有14根植入电极,其中4号(d)为长电极,其他均为短电极。默认情况下,触点直径略大于电极直径,否则有可能会导致触点的显示错误

在iEEG选项卡中选中要标记的电极,在MRI查看器中,将光标置于电极的尖端(颅内),点击选项卡中“Set tip”,再将光标放在头骨中电极的入口处(也可放在其他位置,该点仅用于计算方向,不使用该点的实际位置),点击Set skull entry,电极会出现在MRI查看器中,各触点位置会自动计算标出

在这一步中,您可以同时打开MRI图像与病例对照以获得电极位置的参考,由于MRI图像与CT图像已经配准,因此两个图像电极的坐标完全相同。在MRI查看器中,您可以通过Neurological/Radiological选项,在神经科医生和影像科医生的影像查看习惯中切换,以符合病例的左右方向

完成所有标记后,点击save,保存电极位置信息

右键SEEG/ECOG,选择在皮层模型下查看电极,以大致检查电极位置是否正确,在查看器中,每个触点都可被选中,右键触点可查看其详细信息

右键SEEG/ECOG-File-Export to file,将电极触点信息导出,在文件类型中选择EEG:ASCII:Name,XYZ_MNI导出MNI坐标,选择EEG:ASCII:Name,XYZ_World导出World坐标,由于我们使用的是已经过acpc配准的图像,因此World坐标也就是acpc坐标。导出的坐标是两个txt文件

6. 确定电极触点分区

在brainstorm里,可以很方便地查看各类大脑分区信息,freesurfer已为我们提供Deskian-Killiany、Destrieux等多种native坐标系下的分区信息,右键Subject-Add MNI parcellation,下载基于MNI空间的分区信息,在之前的步骤中我们已将图像注册到标准MNI空间,这些分区信息将基于此与图像对应,在这里我们下载AAL3分区

切换回电极信息页面,右键SEEG/ECOG,在MRI查看器中查看电极,可以看到,MRI查看器右上角已经显示光标位置所在区域,右键空白处-Anatomical atlas可以切换当前所使用的分区,“Show atlas“选项可以在MRI图像上以不同颜色显示各个分区。通过iEEG选项卡切换电极,可以很方便地查看各个位置所在的区域

使用示例代码,导入电极的MNI、World(acpc)坐标,调用generate_electable_v3方法,计算各电极触点分区信息并储存在excel表格中

、示例代码

Brainstorm+FieldtripiEEG分区定位示例代码.rar_fieldtrip的MNI模板-其它文档类资源-CSDN下载

四、 资源下载

鼓励您在原本的网站注册下载,与全世界数千名研究者共同讨论,若存在困难,我上传了brainstorm和fieldtrip的最新版本:

brainstorm:brainstorm_210711_fieldtrip的MNI模板-医疗文档类资源-CSDN下载

fieldtrip:fieldtrip-20210614.zip-医疗文档类资源-CSDN下载

Brainstorm + Fieldtrip IEEG定位及分区教程相关推荐

  1. win10一按右键就闪屏_一步一步的详细讲解Win10磁盘分区教程

    大家好,今天给大家带来的是 W10的磁盘分区教程,这个教程已经很常见了. 但是对于电脑小白而言,还是有很多小白不懂得,感觉磁盘分区会很复杂,希望这篇教程可以给予他们一些帮助. 本篇教程,目录如下: 目 ...

  2. linux命令管理GPT分区,Linux磁盘管理GPT分区教程

    Linux内核代码量大.逻辑关系复杂,因此对内核中的错误进行追溯和调试一直以来都是一件既耗费时间又耗费精力的事情.接下来是小编为大家收集的 Linux磁盘管理GPT分区教程,希望能帮到大家. Linu ...

  3. 计算机WIN7动态硬盘分区,win7硬盘分区教程

    win7硬盘分区教程 硬盘分区其实是对硬盘的一种格式化,进行了此步骤后才可以使用硬盘保存各种数据,下面是小编为大家详细介绍win7硬盘分区教程,欢迎大家阅读! 一.进入磁盘管理工具 我们点击要对硬盘进 ...

  4. 美妆短视频的定位分享,教程、技巧都可做,还可尝试变装内容

    美妆短视频的定位分享,教程.技巧都可做,还可尝试变装内容 爱美之心人皆有之,因此,美妆短视频的受众其实还是很广的,除了广大的爱美女孩之外,也有很多精致的男士是美妆短视频的受众,而这也是很多博主选择做美 ...

  5. WIN10解包分区和磁盘分区教程

    From: https://thinkpad.bbs.taobao.com/detail.html?postId=5390151 新机到手第一件事情就是要系统解包和磁盘分区.建议在安装和卸载软件之前, ...

  6. 服务器电脑增加机械硬盘分区,台式电脑加新、旧机械硬盘分区教程 旧机械硬盘可以装新电脑上吗?...

    有不少新电脑用户为了节省几百元的预算,将家里的老电脑的旧机械硬盘卸下安装到新电脑上,来继续使用.那么旧机械硬盘可以装新电脑上吗?下面分享一下台式电脑加新/旧机械硬盘+分区教程,无论是增加新机械硬盘或者 ...

  7. diskgenius linux 分区,DiskGenius怎么分区,DiskGenius分区教程

    DiskGenius分区教程 方法一: 一.首先,打开DiskGenius软件并单击以选择要操作的硬盘驱动器. 二.然后右键单击"快速分区"按钮并单击"确定". ...

  8. 用计算机管理从新分区,电脑如何分区硬盘分区_电脑怎么重新分区教程-win7之家...

    在给电脑安装系统的过程中,一些粗心大意的用户总是会忘记安装过程中的某些操作,比如分区硬盘的步骤,当然也存在着有的用户不知道如何分区硬盘比较合适的问题,因此常常会在后续操作中来研究,那么电脑如何分区硬盘 ...

  9. 7070mt安装ubuntu dell_戴尔(dell)7070MT台式机装win7系统及分区教程(支持usb驱动)

    [文章导读]戴尔OptiPlex 7070MT是一款商用台式机系列,该系列采用Intel 酷睿第八代处理器.默认预装了win10系统,最近有很多网友改了win7后发现进不了系统,由于默认戴尔OptiP ...

最新文章

  1. 这只狗,其实是猫变的:“撸猫神器”StarGAN v2来了!
  2. leetcode_zigzag conversion
  3. Lua 操作 MongoDB 数据库实例
  4. python matplotlib_Python の Matplotlib
  5. SpringBoot实战教程(8)| 整合mybatis-plus
  6. 论文A Tutorial on Graph-Based SLAM 学习笔记
  7. 北京特9内环和外环的区别_2021年2月CFA北京机考考点在哪里?
  8. python 颜色调整(饼状图)
  9. Python中requests上传大文件
  10. 多目标启发式算法(NSGA2, MOEA, MOPSO)python实现
  11. 计算机视觉专业名词中英文对照(转)
  12. 金字塔简单代码(java)
  13. CTU Open Contest 2019 J. Beer Vision
  14. C语言编译预处理#define(x,y)语句中x*y 、(x)*(y)、(x*y)的问题
  15. 【GMS认证】关于GSI和VTS解锁刷镜像命令总结
  16. 美国大学计算机科学gpa,美国大学gpa计算
  17. CDR类CAD制作室内装修平面图
  18. Python计算分位点与逆运算:根据给定的值,计算在序列中的分位水平
  19. YUV420 Planar 和YUV420 Semi-Planar详解
  20. 实训五#5.5猫与餐厅的故事

热门文章

  1. linux修改文件元信息,Linux 文件元数据详细讲解
  2. Gps车辆监控软件平台开发技术精华
  3. 腾讯云轻量应用服务器 CentOS7.6 +Tomcat/apache 搭建个人web项目,并允许外网通过80或8080端口访问【超详细】
  4. python语言案例教程单元测试答案_Python单元测试框架(附例子)
  5. 事件驱动模型基本结构
  6. pymy 监控mysql_用Python对MySQL同步状态进行监控_MySQL
  7. Android 使用svg构造交互式中国地图
  8. 总结matlab2018a安装后遇到matlab license manager error -8
  9. 操作系统三十年的兴衰史
  10. 【转】手机软件商店(wiki)