Rectification for cone-beam projection and backprojection

The purpose of this paper is to derive a technique for accelerating the computation of cone-beam forward and backward projections that are the basic steps of tomographic reconstruction. The cone-beam geometry of C-arm systems is commonly described with projection matrices. Such matrices provide a continuous framework for analyzing the flow of operations needed to compute backprojection for analytical reconstruction, as well as the combination of forward and backward projections for iterative reconstruction. The proposed rectification technique resamples the original data to planes that are aligned with two of the reconstructed volume main axes, so that the original cone-beam geometry can be replaced by a simpler geometry, where succession of plane magnifications are involved only. Rectification generalizes previous independent results to the cone-beam backprojection of preprocessed data as well as to cone-beam iterative reconstruction. The memory access pattern of simple magnifications provides superior predictability and is, therefore, easier to optimize, independently of the choice of the interpolation technique. Rectification is also shown to provide control over interpolation errors through oversampling, allowing tradeoffs between computation speed and precision to be set. Experimental results are provided for linear and nearest neighbor interpolations, based on simulations, as well as phantom and patient data acquired on a digital C-arm system