Fast and exact 2D image reconstruction by means of Chebyshev decomposition and backprojection.

A new algorithm for the reconstruction of two-dimensional (2D) images from projections is described. The algorithm is based on the decomposition of the projections into Chebyshev polynomials of the second kind, which are the ideal basis functions for this application. The Chebyshev decomposition is done via the fast discrete sine transform. A discrete reconstruction filter is applied that corresponds to the ramp filter used in standard filtered backprojection (FBP) reconstruction. In contrast to FBP, the filter is applied to the Chebyshev coefficients and not to the Fourier coefficients of the projections. Then the reconstructed image is simply obtained by means of backprojection. Consequently, the method can be considered as a Chebyshev-domain filtered backprojection (CD-FBP). The total calculation time is dominated by the backprojection step only and is comparable to FBP. The merits of CD-FBP as compared with standard FBP are that: (a) The result is exact if the 2D function to be reconstructed can be decomposed into polynomials of finite degree, and if the sampling is adequate. Otherwise a polynomial approximation results. (b) The algorithm is inherently discrete. (c) It is particularly well suited for reconstructions from projections with non-equidistant samples that occur for instance in 2D PET (positron emission tomography) imaging and in a special form of fan beam scanning. Examples of applications comprise reconstructions of the Shepp and Logan head phantom in various sampling geometries, and a real PET test object. In the PET example an increased resolution is observed in comparison with standard FBP.

[1]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .

[2]  R M Lewitt,et al.  Performance of the Fourier rebinning algorithm for PET with large acceptance angles. , 1998, Physics in medicine and biology.

[3]  R. Marr,et al.  On the reconstruction of a function on a circular domain from a sampling of its line integrals , 1974 .

[4]  W. J. Lorenz,et al.  Performance evaluation of the whole-body PET scanner ECAT EXACT HR + , 1997 .

[5]  Steven G. Johnson,et al.  The Fastest Fourier Transform in the West , 1997 .

[6]  G. N. Ramachandran,et al.  Three-dimensional reconstruction from radiographs and electron micrographs: application of convolutions instead of Fourier transforms. , 1971, Proceedings of the National Academy of Sciences of the United States of America.

[7]  A. Cormack.  Representation of a Function by Its Line Integrals, with Some Radiological Applications , 1963 .

[8]  William H. Press,et al.  Numerical Recipes in C, 2nd Edition , 1992 .

[9]  Michel Defrise,et al.  Exact and approximate rebinning algorithms for 3-D PET data , 1997, IEEE Transactions on Medical Imaging.

[10]  S J Riederer,et al.  Ripple suppression during reconstruction in transverse tomography. , 1975, Physics in medicine and biology.

[11]  Ronald N. Bracewell,et al.  The Fourier Transform and Its Applications , 1966 .

[12]  Alfred K. Louis.  Orthogonal Function Series Expansions and the Null Space of the Radon Transform , 1984 .

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

[14]  G Brix,et al.  Performance evaluation of a whole-body PET scanner using the NEMA protocol. National Electrical Manufacturers Association. , 1997, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[15]  Thomas Bortfeld,et al.  The exponential radon transform and projection filtering in radiotherapy planning , 1995, Int. J. Imaging Syst. Technol..

[16]  A K Louis,et al.  A novel filter design technique in 2D computerized tomography , 1996 .