[发明专利]一种低剂量X射线CT图像统计迭代重建方法

摘要

本发明公开一种低剂量X射线CT图像统计迭代重建方法，包括步骤：对CT设备的低剂量X射线CT图像的投影数据y_raw进行图像重建，得到迭代初始图像μ^init；对投影数据y_raw进行数据恢复处理，得到恢复后的投影数据y_restored，对恢复后的投影数据y_restored进行图像重建，得到参考图像μ^ref；以参考图像μ^ref和迭代初始图像_μ^init，按R(μ^init)=φ(μ^init-SRNLM(μ^init))构建边缘保持先验R(μ^init)，其中，φ()为能量势函数，SRNLM(μ^init)为由参考图像μ^ref引导的非局部均值滤波；根据迭代初始图像μ^init边缘保持先验R(μ^init)，利用统计迭代公式进行迭代计算，得到迭代重建图像μ^iter；当对重建图像μ^iter进行迭代计算的结果满足迭代终止条件时，迭代终止，获得低剂量X射线CT图像最终的重建图像。本发明能够有效地去除图像噪声并抑制条形伪影，同时较好地保持图像细节信息。

主权项

一种低剂量X射线CT图像统计迭代重建方法，其特征在于，包括步骤：获取CT设备的系统参数和低剂量X射线CT图像的投影数据y_raw，对投影数据y_raw进行图像重建，得到迭代初始图像μ^init；对投影数据y_raw进行数据恢复处理，得到恢复后的投影数据y_restored，对恢复后的投影数据y_restored进行图像重建，得到参考图像μ^ref；以参考图像μ^ref和迭代初始图像μ^init，按R(μ^init)=φ(μ^init‑SRNLM(μ^init))构建边缘保持先验R(μ^init)，其中，φ()为能量势函数，SRNLM(μ^init)为由参考图像μ^ref引导的非局部均值滤波；根据迭代初始图像μ^init边缘保持先验R(μ^init)，利用统计迭代公式进行迭代计算，得到迭代重建图像μ^iter；当对重建图像μ^iter进行迭代计算的结果满足迭代终止条件时，迭代终止，获得低剂量X射线CT图像最终的重建图像。