[发明专利]一种探测fMRI脑网络动态协变的方法

权利要求书

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的均值，X_i、Y_i分别表示在第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≤LN，滑动步长n，1≤nL，将每个ROI中的被试度分布序列1×N向量划分为段的1×L向量，其中S为滑动窗的个数，N为样本中被试人数，L、n分别为以时间点为单位滑动窗的长度以及步长，为向下取整函数；

(7)在一个滑动窗内计算每个ROI脑区度序列之间的皮尔逊相关系数，得到一个M×M的ROI脑区间协变矩阵，对所有滑动窗重复此操作从而得到一M×M×S的协变矩阵序列；

(8)得到协变矩阵序列后分别计算每个ROI脑区在这S个协变矩阵内的变异情况。

2.根据权利要求1所述的探测fMRI脑网络动态协变的方法，其特征在于，步骤(8)所述变异情况通过方差、变异系数来刻画。

3.根据权利要求1所述的探测fMRI脑网络动态协变的方法，其特征在于，当有多组被试时，通过组间统计检验比较组别之间的差异。

4.根据权利要求1所述的探测fMRI脑网络动态协变的方法，其特征在于，所述步骤(4)中所有被试的网络矩阵利用行为学数据进行排序得到跨被试的网络矩阵序列。

5.根据权利要求4所述的探测fMRI脑网络动态协变的方法，其特征在于，步骤(4)中所述行为学数据为按照被试的年龄大小、病程长短、教育程度高低进行排列。

6.根据权利要求1所述的探测fMRI脑网络动态协变的方法，其特征在于，所述滑动窗为：矩形窗、高斯窗或者汉明窗，窗长以及步长的选择可以灵活选取。