一种采集大脑功能连接自发波动变异性方法技术

本发明专利技术公开了一种采集大脑功能连接自发波动变异性方法，包括如下步骤：步骤一，利用磁共振扫描仪记录人脑静息态过程中的血氧水平依赖信号，获得fMRI数据；步骤二，对fMRI数据进行滤波处理选取低频信号；步骤三，对低频信号进行相关系数的度量算法获得功能连接值FC；步骤四，在选取的低频信号以及功能连接值FC的基础上通过滑动窗法构建动态功能连接序列dFC；步骤五，对动态功能连接序列dFC进行低通滤波处理后按照公式(1)获得动态功能连接序列dFC＝[FC

A method of collecting spontaneous fluctuations in functional connectivity of the brain

The invention discloses a collection of functional connectivity in the brain of spontaneous fluctuations of variation method, which comprises the following steps: step one, using blood oxygen level magnetic resonance scanner records the brain resting state dependence in the process of signal acquisition of fMRI data; step two, the fMRI data filtering processing selected frequency signal; step three, the low frequency signal the correlation coefficient measurement algorithm to obtain the functional connectivity value FC; step four, in the selection of the low frequency signal and the function based on FC connection value by sliding window method to construct a dynamic functional connection sequence dFC; step five, the connection sequence dFC low-pass filtering dynamic function according to the formula (1) connected sequence dFC = [FC dynamic function

【技术实现步骤摘要】
一种采集大脑功能连接自发波动变异性方法
本专利技术属于静息态功能磁共振技术的动态功能连接
，特别涉及一种采集大脑功能连接自发波动变异性方法。
技术介绍
人脑是一个极其复杂且高效的网络系统，脑区间的协调工作将空间分离的脑区联系起来，使其构成一个拓扑网络。基于功能连接的研究表明，利用fMRI技术采集到的一些脑区之间的BOLD信号存在高度的时域同步性。研究人员通过采集人脑在特定生理状态下的fMRI数据，计算不同脑区的BOLD信号之间的相关性，构建大尺度的脑区功能连接网络，从而来估计在这种生理状态下不同脑区之间的联系。尽管基于静息态fMRI技术的大脑功能连接的研究引起越来越多研究者的关注，并且在研究大脑认知功能领域取得了很大的进步。但是，静息态功能连接(resting-statefunctionconnectivity,RSFC)的估计在很大程度上是受限的，因为以往的研究中，fMRI采集的时间通常是5到30分钟，这隐含的假设在采集时间内大脑的信号之间的关联关系在时间和空间上是恒定的，也就是，扫描时间内脑区之间静息态功能连接为一个固定不变的值。这个假设避免了整个大脑连接极其复杂的分析，为研究带来了很多方便。但是，它太过于简单，这个假设表示的是复杂的时间和空间现象的平均值。因此，为了更深入的了解大脑网络的基本性能，在最近几年有人提出并证明了动态功能连接的概念：静息态下大脑的功能连接随着时间的变化存在一定的波动性。但是怎么描述功能连接的波动性，就成了一个亟待解决的问题。
技术实现思路
针对现有技术存在不足，本专利技术提出一种采集大脑功能连接自发波动变异性方法，本专利技术旨在利用图像数据，进行特征提取，进而用于数据分析，填补了描述动态功能连接变异性的指标缺乏的空白。为了解决现有技术中存在技术问题，本专利技术采用如下技术方案：一种采集大脑功能连接自发波动变异性方法，包括如下步骤：步骤一，利用磁共振扫描仪记录人脑静息态过程中的血氧水平依赖信号，获得fMRI数据；步骤二，对fMRI数据进行滤波处理选取低频信号；步骤三，对低频信号进行相关系数的度量算法获得功能连接值FC；步骤四，在选取的低频信号以及功能连接值FC的基础上通过滑动窗法构建动态功能连接序列dFC；步骤五，对动态功能连接序列dFC进行低通滤波处理后按照公式(1)获得动态功能连接序列dFC＝[FC1,FC2,...FCn]的平均值；步骤六，对步骤五中动态功能连接序列dFC＝[FC1,FC2,...FCn]按照公式(2)获得功能连接变异性模型。所述步骤二中的低频信号范围为0.01Hz-0.08Hz。本专利技术有益效果：第一，本专利技术提出了一种定量检测两个脑区之间功能连接动态变化的方法，量化了功能连接自发波动的幅度，反映了功能连接自发波动或时变的能力。第二，本专利技术证明了动态功能连接的存在，填补了描述动态功能连接变异性的指标缺乏的空白，为大脑的功能连接的研究提供了方法学的支持。第三，本专利技术实现一体功能连接自发波动的幅度，反映功能连接自发波动或时变的能力，进而为后续研究的分析提供基础数据。附图说明图1是本专利技术一种采集大脑功能连接自发波动变异性的方法流程图；图2是本专利技术涉及动态功能连接曲线图；图3是本专利技术涉及功能连接变异性矩阵图；图4是本专利技术涉及青年组与老年组的功能连接变异性结果比较图。具体实施方式下面结合附图对本专利技术做出详细地说明：本专利技术设计基于一种定量测量大脑功能连接自发波动的量化指标和技术方法，可应用到基于静息态功能磁共振技术的动态功能连接研究。其技术流程是：用磁共振扫描仪记录人脑静息态过程中的血氧水平依赖(BloodOxygenLevelDependence,BOLD)信号，提取时间序列，按照一定的方法，构建动态功能连接序列，然后计算其动态功能连接序列的标准差，即功能连接变异性(FunctionalConnectivityVariability,FCV)，用来描述大脑功能连接的波动性。基于fMRI数据设计定量描述大脑功能连接自发波动指标方法的流程图如图1所示。流程包括，数据采集，数据预处理(滤波处理)，动态功能连接序列的构建以及特征提取。具体内容如下：如图1所示，本专利技术为一种采集大脑功能连接自发波动变异性方法，包括如下步骤：步骤一101，利用磁共振扫描仪记录人脑静息态过程中的血氧水平依赖信号，获得fMRI数据；步骤二102，对fMRI数据进行滤波处理选取低频信号；所述步骤二中的低频信号范围为0.01Hz-0.08Hz。本专利技术的数据描述是基于实现的过程中所用到的数据为在磁共振扫描仪中用功能磁共振成像序列检测样品得到的时间序列信号。取样本中任意两个空间体素或者区域平均(记为N1,N2)的时间序列，设信号为S1和S2，数据长度为L个时间点，采样周期为TR。数据采集之后，需要进行一个滤波处理，即将某个频段的信号去除，功能磁共振成像中数据的预处理排除高频信号，因为只有低频信号(0.01Hz-0.08Hz)能够反映大脑的自发神经活动。步骤三103，对低频信号进行相关系数的度量算法获得功能连接值FC；本专利技术中功能连接(functionconnectivity,FC)的计算分为两个步骤，首先计算不同脑区信号之间的相关系数corr，然后相关系数经过fisherZ变换，得到功能连接值FC；本研究使用的相关系数的度量方法是皮尔森相关分析方法。皮尔森相关系数是Pearson在19世纪80年代提出的，用于度量两个变量X与Y之间的线性相关，其值介于-1与1之间。符号的变化反映变量之间相关性的方向，正数表示正相关，负数表示负相关。其中，1表示X和Y具有很好的线性关系，所有的数据点可以很好地拟合成一条直线，且Y随着X的增加而增加；-1表示X与Y具有很好的线性相关，然而Y随X的增加而减少；0表示不存在相关性。计算公式如下：fisherZ变换的作用是使功能连接值符合正态分布。corr经过fisherZ变换可以得到功能连接FC。步骤四104，在选取的低频信号以及功能连接值FC的基础上通过滑动窗法构建动态功能连接序列dFC；本专利技术中动态功能连接构建采用滑动窗法。首先选取一定宽度的时间窗，设宽度为w，然后按照一定的步长，记为s，依次在采集时间序列上滑动，并计算每一个滑动窗下的时间序列两两之间的相关性，然后经过fisherZ变换，得到该时刻的功能连接值，记为FCi(i为时刻)，逐步滑动时间窗，即可得到一个随时间变化的功能连接序列[FC1,FC2,...FCn]，其中，n＝L-w+1，这个序列即为动态功能连接(dynamicfunctionalconnectivity,dFC)，即dFC＝[FC1,FC2,...FCn]，具体流程如图2中显示的是24TR的时间窗，通过按照一定的步长依次滑动，得到172个FC值，形成一个随时间变化的功能连接序列，即为动态功能连接序列，这个序列能反应出两个节点之间的功能连接随时间动态变化的情况。步骤五105，对动态功能连接序列dFC进行低通滤波处理后按照公式(1)获得动态功能连接序列dFC＝[FC1,FC2,...FCn]的平均值；步骤六106，对步骤五中动态功能连接序列dFC＝[FC1,FC2,...FCn]按照公式(2)获得功能连接变异性模型。本专利技术构建完成了动态功能连接序列，即需要一本文档来自技高网...

【技术保护点】
一种采集大脑功能连接自发波动变异性方法，其特征在于：包括如下步骤：步骤一，利用磁共振扫描仪记录人脑静息态过程中的血氧水平依赖信号,获得fMRI数据；步骤二，对fMRI数据进行滤波处理选取低频信号；步骤三，对选取的低频信号进行相关系数的度量算法获得功能连接值FC；步骤四，在选取的低频信号以及功能连接值FC的基础上通过滑动窗法构建动态功能连接序列dFC；步骤五，对动态功能连接序列dFC进行低通滤波处理后按照公式(1)获得动态功能连接序列dFC＝[FC

【技术特征摘要】
1.一种采集大脑功能连接自发波动变异性方法，其特征在于：包括如下步骤：步骤一，利用磁共振扫描仪记录人脑静息态过程中的血氧水平依赖信号,获得fMRI数据；步骤二，对fMRI数据进行滤波处理选取低频信号；步骤三，对选取的低频信号进行相关系数的度量算法获得功能连接值FC；步骤四，在选取的低频信号以及功能连接值FC的基础上通过滑动窗法构建动态功能连接序列dFC；步骤五，对动态功能连接序列dFC进行低通滤波处理后按照公式(1)获得动态功能连接序列dFC＝[FC...

【专利技术属性】
技术研发人员：赵欣，张雄，陈元园，周鹏，綦宏志，明东，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津,12