O(N3 log N) backprojection algorithm for the 3-D radon transform.

We present a novel backprojection algorithm for three-dimensional (3-D) radon transform data that requires O(N3 log2 N) operations for reconstruction of an N x N x N volume from O(N2) plane-integral projections. Our algorithm uses a hierarchical decomposition of the 3-D radon transform to recursively decompose the backprojection operation. Simulations are presented demonstrating reconstruction quality comparable to the standard filtered backprojection, which requires O(N5) computations under the same circumstances.

[1]  Xiaochuan Pan,et al.  Quasi Band-Limited Properties of Radon Transforms and Their Implications for Increasing Angular Sampling Densities , 1998, IEEE Trans. Medical Imaging.

[2]  Peter Steffen,et al.  An efficient Fourier method for 3-D radon inversion in exact cone-beam CT reconstruction , 1998, IEEE Transactions on Medical Imaging.

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

[4]  R. Van de Walle,et al.  The motion sensitivity of magnetic resonance imaging: a comparison between Fourier transform imaging and projection reconstruction , 1995, Proceedings of 17th International Conference of the Engineering in Medicine and Biology Society.

[5]  Charles C. Peck,et al.  Implementations, comparisons, and an investigation of heuristic techniques for cone-beam tomography , 1996, IEEE Trans. Medical Imaging.

[6]  V. Palamodov,et al.  Localization of harmonic decomposition of the Radon transform , 1995 .

[7]  Yoram Bresler,et al.  A multilevel domain decomposition algorithm for fast O(N2logN) reprojection of tomographic images , 2000, IEEE Trans. Image Process..

[8]  C Axelsson,et al.  Three-dimensional reconstruction from cone-beam data in O(N3logN) time , 1994, Physics in medicine and biology.

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

[10]  Yoram Bresler,et al.  O(N2log2N) filtered backprojection reconstruction algorithm for tomography , 2000, IEEE Trans. Image Process..

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

[12]  Alan V. Oppenheim,et al.  Discrete-Time Signal Pro-cessing , 1989 .