基于多模态磁共振影像的大脑灰质及白质追踪方法及装置
阅读：441发布：2020-05-27
专利汇可以提供基于多模态磁共振影像的大脑灰质及白质追踪方法及装置专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种基于多模态磁共振影像的大脑灰质及白质追踪方法及装置。方法包括：获取样本人群的脑影像数据，得到结构磁共振图像和弥散张量图像；对结构磁共振图像、弥散张量图像进行预处理，得到大脑灰质图像、部分 各向异性 系数图像；对大脑灰质图像和部分各向异性系数图像进行融合，实现大脑灰质及白质追踪。本发明采用多变量分析方法对多模态磁共振图像(结构磁共振图像以及弥散张量图像)进行融合，解决了 现有技术 中使用单变量方法对结构磁共振影像或 弥散张量成像 数据这些单一模态磁共振数据进行分析导致一些隐藏的联合信息被忽略，分析不全面的技术问题，从而提供一种数据更丰富、分析更全面、结果更准确的大脑灰质及白质追踪方法和装置。，下面是基于多模态磁共振影像的大脑灰质及白质追踪方法及装置专利的具体信息内容。
1.一种基于多模态磁共振影像的大脑灰质及白质追踪方法，其特征在于，包括以下步骤：
获取样本人群的脑影像数据，得到多模态磁共振图像，所述多模态磁共振图像包括结构磁共振图像和弥散张量图像；
对所述结构磁共振图像进行预处理，得到大脑灰质图像；
对所述弥散张量图像进行预处理，得到部分各向异性系数图像，所述部分各向异性系数图像中包含有大脑白质纤维束信息；
对所述大脑灰质图像和所述部分各向异性系数图像进行融合，实现大脑灰质及白质追踪。
2.根据权利要求1所述的一种基于多模态磁共振影像的大脑灰质及白质追踪方法，其特征在于，所述步骤对所述结构磁共振图像进行预处理，得到大脑灰质图像具体包括：
对所述结构磁共振图像进行分割，得到初步大脑灰质图像、大脑白质图像和脑脊液图像；
对所述初步大脑灰质图像进行标准化；
将标准化后的大脑灰质图像利用半高全宽高斯核进行平滑，得到大脑灰质图像。
3.根据权利要求1所述的一种基于多模态磁共振影像的大脑灰质及白质追踪方法，其特征在于，所述步骤对所述弥散张量图像进行预处理，得到部分各向异性系数图像具体包括：
对所述弥散张量图像进行涡流校正和头动校正；
移除涡流校正和头动校正后的弥散张量图像中大脑以外的组织，提取大脑组织图像；
对所述大脑组织图像进行张量拟合，得到初步部分各向异性系数图像；
对所述初步部分各向异性系数图像进行标准化，得到部分各向异性系数图像。
4.根据权利要求1所述的一种基于多模态磁共振影像的大脑灰质及白质追踪方法，其特征在于，所述步骤获取样本人群的脑影像数据之前还包括步骤：
对被试人群进行头颅磁共振扫描，筛选无器质性病变的被试人群的脑影像数据，作为样本人群的脑影像数据。
5.根据权利要求4所述的一种基于多模态磁共振影像的大脑灰质及白质追踪方法，其特征在于，所述被试人群的头部与扫描设备的间隙具有弹性物质。
6.根据权利要求5所述的一种基于多模态磁共振影像的大脑灰质及白质追踪方法，其特征在于，所述被试人群的耳部佩戴有隔音耳塞。
7.一种基于多模态磁共振影像的大脑灰质及白质追踪装置，其特征在于，包括：
获取模块，用于获取样本人群的脑影像数据，得到多模态磁共振图像，所述多模态磁共振图像包括结构磁共振图像和弥散张量图像；
第一处理模块，用于对所述结构磁共振图像进行预处理，得到大脑灰质图像；
第二处理模块，用于对所述弥散张量图像进行预处理，得到部分各向异性系数图像，所述部分各向异性系数图像中包含有大脑白质纤维束信息；
融合模块，用于对所述大脑灰质图像和所述部分各向异性系数图像进行融合，实现大脑灰质及白质追踪。
8.根据权利要求7所述的一种基于多模态磁共振影像的大脑灰质及白质追踪装置，其特征在于，所述装置还包括：
筛选模块，用于对被试人群进行头颅磁共振扫描，筛选无器质性病变的被试人群的脑影像数据，作为样本人群的脑影像数据。
基于多模态磁共振影像的大脑灰质及白质追踪方法及装置
技术领域
背景技术
[0002] MCI：Mild Cognitive Impairment，轻度认知障碍。
[0003] HC：Healthy Control，健康对照。
[0004] GM：Grey Matter，灰质，是中枢神经系统的重要组成部分，对信息进行深入处理。
[0005] GMV：Grey Matter Volume，灰质体积。
[0006] WM：White Matter，白质，控制着神经元共享的讯号，协调脑区之间的正常运作。
[0008] sMRI：Structural Magnetic Resonance Imaging，结构磁共振成像，可以用来评估每个体素上灰质的局部浓度或体积的变化，从而反映相应解剖结构的变化。
[0010] FA：Fractional Anisotropy，部分各向异性系数。
[0011] CCA：Canonical Correlation Analysis，典型相关分析，一种多变量分析方法。
[0012] ICA：Independent Component Analysis，独立成分分析，一种多变量分析方法。
[0013] mCCA：multimodal Canonical Correlation Analysis，多模态典型相关分析，一种多变量分析方法，可用于多种模态数据进行融合分析。
[0014] jICA：joint Independent Component Analysis，联合独立成分分析，一种多变量分析方法，可用于对多种模态数据进行融合分析。
[0015] mCCA-jICA：multimodal Canonical Correlation Analysis-joint Independent Component Analysis，多模态典型相关分析-联合独立成分分析，一种多变量分析方法，结合了典型相关分析以及独立成分分析两种方法，对多模态数据进行融合分析。
[0016] 轻度认知障碍是介于正常衰老和痴呆之间的一种中间状态，会导致患者认知能力轻微而显著的下降，其中包括记忆和思维能力。结构核磁共振是一种神经影像类型，可以显示大脑的解剖结构，常用于分析大脑灰质结构信息；弥散张量成像则可以反映大脑的纤维束弥散相关信息，常用于研究大脑白质的纤维束。利用结构磁共振影像以及弥散张量成像技术，这些不同模态的数据可以表征轻度认知障碍引起的大脑灰质与白质的形态差异变化。
[0017] 如今，关于轻度认知障碍相关的灰质以及白质研究大多是使用单变量方法来分析单一模态的影像数据，这些传统方法有感兴趣区域、基于体素的形态测量法和基于纤维束示踪的空间统计方法。利用这些方法，可以对结构磁共振影像或弥散张量成像数据单独进行分析，确定灰质体积变化脑区以及大脑白质纤维束的萎缩状况。然而针对单一数据模态的分析方法无法揭示结构磁共振影像以及弥散张量成像联合提供的隐藏信息，忽略了多个模态之间隐藏的联合信息。这些隐藏的联合信息非常重要，它可以揭示灰质与白质的变化是如何相互作用的，从而更全面地展示轻度认知障碍患者大脑的结构异常，确定这一疾病的大脑神经标记物，为临床上这一疾病的诊断评估做出贡献。
发明内容
[0018] 本发明旨在至少在一定程度上解决相关技术中的技术问题之一。为此，本发明的一个目的是提供一种基于多模态磁共振影像的大脑灰质及白质追踪方法及装置，能够对多模态磁共振图像(结构磁共振图像以及弥散张量图像)进行处理、融合，实现大脑灰质及白质追踪。
[0019] 本发明所采用的技术方案是：
[0020] 第一方面，本发明提供一种基于多模态磁共振影像的大脑灰质及白质追踪方法，包括以下步骤：
[0021] 获取样本人群的脑影像数据，得到多模态磁共振图像，所述多模态磁共振图像包括结构磁共振图像和弥散张量图像；
[0022] 对所述结构磁共振图像进行预处理，得到大脑灰质图像；
[0023] 对所述弥散张量图像进行预处理，得到部分各向异性系数图像，所述部分各向异性系数图像中包含有大脑白质纤维束信息；
[0024] 对所述大脑灰质图像和所述部分各向异性系数图像进行融合，实现大脑灰质及白质追踪。
[0025] 进一步地，所述步骤对所述结构磁共振图像进行预处理，得到大脑灰质图像具体包括：
[0026] 对所述结构磁共振图像进行分割，得到初步大脑灰质图像、大脑白质图像和脑脊液图像；
[0027] 对所述初步大脑灰质图像进行标准化；
[0028] 将标准化后的大脑灰质图像利用半高全宽高斯核进行平滑，得到大脑灰质图像。
[0029] 进一步地，所述步骤对所述弥散张量图像进行预处理，得到部分各向异性系数图像具体包括：
[0030] 对所述弥散张量图像进行涡流校正和头动校正；
[0031] 移除涡流校正和头动校正后的弥散张量图像中大脑以外的组织，提取大脑组织图像；
[0032] 对所述大脑组织图像进行张量拟合，得到初步部分各向异性系数图像；
[0033] 对所述初步部分各向异性系数图像进行标准化，得到部分各向异性系数图像。
[0034] 进一步地，所述步骤获取样本人群的脑影像数据之前还包括步骤：
[0035] 对被试人群进行头颅磁共振扫描，筛选无器质性病变的被试人群的脑影像数据，作为样本人群的脑影像数据。
[0036] 进一步地，所述被试人群的头部与扫描设备的间隙具有弹性物质。
[0038] 第二方面，本发明提供一种基于多模态磁共振影像的大脑灰质及白质追踪装置，包括：
[0039] 获取模块，用于获取样本人群的脑影像数据，得到多模态磁共振图像，所述多模态磁共振图像包括结构磁共振图像和弥散张量图像；
[0040] 第一处理模块，用于对所述结构磁共振图像进行预处理，得到大脑灰质图像；
[0041] 第二处理模块，用于对所述弥散张量图像进行预处理，得到部分各向异性系数图像，所述部分各向异性系数图像中包含有大脑白质纤维束信息；
[0042] 融合模块，用于对所述大脑灰质图像和所述部分各向异性系数图像进行融合，实现大脑灰质及白质追踪。
[0043] 进一步地，所述装置还包括：
[0044] 筛选模块，用于对被试人群进行头颅磁共振扫描，筛选无器质性病变的被试人群的脑影像数据，作为样本人群的脑影像数据。
[0045] 本发明的有益效果是：
[0046] 本发明通过采用一种多变量分析方法对多模态磁共振图像(结构磁共振图像以及弥散张量图像)进行融合，解决了现有技术中使用单变量方法对结构磁共振影像或弥散张量成像数据这些单一模态磁共振数据进行分析导致一些隐藏的联合信息被忽略，分析不全面的技术问题，从而提供了一种数据更丰富、分析更全面、结果更准确的大脑灰质及白质追踪方法和装置。附图说明
[0047] 图1是本发明中基于多模态磁共振影像的大脑灰质及白质追踪方法的一实施例的流程示意图；
[0048] 图2是本发明中存在显著组间差异的灰质成分对应的混合系数散点图；
[0049] 图3是本发明中存在显著组间差异的白质成分对应的混合系数散点图。
具体实施方式
[0050] 需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。
[0051] 实施例一
[0052] 本实施例提供了一种基于多模态磁共振影像的大脑灰质及白质追踪方法，如图1所示，该方法包括以下步骤：
[0053] S1.对被试人群进行头颅磁共振扫描，筛选无器质性病变的被试人群的脑影像数据，作为样本人群的脑影像数据。
[0054] S2.获取样本人群的脑影像数据，得到多模态磁共振图像，多模态磁共振图像包括结构磁共振图像和弥散张量图像；
[0055] S3.对结构磁共振图像进行预处理，得到大脑灰质图像；
[0056] S4.对弥散张量图像进行预处理，得到部分各向异性系数图像，部分各向异性系数图像中包含有大脑白质纤维束信息；
[0057] S5.对大脑灰质图像和部分各向异性系数图像进行融合，实现大脑灰质及白质追踪。
[0058] 具体地，
[0059] 步骤S1：采用Magnetom Verio3.0T超导型磁共振成像仪和标准EPI功能头线圈进行脑影像数据采集。采集数据前对被试人群进行常规头颅磁共振平扫，排除有器质性病变的被试人群，采集剩余的被试人群的脑影像数据作为样本人群的脑影像数据。
[0060] 优选地，采集过程中，为了降低被试人群头部移动所产生的影响，被试人群的头部和头部扫描线圈的间隙添加了弹性物质，例如海绵软垫；为了降低噪音对被试人群的影响，用隔音耳塞填充被试人群的耳朵，并告知其闭眼，放松心情，不要进行任何思维活动并保持安静、清醒的状态(且不可睡着)。
[0061] 具体扫描参数如下：(1)结构磁共振图像：采用三维解剖像，磁化准备快速梯度回波序列，矢状位扫描，TR/TE＝1900ms/2.22ms，FOV＝250mm×250mm，层厚＝1mm，分辨率＝256×256×168，翻转角度＝9°，层数＝176；(2)弥散张量成像：弥散敏感梯度沿30个非共线方向(b＝1000s/mm2)分布，另有一次采集无弥散权重(b＝0s/mm2)，TR/TE＝6800ms/93ms，FOV＝240mm×240m，层厚＝3mm，层数＝46。
[0062] 步骤S2：获取步骤S1中采集的样本人群的脑影像数据，得到多个sMRI图像(即结构磁共振图像)和多个DTI图像(即弥散张量图像)。
[0063] 步骤S3：对多个sMRI图像进行预处理，得到灰质图像，具体包括：
[0064] 步骤S31：将每个被试者的结构磁共振图像分割为初步GM图像即初步大脑灰质图像、大脑白质图像以及脑脊液图像；
[0065] 步骤S32：对分割出来的初步GM图像进行标准化；
[0066] 步骤S33：将标准化后的GM图像利用半高全宽高斯核进行平滑，得到GM图像。
[0067] 步骤S4：对DTI图像进行预处理，得到部分各向异性系数图像，具体包括：
[0068] 步骤S41：对弥散张量图像进行涡流校正和头动校正；
[0069] 步骤S42：移除涡流校正和头动校正后的弥散张量图像中大脑以外的组织，提取大脑组织图像；
[0070] 步骤S43：对大脑组织图像进行张量拟合，得到初步FA图像，即初步部分各向异性系数图像；
[0071] 步骤S44：对初步FA图像进行标准化(与MNI空间对齐)，得到FA图像。
[0072] 步骤S5：利用mCCA-jICA对GM图像和FA图像进行融合，实现大脑灰质及白质追踪。
[0073] 具体的流程如下：
[0074] 首先，以纳入分析的有两组人群为例，第一组：MCI患者即轻度认知障碍患者共99例，其中男性28例，女性71例，平均年龄为65.11±6.55岁；第二组：HC即健康对照有70例，其中男性26例，女性44例，平均年龄为64.64±5.76岁。对于采集到的上述两组人群的sMRI图像(即结构磁共振图像)以及DTI图像(即弥散张量图像)，将sMRI这一模态的数据堆叠成一个矩阵X1，矩阵的维度为：被试者数量×体素个数。其中，体素为构成大脑三维图像的最小单位，因此矩阵X1的每一行对应一个被试者的全脑三维图像数据的展开。同理，所有被试者的DTI这一模态的数据被堆叠成为另一个矩阵X2。
[0075] 利用mCCA即多模态典型相关分析法，分别对两个矩阵X1和X2进行源分离，得到有代表性的典型变量B1、B2以及对应的源构成的矩阵C1、C2。其中，典型变量B1、B2间的相关关系可以反映X1以及X2之间的整体相关性，其维度为：被试者数量×成分个数；与典型变量B1、B2对应的C1、C2为全脑源信号构成的矩阵，维度为：成分个数×体素个数。由于mCCA分离出来的源可能仍包含一些真正的源的线性组合，因此需要将C1以及C2级联成为一个矩阵[C1，C2]后，利用jICA即联合独立成分分析法，再度进行源分离，得到系数矩阵D以及最终的源矩阵[S1，S2]。两个模态数据的混合系数矩阵A1以及A2由mCCA得到的典型变量B1、B2以及jICA得到的系数矩阵D计算得到。混合系数矩阵A1以及A2的维度为：被试者数量×成分个数，即对于每个成分而言，每个被试者包含一个混合系数。矩阵X1、X2、A1、A2、S1、S2、D以及典型变量B1、B2之间的关系如下：
[0076] X1＝A1×S1，X2＝A2×S2，A1＝B1×D，A2＝B2×D
[0077] 后续的统计分析基于上述得到的混合系数矩阵A1以及A2进行。对于每个模态的每个成分而言，对两组被试(MCI以及HC)的混合系数进行双样本t检验分析。例如：对于sMRI的第一个成分，每个被试者会有一个混合系数，这些被试者分为MCI患者(99人)及HC(70人)，将这两组被试者的混合系数进行双样本t检验。在依次对sMRI以及DTI两个模态所有成分的混合系数都进行双样本t检验后，将得到的检验值p进行多重比较校正。以每一模态数据源分离得到的成分个数15为例，若某一成分在双样本t检验得到的p小于3.33×10-3(0.05/成分个数，即0.05/15)，则视这一成分为可以通过Bonferroni校正的成分，存在显著组间差异；若p大于3.33×10-3且小于0.05，则视这一成分在未校正的水平上存在组间差异。混合系数存在组间差异的成分称为组间判别成分，认为这一成分对应的脑区即为MCI患者与HC相比的差异脑区。若在同一个序数的成份中，两个模态的混合系数皆存在显著组间差异的话，这一成分称为联合组间判别成分，即两个模态对应这一序数的成分内的变化脑区是相关的，由此可以确定哪些脑区的灰质体积与白质纤维束的变化是相关的。
[0078] 下面，仍以每一模态数据源分离得到的成分个数15为例，在对每个模态的所有成分对应的混合系数逐一进行双样本t检验后，如图2和图3所示，假设得到大脑灰质存在显著组间差异的成分为IC03、IC05、IC07和IC09，得到大脑白质存在显著组间差异的成分为IC03和IC09，这些成分对应的散点图分别如图2和图3所示，图2是存在显著组间差异的GM(灰质)成分对应的混合系数散点图，图3是存在显著组间差异的WM(白质)成分对应的混合系数散点图。
[0079] 由上述对于混合系数的比较检验结果可知，大脑灰质存在显著组间差异的成分为4个，而大脑白质存在显著组间差异的成分为2个。这些成分对应的脑区即为MCI患者与HC健康对照相比的差异脑区，即MCI患者在这些脑区的灰质体积或白质纤维束与HC健康对照有显著变化。
[0080] 一般情况下，轻度认知障碍患者的灰质体积较健康对照相比的差异区域主要有：丘脑、脑岛、中央前回、颞下回、颞中回、颞上回、梭状回、额中回、中央前回、尾状核、海马以及海马旁回。轻度认知障碍患者的白质纤维束较健康对照相比的差异区域主要有：胼胝体辐射线额部、皮质脊髓束、丘脑前束、上纵束以及扣带(海马)。除了可以查找轻度认知障碍患者与健康对照相比的差异区域外，还可以判断联合组间判别成分。判断方法为：结构磁共振图像以及弥散张量图像两种模态的成分(IC)有15个，如果在两种模态内存在显著组间差异的成分序号相同(参照图2和图3中的GM_IC03和WM_IC03，以及GM_IC09以及WM_IC09)，则称这个成分为联合组间判别成分，即IC03以及IC09为联合组间判别成分。联合组间判别成分中对应的脑区的灰质体积与白质纤维束的变化是相关的。
[0081] 通过判断上述这些脑区的灰质体积或白质纤维束的变化，可以帮助医生对轻度认知障碍患者及早进行无创诊断，制定合适的治疗方案，防止疾病的进一步发展以及阿尔茨海默症的发生。
[0082] 实施例二
[0083] 本实施例提供了一种基于多模态磁共振影像的大脑灰质及白质追踪装置，包括：
[0084] 筛选模块，用于对被试人群进行头颅磁共振扫描，筛选无器质性病变的被试人群的脑影像数据，作为样本人群的脑影像数据；
[0085] 获取模块，用于获取样本人群的脑影像数据，得到多模态磁共振图像，多模态磁共振图像包括结构磁共振图像和弥散张量图像；
[0086] 第一处理模块，用于对结构磁共振图像进行预处理，得到大脑灰质图像；
[0087] 第二处理模块，用于对弥散张量图像进行预处理，得到部分各向异性系数图像，部分各向异性系数图像中包含有大脑白质纤维束信息；
[0088] 融合模块，用于对大脑灰质图像和部分各向异性系数图像进行融合，实现大脑灰质及白质追踪。
[0089] 基于多模态磁共振影像的大脑灰质及白质追踪装置的具体工作过程参照实施例一中对基于多模态磁共振影像的大脑灰质及白质追踪方法的描述，此处不再赘述。
[0090] 本发明的基于多模态磁共振影像的大脑灰质及白质追踪方法及装置，同时利用了结构磁共振图像以及弥散张量图像两种模态的数据来研究分析大脑灰质和大脑白质纤维束，以及大脑灰质与大脑白质的变化是如何相互作用的。将该方法运用到轻度认知障碍患者的诊断评估中，可以更全面地展示轻度认知障碍患者大脑的结构异常，从而更容易确定轻度认知障碍患者大脑灰质体积变化区域以及大脑白质纤维束的萎缩区域，即确定这一疾病的大脑神经标记物，为临床上这一疾病的诊断评估做出贡献。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种基于弥散张量的神经成像方法、装置和磁共振成像设备 | 2020-05-26 | 505 |
| 弥散张量成像方法及系统 | 2020-05-12 | 898 |
| 磁共振弥散张量成像方法和系统 | 2020-05-13 | 38 |
| 一种基于弥散张量成像实现经颅磁刺激深部脑区的方法 | 2020-05-18 | 146 |
| 一种用于磁共振弥散张量成像重建算法验证的立体交叉仿生模型及其制备方法 | 2020-05-23 | 546 |
| 一种基于弥散张量成像的抑郁症患者疾病恢复的评估方法 | 2020-05-20 | 896 |
| 一种基于弥散张量成像技术的疾病诊断辅助系统 | 2020-05-20 | 709 |
| 弥散张量成像方法及系统 | 2020-05-11 | 846 |
| 多b值弥散张量成像采样方法 | 2020-05-11 | 232 |
| 确定弥散张量成像感兴趣区的方法 | 2020-05-12 | 530 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈