基于双残差网络的低剂量CT牙齿图像去噪方法技术

本发明专利技术公开一种基于双残差网络的低剂量CT牙齿图像去噪方法，首先训练辅助网络，并保存模型；利用经过双残差网络提取的特征空间图之间的相似关系，计算损失函数值，通过合并去噪网络和双残差网络之间的损失函数值，进而更新参数以辅助去噪网络的训练。该方法主要包括数据加载模块、辅助网络训练模块、双残差网络去噪模块、测试模块。本发明专利技术对低剂量CT牙齿图像去噪有着显著的效果，能够保留图像的一些细节特征，具有重要的应用价值。

Denoising method of low dose CT tooth image based on double residual network

【技术实现步骤摘要】
基于双残差网络的低剂量CT牙齿图像去噪方法
本专利技术方法涉及一种图像去噪方法，尤其是一种基于双残差网络的低剂量CT牙齿图像去噪方法。
技术介绍
近年来，口腔CT成像技术越来越广泛的应用在口腔、牙齿的诊断与治疗过程中。口腔CT利用采集仪围绕待检测人体组织区域螺旋式地采集数据，可以获得各方位断层图像，进而被用于辅助牙齿手术前测量、口腔炎症、肿瘤与其他口腔疾病的诊断。采用低剂量CT设备对人体的辐射剂量约是传统CT扫描仪的1/5，因此采用低剂量设备可以降低辐射剂量，让广大患者减少不必要的心理负担，更加安心地接受CT扫描。但是，在选择降低剂量CT扫描的同时，图像的清晰度也有所下降，同时射线穿透性降低会因散射、干涉现象导致大量的噪声产生。近年来图像处理、深度学习技术日益成熟，基于深度学习的图像去噪技术已经可以获得比较优异的去噪效果。然而大部分基于深度学习的去噪方法是采用训练单神经网络来实现去噪任务的，但是恢复微小细节的能力尚有待提高，而牙齿图像细小纹理正是医师判断牙齿是否存在内部裂纹的重要依据。为此，低清晰度，高噪声的低剂量CT图像，只能用于疾病的早期预防。对于诸如牙齿断裂、部分病灶的进一步确诊还是需要高分辨CT甚至是增强CT来确认的。
技术实现思路
本专利技术是为了解决现有技术所存在的上述技术问题，提供一种基于双残差网络的低剂量CT牙齿图像去噪方法。本专利技术的技术解决方案是：一种基于双残差网络的低剂量CT牙齿图像去噪方法，按照如下步骤进行：01部分，取BSD500数据集的灰度图像，在进入网络模型之前进行预处理操作，具体步骤如下：步骤C011：导入BSD500数据集的灰度图像，分别为500张预处理图像和500张真实图像，各500张灰度图像中432张做训练集，68张做测试集，预处理图像数据集记为Pre_Image，用于验证的真实图像数据集记为Real_Image；500张预处理图像分别记为Pre_Image1，Pre_Image2，…，Pre_Image500；500张真实图像分别记为Real_Image1，Real_Image2，…，Real_Image500；步骤C012：将500张预处理图像和500张真实图像切块，块的大小为48*48像素，将预处理图像与真实图像所切的块分别记为预处理图像块Pre_P1，Pre_P2,…,Pre_Pn和真实图像块Real_P1，Real_P2,…,Real_Pn；步骤C013：将预处理图像块Pre_P1，Pre_P2,…,Pre_Pn分别加入σ值为50的高斯噪声，记为噪声图像块Noise_P1，Noise_P2,…,Noise_Pn；步骤C014：每次随机选取16000个噪声图像块和16000个真实图像块输入网络，其中噪声图像块记作Noise_Pi，真实图像块记作Real_Pi，1≤i≤16000；02部分，训练辅助网络20次，为双残差网络去噪提供数据，具体步骤如下：步骤C021：取单通道的真实图像块Real_Pi，输进网络；步骤C022：开始训练M＝16000对应的线性关系模型Model1，约定变量epoch_pro为辅助网络循环训练的次数，初始为0；所述模型Model1的目标函数定义为：式中为真实图像块Real_Pi，L(·)为损失函数，f(·|Φ)表示带有Φ参数化的辅助网络，λR(Φ)是一个带有参数Φ和λ的正则化项，所述λ>0；步骤C0221：进入浅层特征提取层，记为Conv_pro1层；Conv_pro1层由64个3*3大小的滤波器组成，真实图像块Real_Pi经过Conv_pro1层后，得到64个通道的浅层特征空间图，记为Real_proPi_fm1；步骤C0222：进入深层特征提取层，记为Deep_proRes2层；Deep_proRes2层中包含5个残差组和一个卷积层以及一个长跳跃连接；每一个残差组中包含10个带有通道注意机制的残差块和一个短跳连接；其中残差块是由两个卷积层和Relu激活层以及跳跃连接构成；浅层特征空间图Real_proPi_fm1经过Deep_proRes2层后，得到64个通道的深层特征空间图，记为Real_proPi_fm2；步骤C0223：进入重构层，记为Conv_pro3层；Conv_pro3层由1个3*3大小的滤波器组成，深层特征空间图Real_proPi_fm2经过Conv_pro3层后，得到重构特征空间图，记为Real_proPi_fm3；步骤C023：求重构特征空间图Real_proPi_fm3和真实灰度图像块Real_Pi之间的差异，通过L1Loss计算方法可得相应的损失函数值，记为loss_pro，置epoch_pro＝epoch_pro+1；若loss_pro值小于0.01或者训练次数epoch_pro达到20，网络停止训练，保存训练模型model_pro.pt文件；否则将loss_pro值反向传播，通过ADAM优化算法重新更新参数，再一次进入步骤C022，循环重复训练；03部分，利用双残差网络中得到的特征空间图的差异，进行去噪处理，具体步骤如下：步骤C031：取单通道的噪声图像块Noise_Pi和单通道的真实图像块Real_Pi；步骤C032：构建与训练辅助网络Model1结构一致的新辅助网络模型Model1*，将新辅助网络模型Model1*的参数替换为020部分训练辅助网络Model1得到的model_pro.pt模型和参数，并且固定该新辅助网络模型Model1*的参数不再更新；步骤C033：开始训练去噪网络M＝16000对应的线性关系模型Model，约定变量epoch为双残差网络循环训练的次数，初始化为0；所述模型Model的目标函数定义为：式中为真实图像块Real_Pi，为噪声图像块Noise_Pi，Ω为层索引，L(·)为损失函数，为通过第j个隐藏层得到的特征空间图的损失函数，fj(·)表示第j个隐藏层之前的网络，λR(θ)是一个带有参数θ和λ的正则化项，所述λ>0；步骤C0331：将噪声图像块Noise_Pi输入去噪网络Model；步骤0332：进入浅层特征提取层，记为Conv1层；Conv1层由64个3*3大小的滤波器组成，噪声图像块Noise_Pi经过Conv1层后，得到64个通道的浅层特征空间图，记为Noise_Pi_fm1；步骤0333：进入深层特征提取层，记为Deep_Res2层；Deep_Res2层中包含5个残差组和一个卷积层，以及一个长跳跃连接；每一个残差组中包含10个带有通道注意机制的残差块和一个短跳连接；其中残差块是由两个卷积层和Relu激活层以及跳跃连接构成；浅层特征空间图Noise_Pi_fm1经过Deep_Res2层后，得到64个通道的深层特征空间图，记为Noise_Pi_fm2；步骤0334：进入重构层，记为Conv3层；Conv3层由1个3*3大小的滤波器组成，深层特征空间图Noise_Pi_fm2经过Conv3层后本文档来自技高网...

【技术保护点】
1.一种基于双残差网络的低剂量CT牙齿图像去噪方法，其特征在于按照如下步骤进行：/n01部分，取BSD500数据集的灰度图像，在进入网络模型之前进行预处理操作，具体步骤如下：/n步骤C011：导入BSD500数据集的灰度图像，分别为500张预处理图像和500张真实图像，各500张灰度图像中432张做训练集，68张做测试集，预处理图像数据集记为Pre_Image，用于验证的真实图像数据集记为Real_Image；500张预处理图像分别记为Pre_Image1，Pre_Image2，…，Pre_Image500；500张真实图像分别记为Real_Image1，Real_Image2，…，Real_Image500；/n步骤C012：将500张预处理图像和500张真实图像切块，块的大小为48*48像素，将预处理图像与真实图像所切的块分别记为预处理图像块Pre_P

【技术特征摘要】
1.一种基于双残差网络的低剂量CT牙齿图像去噪方法，其特征在于按照如下步骤进行：
01部分，取BSD500数据集的灰度图像，在进入网络模型之前进行预处理操作，具体步骤如下：
步骤C011：导入BSD500数据集的灰度图像，分别为500张预处理图像和500张真实图像，各500张灰度图像中432张做训练集，68张做测试集，预处理图像数据集记为Pre_Image，用于验证的真实图像数据集记为Real_Image；500张预处理图像分别记为Pre_Image1，Pre_Image2，…，Pre_Image500；500张真实图像分别记为Real_Image1，Real_Image2，…，Real_Image500；
步骤C012：将500张预处理图像和500张真实图像切块，块的大小为48*48像素，将预处理图像与真实图像所切的块分别记为预处理图像块Pre_P1，Pre_P2,…,Pre_Pn和真实图像块Real_P1，Real_P2,…,Real_Pn；
步骤C013：将预处理图像块Pre_P1，Pre_P2,…,Pre_Pn分别加入σ值为50的高斯噪声，记为噪声图像块Noise_P1，Noise_P2,…,Noise_Pn；
步骤C014：每次随机选取16000个噪声图像块和16000个真实图像块输入网络，其中噪声图像块记作Noise_Pi，真实图像块记作Real_Pi，1≤i≤16000；
02部分，训练辅助网络20次，为双残差网络去噪提供数据，具体步骤如下：
步骤C021：取单通道的真实图像块Real_Pi，输进网络；
步骤C022：开始训练M＝16000对应的线性关系模型Model1，约定变量epoch_pro为辅助网络循环训练的次数，初始为0；
所述模型Model1的目标函数定义为：式中为真实图像块Real_Pi，L(·)为损失函数，f(·|Φ)表示带有Φ参数化的辅助网络，λR(Φ)是一个带有参数Φ和λ的正则化项，所述λ>0；
步骤C0221：进入浅层特征提取层，记为Conv_pro1层；Conv_pro1层由64个3*3大小的滤波器组成，真实图像块Real_Pi经过Conv_pro1层后，得到64个通道的浅层特征空间图，记为Real_proPi_fm1；
步骤C0222：进入深层特征提取层，记为Deep_proRes2层；Deep_proRes2层中包含5个残差组和一个卷积层以及一个长跳跃连接；每一个残差组中包含10个带有通道注意机制的残差块和一个短跳连接；其中残差块是由两个卷积层和Relu激活层以及跳跃连接构成；浅层特征空间图Real_proPi_fm1经过Deep_proRes2层后，得到64个通道的深层特征空间图，记为Real_proPi_fm2；
步骤C0223：进入重构层，记为Conv_pro3层；Conv_pro3层由1个3*3大小的滤波器组成，深层特征空间图Real_proPi_fm2经过Conv_pro3层后，得到重构特征空间图，记为Real_proPi_fm3；
步骤C023：求重构特征空间图Real_proPi_fm3和真实灰度图像块Real_Pi之间的差异，通过L1Loss计算方法可得相应的损失函数值，记为loss_pro，置epoch_pro＝epoch_pro+1；若loss_pro值小于0.01或者训练次数epoch_pro达到20，网络停止训练，保存训练模型model_pro.pt文件；否则将loss_pro值反向传播，通过ADAM优化算法重新更新参数，再一次进入步骤C022，循环重复训练；
03部分，利用双残差网络中得到的特征空间图的差异，进行去噪处理，具体步骤如下：
步骤C031：取单通道的噪声图像块Noise_Pi和单通道的真实图像块Rea...

【专利技术属性】
技术研发人员：傅博，王丽妍，杜飞飞，刘芳菲，
申请(专利权)人：辽宁师范大学，
类型：发明
国别省市：辽宁;21