基于多尺度脑网络特征的功能磁共振影像数据分类方法技术

本发明专利技术涉及图像处理技术，具体是一种基于多尺度脑网络特征的功能磁共振影像数据分类方法。本发明专利技术解决了传统磁共振影像数据分类方法分类准确率低的问题。基于多尺度脑网络特征的功能磁共振影像数据分类方法，该方法是采用如下步骤实现的：步骤S1：对静息态功能磁共振影像进行预处理；步骤S2：采用动态随机种子方法对影像进行区域分割，并对所分割的脑区进行平均时间序列提取；步骤S3：计算各脑区平均时间序列两两间的关联程度；步骤S4：对关联矩阵进行二值化处理；步骤S5：计算静息态功能脑网络的局部属性及其在特定阈值空间的AUC值；步骤S6：进行分类器构建；步骤S7：量化所选特征在分类器中的重要度和冗余度。本发明专利技术适用于磁共振影像数据分类。

【技术实现步骤摘要】

本专利技术涉及图像处理技术，具体是一种基于多尺度脑网络特征的功能磁共振影像数据分类方法。
技术介绍
作为功能磁共振成像(functionalMagneticResonanceImaging,fMRI)技术与复杂网络理论的结合，磁共振影像数据分类方法当前已经成为脑科学领域的热点之一。然而，传统磁共振影像数据分类方法由于自身原理所限，仅能够对单一尺度的脑网络进行描述，由此导致其分类准确率低，从而严重影响其应用价值。基于此，有必要专利技术一种全新的磁共振影像数据分类方法，以解决传统磁共振影像数据分类方法存在的上述问题。
技术实现思路
本专利技术为了解决传统磁共振影像数据分类方法分类准确率低的问题，提供了一种基于多尺度脑网络特征的功能磁共振影像数据分类方法。本专利技术是采用如下技术方案实现的：基于多尺度脑网络特征的功能磁共振影像数据分类方法，该方法是采用如下步骤实现的：步骤S1：对静息态功能磁共振影像进行预处理；步骤S2：根据选定的标准化脑图谱，采用动态随机种子方法对预处理后的静息态功能磁共振影像进行区域分割，分割尺度分别为90，256，497，1003，1501，然后对所分割的各脑区进行平均时间序列的提取；步骤S3：采用皮尔逊相关方法，计算各脑区的平均时间序列两两间的关联程度，由此得到关联矩阵；步骤S4：设定阈值，然后根据阈值对关联矩阵进行二值化处理，由此得到静息态功能脑网络模型；步骤S5：计算静息态功能脑网络模型的局部属性及局部属性在特定阈值空间内的AUC值；所述局部属性包括：静息态功能脑网络模型中各节点的度、节点效率、中间中心度；步骤S6：采用支持向量机分类算法，选择静息态功能脑网络模型的局部属性作为分类特征，由此进行分类器的构建，然后采用交叉验证方法对构建的分类器进行检验；步骤S7：采用互信息分析方法，对所选特征在分类器中的重要度和冗余度进行量化，然后根据量化结果对所选特征进行二次筛选，由此对静息态功能脑网络模型进行优化。与传统磁共振影像数据分类方法相比，本专利技术所述的基于多尺度脑网络特征的功能磁共振影像数据分类方法通过采用动态随机种子方法、皮尔逊相关方法、支持向量机分类算法、交叉验证方法、互信息分析方法，实现了对多尺度的脑网络进行描述，由此大幅提高了分类准确率(如图1所示，本专利技术的分类准确率明显高于传统磁共振影像数据分类方法的分类准确率)，从而使得应用价值更高。本专利技术有效解决了传统磁共振影像数据分类方法分类准确率低的问题，适用于磁共振影像数据分类。附图说明图1是本专利技术与传统磁共振影像数据分类方法的对比示意图。具体实施方式基于多尺度脑网络特征的功能磁共振影像数据分类方法，该方法是采用如下步骤实现的：步骤S1：对静息态功能磁共振影像进行预处理；步骤S2：根据选定的标准化脑图谱，采用动态随机种子方法对预处理后的静息态功能磁共振影像进行区域分割，分割尺度分别为90，256，497，1003，1501，然后对所分割的各脑区进行平均时间序列的提取；步骤S3：采用皮尔逊相关方法，计算各脑区的平均时间序列两两间的关联程度，由此得到关联矩阵；步骤S4：设定阈值，然后根据阈值对关联矩阵进行二值化处理，由此得到静息态功能脑网络模型；步骤S5：计算静息态功能脑网络模型的局部属性及局部属性在特定阈值空间内的AUC值；所述局部属性包括：静息态功能脑网络模型中各节点的度、节点效率、中间中心度；步骤S6：采用支持向量机分类算法，选择静息态功能脑网络模型的局部属性作为分类特征，由此进行分类器的构建，然后采用交叉验证方法对构建的分类器进行检验；步骤S7：采用互信息分析方法，对所选特征在分类器中的重要度和冗余度进行量化，然后根据量化结果对所选特征进行二次筛选，由此对静息态功能脑网络模型进行优化。所述步骤S1中，预处理采用SPM软件进行，预处理步骤具体包括：时间层校正、头动校正、联合配准、空间标准化、低频滤波。所述步骤S2中，标准化脑图谱采用AAL模板；区域分割步骤具体包括：首先，计算AAL模板中每个脑区占所有脑区的体素比例V；然后，计算AAL模板中原有脑区在节点规模N下可细化的子区域个数k，k＝VN；然后，对相应脑区设置k个随机种子体素，并依次计算所有剩余体素与种子体素的距离；然后，采用动态随机种子点的设置方法，将当前体素与距离最近的体素组合形成新的子区域，并将新的物理中心设置为新的种子体素；依次循环，直至脑区内所有体素均分割完成；平均时间序列的提取步骤具体包括：提取AAL模板中每个脑区所包含的所有体素在不同时间点上的BOLD强度，并将各体素在不同时间点上的BOLD强度进行算术平均，由此得到各脑区的平均时间序列。所述步骤S3中，计算公式具体表示如下：rij=Σi,j=1n(xi(t)-xi‾)(xj(t)-xj‾)Σi,j=1n(xi(t)-xi‾)2Σi,j=1n(xj(t)-xj‾)2---(1);]]>公式(1)中：rij表示关联矩阵中第i行第j列的元素；n表示时间点个数；xi(t)表示第i个脑区的时间序列；表示第i个脑区的时间序列的平均值；xj(t)表示第j个脑区的时间序列；表示第j个脑区的时间序列的平均值；关联矩阵的维度分别为90×90，256×256，497×497，1003×1003，1501×1501。所述步骤S4中，二值化处理公式具体表示如下：bij=1,rij≥τ0,rij<τ---(2);]]>公式(2)中：bij表示静息态功能脑网络模型中第i行第j列的元素；rij表示关联矩阵中第i行第j列的元素；τ表示阈值；静息态功能脑网络模型的维度分别为90×90，256×256，497×497，1003×1003，1501×1501。所述步骤S5中，计算公式具体表示如下：ki=Σj∈Gaij---(3);]]>公式(3)中：ki表示静息态功能脑网络模型中任意一节点i的度；aij表示静息态功能脑网络模型中节点i与节点j之间的连接；ei=1N-1Σj≠i∈G1dij---(4);]]>公式(4)中：ei表示静息态功能脑网络模型中任意一节点i的节点效率；dij表示静息态功能脑网络模型中节点i与节点j之间的最短路径长度；bi=Σm≠i≠n∈Gσmn(i)σmn---(5);]]>公式(5)中：bi表示静息态功能脑网络模型中任意一节点i的中间中心度；σmn表示从节点m到节点n的最短路径的数量；σmn(i)表示从节点m到节点n经过节点i的最短路径的数量；YAUC=Σk=1n-1[Y(Sk)+Y(Sk-1)]×ΔS/2---(6);]]>公式(6)中：YAUC表示各局部属性值在特定阈值空间内的AUC值；Y(Sk)表示阈值Sk对应的局部属性值；Y(Sk-1)表示阈值Sk-1对应的局部属性值；ΔS表示两个阈值之间的间隔。所述步骤S6中，分类器的构建步骤具体包括：采用RBF核函数，选择双样本文档来自技高网...

【技术保护点】
一种基于多尺度脑网络特征的功能磁共振影像数据分类方法，其特征在于：该方法是采用如下步骤实现的：步骤S1：对静息态功能磁共振影像进行预处理；步骤S2：根据选定的标准化脑图谱，采用动态随机种子方法对预处理后的静息态功能磁共振影像进行区域分割，分割尺度分别为90，256，497，1003，1501，然后对所分割的各脑区进行平均时间序列的提取；步骤S3：采用皮尔逊相关方法，计算各脑区的平均时间序列两两间的关联程度，由此得到关联矩阵；步骤S4：设定阈值，然后根据阈值对关联矩阵进行二值化处理，由此得到静息态功能脑网络模型；步骤S5：计算静息态功能脑网络模型的局部属性及局部属性在特定阈值空间内的AUC值；所述局部属性包括：静息态功能脑网络模型中各节点的度、节点效率、中间中心度；步骤S6：采用支持向量机分类算法，选择静息态功能脑网络模型的局部属性作为分类特征，由此进行分类器的构建，然后采用交叉验证方法对构建的分类器进行检验；步骤S7：采用互信息分析方法，对所选特征在分类器中的重要度和冗余度进行量化，然后根据量化结果对所选特征进行二次筛选，由此对静息态功能脑网络模型进行优化。

【技术特征摘要】
1.一种基于多尺度脑网络特征的功能磁共振影像数据分类方法，其特征在于：该方法是采用如下步骤实现的：步骤S1：对静息态功能磁共振影像进行预处理；步骤S2：根据选定的标准化脑图谱，采用动态随机种子方法对预处理后的静息态功能磁共振影像进行区域分割，分割尺度分别为90，256，497，1003，1501，然后对所分割的各脑区进行平均时间序列的提取；步骤S3：采用皮尔逊相关方法，计算各脑区的平均时间序列两两间的关联程度，由此得到关联矩阵；步骤S4：设定阈值，然后根据阈值对关联矩阵进行二值化处理，由此得到静息态功能脑网络模型；步骤S5：计算静息态功能脑网络模型的局部属性及局部属性在特定阈值空间内的AUC值；所述局部属性包括：静息态功能脑网络模型中各节点的度、节点效率、中间中心度；步骤S6：采用支持向量机分类算法，选择静息态功能脑网络模型的局部属性作为分类特征，由此进行分类器的构建，然后采用交叉验证方法对构建的分类器进行检验；步骤S7：采用互信息分析方法，对所选特征在分类器中的重要度和冗余度进行量化，然后根据量化结果对所选特征进行二次筛选，由此对静息态功能脑网络模型进行优化。2.根据权利要求1所述的基于多尺度脑网络特征的功能磁共振影像数据分类方法，其特征在于：所述步骤S1中，预处理采用SPM软件进行，预处理步骤具体包括：时间层校正、头动校正、联合配准、空间标准化、低频滤波。3.根据权利要求1所述的基于多尺度脑网络特征的功能磁共振影像数据分类方法，其特征在于：所述步骤S2中，标准化脑图谱采用AAL模板；区域分割步骤具体包括：首先，计算AAL模板中每个脑区占所有脑区的体素比例V；然后，计算AAL模板中原有脑区在节点规模N下可细化的子区域个数k，k＝VN；然后，对相应脑区设置k个随机种子体素，并依次计算所有剩余体素与种子体素的距离；然后，采用动态随机种子点的设置方法，将当前体素与距离最近的体素组合形成新的子区域，并将新的物理中心设置为新的种子体素；依次循环，直至脑区内所有体素均分割完成；平均时间序列的提取步骤具体包括：提取AAL模板中每个脑区所包含的所有体素在不同时间点上的BOLD强度，并将各体素在不同时间点上的BOLD强度进行算术平均，由此得到各脑区的平均时间序列。4.根据权利要求1所述的基于多尺度脑网络特征的功能磁共振影像数据分类方法，其特征在于：所述步骤S3中，计算公式具体表示如下：rij=Σi,j=1n(xi(t)-xi‾)(xj(t)-xj‾)Σi=1n(xi(t)-xi‾)2Σj=1n(xj(t)-xj‾)2---(1);]]>公式(1)中：rij表示关联矩阵中第i行第j列的元素；n表示时间点个数；xi(t)表示第i个脑区的时间序列；表示第i个脑区的时间序列的平均值；xj(t)表示第j个脑区的时间序列；表示第j个脑区的时间序列的平均值；关联矩阵的维度分别为90×90，256×256，497×497，1003×1003，1501×1501。5.根据权利要求1所述的基于多尺度脑网络特征的功能磁共振影像数据分类方法，其特征在于：所述步骤S4中，二值化处理公式具体表示如下：bij=1,rij&Greater...

【专利技术属性】
技术研发人员：郭浩，曹锐，王彬，相洁，李海芳，陈俊杰，
申请(专利权)人：太原理工大学，
类型：发明
国别省市：山西;14