Fourier-based reconstruction for fully 3-D PET: Optimization of interpolation parameters

被引:16
作者
Matej, Samuel [1 ]
Kazantsev, Ivan G. [1 ]
机构
[1] Univ Penn, Dept Radiol, Philadelphia, PA 19104 USA
关键词
direct Fourier method (DFM); fully three-dimensional (3-D) image reconstruction; gridding interpolation; positron emission tomography (PET); reprojection; spectrum density weighting;
D O I
10.1109/TMI.2006.873219
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Fourier-based approaches for three-dimensional (3-D) reconstruction are based on the relationship between the 3-D Fourier transform (FT) of the volume and the two-dimensional (2-D) FT of a parallel-ray projection of the volume. The critical step in the Fourier-based methods is the estimation of the samples of the 3-D transform of the image from the samples of the 2-D transforms of the projections on the planes through the origin of Fourier space, and vice versa for forward-projection (reprojection). The Fourier-based approaches have the potential for very fast reconstruction, but their straightforward implementation might lead to unsatisfactory results if careful attention is not paid to interpolation and weighting functions. In our previous work, we have investigated optimal interpolation parameters for the Fourier-based forward and back-projectors for iterative image reconstruction. The optimized interpolation kernels were shown to provide excellent quality comparable to the ideal sine interpolator. This work presents an optimization of interpolation parameters of the 3-D direct Fourier method with Fourier reprojection (3D-FRP) for fully 3-D positron emission tomography (PET) data with incomplete oblique projections. The reprojection step is needed for the estimation (from an initial image) of the missing portions of the oblique data. In the 3D-FRP implementation, we use the gridding interpolation strategy, combined with proper weighting approaches in the transform and image domains. We have found that while the 3-D reprojection step requires similar optimal interpolation parameters as found in our previous studies on Fourier-based iterative approaches, the optimal interpolation parameters for the main 3D-FRP reconstruction stage are quite different. Our experimental results confirm that for the optimal interpolation parameters a very good image accuracy can be achieved even without any extra spectral oversampling, which is a common practice to decrease errors caused by interpolation in Fourier reconstruction.
引用
收藏
页码:845 / 854
页数:10
相关论文
共 44 条
[1]  
BAL G, 2002, P 2002 IEEE NUCL SCI
[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]   Rapid gridding reconstruction with a minimal oversampling ratio [J].
Beatty, PJ ;
Nishimura, DG ;
Pauly, JM .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 2005, 24 (06) :799-808
[4]  
Bendriem B., 1998, THEORY PRACTICE 3D P, DOI [10.1007/978-94-017-3475-2, DOI 10.1007/978-94-017-3475-2]
[5]  
Braddick O, 1998, PERCEPTION, V27, P1
[6]   Fast fully 3-D image reconstruction in PET using planograms [J].
Brasse, D ;
Kinahan, PE ;
Clackdoyle, R ;
Defrise, M ;
Comtat, C ;
Townsend, DW .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 2004, 23 (04) :413-425
[7]   FULLY 3-DIMENSIONAL POSITRON EMISSION TOMOGRAPHY [J].
COLSHER, JG .
PHYSICS IN MEDICINE AND BIOLOGY, 1980, 25 (01) :103-115
[8]  
DANIELSSON P, 1997, P 1997 IEEE NUCL SCI
[9]   Distance-driven projection and backprojection in three dimensions [J].
De Man, B ;
Basu, S .
PHYSICS IN MEDICINE AND BIOLOGY, 2004, 49 (11) :2463-2475
[10]  
Dunne S., 1990, Proceedings of the First Conference on Visualization in Biomedical Computing (Cat. No.90TH0311-1), P11, DOI 10.1109/VBC.1990.109296