一种阿尔兹海默病MRI图像分类方法技术

本发明专利技术公开了一种阿尔兹海默病MRI图像的分类方法。发明专利技术使用度量学习的思想构造损失函数并用于训练一组卷积神经网络。一方面，对于每一个3D数据在Axial、Sagittal、coronal三个方向上得到多个2D形式切片，同时输入到一组卷积神经网络中，既在一定程度上保留了3D的数据信息，又避免了使用3D卷积神经网络难以找到合适的数据集进行预训练而过拟合严重的问题；另一方面，以度量学习的思想构造损失函数训练得到的深度学习模型提取出的特征更具有可分性。最终可以得到更好的分类结果。

【技术实现步骤摘要】
一种阿尔兹海默病MRI图像分类方法
本专利技术涉及图像分类
，具体涉及一种阿尔兹海默病MRI图像分类方法。
技术介绍
医学图像处理作为图像处理和机器学习的一个发展方向，是关系人类生活最为密切的领域之一。随着人口老龄化的发展，阿尔兹海默病作为最常见的老年痴呆症之一，给人们的生活尤其是患者及其家人的生活带来了极大的影响。阿尔兹海默病(Alzheimer’sDisease,AD)的分类作为医学图像分类领域的一个重要分支，对AD的计算机辅助诊断具有重要意义，特别是对于病情的前期诊断及病情恶化的及时控制尤为重要。本专利技术采用的影像数据为MRI(MagneticResonanceImaging)图像。MRI图像能够提供丰富的脑组织形态特征，特别是与阿尔兹海默病相关的海马体、颞叶、额叶等部位的特征，同时阿尔兹海默病患者脑部会出现大量的萎缩，在MRI图像中可以直观看出。传统的阿尔茨海默病分类方法主要分为特征提取和分类两个部分：(1)首先通过各种人工或半人工的手段从采集的原始影像数据(MRI，PET等)中提取特征例如灰质体积、皮层厚度、海马体积等；(2)将(1)中提取到的特征送入传统的分类器——支持向量机(SVM)、随机森林(RF)等中进行分类，但该分类方法存在特征提取困难、分类精度低的缺点。
技术实现思路
针对现有技术中的上述不足，本专利技术提供的一种阿尔兹海默病MRI图像分类方法解决了现有分类方法特征提取困难、分类精度低的问题。为了达到上述专利技术目的，本专利技术采用的技术方案为：一种阿尔兹海默病MRI图像分类方法，包括以下步骤：S1、采集足够数量的AD、MCI和NC三类人群的MRI图像；S2、对MRI图像进行预处理，得到Axial、Sagittal和coronal三个方向上的多个2D形式的切片，并将图像数据集分为训练集和测试集；S3、构建深度学习模型，并对深度学习模型中每个子网络进行预训练；所述深度学习模型中的每个子网络都是一个卷积神经网络，子网络个数与每个样本的切片数量相等；S4、构建基于度量学习的损失函数，将训练集中的训练数据作为深度学习模型的输入，并通过反向传播及随机梯度下降法优化网络参数，得到测试模型；S5、将测试集中的测试数据输入测试模型中，通过前向传播得到输出数据，将输出数据作为分类特征输入KNN分类器得到分类结果。进一步地：所述步骤S3中的预训练具体为：在每个子网络的最后一层增加softmax层，通过ImageNet数据集对每个子网络单独进行预训练，预训练完成后，去除softmax层，并用各个子网络的参数预训练结果初始化此深度学习模型中的相应参数。进一步地：所述步骤S4中的损失函数J具体为：J＝J1+α·J2+regularization上式中，J1为第一损失函数，J2为第二损失函数，α为大于0的常数，regularization为网络参数的正则化项。其中，第一损失函数为：上式中，N为训练样本的个数，训练集D＝{X1,X2,…,Xn,…,XN}；i,j＝1,2,…,M，M为切片数量，即为子网络个数，为深度学习模型第i个子网络的输出，其输入为为样本Xn的第i个切片，d(x,y)为向量x,y之间的相似度度量；第二损失函数为：上式中，m＝1,2,…,M，从训练集D得到两个样本对集P和N，即对于若Xi,Xj属于同一个类别，则(Xi,Xj)∈P，否则(Xi,Xj)∈N；为模型第m个子网络的输出，其输入为为样本Xi的第m个切片；C为大于0的常数。进一步地：步骤S5的具体步骤为：S51、把训练数据集D＝{X1,X2,…,Xn,…,XN}的输出Z＝{Z1,Z2,…,Zn,…,ZN}作为KNN分类器的训练数据，并把Xn的标签作为Zn的标签，为KNN分类器设置K值；S52、将测试数据Xtest输入到测试模型中输出Ztest，通过KNN分类器对Ztest进行分类，得到测试数据Xtest的类别；Ztest与训练数据Zn的距离计算公式为：上式中，为Ztest与训练数据集中第n个数据Xn的输出Zn的距离；为深度学习模型第m个子网络的输出，ωm为第m个网络的权重。本专利技术的有益效果为：本专利技术使用度量学习的思想构造损失函数并用于训练一组卷积神经网络。一方面，对于每一个3D数据在Axial、Sagittal、coronal三个方向上得到多个2D形式切片，同时输入到一组卷积神经网络中，既在一定程度上保留了3D的数据信息，又避免了使用3D卷积神经网络难以找到合适的数据集进行预训练而过拟合严重的问题；另一方面，以度量学习的思想构造损失函数训练得到的深度学习模型提取出的特征更具有可分性。最终可以得到更好的分类结果。附图说明图1为本专利技术流程图；图2为本专利技术的模型结构图。具体实施方式下面对本专利技术的具体实施方式进行描述，以便于本
的技术人员理解本专利技术，但应该清楚，本专利技术不限于具体实施方式的范围，对本
的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本专利技术的精神和范围内，这些变化是显而易见的，一切利用本专利技术构思的专利技术创造均在保护之列。如图1所示，一种阿尔兹海默病MRI图像的分类方法，包括以下步骤：S1、采集足够数量的AD(阿尔兹海默病)、MCI(轻度认知障碍)和NC(正常)三类人群的MRI(磁共振成像)图像。S2、对MRI图像进行预处理，得到Axial(轴向)、Sagittal(径向)、coronal(冠状)三个方向上的多个2D形式的切片，并将图像数据集分为训练集和测试集。对3个方向上对每个MRI图像切片分别得到尺寸为79×95×1、95×79×1、79×95×1的3组2D图像，并把每个切片图像复制2个通道变成3个通道的图像，即79×95×3、95×79×3、79×95×3。S3、构建深度学习模型，并对深度学习模型中每个子网络进行预训练；所述深度学习模型中的每个子网络都是一个卷积神经网络，子网络个数与每个样本的切片数量相等。在每个子网络的最后一层增加softmax层，通过ImageNet数据集对每个子网络单独进行预训练，预训练完成后，去除softmax层，并用各个子网络的参数预训练结果初始化此深度学习模型中的相应参数。如图2所示，整个深度学习模型包含M个子网络，M为每个受试者的MRI图像在Axial(轴向)、Sagittal(径向)、coronal(冠状)三个方向上进行切片获得的2D图像数，即M＝60。每个子网络的结构根据实际效果可以不同也可以相同，但需保证输出尺寸相同，其中输入数据为Axial方向上的切片的子网络的输入尺寸为79×95×3，输入数据为Sagittal方向上的切片的子网络的输入尺寸为79×79×3，输入数据为coronal方向上的切片的子网络的输入尺寸为95×79×3。S4、构建基于度量学习的损失函数，将训练集中的训练数据作为深度学习模型的输入，并通过反向传播及随机梯度下降法优化网络参数，得到测试模型。在训练过程中，每次迭代从训练集中随机选取一个batch的数据输入到网络中，batch的大小设为N，所以每个batch共有N×M张2D图像，即X＝{X1,X2,…,Xn,…,XN},其中每个数据Xn包括M张图像，即对于其中的每个子网络，一个batch包括N张2D图像，即第m个子网络输入的是一个batch中每个数据的第本文档来自技高网...

【技术保护点】
1.一种阿尔兹海默病MRI图像的分类方法，其特征在于，包括以下步骤：S1、采集足够数量的AD、MCI和NC三类人群的MRI图像；S2、对MRI图像进行预处理，得到Axial、Sagittal和coronal三个方向上的多个2D形式的切片，并将图像数据集分为训练集和测试集；S3、构建深度学习模型，并对深度学习模型中每个子网络进行预训练；所述深度学习模型中的每个子网络都是一个卷积神经网络，子网络个数与每个样本的切片数量相等；S4、将训练集中的训练数据作为深度学习模型的输入，构建基于度量学习的损失函数，并通过反向传播及随机梯度下降法优化网络参数，得到测试模型；S5、将测试集中的测试数据输入测试模型中，通过前向传播得到输出数据，将输出数据作为分类特征输入KNN分类器得到分类结果。

【技术特征摘要】
1.一种阿尔兹海默病MRI图像的分类方法，其特征在于，包括以下步骤：S1、采集足够数量的AD、MCI和NC三类人群的MRI图像；S2、对MRI图像进行预处理，得到Axial、Sagittal和coronal三个方向上的多个2D形式的切片，并将图像数据集分为训练集和测试集；S3、构建深度学习模型，并对深度学习模型中每个子网络进行预训练；所述深度学习模型中的每个子网络都是一个卷积神经网络，子网络个数与每个样本的切片数量相等；S4、将训练集中的训练数据作为深度学习模型的输入，构建基于度量学习的损失函数，并通过反向传播及随机梯度下降法优化网络参数，得到测试模型；S5、将测试集中的测试数据输入测试模型中，通过前向传播得到输出数据，将输出数据作为分类特征输入KNN分类器得到分类结果。2.根据权利要求1所述的阿尔兹海默病MRI图像的分类方法，其特征在于，所述步骤S3中的预训练具体为：在每个子网络的最后一层增加softmax层，通过ImageNet数据集对每个子网络单独进行预训练，预训练完成后，去除softmax层，并用各个子网络的参数预训练结果初始化此深度学习模型中的相应参数。3.根据权利要求1所述的阿尔兹海默病MRI图像的分类方法，其特征在于，所述步骤S4中的损失函数J具体为：J＝J1+α·J2+regularization上式中，J1为第一损失函数，J2为第二...

【专利技术属性】
技术研发人员：程建，周娇，苏炎洲，郭桦，林莉，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51