Fourier-based forward and back-projectors in iterative fan-beam tomographic image reconstruction

Fourier-based forward and back-projection methods have the potential to reduce computation demands in iterative tomographic image reconstruction. Interpolation errors are a limitation of conventional Fourier-based projectors. Recently, the min-max optimized Kaiser-Bessel interpolation within the nonuniform fast Fourier transform (NUFFT) approach has been applied in parallel-beam image reconstruction, whose results show lower approximation errors than conventional interpolation methods. However, the extension of min-max NUFFT approach to fan-beam data has not been investigated. We have extended the min-max NUFFT framework to the fan-beam tomography case, using the relationship between the fan-beam projections and corresponding projections in parallel-beam geometry. Our studies show that the fan-beam Fourier-based forward and back-projection methods can significantly reduce the computation time while still providing comparable accuracy as their space-based counterparts.

[1]  Yoram Bresler,et al.  NATIVE FAN-BEAM TOMOGRAPHIC RECONSTRUCTION , .

[2]  Augusto Marques Ferreira da Silva,et al.  Efficient NUFFT-based direct Fourier algorithm for fan beam CT reconstruction , 2004, Medical Imaging: Image Processing.

[3]  Jeffrey A. Fessler,et al.  Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction , 1997, IEEE Transactions on Medical Imaging.

[4]  Jeffrey A. Fessler,et al.  Nonuniform fast Fourier transforms using min-max interpolation , 2003, IEEE Trans. Signal Process..

[5]  Jeffrey A. Fessler,et al.  Conjugate-gradient preconditioning methods for shift-variant PET image reconstruction , 1999, IEEE Trans. Image Process..

[6]  P L Bellon,et al.  Fast direct Fourier methods, based on one- and two-pass coordinate transformations, yield accurate reconstructions of x-ray CT clinical images. , 1997, Physics in medicine and biology.

[7]  Shu Xiao,et al.  O(N/sup 2/ log N) native fan-beam tomographic reconstruction , 2002, Proceedings IEEE International Symposium on Biomedical Imaging.

[8]  B. De Man,et al.  Distance-driven projection and backprojection in three dimensions. , 2004, Physics in medicine and biology.

[9]  Donald Fraser,et al.  Interpolation by the FFT revisited-an experimental investigation , 1989, IEEE Trans. Acoust. Speech Signal Process..

[10]  Jeffrey A. Fessler,et al.  Iterative tomographic image reconstruction using Fourier-based forward and back-projectors , 2004, IEEE Transactions on Medical Imaging.

[11]  Henry Stark,et al.  Direct Fourier Reconstruction in Fan-Beam Tomography , 1987, IEEE Transactions on Medical Imaging.

[12]  Jeffrey A. Fessler,et al.  Spatial resolution properties of penalized-likelihood image reconstruction: space-invariant tomographs , 1996, IEEE Trans. Image Process..

[13]  Jeffrey A. Fessler,et al.  A min-max approach to the multidimensional nonuniform FFT: application to tomographic image reconstruction , 2001, Proceedings 2001 International Conference on Image Processing (Cat. No.01CH37205).

[14]  J. Fessler,et al.  Fourier-based forward and back-projectors in iterative fan-beam tomographic image reconstruction , 2006, IEEE Transactions on Medical Imaging.

[15]  Guy M. Besson,et al.  CT fan-beam parametrizations leading to shift-invariant filtering , 1996 .

[16]  A. Macovski,et al.  Selection of a convolution function for Fourier inversion using gridding [computerised tomography application]. , 1991, IEEE transactions on medical imaging.

[17]  J. D. O'Sullivan,et al.  A Fast Sinc Function Gridding Algorithm for Fourier Inversion in Computer Tomography , 1985, IEEE Transactions on Medical Imaging.

[18]  Yoram Bresler,et al.  Error analysis and performance optimization of fast hierarchical backprojection algorithms , 2001, IEEE Trans. Image Process..

[19]  James J. Clark,et al.  A transformation method for the reconstruction of functions from nonuniformly spaced samples , 1985, IEEE Trans. Acoust. Speech Signal Process..

[20]  Wayne Lawton,et al.  A new polar Fourier transform for computer-aided tomography and spotlight synthetic aperture radar , 1988, IEEE Trans. Acoust. Speech Signal Process..

[21]  Ken D. Sauer,et al.  A local update strategy for iterative reconstruction from projections , 1993, IEEE Trans. Signal Process..

[22]  C. Mistretta,et al.  A novel extension of the parallel-beam projection-slice theorem to divergent fan-beam and cone-beam projections. , 2005, Medical physics.

[23]  A. Averbuch,et al.  3D Fourier based discrete Radon transform , 2003 .

[24]  Johan Waldén,et al.  Analysis of the direct Fourier method for computer tomography , 2000, IEEE Transactions on Medical Imaging.

[25]  J. Fessler.  Iterative Tomographic Image Reconstruction Using Nonuniform Fast Fourier Transforms , 2002 .

[26]  Karsten Fourmont.  Non-Equispaced Fast Fourier Transforms with Applications to Tomography , 2003 .

[27]  Gabriele Steidl,et al.  A new linogram algorithm for computerized tomography , 2001 .

[28]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[29]  P M Joseph,et al.  The exponential edge-gradient effect in x-ray computed tomography. , 1981, Physics in medicine and biology.

[30]  Philippe Moireau,et al.  Fast numerical inversion of the attenuated Radon transform with full and partial measurements , 2004 .