[发明专利]一种基于边界信息扩散的X射线CL图像重建算法

摘要

本发明公开一种基于边界信息扩散的X射线CL图像重建算法，步骤如下：输入变量；初始化；假设已得到估计图像u(k)，以u(k)为初值，利用CL扫描数据集p，更新估计图像u(k+1/2)＝RG(p,u(k))；利用图像u(k+1/2)的边界信息扩散修正估计图像，得到u(k+1)＝P(u(k+1/2))；判断相邻两次迭代图像间的差别是否小于给定阈值，或是否达到迭代次数上限N，满足则终止迭代；否则，以u(k+1)为初值，开始新一轮迭代，直至满足则终止，基于边界信息扩散的X射线CL图像重建算法结束。本发明算法利用重建出的图像边界信息，通过迭代的方式，逐步恢复出图像内部灰度值，能够有效降低重建图像层间信息混叠，提高重建图像的层间分辨率，进而提高CL系统的实用性。

主权项

一种基于边界信息扩散的X射线CL图像重建算法，其特征在于，步骤如下：步骤1、输入变量：CL扫描数据集p，CL扫描几何参数集G；步骤2、初始化：初始估计图像u(0)，迭代终止阈值ε或迭代次数上限N；步骤3、假设已得到估计图像u(k)，以u(k)为初值，利用CL扫描数据集p，更新估计图像u(k+1/2)＝RG(p,u(k))，其中RG表示扫描几何参数集G相关的图像重建算子；步骤4、利用图像u(k+1/2)的边界信息扩散修正估计图像，得到u(k+1)＝P(u(k+1/2))，其中P表示边界信息扩散修正算子；步骤5、判断相邻两次迭代图像间的差别是否小于给定阈值，即||u(k+1)‑u(k)||≤ε，或是否达到迭代次数上限N，满足则终止迭代；否则，转步骤3，直至相邻两次迭代图像间的差别小于给定阈值，即||u(k+1)‑u(k)||≤ε，或达到迭代次数上限N，则终止，基于边界信息扩散的X射线CL图像重建算法结束。