基于深度多重解析网络的缺失CT投影数据估计方法技术

本发明专利技术公开了一种基于深度多重解析网络的缺失CT投影数据估计方法，属于医学图像重建技术领域。方法包括：获取一定数量的完全CT投影数据图像，对每个图像进行部分遮挡，得到缺失部分的CT投影数据原图和具有缺失区域的缺失CT投影数据图像；将缺失部分的CT投影数据原图和具有缺失区域的缺失CT投影数据图像作为训练数据，输入预设的深度多重解析网络模型进行训练，深度多重解析网络模型使用局部鉴别器和全局鉴别器两个鉴别器；利用S2中已训练的深度多重解析网络模型进行预测得到缺失部分的CT投影数据图像；根据预测的缺失部分CT投影数据图像，重建出CT投影数据图像。本发明专利技术的方法使得生成的CT投影数据图像的缺失区域边界更具连贯性。

【技术实现步骤摘要】
基于深度多重解析网络的缺失CT投影数据估计方法
本专利技术涉及医学图像重建
，具体涉及一种基于深度多重解析网络的缺失CT投影数据估计方法。
技术介绍
作为目前一种常规有效的临床医学诊断工具，X射线计算机断层摄影术为临床医生的诊断提供了丰富的人体器官组织信息。但是由相关研究表明：一次完全的CT扫描通常伴随着较高程度的电离辐射，而高剂量电离辐射可诱发人体新陈代谢异常乃至癌症、白血病等疾病。临床上减少病患辐射量的重要方法之一就是减小CT扫描范围，即将探测器的旋转角度范围限制在某个小于标准的区间内，从而在总体上大幅减少了患者所受X射线辐射量。虽然限制CT设备扫描范围能够降低患者所受X射线辐射量，但同时会造成所获CT投影数据部分缺失，即获得的是不完全投影数据，使重建CT图像质量明显下降，以至于无法满足临床诊断的需要。同样在多排CT成像中，X射线辐射量的降低会造成重建图像质量的明显下降。随着扫描范围的减小，虽然患者受到的辐射剂量大幅降低，但是重建图像出现大量星条状伪影和噪声，严重影响了对特征点的分辨。因此，如何在减小扫描范围，即投影数据不完全条件下重建出符合临床诊断要求、高质量的CT图像具有重要的科学意义和临床实用价值，也引起了国内外学者越来越多的重视。2014年美国北卡罗来纳大学教堂山分校IDEA研究团队利用随机森林、卷积神经网络等算法，结合自动上下文模型从MRI和低剂量PET图像中估计正常剂量PET图像或CT图像。Boublil等人提出了使用人工神经网络来提升常用CT图像重建算法性能的理论框架，并将其成功应用于低剂量医学图像重建。Dosovitskiy等人证明了可以通过通过解码器网络反转深度卷积网络特征来重建目标图像。Kingma等人提出变分自动编码器（VAE），其通过在潜在单元上施加先验来使编码器正规化，使得可以通过从潜在单元采样或插入潜在单元来生成图像。然而，由于其基于像素方式高斯似然的训练目标，VAE生成的图像通常是模糊的。Wright等人将图像完成作为从输入中恢复稀疏信号的任务，通过求解稀疏线性系统，可以从一些损坏的输入中恢复图像。然而，该算法要求图像高度结构化（即，假设数据点位于低维子空间中），例如：良好对齐的面部图像。Pathak等人提出用上下文编码器模型来对图像进行重建，但仍存在生成图像缺失边界的像素值不一致性的问题。
技术实现思路
本专利技术的目的在于克服现有技术中的不足，弥补生成图像缺失部分边界的像素值不一致的问题，提出了一种基于深度多重解析网络的缺失CT投影数据估计方法，提高生成的完全CT图像质量。为解决上述技术问题，本专利技术提供了一种基于深度多重解析网络的缺失CT投影数据估计方法，其特征是，包括以下步骤：S1，获取一定数量的完全CT投影数据图像，对每个图像进行部分遮挡，得到缺失部分的CT投影数据原图和具有缺失区域的缺失CT投影数据图像；S2，将缺失部分的CT投影数据原图和具有缺失区域的缺失CT投影数据图像作为训练数据，输入预设的深度多重解析网络模型，以训练该模型预测缺失部分CT投影数据图像；S3，对待测的缺失CT投影数据图像，利用S2中已训练的深度多重解析网络模型进行预测得到缺失部分的CT投影数据图像；S4，根据预测的缺失部分CT投影数据图像，重建出CT投影数据图像。进一步的，深度多重解析网络模型包括生成器、局部鉴别器和全局鉴别器。进一步的，深度多重解析网络模型的训练过程为：将缺失CT投影数据图像作为生成器输入，生成器输出缺失部分的CT投影数据图像；将生成器输出的缺失部分的CT投影数据图像与缺失部分的CT投影数据原图输入到局部鉴别器中，判断生成器输出的缺失部分的CT投影数据图像是真实或是合成；将生成器输出的缺失部分的CT投影数据图像与完全CT投影数据图像拼接成新的完全投影数据图像，和完全CT投影数据原图均作为全局鉴别器的输入，判断生成器输出的缺失部分的CT投影数据图像是真实或是合成。进一步的，生成器网络中包含五层卷积层、两层全连接层和五层反卷积层；卷积层共5层，卷积核尺寸均为4×4，5层步长分别为(2,2)、(3,2)、(3,4)、(4,4)、(2,2)，填充采用补零填充，5层卷积核数量分别为32，32，64，128，256；5层反卷积层，卷积核尺寸均为4×4，步长分别为(2,2)、(2,2)、(3,4)、(3,4)、(2,2)，填充采用补零填充，5层卷积核数量分别为256，128，64，32，32，每一层均进行批量归一化，激活函数采用修正线性单元函数，最后一层激活函数采用双曲正切函数。进一步的，局部鉴别器网络具有五层卷积层，包括卷积、批量归一化和激活函数操作，以及一个输出层；卷积层共5层，卷积核尺寸均为4×4，5层步长分别为(2,2)、(2,2)、(3,4)、(3,4)、(2,2)，填充采用补零填充，5层卷积核数量分别为32，64，128，128，256；每一层均进行批量归一化，激活函数采用带泄露线性整流函数。进一步的，全局鉴别器网络具有五层卷积层，包括卷积、批量归一化和激活函数操作，以及一个输出层；卷积层共5层，卷积核尺寸均为4×4，5层步长分别为(2,2)、(2,2)、(3,4)、(4,4)、(3,2)，填充采用补零填充，5层卷积核数量分别为32，64，128，128，256；每一层均进行批量归一化，激活函数采用带泄露线性整流函数。进一步的，采用卷积滤波反投影方法从预测的缺失部分CT投影数据图像中重建出CT投影数据图像。与现有技术相比，本专利技术所达到的有益效果是：本专利技术的方法使得生成的CT投影数据图像的缺失区域边界更具连贯性，psnr（PeakSingnaltoNoise，峰值信噪比）值和ssim（StructureSimilarityIndex，结构相似性）值均体现了本专利技术方法的优越性。附图说明图1是本专利技术方法的流程示意图；图2是深度多重解析网络模型。具体实施方式下面结合附图对本专利技术作进一步描述。以下实施例仅用于更加清楚地说明本专利技术的技术方案，而不能以此来限制本专利技术的保护范围。本专利技术的一种基于深度多重解析网络的缺失CT投影数据估计方法，参见图1所示，包括以下过程：步骤（1），将一定数量的完全CT投影数据图像作为训练数据，对每个图像进行部分遮挡，得到缺失（遮挡）部分的CT投影数据原图和具有缺失区域的缺失CT投影数据图像。本专利技术中训练数据图像数量是124张，图像尺寸是720×1024像素。步骤（2），构建深度多重解析网络模型，如图2所示，该模型包括生成器、局部鉴别器和全局鉴别器。其中生成器网络，将缺失CT投影数据图像作为生成器输入，输出缺失部分的CT投影数据图像。生成器网络中包含五层卷积层（编码器）、两层全连接层和五层反卷积层（解码器）。卷积层共5层，卷积核尺寸均为4×4，5层步长分别为(2,2)、(3,2)、(3,4)、(4,4)、(2,2)，填充采用补零填充，5层卷积核数量分别为32，32，64，128，256；5层反卷积层，卷积核尺寸均为4×4，步长分别为(2,2)、(2,2)、(3,4)、(3,4)、(2,2)，填充采用补零填充，5层卷积核数量分别为256，128，64，32，32，每一层均进行批量归一化，激活函数采用修正线性单元函数，最后一层激活函数采用双曲正切函数。步本文档来自技高网...

【技术保护点】
1.一种基于深度多重解析网络的缺失CT投影数据估计方法，其特征是，包括以下步骤：S1，获取一定数量的完全CT投影数据图像，对每个图像进行部分遮挡，得到缺失部分的CT投影数据原图和具有缺失区域的缺失CT投影数据图像；S2，将缺失部分的CT投影数据原图和具有缺失区域的缺失CT投影数据图像作为训练数据，输入预设的深度多重解析网络模型，以训练该模型预测缺失部分CT投影数据图像；S3，对待测的缺失CT投影数据图像，利用S2中已训练的深度多重解析网络模型进行预测得到缺失部分的CT投影数据图像；S4，根据预测的缺失部分CT投影数据图像，重建出CT投影数据图像。

【技术特征摘要】
1.一种基于深度多重解析网络的缺失CT投影数据估计方法，其特征是，包括以下步骤：S1，获取一定数量的完全CT投影数据图像，对每个图像进行部分遮挡，得到缺失部分的CT投影数据原图和具有缺失区域的缺失CT投影数据图像；S2，将缺失部分的CT投影数据原图和具有缺失区域的缺失CT投影数据图像作为训练数据，输入预设的深度多重解析网络模型，以训练该模型预测缺失部分CT投影数据图像；S3，对待测的缺失CT投影数据图像，利用S2中已训练的深度多重解析网络模型进行预测得到缺失部分的CT投影数据图像；S4，根据预测的缺失部分CT投影数据图像，重建出CT投影数据图像。2.根据权利要求1所述的一种基于深度多重解析网络的缺失CT投影数据估计方法，其特征是，深度多重解析网络模型包括生成器、局部鉴别器和全局鉴别器。3.根据权利要求2所述的一种基于深度多重解析网络的缺失CT投影数据估计方法，其特征是，深度多重解析网络模型的训练过程为：将缺失CT投影数据图像作为生成器输入，生成器输出缺失部分的CT投影数据图像；将生成器输出的缺失部分的CT投影数据图像与缺失部分的CT投影数据原图输入到局部鉴别器中，判断生成器输出的缺失部分的CT投影数据图像是真实或是合成；将生成器输出的缺失部分的CT投影数据图像与完全CT投影数据图像拼接成新的完全投影数据图像，和完全CT投影数据原图均作为全局鉴别器的输入，判断生成器输出的缺失部分的CT投影数据图像是真实或是合成。4.根据权利要求2所述的一种基于深度多重解析网络的缺失CT投影数据估计方法，其特征是，生成器网络中包含五层卷积层、两层全连接层和五层反卷积层；卷积层共5层，卷积核尺寸均为...

【专利技术属性】
技术研发人员：戴修斌，林语萱，刘天亮，晏善成，
申请(专利权)人：南京邮电大学，
类型：发明
国别省市：江苏,32