一种基于多步学习的三维核磁共振胰腺图像分割方法技术

本发明专利技术公开了一种基于多步学习的三维核磁共振胰腺图像分割方法，包括：(1)归一化预处理三维核磁共振图像；(2)随机将预处理后的图像分为训练集、验证集和测试集；(3)利用分割标注将训练集中包含胰腺的块和对应的标注切出备用；(4)压缩原始图像训练一个Q‑net模型计算胰腺大致位置；(5)使用步骤3中的配对数据预训练另一个P‑net模型；(6)使用预训练的Q‑net产生胰腺的3D位置预测图，选取概率高的图块，映射回原图，分块输入预训练的P‑net中合并训练预测胰腺位置；(7)使用训练好的分割模型在测试集上进行预测检测效果。利用本发明专利技术，能够精确地从三维的核磁共振图像中分割出胰腺，可为胰腺的放射治疗提供依据和指导。

【技术实现步骤摘要】
一种基于多步学习的三维核磁共振胰腺图像分割方法
本专利技术属于医学图像处理领域，尤其是涉及一种基于多步学习的三维核磁共振胰腺图像分割方法。
技术介绍
胰腺癌是长期以来严重威胁人类健康的一种出现在胰腺上的癌症，是常见恶性肿瘤中恶性程度最高，死亡率最高的癌症，其发病率和死亡率几乎相等。约90％以上的患者无法通过手术治疗根治，转移极快，五年存活率仅有3％。虽然随着科学技术的不断进步和发展，癌症的总体治疗手段和方法都有了很大突破，但胰腺癌的治疗并没有实质性的进展，患病死亡率也长期位居第一。因此在胰腺癌的诊断和治疗上，必须加强科学技术的研究，争取早日改善这一现状。一方面，在传统的胰腺癌诊断中，主要依靠放射科医生的经验和观察阅片进行，有一定的主观性，随着工作量的增加，诊断的准确率和效率也会有所下降。准确分割出胰腺的三维结构对于放射科医生的诊断有重要的促进作用。另一方面，对于无法手术根治的胰腺癌患者，放射治疗是一种重要的治疗方法，该方法利用强辐射射线杀死肿瘤细胞，但是也会对路径上的正常组织和细胞造成不可逆的破坏。所以在实际操作中，肿瘤组织的精准定位、放射治疗的射线方向、角度，放射使用剂量等放疗方案制定成为了放射治疗的重要因素，在精准医疗时代，这一点更被多次强调。为了制定出更科学高效的胰腺癌放疗方案，胰腺组织的三维分割十分重要，它能够有效帮助医生分析病情，精准定位和制定放疗手段。在现有的医学影像检测手段中，主要的方法为三维核磁共振成像，这种方法具有强大的成像能力，尤其对于软组织的成像更加清晰。比起其他同类检查手段而言，对人体的损害相对较小，已经成为了胰腺疾病的检查、诊断、治疗的重要医学影像手段。而针对三维核磁同振图像的胰腺分割，也就是实现以上目标的实际手段。在现有的胰腺分割方法中，如基于区域分割、基于边缘分割、基于图谱分割等方法，都是基于二维的序列图像分割的方法，利用三维核磁共振图的切片信息进行分割，无法有效利用到不同切面之间的相关信息。另一方面，由于胰腺的固有特性，比起腹腔来说体积过小，正负例样本数目差距过大，利用现有的直接的分割方法很难精确的分割出胰腺的位置，与此同时，核磁共振图像又有噪声偏大、对比度较低、临近组织边界不清晰等特点，使得这一任务难以通过普通的3D分割神经网络进行，多会出现边界模糊，分割不足或过分割的情况。
技术实现思路
本专利技术提供了一种基于多步学习的三维核磁共振胰腺图像分割方法，通过多步学习的方法，提高三维核磁同振图像中胰腺分割的准确率和精细程度。一种基于多步学习的三维核磁共振胰腺图像分割方法，包括以下步骤：(1)输入三维核磁共振序列图像，对其进行归一化处理，得到归一化后的序列图像集；(2)将序列图像集人工标注后按一定比例随机分成训练集、验证集和测试集；(3)将训练集中图像预处理成相同大小尺寸后，将含有胰腺的图块和对应的标注切出进行配对；(4)搭建分割模型，所述分割模型包含Q-net模型和P-net模型，所述Q-net模型用于对图像粗分割生成候选区域，所述P-net模型用于对候选区域的胰腺进行精准分割；(5)将训练集的图像压缩后输入Q-net模型进行预训练，将步骤(3)中切出的含有胰腺的图块和对应的标注图块输入P-net模型进行预训练；(6)将预训练好的Q-net模型和P-net模型进行合并训练，将训练集图像压缩后输入预训练好的Q-net模型进行粗分割，产生胰腺的3D位置预测图，选取高于阈值的图块映射回原图，取出这些位置的图块和对应标注输入预训练好的P-net模型中进行分割预测，达到预设的训练次数后结束训练；(7)将需要分割的图像归一化后输入训练好的分割模型中，输出最终的分割结果。本专利技术使用了多步学习的方法，使用压缩的原始图像训练Q-net来粗分割生成候选区域，再使用预训练的P-net来对候选区域中的胰腺进行精准分割，可以解决现有胰腺分割技术中胰腺分割边界粗糙，假阳性率高的问题，提高了胰腺分割的准确率和精细程度。步骤(1)和步骤(7)中，所述归一化处理的步骤如下：(1-1)计算全部图像像素值的平均值和标准差；(1-2)对于全部图像，图像像素值减去上步计算的平均值，再除以上步计算的标准差；(1-3)对于上步处理后的图像，每张图减去该图中像素最小值，再除以该图中像素值的极值差。步骤(3)中，所述预处理的步骤如下：对于长、宽、高任何一边小于设定像素N的图块，利用双线性插值扩大到N像素；对于长、宽、高任何一边大于N像素的图块，按照N×N×N进行裁剪，剩下部分如果大于则同样利用双线性插值扩大到N×N×N，否则舍弃。在对训练集中含有胰腺的图块和对应的标注图块切出配对时，需按照步骤(5)中所述的P-net的输入大小进行图块的切割、缩放等操作。步骤(4)中，所述Q-net模型和P-net模型的网络结构相同，为3D-Unet模型，模型均包含编码和解码两条路径，两条路径分别有4个卷积层；在编码路径，每个卷积层包含两个3×3×3的卷积，后接一个ReLu激活函数，接着连接一个2×2×2的最大值池化，在每个维度的步长都为2；在解码路径，每一个卷积层包括一个2×2×2的反卷积，在每个维度的步长为2，后接两个3×3×3的卷积层，最后连接ReLu激活函数；相同维度的编码和解码路径上的卷积层之间设有快捷连接。步骤(5)中，预训练Q-net模型的步骤如下：(5-1-1)对训练集的图像压缩至N×N×N，输入Q-net模型进行训练；(5-1-2)使用1-dice损失作为损失函数，使用随机梯度下降的优化方法，每个epoch训练结束后在验证集上计算dice值，当Q-net模型在验证集上的预测的dice值达到0.6时结束预训练。在Q-net模型预训练结束后，通过验证集的预测结果和标注进行阈值选择，选取能够使得验证集召回率大于90％的最大阈值作为挑选的阈值。阈值是用来界定哪些Q-net的输出作为阳性的预测结果。当阈值为1，输出结果为0个像素，阈值0时，输出结果为全部输入像素。在Q-net模型预训练完后，通过验证集，调整阈值使得输出的像素点能够包括至少90％的标注像素。步骤(5)中，预训练P-net模型的步骤如下：(5-2-1)使用步骤(3)中的切出图块和对应的标注作为训练数据，并且按8:2的比例随机分为训练集和验证集；(5-2-2)使用1-dice损失作为损失函数，使用随机梯度下降的优化方法，每个epoch训练结束后在验证集上计算dice值，当在P-net在验证集上的预测的dice值达到0.8时结束预训练。步骤(6)中，使用Q-net模型预训练后挑选的阈值作为标准，高于该阈值的预测区域作为Q-net模型的预测结果；并将这些区域根据Q-net的压缩比例映射回原图；这些图块根据P-net的输入大小进行裁剪和缩放操作后输入P-net进行端对端训练。每经过5个epoch训练后就检验阈值是否依然能够覆盖90％的标记对照，如果不能，则自动调整为符合要求的阈值继续训练，直到达到预设的训练次数。在模型训练完毕后，将训练集输入该模型，使用dice值作为最终的模型预测结果评价指标。与现有技术相比，本专利技术具有以下有益效果：1、本专利技术针对比较困难的胰腺分割任务，使用了多步学习的方法，使用压缩的原始图像训练Q-net来粗分割生成候选区域，再使用预训练的P-net来本文档来自技高网...

【技术保护点】
1.一种基于多步学习的三维核磁共振胰腺图像分割方法，其特征在于，包括以下步骤：(1)输入三维核磁共振序列图像，对其进行归一化处理，得到归一化后的序列图像集；(2)将序列图像集人工标注后按一定比例随机分成训练集、验证集和测试集；(3)将训练集中图像预处理成相同大小尺寸后，将含有胰腺的图块和对应的标注切出；(4)搭建分割模型，所述分割模型包含Q‑net模型和P‑net模型，所述Q‑net模型用于对图像粗分割生成候选区域，所述P‑net模型用于对候选区域的胰腺进行精准分割；(5)将训练集的图像压缩后输入Q‑net模型进行预训练，将步骤(3)中切出的含有胰腺的图块和对应的标注输入P‑net模型进行预训练；(6)将预训练好的Q‑net模型和P‑net模型进行合并训练，将训练集图像压缩后输入预训练好的Q‑net模型进行粗分割，产生胰腺的3D位置预测图，选取高于阈值的图块映射回原图，取出这些位置的图块和对应标注输入预训练好的P‑net模型中进行分割预测，达到预设的训练次数后结束训练；(7)将需要分割的图像归一化后输入训练好的分割模型中，输出最终的分割结果。

【技术特征摘要】
1.一种基于多步学习的三维核磁共振胰腺图像分割方法，其特征在于，包括以下步骤：(1)输入三维核磁共振序列图像，对其进行归一化处理，得到归一化后的序列图像集；(2)将序列图像集人工标注后按一定比例随机分成训练集、验证集和测试集；(3)将训练集中图像预处理成相同大小尺寸后，将含有胰腺的图块和对应的标注切出；(4)搭建分割模型，所述分割模型包含Q-net模型和P-net模型，所述Q-net模型用于对图像粗分割生成候选区域，所述P-net模型用于对候选区域的胰腺进行精准分割；(5)将训练集的图像压缩后输入Q-net模型进行预训练，将步骤(3)中切出的含有胰腺的图块和对应的标注输入P-net模型进行预训练；(6)将预训练好的Q-net模型和P-net模型进行合并训练，将训练集图像压缩后输入预训练好的Q-net模型进行粗分割，产生胰腺的3D位置预测图，选取高于阈值的图块映射回原图，取出这些位置的图块和对应标注输入预训练好的P-net模型中进行分割预测，达到预设的训练次数后结束训练；(7)将需要分割的图像归一化后输入训练好的分割模型中，输出最终的分割结果。2.根据权利要求1所述的基于多步学习的三维核磁共振胰腺图像分割方法，其特征在于，归一化处理的步骤如下：(1-1)计算全部图像像素值的平均值和标准差；(1-2)对于全部图像，图像像素值减去上步计算的平均值，再除以上步计算的标准差；(1-3)对于上步处理后的图像，每张图减去该图中像素最小值，再除以该图中像素值的极值差。3.根据权利要求1所述的基于多步学习的三维核磁共振胰腺图像分割方法，其特征在于，步骤(3)中，所述预处理的步骤如下：对于长、宽、高任何一边小于设定像素N的图块，利用双线性插值扩大到N像素；对于长、宽、高任何一边大于N像素的图块，按照N×N×N进行裁剪，剩下部分如果大于则同样利用双线性插值扩大到N×N×N，否则舍弃。4.根据权利要求1所述的基于多步学习的三维核磁共振胰腺图像分割方法，其特征在于，步骤(4)中，所述Q-net模型和P-net模型的网络结构相同...

【专利技术属性】
技术研发人员：吴健，余柏翰，王文哲，冯芮苇，陆逸飞，吴福理，
申请(专利权)人：浙江大学山东工业技术研究院，
类型：发明
国别省市：山东,37