O(N/sup 3/ 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(N/sup 3/ log/sub 2/ N) operations for reconstruction of an N/spl times/N/spl times/N volume from O(N/sup 2/) 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(N/sup 5/) computations under the same circumstances.

[1]  Kevin Barraclough,et al.  I and i , 2001, BMJ : British Medical Journal.

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

[3]  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.

[4]  P. Grangeat.  Mathematical framework of cone beam 3D reconstruction via the first derivative of the radon transform , 1991 .

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

[6]  R. Marr,et al.  On Two Approaches to 3D Reconstruction in NMR Zeugmatography , 1981 .

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

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

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

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

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

[12]  A. W. M. van den Enden,et al.  Discrete Time Signal Processing , 1989 .

[13]  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.

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

[15]  Xiaochuan Pan.  Quasi-bandlimited properties of Radon transforms and their implications for increasing angular sampling densities , 1998, IEEE Transactions on Medical Imaging.