一种基于多模态磁共振成像的脑区功能定位方法

摘要

本发明公开了一种基于多模态磁共振成像的大脑脑区功能个体化定位的方法，本发明利用多种模态的磁共振成像方式，包括结构磁共振成像、功能磁共振成像以及扩散磁共振成像，综合考察在脑皮层沟回形态、功能任务参与脑区，以及与之相连的白质纤维走向等三个方面，准确对特定脑区进行个体化的功能区定位。本发明可以用于许多临床神经精神科学的研究，包括临床经颅磁刺激的靶区定位、脑卒中(肿瘤)患者的手术计划及预后判断，以及脑神经精神科学研究中感兴趣区域标定，等等。由于本发明针对被试个体综合考察了脑皮层结构、功能任务脑区粗定位以及与之相连的脑白质纤维分布，因此，从根本上克服了传统基于统一结构脑区模板的局限性。

说明书

技术领域

本发明属于医疗器械及医学图像处理和神经影像领域，尤其是一种基 于多模态磁共振成像的脑区功能定位方法。

背景技术

神经科学以及临床医学的神经精神疾病诊疗中，常常会涉及针对某个 功能脑区的定位，然后对该特定脑区进行手术治疗，或者物理器械干预等。 例如经颅磁刺激(Transcranial Magnetic Stimulation，TMS)设备对脑部背外 侧前额叶(Dorsolateral Prefrontal Cortex，DLPFC)的磁场干预，必须事先对 该手术患者的头部背外侧前额叶有准确的定位，这样才能保证刺激的有效 性，同时也避免对其周围其它的功能脑区造成不必要的伤害。目前临床常 用的脑功能区定位方式有两种：第一，选用适度的频率的磁场(一般由低到 高试探，直至食指有无控制意识的抽动)刺激大脑中央沟两侧的顶叶和额叶 边界靠近额叶的地方，在刺激下引起食指小幅度的抽搐对应的磁刺激的位 置为背外侧前额叶；第二，选用一个标准的通用脑皮层结构图像模板，比 如MNI(Montreal Neurological Institute)模板，以此作为该被测试者的结构 图像，然后在通用的背外侧前额叶的地方作为该被测试者的刺激脑区。上 述两种方式中，前者是以比较粗糙的方式对感兴趣区域进行个体化定位， 而后者则完全没有个体化的考虑，直接选用通用的模板代替个体被测试者， 现有技术的这两种方式都没有真正实现个体化的准确定位，从而降低了后 期治疗效果的临床价值，同时也无法定量的评估实际疗效。

发明内容

为了克服上述现有技术所存在的缺陷，本发明提出一种基于多模态磁 共振成像的脑区功能定位方法，本发明方法针对单个个体被试，分别用结 构磁共振成像、功能磁共振成像和扩散磁共振成像，分析测量该个体针对 某特定功能的脑区的沟回形态、功能任务参与脑区，以及与这些脑区相连 的白质纤维走向等三个特征值，共同从形态、功能及解剖结构三个方面对 个体脑区进行个体化功能定位，以消除传统的统一模板定位方式忽视个体 差异的弊端。

本发明所提出的一种基于多模态磁共振成像的脑功能区定位方法包 括包括以下步骤：

步骤S1，对被测试个体头部进行结构磁共振成像，得到该个体头部灰 质皮层沟回折叠的结构特征；

步骤S2，结合已有的脑解剖学模板，根据步骤S1得到的灰质结构图 像初步确定感兴趣脑功能区的位置；

步骤S3，在对该个体进行功能磁共振成像的同时进行一特定功能任务 测试，并基于功能磁共振成像结果和功能任务测试结果结合所述步骤S2 确定的脑功能区的位置进行筛选优化；

步骤S4，对所述个体进行扩散磁共振成像，对于所述步骤S3得到的 每一优化脑功能区，基于采集得到的扩散磁共振成像数据，得到与所述优 化脑功能区相连的脑白质纤维在成像体素尺度水平的宏观走向分布图，即 神经纤维分布图；

步骤S5，基于所述步骤S3得到的筛选优化后的脑功能区和所述步骤 S4中得到的神经纤维分布图，从筛选优化后的脑功能区中挑选出所有有神 经纤维相互连接的脑功能区，这些脑功能区即为对于所述特定功能任务来 说准确的感兴趣脑功能区。

本发明利用多模态磁共振成像共同从形态、功能及解剖结构三个方面 对个体脑区进行个体化功能定位，最终得到感兴趣区域及与该感兴趣区域 有神经纤维连接的脑功能区。因此，本发明在临床神经精神疾病诊断，手 术计划、导航和定位，治疗的多靶区优化设计，疾病康复的预后判断，以 及神经科学研究上都将有奠基性的应用。

附图说明

图1是本发明基于多模态磁共振成像的脑区功能定位方法流程图；

图2是根据本发明一实施例的神经纤维分布图；

图3是对于抑郁症的标准治疗流程图；

图4是根据本发明方法对于抑郁症的治疗流程图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实 施例，并参照附图，对本发明进一步详细说明。

图1是本发明基于多模态磁共振成像的脑区功能定位方法流程图，如 图1所示，本发明所提出的基于多模态磁共振成像进行脑功能区定位的方 法包括以下步骤：

步骤S1，对被测试个体头部进行结构磁共振成像，得到该个体头部灰 质皮层沟回折叠的结构特征，所述结构特征包括沟回褶皱的深浅，褶皱曲 率，皮层厚度，等等。

所述结构磁共振成像为核磁共振T1加权成像。

步骤S2，结合已有的脑解剖学模板，根据步骤S1得到的灰质结构图 像初步确定感兴趣脑功能区的位置；

所述脑解剖学模板，可以同时结合临床解剖学知识和目前广泛采用的 数字化大脑结构模板，包括加拿大McGill大学主导开发的MNI模板和德 国神经科医生科比尼安·布洛德曼提出的Brodmann图谱。基于所述脑解剖 学模板和已经采集得到的灰质皮层结构图像可以确定该个体感兴趣脑功 能区的大概物理位置。

步骤S3，在对该个体进行功能磁共振成像的同时进行一特定功能任务 测试，并基于功能磁共振成像结果和功能任务测试结果结合所述步骤S2 确定的脑功能区的位置进行筛选优化。

在对该个体进行功能磁共振成像的同时进行一特定功能任务测试，是 指针对某脑功能区，让被试个体在进行功能磁共振扫描的同时参与能特定 激活该脑功能区及相关活动脑功能区的功能任务，使得在功能磁共振成像 中，能显著性的获取该被试个体针对上述特定功能任务的高表达信号。

所述特定功能任务的选择是基于现有技术中的荟萃分析方法 (Meta-Analysis)进行的，具体来说，在使用所述荟萃分析方法选择功能任 务时，首先统计分析已有的文献报道中不同的功能任务与相应的脑功能区 之间已有的反应关系，比如针对经颅磁刺激治疗抑郁症所关注的背外侧前 额叶这一脑功能区，统计有哪些可能的功能任务能有效的激活背外侧前额 叶区域及其他相关联的脑功能区；然后利用文献数量投票方法依据脑功能 区激活的特异性、相关性和易操作性选出对于该脑功能区来说表现最为强 烈的功能任务作为所述特定功能任务。

基于功能磁共振成像结果和功能任务测试结果对所述步骤S2确定的 多个脑功能区的位置进行筛选优化进一步为，首先使用现有技术中常用的 激活区检测的分析方法，确定与所述特定功能任务对应的具体的激活脑功 能区；然后以步骤S2中得到的相应感兴趣脑功能区为参考脑区，剔除掉 那些与参考脑区没有交集的激活脑功能区，剩下的激活脑功能区即为筛选 优化的结果。

步骤S4，对所述个体进行扩散磁共振成像，对于所述步骤S3得到的 每一优化脑功能区，基于采集得到的扩散磁共振成像数据，根据现有技术 中的高角度分辨率成像(High Angular Resolution Diffusion Imaging，HARDI) 的扩散磁共振成像计算方法得到与所述优化脑功能区相连的脑白质纤维 在成像体素尺度水平的宏观走向分布图，即神经纤维分布图，从而得到与 所述优化脑功能区有直接神经纤维连接的每一个脑功能区，从而从生物学 上神经元连接的角度进一步优化所述步骤S3得到的优化脑功能区；

根据高角度分辨率成像的扩散磁共振成像计算方法得到神经纤维分 布图的步骤进一步包括：

首先，获得每一单个体素内表征神经纤维走向分布的水分子扩散方向 分布函数(Orientation Distribution Functions，ODF)；

然后，计算体素与体素之间ODF的相关性；

最后，连接相关性超过特定阈值(比如前后两体素的ODF代表的主 方向夹角小于70度)的相邻体素，得到脑白质纤维在成像体素尺度水平 的宏观走向分布图，即神经纤维分布图，如图2所示，图中包含很多根体 素水平的纤维显示。

步骤S5，基于所述步骤S3得到的筛选优化后的脑功能区和所述步骤 S4中得到的神经纤维分布图，从筛选优化后的脑功能区中挑选出所有有神 经纤维相互连接的脑功能区，这些脑功能区即为对于所述特定功能任务来 说准确的感兴趣脑功能区。

该步骤有两个目的：1，步骤S3得到的脑功能区不是所有的体素都有 神经纤维相互连接，可以从这些脑功能区中剔除没有神经纤维相互连接的 体素，从而进一步缩小优化步骤S3得到的脑功能区(筛选优化得到的脑功 能区命名为A1)；2，寻找与脑功能区A1有神经纤维直接相连接的其它的 脑灰质皮层区域(命名为B1)，B1可以作为经颅磁刺激中针对A1靶区的代 偿区域。

通过上述计算，就可以得到准确的感兴趣的脑功能区。

例如在经颅磁刺激治疗情绪障碍的疾病中，如抑郁症，根据已有的临 床神经精神科学的研究，通常会选择针对个体患者的背外侧前额叶，用低 频的磁场刺激患者的左侧背外侧前额叶，高频的磁场刺激右侧的背外侧前 额叶。背外侧前额叶，其位置大约包括额中回的中间三分之一和额下回的 最前侧的一部分，相当于Brodmann模板的BA46区，也包括Brodmann 9 区的一部分。标准的治疗流程是，先通过一些手段知道该疾病患者的背外 侧额叶的位置，然后分别把刺激线圈分别移到左侧和右侧施加低频和高频 的磁场，低频一般在0-1HZ，高频一般大于8HZ，实际的大小以被试患者 的体验为准(无痛苦但同时有较强的刺激反应)，流程如图4所示。但究竟 通过怎样的手段定位该疾病患者个体的背外侧额叶的位置，本发明就是为 了准确解决这个问题，在该示例中，本发明具体的操作流程如图4所示， 即先对该疾病患者进行结构磁共振扫描，得到头部灰质皮层的脑沟回分布 图，结合临床解剖学知识和已有的脑结构图谱，比如上述描述的MNI图 谱或者Brodmann图谱，初步定位背外侧前额叶所在的沟回位置；然后通 过荟萃分析(Meta-Analysis)综合统计已有关于背外侧前额叶的功能磁共振 成像中功能任务设计，找出强激活或强关联背外侧前额叶的功能任务，并 在被试患者进行功能磁共振扫描过程中执行该测试任务，用激活区检测的 方法，进一步喜欢和筛查上述结构磁共振定位的脑区；最后通过扩散磁共 振成像，选用高角度分辨率成像进行脑白质纤维在成像体素水平的建模和 纤维跟踪，结合上述功能磁共振成像得到的脑区，综合得到更优化的目标 脑区，同时选定脑区的密切关联的脑功能区。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行 了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而 已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修 改、等同替换、改进等，均应包含在本发明的保护范围之内。