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

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

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


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

技术介绍

[0002]计算机断层成像(Computed Tomography，CT)技术已被广泛应用于临床的诊断，但大量辐射会对人体造成不可逆的伤害，因此低剂量CT图像重建在临床上具有重要的研究价值。现有的重建算法包括滤波反投影算法、迭代重建算法和统计迭代重建算法。滤波反投影算法(FBP算法)对数据完备性要求较高，在低剂量的情况下无法产生理想的诊断图像质量；代数迭代重建和同时代数迭代重建等的迭代算法能够在一定程度上改善这一问题，但由于缺乏额外的先验知识，很难获得满意的结果；在重建过程中引入图像先验信息的统计迭代算法已被证明能够更好的去除因降低射线剂量而产生的噪声和伪影，其相关的目标函数一般由数据保真项和正则项两部分组成。Candes等人提出的压缩感知(Compressed Sensing,CS)理论，已经证明了可以使用远低于奈奎斯特采样速率的数据来精确恢复原始信号。随着CS的兴起，基于稀疏表示的图像重建方法得到了快速的发展，全变分算法利用梯度最小化来达到去噪的目的，但它假设图像是分段平滑的，导致重建图像出现阶梯伪影；基于字典学习的图像重建方法从图像中学习字典，能够得到比一般方法更好的稀疏表示，但这种基于图像块的方法会使特征表示不具备平移不变性，且字典原子高度冗余会导致重建图像的边界出现块状聚合伪影。
[0003]与传统字典学习相同，卷积字典学习具有明确的物理意义，且将平移不变性模型化于目标函数中，使用卷积滤波器和稀疏特征图对整幅图像进行稀疏分解，能够快速有效的解决基于图像分块的字典学习中的边界聚合伪影。
[0004]受深度学习的启发，一些基于卷积神经网络的低剂量CT图像重建方法被相继提出。卷积神经网络通过训练过程使得重建图像和标签图像之间的误差最小化，网络模型具有强大的特征提取能力，但缺点在于所需的数据量较大且难以让人理解内部参数优化的逻辑性。基于此，本专利技术提出了一种结合卷积神经网络和卷积字典学习的低剂量CT重建方法，凭借卷积神经网络强大的特征提取能力和卷积字典学习的可解释性，来解决低剂量下CT重建图像质量较差的问题。所提方法能够有效的去除了噪声和伪影，且能够更好的保留图像细节。
[0005]传统的基于字典学习和稀疏表示的低剂量CT图像重建方法，是对图像进行分块并逐块处理，这导致重建图像出现边界聚合伪影，且大多图像重建方法都使用手工先验(如L0先验)，这也导致学习到的字典对图像结构的表达能力有限。

技术实现思路

[0006]本专利技术的目的在于提供一种结合卷积神经网络和卷积字典学习的低剂量CT图像重建方法。该方法利用网络自适应的学习模型的先验，使模型能够更好的适用于重建图像。
卷积字典学习是基于整幅图像的，能有效的解决边界聚合伪影问题，且可解释性使得网络优化更直观。两者的结合能够更好的去除低剂量图像中的噪声和伪影。
[0007]本专利技术提供的一种结合卷积神经网络和卷积字典学习的低剂量CT图像重建方法包括如下步骤：
[0008]步骤一：选取一定数量的正常剂量CT数据集，该数据集中包括训练集和测试集；
[0009]步骤二：数据初始化，采集待重建CT图像的投影数据，通过在正常剂量CT投影数据上添加噪声来模拟低剂量CT图像的投影数据，再对噪声投影进行滤波反投影处理，得到初始的重建图像，通过稀疏特征图初始化网络来得到初始卷积稀疏特征图，并将与所述卷积稀疏特征图相对应的卷积字典原子初始化为全零；
[0010]步骤三：模型训练，将总的目标函数分解为三个子问题，通过在每次迭代中交替的更新重建图像、卷积稀疏特征图和卷积字典来训练网络模型；
[0011]步骤四：重复步骤三，直至达到设定的迭代次数，即stage数；
[0012]步骤五：对重建图像计算损失，更新模型，待损失稳定后保存模型；
[0013]步骤六：将测试集输入训练好的模型中，获得重建后的CT图像。
[0014]本专利技术提出的一种结合卷积神经网络和卷积字典学习的低剂量CT图像重建方法，利用了卷积神经网络和卷积字典学习的优势，将网络强大的特征提取能力与卷积字典学习的可解释性有效结合，通过网络为每个图像自适应的学习卷积字典以及先验。
[0015]该方法凭借卷积字典学习和CT重建原理构建目标函数，并基于数学优化理论求解和构建网络模块，在模型性能维护时逻辑性更强，在重建效果方面不仅能够去除图像中的噪声和伪影，还能够保留更多的图像细节。
附图说明
[0016]图1是本专利技术结合卷积神经网络和卷积字典学习的低剂量CT图像重建方法的流程图。
[0017]图2是本专利技术每轮迭代时更新重建图像，更新卷积字典以及稀疏特征图的结构示意图。
[0018]图3是本专利技术中各个模块以及网络的结构示意图。
具体实施方式
[0019]为了更好地理解本专利技术的技术方案，以下结合附图对本专利技术技术方案进行详细说明。
[0020]步骤一：选取合适的CT数据集。
[0021]选取一定数量的正常剂量CT数据集，本实验选取的数据集共有800张，尺寸大小为512
×
512，其中765张作为训练集，35张作为测试集。
[0022]步骤二：数据初始化：采集待重建图像的投影数据。通过在正常剂量的投影数据上添加噪声来模拟低剂量的投影数据，再对添加噪声后的投影数据进行滤波反投影处理，来得到初始的重建图像。通过“稀疏特征图初始化网络”来得到初始的卷积稀疏特征图，并将卷积字典初始化为全零。
[0023]步骤201)：研究表明，CT投影数据噪声产生的两个主要来源是X射线量子噪声和系
统电子噪声。探测器接收到的对数变化前的投影可以用统计独立的泊松分布加上统计独立的高斯分布来描述，即：
[0024][0025]因此，通过给正常剂量的投影添加服从泊松分布加上高斯分布的噪声来模拟低剂量的投影数据P。其中，λ为光子数量的期望值，m
e
是探测器收集的辐射强度的测量值，是电子噪声的均值(通常m
e
为0)，为电子噪声的方差。
[0026]步骤202)：数据初始化。通过滤波反投影(FBP)算法得到初始化的重建图像Y0。将初始图像Y0和估计的噪声方差σ
i
(即图像中加入两种类型的噪声后，整体的噪声方差)带入卷积稀疏编码的初始模块中(对应图3中卷积稀疏特征图初始化网络)，得到m个大小为N
×
N的初始卷积稀疏特征图X＝{X1,X2,
…
,X
m
}，将X
i
与相对应的卷积字典原子D
i
用全0进行初始化，通常卷积字典原子的尺寸远小于稀疏特征图。其中，噪声方差σ
i
可被定义为：
[0027][0028]其中，I0为入射光子数(表示射线到达身体前的本文档来自技高网...

【技术保护点】

【技术特征摘要】
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
...

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