Abstract:The Radon transform and its adjoint, the back-projection operator, can both be expressed as convolutions in log-polar coordinates. Hence, fast algorithms for the application of these operators can be constructed by using the FFT, if data is resampled at log-polar coordinates. Radon data is typically measured on an equally spaced grid in polar coordinates, and reconstructions are represented (as images) in Cartesian coordinates. Therefore, in addition to FFT, several steps of interpolation have to be conducted in order to apply the Radon transform and the back-projection operator by means of convolutions. However, in comparison to the interpolation conducted in Fourier-based gridding methods, the interpolation performed in the Radon and image domains will typically deal with functions that are substantially less oscillatory. Reasonable reconstruction results can thus be expected using interpolation schemes of moderate order. The approach also provides better control over the artifacts that can appear due to measurement errors. Both the interpolation and the FFT operations can be efficiently implemented on graphical processor units (GPUs). For the interpolation, it is possible to make use of the fact that linear interpolation is hard-wired on GPUs, meaning that it has the same computational cost as direct memory access. Cubic order interpolation schemes can be constructed by combining linear interpolation steps, and this provides important computation speedup. We provide details about how the Radon transform and the back-projection can be implemented efficiently as convolution operators on GPUs. For large data sizes, these algorithms are several times faster than those of other software packages based on GPU implementations of the Radon transform and the back-projection operator. Moreover, the gain in computational speed is substantially higher when comparing against other CPU-based algorithms.

Fast Computation of the Multidimensional Discrete Fourier Transform and Discrete Backward Fourier Transform on Sparse Grids

Fast Discrete Algorithms for Sparse Fourier Expansions of High Dimensional Functions

A Pre-splitting Green's Function Based Hybrid Fast Algorithm for Multiscale Problems

Minimizing communication in the multidimensional FFT

Orthogonal Polynomial Expansions on Sparse Grids.

Fast Fourier transforms for the evaluation of convolution products: CPU versus GPU implementation

Fast Kernel Summation in High Dimensions via Slicing and Fourier Transforms

Dimension-independent Sparse Fourier Transform

Flexible Multi-Dimensional FFTs for Plane Wave Density Functional Theory Codes

A Fast Algorithm for Computing Multidimensional DCT on Certain Small Sizes

Fast Algorithms and Efficient GPU Implementations for the Radon Transform and the Back-Projection Operator Represented as Convolution Operators

Fast computation of general Fourier Transforms on GPUS

Hermite Spectral Method with Hyperbolic Cross Approximations to High-Dimensional Parabolic Pdes

An adaptive fast Gauss transform in two dimensions

A Dual-space Multilevel Kernel-splitting Framework for Discrete and Continuous Convolution

Improved Reconstruction for Fourier-Sparse Signals

Gohberg-Semencul Factorization-Based Fast Implementation of Sparse Bayesian Learning With a Fourier Dictionary

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High Frequency Wave Equations

Matrix Factorization for Fast DCT Algorithms

Super-resolution and Robust Sparse Continuous Fourier Transform in Any Constant Dimension: Nearly Linear Time and Sample Complexity

Fast and Scalable FFT-Based GPU-Accelerated Algorithms for Hessian Actions Arising in Linear Inverse Problems Governed by Autonomous Dynamical Systems