[发明专利]一种免阈值设定的静息态功能脑网络构建方法

说明书

本发明公开了一种免阈值设定的静息态功能脑网络构建方法，包括以下步骤：步骤S1：获取fMRI数据并提取各个脑区的BOLD信号矩阵；步骤S2：对BOLD信号进行信号处理并计算各个脑区间的皮尔森相关系数以此构建实数脑网络矩阵；步骤S3：将步骤S2中构建的实数脑网络矩阵输入神经网络并输出二值脑网络矩阵；步骤S4：采用遗传优化算法优化步骤S3输出的二值脑网络矩阵；步骤S5：如果步骤S4输出的二值脑网络矩阵符合优化条件，则获取最优二值脑网络矩阵；否则，重新配置神经网络参数，并重复步骤S2‑S5；步骤S6：根据步骤S5中最优二值脑网络矩阵产生免阈值静息态脑功能网络。

技术领域

本发明涉及静息态下功能脑网络构建技术领域，尤其涉及一种免阈值设定的静息态功能脑网络构建方法。

背景技术

大量来自功能磁共振成像(fMRI)的研究表明人处于静息状态(即大脑无特定任务状态)下，大脑区域间存在着低频自发神经元活动的同步性，这种同步活动具有重要的生理意义，同步的脑区构成了功能脑网络。基于图论的复杂网络分析技术使得关于大脑同步活动网络的拓扑结构研究成为可能。

目前三种传统的网络相似度比较方法，即基于退火算法的属性比较法，皮尔逊比较法和Jaccard比较法。基于退火算法的属性比较法是通过对两组网络的属性的统计检验然后利用模拟退火算法对其进行分析。Jaccard的比较方法主要是通过衡量网络重合度来比较两个网络的相似性，重合度越高相似度越强。皮尔逊相关的比较方法是将向量相关的比较方法扩展到矩阵相关，比较两个网络矩阵的相关性，若网络矩阵的相关性好，则网络的相似度高。

相关系数是由统计学家卡尔·皮尔逊设计的一种统计指标，用来反映变量之间相关关系密切程度。相关系数是通过积差方法计算，在两变量与各自平均值的离差基础上将两个离差相乘，以达到反映两变量之间相关程度的目的。这个相关系数就是皮尔逊相关。我们将这种对向量之间进行相关的方法引用到对矩阵的相关性计算中就得到了皮尔逊比较方法。具体而言，在脑网络研究中一个网络的载体是一个二值矩阵，比较网络之间的相似程度可以转化为比较代表这两个网络的两个二值矩阵的相似程度，而比较矩阵的相似程度可以用到皮尔逊相关系数，所以可以得到这样的结论：代表两个网络的二值矩阵的皮尔逊相关系数越大，那么两个网络越相似。当前主流的方法是人为设定一个阈值，比较皮尔逊相关系数与阈值的大小，当皮尔逊相关系数大于阈值时，将皮尔逊相关系数设置为“1”，认为脑区相关；反之将皮尔逊相关系数设置为“0”，,认为脑区不相关。由于阈值的选择受人为因素影响较大，如何选择合理的阈值从而构建真正的脑网络方法一致困扰着业界。

故，针对现有技术的缺陷，实有必要提出一种技术方案以解决现有技术存在的技术问题。

发明内容

有鉴于此，确有必要提供一种免阈值设定的静息态功能脑网络构建方法，采用预处理过的BOLD(血氧饱和度)信号，依次采用滑动窗口采样技术、皮尔森相关性检验技术、遗传算法等自主重构静息态功能脑网络作为认知大脑内部的工作机制，从而构建合理的脑网络用于脑疾病分析，为研究人脑内部的神经活动规律以及脑部疾病诊断提供有效的手段和依据。

为了解决现有技术存在的技术问题，本发明的技术方案如下：

一种免阈值设定的静息态功能脑网络构建方法，包括以下步骤：

步骤S1：获取fMRI数据并提取各个脑区的BOLD信号矩阵；

步骤S2：对BOLD信号进行信号处理并计算各个脑区间的皮尔森相关系数以此构建实数脑网络矩阵；

步骤S3：将步骤S2中构建的实数脑网络矩阵输入神经网络并输出二值脑网络矩阵；

步骤S4：采用遗传优化算法优化步骤S3输出的二值脑网络矩阵；

步骤S5：如果步骤S4输出的二值脑网络矩阵符合优化条件，则获取最优二值脑网络矩阵；否则，重新配置神经网络参数，并重复步骤S3-S5；