Quasi Band-Limited Properties of Radon Transforms and Their Implications for Increasing Angular Sampling Densities

The n-dimensional (n-D) radon transform, which forms the mathematical basis for a broad variety of tomographic imaging applications, can be viewed as an n-D function in n-D sinogram space. Accurate reconstruction of continuous or discrete tomographic images requires full knowledge of the radon transform in the corresponding n-D sinogram space. In practice, however, one can have only a finite set of discrete samples of the radon transform in the sinogram space. One often derives the desired full knowledge of the radon transform from its discrete samples by invoking various interpolation algorithms. According to the Wittaker-Shannon sampling theorem, a necessary condition for a full and unique recovery of the radon transform from its discrete samples is that the radon transform itself be bandlimited. Therefore, it is necessary to analyze the bandlimited properties of the radon transform. In this work, we analyze explicitly the bandlimited properties of the radon transform and show that the radon transform is mathematically quasi-bandlimited [or essentially bandlimited] in two quantitative senses and can essentially be treated as bandlimited in practice. The quasi-bandlimited properties can be used for increasing the angular sampling density of the radon transform.

[1]  G. A. Watson A treatise on the theory of Bessel functions , 1944 .

[2]  Stephen J. Norton,et al.  Reconstruction of a reflectivity field from line integrals over circular paths , 1980 .

[3]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

[4]  R M Mersereau,et al.  Direct Fourier transform techniques in 3-D image reconstruction. , 1976, Computers in biology and medicine.

[5]  Robert M. Lewitt,et al.  Fourier correction for spatially variant collimator blurring in SPECT , 1995, IEEE Trans. Medical Imaging.

[6]  S. Deans The Radon Transform and Some of Its Applications , 1983 .

[7]  R. Bracewell The Fourier Transform and Its Applications , 1966 .

[8]  A. J. Jerri The Shannon sampling theorem—Its various extensions and applications: A tutorial review , 1977, Proceedings of the IEEE.

[9]  H Guan,et al.  Computed tomography using algebraic reconstruction techniques (ARTs) with different projection access schemes: a comparison study under practical situations , 1996, Physics in medicine and biology.

[10]  A. Erdélyi,et al.  Higher Transcendental Functions , 1954 .

[11]  R. Lewitt Reconstruction algorithms: Transform methods , 1983, Proceedings of the IEEE.

[12]  A. Kak,et al.  A computational study of reconstruction algorithms for diffraction tomography: Interpolation versus filtered-backpropagation , 1983 .

[13]  A. Lindgren,et al.  Sampling the 2-D Radon transform , 1981 .

[14]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .

[15]  R. Hingorani,et al.  Direct Fourier reconstruction in computer tomography , 1981 .

[16]  Gabor T. Herman,et al.  Image reconstruction from projections : the fundamentals of computerized tomography , 1980 .

[17]  M. Kaveh,et al.  Reconstructive tomography and applications to ultrasonics , 1979, Proceedings of the IEEE.

[18]  Gabor T. Herman Image Reconstruction from Projections , 1979 .

[19]  Harrison H. Barrett,et al.  III The Radon Transform and Its Applications , 1984 .

[20]  A. Papoulis Systems and transforms with applications in optics , 1981 .

[21]  N C Yang,et al.  The circular harmonic transform for SPECT reconstruction and boundary conditions on the Fourier transform of the sinogram. , 1988, IEEE transactions on medical imaging.