基于最小生成树的多模态磁共振影像数据分类方法技术

本发明专利技术涉及图像处理技术，具体是一种基于最小生成树的多模态磁共振影像数据分类方法。本发明专利技术解决了传统磁共振影像数据分类方法分类准确率低的问题。基于最小生成树的多模态磁共振影像数据分类方法，该方法是采用如下步骤实现的：步骤S1：对静息态功能磁共振影像进行预处理；步骤S2：计算各脑区的平均时间序列两两之间的皮尔逊相关系数；步骤S3：对磁共振弥散张量成像影像进行预处理；步骤S4：计算各脑区两两之间的权重；步骤S5：构建多模态最小生成树脑网络；步骤S6：计算多模态最小生成树脑网络的局部属性；步骤S7：进行分类器的构建；步骤S8：对所选特征在分类器中的重要度和冗余度进行量化。本发明专利技术适用于磁共振影像数据分类。

【技术实现步骤摘要】

本专利技术涉及图像处理技术，具体是一种基于最小生成树的多模态磁共振影像数据分类方法。
技术介绍
作为多模态磁共振成像技术与复杂网络理论的结合，磁共振影像数据分类方法当前已经成为脑科学领域的热点之一。其广泛应用在各类研究中，特别是脑疾病的研究，并取得许多令人惊喜的成果。功能磁共振(functionalmagneticresonanceimaging，fMRI)因其无创、高分辨率及可重复性等优点，在对以功能改变为特征的疾病研究中具有良好的应用前景。功能磁共振是在较短的时间内连续抓取一个时间序列内相同部位的成像，时间分辨率较高。但由于功能磁共振技术的天生弱点，相应的代价就是空间分辨率会下降。结构磁共振成像(structuralmagneticresonanceimaging，sMRI)能够客观记录下从疾病潜伏期到发作期整个过程中患者脑结构生物标记的变化，这些数据能够从根本上改变人们对这种疾病的认识，并且能够影响和引导疾病的后续诊断和治疗。结构磁共振成像可以清楚的分辨病人的病灶位置，空间分辨率高，可以弥补功能磁共振成像的缺点。因此，本专利技术采用基于功能磁共振和结构磁共振的多模态磁共振成像技术，来提高测量结果的准确性。传统影像数据分类方法由于其自身原理所限，普遍存在方法论的限制。比如，图论方法受网络大小(即节点的数目)、网络稀疏度(即存在连接的百分比)和平均度(即每个节点的连接数)的影响。而在脑网络的构建过程中，阈值、冗余链接以及虚假链接会影响脑网络的拓扑结构。由此导致其分类准确率低，从而严重影响其应用价值。虽然一些技术可以在一定程度上控制上述因素所带来的误差，却很难做到完全避免这些因素对结果的影响。基于此，有必要专利技术一种全新的磁共振影像数据分类方法，以解决传统磁共振影像数据分类方法存在的上述问题。最小生成树是一个数学定义的、无偏的子网，计算结果不会因为网络大小、平均度或密度的影响而出现误差，可以避免一般的方法论问题。因此，本专利技术采用了一种基于最小生成树的多模态磁共振影像数据分类方法，以提高分类结果的准确性。
技术实现思路
本专利技术为了解决传统磁共振影像数据分类方法分类准确率低的问题，提供了一种基于最小生成树的多模态磁共振影像数据分类方法。本专利技术是采用如下技术方案实现的：基于最小生成树的多模态磁共振影像数据分类方法，该方法是采用如下步骤实现的：步骤S1：对静息态功能磁共振影像进行预处理，并根据选定的标准化脑图谱对预处理后的静息态功能磁共振影像进行区域分割，然后对所分割的各脑区进行平均时间序列的提取；步骤S2：计算各脑区的平均时间序列两两之间的协方差系数，由此得到协方差矩阵；然后根据协方差矩阵计算各脑区的平均时间序列两两之间的皮尔逊相关系数，由此得到皮尔逊相关矩阵；步骤S3：对磁共振弥散张量成像影像进行预处理，并根据选定的标准化脑图谱对预处理后的磁共振弥散张量成像影像进行区域分割，然后对所分割的各脑区进行全脑神经纤维追踪，由此得到各脑区的平均局部各向异性及各脑区之间的神经纤维数量；步骤S4：计算各脑区两两之间的权重，由此得到神经纤维矩阵；步骤S5：结合皮尔逊相关矩阵和神经纤维矩阵，利用克鲁斯卡尔算法构建多模态最小生成树脑网络；步骤S6：计算多模态最小生成树脑网络的局部属性；所述局部属性包括：多模态最小生成树脑网络中各节点的度、介数中心度、离心率；步骤S7：采用支持向量机分类算法，选择多模态最小生成树脑网络的局部属性作为分类特征，由此进行分类器的构建，然后采用交叉验证方法对构建的分类器进行检验；步骤S8：采用互信息分析方法，对所选特征在分类器中的重要度和冗余度进行量化，然后根据量化结果对所选特征进行二次筛选，由此对多模态最小生成树脑网络进行优化。与传统磁共振影像数据分类方法相比，本专利技术所述的基于最小生成树的多模态磁共振影像数据分类方法通过采用皮尔逊相关方法、克鲁斯卡尔算法、支持向量机分类算法、交叉验证方法、互信息分析方法，实现了多模态最小生成树脑网络的描述，由此大幅提高了分类准确率(如图1所示，本专利技术的分类准确率明显高于传统磁共振影像数据分类方法的分类准确率)，从而使得应用价值更高。本专利技术有效解决了传统磁共振影像数据分类方法分类准确率低的问题，适用于磁共振影像数据分类。附图说明图1是本专利技术与传统磁共振影像数据分类方法的对比示意图。具体实施方式基于最小生成树的多模态磁共振影像数据分类方法，该方法是采用如下步骤实现的：步骤S1：对静息态功能磁共振影像进行预处理，并根据选定的标准化脑图谱对预处理后的静息态功能磁共振影像进行区域分割，然后对所分割的各脑区进行平均时间序列的提取；步骤S2：计算各脑区的平均时间序列两两之间的协方差系数，由此得到协方差矩阵；然后根据协方差矩阵计算各脑区的平均时间序列两两之间的皮尔逊相关系数，由此得到皮尔逊相关矩阵；步骤S3：对磁共振弥散张量成像影像进行预处理，并根据选定的标准化脑图谱对预处理后的磁共振弥散张量成像影像进行区域分割，然后对所分割的各脑区进行全脑神经纤维追踪，由此得到各脑区的平均局部各向异性及各脑区之间的神经纤维数量；步骤S4：计算各脑区两两之间的权重，由此得到神经纤维矩阵；步骤S5：结合皮尔逊相关矩阵和神经纤维矩阵，利用克鲁斯卡尔算法构建多模态最小生成树脑网络；步骤S6：计算多模态最小生成树脑网络的局部属性；所述局部属性包括：多模态最小生成树脑网络中各节点的度、介数中心度、离心率；步骤S7：采用支持向量机分类算法，选择多模态最小生成树脑网络的局部属性作为分类特征，由此进行分类器的构建，然后采用交叉验证方法对构建的分类器进行检验；步骤S8：采用互信息分析方法，对所选特征在分类器中的重要度和冗余度进行量化，然后根据量化结果对所选特征进行二次筛选，由此对多模态最小生成树脑网络进行优化。所述步骤S1中，预处理步骤具体包括：时间层校正、头动校正、联合配准、空间标准化、低频滤波；标准化脑图谱采用AAL模板；平均时间序列的提取步骤具体包括：提取AAL模板中每个脑区所包含的所有体素在不同时间点上的BOLD强度，并将各体素在不同时间点上的BOLD强度进行算术平均，由此得到各脑区的平均时间序列。所述步骤S2中，计算公式具体表示如下：公式(1)中：Sij表示协方差矩阵中第i行第j列的元素，即第i个脑区和第j个脑区之间的协方差系数；M表示时间点个数；xi(t)表示第i个脑区的时间序列；表示第i个脑区的时间序列的平均值；xj(t)表示第j个脑区的时间序列；表示第j个脑区的时间序列的平均值；协方差矩阵的维度为90×90；公式(2)中：Rij表示皮尔逊相关矩阵中第i行第j列的元素，即第i个脑区和第j个脑区之间的皮尔逊相关系数；xi表示第i个脑区的协方差序列；yj表示第j个脑区的协方差序列；表示第i个脑区和第j个脑区的协方差；表示第i个脑区的协方差序列的标准差；表示第j个脑区的协方差序列的标准差；皮尔逊相关矩阵的维度为90×90。所述步骤S3中，预处理采用PANDA软件进行；预处理步骤具体包括：配准、标准化、涡流矫正；标准化脑图谱采用AAL模板；全脑神经纤维追踪采用FACT算法进行。所述步骤S4中，计算公式具体表示如下：公式(3)中：Dij表示神经纤维矩阵中第i行第j列的元素，即第i个脑区和第j个脑区之间的权本文档来自技高网...

【技术保护点】
一种基于最小生成树的多模态磁共振影像数据分类方法，其特征在于：该方法是采用如下步骤实现的：步骤S1：对静息态功能磁共振影像进行预处理，并根据选定的标准化脑图谱对预处理后的静息态功能磁共振影像进行区域分割，然后对所分割的各脑区进行平均时间序列的提取；步骤S2：计算各脑区的平均时间序列两两之间的协方差系数，由此得到协方差矩阵；然后根据协方差矩阵计算各脑区的平均时间序列两两之间的皮尔逊相关系数，由此得到皮尔逊相关矩阵；步骤S3：对磁共振弥散张量成像影像进行预处理，并根据选定的标准化脑图谱对预处理后的磁共振弥散张量成像影像进行区域分割，然后对所分割的各脑区进行全脑神经纤维追踪，由此得到各脑区的平均局部各向异性及各脑区之间的神经纤维数量；步骤S4：计算各脑区两两之间的权重，由此得到神经纤维矩阵；步骤S5：结合皮尔逊相关矩阵和神经纤维矩阵，利用克鲁斯卡尔算法构建多模态最小生成树脑网络；步骤S6：计算多模态最小生成树脑网络的局部属性；所述局部属性包括：多模态最小生成树脑网络中各节点的度、介数中心度、离心率；步骤S7：采用支持向量机分类算法，选择多模态最小生成树脑网络的局部属性作为分类特征，由此进行分类器的构建，然后采用交叉验证方法对构建的分类器进行检验；步骤S8：采用互信息分析方法，对所选特征在分类器中的重要度和冗余度进行量化，然后根据量化结果对所选特征进行二次筛选，由此对多模态最小生成树脑网络进行优化。...

【技术特征摘要】
1.一种基于最小生成树的多模态磁共振影像数据分类方法，其特征在于：该方法是采用如下步骤实现的：步骤S1：对静息态功能磁共振影像进行预处理，并根据选定的标准化脑图谱对预处理后的静息态功能磁共振影像进行区域分割，然后对所分割的各脑区进行平均时间序列的提取；步骤S2：计算各脑区的平均时间序列两两之间的协方差系数，由此得到协方差矩阵；然后根据协方差矩阵计算各脑区的平均时间序列两两之间的皮尔逊相关系数，由此得到皮尔逊相关矩阵；步骤S3：对磁共振弥散张量成像影像进行预处理，并根据选定的标准化脑图谱对预处理后的磁共振弥散张量成像影像进行区域分割，然后对所分割的各脑区进行全脑神经纤维追踪，由此得到各脑区的平均局部各向异性及各脑区之间的神经纤维数量；步骤S4：计算各脑区两两之间的权重，由此得到神经纤维矩阵；步骤S5：结合皮尔逊相关矩阵和神经纤维矩阵，利用克鲁斯卡尔算法构建多模态最小生成树脑网络；步骤S6：计算多模态最小生成树脑网络的局部属性；所述局部属性包括：多模态最小生成树脑网络中各节点的度、介数中心度、离心率；步骤S7：采用支持向量机分类算法，选择多模态最小生成树脑网络的局部属性作为分类特征，由此进行分类器的构建，然后采用交叉验证方法对构建的分类器进行检验；步骤S8：采用互信息分析方法，对所选特征在分类器中的重要度和冗余度进行量化，然后根据量化结果对所选特征进行二次筛选，由此对多模态最小生成树脑网络进行优化。2.根据权利要求1所述的基于最小生成树的多模态磁共振影像数据分类方法，其特征在于：所述步骤S1中，预处理步骤具体包括：时间层校正、头动校正、联合配准、空间标准化、低频滤波；标准化脑图谱采用AAL模板；平均时间序列的提取步骤具体包括：提取AAL模板中每个脑区所包含的所有体素在不同时间点上的BOLD强度，并将各体素在不同时间点上的BOLD强度进行算术平均，由此得到各脑区的平均时间序列。3.根据权利要求1所述的基于最小生成树的多模态磁共振影像数据分类方法，其特征在于：所述步骤S2中，计算公式具体表示如下：Sij=1MΣt=1M(xi(t)-xi‾)(xj(t)-xj‾)---(1);]]>公式(1)中：Sij表示协方差矩阵中第i行第j列的元素，即第i个脑区和第j个脑区之间的协方差系数；M表示时间点个数；xi(t)表示第i个脑区的时间序列；表示第i个脑区的时间序列的平均值；xj(t)表示第j个脑区的时间序列；表示第j个脑区的时间序列的平均值；协方差矩阵的维度为90×90；Rij=Σi,j=1n(xi-x‾)(yj-y‾)Σi=1n(xi-x‾)2Σj=1n(yj-y‾)2---(2);]]>公式(2)中：Rij表示皮尔逊相关矩阵中第i行第j列的元素，即第i个脑区和第j个脑区之间的皮尔逊相关系数；xi表示第i个脑区的协方差序列；yj表示第j个脑区的协方差序列；表示第i个脑区和第j个脑区的协方差；表示第i个脑区的协方差序列的标准差；表示第j个脑区的协方差序列的标准差；皮尔逊相关矩阵的维度为90×90。4.根据权利要求1所述的基于最小生成树的多模态磁共振影像数据分类方法，其特征在于：所述步骤S3中，预处理采用PANDA软件进行；预处理步骤具体包括：配准、标准化、涡流矫正；标准化脑图谱采用AAL模板；全脑神经纤维追踪采用FACT算法进行。5.根据权利要求1所述的基于最小生成树的多模态磁共振影像数据分类方法，其特征在于：所述步骤S4中，计算公式具体表示如下：Dij=FNij×(FAi+FAj2)---(3);]]>公式(3)中：Dij表示神经纤维矩阵中第i行第j列的元素，即第i个脑区和第j个脑区之间的权重；FNij表示第i个脑区和第j个脑区之间的神经纤维数量；FAi表示第i个脑区的平均局部各向异性；FAj表示第j个脑区的平均局部各向异...

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