基于截断数据外推网络的CT局部重建方法
文献发布时间:2023-06-19 10:16:30
技术领域
本发明涉及CT图像局部重建技术领域,具体涉及一种基于截断数据外推网络的CT局部重建方法。
背景技术
计算机断层成像技术(CT)是一种X光显影技术,以其无损、可视化等优点广泛应用于医疗诊断、工业无损探伤等领域。虽然CT技术已经发展得相当成熟,但受到现有成像方法的制约,在高分辨率成像、大物体成像和减少辐射剂量等方面仍然存在一定的问题。
一方面,医疗领域中通常只对人体的某个发生病变的器官或部位感兴趣,另一方面,实际的工程应用中,由于物体结构、尺寸的约束,不能进行全角度的扫描。以上情况均会造成投影数据的截断,若使用传统重建算法如滤波反投影算法(FBP)进行重建,图像就会出现明显的截断伪影。为消除截断伪影,可以采用投影数据外推延拓的方法对投影数据进行补全,再用FBP算法进行重建。
现有数据外推方法主要是通过插值或空间变换迭代来实现。例如,常量外推法:以边界点的投影值作为缺失的数据值向两侧外延数据;曲线拟合法:利用最小二乘法对采样点进行曲线拟合,再利用这条曲线外推出其他投影值。这些外推方法只利用投影边缘数据进行补全,没有有效利用投影数据间的冗余性,其精确度较低。
发明内容
本发明为解决CT图像局部重建出现截断伪影的问题,提供一种基于截断数据外推网络的CT局部重建方法,通过N张同一物体的CT图像及截断图像获得训练样本数据集,利用基于U-net网络的截断数据外推网络通过训练样本数据集进行深度学习,得到一个能够较好的补全截断投影数据的神经网络,将待重建截断图像输入截断数据外推网络,对待重建截断图像进行补全,最终在得到补全后的完整数据后进行FBP重建,从而解决截断伪影的问题。
本发明的技术方案是:
基于截断数据外推网络的CT局部重建方法,所述方法包括:
步骤1:基于投影图像特征,制作匹配的训练样本数据集;
步骤2:基于投影图像特点,计算截断数据外推网络损失;
步骤3:利用截断数据外推网络损失作为约束对截断数据外推网络的构建进行优化、并通过训练样本数据集对截断数据外推网络进行训练;
对特征图进行融合完成对截断区域信息的补全,获得修补后的投影正弦图;
步骤4:对截断数据外推网络补全后的投影正弦图进行FBP重建。
进一步地,所述步骤1具体包括:
步骤1.1:获取N张同类物体的不同角度的完整CT图像、并对N张所述完整CT图像分别进行随机像素的上下截断处理,得到N张与完整CT图像对应的截断图像;
步骤1.2:将截断图像的被截断部分投影数据置零,获得截断图像的投影正弦图,同时获得完整CT图像的投影正弦图;
步骤1.3:构建完整CT图像的投影正弦图与截断图像的投影正弦图一一对应的图像对、并将所述图像对整合建立训练样本数据集。
进一步地,所述步骤2具体包括:
步骤2.1:在截断处理区域将L2损失函数作为外推网络的损失约束条件,按照公式(1)进行截断数据外推网络损失计算:
其中f和
进一步地,所述步骤3具体包括:
步骤3.1基于公式(1)约束条件,对应调整截断数据外推网络的层数设置和模块组成,进行截断数据外推网络结构优化;
步骤3.2:将训练样本数据集上传截断数据外推网络,利用训练样本数据集进行网络的训练,通过深度学习完善截断数据外推网络;
步骤3.3:获取待重建的截断图像投影正弦图,并将其输入截断数据外推网络进行数据外推,通过截断数据补全得到补全投影正弦图。
进一步地,步骤3.1中所述层数设置包括将原始的池化层替换为步长为2的卷积操作来进行数据降维处理,提升截断数据外推网络细节特征的提取能力。
进一步地,步骤3.1中所述模块组成包括多个降维模块及升维模块,多个所述降维模块及升维模块均包含两个卷积核大小为3×3且步长为1的卷积操作。
进一步地,所述步骤3.3具体包括:
步骤3.3.1:在降维模块中,通过步长为2,卷积核大小为3×3的卷积降维过程,实现图像特征的提取,截断数据外推网络对训练样本数据集进行映射,通过多层次的非线性映射逐步提取非截断区域信息对应的特征图;
步骤3.3.2:利用Relu激活函数,通过升维模块完成特征图融合,对截断区域信息进行补全,获得补全投影正弦图。
通过上述技术方案,本发明的有益效果为:
本发明提出了一种基于截断数据外推网络的CT局部重建方法,针对投影图像的特点,将截断数据外推网络进行了优化,通过深度学习将截断数据补全,最终将网络输出结果用FBP算法进行重建得到完整的CT图像。
与直接利用截断投影数据重建得到的图像相比,利用截断数据外推网络补全,再进行FBP算法重建的CT图像,其截断伪影得到有效的抑制,图像更加清晰。从数据指标上看,RMSE(均方根误差)与NMAD(归一化平均绝对距离)都大大降低,NMAD和RMSE分别下降44.24%、43.18%左右。重建结果和上述指标表明,该方法可以有效处理数据截断引起的伪影,提高重建图像的质量。
附图说明
图1基于截断数据外推网络的CT局部重建的截断数据外推网络结构图。
图2基于截断数据外推网络的CT局部重建的流程框图。
图3基于截断数据外推网络的CT局部重建的结果对比图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步说明:
如图1~2所示,基于截断数据外推网络的CT局部重建方法,所述方法包括:
步骤1:基于投影图像特征,制作匹配的训练样本数据集;
步骤2:基于投影图像特点,计算截断数据外推网络损失;
步骤3:利用截断数据外推网络损失作为约束对截断数据外推网络的构建进行优化、并通过训练样本数据集对截断数据外推网络进行训练;
对特征图进行融合完成对截断区域信息的补全,获得修补后的投影正弦图;
步骤4:对截断数据外推网络补全后的投影正弦图进行FBP(Filter back-projection滤波反投影)重建。
为使截断数据外推网络得到足够的学习数据以保证截断数据外推网的补全效果,所述步骤1具体包括:
步骤1.1:利用CT扫描设备进行360°扫描获取同类物体的N个完整CT图像、并对N个完整CT图像分别进行随机像素的上下截断处理得到N个与完整CT图像对应的截断图像;
步骤1.2:将截断图像的被截断部分投影数据置零,获得截断图像的投影正弦图,同时获取完整CT图像的投影正弦图;
步骤1.3:构建完整CT图像的投影正弦图与截断图像的投影正弦图一一对应的图像对、并将所述图像对整合建立训练样本数据集。
为提高补全图像匹配效果,所述步骤2具体包括:
步骤2.1:在截断处理区域将L2损失函数作为外推网络的损失约束条件,按照公式(1)进行截断数据外推网络损失计算:(由于投影正弦图像分布均匀,不存在较大的偏差值,在本方法中选择L2损失函数作为网络的损失约束条件。L2损失函数可以在整个网络补偿过程中有效防止偏差较大的离群值出现,以确保外推数据的可靠性。)
其中f和
为实现基于截断数据外推网络的截断投影正弦图外推补全,所述步骤3具体包括:
步骤3.1基于公式(1)约束条件,对应调整截断数据外推网络的层数设置和模块组成,进行截断数据外推网络结构优化;
步骤3.2:将训练样本数据集上传截断数据外推网络,利用训练样本数据集进行网络的训练,通过深度学习完善截断数据外推网络;
步骤3.3:获取待重建的截断图像投影正弦图,并将其输入截断数据外推网络进行数据外推,通过截断数据补全得到补全投影正弦图。
为优化补全效果,步骤3.1中所述层数设置包括将原始的池化层替换为步长为2的卷积操作来进行数据降维处理,提升截断数据外推网络细节特征的提取能力。步骤3.1中所述模块组成包括多个降维模块及升维模块,多个所述降维模块及升维模块均包含两个卷积核大小为3×3且步长为1的卷积操作。
为进一步优化补全效果,所述步骤3.3具体包括:
步骤3.3.1:在降维模块中,通过步长为2,卷积核大小为3×3的卷积降维过程,实现图像特征的提取,截断数据外推网络对训练样本数据集进行映射,通过多层次的非线性映射逐步提取非截断区域信息对应的特征图;
步骤3.3.2:利用Relu激活函数,通过升维模块完成特征图融合,对截断区域信息进行补全,获得补全投影正弦图。
为证明本发明效果,进行如下实验设置:
在本实施例中,选择临床上的腹部CT图像作为模型,实验时选取了2121张完整的不同角度的腹部CT图像,然后将2121张完整的不同角度的腹部CT图像做投影数据处理,得到完整CT图像的投影正弦图(正弦图大小为512×360),然后对投影数据的正弦图模拟实际截断情况随机上下截取80-150,得到2121张不同截断部位的截断图像的投影正弦图,将其中20张不同截断部位的截断图像的投影正弦图取出作为测试样本,并将与之对应的20张完整CT图像的投影正弦图作为最终比对样本;
将其余的2101张完整CT图像的投影正弦图及2101张截断图像的投影正弦图按步骤1形成训练样本数据集,输入截断数据外推网络进行训练;
在本实施例中,所述截断数据外推网络采用8层的U-net模型的神经网络,选择L2损失函数作为截断数据外推网络的约束损失,保证补全数据的有效性和可靠性。根据损失函数的计算值,利用Adam优化器不断优化神经网络参数,使其逼近或达到最优值,从而使损失函数值最小,最终得到一个能够较好的补全截断投影数据的神经网络。
试验时运用AMAX工作站上的TensorFlow(版本1.4.0)进行截断数据外推网络的训练和测试。该工作站具有两个Intel Xeon E5-2640 v4 CPU 2.4GHz和64GB内存。
在完成截断数据外推网络的构建、学习及优化后,将20张测试样本输入截断数据外推网络,运用截断数据外推网络对20张测试样本进行外推补全,之后进行FBP重建。
其补全后的腹部CT图像与比对样本的比对图像如图3所示。图中(a)表示完整CT图像的投影正弦图,(b)表示截断图像的投影正弦图,(c)表示补全后的CT图像的投影正弦图,(d)表示完整CT图像,(e)表示截断CT图像,(f)表示补全后的CT图像。
以上所述之实施例,只是本发明的较佳实施例而已,并非限制本发明的实施范围,故凡依本发明专利范围所述的构造、特征及原理所做的等效变化或修饰,均应包括于本发明申请专利范围内。
- 基于截断数据外推网络的CT局部重建方法
- 基于非对称横向双边截断投影数据的扇束CT重建方法