A fast and accurate multilevel inversion of the radon transform

A number of imaging technologies reconstruct an image function from its Radon projection using the convolution backprojection method. The convolution is an O(N-2 log N) algorithm, where the image consists of N x N pixels, while the backprojection is an O(N-3) algorithm, thus constituting the major computational burden of the convolution backprojection method. An O(N-2 log N) multilevel backprojection method is presented here. When implemented with a Fourier-domain postprocessing technique, also presented here, the resulting image quality is similar or superior to the image quality of the classical backprojection technique.