基于独立成分高序不确定脑网络的磁共振影像分类方法技术

本发明专利技术公开了基于独立成分高序不确定脑网络的磁共振影像分类方法，涉及图像处理技术。先使用独立成分分析提取默认网络独立成分，接着在独立成分的时间序列上利用时间窗滑动的方法构建高阶功能连接网络，然后使用基于加权图的频繁子图挖掘和特征选择方法，最后以判别性子图作为特征对进行分类识别。该方法不需要依赖先验的脑图谱模板，充分考虑了扫描时间内的时变特性，且将频繁子图挖掘应用到了加权图上，呈现出了更高级、更复杂的大脑区域之间的相互作用。解决了传统的磁共振成像数据分类方法准确率低的问题。

【技术实现步骤摘要】
基于独立成分高序不确定脑网络的磁共振影像分类方法
本专利技术属于图像处理
，特别是涉及一种基于独立成分高序不确定脑网络的磁共振影像分类方法。
技术介绍
人脑是一个高度复杂的系统，探索其内在的结构和功能是一项巨大的挑战。功能磁共振成像技术与复杂网络理论的结合，已成为脑科学领域的研究热点之一，并广泛应用在各项研究中。该方法在探索大脑区域间结构和功能的相互作用上取得许多令人惊喜的成果。然而，由于其自身原理和特点及方法论的限制，导致分类准确率低，严重影响了应用价值。近年来大多数基于脑网络的研究，在定义网络节点时都是按照脑图谱模板定义的。只是这种节点的定义方法过于依赖先验的脑图谱模板，不同模板建立起来的网络可能会存在显著差异，从而严重影响了功能连接网络的分析结果。其次，传统的功能连接网络分析中忽略了大脑区域之间功能连接的动态变化，即边的高阶关系。然而已有研究发现这些高阶信息的损失对于疾病诊断可能是重要的。此外，目前大多数频繁子图挖掘的研究集中在确定图上，判别性子图特征的选择和各种属性的计算也都是基于确定图的。虽然确定图分析起来比较简单，但由于确定图忽略了脑网络功能的不确定性，不能全面准确地反映脑区之间的交互。为了解决这些问题，本专利技术先使用独立成分分析选择静息态脑网络成分，然后在其上构建高序功能连接网络，最后对高序不确定功能脑网络进行频繁子图挖掘。这种方法可以更好地反映了全脑功能连接的时频特性，呈现了更高级、更复杂的大脑区域之间的相互作用，由此提高了分类准确率。
技术实现思路
本专利技术的目的在于提供一种基于独立成分高序不确定脑网络的磁共振影像分类方法，分类准确率高。为解决上述技术问题，本专利技术所采用的技术方案是，基于独立成分高序不确定脑网络的磁共振影像分类方法，按照以下步骤进行：步骤S1：对静息态功能磁共振成像数据进行预处理，然后使用独立成分分析方法提取出各独立成分；步骤S2：筛选出属于默认网络的独立成分，提取时间序列，并进行后处理；步骤S3：选定长度固定的滑动窗口，按照1s的移动步长，划分默认网络独立成分的时间序列；步骤S4：计算每个时间窗下默认网络成分的时间序列两两间的皮尔逊相关系数，得到低阶关联矩阵；步骤S5：提取每个时间窗下低阶关联矩阵中的元素值，生成高阶关联矩阵；步骤S6：构建高序不确定功能脑网络，节点之间全连接，节点之间边的权值即高阶关联矩阵中的元素；步骤S7：对构建的全连接加权的脑网络进行频繁子图挖掘，得到加权的频繁子图模式；步骤S8：计算基于动态规划的判别性函数的相关统计指标；步骤S9：采用支持向量机分类算法，选择判别性子图特征作为分类依据，构建分类器，然后采用交叉验证方法进行检验；步骤S10：采用互信息分析方法，量化所选判别性子图特征在分类器中的重要度和冗余度，然后根据量化结果对所选判别子图进行二次筛选，优化特征子集。进一步的，所述步骤S1中，预处理步骤包括时间层校正，头部校正，空间标准化，平滑；然后，使用基于负熵最大的独立成分分析方法提取独立成分；其中，负熵的定义为：Nh(Y)＝{E[h(Y)]-E[h(YGauss)]}2；(1)公式(1)中，Y是源信号的随机变量，Nh(Y)为随机变量Y的负熵；YGauss是一与Y具有相同方差的高斯随机变量；E为均值运算；h为非线性函数，这里为使Y具有最大的非高斯性，即找到Y负熵的最大近似值，可通过对E{h(Y)}进行优化来获得，而E{h(Y)}的最适条件可通过下式获得：Y＝WTXE{Xh(WTX)}+βW＝0；(2)公式(2)中，Y是源信号的随机变量；X是观测信号的随机变量；W是解混矩阵，WT是W的转置；E为均值运算；h为非线性函数；β是一个恒定值；用牛顿迭代方法解此方程后，可以获得近似牛顿迭代公式，简化后就是独立成分分析方法的迭代公式：W+＝E{Xh(WTX)}-E{h′(WTX)}W；(3)公式(3)中，h′是h的导数；W+是W的新值，W＝W+/||W+||。进一步的，所述步骤S2中，需要使用默认网络模板匹配的方法对独立成分进行筛选，接着对每个被试默认网络成分的时间序列进行后处理，包括去线性漂移和低通滤波，以去除线性漂移和高频生理噪声的影响。进一步的，所述步骤S4中，计算每个时间窗下默认网络成分的时间序列两两间的皮尔逊相关系数，得到低阶关联矩阵的具体公式是：公式(4)中，rij表示每个时间窗下关联矩阵中第i行第j列的元素，是默认网络成分成分i与j之间的皮尔逊相关系数；n表示时间点个数；xi(t)表示第i个默认网络成分的时间序列；表示第i个默认网络成分的平均时间序列。进一步的，所述步骤S6中，在高阶关联矩阵上，将其边的权值按照如下公式进行处理，得到的即为高序不确定功能脑网络：公式(5)中，bij表示高序不确定功能脑网络模型中第i行第j列的元素；rij表示高阶关联矩阵中第i行第j列的元素。进一步的，所述步骤S7中，对构建的全连接加权的脑网络进行频繁子图挖掘，方法如下：不确定图的频繁子图挖掘算法：输入：不确定图数据集和阈值minsup，输出：中的全部频繁子图模式集R，首先，初始化结果集R为空，遍历得到中所有的1-子图模式g，新建空集M，接着，对中的每个1-子图模式g，调用算法1，计算期望支持度如果子图的期望支持度就把它加入到集合M,R中，然后，对M中的每一个1-子图模式，调用算法2，得到新的子图模式集合N，并使用算法1计算新子图模式的期望支持度F，如果新子图模式的期望支持度F≥minsup，就把它加入到集合R中；最后，重复上步，直到没有满足条件的子图生成为止；其中，算法1：计算期望支持度输入：子图模式g、不确定图数据集输出：g在中的期望支持度；首先，对中的每个不确定图按公式(6)计算g在中出现的概率然后，按公式(7)计算g在中的期望支持度并输出算法2：模式增长输入：k-子图模式，1-子图模式集合M；输出：(k+1)-子图模式结果集N；首先，将k-子图模式所在子搜索空间标记为i；然后，在1-子图模式集合M中找出标号>i的子图，将其添加到k-子图上，得到(k+1)-子图，如果(k+1)-子图模式连通，则将其加入结果集N；计算g在中出现的概率公式如下：公式(6)中，表示不确定图包含一个子图g的概率；表示一个确定图G被不确定图所蕴含的概率；是示性函数，如果则反之，若则计算g在中的期望支持度公式如下：公式(7)中，表示g在不确定图集中的期望支持度；表示子图模式g在不确定图中的概率；n表示不确定图个数。进一步的，所述步骤S8中，统计指标偏度的计算公式具体表示如下：公式(8)中，表示的偏度；n+表示正样本的数目；表示子图模式g在确定图正样本中出现的次数；表示子图模式g在不确定图中的概率；表示子图模式g在不确定图集中概率的平均值；判别性函数G-test的计算公式具体表示如下：公式(9)中，n+和n_分别表示正负样本的数目；与分别表示子图模式g在确定图正负样本中出现的次数；在不确定图上计算时，需要将公式中的与替换成与进一步的，所述步骤S9中，构建分类器的具体步骤包括：采用线性核函数，选择最具有判别性的子图特征出现的概率值作为分类特征，构建分类器；然后采用交叉验证方法进行检验，具体步骤为：从样本集中随机选择90％的样本作为训练集，剩余10％的样本作为测试集本文档来自技高网...

【技术保护点】
1.基于独立成分高序不确定脑网络的磁共振影像分类方法，其特征在于，按照以下步骤进行：步骤S1：对静息态功能磁共振成像数据进行预处理，然后使用独立成分分析方法提取出各独立成分；步骤S2：筛选出属于默认网络的独立成分，提取时间序列，并进行后处理；步骤S3：选定长度固定的滑动窗口，按照1s的移动步长，划分默认网络独立成分的时间序列；步骤S4：计算每个时间窗下默认网络成分的时间序列两两间的皮尔逊相关系数，得到低阶关联矩阵；步骤S5：提取每个时间窗下低阶关联矩阵中的元素值，生成高阶关联矩阵；步骤S6：构建高序不确定功能脑网络，节点之间全连接，节点之间边的权值即高阶关联矩阵中的元素；步骤S7：对构建的全连接加权的脑网络进行频繁子图挖掘，得到加权的频繁子图模式；步骤S8：计算基于动态规划的判别性函数的相关统计指标；步骤S9：采用支持向量机分类算法，选择判别性子图特征作为分类依据，构建分类器，然后采用交叉验证方法进行检验；步骤S10：采用互信息分析方法，量化所选判别性子图特征在分类器中的重要度和冗余度，然后根据量化结果对所选判别子图进行二次筛选，优化特征子集。

【技术特征摘要】
1.基于独立成分高序不确定脑网络的磁共振影像分类方法，其特征在于，按照以下步骤进行：步骤S1：对静息态功能磁共振成像数据进行预处理，然后使用独立成分分析方法提取出各独立成分；步骤S2：筛选出属于默认网络的独立成分，提取时间序列，并进行后处理；步骤S3：选定长度固定的滑动窗口，按照1s的移动步长，划分默认网络独立成分的时间序列；步骤S4：计算每个时间窗下默认网络成分的时间序列两两间的皮尔逊相关系数，得到低阶关联矩阵；步骤S5：提取每个时间窗下低阶关联矩阵中的元素值，生成高阶关联矩阵；步骤S6：构建高序不确定功能脑网络，节点之间全连接，节点之间边的权值即高阶关联矩阵中的元素；步骤S7：对构建的全连接加权的脑网络进行频繁子图挖掘，得到加权的频繁子图模式；步骤S8：计算基于动态规划的判别性函数的相关统计指标；步骤S9：采用支持向量机分类算法，选择判别性子图特征作为分类依据，构建分类器，然后采用交叉验证方法进行检验；步骤S10：采用互信息分析方法，量化所选判别性子图特征在分类器中的重要度和冗余度，然后根据量化结果对所选判别子图进行二次筛选，优化特征子集。2.根据权利要求1所述的基于独立成分高序不确定脑网络的磁共振影像分类方法，其特征在于，所述步骤S1中，预处理步骤包括时间层校正，头部校正，空间标准化，平滑；然后，使用基于负熵最大的独立成分分析方法提取独立成分；其中，负熵的定义为：Nh(Y)＝{E[h(Y)]-E[h(YGauss)]}2；(1)公式(1)中，Y是源信号的随机变量，Nh(Y)为随机变量Y的负熵；YGauss是一与Y具有相同方差的高斯随机变量；E为均值运算；h为非线性函数，这里为使Y具有最大的非高斯性，即找到Y负熵的最大近似值，可通过对E{h(Y)}进行优化来获得，而E{h(Y)}的最适条件可通过下式获得：公式(2)中，Y是源信号的随机变量；X是观测信号的随机变量；W是解混矩阵，WT是W的转置；E为均值运算；h为非线性函数；β是一个恒定值；用牛顿迭代方法解此方程后，可以获得近似牛顿迭代公式，简化后就是独立成分分析方法的迭代公式：W+＝E{Xh(WTX)}-E{h′(WTX)}W；(3)公式(3)中，h′是h的导数；W+是W的新值，W＝W+/||W+||。3.根据权利要求1所述的基于独立成分高序不确定脑网络的磁共振影像分类方法，其特征在于：所述步骤S2中，需要使用默认网络模板匹配的方法对独立成分进行筛选，接着对每个被试默认网络成分的时间序列进行后处理，包括去线性漂移和低通滤波，以去除线性漂移和高频生理噪声的影响。4.根据权利要求1所述的基于独立成分高序不确定脑网络的磁共振影像分类方法，其特征在于：所述步骤S4中，计算每个时间窗下默认网络成分的时间序列两两间的皮尔逊相关系数，得到低阶关联矩阵的具体公式是：公式(4)中，rij表示每个时间窗下关联矩阵中第i行第j列的元素，是默认网络成分成分i与j之间的皮尔逊相关系数；n表示时间点个数；xi(t)表示第i个默认网络成分的时间序列；表示第i个默认网络成分的平均时间序列。5.根据权利要求1所述的基于独立成分高序不确定脑网络的磁共振影像分类方法，其特征在于：所述步骤S6中，在高阶关联矩阵上，将其边的权值按照如下公式进行处理，得到的即为高序不确定功能脑网络：公式(5)中，bij表示高序不确定功能脑网络模型中第i行第j列的元素；rij表示高阶关联矩阵中第i行第j列的元素。6.根据权利要求1所述的基于独立成分高序不确定脑网络的磁共振影像分类方法，其特征在于：所述步骤S7中，...

【专利技术属性】
技术研发人员：郭浩，郭涛，程忱，雷波，王恁，李瑶，李欣芸，孙超，
申请(专利权)人：太原理工大学，
类型：发明
国别省市：山西,14