一种结合神经网络和卷积字典学习的低剂量CT图像重建方法技术

本发明专利技术涉及一种结合卷积神经网络和卷积字典学习的CT图像重建方法。该方法利用网络自适应的学习模型的先验，使模型能够更好的适用于重建图像，卷积字典学习是基于整幅图像的，能有效的解决边界聚合伪影问题，且可解释性使得网络优化更直观。两者的结合能够有效的去除低剂量图像中的噪声和伪影，能够更好的保留图像细节。像细节。像细节。

【技术实现步骤摘要】
一种结合神经网络和卷积字典学习的低剂量CT图像重建方法


[0001]本专利技术涉及一种低剂量CT图像重建方法，尤其是一种结合了神经网络和卷积字典学习的低剂量CT图像重建方法。

技术介绍

[0002]计算机断层成像(Computed Tomography，CT)技术已被广泛应用于临床的诊断，但大量辐射会对人体造成不可逆的伤害，因此低剂量CT图像重建在临床上具有重要的研究价值。现有的重建算法包括滤波反投影算法、迭代重建算法和统计迭代重建算法。滤波反投影算法(FBP算法)对数据完备性要求较高，在低剂量的情况下无法产生理想的诊断图像质量；代数迭代重建和同时代数迭代重建等的迭代算法能够在一定程度上改善这一问题，但由于缺乏额外的先验知识，很难获得满意的结果；在重建过程中引入图像先验信息的统计迭代算法已被证明能够更好的去除因降低射线剂量而产生的噪声和伪影，其相关的目标函数一般由数据保真项和正则项两部分组成。Candes等人提出的压缩感知(Compressed Sensing,CS)理论，已经证明了可以使用远低于奈奎斯特采样速率的数据来精确恢复原始信号。随着C本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种结合卷积神经网络和卷积字典学习的低剂量CT图像重建方法，该方法包括如下步骤：步骤一：选取一定数量的正常剂量CT数据集，该数据集中包括训练集和测试集；步骤二：数据初始化，采集待重建CT图像的投影数据，通过在正常剂量CT投影数据上添加噪声来模拟低剂量CT图像的投影数据，再对噪声投影进行滤波反投影处理，得到初始的重建图像，通过稀疏特征图初始化网络来得到初始卷积稀疏特征图，并将与所述卷积稀疏特征图相对应的卷积字典原子初始化为全零；步骤三：模型训练，将总的目标函数分解为三个子问题，通过在每次迭代中交替的更新重建图像Y、卷积稀疏特征图X和卷积字典D来训练网络模型；步骤四：重复步骤三，直至达到设定的迭代次数，即stage数；步骤五：对重建图像计算损失，更新模型，待损失稳定后保存模型；步骤六：将测试集输入训练好的模型中，获得重建后的CT图像。2.根据权利要求1所述的低剂量CT图像重建方法，其特征在于CT投影数据噪声的两个主要来源是X射线量子噪声和系统电子噪声，探测器接收到的对数变化前的投影可以用统计独立的泊松分布加上统计独立的高斯分布来描述，即：通过给正常剂量的投影添加服从泊松分布加上高斯分布的噪声来模拟低剂量的投影数据P，其中，λ为光子数量的期望值，是探测器收集的辐射强度的测量值，m
e
是电子噪声的均值(通常m
e
为0)，为电子噪声的方差。3.根据权利要求1或2所述的低剂量CT图像重建方法，其特征在于通过滤波反投影(FBP)算法得到初始化的重建图像Y0，将初始图像Y0和估计的噪声方差σ
i
带入卷积稀疏编码的初始模块中，得到m个大小为N
×
N的初始卷积稀疏特征图X＝{X1,X2,
…
,X
m
}，将与相对应的卷积字典原子用全0进行初始化，通常卷积字典原子的尺寸远小于稀疏特征图。其中，噪声方差可被定义为：其中，I0为入射光子数，为投影数据在第i个探测器上的期望值。4.根据权利要求1所述的低剂量CT图像重建方法，其特征在于借助可分离抛物线代理以及半二次方分裂算法来交替最小化所述子问题，通过在每次迭代中交替的更新重建图像Y、卷积稀疏特征图X和卷积字典D来训练网络模型，所述目标函数可表示为：其中，A为系统矩阵，P为投影数据，Y为待重建图像，Σ是一个对角矩阵，其矩阵元素为D
i
表示第i个卷积字典原子，X
i
为字典原子对应的稀疏特征图，φ(X)表示稀疏特征图X的先验，表示字典D上的正则化项，β
X
和β
D
分别为X和D的正则化参数。5.根据权利要求1或4所述的低剂量CT图像重建方法，其特征在于用迭代重建算法对图像进行重建，在第t(t＝1,2,...M)次迭代过程中，将上一次迭代后的结果，即图像Y
(t
‑
1)
，卷积稀疏特征图X
(t
‑
1)
和卷积稀疏编码D
...

【专利技术属性】
技术研发人员：桂蕴嘉，刘祎，谷亚男，颜溶標，
申请(专利权)人：中北大学，
类型：发明
国别省市：