一种基于多数据中心的多模态脑网络特征分类方法技术

本发明专利技术公开了一种基于多数据中心的多模态脑网络特征分类方法，对多个不同数据中心的各模态脑影像数据进行预处理；提取预处理图像中的参数，构建各脑网络；将每个脑网络的上三角元素拼接并聚合为各连接聚合矩阵；将多个数据中心设定为目标域和源域；使用基于超图的多源领域自适应算法将源域中的样本分布转换为与目标域一致，并将其与目标域中的样本组合；使用基于局部有序结构的多任务特征选择算法进行特征选择；对选择出的判别特征采用基于多核学习的支持向量机算法进行分类。该方法降低了多中心脑网络特征分布的异质性；使用多核学习融合功能网络和结构网络以更好地辅助脑网络特征分类，对研究大脑的认知功能和结构的障碍具有重要的应用价值。碍具有重要的应用价值。碍具有重要的应用价值。

【技术实现步骤摘要】
一种基于多数据中心的多模态脑网络特征分类方法


[0001]本专利技术涉及生物医学影像信息处理，具体是涉及一种基于多数据中心的多模态脑网络特征分类方法。

技术介绍

[0002]人类大脑是自然界存在的一种极具复杂的系统，各类神经元细胞通过突触连接在一起，形成了一个非常复杂的网络，因此脑网络是大脑进行各种生理和认知活动的结构基础。随着科学的进步，脑成像技术的发展为脑网络研究提供了广阔的平台。目前人们主要利用功能性磁共振成像技术(Functional magnetic resonance imaging,fMRI)和结构磁共振成像技术(Structural magnetic resonance imaging,sMRI)技术分别构建功能网络和结构网络对大脑进行研究。fMRI具有较高的空间分辨率和非侵入式等特点，而sMRI通过人脑内部的质子密度的分布现象对大脑组织结构部位进行成像，从而获取大脑组织的形态学数据。两种技术在脑网络研究和脑疾病辅助诊断中应用非常广泛，对于人类脑科学研究的发展有着十分重要的意义。
[0003]功能网络是指将不同的神经元集合或者是将某些脑区之间的功能利用一种相关性关系进行整合，并通过功能连接来记录不同脑区之间信号的动态协调性。结构网络是根据不同脑区之间的形态学相关性寻找脑区之间的关联，并且结构连接是大脑中网络连接的物质基础，神经元的各种活动通过结构网络传导至全脑，产生协同作用。研究功能网络可以探索人脑的功能组织模式，有助于理解人脑在正常或疾病状态下的功能活动规律；研究结构网络有助于揭示人脑结构网络的拓扑性质。对两种脑网络进行融合研究可以进一步加速人类探索大脑的脚步，可以更全面地分析和理解脑疾病地病理生理机制。
[0004]近年来，共享数据呈现出增长的趋势，世界各地先后公开发布了多个医学数据库来解决不同的健康问题。尽管多中心医学影像数据库的发展在一定程度上方便了对数据访问权限的获取，但是对多中心数据进行分析是充满挑战性的。多中心数据来源复杂，大多数方法都假设多中心数据来自相同的分布，即现有的方法忽略了不同中心之间扫描设备或成像参数间的差异所造成的数据异构问题，所以无法很好地处理多中心异构数据，这不仅降低了数据的统计能力，也会削弱多中心数据分析研究的可信度。
[0005]现有的多中心异构数据研究可大致分为两类，一类是为每个中心分别学习模型并最终整合，但是这类方法忽略了多中心数据间的异质性；另一类是多源领域自适应问题，即首先基于不同的源域数据训练模型，然后在目标域上测试所学习的模型。无监督低秩表示被提出用于领域自适应，源域数据通过低秩表示变换到目标域以进行后续的学习任务，但是现有的方法忽略了不同源域中数据之间的复杂关系，即高阶关系，而这种关系也许包含有用的信息来提高分类或回归的性能。
[0006]在机器学习领域，多模态学习已经引起了研究人员广泛兴趣。在医学影像领域研究人员开始利用多模态数据进行疾病诊断和分析，整合多种模态数据能够提供必要的互补信息，从而帮助提高脑疾病的诊断精度。特征选择经常被用于去除多模态数据中多余或不
相关的特征，然而由于大脑和脑疾病的复杂性，从单模态数据中检测出所有与疾病相关的特征是非常困难的。由于不同模态数据反映相同的基础病理，也许能互相提供互补的信息有助于寻找异常的脑区。目前已有相关研究探索从多模态数据中联合地选择特征，但是由于患有相同脑疾病的不同被试具有异质性，仅考虑数据和被试标签之间的关系用于特征选择并不能很好地挖掘数据中复杂的内在关系，并且在数据采集过程中往往同时使用多个成像设备采集数据，使得不同类别的样本在某些模态下可能会非常相似，通过一般方法提取出的特征容易重复或相似，判别性较低。

技术实现思路

[0007]专利技术目的：针对以上缺点，本专利技术提供一种降低数据异质性的基于多数据中心的多模态脑网络特征分类方法。
[0008]技术方案：为解决上述问题，本专利技术采用一种基于多数据中心的多模态脑网络特征分类方法，包括以下步骤：
[0009](1)采集多个不同数据中心的各模态脑影像数据，并分别进行预处理，得到各数据中心各模态的预处理图像；
[0010](2)将得到的预处理图像基于标准化大脑分区模板进行脑区划分，均划分为N个脑区，每个脑区分别对应脑网络中的一个节点；
[0011](3)分别提取各数据中心各模态预处理图像中的参数，分别构建各模态的脑网络；
[0012](4)将每个脑网络的上三角元素拼接成一个维数为的连接向量；并将每个数据中心的相同模态的所有连接向量组合成连接聚合矩阵；
[0013](5)选择一个数据中心为目标域，其他数据中心为源域，对目标域和源域中连接聚合矩阵的连接向量分别划分训练集和测试集；使用基于超图的多源领域自适应算法为每个源域中的训练集求解得到一个转换矩阵，通过转换矩阵将多个源域中的连接向量分布转换为与目标域中的连接向量分布一致，并将转换后源域中的训练集与目标域中的训练集组合为最终的各模态连接训练集；
[0014](6)使用基于局部有序结构的多任务特征选择算法在最终的各模态连接训练集中训练模型和选择特征；基于训练集训练的模型以及选择的特征，在目标域测试集中选出判别特征；
[0015](7)对选择出的判别特征采用基于多核学习的支持向量机算法进行分类，最终得到分类结果。
[0016]进一步的，所述步骤(1)中采集多个不同数据中心的大脑功能磁共振成像和大脑结构磁共振成像；对大脑结构核磁共振成像预处理后得到大脑灰质图像。
[0017]所述步骤(2)中标准化大脑分区模板包括AAL分区模板和Brodmann分区模板。
[0018]所述步骤(3)中各模态的脑网络包括大脑功能网络和大脑结构网络，对于不同数据中心的大脑功能磁共振图像，提取不同脑区的平均时间序列，计算两两脑区平均时间序列之间的相关系数，构建每个被试的大脑功能网络；对于不同数据中心的大脑灰质图像，使用核密度估计方法估计每个脑区所有体素的灰质体积值的概率密度函数，然后计算两两脑区灰质体积的概率密度函数之间的Jensen散度，将Jensen散度定义为结构网络的连接权重
以评估不同脑区灰质体积的相似性，从而构建每个被试的大脑结构网络。
[0019]所述Jensen散度的计算公式为：
[0020][0021]式中，和分别表示脑区P和脑区Q灰质体积的概率分布，和对应的概率密度函数为和D
KL
(.||.)为计算两两脑区灰质体积的概率密度函数之间的KL散度；(.||.)为计算两两脑区灰质体积的概率密度函数之间的KL散度；和之间KL散度的计算公式为：
[0022][0023]进一步的，所述步骤(5)中基于超图的多源领域自适应算法中，第m个模态的目标函数如下所示：
[0024][0025]式中，||.||
*
为矩阵的核范数，||.||
2,1
为矩阵的l
2,1
范数，tr(.)为矩阵的迹运算，为目标域训练集的特征矩阵，D
m本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于多数据中心的多模态脑网络特征分类方法，其特征在于，包括以下步骤：(1)采集多个不同数据中心的各模态脑影像数据，并分别进行预处理，得到各数据中心各模态的预处理图像；(2)将得到的预处理图像基于标准化大脑分区模板进行脑区划分，均划分为N个脑区，每个脑区分别对应脑网络中的一个节点；(3)分别提取各数据中心各模态预处理图像中的参数，分别构建各模态的脑网络；(4)将每个脑网络的上三角元素拼接成一个维数为的连接向量；并将每个数据中心的相同模态的所有连接向量组合成连接聚合矩阵；(5)选择一个数据中心为目标域，其他数据中心为源域，对目标域和源域中连接聚合矩阵的连接向量分别划分训练集和测试集；使用基于超图的多源领域自适应算法为每个源域中的训练集求解得到一个转换矩阵，通过转换矩阵将多个源域中的连接向量分布转换为与目标域中的连接向量分布一致，并将转换后源域中的训练集与目标域中的训练集组合为最终的各模态连接训练集；(6)使用基于局部有序结构的多任务特征选择算法在最终的各模态连接训练集中训练模型和选择特征；基于训练集训练的模型以及选择的特征，在目标域测试集中选出判别特征；(7)对选择出的判别特征采用基于多核学习的支持向量机算法进行分类，最终得到分类结果。2.根据权利要求1所述的多模态脑网络特征分类方法，其特征在于，所述步骤(1)中采集多个不同数据中心的大脑功能磁共振成像和大脑结构磁共振成像；对大脑结构核磁共振成像预处理后得到大脑灰质图像。3.根据权利要求1所述的多模态脑网络特征分类方法，其特征在于，所述步骤(2)中标准化大脑分区模板包括AAL分区模板和Brodmann分区模板。4.根据权利要求2所述的多模态脑网络特征分类方法，其特征在于，所述步骤(3)中各模态的脑网络包括大脑功能网络和大脑结构网络，对于不同数据中心的大脑功能磁共振图像，提取不同脑区的平均时间序列，计算两两脑区平均时间序列之间的相关系数，构建每个被试的大脑功能网络；对于不同数据中心的大脑灰质图像，使用核密度估计方法估计每个脑区所有体素的灰质体积值的概率密度函数，然后计算两两脑区灰质体积的概率密度函数之间的Jensen散度，将Jensen散度定义为结构网络的连接权重以评估不同脑区灰质体积的相似性，从而构建每个被试的大脑结构网络。5.根据权利要求4所述的多模态脑网络特征分类方法，其特征在于，所述Jensen散度的计算公式为：式中，和分别表示脑区P和脑区Q灰质体积的概率分布，和对应的概率密度函数为和D
KL
(.||.)为计算两两脑区灰质体积的概率密度函数之间的KL散度；(.||.)为计算两两脑区灰质体积的概率密度函数之间的KL散度；和之间KL散度的计算公式为：
6.根据权利要求1所述的多模态脑网络特征分类方法，其特征在于，所述步骤(5)中基于超图的多源领域自适应算法中，第m个模态的目标函数如下所示：于超图的多源领域自适应算法中，第m个模态的目标函数如下所示：式中，||.||
*
为矩阵的核范数，||.||
2,1
为矩阵的l
2,1
范数，tr(.)为矩阵的迹运算，为目标域训练集...

【专利技术属性】
技术研发人员：戚世乐，季一新，张道强，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：