基于相对熵损失函数卷积神经网络的脑龄预测方法技术

本发明专利技术公开了一种基于相对熵损失函数卷积神经网络的脑龄预测方法，包括如下步骤：步骤1.利用采集被试的功能型磁共振影像数据形成原始样本集；步骤2.原始样本集中功能型磁共振影像数据预处理形成三维的T1图像数据并形成样本集；步骤3.将样本集划分为训练样本集和测试样本集；步骤4.利用训练样本集训练3DCNN形成脑龄预测模型，所述3DNN在训练过程中采用分类网络的相对熵损失函数反向更新3DCNN的网络参数；步骤5.将测试集输入脑龄预测模型得到预测的脑龄；用相对熵损失函数反向更新3DCNN的网络参数，同时采用中心化功能型磁共振影像数据中剔除周围无用信息，提高了脑龄预测模型的鲁棒性及预测精度。的鲁棒性及预测精度。的鲁棒性及预测精度。

【技术实现步骤摘要】
基于相对熵损失函数卷积神经网络的脑龄预测方法


[0001]本专利技术涉及到及脑磁共振图像处理领域，具体来说涉及一种基于相对熵损失函数卷积神经网络的脑龄预测方法。

技术介绍

[0002]大脑在发育和老化的过程中表现出结构和机能上的有规律的改变，这种改变的模式非常复杂，临床上不可能用肉眼得出结果。在实际应用中，采用脑龄作为度量该过程中的变化模式的指标。
[0003]目前关于大脑的发育、老化轨迹的研究并没有十分的深入，我们迫切地需要一个高精确度的脑龄预测模型去拟合大脑在发育、老化过程中的变化轨迹。一方面可以帮助我们理解大脑的发育、老化机制，加强我们对于人脑的理解和认识；另一方面，可以将脑龄与实际年龄的差值作为多种精神疾病的疾病标记物，对于精神疾病的早期诊断具有非常重要的临床意义。
[0004]功能型磁共振成像的应用分为三种情况：(1)扩散成像，人体内的水分子存在布朗运动形式的随机扩散。这种扩散信息与弛豫时间T1、T2是无关的，它能在分子水平上提供功能型的信息。(2)灌注成像，在显微毛细血管层次上的血液动力学成像，传统上是用同位素成像的方法来解决的。在磁共振成像中的平面回波成像方法不仅能同样提供有关的区域脑血流及脑血流量的信息，而且比传统方法具有更高的空间分辨率。(3)任务急活的图像，人体在做某项活动时，大脑皮层特殊的区域中会有相应的反映。用fMRI测定大脑血液的氧合水平就能直接进行脑功能的研究。功能磁共振成像的新技术，将正电子发射断层扫描技术和磁共振成像两项技术优势结合起来，通过检验血流进入脑细胞的磁场变化而实现脑功能成像，它给出更精确的结构与功能关系。
[0005]在现有的基于功能型的脑龄预测模型中，大多采用传统的机器学习方法，这些方法需要经过特征提取、特征选择、训练模型等过程，大量的依赖第三方工具，缺乏端到端的便捷性，且预测精度无法满足临床应用需求。
[0006]随着深度学习技术的发展，尤其是卷积神经网络，为脑龄预测提供了新的方法。卷积神经网络是一种常见的深度学习网络架构，受生物自然视觉认知机制启发而来。卷积神经网络具备非常优异的特征提取能力，在计算机视觉领域中的各种子任务中都表现出了绝对的优势。

技术实现思路

[0007]为了解决上述问题，本专利技术提供一种基于相对熵损失函数卷积神经网络的脑龄预测方法，所述脑龄预测方法精度高。
[0008]为了实现以上目的，本专利技术采用的一种技术方案如下：
[0009]一种基于相对熵损失函数卷积神经网络的脑龄预测方法，包括以下步骤：
[0010]步骤1.采集被试的功能型磁共振影像数据形成原始样本集，所述原始样本集的样
本包括被试的功能型磁共振影像数据及对应的被试的实际年龄；
[0011]步骤2.对原始样本集中功能型磁共振影像数据预处理形成三维的T1图像数据，并将T1图像数据和实际年龄一一对应形成样本集；
[0012]步骤3.将样本集划分为训练样本集和测试样本集；
[0013]步骤4.利用训练样本集训练3DCNN形成脑龄预测模型，所述3DNN在训练过程中采用分类网络的相对熵损失函数反向更新3DCNN的网络参数；
[0014]步骤5.将测试集输入脑龄预测模型得到预测的脑龄。
[0015]进一步地，所述步骤2中包含以下步骤：
[0016]步骤21.对原始样本集中功能型磁共振影像数据利用切片时间矫正、头动矫正和空间正则化预处理将噪声信息去除并形成三维图像数据；
[0017]步骤22.对所述三维图像数据中心化处理形成三维的T1图像数据。
[0018]进一步地，所述步骤3包含以下步骤：
[0019]步骤31.将样本训练集中的T1图像数据送入3DCNN的特征提取网络提取特征图；
[0020]步骤32.将特征图送至3DCNN的分类网络得到各实际年龄标签age
c
的分类概率p
c
，利用得到预测脑龄；
[0021]步骤33.利用预测脑龄pred和实际年龄age计算相对熵损失函数：其中，n为样本训练集中T1图像数据的个数，利用相对熵损失函数反向更新特征提取网络及分类网络的网络参数，直至达到预设的迭代停止条件。
[0022]进一步地，所述步骤33中预设的迭代停止条件包含迭代次数或者相对熵损失函数最大阈值。
[0023]进一步地，所述特征提取网络包含五个特征提取模块，前四个特征提取模块均包含3*3*3的卷积层、批归一化层、Relu激活函数层、步长为2的池化层，第五个特征提取模块包含1*1*1的卷积层、批归一化层、Relu激活函数层、步长为2的池化层，所述五个特征提取模块的通道数分别为32、64、128、128及32。
[0024]进一步地，所述批归一化层中的批归一化处理具体为：
[0025]步骤a.计算经过卷积层后批量T1图像数据的均值μ
β
：
[0026]步骤b.计算经过卷积层后批量T1图像数据的方差
[0027]步骤c.规划范处理：
[0028]本专利技术的有益效果在于：利用相对熵损失函数反向更新3DCNN的网络参数，同时采用中心化功能型磁共振影像数据中剔除周围无用信息，提高了脑龄预测模型的鲁棒性及预测精度；利用平均绝对误差作为模型精确度的指标，可以发现本专利技术的脑龄预测方法相对于支持向量回归算法(SVR)及深度学习法(DNN)具有更高的精度。
附图说明
[0029]图1为本申请实施例中一种基于相对熵损失函数卷积神经网络的脑龄预测方法的流程图；
[0030]图2为本申请实施例中脑龄预测模型的网络架构图。
具体实施方式
[0031]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例，本专利技术中的第五、五个均是为了说明算法模型中的不同卷积层，没有限定意义。基于本专利技术中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0032]如图2所示，脑龄预测模型对应的3DCNN网络包含特征提取网络及分类网络，所述特征提取网络包含五个特征提取模块，前四个特征提取模块均包含3*3*3的卷积层、批归一化层、Relu激活函数层、步长为2的池化层，第五个特征提取模块包含1*1*1的卷积层、批归一化层、Relu激活函数层、步长为2的池化层，所述五个特征提取模块的通道数分别为32、64、128、128及32；分类网络包含全连接层及softmax分类器。
[0033]如图1所述，本实施例中一种基于相对熵损失函数卷积神经网络的脑龄预测方法包括以下步骤：
[0034]步骤1.采集被试的功能型磁共振影像数据形成原始样本集，所述原始样本集的样本包括被试的功能型磁共振影像数据及对应的被试的实际年龄；
[0035]步骤2.对原始样本集中功能型磁共振影像数据预处理形成三维的T1图像数据，并将T1图像数据和实际年龄本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于相对熵损失函数卷积神经网络的脑龄预测方法，其特征在于，包括以下步骤：步骤1.采集被试的功能型磁共振影像数据形成原始样本集，所述原始样本集的样本包括被试的功能型磁共振影像数据及对应的被试的实际年龄；步骤2.对原始样本集中功能型磁共振影像数据预处理形成三维的T1图像数据，并将T1图像数据和实际年龄一一对应形成样本集；步骤3.将样本集划分为训练样本集和测试样本集；步骤4.利用训练样本集训练3DCNN形成脑龄预测模型，所述3DNN在训练过程中采用分类网络的相对熵损失函数反向更新3DCNN的网络参数；步骤5.将测试集输入脑龄预测模型得到预测的脑龄。2.根据权利要求1所述的基于相对熵损失函数卷积神经网络的脑龄预测方法，其特征在于，所述步骤2中包含以下步骤：步骤21.对原始样本集中功能型磁共振影像数据利用切片时间矫正、头动矫正和空间正则化预处理将噪声信息去除并形成三维图像数据；步骤22.对所述三维图像数据中心化处理形成三维的T1图像数据。3.根据权利要求1所述的基于相对熵损失函数卷积神经网络的脑龄预测方法，其特征在于，所述步骤3包含以下步骤：步骤31.将样本训练集中的T1图像数据送入3DCNN的特征提取网络提取特征图；步骤32.将特征图送至3DCNN的分类网络得到各实际年龄标签age
c
...

【专利技术属性】
技术研发人员：丁卫平，沈鑫杰，黄嘉爽，李铭，秦廷桢，孙颖，潘柏儒，冯志豪，鞠恒荣，程纯，曹金鑫，
申请(专利权)人：南通大学，
类型：发明
国别省市：