本发明专利技术提供了一种基于复杂网络的偏头痛辅助诊断测评方法。该方法首先对预处理后的图像进行分割，定义感兴趣区域，然后提取各个感兴趣区域的时间序列，再进行相关矩阵的计算，然后通过对相关矩阵进行二值化处理，来构建被试的静息态脑功能复杂网络。并进行复杂网络的拓扑结构性质计算，最后再将所得到被试的拓扑性质进行双样本T检验，并分析得到的显著异常。该方法对偏头痛患者进行测评之后，能很好地辅助对偏头痛患者的诊断。
【技术实现步骤摘要】
本专利技术属于医学影像学图像处理
，尤其是涉及一种基于复杂网络的偏头 痛患者大脑皮层的辅助分类分析方法。
技术介绍
偏头痛是一种常见的反复发作的慢性神经综合征，常伴有恶心、呕吐及畏光、畏声 等。随着神经影像技术在偏头痛研究上的深入，人们对偏头痛的生理病理机制的认识已逐 渐从最初的血管学说发展为神经血管学说，直至中央神经系统紊乱学说。并认为皮层扩散 性抑制是先兆性偏头痛重要病理生理机制。进一步研究发现，偏头痛患者局部大脑皮层的 结构和功能较之于正常人存在显著异常，这些区域包括：前额皮层，前喙扣带皮层，眶额皮 层以及脑岛等，而这些皮层在大脑中都参与到疼痛处理，这说明偏头痛的产生与大脑疼痛 处理机制的异常有关。 作为一种新型的、非侵入式的测量神经元自发活动的方法，静息态功能磁共振成 像技术越来越受到关注。有学者采用局部一致性以及低频振幅比率等方法分析偏头痛患者 的神经元自发性活动，发现在左侧前喙扣带皮层，双侧前额叶和右丘脑等区域，偏头痛患者 较之正常人会发生改变。 近年来，将图论知识应用于功能磁共振成像数据分析成为了研究热点，研究人员 通过定义人的大脑区域为节点、区域之间的功能连接为边，构建人脑网络。用图论的方法研 究发现许多神经疾病患者的（如：精神分裂症，阿尔兹海默症等）大脑静息态脑功能网络的 拓扑结构（如：聚类系数、特征路径长度等）较之于正常人会呈现显著异常，有学者认为这 些异常可以作为辅助相关疾病诊断的生理病理标志。
技术实现思路
本专利技术的目的是要针对目前技术上的不足，提出基于复杂网络的智能测评方法， 通过构建偏头痛患者大脑皮层的静息态脑功能复杂网络，并分析相关的网络拓扑结构测度 在偏头痛患者（MP)与正常人（NC)之间的统计学差异，实现偏头痛患者样本的分类。 为解决上述技术问题，本专利技术提供的技术步骤如下： 步骤1:按照临床医学诊断标准，分别采集正常人对照组（NC)和偏头痛患者组 (MP)的大脑静息态功能磁共振图像，两组样本的数目相等或者接近，获取过程中被试要求 保持大脑清醒，平躺于磁共振仪器内，不做任何定性思考。 步骤2 :对所采集的被试数据进行预处理。 所有数据均去除了前10个时间点，以消除被试在习惯实验环境过程中受到不均 衡磁场的影响。预处理流程主要包括时间层矫正（slice timing),头动矫正（realign)、图 像标准化（normalization)以及平滑处理（smooth)等操作。首先，所有的数据集采用正弦 内插方法消除时间偏移，进而再用六度间变换方法的方法消除空间偏移；其次，为了尽量减 少伪影，在任何方向上位移大于1. 5mm或头部转动大于1. 5°的时间片都会被丢弃；第三， 所有数据集进行空间标准化时，都以蒙特利尔神经学研究所（MNI)平面回波成像图像为模 进行重采样；最后，通过带通滤波器（〇. 01Hz < f < 0. 1Hz)，以除去低频漂移和高频生理噪 声。 步骤3 :计算正常人对照组和偏头痛患者组两组大脑静息态功能磁共振图像的大 脑功能连接矩阵。 本专利技术采用的是自动解剖标签模板，该模板将大脑划分为90个感兴趣区域，在网 络构建中，我们将这90个大脑区域定义为一系列节点；然后，将所有脑区的时间序列都通 过带通滤波（〇. 01-0. 〇8Hz)进行过滤，以去除磁场漂移和高频噪音的影响；随后，将大脑白 质深部的平均时间过程和6个刚体运动参数从已滤波后的时间序列中进行回归分析；最 后，再计算每个种子区域的平均时间序列，并获得由所有可能连接节点之间皮尔森偏相关 系数所组成的90*90的矩阵。 步骤4 :构建两组大脑的复杂网络。 正如先前研究所提到的，我们首先将得到的90*90的矩阵进行阈值化处理转化为 连接矩阵，即：如果两个节点之间的相关系数大于特定的阈值，我们就认为这两个节点之间 存在功能连接；否则，我们就认为两个节点之间不存在功能连接。大脑的功能连接可以用无 向边表示。随后，我们再对得到的连接矩阵进行二值化处理，即：在连接矩阵中，用"1"表示 两个节点之间有连接，"〇"表示两个节点之间没有连接。经过这两步处理之后，我们便获得 了大脑的静息态脑功能复杂网络的图。通过这种方式，我们获得的图是由相同数量的连接 组成的，这可以使得正常人和偏头痛患者这两组被试的皮层网络具有相同的布线成本。我 们将相关矩阵通过一个较宽的阈值范围（15% -25% )，然后估算每一个对应阈值的网络性 质。考虑到正常人和患者的异常连接和皮层中枢，阈值为（S) = 15%在随后分析中最为典 型。这个值是能够避免网络中产生孤立点（皮质网络完全连接）的最低阈值，同时，它还可 以尽量减少网络中杂散边的数量。这个阈值在先前的很多图论分析研究中也有应用。 步骤5 :计算两组被试的复杂网络拓扑结构 1?聚类系数 聚类系数是表示一个网络中节点的集群连接程度的系数，它反应的是大脑网络中 功能分割的情况。通常，聚类系数Ci被定义为：【主权项】1. 一种基于复杂网络的偏头痛患者大脑皮层辅助分类分析方法，包括w下几个步骤： 步骤1;按照临床医学诊断标准，分别采集正常人对照组和偏头痛患者组的大脑皮层 静息态功能磁共振图像，两组样本的数目相等或者接近； 步骤2 ;对所采集的所述正常人对照组和偏头痛患者组的大脑皮层静息态功能磁共振 图像的数据进行预处理，预处理流程主要包括时间层矫正、头动矫正、图像标准化W及平滑 处理等操作； 步骤3 ;分别计算所述正常人对照组和所述偏头痛患者组大脑皮层静息态功能磁共振 图像两组数据的大脑功能连接矩阵，采用自动解剖标签模板，该模板将大脑划分为90个感 兴趣区域，通过提取所述感兴趣区域的时间序列，进而计算每个种子区域的平均时间序列， 并获得由所有可能连接节点之间皮尔森偏相关系数所组成的分别对应于正常人对照组和 偏头痛患者组的两组90*90的矩阵； 步骤4 ;将所述的两组90*90的矩阵分别进行阔值化处理转化为连接矩阵，随后，再对 得到的连接矩阵进行二值化处理，获得两组大脑静息态脑功能复杂网络图； 步骤5 ;计算所述所述两组大脑静息态脑功能复杂网络图的复杂网络拓扑结构，即计 算其聚类系数、特征路径长度、介数中也度、同配性、和小世界性； 步骤6 ;对两组所述正常人对照组和偏头痛患者组的大脑皮层静息态功能磁共振图像 数据进行双样本T检验统计分析，并采用假发现率校正方法对多重比较的P值进行校正； 步骤7 ;通过对步骤5的所述复杂网络拓扑结构进行统计分析及特征提取，提取所述两 组大脑静息态脑功能复杂网络图中的尾状核、豆状壳核等节点的特征路径长度，提取网络 中丘脑、左侧枕下回、右侧枕中回等节点的介数中也度，并统计所述小世界性和同配性。2. 根据权利要求1所述的基于复杂网络的偏头痛患者大脑皮层辅助分类分析方法，其 特征是：步骤5中所述聚类系数、特征路径长度、介数中也度、同配性、和小世界性是指： (1) 聚类系数，聚类系数是表示一个网络中节点的集群连接程度的系数，它反应的是大 脑网络中功能分割的情况； (2) 特征路径长度，特征路径长度也称为平均最短路径长度，它反应的是大脑的功能集 成的情况； (3) 介数中也度，介数中也度是反应的节点在网络中重要性的重要本文档来自技高网...
【技术保护点】
一种基于复杂网络的偏头痛患者大脑皮层辅助分类分析方法，包括以下几个步骤：步骤1：按照临床医学诊断标准，分别采集正常人对照组和偏头痛患者组的大脑皮层静息态功能磁共振图像，两组样本的数目相等或者接近； 步骤2：对所采集的所述正常人对照组和偏头痛患者组的大脑皮层静息态功能磁共振图像的数据进行预处理，预处理流程主要包括时间层矫正、头动矫正、图像标准化以及平滑处理等操作； 步骤3：分别计算所述正常人对照组和所述偏头痛患者组大脑皮层静息态功能磁共振图像两组数据的大脑功能连接矩阵，采用自动解剖标签模板，该模板将大脑划分为90个感兴趣区域，通过提取所述感兴趣区域的时间序列，进而计算每个种子区域的平均时间序列，并获得由所有可能连接节点之间皮尔森偏相关系数所组成的分别对应于正常人对照组和偏头痛患者组的两组90*90的矩阵；步骤4：将所述的两组90*90的矩阵分别进行阈值化处理转化为连接矩阵，随后，再对得到的连接矩阵进行二值化处理，获得两组大脑静息态脑功能复杂网络图； 步骤5：计算所述所述两组大脑静息态脑功能复杂网络图的复杂网络拓扑结构，即计算其聚类系数、特征路径长度、介数中心度、同配性、和小世界性； 步骤6：对两组所述正常人对照组和偏头痛患者组的大脑皮层静息态功能磁共振图像数据进行双样本T检验统计分析，并采用假发现率校正方法对多重比较的P值进行校正； 步骤7：通过对步骤5的所述复杂网络拓扑结构进行统计分析及特征提取，提取所述两组大脑静息态脑功能复杂网络图中的尾状核、豆状壳核等节点的特征路径长度，提取网络中丘脑、左侧枕下回、右侧枕中回等节点的介数中心度，并统计所述小世界性和同配性。...
【技术特征摘要】
【专利技术属性】
技术研发人员：陈敦耀，曾卫明，
申请(专利权)人：上海海事大学，
类型：发明
国别省市：上海;31
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。