一种探测fMRI脑网络动态协变的方法技术

本发明专利技术提供一种探测fMRI脑网络动态协变的方法，该方法包括：首先利用功能磁共振数据构建大脑功能网络，再按照一定的依据(年龄大小、病程长短、教育程度高低等因素)分别将不同的网络矩阵顺序排列起来，得到跨被试的网络矩阵序列，而后利用滑动窗方法进行跨被试的ROI脑区之间的协变关系研究，得到一系列协变矩阵，再通过计算刻画协变关系变异性的数学指标或通过统计检验对跨被试尺度下脑区间协同活动模式的差异，从而更好的理解在网络中脑区之间的交互关系。

【技术实现步骤摘要】
一种探测fMRI脑网络动态协变的方法
本方法属于生物医学图像模式识别
，具体涉及一种探测fMRI脑网络动态协变的方法。
技术介绍
功能磁共振成像(functionalmagneticresonanceimaging,fMRI)技术为我们提供了一种无创的观测大脑的方法，其利用血氧水平依赖(BOLD)信号可以很好的反应大脑的功能活动。人的大脑各个部分之间的联系是非常错综复杂而且是具有一定特殊的模式结构的，近年来脑连接组(brainconnectome)概念的提出以及快速的发展，越来越多的研究都从脑网络的角度来考量大脑的功能活动、结构特征等。现在的脑科学已经从以前的将大脑划分为独立的功能分区转变为探索各脑区之间是如何互相影响以及协同关系的变化情况。目前，已有研究利用滑动窗的方式来构建脑动态网络,但是，目前的方式纵向数据获取难度大，不能很好的探究疾病以及人脑认知功能等与脑活动模式间的关系。
技术实现思路
本专利技术的目的在于解决上述现有技术存在的缺陷，利用滑动窗方法度量跨被试尺度上各脑区之间的协同变化关系，提出一种基于fMRI探测脑功能网络动态协变的方法。一种探测fMRI脑网络动态协变的方法，具体步骤如下：(1)获取原始血氧依赖水平的功能磁共振(BOLD-fMRI)数据图像，将原始四维数据矩阵转换为二维V×T矩阵，其中，四维数据矩阵为三维空间×一维时间，V为磁共振图像中的体素个数，T为时间点个数；并且根据事先定义好的空间感兴趣区ROI模板，在所定义的同一ROI区域内将所有体素的BOLD-fMRI信号求取均值，并以此均值时间信号作为该ROI的BOLD-fMRI信号值，从而将V×T矩阵压缩为M×T矩阵，其中M为空间ROI模板所定义的大脑区域数量，T为时间点长度；(2)计算所有ROI脑区之间两两配对得到的皮尔逊相关系数，其中两个ROI脑区之间的皮尔逊相关系数，计算如下式：其中，r为两个ROI之间BOLD-fMRI信号之间的相关系数，X、Y为一个内任意两个ROI的时间序列，T为时间序列中时间点的个数，分别为X、Y的均值，Xi、Yi分别表示在第i个时间点A、B两个ROI的信号值；(3)对所有被试重复(2)中的计算过程，从而在每个被试中都可以得到一个网络连接矩阵M×M，定义网络稀疏度将每个网络矩阵作稀疏化处理；其中M为脑区划分ROI的个数；(4)将N名被试的网络矩阵按照一定的依据排列起来，得到网络矩阵序列M×M×N；(5)针对所述网络矩阵序列M×M×N，求取其节点的度分布，从而得到网络度分布矩阵序列M×N，其中M为ROI脑区个数，N为被试人数；(6)利用滑动窗思想，对于网络度分布矩阵序列中每个ROI脑区，均产生一个1×N序列，N为样本中被试人数，针对这一序列根据预先设定好的滑动窗长度L，1≤L<N，，滑动步长n，1≤n<L，将每个ROI中的被试度分布序列1×N向量划分为段的1×L向量，其中S为滑动窗的个数，N为样本中被试人数，L、n分别为以时间点为单位滑动窗的长度以及步长，为向下取整函数；(7)在一个滑动窗内计算每个ROI脑区度序列之间的皮尔逊相关系数，得到一个M×M的ROI脑区间协变矩阵，对所有滑动窗重复此操作从而得到一M×M×S的协变矩阵序列；(8)得到协变矩阵序列后分别计算每个ROI脑区在这S个协变矩阵内的变异情况。进一步地，如上所述的探测fMRI脑网络动态协变的方法，步骤(8)所述变异情况通过方差、变异系数来刻画。进一步地，如上所述的探测fMRI脑网络动态协变的方法，当有多组被试时，通过组间统计检验比较组别之间的差异。进一步地，如上所述的探测fMRI脑网络动态协变的方法，所述步骤(4)中所有被试的网络矩阵利用行为学数据进行排序得到跨被试的网络矩阵序列。进一步地，如上所述的探测fMRI脑网络动态协变的方法，步骤(4)中所述行为学数据为按照被试的年龄大小、病程长短、教育程度高低进行排列。进一步地，如上所述的探测fMRI脑网络动态协变的方法，所述滑动窗为：矩形窗、高斯窗或者汉明窗，窗长以及步长的选择可以灵活选取。有益效果：该方法是按照一定的人口学信息以及行为信息、疾病症状信息等为依据顺序将被试排列，从而实现在跨被试尺度上研究排列依据因素对大脑活动的影响以及两者之间的关系，可以作为纵向研究的一种替代方案，解决纵向数据获取难度大的问题，更好的探究疾病以及人脑认知功能等与脑活动模式间的关系。附图说明图1为ROI脑区时间序列信号图；图2为稀疏化后的网络矩阵及连接图；图3为网络矩阵序列图；图4为度分布序列滑窗示意图；图5为协变矩阵序列以及所得各ROI脑区方差示意图。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚，下面本专利技术中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。本专利技术提供一种探测fMRI脑网络动态协变的方法，属于生物医学图像模式识别
，具体涉及基于功能磁共振时间序列的脑功能网络动态协变探测方法，该方法包括：首先利用功能磁共振(fMRI)数据构建大脑功能网络，再按照一定的依据(年龄大小、病程长短、教育程度高低等因素)分别将不同的网络矩阵顺序排列起来，得到跨被试的网络矩阵序列，而后利用滑动窗方法进行跨被试的ROI脑区之间的协变关系研究，得到一系列协变矩阵，再通过计算刻画协变关系变异性的数学指标或通过统计检验对跨被试尺度下脑区间协同活动模式的差异，从而更好的理解在网络中脑区之间的交互关系。该方法具体实施步骤如下：(1)将一组正常被试共计30人的原始功能磁共振数据，经过一系列预处理步骤消除存在的噪声影响，其主要操作包括：删除前5个时间点；以第30层为参考层进行时间层校正；将所有时间点的数据空间对其到第一个时间点，并估算包括3个平动、3个转动头动参数、以及层间位移等参数，并依据平动不超过1毫米，转动不超过1°，平均层间位移不超过0.2毫米的标准筛查被试数据；依照EPI模板将数据标准化到MNI(蒙特利尔神经学研究所)人脑标准空间去；将24个头动参数，平均白质、脑脊液信号，以及全脑信号通过广义线性模型(GLM)作回归处理；去除线性漂移；以0.01-0.08赫兹的频段进行带通滤波；从而得到预处理好的原始数据。具体地，原始功能磁共振数据中的噪声信号非常多，如果不做处理的话会带来严重的影响，极大削弱结果的可信度。删除时间点是消除核磁共振仪器开始时磁场不稳定的影响；时间层校正以及头动参数的估计都是为了消除头动的影响，头动在该领域是公认的极大的数据污染因素；标准化到MNI空间是为了消除个体大脑结构差异造成的不可比性，将所有个体的大脑都标准化到相同的空间上去，从而才可以实现后续的统计比较分析；带通滤波是为了消除其他频段的噪声影响大脑活动信号的频段集中在0.01-0.08Hz频段上，这样做可以最大程度上的去除其他频段噪声的影响，保留有用目标信号。(2)将原始四维数据矩阵(三维空间×一维时间)转换为二维V×T矩阵，其中V＝61×73×61＝271633为磁共振图像中的体素个数，T＝240为时间点长度；并且根据事先定义好的空间感兴趣区(ROI)模板，在所定义本文档来自技高网...

【技术保护点】
1.一种探测fMRI脑网络动态协变的方法，其特征在于，包括以下步骤：(1)获取原始血氧依赖水平的功能磁共振(BOLD‑fMRI)数据图像，将原始四维数据矩阵转换为二维V×T矩阵，其中，四维数据矩阵为三维空间×一维时间，V为磁共振图像中的体素个数，T为时间点个数；并且根据事先定义好的空间感兴趣区ROI模板，在所定义的同一ROI区域内将所有体素的BOLD‑fMRI信号求取均值，并以此均值时间信号作为该ROI的BOLD‑fMRI信号值，从而将V×T矩阵压缩为M×T矩阵，其中M为空间ROI模板所定义的大脑区域数量，T为时间点长度；(2)计算所有ROI脑区之间两两配对得到的皮尔逊相关系数，其中两个ROI脑区之间的皮尔逊相关系数，计算如下式：

【技术特征摘要】
1.一种探测fMRI脑网络动态协变的方法，其特征在于，包括以下步骤：(1)获取原始血氧依赖水平的功能磁共振(BOLD-fMRI)数据图像，将原始四维数据矩阵转换为二维V×T矩阵，其中，四维数据矩阵为三维空间×一维时间，V为磁共振图像中的体素个数，T为时间点个数；并且根据事先定义好的空间感兴趣区ROI模板，在所定义的同一ROI区域内将所有体素的BOLD-fMRI信号求取均值，并以此均值时间信号作为该ROI的BOLD-fMRI信号值，从而将V×T矩阵压缩为M×T矩阵，其中M为空间ROI模板所定义的大脑区域数量，T为时间点长度；(2)计算所有ROI脑区之间两两配对得到的皮尔逊相关系数，其中两个ROI脑区之间的皮尔逊相关系数，计算如下式：其中，r为两个ROI之间BOLD-fMRI信号之间的相关系数，X、Y为一个内任意两个ROI的时间序列，T为时间序列中时间点的个数，分别为X、Y的均值，Xi、Yi分别表示在第i个时间点A、B两个ROI的信号值；(3)对所有被试重复(2)中的计算过程，从而在每个被试中都可以得到一个网络连接矩阵M×M，定义网络稀疏度将每个网络矩阵作稀疏化处理；其中M为脑区划分ROI的个数；(4)将N名被试的网络矩阵按照一定的依据排列起来，得到网络矩阵序列M×M×N；(5)针对所述网络矩阵序列M×M×N，求取其节点的度分布，从而得到网络度分布矩阵序列M×N，其中M为ROI脑区个数，N为被试人数；(6)利用滑动窗思想，对于网络度分...

【专利技术属性】
技术研发人员：廖伟，孟耀，陈华富，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51