割线体投影器:三维锥束CT算子构建的新方法

Cutting Voxel Projector a New Approach to Construct 3D Cone Beam CT Operator

摘要 Abstract

本文介绍了一类新型的三维锥束断层扫描重建投影器。推导出了计算被投影到探测器像素上的体素体积与其对探测器记录的衰减线积分贡献之间关系的解析公式。基于这些公式,我们构建了一个接近精确的投影器和反投影器,特别适合代数重建技术和具有非均匀体素网格的分层重建方法。与传统假设均匀网格且固定体素大小的投影器不同,我们的方法能够实现体素的局部细化,从而支持自适应网格分辨率,并在感兴趣区域提高重建质量。我们实现了这种割线体投影器以及一个放松版本(速度优化版),并与两种已建立的投影器进行了比较:基于Siddon算法的射线追踪投影器和TT足迹投影器。结果表明,割线体投影器在大锥束角情况下比TT投影器具有更高的准确性。此外,割线体投影器的放松版本在保持相当精度的同时提供了显著的速度优势。相比之下,经过调整以达到相同精度的Siddon算法明显慢于割线体投影器。所有算法均在GPU优化的开源代数重建框架中实现。项目GitHub仓库为https://github.com/kulvait/KCT_cbct。

We introduce a novel class of projectors for 3D cone beam tomographic reconstruction. Analytical formulas are derived to compute the relationship between the volume of a voxel projected onto a detector pixel and its contribution to the line integral of attenuation recorded by that pixel. Based on these formulas, we construct a near-exact projector and backprojector, particularly suited for algebraic reconstruction techniques and hierarchical reconstruction approaches with nonuniform voxel grids. Unlike traditional projectors, which assume a uniform grid with fixed voxel sizes, our method enables local refinement of voxels, allowing for adaptive grid resolution and improved reconstruction quality in regions of interest. We have implemented this cutting voxel projector along with a relaxed, speed-optimized version and compared them to two established projectors: a ray-tracing projector based on Siddon's algorithm and a TT footprint projector. Our results demonstrate that the cutting voxel projector achieves higher accuracy than the TT projector, especially for large cone beam angles. Furthermore, the relaxed version of the cutting voxel projector offers a significant speed advantage, while maintaining comparable accuracy. In contrast, Siddon's algorithm, tuned to achieve the same accuracy, is considerably slower than the cutting voxel projector. All algorithms are implemented in a GPU optimized open-source framework for algebraic reconstruction. GitHub repository of the project https://github.com/kulvait/KCT_cbct.