用于图像重建的FDK源码

资 源 简 介

用于图像重建的FDK源码

详 情 说 明

FDK算法是滤波反投影(FBP)类重建方法在锥束CT中的经典实现，特别适用于直线扫描轨迹的有限角度重建场景。该算法的核心思想是通过三维空间的加权滤波和反投影操作，从二维投影数据中重建出三维物体图像。

在直线扫描的有限角度场景下，FDK通过以下关键步骤实现重建：首先对采集到的二维投影数据进行加权预处理，补偿不同探测器单元接收的X射线路径差异。接着进行沿探测器行方向的一维滤波，常用的滤波函数包括Ram-Lak、Shepp-Logan等，这一步能有效抑制反投影过程中产生的星状伪影。滤波后的投影数据会经过锥束几何特有的加权因子调整，最后通过三维反投影运算重建出物体图像。

值得注意的是，有限角度重建会导致数据不完备问题，FDK通过引入Parker加权等方法来缓解由此产生的重建伪影。算法实现时需要特别注意锥束几何中的射线采样均匀性，这直接影响最终重建图像的质量。现代改进算法常在经典FDK基础上加入正则化约束或深度学习方法，以进一步提升有限角度下的重建精度。