Fourier-based reprojection methods have the potential to reduce the computation time in iterative tomographic image reconstruction. Interpolation errors are a limitation of Fourier-based reprojection methods. We apply a min-max interpolation method for the nonuniform fast Fourier transform (NUFFT) to minimize the interpolation errors. Numerical results show that the min-max NUFFT approach provides substantially lower approximation errors in tomographic reprojection and backprojection than conventional interpolation methods.

We investigate the use of direct-Fourier (DF) image reconstruction in computerized tomography and synthetic aperture radar (SAR). One of our aims is to determine why the convolutionbackprojection (CBP) method is favored over DF methods in tomography, while DF methods are virtually always used in SAR. We show that the CBP algorithm is equivalent to DF reconstruction using a Jacobian-weighted 2-D periodic sinc-kernel interpolator. This interpolation is not optimal in any sense, which suggests that DF algorithms utilizing optimal interpolators may surpass CBP in image quality. We consider use of two types of DF interpolation: a windowed sinc kernel, and the least-squares optimal Yen interpolator. Simulations show that reconstructions using the Yen interpolator do not possess the expected visual quality, because of regularization needed to preserve numerical stability. Next, we show that with a concentric-squares sampling scheme, DF interpolation can be performed accurately and efficiently...

This paper also proposed another reconstruction method based on a direct approximation of the Fourier inversion formula using a twodimensional (2-D) trapezoidal rule. In addition, the possibility of reconstruction from a concentric-squares raster was discussed. Numerous simple interpolators have been tried in DF reconstruction with the results compared with CBP [33]. In [34] and [35], the concept of angular bandlimiting was used to interpolate the polar data onto a Cartesian grid. In [36], a DF reconstruction using bilinear interpolation for diffraction tomography provided image quality that was comparable to that produced by the CBP algorithm. Very good reconstruction quality was obtained in [37] and [38] using a spline interpolator, or a hybrid type of spline interpolator. The notion of "gridding" was introduced in [39] as a method of obtaining optimal inversion of Fourier data. An optimal gridding function was proposed, and successful results were obtained when applied to the tomographic reconstruction problem. In [40], several different gridding functions were tried for DF reconstruction, and the performances were compared. In [41, 42], the linogram reconstruction method was proposed as a form of DF reconstruction. The data collection grid in the linogram method is the same as in the concentric-squares sampling scheme. The inversion of the Fourier data in [41, 42] was accomplished by first applying the chirp-z transform in one direction and then computing FFTs in the other direction. In CT, many of these attempts at DF reconstruction have given a poorer result than the CBP algorithm, due to the error incurred in the process of the polar-to-Cartesian interpolation. The attraction of DF reconstruction, however, is that it is thought to require less computation than ...

All errors in this manual and all bugs in the program

1.1 Statement of purpose and history................................. 5