一种基于曲面全变差的低剂量CT投影域去噪及重建方法技术

本发明专利技术涉及一种基于曲面全变差的低剂量CT投影域去噪及重建方法，步骤一：根据投影域数据噪声的产生机制，建立优化泛函数据项Ψ；建立包括基于l

【技术实现步骤摘要】
一种基于曲面全变差的低剂量CT投影域去噪及重建方法


[0001]本专利技术涉及机CT图像处理方法领域，更具体地，涉及一种基于曲面全变差的低剂量CT投影域去噪及重建方法。

技术介绍

[0002]计算机断层扫描(ComputedTomography,CT)能够以非侵入的方式得到人体内部的解剖图像，清晰地显示各部位的结构组织及病灶，帮助医生进行临床诊断，是一种广泛应用于临床检查的成像技术。然而，CT扫描过程中需要利用X射线穿透人体来采集各个角度的投影数据。许多研究表明，CT扫描过程中产生的X射线辐射会诱发人体的新陈代谢异常甚至有致癌风险，其对人体健康的影响不容忽视。为了应对标准CT扫描所产生的医源性辐射伤害问题，Naidich等人在20世纪90年代提出了低剂量CT(LowDoseCT,LDCT)。LDCT通过降低CT扫描设置中的毫安秒(mAs)设置，即降低通过X射线管的电流或缩短曝光时间来减少X射线的辐射剂量。然而，降低mAs设置会使到达探测器的光子所服从的分布由高斯分布退化为泊松分布。与标准CT相比，LDCT的原始数据生成过程中会受到相对更高水平噪声的污染。
[0003]目前在图像降噪领域被广泛采用的全变差(Total Variation，TV)正则化假设真实图像是分片常数的，如公开号为“CN104574416A”，公开日为2015年4月29日的专利公布文件中公开了一种低剂量能谱CT图像去噪方法。但CT投影域数据通常并不满足这一假设，因此直接将全变差正则化方法应用到投影域降噪问题中容易产生阶梯伪影，影响了去噪的精度。同时，在图像降噪后需要对图像进行重建，目前图像重建方法大部分是基于离散模型的算法，其所求解的图像函数解属于分片常数函数的范畴，所需要的重建时间也比较长。而且由于离散化的过程会引入误差，使用这种方法导出的优化模型天然地存在模型误差。

技术实现思路

[0004]本专利技术为克服上述现有技术中低剂量CT去噪存在阶梯伪影导致去噪精度低和重建时间长的问题，提供一种基于曲面全变差的低剂量CT投影域去噪及重建方法，实现高精度去噪和加快图像重建速度。
[0005]为解决上述技术问题，本专利技术采用的技术方案是：一种基于曲面全变差的低剂量CT投影域去噪及重建方法，包括如下步骤：
[0006]步骤一：根据投影域数据噪声的产生机制，建立优化泛函数据项Ψ；建立包括基于范数的曲面全变差正则项的优化泛函正则项Υ，将优化泛函数据项和优化泛函正则项Υ之和作为优化泛函；
[0007]步骤二：对步骤一中的优化泛函进行求解，得到降噪后的CT投影域数据；
[0008]步骤三：基于连续积分方程模型，根据步骤二中降噪后的CT投影域数据进行图像重建。
[0009]在上述的技术方案中，通过建立包括正则项和数据项的优化泛函并对其进行优化求解，将输入的投影域数据中的噪声去除，实现CT投影域的去噪。其中，曲面全变差正则项
引入曲面的主曲率方向，用主曲率方向代替梯度方向，并结合投影域本身的性质，主曲率方向同时包含了一阶梯度和二阶梯度信息，判断图像某点的变化强度大的方向和平滑方向比一阶梯度方向更为精确，避免图像中出现伪影，令图像去噪更加准确。在图像重建的时候采用连续积分方程模型来设计CT图像重建算法，连续积分方程基于实际物理模型，具有更高的精度和重建速度。
[0010]优选的，在所述步骤一中，优化泛函正则项Υ还包括基于范数的曲波变换系数正则项。根据CT投影域数据曲波系数的稀疏性，基于范数的曲波变换系数正则项能够改善去噪效果。
[0011]优选的，曲面全变差正则项建立方法为：将CT投影域数据视为空间中的曲面，约束曲面在主曲率方向上的跳变，将主曲率方向上梯度的范数作为曲面全变差正则项，具体表示为：
[0012][0013][0014]式中，Y为带噪声的输入投影域数据；B为交错差分算子；S为主曲率方向矩阵；为横向差分算子；为纵向差分算子；L是lipschitz连续系数；M为Y的行数。
[0015]具体的，定义主曲率方向矩阵表示为：
[0016][0017]式中，v1和v2分别为Y的最大主曲率和最小主曲率方向；D
Yρ
为主曲率绝对差。
[0018]优选的，曲波变换系数正则项建立方法为对CT投影域数据Y进行离散曲波变换，得到系数CY，具体表示为：Υ2(Y)＝||CY||1。
[0019]优选的，步骤一中的优化泛函数据项Ψ具体为：
[0020][0021]其中，Y为带噪声的输入投影域数据；Q为到达探测板的光子数；P为探测板所接收到的光子数；I0为X射线入射强度。
[0022]优选的，步骤一中的优化泛函具体为：
[0023][0024]式中，I0为X射线入射的强度；λ1和λ2为两种正则化项的惩罚系数，取为正实数；C为曲波变换算子；S为主曲率方向算子；B为交替差分算子。
[0025]优选的，所述步骤三中的具体流程为：
[0026]S3.1：建立图像重建的优化泛函如下：
[0027][0028]式中，f是图像重建函数；g为平行束投影域函数；为Radon变换；λ是正则化项系数；Φ是恒大于零且L
‑
lipschitz连续的函数；B表示线性变换；
[0029]S3.2：求解图像重建的优化泛函；
[0030]S3.3：对图像重建的优化泛函进行迭代求解，求得图像重建函数f。
[0031]在连续积分方程模型的基础上，建立图像重建的优化泛函中引入了Radon变换，通过Radon变换的卷积性质对后续的迭代算法进行加速，使得图像重建的速度更加快，重建时间更短。
[0032]优选的，在所述步骤S3.2中，采用邻近梯度法求解图像重建的优化泛函，具体的流程为：
[0033]首先计算的梯度。设则其中是Radon变换的共轭算子。利用的卷积性质与傅里叶变换完成加速计算。
[0034]计算计算计算是f在L2空间上某组基函数，z
i
是基函数的系数，是探测板的探测范围，由此可以解得其中利用卷积性质可以加速的计算，设矩阵其元素由于W是循环矩阵，因此w为矩阵W的第一行，z＝[z1,z2,
…
,z
L
]T
。
[0035]计算计算其中其中
[0036]表示g(θ,u)关于u分量的傅里叶变换。
[0037]优选的，在所述步骤S3.3中，采用迭代公式对图像重建的优化泛函进行迭代求解，迭代公式具体为：
[0038][0039]其中，表示傅里叶变换；表示反傅里叶变换；prox
Φ
是Φ的临近梯度算子；L是Φ的lipschitz连续系数；e表示向量按元素相乘；b
(k+1)
表示第k+1步迭代的中间结果，z
(k+1)
表示第k+1步的图像函数。
[0040]优选的，在所述步骤二中，使用带有随机步长的ADMM优化算法对步骤一中的优化泛函进行求解。采用随机步长可以加快优化算法的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于曲面全变差的低剂量CT投影域去噪及重建方法，其特征在于，包括如下步骤：步骤一：根据投影域数据噪声的产生机制，建立优化泛函数据项Ψ；建立包括基于l
1/2
范数的曲面全变差正则项的优化泛函正则项Υ，将优化泛函数据项和优化泛函正则项Υ之和作为去噪优化泛函；步骤二：对步骤一中的去噪优化泛函进行求解，得到降噪后的CT投影域数据；步骤三：基于连续积分方程模型，根据步骤二中降噪后的CT投影域数据进行图像重建。2.根据权利要求1所述的一种基于曲面全变差的低剂量CT投影域去噪及重建方法，其特征在于，在所述步骤一中，优化泛函正则项γ还包括基于l1范数的曲波变换系数正则项。3.根据权利要求2所述的一种基于曲面全变差的低剂量CT投影域去噪及重建方法，其特征在于，建立曲面全变差正则项的方法为将CT投影域数据视为空间中的曲面，约束曲面在主曲率方向上的跳变，将主曲率方向上梯度的范数作为曲面全变差正则项；具体表示为：示为：式中，Y为带噪声的输入投影域数据；B为交错差分算子；S为主曲率方向矩阵；为横向差分算子；为纵向差分算子；L是lipschitz连续系数；M为Y的行数。4.根据权利要求3所述的一种基于曲面全变差的低剂量CT投影域去噪及重建方法，其特征在于，建立曲波变换系数正则项的方法为对CT投影域数据Y进行离散曲波变换，得到系数CY，将曲波变换系数正则项表示为：γ2(Y)＝||CY||1。5.根据权利要求4所述的一种基于曲面全变差的低剂量CT投影域去噪及重建方法，其特征在于，步骤一中的优化泛函数据项Ψ具体为：其中，Y为带噪声的输入投影域数据；Q为到达探测板的光子数；P为探测板所接收到的光子数；I0为X射线入射强度。6.根据权利要求5所述的一种基于曲面全...

【专利技术属性】
技术研发人员：江颖，刘婷，刘伟锋，吴锐帆，
申请(专利权)人：中山大学，
类型：发明
国别省市：