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

TL;DR

This paper presents a new adaptive voxel projector for 3D cone beam CT reconstruction that improves accuracy and speed, especially for large angles, by enabling local voxel refinement and is implemented in a GPU framework.

Contribution

The paper introduces a novel cutting voxel projector with analytical formulas, allowing adaptive grid refinement and superior accuracy and speed compared to traditional methods.

Findings

Higher accuracy than TT projector for large cone angles

Relaxed version offers significant speed improvements

Outperforms Siddon's algorithm in accuracy and speed

Abstract

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.