本发明属于脑影像与计算机科学的交叉领域，涉及图像分类技术领域，尤其涉及一种基于融合网络的阿尔兹海默病分类方法。
背景技术：
随着人工智能和神经影像技术的快速发展，医学图像处理作为与医疗科学技术相关最为密切的领域之一正在蓬勃发展。阿尔兹海默病的分类是医学图像分类领域的一个重要研究方向，在计算机辅助诊断中具有重要意义。通过机器学习对磁共振图像进行分析，能够实现阿尔茨海默症(alzheimer'sdisease，ad)的快速、准确分类。在某些情况下，这种方法表现出比临床医生更好的分类准确率。基于机器学习的计算机辅助诊断发展成为一个重要的研究领域，表现出巨大的优势和发展的前景。
早期ad的分类方法依靠领域知识结合经典的机器学习方法，揭示阿尔兹海默病发生的病理细节。但是这些方法基于人为的先验假设，高度依赖于手工提取的特征，分类的效果受图像预处理步骤，如分割和配准等的影响。近年来，深度学习的成功进一步推动了ad分类的研究。深度学习作为一种先进的机器学习方法，避免了图像预处理步骤带来的负面影响，能够在大规模且复杂的高维数据中提取重要的信息，从而获得数据中隐含的层次化特征和数据的复杂结构。这些利用基于深度学习的基网络实现ad分类的模型，相较于传统的机器学习方法表现出良好的效果，但存在模型泛化性差的问题，以及准确率还有待进一步提升。
技术实现要素：
针对以上问题，本发明提出一种基于融合网络的阿尔兹海默病分类方法，有效使用被试的结构核磁共振图像数据进行自动辅助疾病分类。
本发明的技术方案：一种基于融合网络的阿尔兹海默病分类方法，包括以下步骤：
步骤s1：给定被试的磁共振图像数据集，数据集包含磁共振图像和被试的标签，磁共振图像数据集记为c＝{(ai,bi)|i∈[1,r]}，式中，a表示磁共振图像，b表示被试的对应标签，r表示样本数量，i表示样本下标。将磁共振图像全集c拷贝并划分成磁共振图像训练集e＝{(ai,bi)|i∈[1,m]}和磁共振图像测试集e＝{(ai,bi)|i∈[1,n]}，其中n m＝r，且m和n分别表示划分后的样本数量；
步骤s2：对磁共振图像训练集依次进行数据预处理、水平翻转扩增操作和数据归一化操作，得到处理后的磁共振图像训练集e'＝{(ai,bi)|i∈[1,m']}，经过扩增后，其数量相应的增多至m'；对磁共振图像测试集f依次进行数据预处理、数据归一化操作得到归一化的磁共振图像测试集f'＝{(ai,bi)|i∈[1,n]}；
步骤s3：将磁共振图像训练集e'作为输入，载入到计算机内存中，并将磁共振图像训练集e'＝{(ai,bi)|i∈[1,m']}随机打乱，为接下来的训练阶段做准备；
步骤s4：分别构建基网络g1和g2。在pytorch框架中使用sgd优化方法，对基网络g1和g2的网络参数进行初始化，得到基网络g1和g2的初始参数，并设定相关训练超参数，用于网络模型的训练优化。通过基网络g1和g2分别对训练磁共振图像进行特征提取，得到特征h1和h2，并对分别提取的特征进行特征融合得到融合特征h3；
步骤s5：构建融合特征分类决策层p，将基网络g1和g2提取到的特征h1和h2以及融合特征h3分别输入到各网络的决策层，分别得到分类决策层的标签p1、p2和p3。对p1、p2和p3进行融合，得到对应输入样本最终的分类标签b'；
步骤s6：通过损失函数计算，计算训练数据通过网络分类得到的标签b'和训练真实标签b之间的误差，得出损失值。利用得出的损失值进行反向传播，对基网络g1和g2，以及融合网络分别进行网络参数调整；
步骤s7：利用测试集f'对融合网络模型进行评估，保留网络模型g。同时在参数更新结束阶段进行判断，判断训练迭代次数是否已达到最大迭代次数，若已经达到最大迭代次数，则训练阶段结束，进入下一步骤。反之，将进入步骤s3进行循环迭代训练；
步骤s8：训练完成的网络模型，输入一个通过数据预处理的被试磁共振图像数据，得到被试的输出标签。
可选地，所述步骤s1中，对给定被试的磁共振图像数据集进行样本划分的具体方式为：根据被试的样本标签进行样本的随机划分，将被试对应的样本数据分为训练集和测试集。
可选地，所述步骤s2中，对被试的磁共振图像数据预处理、水平翻转扩增操作和数据归一化操作具体为：使用基于体素的形态测量(vbm)对所有数据进行预处理，将图像归一化到一个标准模版空间中，然后提取数据的灰质信息，并对其进行平滑操作，最后通过水平翻转进行数据扩增。
可选地，所述步骤s3中，输入模型训练的数据进行随机打乱。
可选地，所述步骤s4中，在pytorch框架中采用sgd优化方法对训练超参数动态调整，优化训练。基网络g1和g2分别对训练磁共振图像进行特征提取，得到的特征h1和h2，通过对特征h1和h2进行通道拼接，得到一个融合的特征h3。
可选地，所述步骤s5中，对于基网络提取的特征h1和h2，以及融合的特征h3，在通过对应的决策层后，分别得到分类决策层的标签p1、p2和p3。对p1、p2和p3进行拼接，得到对应输入样本最终的分类标签b'。
可选地，所述步骤s6中，对分类得到的标签b'和训练真实标签b进行误差计算，得出损失值，利用网络的反向传播，对模型进行参数调整。
可选地，所述步骤s7中，通过输入磁共振图像测试样本计算模型分类的准确率，对训练的分类模型进行评估，保留最优模型参数。
有益效果：本发明使用了新颖的行卷积神经网络，定义了行卷积操作，不同于传统方格卷积，该卷积操作可以有效获取图结构数据上的局部信息；定义了特征融合方法和决策融合方法，构造的基于深度融合网络的分类模型，能够有效利用多层次多分类器不同尺度的特征信息，从而获得更好的阿尔兹海默病分类效果。
附图说明
图1是本发明的总体流程图；
图2是本发明的预处理流程图；
图3是本发明的总体框架图；
图4是本发明的网络结构框架图。
具体实施方式
为了将本发明实施例中的技术方案描述的完整、清楚，下面将结合本发明实施例中的附图展开进一步详细说明。应当理解，此处所描述的具体实施例仅用于解释本发明，并不限定本发明。
请参阅图1，本发明提供了一种技术方案：一种基于融合网络的阿尔兹海默病分类方法，包括以下步骤：
步骤s1：从adni(http://adni.loni.usc.edu/)公开数据集获取被试的磁共振图像数据集，该数据集为t1加权的结构mri扫描，经过特定的图像预处理步骤，包括多平面重建(mpr)、灰度偏差(gradwarp)、b1非均匀性校正和n3强度归一化。数据集包含磁共振图像和被试的标签，其中标签是根据受教育程度和其他标准修正的记忆任务的得分来提取的。这些标签分别是ad(诊断为ad的轻度痴呆患者)、mci(ad前驱期轻度认知障碍患者)和nc(老年对照组参与者)。磁共振图像数据集记为c＝{(ai,bi)|i∈[1,r]}，式中，a表示磁共振图像，b表示被试的对应标签，r表示样本数量，i表示样本下标。将磁共振图像全集c拷贝并划分成磁共振图像训练集e＝{(ai,bi)|i∈[1,m]}和磁共振图像测试集e＝{(ai,bi)|i∈[1,n]}，其中n m＝r，且m和n分别表示划分后的样本数量；
步骤s2：使用基于体素的形态测量(vbm)，对磁共振图像集进行数据预处理和数据归一化操作，数据处理流程如图2所示。利用matlab平台的医学影像预处理工具spm12中的cat12包，将图像归一化到一个标准模版空间中，然后将数据分割为灰质、白质和脑脊液，提取包含丰富信息的灰质图像作为后续神经网络训练的输入图像。灰质图像的输出维度为121*145*121。在得到灰质图像后，使用spm12中的smooth模块对灰质图像进行平滑(fmwh＝4mm)以得到更高的信噪比和减小个体间的解剖差异。得到处理后的磁共振图像训练集e'＝{(ai,bi)|i∈[1,m']}，经过水平翻转的数据扩增后，其数量相应的增多至m'；得到磁共振图像测试集f'＝{(ai,bi)|i∈[1,n]}；
以vgg11和resnet18的网络结构为基础，构建基于特征融合和决策融合的深度融合网络架构，网络结构如图3所示。编写深度卷积神经网络代码，并利用深度学习框架工具pytorch实现算法。
步骤s3：将磁共振图像训练集e'作为输入，载入到计算机内存中，并使用python编程语言中的shuffle语句，将磁共振图像训练集e'＝{(ai,bi)|i∈[1,m']}随机打乱；
步骤s4：分别以vgg11和resnet18为基础，构建基网络g1和g2，网络模型具体结构如图4，网络具体参数如表1。
表1本发明的实例网络参数表
在pytorch框架中使用sgd优化方法，对基网络g1和g2的网络参数进行初始化，得到基网络g1和g2的初始参数，并设定相关训练超参数(迭代次数、批大小、学习率、衰减动量等)，用于网络模型的训练优化。通过基网络g1和g2分别对训练磁共振图像进行特征提取，得到特征h1和h2，并对提取的特征进行特征通道拼接得到融合特征h3；
步骤s5：构建融合特征分类决策层p，将基网络g1和g2提取到的特征以及融合特征h3分别输入到各网络的决策层，分别得到分类决策层的标签p1、p2和p3。对p1、p2和p3进行决策通道拼接融合，得到对应输入样本最终的分类标签b'；
步骤s6：通过损失函数计算，计算训练数据通过网络分类得到的标签b'和训练真实标签b之间的误差，得出损失值。利用得出的损失值进行反向传播，对基网络g1和g2，以及融合网络分别进行网络参数调整；
步骤s7：利用测试集f'对融合网络模型进行评估，保留网络模型g。同时在参数更新结束阶段进行判断，判断训练迭代次数是否已达到最大迭代次数，若已经达到最大迭代次数，则训练阶段结束，进入下一步骤。反之，将进入步骤s3进行循环迭代训练；
融合网络结合vgg11和resnet18结构的优点，使得网络在上采样的同时结合图像不同深度的特征，从而更加全面地利用了特征图信息，也防止了深层网络中的网络退化问题，并且融合的网络避免了单个基网络的误差领域，提升网络的分类精度和泛化能力。
步骤s8：训练完成的网络模型，输入一个通过数据预处理的被试磁共振图像数据，得到被试的输出标签，是ad、mci或是nc。
综上所述，本发明采用一种基于融合网络的阿尔兹海默病分类方法，融合网络结构在多个深度卷积神经网络的基础上加入特征融合和决策融合操作，使得网络更加全面地捕捉到图像的特征，并且避免了单个基网络的误差领域，提升网络的分类精度和泛化能力，可以更好地对阿尔兹海默病进行快速、准确的分类。
对于本领域的技术人员，本发明不限于上述实施例的细节，在不背离本发明的精神和范围的情况下，能够以其它具体形式实现本发明。此外，本领域的技术人员可以对本发明进行各种改动和变形而不脱离本发明的精神和范围，这些改进和变形也应视为本发明的保护范围。
1.一种基于融合网络的阿尔兹海默病分类方法，其特征在于，所述方法包括以下步骤：
步骤s1：给定被试的磁共振图像数据集，数据集包含磁共振图像和被试的标签，磁共振图像数据集记为c＝{(ai,bi)|i∈[1,r]}，式中，a表示磁共振图像，b表示被试的对应标签，r表示样本数量，i表示样本下标，将磁共振图像全集c拷贝并划分成磁共振图像训练集e＝{(ai,bi)|i∈[1,m]}和磁共振图像测试集e＝{(ai,bi)|i∈[1,n]}，其中n m＝r，且m和n分别表示划分后的样本数量；
步骤s2：对磁共振图像训练集依次进行数据预处理、水平翻转扩增操作和数据归一化操作，得到处理后的磁共振图像训练集e'＝{(ai,bi)|i∈[1,m']}，经过扩增后，其数量相应的增多至m'；对磁共振图像测试集f依次进行数据预处理、数据归一化操作得到归一化的磁共振图像测试集f'＝{(ai,bi)|i∈[1,n]}；
步骤s3：将磁共振图像训练集e'作为输入，载入到计算机内存中，并将磁共振图像训练集e'＝{(ai,bi)|i∈[1,m']}随机打乱，为接下来的训练阶段做准备；
步骤s4：分别构建基网络g1和g2，在pytorch框架中使用sgd优化方法，对基网络g1和g2的网络参数进行初始化，得到基网络g1和g2的初始参数，并设定相关训练超参数，用于网络模型的训练优化，通过基网络g1和g2分别对训练磁共振图像进行特征提取，得到特征h1和h2，并对分别提取的特征进行特征融合得到融合特征h3；
步骤s5：构建融合特征分类决策层p，将基网络g1和g2提取到的特征h1和h2以及融合特征h3分别输入到各网络的决策层，分别得到分类决策层的标签p1、p2和p3，对p1、p2和p3进行融合，得到对应输入样本最终的分类标签b'；
步骤s6：通过损失函数计算，计算训练数据通过网络分类得到的标签b'和训练真实标签b之间的误差，得出损失值，利用得出的损失值进行反向传播，对基网络g1和g2，以及融合网络分别进行网络参数调整；
步骤s7：利用测试集f'对融合网络模型进行评估，保留网络模型g，同时在参数更新结束阶段进行判断，判断训练迭代次数是否已达到最大迭代次数，若已经达到最大迭代次数，则训练阶段结束，进入下一步骤，反之，将进入步骤s3进行循环迭代训练；
步骤s8：训练完成的网络模型，输入一个通过数据预处理的被试磁共振图像数据，得到被试的输出标签。
2.根据权利要求书1所述的一种基于融合网络的阿尔兹海默病分类方法，其特征在于，所述步骤s1中，对给定被试的磁共振图像数据集进行样本划分的具体方式为：根据被试的样本标签进行样本的随机划分，将被试对应的样本数据分为训练集和测试集。
3.根据权利要求书1所述的一种基于融合网络的阿尔兹海默病分类方法，其特征在于，所述步骤s2中，对被试的磁共振图像数据预处理、水平翻转扩增操作和数据归一化操作具体为：使用基于体素的形态测量(vbm)对所有数据进行预处理，将图像归一化到一个标准模版空间中，然后提取数据的灰质信息，并对其进行平滑操作，最后通过水平翻转进行数据扩增。
4.根据权利要求书1所述的一种基于融合网络的阿尔兹海默病分类方法，其特征在于，所述步骤s3中，输入模型训练的数据进行随机打乱。
5.根据权利要求书1所述的一种基于融合网络的阿尔兹海默病分类方法，其特征在于，所述步骤s4中，在pytorch框架中采用sgd优化方法对训练超参数动态调整，优化训练，基网络g1和g2分别对训练磁共振图像进行特征提取，得到的特征h1和h2，通过对特征h1和h2进行通道拼接，得到一个融合的特征h3。
6.根据权利要求书1所述的一种基于融合网络的阿尔兹海默病分类方法，其特征在于，所述步骤s5中，对于基网络提取的特征h1和h2，以及融合的特征h3，在通过对应的决策层后，分别得到分类决策层的标签p1、p2和p3，对p1、p2和p3进行拼接，得到对应输入样本最终的分类标签b'。
7.根据权利要求书1所述的一种基于融合网络的阿尔兹海默病分类方法，其特征在于，所述步骤s6中，对分类得到的标签b'和训练真实标签b进行误差计算，得出损失值，利用网络的反向传播，对模型进行参数调整。
8.根据权利要求书1所述的一种基于融合网络的阿尔兹海默病分类方法，其特征在于，所述步骤s7中，通过输入磁共振图像测试样本计算模型分类的准确率，对训练的分类模型进行评估，保留最优模型参数。
技术总结