一种基于fMRI与DTI融合的穹隆白质分割的方法专利检索-弥散张量成像磁共振成像分子成像医学影像学诊断设备和程序专利检索查询-专利查询网

一种基于fMRI与DTI融合的穹隆白质分割的方法

阅读：970发布：2020-07-05

专利汇可以提供一种基于fMRI与DTI融合的穹隆白质分割的方法专利检索，专利查询，专利分析的服务。并且一种基于fMRI与DTI融合的穹隆白质分割的方法，包括对fMRI和DTI磁共振影像数据进行预处理、基于静息态功能连接对海马脑区进行区域分割、ROI选取、基于海马子区域和丘脑间白质纤维束示踪技术的穹隆白质分割四个步骤，本发明的创新之处在于提出了通过fMRI功能连接对海马脑区进行分割并以分割子区域作为种子点，以丘脑作为目标脑区进行白质纤维示踪对穹隆白质进行分割的新方法；本发明结合了人脑功能与结构的信息对穹隆白质进行了分割，思想新颖，过程简单，具有重要的科研和临床价值。，下面是一种基于fMRI与DTI融合的穹隆白质分割的方法专利的具体信息内容。

权利要求

1.一种基于fMRI与DTI融合的穹隆白质分割的方法，其特征在于，具体步骤如下：
(1)、对fMRI和DTI磁共振影像数据进行预处理，预处理包括功能磁共振成像和弥散张量成像两个方面，fMRI预处理分别有头颅分割、时间校正、头动校正、空间平滑、时域带通滤波、去线性漂移和冗余去除；DTI预处理包括涡流变形校正、去头皮操作、计算每一个体素的弥散张量和计算FA，MD弥散参数步骤；
(2)、基于静息态功能连接对海马脑区进行区域分割，通过计算海马脑区每个体素的全脑连接模式，求得连接模式的相似度矩阵，运用谱聚类对全脑连接模式相似度矩阵进行聚类分割，实现海马的子区域分割；通过计算海马区域中每个体素的全脑连接模式，进而得到每两个体素之间的全脑连接模式相似度矩阵，以此矩阵进行子区划分；具体做法是提取海马每个体素的时间序列信号，对其他88个脑区的静息态信号，采用脑区内所有体素信号平均的方式提取；其次，通过计算每个海马内部的体素信号与其他88个脑区的信号之间的相关系数，得到海马体素的全脑连接矩阵M，其中元素(i，j)表示第i个海马体素到第j个脑区的功能连接强度；该矩阵的每一行描述了一个海马体素在静息态下和全脑其他脑区之间的连接模式，值越接近1，相关性越大；接着，计算每两行连接模式之间的相关系数就可以度量两个体素在全脑功能连接模式上的相似度，记为相似度矩阵N，其中每个元素表示海马第i个体素到第j个体素在全脑功能连接模式上的相似度；
(3)、ROI选取，采用AAL模板将皮层和皮层下的灰质划分成90个脑区，左右丘脑分别对应为77号脑区和78号脑区，将左右丘脑配准到MNI空间作为目标ROI，将步骤(2)中分割得到的海马子区域作为种子ROI；
(4)、基于海马子区域和丘脑间白质纤维束示踪技术的穹隆白质分割，在MNI空间以左右丘脑作为目标ROI，以海马子区域作为种子ROI，采用白质纤维示踪技术得到海马各子脑区间与丘脑的白质连接，从而实现穹隆白质的分割。
2.根据权利要求1所述的一种基于fMRI与DTI融合的穹隆白质分割的方法，其特征在于，步骤(2)所述的基于静息态功能连接对海马脑区进行区域分割，为了准确定位海马脑区，首先将静息态磁共振图像配准至MNI标准空间；为提高配准精度，采用两步配准：首先将功能像配准到结构像，其次将结构像配准到标准空间，从而利用所得到的变换矩阵将功能像配准到标准空间并将体素重采样为3mm×3mm×3mm；配准模板采用MNI模板；
采用AAL模板将皮层和皮层下的灰质划分成90个脑区，左右海马分别为37和38号脑区，将海马ROI模板以及静息态fMRI数据都配准到MNI空间；在此相似度矩阵的基础上，采用谱聚类算法对海马的体素进行分类。
3.根据权利要求1所述的一种基于fMRI与DTI融合的穹隆白质分割的方法，其特征在于，
所述的步骤(4)具体原理是通过施加多个方向的梯度场，通过计算得到扩散张量矩阵，通过该矩阵的3个特征向量得到包括各向异性参数FA和平均弥散系数的衡量人体内水分子扩散情况的标量，而主特征向量则代表着该体素内水分子主要的扩散方向，首先从白质纤维上选定的某点开始进行白质追踪，确定该点的扩散张量方向，然后沿着该方向前进一定步长，以得到的新体素作为开始点重复进行白质方向确定，直到满足白质追踪停止的条件，最后将所有点连接起来即得到所追踪的白质曲线，用公式表示为：
其中xt+1代表白质纤维追踪下一个体素的坐标位置，而xt表示当前体素的坐标点，α代表每次迭代的步长，则表示当前体素内水分子扩散方向；
白质路径示踪的终止条件设定为角度大于35°或者各向异性参数FA小于0.15，质示踪的纤维束长度如果小于20mm则自动舍弃；最终得到的白质纤维示踪结果，从而实现穹隆白质的分割。

说明书全文

一种基于fMRI与DTI融合的穹隆白质分割的方法

技术领域

[0001] 本发明属于图像处理技术领域，具体涉及一种基于fMRI与DTI融合的穹隆白质分割的方法，尤其涉及通过fMRI功能连接对海马脑区进行分割并以分割子区域作为种子点，以丘脑作为目标脑区进行白质纤维示踪对穹隆白质进行分割的方法。

背景技术

[0002] fMRI是研究脑活动、脑功能的主要的无创方法之一，具有毫米级的空间分辨率。fMRI方法的发展对脑认知功能的研究具有突破性的进展，它已成为神经科学探索人类大脑神经机制的重要工具。DTI利用在不均质组织中水分子扩散表现出各向异性的特性，通过施加多个方向的弥散敏感梯度场可以测得体素内部水分子在每个方向上的弥散强度，计算相关弥散指标，利用所得的多种参数值成像可以描述水分子在组织内的弥散运动。DTI使磁共振(Magnetic Resonance Imaging,MRI)对人体的研究深入到了更微观的水平，它是目前在活体上进行水分子弥散测量和成像的唯一方法。它对于研究正常人以及病理状态下脑白质的完整性和方向性具有很大优越性。fMRI主要反映脑皮质情况，而对于脑白质相关情况则不能显示，而DTI也有其局限性即不能提供皮层功能情况信息，在这方面，fMRI恰好可以弥补其不足。DTI与fMRI结合的多模态成像技术提供了人体神经系统的大量结构和功能信息，在神经发育和衰老、精神疾病研究、神经外科术前准备和术后评价等领域具有广泛的应用。

[0003] 穹隆是海马脑区主要的传入和传出白质纤维，另外也是海马信息传递至丘脑最主要的途径。该白质是边缘系统的关键结构，它涉及很多重要的大脑功能，比如空间记忆、记忆检索、文字记忆，相关解剖研究也证实当沿着穹隆体切断该白质会造成记忆缺失的现象。所以，穹隆的完整性可以作为记忆缺陷障碍的判断，最近的一些研究报道一些具有认知障碍尤其是情景记忆方面的病人穹隆都出现受损情况。作为边缘系统的一部分，穹隆的异常可以比海马体积的萎缩更早地预示认知受损。但目前对穹隆白质的子分支的研究和了解很少涉及，穹隆白质子分支的分割对于记忆相关的精神疾病的预防和靶向治疗有着重要意义，

发明内容

[0004] 为了克服上述现有技术的缺陷，本发明的目的在于提供一种基于fMRI与DTI融合的穹隆白质分割的方法，能够简单、高效地对穹隆白质进行分割，穹隆白质子分支的分割对于记忆相关的精神疾病的预防和靶向治疗有着重要意义。

[0005] 为了达到上述目的，本发明的技术方案为：

[0006] 一种基于fMRI与DTI融合的穹隆白质分割的方法，具体步骤如下：

[0007] (1)、对fMRI和DTI磁共振影像数据进行预处理，预处理包括功能磁共振成像和弥散张量成像两个方面，fMRI预处理分别有头颅分割、时间校正、头动校正、空间平滑、时域带通滤波、去线性漂移和冗余去除；DTI预处理包括涡流变形校正、去头皮操作、计算每一个体素的弥散张量和计算FA，MD弥散参数步骤；

[0008] (2)、基于静息态功能连接对海马脑区进行区域分割，通过计算海马脑区每个体素的全脑连接模式，求得连接模式的相似度矩阵，运用谱聚类对全脑连接模式相似度矩阵进行聚类分割，实现海马的子区域分割；

[0009] (3)、ROI选取，采用AAL模板将皮层和皮层下的灰质划分成90个脑区，左右丘脑分别对应为77号脑区和78号脑区，将左右丘脑配准到MNI空间作为目标ROI，将步骤2中分割得到的海马子区域作为种子ROI；

[0010] (4)、基于海马子区域和丘脑间白质纤维束示踪技术的穹隆白质分割，在MNI空间以左右丘脑作为目标ROI，以海马子区域作为种子ROI，采用白质纤维示踪技术得到海马各子脑区间与丘脑的白质连接，从而实现穹隆白质的分割。

[0011] 本发明的创新之处在于提出了通过fMRI功能连接对海马脑区进行分割并以分割子区域作为种子点，以丘脑作为目标脑区进行白质纤维示踪对穹隆白质进行分割的新方法。传统白质分割需要临床医生进行生理解剖，无法满足临床无创有效的诊断和治疗白质子分支出现病变造成的精神疾病问题，弥散张量成像是目前在体白质纤维成像的唯一方法，本发明结合了人脑功能与结构的信息对穹隆白质进行了分割，思想新颖，过程简单，具有重要的科研和临床价值。附图说明

[0012] 图1是静息态功能图像预处理流程图。

[0013] 图2是头颅分割前后对比图。

[0014] 图3是头动校正结果。

[0015] 图4是功能像配准前后对比，左侧为配准前功能像，中间为MNI152标准模板，右侧为配准后功能像，经过配准图像与模板之间的差异减小。

[0016] 图5是基于静息态功能连接的海马脑区分割结果。

[0017] 图6是白质纤维示踪算法流程示意图。

[0018] 图7是穹隆白质分割结果图。

具体实施方式

[0019] 下面结合附图对本发明做详细叙述。

[0020] 本发明是基于fMRI与DTI融合的穹隆白质分割的方法，具体步骤如下：

[0021] (1)、对fMRI和DTI磁共振影像数据进行预处理。预处理包括功能磁共振成像和弥散张量成像两个方面，fMRI预处理分别有头颅分割、时间校正、头动校正、空间平滑、时域带通滤波、去线性漂移和冗余去除。DTI预处理包括涡流变形校正、去头皮操作、计算每一个体素的弥散张量和计算FA，MD等弥散参数等步骤。

[0022] fMRI数据预处理是在Linux系统Ubuntu14.04下基于AFNI和FSL 软件编写Batch文件(数据批处理)完成。预处理流程见图1，主要包含以下几个方面：

[0023] 1)头颅分割

[0024] 采集得到的结构像通常包含头颅信息，需要进行头颅分割，从而消除颅骨部位引入的伪影对后续数据分析的影响。基于FSL软件提供的3drefit、3dresample和fast segment实现头部颅骨和大脑内部组织结构的分割，分割结果如图2所示。

[0025] 2)时间校正

[0026] 血液动力学函数表明血液对刺激的响应有一定的时间延迟，由于一个TR期间采集到全脑图像，这导致每一层的图像并非在同一时刻采集，而是发生在整个扫描时间段内，时间校正就是通过类似于插值的方法对每一层图像进行层时间处理，使得一个TR周期内各层图像近似在同一时刻获取。

[0027] 3)头动校正

[0028] 进行头动校正时一般将被试的大脑看做一个刚体，因此在fMRI实验中被试头部的运动可以近似成一种刚体运动，即只有平移变换与旋转变换的组合。选择单个被试的第一帧图像作为参考图像，通过AFNI的3dvolreg函数使得其余的所有图像与参考图像配准，结果如图3所示。如果头动超过一个体素则去除该被试。

[0029] 4)空间平滑

[0030] 空间平滑采用高斯函数进行高斯平滑，能有效地消弱随机噪声对fMRI 信号的影响，提高数据的信噪比。三维高斯函数是比较常用的空间平滑方法，其半高全宽决定了空间平滑的力度，本文选用半高全宽(Full Width at Half Maximum,FWHM)为6mm的高斯核函数进行数据平滑。

[0031] 5)时域带通滤波

[0032] 静息态fMRI信号是一种低频波动，频率主要集中在0.01～0.1Hz，这种低频波动反映了自发的神经活动。因此采用频率范围为0.01～0.1Hz的带通滤波器去除与呼吸、心跳等有关的生理噪声。

[0033] 6)去线性漂移

[0034] 由于机器的长期工作导致温度升高或者被试的不适应，随着时间的积累会存在线性漂移，需要进行去线性。

[0035] 7)冗余去除

[0036] 为了去除脑脊液(Cerebro-Spinal Fluid,CSF)、白质(White Matter,WM)等冗余信号，需要对结构像进行分割，利用分割得到的信息制作脑脊液与白质模板，然后去除白质、脑脊液、全脑信号以及头动伪迹等冗余信号。

[0037] DTI数据预处理为使用FMRIB线性配准工具(FLIRT)对DTI数据进行涡流变形校正，将弥散图像都配准至第一幅无弥散加权图像上。接着使用FMRIB中提供的工具(BET)进行去头皮操作。最后使用FMRIB的弥散工具箱计算每一个体素的弥散张量，从而得到单个被试的FA和MD标量图。然后将单个被试的FA和MD标量图配准到结构像，其次将结构像配准到标准空间，从而利用所得到的变换矩阵将功能像配准到标准空间。

[0038] (2)、基于静息态功能连接对海马脑区进行区域分割。为了准确定位海马脑区，首先将静息态磁共振图像配准至MNI标准空间；为提高配准精度，采用两步配准：首先将功能像配准到结构像，其次将结构像配准到标准空间，从而利用所得到的变换矩阵将功能像配准到标准空间并将体素重采样为3mm×3mm×3mm；配准模板采用MNI模板，该模板是由加拿大蒙特利尔神经科学研究所(Montreal Neurological Institute,MNI)研发而成；整个配准过程使用FMRIB提供的线性配准工具实现，配准前后图像对比如图4所示。

[0039] 采用AAL模板将皮层和皮层下的灰质划分成90个脑区。左右海马分别为37和38号脑区，将海马ROI模板以及静息态fMRI数据都配准到MNI空间。通过计算海马区域中每个体素的全脑连接模式，进而得到每两个体素之间的全脑连接模式相似度矩阵，以此矩阵进行子区划分。具体做法是提取海马每个体素的时间序列信号，对其他88个脑区的静息态信号，采用脑区内所有体素信号平均的方式提取。其次，通过计算每个海马内部的体素信号与其他88个脑区的信号之间的相关系数，得到海马体素的全脑连接矩阵M，其中元素(i，j)表示第i个海马体素到第j个脑区的功能连接强度。该矩阵的每一行描述了一个海马体素在静息态下和全脑其他脑区之间的连接模式，值越接近1，相关性越大。接着，计算每两行连接模式之间的相关系数就可以度量两个体素在全脑功能连接模式上的相似度，记为相似度矩阵N，其中每个元素表示海马第i个体素到第j个体素在全脑功能连接模式上的相似度。在此相似度矩阵的基础上，本发明采用谱聚类算法对海马的体素进行分类，分割结果如图5所示；

[0040] (3)、ROI选取。采用AAL(automated anatomical labeling)模板将皮层和皮层下的灰质划分成90个脑区，左右丘脑分别对应为77号脑区和78号脑区，将左右丘脑配准到MNI空间作为目标ROI，将步骤2中分割得到的海马子区域作为种子ROI；

[0041] (4)、基于海马子区域和丘脑间白质纤维束示踪技术的穹隆白质分割。在MNI空间以左右丘脑作为目标ROI，以的海马子区域作为种子ROI，采用白质纤维示踪技术得到海马各子脑区间与丘脑的白质连接。具体原理是通过施加多个方向的梯度场，通过计算得到扩散张量矩阵，通过该矩阵的3个特征向量得到包括各向异性参数FA和平均弥散系数的衡量人体内水分子扩散情况的标量，而主特征向量则代表着该体素内水分子主要的扩散方向，该算法主要就是基于每个体素的扩散向量确定白质走向，首先从白质纤维上选定的某点开始进行白质追踪，确定该点的扩散张量方向，然后沿着该方向前进一定步长，以得到的新体素作为开始点重复进行白质方向确定，直到满足白质追踪停止的条件，最后将所有点连接起来即得到所追踪的白质曲线，具体流程如图6所示。用公式表示为：

[0042]

[0043] 其中xt+1代表白质纤维追踪下一个体素的坐标位置，而xt表示当前体素的坐标点，a代表每次迭代的步长，则表示当前体素内水分子扩散方向。

[0044] 白质路径示踪的终止条件设定为角度大于35°或者FA小于0.15，质示踪的纤维束长度如果小于20mm则自动舍弃。最终得到的白质纤维示踪结果，从而实现穹隆白质的分割，具体结果如图7；

[0045] 综上所述，本发明基于fMRI与DTI融合的多模态成像对穹隆白质进行分割有着实际的临床应用价值。提出了通过fMRI功能连接对海马脑区进行分割并以分割子区域作为种子点，以丘脑作为目标脑区进行白质纤维示踪对穹隆白质进行分割的新思路。本发明结合了人脑功能与结构的信息对穹隆白质进行了分割，思想新颖，过程简单，具有重要的科研和临床价值。

标题	发布/更新时间	阅读量
去除磁共振弥散张量成像噪声的方法和系统	2020-05-17	243
基于脑磁图和弥散张量成像的多模态脑功能重建评估方法	2020-05-21	810
一种汉语语言功能区测查用功能磁共振成像数据的处理方法	2020-05-24	156
一种基于弥散张量的神经成像方法、装置和磁共振成像设备	2020-05-26	505
去除磁共振弥散张量成像噪声的方法和系统	2020-05-16	57
一种基于弥散张量成像技术的疾病诊断辅助系统及方法	2020-05-19	657
一种基于弥散张量成像实现经颅磁刺激深部脑区的方法	2020-05-18	146
磁共振弥散张量成像质量控制综合测试体模	2020-05-17	645
弥散张量成像方法及系统	2020-05-15	641
基于多模态磁共振影像的大脑灰质及白质追踪方法及装置	2020-05-27	441

一种基于fMRI与DTI融合的穹隆白质分割的方法

一种基于fMRI与DTI融合的穹隆白质分割的方法

技术领域

背景技术

发明内容

具体实施方式

该功能需要专业版企业版VIP权限，您可以：