CN111340821A - 一种基于模块连接的大脑结构网络的偏测性检测方法 - Google Patents

一种基于模块连接的大脑结构网络的偏测性检测方法 Download PDF

Info

Publication number
CN111340821A
CN111340821A CN202010105253.0A CN202010105253A CN111340821A CN 111340821 A CN111340821 A CN 111340821A CN 202010105253 A CN202010105253 A CN 202010105253A CN 111340821 A CN111340821 A CN 111340821A
Authority
CN
China
Prior art keywords
brain
fiber bundle
module
network
matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010105253.0A
Other languages
English (en)
Other versions
CN111340821B (zh
Inventor
李丹丹
王彬
崔晓红
相洁
曹锐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Taiyuan University of Technology
Original Assignee
Taiyuan University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CN202010105253.0A priority Critical patent/CN111340821B/zh
Publication of CN111340821A publication Critical patent/CN111340821A/zh
Application granted granted Critical
Publication of CN111340821B publication Critical patent/CN111340821B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/11—Region-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/136—Segmentation; Edge detection involving thresholding
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30016—Brain

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明公开一种基于模块连接的大脑结构网络的偏测性检测方法，包括以下步骤：对弥散张量成像进行预处理，根据选定的标准化脑图谱，对预处理后的弥散张量成像进行区域分割；采用确定性纤维束追踪算法，基于追踪的结束条件对全脑的纤维进行重建，并计算两两脑区间的纤维束数量与部分各向异性指数、以及每个脑区的表面积，得到脑区间的纤维束数量矩阵FN，部分各向异性指数矩阵FA和脑区表面积矩阵Surface；与传统的偏侧性检测方法相比，本发明所述的一种基于模块连接的大脑结构网络的偏侧性检测方法，将网络模块性与偏侧性进行融合，能够更有效、精确地挖掘大脑局部网络指标的偏侧性，从而为探究大脑工作机制提供一定的帮助。

Description

一种基于模块连接的大脑结构网络的偏测性检测方法
技术领域
本发明属于脑结构成像领域以及大脑网络拓扑结构分析技术领域，尤其涉及一种基于模块连接的大脑结构网络的偏测性检测方法。
背景技术
偏侧性是指大脑左右半球在结构和功能上存在的不对称性，其源于胚胎时期且随着年龄、环境和经验的积累而发生变化，是人脑发育的一个基本特征。大脑结构偏侧性主要包括大脑白质体积、白质整合性等。大脑功能偏侧性是指时空处理、语言、运动以及认知控制功能在左右半球存在的差异。许多神经学研究揭示正常人存在白质偏侧性的区域往往伴随着大脑功能的偏侧性。如研究表明正常人在枕顶区和眶额区的白质呈右偏侧性，这与该区域对应的视空以及注意功能的右偏侧性一致。因此，研究左右半球白质的差异性对我们探究大脑认知功能的加工机制非常重要。
基于神经影像技术和复杂网络理论的脑网络方法，为我们研究脑区间的连接模式提供了有效途径，已被成功用于大脑偏侧性研究中。结合弥散张量成像(diffusion tensorimaging,DTI)技术和图理论，通过研究左右两个半球内部的白质网络，一些研究已经证实大脑两个半球网络之间存在显著的拓扑不对称性。例如，健康成年人的右半球网络在全局效率显著优于左半球网络，表明右半球在视空注意功能显著优于左半球。
近来，大量研究表明脑网络具备模块化结构的拓扑性质。模块化结构表现了网络的聚集程度，是复杂网络的基本属性之一。每一个模块都由网络中的部分节点组成，人脑可以划分为若干个模块(主要包括运动模块、视觉模块、注意模块、默认模块和边缘系统)，每个模块对应不同的大脑功能，其解释了大脑结构和功能之间的潜在关系，为探究大脑工作机制提供了新的视角。传统的关于大脑结构网络的偏侧性研究主要是探究大脑左右半球网络在拓扑结构上存在的差异性，忽视了网络模块内以及模块间拓扑性质的偏侧性。因此，基于模块连接的偏侧性研究，更有助于从局部脑区连接的角度，探究大脑拓扑结构的偏侧性，进一步探究大脑工作机制。
发明内容
针对现有技术的不足，本发明提出了采用一种基于社区划分算法Louvain对大脑结构网络进行模块划分，基于划分的模块，利用偏侧性指数对模块内部与模块间结构连接的不对称性进行检测。
一种基于模块连接的大脑结构网络的偏测性检测方法，采用如下步骤实现的：
步骤S1：对弥散张量成像进行预处理，然后根据选定的标准化脑图谱，对预处理后的弥散张量成像进行区域分割；
步骤S2：采用确定性纤维束追踪算法，基于追踪的结束条件对全脑的纤维进行重建，并计算两两脑区间的纤维束数量与部分各向异性指数以及每个脑区的表面积，由此得到脑区间的纤维束数量矩阵FN，部分各向异性指数矩阵FA和脑区表面积矩阵Surface；
步骤S3：将矩阵FN与FA进行耦合，并基于Surface矩阵进行归一化，构建大脑结构网络矩阵SC；
步骤S4：将大脑结构网络矩阵SC划分为左半球网络SCL和右半球网络SCR；
步骤S5：基于社区划分算法Louvain分别对两个网络SCL和SCR进行模块划分；
步骤S6：基于图理论，从全局效率、局部效率、连接密度、连接强度的角度分析各模块内部、以及模块间的结构连接；
步骤S7：通过计算各模块内部与两两模块间连接指标的偏侧性指数，得到基于模块连接的大脑结构网络的偏侧性分析结果。
在步骤S1中，首先利用FSL中的BET工具箱对DTI原始图像进行去头骨操作；然后利用FSL中的FLIRT工具进行头动矫正，以及利用FSL中的eddy_correct工具进行涡流矫正；最后利用PANDA软件对DTI图像进行张量重建；根据选定的标准化脑图谱BN模板，对预处理后的DTI图像进行分割，将大脑分为246个脑区；
其中，PANDA软件是由北京师范大学研发的用于分析大脑弥散张量成像的工具箱；
在步骤S2中，纤维束追踪的结束条件具体包括：1)在纤维束追踪过程中，若某条纤维束追踪到达某一体素时，该条纤维束的部分各向异性指数小于0.1，则该条纤维束的追踪终止；2)在纤维束追踪过程中，若某条纤维束追踪到达某一体素时，该条纤维束位于大脑皮层的边界，则该条纤维束的追踪终止；3)在纤维束追踪过程中，若某条纤维束追踪到达某一体素时，该条纤维束的偏转角度大于35°，则该条纤维束的追踪终止。基于纤维束追踪的结果，计算两两脑区间的纤维束数量、部分各向异性指数以及脑区表面积，由此得到脑区间的纤维束数量矩阵FN，部分各向异性指数矩阵FA和脑区表面积矩阵Surface；
在步骤S3中，构建结构脑网络SC的处理公式如下所示：
Figure BDA0002388332000000031
公式(1)中，SCi,j表示大脑结构网络矩阵SC中第i行第j列的元素；FAi,j表示部分各向异性指数矩阵FA中第i行第j列的元素；FNi,j表示纤维束数量矩阵FN中第i行第j列的元素；Surfacei表示脑区i的表面积；τ表示阈值，取值为3；大脑结构网络矩阵SC的维度为246×246。
在步骤S4中，将大脑结构网络SC划分为左右两个半球网络的处理公式如下；
Figure BDA0002388332000000041
Figure BDA0002388332000000042
其中，
Figure BDA0002388332000000043
公式(2)和(3)中，SCi,j表示大脑结构网络矩阵SC中第i行第j列的元素；SCLiL,jL表示大脑左半球结构网络矩阵SCL中第iL行第jL列的元素，其中；SCRiR,jR表示大脑右半球结构网络矩阵SCR中第iR行第jR列的元素。
在步骤S5中，Louvain算法是基于模块度的社区发现算法，其优化目标是最大化整个社区网络的模块度。基于社区发现算法Louvain对大脑结构网络进行模块划分包括以下步骤：
步骤S51：社区网络的模块度是评估一个社区网络划分好坏的度量方法，它的含义是社区内节点的连边数与随机情况下的边数之差，它的取值范围是(0,1)，其定义如下：
Figure BDA0002388332000000044
公式中，Wi,j表示节点i和j之间的权值，ki表示节点的度，m为网络中的总连接边数，ci为节点所在的社团编号。∑in表示社区c内的所有边的权重之和。∑tot表示与社区∑tot内的节点相连的边的权重之和。
步骤S52:社区划分算法Louvain划分模块的方法如下：①初始时将每个顶点当作一个社区，社区个数与顶点个数相同；②依次将每个顶点与之相邻顶点合并在一起，计算它们的模块度增益ΔQ是否大于0，如果大于0，就将该结点放入该相邻结点所在社区；③迭代第二步，直至算法稳定，即所有顶点所属社区不再变化；④将各个社区所有节点压缩成为一个结点，社区内点的权重转化为新结点环的权重，社区间权重转化为新结点边的权重；⑤重复步骤1-3，直至算法稳定。
Figure BDA0002388332000000045
在步骤S6中，利用图理论分析各模块内部的结构连接，包括全局效率、局部效率、连接密度、连接强度；各指标的具体计算公式如下：
全局效率(Eg):
Figure BDA0002388332000000051
局部效率(Eloc):
Figure BDA0002388332000000052
连接密度(D):
Figure BDA0002388332000000053
连接强度(S):
Figure BDA0002388332000000054
公式(6)-(9)中，H表示半球网络；c表示模块网络；
Figure BDA0002388332000000055
表示半球H内模块c的节点个数；
Figure BDA0002388332000000056
表示半球H内模块c的结构连接个数；
Figure BDA0002388332000000057
表示半球H内模块c的连接密度；
Figure BDA0002388332000000058
表示半球H内模块c的所有连接权重之和；
Figure BDA0002388332000000059
表示半球H内模块c的连接强度；dij表示节点i和节点j间的最短路径长度；
Figure BDA00023883320000000510
表示半球H内模块c的全局效率，
Figure BDA00023883320000000511
表示半球H内模块c的局部效率。
在步骤S7中:不对称性指数的计算方法如下：
Figure BDA00023883320000000512
公式(10)中，GR表示右半球结构网络的图论连接指标，GL表示左半球结构网络的图论连接指标，
Figure BDA00023883320000000513
表示右半球内模块c的图论连接指标，
Figure BDA00023883320000000514
表示左半球内模块c的图论连接指标。
与现有技术相比，本发明的有益效果是：与传统的偏侧性检测方法相比，本发明所述的一种基于模块连接的大脑结构网络的偏侧性检测方法，将网络模块性与偏侧性进行融合，能够更有效、精确地挖掘大脑局部网络指标的偏侧性，从而为探究大脑工作机制提供一定的帮助。
附图说明
图1为本发明方法流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。
实施例
(1)实验数据
实验选取了52名正常人(Normal control,NC)作为实验对象。
(2)实验过程
请参阅图1，按照图1所示的流程，采用如下步骤对52名被试的大脑结构网络的偏侧性进行分析：
步骤S1：对弥散张量成像进行预处理，然后根据选定的标准化脑图谱，对预处理后的弥散张量成像进行区域分割；
在步骤S1中，采用脑功能磁共振成像软件(FMRIB Software Library,FSL)进行，预处理的步骤具体包括：磁化系数修正、涡流失真修正、头动矫正；标准化脑图谱采用脑成像国际联盟(the International Consortium of Brain Mapping，ICBM)图谱；对预处理后的弥散张量图像进行分割，将大脑分为246个脑区；
步骤S2：采用确定性纤维束追踪算法，基于追踪的结束条件对全脑的纤维进行重建，并计算两两脑区间的纤维束数量、部分各向异性指数以及每个脑区的表面积，由此得到脑区间的纤维束数量矩阵FN、部分各向异性指数矩阵FA以及表面积矩阵Surface；
在步骤S2中，确定性纤维束追踪算法采用以下四种算法的任一种：纤维联络连续追踪(Fiber Assignment by Continuous Tracking,FACT)、2nd order runge-kutta(2ndorderRK)、Tensoline与Interplated Streamline；纤维束追踪的结束条件具体包括：1)在纤维束追踪过程中，若某条纤维束追踪到达某一体素时，该条纤维束的部分各向异性指数小于0.1，则该条纤维束的追踪终止；2)在纤维束追踪过程中，若某条纤维束追踪到达某一体素时，该条纤维束位于大脑皮层的边界，则该条纤维束的追踪终止；3)在纤维束追踪过程中，若某条纤维束追踪到达某一体素时，该条纤维束的偏转角度大于35°，则该条纤维束的追踪终止。
步骤S3：将矩阵FN与FA进行耦合，并基于Surface矩阵进行归一化，构建大脑结构网络矩阵SC；
在步骤S3中，利用公式(1)计算两个脑区之间的连接：
Figure BDA0002388332000000071
公式(1)中，SCi,j表示大脑结构网络矩阵SC中第i行第j列的元素；FAi,j表示部分各向异性指数矩阵FA中第i行第j列的元素；FNi,j表示纤维束数量矩阵FN中第i行第j列的元素；Surfacei表示脑区i的表面积；τ表示阈值，取值为3；大脑结构网络矩阵SC的维度为246×246。
步骤S4：将大脑结构网络SC划分为左半球网络SCL和右半球网络SCR；
在步骤S4中，利用公式(2)构建左半球网络SCL，利用公式(3)构建右半球网络SCR：
Figure BDA0002388332000000072
Figure BDA0002388332000000073
其中，
Figure BDA0002388332000000074
公式(2)和(3)中，SCi,j表示大脑结构网络矩阵SC中第i行第j列的元素；SCLiL,jL表示大脑左半球结构网络矩阵SCL中第iL行第jL列的元素，其中；SCRiR,jR表示大脑右半球结构网络矩阵SCR中第iR行第jR列的元素。
步骤S5：基于社区划分算法Louvain分别对两个网络SCL和SCR进行模块划分；
在步骤S5中，Louvain算法是基于模块度的社区发现算法，其优化目标是最大化整个社区网络的模块度。基于社区发现算法Louvain大脑结构网络进行模块划分包括以下步骤：
步骤S51：社区网络的模块度是评估一个社区网络划分好坏的度量方法，它的含义是社区内节点的连边数与随机情况下的边数之差，它的取值范围是(0,1)，其定义如下：
Figure BDA0002388332000000081
公式中，Wi,j表示节点i和j之间的权值，ki表示节点的度，m为网络中的总连接边数，ci为节点所在的社团编号。∑in表示社区c内的所有边的权重之和。∑tot表示与社区∑tot内的节点相连的边的权重之和。
步骤S52:社区划分算法Louvain划分模块的方法如下：①初始时将每个顶点当作一个社区，社区个数与顶点个数相同；②依次将每个顶点与之相邻顶点合并在一起，计算它们的模块度增益ΔQ是否大于0，如果大于0，就将该结点放入该相邻结点所在社区；③迭代第二步，直至算法稳定，即所有顶点所属社区不再变化；④将各个社区所有节点压缩成为一个结点，社区内点的权重转化为新结点环的权重，社区间权重转化为新结点边的权重；⑤重复步骤1-3，直至算法稳定。
Figure BDA0002388332000000082
步骤S6：利用图理论分析各模块内部的结构连接，包括全局效率、局部效率、连接密度、连接强度；
在步骤S6中，利用图理论分析各模块内部的结构连接，包括全局效率、局部效率、连接密度、连接强度；各指标的具体计算公式如下：
全局效率(Eg):
Figure BDA0002388332000000083
局部效率(Eloc):
Figure BDA0002388332000000084
连接密度(D):
Figure BDA0002388332000000085
连接强度(S):
Figure BDA0002388332000000086
公式(6)-(9)中，H表示半球网络；c表示模块网络；
Figure BDA0002388332000000087
表示半球H内模块c的节点个数；
Figure BDA0002388332000000091
表示半球H内模块c的结构连接个数；
Figure BDA0002388332000000092
表示半球H内模块c的连接密度；
Figure BDA0002388332000000093
表示半球H内模块c的所有连接权重之和；
Figure BDA0002388332000000094
表示半球H内模块c的连接强度；dij表示节点i和节点j间的最短路径长度；
Figure BDA0002388332000000095
表示半球H内模块c的全局效率，
Figure BDA0002388332000000096
表示半球H内模块c的局部效率。
步骤S7：通过计算各模块内部与两两模块间连接指标的偏侧性指数，得到基于模块连接的大脑结构网络的偏侧性分析结果。
在步骤S7中:利用公式(10)计算偏侧性指数：
Figure BDA0002388332000000097
公式(10)中，GR表示右半球结构网络的图论连接指标，GL表示左半球结构网络的图论连接指标，
Figure BDA0002388332000000098
表示右半球内模块c的图论连接指标，
Figure BDA0002388332000000099
表示左半球内模块c的图论连接指标。
(3)实验结果
表1半球网络内部与模块内部结构连接存在显著偏侧性的指标
Figure BDA00023883320000000910
表2模块间的结构连接存在显著偏侧性指标
Figure BDA00023883320000000911
Figure BDA0002388332000000101
与传统的偏侧性检测方法相比，本发明所述的一种基于模块连接的大脑结构网络的偏侧性检测方法，将网络模块性与偏侧性进行融合，能够更有效、精确地挖掘大脑局部网络指标的偏侧性，从而为探究大脑工作机制提供一定的帮助。
以上仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种基于模块连接的大脑结构网络的偏测性检测方法，其特征在于，包括以下步骤：
步骤S1：对弥散张量成像进行预处理，根据选定的标准化脑图谱，对预处理后的弥散张量成像进行区域分割；
步骤S2：采用确定性纤维束追踪算法，基于追踪的结束条件对全脑的纤维进行重建，并计算两两脑区间的纤维束数量与部分各向异性指数、以及每个脑区的表面积，得到脑区间的纤维束数量矩阵FN，部分各向异性指数矩阵FA和脑区表面积矩阵Surface；
步骤S3：将矩阵FN与FA进行耦合，并基于Surface矩阵进行归一化，构建大脑结构网络矩阵SC；
步骤S4：将大脑结构网络矩阵SC划分为左半球网络SCL和右半球网络SCR；
步骤S5：基于社区划分算法Louvain分别对两个网络SCL和SCR进行模块划分；
步骤S6：基于图理论，从全局效率、局部效率、连接密度、连接强度的角度分析各模块内部、以及模块间的结构连接；
步骤S7：通过计算各模块内部与两两模块间连接指标的偏侧性指数，得到基于模块连接的大脑结构网络的偏侧性分析结果。
2.根据权利要求1所述的检测方法，其特征在于，在步骤S1中，具体包括：首先利用FSL中的BET工具箱对弥散张量成像的原始图像进行去头骨操作；然后利用FSL中的FLIRT工具进行头动矫正，以及利用FSL中的eddy_correct工具进行涡流矫正；最后利用PANDA软件对弥散张量成像图像进行张量重建；根据选定的标准化脑图谱，对预处理后的弥散张量成像图像进行分割，将大脑分为246个脑区。
3.根据权利要求2所述的检测方法，其特征在于，在步骤S2中，纤维束追踪的结束条件具体包括：1)在纤维束追踪过程中，若某条纤维束追踪到达某一体素时，该条纤维束的部分各向异性指数小于0.1，则该条纤维束的追踪终止；2)在纤维束追踪过程中，若某条纤维束追踪到达某一体素时，该条纤维束位于大脑皮层的边界，则该条纤维束的追踪终止；3)在纤维束追踪过程中，若某条纤维束追踪到达某一体素时，该条纤维束的偏转角度大于35°，则该条纤维束的追踪终止。
4.根据权利要求3所述的检测方法，其特征在于，在步骤S3中，构建结构脑网络SC的处理公式如下所示：
Figure FDA0002388331990000021
公式(1)中，SCi,j表示大脑结构网络矩阵SC中第i行第j列的元素；FAi,j表示部分各向异性指数矩阵FA中第i行第j列的元素；FNi,j表示纤维束数量矩阵FN中第i行第j列的元素；Surfacei表示脑区i的表面积；τ表示阈值，取值为3；大脑结构网络矩阵SC的维度为246×246。
5.根据权利要求4所述的检测方法，其特征在于，在步骤S4中，将大脑结构网络SC划分为左右两个半球网络的处理公式如下；
Figure FDA0002388331990000022
Figure FDA0002388331990000023
其中，
Figure FDA0002388331990000024
公式(2)和(3)中，SCi,j表示大脑结构网络矩阵SC中第i行第j列的元素；SCLiL,jL表示大脑左半球结构网络矩阵SCL中第iL行第jL列的元素，其中；SCRiR,jR表示大脑右半球结构网络矩阵SCR中第iR行第jR列的元素。
6.根据权利要求5所述的检测方法，其特征在于，在步骤S5中，基于社区发现算法Louvain对大脑结构网络进行模块划分包括以下步骤：
步骤S51：社区网络的模块度是社区内节点的连边数与随机情况下的边数之差，它的取值范围是(0,1)，其定义如下：
Figure FDA0002388331990000031
公式中，Wi,j表示节点i和j之间的权值，ki表示节点的度，m为网络中的总连接边数，ci为节点所在的社团编号，∑in表示社区c内的所有边的权重之和，∑tot表示与社区∑tot内的节点相连的边的权重之和；
步骤S52:社区划分算法Louvain划分模块的方法包括：
①初始时将每个顶点当作一个社区，社区个数与顶点个数相同；
②依次将每个顶点与之相邻顶点合并在一起，计算它们的模块度增益ΔQ是否大于0，如果大于0，就将该结点放入该相邻结点所在社区，计算公式如下：
Figure FDA0002388331990000032
③迭代第二步，直至算法稳定，即所有顶点所属社区不再变化；
④将各个社区所有节点压缩成为一个结点，社区内点的权重转化为新结点环的权重，社区间权重转化为新结点边的权重；
⑤重复步骤①-③，直至算法稳定。
7.根据权利要求6所述的检测方法，其特征在于，在步骤S6中，利用图理论分析各模块内部的结构连接，包括全局效率、局部效率、连接密度、连接强度；各指标的具体计算公式如下：
全局效率(Eg):
Figure FDA0002388331990000033
局部效率(Eloc):
Figure FDA0002388331990000034
连接密度(D):
Figure FDA0002388331990000035
连接强度(S):
Figure FDA0002388331990000036
公式(6)-(9)中，H表示半球网络；c表示模块网络；
Figure FDA0002388331990000041
表示半球H内模块c的节点个数；
Figure FDA0002388331990000042
表示半球H内模块c的结构连接个数；
Figure FDA0002388331990000043
表示半球H内模块c的连接密度；
Figure FDA0002388331990000044
表示半球H内模块c的所有连接权重之和；
Figure FDA0002388331990000045
表示半球H内模块c的连接强度；dij表示节点i和节点j间的最短路径长度；
Figure FDA0002388331990000046
表示半球H内模块c的全局效率，
Figure FDA0002388331990000047
表示半球H内模块c的局部效率。
8.根据权利要求7所述的检测方法，其特征在于，在步骤S7中，偏侧性指数的计算方法如下：
Figure FDA0002388331990000048
公式(10)中，GR表示右半球结构网络的图论连接指标，GL表示左半球结构网络的图论连接指标，
Figure FDA0002388331990000049
表示右半球内模块c的图论连接指标，
Figure FDA00023883319900000410
表示左半球内模块c的图论连接指标。
CN202010105253.0A 2020-02-20 2020-02-20 一种基于模块连接的大脑结构网络的偏测性检测方法 Active CN111340821B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010105253.0A CN111340821B (zh) 2020-02-20 2020-02-20 一种基于模块连接的大脑结构网络的偏测性检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010105253.0A CN111340821B (zh) 2020-02-20 2020-02-20 一种基于模块连接的大脑结构网络的偏测性检测方法

Publications (2)

Publication Number Publication Date
CN111340821A true CN111340821A (zh) 2020-06-26
CN111340821B CN111340821B (zh) 2020-12-08

Family

ID=71183555

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010105253.0A Active CN111340821B (zh) 2020-02-20 2020-02-20 一种基于模块连接的大脑结构网络的偏测性检测方法

Country Status (1)

Country Link
CN (1) CN111340821B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111784206A (zh) * 2020-07-29 2020-10-16 南昌航空大学 采用LeaderRank算法评估社交网络关键节点的方法
CN112741613A (zh) * 2021-01-13 2021-05-04 武汉大学 一种静息态人脑默认网络功能与结构耦合分析方法
CN113283465A (zh) * 2021-04-02 2021-08-20 电子科技大学 一种弥散张量成像数据分析方法及装置

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090129649A1 (en) * 2007-11-20 2009-05-21 Faycal Djeridane Method and system for processing multiple series of biological images obtained from a patient
CN102609946A (zh) * 2012-02-08 2012-07-25 中国科学院自动化研究所 一种基于黎曼流形的脑白质纤维束跟踪的组间处理方法
CN103345749A (zh) * 2013-06-27 2013-10-09 中国科学院自动化研究所 一种基于模态融合的大脑网络功能连接偏侧性检测方法
CN104346530A (zh) * 2014-10-29 2015-02-11 中国科学院深圳先进技术研究院 大脑异常参数的提取方法和系统
CN104573742A (zh) * 2014-12-30 2015-04-29 中国科学院深圳先进技术研究院 医学图像分类方法和系统
CN105117731A (zh) * 2015-07-17 2015-12-02 常州大学 一种大脑功能网络的社团划分方法
KR101621849B1 (ko) * 2014-12-02 2016-05-31 삼성전자주식회사 뇌 네트워크의 분석을 위한 노드를 결정하는 방법 및 장치
CN106548206A (zh) * 2016-10-27 2017-03-29 太原理工大学 基于最小生成树的多模态磁共振影像数据分类方法
CN107392907A (zh) * 2017-09-01 2017-11-24 上海理工大学 基于静息态功能性磁共振成像的海马旁回功能分区方法
CN107658018A (zh) * 2017-10-12 2018-02-02 太原理工大学 一种基于结构连接和功能连接的融合脑网络构建方法
CN109147941A (zh) * 2018-10-17 2019-01-04 上海交通大学 基于结构磁共振影像数据的大脑健壮性评估方法
CN109921921A (zh) * 2019-01-26 2019-06-21 复旦大学 一种时变网络中时效稳定社团的检测方法和装置
CN110458832A (zh) * 2019-08-14 2019-11-15 电子科技大学 静息态大脑功能对称性分析方法
US10610124B2 (en) * 2012-08-09 2020-04-07 Brainlab Ag Localization of fibrous neural structures
US10687897B2 (en) * 2013-03-15 2020-06-23 Synaptive Medical (Barbados) Inc. System and method for health imaging informatics

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090129649A1 (en) * 2007-11-20 2009-05-21 Faycal Djeridane Method and system for processing multiple series of biological images obtained from a patient
CN102609946A (zh) * 2012-02-08 2012-07-25 中国科学院自动化研究所 一种基于黎曼流形的脑白质纤维束跟踪的组间处理方法
US10610124B2 (en) * 2012-08-09 2020-04-07 Brainlab Ag Localization of fibrous neural structures
US10687897B2 (en) * 2013-03-15 2020-06-23 Synaptive Medical (Barbados) Inc. System and method for health imaging informatics
CN103345749A (zh) * 2013-06-27 2013-10-09 中国科学院自动化研究所 一种基于模态融合的大脑网络功能连接偏侧性检测方法
CN104346530A (zh) * 2014-10-29 2015-02-11 中国科学院深圳先进技术研究院 大脑异常参数的提取方法和系统
KR101621849B1 (ko) * 2014-12-02 2016-05-31 삼성전자주식회사 뇌 네트워크의 분석을 위한 노드를 결정하는 방법 및 장치
CN104573742A (zh) * 2014-12-30 2015-04-29 中国科学院深圳先进技术研究院 医学图像分类方法和系统
CN105117731A (zh) * 2015-07-17 2015-12-02 常州大学 一种大脑功能网络的社团划分方法
CN106548206A (zh) * 2016-10-27 2017-03-29 太原理工大学 基于最小生成树的多模态磁共振影像数据分类方法
CN107392907A (zh) * 2017-09-01 2017-11-24 上海理工大学 基于静息态功能性磁共振成像的海马旁回功能分区方法
CN107658018A (zh) * 2017-10-12 2018-02-02 太原理工大学 一种基于结构连接和功能连接的融合脑网络构建方法
CN109147941A (zh) * 2018-10-17 2019-01-04 上海交通大学 基于结构磁共振影像数据的大脑健壮性评估方法
CN109921921A (zh) * 2019-01-26 2019-06-21 复旦大学 一种时变网络中时效稳定社团的检测方法和装置
CN110458832A (zh) * 2019-08-14 2019-11-15 电子科技大学 静息态大脑功能对称性分析方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DANDAN LI等: ""Reduced hemispheric asymmetry of brain anatomical networks"", 《BRAIN IMAGING AND BEHAVIOR (2019)》 *
ED BULLMORE等: ""The economy of brain network organization"", 《NATURE REVIEWS NEUROSCIENCE》 *
隆晓菁等: ""阿尔茨海默氏症患者大脑的结构偏侧性研究"", 《集成技术》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111784206A (zh) * 2020-07-29 2020-10-16 南昌航空大学 采用LeaderRank算法评估社交网络关键节点的方法
CN111784206B (zh) * 2020-07-29 2021-03-19 南昌航空大学 采用LeaderRank算法评估社交网络关键节点的方法
CN112741613A (zh) * 2021-01-13 2021-05-04 武汉大学 一种静息态人脑默认网络功能与结构耦合分析方法
CN113283465A (zh) * 2021-04-02 2021-08-20 电子科技大学 一种弥散张量成像数据分析方法及装置
CN113283465B (zh) * 2021-04-02 2022-04-29 电子科技大学 一种弥散张量成像数据分析方法及装置

Also Published As

Publication number Publication date
CN111340821B (zh) 2020-12-08

Similar Documents

Publication Publication Date Title
CN111340821B (zh) 一种基于模块连接的大脑结构网络的偏测性检测方法
CN109935336B (zh) 一种儿童呼吸科疾病的智能辅助诊断系统
Roncal et al. MIGRAINE: MRI graph reliability analysis and inference for connectomics
CN111582225B (zh) 一种遥感图像场景分类方法及装置
CN113616184A (zh) 基于多模态磁共振图像的脑网络建模与个体预测方法
WO2024083057A1 (zh) 基于多模态磁共振成像的图卷积神经网络疾病预测系统
CN106548206B (zh) 基于最小生成树的多模态磁共振影像数据分类方法
CN110503635B (zh) 一种基于异构数据融合网络的手骨x光片骨龄评估方法
CN112418337B (zh) 基于脑功能超网络模型的多特征融合数据分类方法
CN110598793A (zh) 一种大脑功能网络特征分类方法
CN103761537B (zh) 基于低秩优化特征字典模型的图像分类方法
CN112348833B (zh) 基于动态连接的脑功能网络变异识别方法及系统
CN117058514B (zh) 基于图神经网络的多模态脑影像数据融合解码方法和装置
Supandi et al. Two step cluster application to classify villages in Kabupaten Madiun based on village potential data
CN116110597B (zh) 一种基于数字孪生的患者病症类别智能分析方法及装置
CN115239674B (zh) 基于多尺度判别的计算机血管造影成像合成方法
CN111861924A (zh) 一种基于进化gan的心脏磁共振图像数据增强方法
CN114299006A (zh) 一种联合图对比学习的自适应多通道图卷积网络
WO2023108873A1 (zh) 一种脑网络和脑成瘾连接计算方法及装置
CN114821299B (zh) 一种遥感图像变化检测方法
Xu et al. Infrared and visible image fusion using a deep unsupervised framework with perceptual loss
CN110192860B (zh) 一种面向网络信息认知的脑成像智能测试分析方法及系统
CN115868923A (zh) 基于展开式循环神经网络的荧光分子断层成像方法及系统
CN117058170A (zh) 基于双分支多尺度交叉融合网络的颈动脉斑块分割方法
CN113283465B (zh) 一种弥散张量成像数据分析方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant