一种卷积神经网络的脑龄预测方法及系统

1.本发明涉及脑龄预测模型技术领域，具体地说，涉及一种卷积神经网络的脑龄预测方法及系统。


背景技术：

2.卷积神经网络是一类包含卷积计算且具有深度结构的前馈神经网络，是深度学习的代表算法之一。卷积神经网络具有表征学习能力，能够按其阶层结构对输入信息进行平移不变分类，因此也被称为“平移不变人工神经网络”。一般地，cnn的基本结构包括两层，其一为特征提取层，每个神经元的输入与前一层的局部接受域相连，并提取该局部的特征。一旦该局部特征被提取后，它与其它特征间的位置关系也随之确定下来;其二是特征映射层，网络的每个计算层由多个特征映射组成，每个特征映射是一个平面，平面上所有神经元的权值相等。特征映射结构采用影响函数核小的sigmoid函数作为卷积网络的激活函数，使得特征映射具有位移不变性。此外，由于一个映射面上的神经元共享权值，因而减少了网络自由参数的个数。卷积神经网络中的每一个卷积层都紧跟着一个用来求局部平均与二次提取的计算层，这种特有的两次特征提取结构减小了特征分辨率。卷积神经网络(cnn)可以通过mri脑部结构扫描准确预测健康个体的实际年龄。这些模型可能会应用于常规临床检查，以检测健康老龄化的偏差，包括早期神经退行性变。这可能对患者护理、药物开发和优化mri数据收集具有重要意义。
3.医院常规核磁共振数据体素不符合mni模板的基本要求，无法清晰的分割脑白质和脑灰质，从而不能获取白质图和灰质图。我们无法利用3d卷积神经网络脑年龄预测模型对脑瘫患儿进行脑年龄预测。2d卷积神经网络脑年龄预测模型的建立为解决这种问题提供了可能性。但由于模型的训练仍是基于白质和灰质图像的输入，故不能满足需求。


技术实现要素：

4.（一）解决的技术问题针对现有技术的不足，本发明的目的在于提供一种卷积神经网络的脑龄预测方法及系统，采用训练模型和预测模型相结合的方式进行脑龄预测，能够让数据得到充分训练，从而得到更准确的数据特征，构建预测模型能够规避白质和灰质图像的输入，使得预测脑龄由特征数据的输入和输出得到结果，满足了脑龄预测的需求。
5.（二）技术方案本为实现上述目的，本发明提供如下技术方案。
6.本发明第一方面提供了一种卷积神经网络的脑龄预测方法，包括：s101，采集健康人群的大脑核磁共振数据并进行预处理，获取采集健康人群的姓名、年龄、性别数据，各测试者的年龄范围为5岁到90岁之间，利用数据集中的弥散张量成像数据及对应测试者的实际年龄形成初始样本集，对所述训练样本集中的各所述历史结构磁共振影像数据预处理，得到对应的处理图像数据集；
s102，在采集的核磁共振数据中随机选取一定量数据，随机选取的数据设为训练集，图像按对应采集人群的实际年龄划分为年龄分布均匀的训练集，剩下的数据设为测试集；s103，将训练集放入训练模型中，训练过程中批归一化层数据依赖与计算访存比在批归一化层前向计算过程中，首先求输入数据的均值与方差，接着使用均值、方差对每个输入数据进行归一化及缩放操作，需要保存每层前向计算时的输出结果，用于反向计算过程中参数误差、输入误差的计算；s104，对训练模型中的数据进行分层采样，分层采样针对训练集中不同年龄段的样本数进行分层采样，对每个年龄段进行标记，得到标记数据值，将标记数据值较少的放入训练模型中，进行多次训练；s105，通过均值绝对误差、相关系数、测定系数、大脑年龄与生理年龄差异评估预测模型的性能。
7.作为优选方案，所述训练模型的过程包括：s106，采用三个卷积核组成2d卷积层连进行特征提取，每个神经元对图像数据局部进行感知，然后在更高层将局部的信息综合起来得到全局信息，给定大脑核磁共振图片数据，用多个卷积核去扫这张图来学习不同的特征来提取原图特征；s107，采用9个跳跃块组成跳跃层捕获更深层次的特征，跳跃层降低了大脑核磁共振图片数据的尺寸，增大了感受野，在某个范围内选择一个数字，让这个数字代表这个范围内的所有的像素得值，增加了鲁棒性，图片中某个目标单纯的位置的移动，不应该影响识别结果，而跳跃层捕捉的恰好是目标的特征，增加了平移不变性，在保留特征信息的前提下，降低了图片的尺寸，从而提升训练速度；s108，通过3组卷积层再次进行特征学习；s109，继续使用同样的跳跃层捕获更深层次的特征；s110，最后通过全连接层回归脑年龄，经多个卷积层和池化层后，连接着1个或1个以上的全连接层，全连接层中的每个神经元与其前一层的所有神经元进行全连接．全连接层可以整合卷积层或者池化层中具有类别区分性的局部信息，为了提升卷积神经网络的性能，全连接层每个神经元的激励函数一般采用relu函数，最后一层全连接层的输出值被传递给一个输出。
8.作为优选方案，所述s101中核磁共振数据预处理是利用fsl脑成像数据分析工具库的线性配准fsllirt命令，将所有核磁共振数据配准到mni152标准空间，并重采样为128x128x18， 体素大小为1.42x1.70x10.11，并对所有数据进行z-score标准化。
9.作为优选方案，所述预测模型的构建包括：s111，提取训练模型中的特征数据，所述特征数据为s110中最后通过全连接层的数据；s112，根据特征数据构建预测模型，在模型构建的过程中卷积层使用步长为32，尺寸为5
×
5的卷积核，池化层使用步长为2，尺寸为2
×
2的窗口，卷积层使用步长为645
×
5的卷积核，池化层使用步长为2，尺寸为2
×
2的窗口，特征处理层使用1024个节点，输出层使用10个节点；s113，将测试集中的数据输入预测模型中；
s114，输出预测的脑龄数据。
10.作为优选方案，所述s101中，大脑核磁共振数据预处理的步骤为：第一步：获取脑部核磁共振图像数据，对脑部核磁共振图像数据进行空间域滤波技术，得到去噪的脑部核磁共振图像；第二步：对脑部核磁共振图像进行图像锐化处理，得到锐化后的脑部核磁共振图像；第三步：将去噪的脑部核磁共振图像和锐化后的脑部核磁共振图像融合在一起，获得到融合后的脑部核磁共振图像，根据灰度值对融合后的脑部核磁共振图像进行分割，得到分割后的脑部核磁共振图像。
11.作为优选方案，所述脑龄数据包括训练集的生理年龄及其预测年龄、测试集的生理年龄及其预测年龄。
12.作为优选方案，在所述s110中全连接层为两层全连接层，并且在全连接层中加入dropout层。
13.作为优选方案，评估预测模型的性能中均值绝对误差，采用公式如下：其中，mae为均值绝对误差，n为测试数据集的样本数，为预测年龄，为生理年龄。
14.本发明第二方面提供了一种卷积神经网络的脑龄预测系统，包括：采集模块，采集健康人群的核磁共振数据，采集模块为存储计算机可执行指令的存储器，存储器采集历史样本集，历史样本集包括多对历史结构磁共振影像数据及对应测试者的生理年龄；选取模块，在采集的核磁共振数据中随机选取一定量数据，选取模块为pythonnumpy随机抽模块；训练模块，对训练集中的数据进行充分训练，训练模块均包含尺寸为5*5*5、步长为32的卷积层、批归一化层、relu激活函数层、步长为2的池化层；构建模块，根据训练模型中提取的特征数据构建预测模型；预测模块，用于对预测集中的数据进行脑龄预测。
15.作为优选方案，通过提取模块提取训练模型中的特征数据，并由输入模块输入到预测模块中，输出模块输出预测模块中的脑龄数据，再由评估模块通过均值绝对误差、相关系数、测定系数、大脑年龄与生理年龄差异进行评估。
16.有益效果与现有技术相比，本发明提供了一种卷积神经网络的脑龄预测方法及系统，具备以下有益效果：本发明采用训练模型和预测模型相结合的方式进行脑龄预测，其中训练模型采用三个卷积核组成2d卷积层连进行特征提取，采用9个跳跃块组成跳跃层捕获更深层次的特征，通过3组卷积层再次进行特征学习，继续使用同样的跳跃层捕获更深层次的特征，能够让数据得到充分训练，从而得到更准确的数据特征，采用构建预测模型能够规避白质和灰质图像的输入，使得预测脑龄由特征数据的输入和输出得到结果，满足了脑龄预测
的需求，并且采用了均值绝对误差、相关系数、测定系数、大脑年龄与生理年龄差异评估预测模型的性能，使得预测模型的精准度得到提高。
附图说明
17.图1为本发明的卷积神经网络的脑龄预测方法流程示意图；图2为本发明的训练模型流程示意图；图3为本发明的预测模型流程示意图；图4为本发明的卷积神经网络的脑龄预测系统示意图；图5为本发明的预测模型系统示意图；图6为本发明的卷积神经网络的结构示意图；图7为本发明的实验数据来源表1。
具体实施方式
18.为了更好地了解本发明的目的、结构及功能，下面结合附图，对本发明一种卷积神经网络的脑龄预测方法及系统做进一步详细的描述，在本说明书的描述中，参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
19.请参阅图1-5，一种卷积神经网络的脑龄预测方法，包括：s101，首先，采集健康人群的核磁共振数据并进行预处理，获取采集健康人群的姓名、年龄、性别数据，各测试者的年龄范围为5岁到90岁之间，利用数据集中的弥散张量成像数据及对应测试者的实际年龄形成初始样本集，对所述训练样本集中的各所述历史结构磁共振影像数据预处理，得到对应的处理图像数据集；s102，在采集的核磁共振数据中随机选取一定量数据，随机选取的数据设为训练集，图像按对应采集人群的实际年龄划分为年龄分布均匀的训练集，剩下的数据设为测试集，随机选取的数据即训练集用于训练模型；s103，将训练集放入训练模型中，训练过程中批归一化层数据依赖与计算访存比在批归一化层前向计算过程中，首先求输入数据的均值与方差，接着使用均值、方差对每个输入数据进行归一化及缩放操作，需要保存每层前向计算时的输出结果，用于反向计算过程中参数误差、输入误差的计算；s104，对训练模型中的数据进行分层采样，分层采样针对训练集中不同年龄段的样本数进行分层采样，对每个年龄段进行标记，得到标记数据值，将标记数据值较少的放入训练模型中，进行多次训练。脑龄数据包括训练集的生理年龄及其预测年龄、测试集的生理年龄及其预测年龄；s105，通过均值绝对误差、相关系数、测定系数、大脑年龄与生理年龄差异评估预测模型的性能。
20.训练模型的过程包括：s106，采用三个卷积核组成2d卷积层连进行特征提取，每个神经元对图像数据局
部进行感知，然后在更高层将局部的信息综合起来得到全局信息，给定大脑核磁共振图片数据，用多个卷积核去扫这张图来学习不同的特征来提取原图特征；s107，采用9个跳跃块组成跳跃层捕获更深层次的特征，跳跃层降低了大脑核磁共振图片数据的尺寸，增大了感受野，在某个范围内选择一个数字，让这个数字代表这个范围内的所有的像素得值，增加了鲁棒性，图片中某个目标单纯的位置的移动，不应该影响识别结果，而跳跃层捕捉的恰好是目标的特征，增加了平移不变性，在保留特征信息的前提下，降低了图片的尺寸，从而提升训练速度；s108，通过3组卷积层再次进行特征学习；s109，继续使用同样的跳跃层捕获更深层次的特征；s110，最后通过全连接层回归脑年龄，经多个卷积层和池化层后，连接着1个或1个以上的全连接层，全连接层中的每个神经元与其前一层的所有神经元进行全连接．全连接层可以整合卷积层或者池化层中具有类别区分性的局部信息，为了提升卷积神经网络的性能，全连接层每个神经元的激励函数一般采用relu函数，最后一层全连接层的输出值被传递给一个输出。
21.s101中核磁共振数据预处理是利用fsl脑成像数据分析工具库的线性配准fsllirt命令，将所有核磁共振数据配准到mni152标准空间，并重采样为128x128x18， 体素大小为1.42x1.70x10.11，并对所有数据进行z-score标准化。
22.预测模型的构建包括：s111，提取训练模型中的特征数据，所述特征数据为s110中最后通过全连接层的数据；s112，根据特征数据构建预测模型，在模型构建的过程中卷积层使用步长为32，尺寸为5
×
5的卷积核，池化层使用步长为2，尺寸为2
×
2的窗口，卷积层使用步长为645
×
5的卷积核，池化层使用步长为2，尺寸为2
×
2的窗口，特征处理层使用1024个节点，输出层使用10个节点；s113，将测试集中的数据输入预测模型中；s114，输出预测的脑龄数据。
23.在所述s101中，大脑核磁共振数据预处理的步骤为：第一步：获取脑部核磁共振图像数据，对脑部核磁共振图像数据进行空间域滤波技术，得到去噪的脑部核磁共振图像；第二步：对脑部核磁共振图像进行图像锐化处理，得到锐化后的脑部核磁共振图像；第三步：将去噪的脑部核磁共振图像和锐化后的脑部核磁共振图像融合在一起，获得到融合后的脑部核磁共振图像，根据灰度值对融合后的脑部核磁共振图像进行分割，得到分割后的脑部核磁共振图像在所述s110中全连接层为两层全连接层，并且在全连接层中加入dropout层。
24.评估预测模型的性能中均值绝对误差，采用公式如下：
其中，mae为均值绝对误差，n为测试数据集的样本数，为预测年龄，为生理年龄。
25.一种卷积神经网络的脑龄预测系统，其特征在于，包括：采集模块，采集健康人群的核磁共振数据，采集模块为存储计算机可执行指令的存储器，存储器采集历史样本集，历史样本集包括多对历史结构磁共振影像数据及对应测试者的生理年龄；选取模块，在采集的核磁共振数据中随机选取一定量数据，选取模块为pythonnumpy随机抽模块；训练模块，对训练集中的数据进行充分训练，训练模块均包含尺寸为5*5*5、步长为32的卷积层、批归一化层、relu激活函数层、步长为2的池化层；构建模块，根据训练模型中提取的特征数据构建预测模型；预测模块，用于对预测集中的数据进行脑龄预测。通过提取模块提取训练模型中的特征数据，并由输入模块输入到预测模块中，输出模块输出预测模块中的脑龄数据，再由评估模块通过均值绝对误差、相关系数、测定系数、大脑年龄与生理年龄差异进行评估。
26.实施例，采集健康人群核磁共振数据共3474例，随机选取2442例核磁共振数据为训练集，用于训练模型中，剩下的1032例核磁共振数据为测试集，由于训练集的健康人群年龄不均匀，为了使样本数量较少的年龄段获得充分训练，对训练集数据采用分层采样，对每个年龄段进行标记，得到标记数据值，将标记数据值较少的放入训练模型中，进行多次训练，再将训练模型中的特征数据输入到预测模型中，构建得到训练集的生理年龄及其预测的脑年龄，其中相关系数r=0.99，均值绝对误差mae=1.85，p＜0.0001，得到测试集的生理年龄及其预测的脑年龄，其中相关系统r=0.95，均值绝对误差mae=3.98，p＜0.05。
27.为了更好的理解本发明，作为示例而言，在一个具体的预测方法中，可以选取年龄范围在5岁-86岁的采集健康人群，以构上述的建训练集和测试集。本技术实验数据来源如表1，但不对本技术保护范围做任何限定。
28.可以理解，本发明是通过一些实施例进行描述的，本领域技术人员知悉的，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。另外，在本发明的教导下，可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此，本发明不受此处所公开的具体实施例的限制，所有落入本技术的权利要求范围内的实施例都属于本发明所保护的范围内。

技术特征：
1.一种卷积神经网络的脑龄预测方法，其特征在于，包括：s101，采集健康人群的大脑核磁共振数据并进行预处理，获取采集健康人群的姓名、年龄、性别数据，各测试者的年龄范围为5岁到90岁之间，利用数据集中的弥散张量成像数据及对应测试者的实际年龄形成初始样本集，对所述训练样本集中的各所述历史结构磁共振影像数据预处理，得到对应的处理图像数据集；s102，在采集的核磁共振数据中随机选取一定量数据，随机选取的数据设为训练集，图像按对应采集人群的实际年龄划分为年龄分布均匀的训练集，剩下的数据设为测试集；s103，将训练集放入训练模型中，训练过程中批归一化层数据依赖与计算访存比在批归一化层前向计算过程中，首先求输入数据的均值与方差，接着使用均值、方差对每个输入数据进行归一化及缩放操作，需要保存每层前向计算时的输出结果，用于反向计算过程中参数误差、输入误差的计算；s104，对训练模型中的数据进行分层采样，分层采样针对训练集中不同年龄段的样本数进行分层采样，对每个年龄段进行标记，得到标记数据值，将标记数据值较少的放入训练模型中，进行多次训练；s105，通过均值绝对误差、相关系数、测定系数、大脑年龄与生理年龄差异评估预测模型的性能。2.根据权利要求1所述的一种卷积神经网络的脑龄预测方法，其特征在于，所述训练模型的过程包括：s106，采用三个卷积核组成2d卷积层连进行特征提取，每个神经元对图像数据局部进行感知，然后在更高层将局部的信息综合起来得到全局信息，给定大脑核磁共振图片数据，用多个卷积核去扫这张图来学习不同的特征来提取原图特征；s107，采用9个跳跃块组成跳跃层捕获更深层次的特征，跳跃层降低了大脑核磁共振图片数据的尺寸，增大了感受野，在某个范围内选择一个数字，让这个数字代表这个范围内的所有的像素得值，增加了鲁棒性，图片中某个目标单纯的位置的移动，不应该影响识别结果，而跳跃层捕捉的恰好是目标的特征，增加了平移不变性，在保留特征信息的前提下，降低了图片的尺寸，从而提升训练速度；s108，通过3组卷积层再次进行特征学习；s109，继续使用同样的跳跃层捕获更深层次的特征；s110，最后通过全连接层回归脑年龄，经多个卷积层和池化层后，连接着1个或1个以上的全连接层，全连接层中的每个神经元与其前一层的所有神经元进行全连接．全连接层可以整合卷积层或者池化层中具有类别区分性的局部信息，为了提升卷积神经网络的性能，全连接层每个神经元的激励函数一般采用relu函数，最后一层全连接层的输出值被传递给一个输出。3.根据权利要求1所述的一种卷积神经网络的脑龄预测方法，其特征在于：所述s101中核磁共振数据预处理是利用fsl脑成像数据分析工具库的线性配准fsllirt命令，将所有核磁共振数据配准到mni152标准空间，并重采样为128x128x18， 体素大小为1.42x1.70x10.11，并对所有数据进行z-score标准化。4.根据权利要求2所述的一种卷积神经网络的脑龄预测方法，其特征在于，所述预测模型的构建包括：
s111，提取训练模型中的特征数据，所述特征数据为s110中最后通过全连接层的数据；s112，根据特征数据构建预测模型，在模型构建的过程中卷积层使用步长为32，尺寸为5
×
5的卷积核，池化层使用步长为2，尺寸为2
×
2的窗口，卷积层使用步长为645
×
5的卷积核，池化层使用步长为2，尺寸为2
×
2的窗口，特征处理层使用1024个节点，输出层使用10个节点；s113，将测试集中的数据输入预测模型中；s114，输出预测的脑龄数据。5.根据权利要求1所述的一种卷积神经网络的脑龄预测方法，其特征在于：所述s101中，大脑核磁共振数据预处理的步骤为：第一步：获取脑部核磁共振图像数据，对脑部核磁共振图像数据进行空间域滤波技术，得到去噪的脑部核磁共振图像；第二步：对脑部核磁共振图像进行图像锐化处理，得到锐化后的脑部核磁共振图像；第三步：将去噪的脑部核磁共振图像和锐化后的脑部核磁共振图像融合在一起，获得到融合后的脑部核磁共振图像，根据灰度值对融合后的脑部核磁共振图像进行分割，得到分割后的脑部核磁共振图像。6.根据权利要求4所述的一种卷积神经网络的脑龄预测方法，其特征在于：所述脑龄数据包括训练集的生理年龄及其预测年龄、测试集的生理年龄及其预测年龄。7.根据权利要求2所述的一种卷积神经网络的脑龄预测方法，其特征在于：在所述s110中全连接层为两层全连接层，并且在全连接层中加入dropout层。8.根据权利要求1所述的一种卷积神经网络的脑龄预测方法，其特征在于：评估预测模型的性能中均值绝对误差，采用公式如下：其中，mae为均值绝对误差，n为测试数据集的样本数，为预测年龄，为生理年龄。9.一种卷积神经网络的脑龄预测系统，其特征在于，包括：采集模块，采集健康人群的核磁共振数据，采集模块为存储计算机可执行指令的存储器，存储器采集历史样本集，历史样本集包括多对历史结构磁共振影像数据及对应测试者的生理年龄；选取模块，在采集的核磁共振数据中随机选取一定量数据，选取模块为pythonnumpy随机抽模块；训练模块，对训练集中的数据进行充分训练，训练模块均包含尺寸为5*5*5、步长为32的卷积层、批归一化层、relu激活函数层、步长为2的池化层；构建模块，根据训练模型中提取的特征数据构建预测模型；预测模块，用于对预测集中的数据进行脑龄预测。10.根据权利要求9所述的一种卷积神经网络的脑龄预测系统，其特征在于：通过提取模块提取训练模型中的特征数据，并由输入模块输入到预测模块中，输出模块输出预测模块中的脑龄数据，再由评估模块通过均值绝对误差、相关系数、测定系数、大脑年龄与生理年龄差异进行评估。

技术总结
本发明公开了一种卷积神经网络的脑龄预测方法及系统，包括：S101，首先，采集健康人群的大脑核磁共振数据并进行预处理；S102，在采集的核磁共振数据中随机选取一定量数据，随机选取的数据设为训练集，剩下的数据设为测试集；S103，将训练集放入训练模型中；S104，对训练模型中的数据进行分层采样，针对不同年龄进行分层采样；S105，通过均值绝对误差、相关系数、测定系数、大脑年龄与生理年龄差异评估预测模型的性能。本发明采用训练模型和预测模型相结合的方式进行脑龄预测，能够规避白质和灰质图像的输入，使得预测脑龄由特征数据的输入和输出得到结果，满足了脑龄预测的需求。满足了脑龄预测的需求。满足了脑龄预测的需求。


技术研发人员：侯昌龙 张建华 李仁峰
受保护的技术使用者：上海市东方医院（同济大学附属东方医院）
技术研发日：2023.03.08
技术公布日：2023/8/14