Mathematical Equivalence of Zero-Padding Interpolation and Circular Sampling Theorem Interpolation with Implications for Direct Fourier Image Reconstruction.

PJ La Riviere,X Pan

DOI: https://doi.org/10.1117/12.310838

1998-01-01

Abstract:The speed and accuracy of Direct Fourier image reconstruction methods have long been hampered by the need to interpolate between the polar grid of Fourier data that is obtained from the measured projection data and the Cartesian grid of Fourier data that is needed to recover an image using the 2D FFT. Fast but crude interpolation schemes such as bilinear interpolation often lead to unacceptable image artifacts, while more sophisticated but computationally intense techniques such as circular sampling theorem (CST) interpolation negate the speed advantages afforded by the use of the 2D FFT. One technique that has been found to yield high-quality images without much computational penalty is a hybrid one in which zero-padding interpolation is first used to increase the density of samples on the polar grid after which bilinear interpolation onto the Cartesian grid is performed. In this work, we attempt to account for the success of the approach relative to the CST approach in three ways. First and most importantly, we establish that zero-padding interpolation of periodic functions that are sampled in accordance with the Nyquist criterion-precisely the sort of function encountered in the angular dimension of the polar grid-is exact and equivalent to circular sampling theorem interpolation. Second, we point out that both approaches make comparable approximations in interpolating in the radial direction. Finally, we indicate that the error introduced by the bilinear interpolation step in the zero-padding approach can be minimized by choosing sufficiently large zero-padding factors.