首页 / 专利库 / 诊断设备和程序 / 磁共振成像 / 弥散张量成像 / 一种基于FMRI及DTI融合的皮质脊髓束纤维追踪方法

一种基于FMRI及DTI融合的皮质脊髓束纤维追踪方法

阅读:741发布:2020-06-19

专利汇可以提供一种基于FMRI及DTI融合的皮质脊髓束纤维追踪方法专利检索,专利查询,专利分析的服务。并且本 发明 公开一种基于FMRI及DTI融合的皮质脊髓束 纤维 追踪方法,包括对fMRI和DTI磁共振图像数据进行预处理、通过任务态功能激活提取最大激活区定义 种子 点来确定重要运动区域、ROI选取、利用多纤维模型的概率追踪法完成纤维追踪,本发明通过选择fMRI的运动激活区作为DTI纤维追踪的种子点,与传统解剖 位置 定义对脑 肿瘤 患者进行CST纤维追踪进行比较,以此来评价CST的不同成分与脑肿瘤边界的空间位置关系,从而为脑肿瘤患者提供有价值的术前信息。,下面是一种基于FMRI及DTI融合的皮质脊髓束纤维追踪方法专利的具体信息内容。

1.一种基于FMRI及DTI融合的皮质脊髓束纤维追踪方法,其特征在于,包括如下步骤:
步骤一、对fMRI和DTI磁共振图像进行预处理,所述预处理过程包括:功能磁共振、弥散张量成像,其中fMRI预处理包括时序校正、头动校正、标准化、空间平滑、去线性漂移、滤波;
DTI预处理包括头动校正、涡流变形校正、去头皮操作、计算每一个体素内的弥散张量和计算FA值;
步骤二、对预处理后的图像通过任务态功能激活、提取最大激活区种子点,所述最大激活区种子点为重要运动区域;
具体过程如下:首先将预处理后的功能图像通过结构像与标准大脑空间进行标准化处理,然后对标准化处理后的图像做广义线性相关模型进行分析,即通过统计学方法选出若干个区域,对所述若干个区域进行任务态功能激活,最后提取出fMRI最大激活区种子点;
步骤三、ROI感兴趣区域的选取:在在DTI的彩色部分各项异性图,即FA图上勾画感兴趣区region of interest,ROI来跟踪corticospinal tract,CST,基于大脑解剖结构将ROI1放置在桥脑延髓区,基于BOLD-fMRI的ROI2放置在BOLD-fMRI激活区的峰值,将ROI1与ROI2两个区域同时作为新的感兴趣区域;
步骤四、利用多纤维模型的概率追踪法完成纤维追踪:以桥脑延髓区的ROI为起点,以fMRI激活区为终点,联合ROI1与ROI2作为新的感兴趣区域,实现皮质脊髓束的追踪。
2.根据权利要求1所述的基于FMRI及DTI融合的皮质脊髓束纤维追踪方法,其特征在于:步骤二中通过统计学方法选出若干个区域过程如下:通过Z统计学的图像阈值选择为Z>
2.3,矫正后的有意义集群设置阈值为P<0.05,再通过任务态功能激活产生选出的3个区域,所述区域包括左侧主要运动皮层区LPMC、右侧主要运动皮层区RPMC、辅助运动皮层SMA,将提取最大激活区定义种子点来确定重要运动区域。
3.根据权利要求1所述的基于FMRI及DTI融合的皮质脊髓束纤维追踪方法,其特征在于:任务态功能激活过程如下:
所有受试者进行双手运动任务激活,所述运动任务采用组设计模式,即静止休息状态及双手运动任务交替进行,当出现每个视觉刺激画面时进行反复的打开手掌和关闭手掌的运动任务。

说明书全文

一种基于FMRI及DTI融合的皮质脊髓束纤维追踪方法

技术领域

[0001] 本发明涉及图像处理技术领域,特别是涉及一种基于FMRI及DTI融合的皮质脊髓束纤维追踪方法。

背景技术

[0002] 血依赖平的功能磁共振成像技术(BOLD-fMRI)是近年来一种常用的无创性术前功能区定位技术,该方法对于功能区的定位被认为与传统的术中经皮层及皮层下电刺激
术得到的功能定位结果较一致,是目前国际上较为流行的一种新的功能区定位技术。然而,与传统的MRI技术相似,该方法无法可视化观察到重要的皮层下白质纤维束的走行情况。皮层下深部脑肿瘤,其重要功能区(如语言及运动功能区)可能尚未被肿瘤组织破坏,然而,负责语言及运动功能区的白质纤维束由于肿瘤的占位效应可能已导致推压、移位甚至破坏。
如果在术前能够明确肿瘤组织与周围重要神经纤维的关系可帮助神经外科医生制定手术
路径并预估手术切除范围。磁共振弥散张量纤维束成像技术(DTT)是目前临床上唯一可用
于活体显示脑内神经白质纤维束走行轨迹的技术。该技术常通过选择白质纤维束通路上的
感兴趣区的种子点及目标区域,三维的观察白质纤维束的走行情况。然而,目前这些感兴趣的选择常常是基于大脑的解剖标志进行人为的勾画,这种勾画的方法在正常人是较为适用
的,但对于脑肿瘤患者而言,肿瘤的占位效应会造成大脑结构的受压移位,因此,基于解剖标志进行的重要功能区的白质纤维束追踪的结果准确性会受到影响。

发明内容

[0003] 本发明的目的是提供一种基于FMRI及DTI融合的皮质脊髓束纤维追踪方法,以解决上述现有技术存在的问题。
[0004] 为实现上述目的,本发明提供了如下方案:本发明提供一种基于FMRI及DTI融合的皮质脊髓束纤维追踪方法,包括如下步骤:
[0005] 步骤一、对fMRI和DTI磁共振图像进行预处理,所述预处理过程包括:功能磁共振、弥散张量成像,其中fMRI预处理包括时序校正、头动校正、标准化、空间平滑、去线性漂移、滤波;DTI预处理包括头动校正、涡流变形校正、去头皮操作、计算每一个体素内的弥散张量和计算FA值;
[0006] 步骤二、对预处理后的图像通过任务态功能激活、提取最大激活区种子点,所述最大激活区种子点为重要运动区域;
[0007] 具体过程如下:首先将预处理后的功能图像通过结构像与标准大脑空间进行标准化处理,然后对标准化处理后的图像做广义线性相关模型进行分析,即通过统计学方法选
出若干个区域,对所述若干个区域进行任务态功能激活,最后提取出fMRI最大激活区种子
点;
[0008] 步骤三、ROI感兴趣区域的选取:在在DTI的彩色部分各项异性图,即FA图上勾画感兴趣区region of interest,ROI来跟踪corticospinal tract,CST,基于大脑解剖结构将ROI1放置在桥脑延髓区,基于BOLD-fMRI的ROI2放置在BOLD-fMRI激活区的峰值,将ROI1与
ROI2两个区域同时作为新的感兴趣区域;
[0009] 步骤四、利用多纤维模型的概率追踪法完成纤维追踪:以桥脑延髓区的ROI为起点,以fMRI激活区为终点,联合ROI1与ROI2作为新的感兴趣区域,实现皮质脊髓束的追踪。
[0010] 优选地,步骤二中通过统计学方法选出若干个区域过程如下:通过Z统计学的图像阈值选择为Z>2.3,矫正后的有意义集群设置阈值为P<0.05,再通过任务态功能激活产生选出的3个区域,所述区域包括左侧主要运动皮层区LPMC、右侧主要运动皮层区RPMC、辅助运动皮层SMA,将提取最大激活区定义种子点来确定重要运动区域。
[0011] 优选地,任务态功能激活过程如下:
[0012] 所有受试者进行双手运动任务激活,所述运动任务采用组设计模式,即静止休息状态及双手运动任务交替进行,当每个视觉刺激画面时进行反复的打开手掌和关闭手掌
的运动任务。
[0013] 本发明公开了以下技术效果:本发明采用多模态MRI成像技术,将每个脑肿瘤患者的功能区的fMRI激活结果作为纤维追踪的种子点进行DTT成像,能够更加完整、精确的观察颅内肿瘤邻近白质纤维束的改变情况。通过这种多模态MRI成像技术可视化的观察大脑重
要功能纤维束的走行,既能帮助患者最大范围的切除肿瘤,又能在有效保护患者皮层功能
区的同时完整的显示重要白质纤维束走行情况,进一步降低由于手术造成的神经功能缺
陷,从而提高患者的术后生活质量
附图说明
[0014] 为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施
例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
[0015] 图1是本发明的流程示意图;
[0016] 图2是在病灶侧大脑半球不同ROI追踪CST结果;
[0017] 图3是受试者进行双手运动任务激活的流程示意图。

具体实施方式

[0018] 下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他
实施例,都属于本发明保护的范围。
[0019] 为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0020] 参照图1-3,本发明公开了一种基于FMRI及DTI融合的皮质脊髓束纤维追踪方法,并具体公开了如下步骤:
[0021] (1)对fMRI和DTI磁共振影像数据进行预处理,所述预处理包括功能磁共振和弥散张量成像两个方面,fMRI预处理分别有时序校正、头动校正、标准化、空间平滑、去线性漂移、滤波;DTI预处理包括头动校正、涡流变形校正、去头皮操作、计算每一个体素内的弥散张量和计算FA值等弥散参数等步骤。
[0022] fMRI数据预处理数据采用FSL软件进行分析。预处理流程见图1,主要包括以下几个方面:
[0023] 1)时序校正
[0024] 三维图像和一维时间两者构成了4D数据的BOLD-fMRI数据。在三维图像中,每层图像在成像时都会产生一个时间差,这个时间差就是TR。由于原始图像的数据点采集会造成
三维图像的各层图像之间产生时间延迟,我们采取图像与一个血流动学函数进行卷积的
方法消除这个延迟,校正时间差,以保证同一个三维数据中的各层图像都处于一个时间点。
[0025] 2)头动校正
[0026] 由于BOLD-MRI扫描时间较长,并且噪声较大,空间密闭,患者在扫描过程中由于烦躁恐惧等心理,可能不可避免地会出现轻微头动。当头动较大时,将会造成图像点的错位,影响处理结果。头动校正采用的方法是以第一个TR的图像为参考,后续的各时间点图像通过刚性旋转到与第一幅图像保持一致,再用内插算法对各个像素点进行重采样
[0027] 3)标准化
[0028] 进行核磁共振实验的被试之间存在着大脑大小和形状的差异,因此在同一个实验中分析不同被试者时,要将不同被试者的大脑统一到一个脑模板中,上述过程为归一化处
理。在BOLD-fRMI的分析中,此处将EPI选定为标准模板。归一化的具体过程为:先用由最小二乘法得到的冠状面和选定的模板的冠状面进行匹配,接着通过直线回归分析自动检测前
连合和后连合(AC-PC)直线,而连接处的细微差别则采用非线性的方法消除。
[0029] 4)空间平滑
[0030] 空间平滑采用高斯函数进行高斯平滑,能有效地消弱随机噪声对fMRI信号的影响,提高数据的信噪比。本文选用半高全宽(Full Width at Half Maximum,FWHM)为6mm的
高斯核函数进行数据平滑。
[0031] 5)去线性漂移
[0032] 由于BOLD效应、大脑中的血流、代谢和电磁信号的综合效应,高时间分辨率和热噪声的噪音始终存在。而另一方面,由于磁共振工作过程中的升温的跟随效应,信号中会存在着随时间偏离基线的趋势。通过去线性漂移可纠正受试者的数据对生理基线的偏离。
[0033] 6)滤波
[0034] 在静息态BOLD-fMRI信号中,0.01~0.08Hz的信号被认为是具有生理意义的,能反映自发的神经元活动。所以静息态BOLD-fMRI数据在分析之前首先需要带通滤波,提取有效信号。
[0035] 对DTI数据进行预处理:对所有的DTI数据均利用线性配准工具(FLIRT)进行涡流变形校正,将弥散图像都配准到第一幅无弥散加权图像上。接着使用FSL软件包中的BET工
具进行去头皮操作,使用掩膜覆盖所有的弥散加权图像。最后使用FMRIB的弥散工具箱计算每一个体素内的弥散张量,从而计算出FA值。
[0036] (2)通过任务态功能激活提取最大激活区定义种子点来确定重要运动区域。
[0037] 首先功能图像通过结构像与MNI152标准大脑空间进行标准化。
[0038] 参照图3,对所有受试者进行双手运动任务激活,所述运动任务采用组块设计模式(如ABAB模式),即30s静止休息状态(A)及30s双手运动任务(B)交替进行,整个任务过程中
共有6个静息状态,5个任务状态组成;所有的病人和健康对照被引导面对每个视觉刺激画
面时进行反复的打开手掌和关闭手掌的运动任务;所有被试者在扫描前都经过专的运动
任务培训,以确保完全理解该任务并能够认真配合完成扫描。
[0039] 再采用FSL FEAT做广义线性相关模型分析:Z统计学的图像阈值选择为Z>2.3,矫正后的有意义集群设置阈值为P<0.05。选定的3个区域包括左侧主要运动皮层区(LPMC)、右侧主要运动皮层区(RPMC)及辅助运动皮层区(SMA)是通过任务态功能激活产生的,用于提
取最大激活区定义种子点来确定重要运动区域。
[0040] (3)ROI种子点的选取。将F的阈值设置为0.2,在DTI的彩色部分各项异性图,即FA图上勾画感兴趣区region of interest,ROI来跟踪corticospinal tract,CST,在横轴位
彩色FA图像上放置ROI来跟踪CST,基于大脑解剖结构将ROI1放置在桥脑延髓区,基于BOLD-fMRI的ROI2放置在BOLD-fMRI激活区的峰值,将ROI1与ROI2两个区域同时作为新的感兴趣
区域;
[0041] (4)最后利用多纤维模型的概率追踪法完成纤维追踪:以桥脑延髓区的ROI为起点,以fMRI激活区为终点,基于新的感兴趣区域实现皮质脊髓束的追踪。
[0042] 图2显示在病灶侧大脑半球不同ROI追踪CST结果(绿色:CST的感兴趣区置于脑桥延髓层面;黄色:CST感兴趣区置于M1功能激活区;紫色:CST感兴趣置于桥脑延髓层面及M1功能激活区);从此结果中可以详尽的观察到整个CST走行情况与脑肿瘤的空间位置关系。
[0043] 本发明基于fMRI及DTI融合的大脑肿瘤术前对皮质脊髓束纤维追踪有着实际的临床应用价值。提出了采用多模态MRI进行扫描,针对患者的肿瘤位置进行个体化任务fMRI刺激,并有效的联合DTT技术进行重要功能区白质纤维束追踪,使重要功能区白质纤维束更加精准的可视化显示出来。本发明思路新颖,过程简单,有重要的临床和科研价值。
[0044] 在本发明的描述中,需要理解的是,术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
[0045] 以上所述的实施例仅是对本发明的优选方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出
的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈