Iterative tomographic image reconstruction using Fourier-based forward and back-projectors

被引:86
作者
Matej, S [1 ]
Fessler, JA
Kazantsev, IG
机构
[1] Univ Penn, Dept Radiol, Philadelphia, PA 19104 USA
[2] Univ Michigan, Dept Elect Engn & Comp Sci, Ann Arbor, MI 48109 USA
关键词
fast forward and back-projectors; gridding; iterative tomographic reconstruction; min-max interpolation; nonuniform FFT;
D O I
10.1109/TMI.2004.824233
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Iterative image reconstruction algorithms play an increasingly important role in modern tomographic systems, especially in emission tomography. With the fast increase of the sizes of the tomographic data, reduction of the computation demands of the reconstruction algorithms is of great importance. Fourier-based forward and back-projection methods have the potential to considerably reduce the computation time in iterative reconstruction. Additional substantial speed-up of those approaches can be obtained utilizing powerful and cheap off-the-shelf fast Fourier transform (FFT) processing hardware. The Fourier reconstruction approaches are based on the relationship between the Fourier transform of the image and Fourier transformation of the parallel-ray projections. The critical two steps are the estimations of the samples of the projection transform, on the central section through the origin of Fourier space, from the samples of the transform of the image, and vice versa for back-projection. Interpolation errors are a limitation of Fourier-based reconstruction methods. We have applied min-max optimized Kaiser-Bessel interpolation within the nonuniform FFT (NUFFT) framework and devised ways of incorporation of resolution models into the Fourier-based iterative approaches. Numerical and computer simulation results show that the min-max NUFFT approach provides substantially lower approximation errors in tomographic forward and back-projection than conventional interpolation methods. Our studies have further confirmed that Fourier-based projectors using the NUFFT approach provide accurate approximations to their space-based counterparts but with about ten times faster computation, and that they are viable candidates for fast iterative image reconstruction.
引用
收藏
页码:401 / 412
页数:12
相关论文
共 40 条
[1]  
AVERBUCH A, 2004, IN PRESS SIAM SCI CO
[2]   O(N2 log2 N) filtered backprojection reconstruction algorithm for tomography [J].
Basu, S ;
Bresler, Y .
IEEE TRANSACTIONS ON IMAGE PROCESSING, 2000, 9 (10) :1760-1773
[3]  
BRASSE D, 2004, IEEE T MED IMAG APR, P413
[4]   A fast and accurate Fourier algorithm for iterative parallel-beam tomography [J].
Delaney, AH ;
Bresler, Y .
IEEE TRANSACTIONS ON IMAGE PROCESSING, 1996, 5 (05) :740-753
[5]  
DEMAN B, P 2002 IEEE NUCL SCI
[6]  
Dunne S., 1990, Proceedings of the First Conference on Visualization in Biomedical Computing (Cat. No.90TH0311-1), P11, DOI 10.1109/VBC.1990.109296
[7]   Ordered subsets algorithms for transmission tomography [J].
Erdogan, H ;
Fessler, JA .
PHYSICS IN MEDICINE AND BIOLOGY, 1999, 44 (11) :2835-2851
[8]   PENALIZED WEIGHTED LEAST-SQUARES IMAGE-RECONSTRUCTION FOR POSITRON EMISSION TOMOGRAPHY [J].
FESSLER, JA .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 1994, 13 (02) :290-300
[9]   Nonuniform fast Fourier transforms using min-max interpolation [J].
Fessler, JA ;
Sutton, BP .
IEEE TRANSACTIONS ON SIGNAL PROCESSING, 2003, 51 (02) :560-574
[10]  
Fessler JA, 2001, IEEE IMAGE PROC, P706, DOI 10.1109/ICIP.2001.959143