Geometrical and Monte Carlo projectors in 3D PET reconstruction.

PURPOSE In the present work, the authors compare geometrical and Monte Carlo projectors in detail. The geometrical projectors considered were the conventional geometrical Siddon ray-tracer (S-RT) and the orthogonal distance-based ray-tracer (OD-RT), based on computing the orthogonal distance from the center of image voxel to the line-of-response. A comparison of these geometrical projectors was performed using different point spread function (PSF) models. The Monte Carlo-based method under consideration involves an extensive model of the system response matrix based on Monte Carlo simulations and is computed off-line and stored on disk. METHODS Comparisons were performed using simulated and experimental data of the commercial small animal PET scanner rPET. RESULTS The results demonstrate that the orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions yield better images in terms of contrast and spatial resolution than those obtained after using the conventional method and the multiray-based S-RT. Furthermore, the Monte Carlo-based method yields slight improvements in terms of contrast and spatial resolution with respect to these geometrical projectors. CONCLUSIONS The orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions represent satisfactory alternatives to factorizing the system matrix or to the conventional on-the-fly ray-tracing methods for list-mode reconstruction, where an extensive modeling based on Monte Carlo simulations is unfeasible.

[1]  M J Welch,et al.  Small animal imaging: current technology and perspectives for oncological imaging , 2002 .

[2]  E. Hoffman,et al.  Calculation of positron range and its effect on the fundamental limit of positron emission tomography system spatial resolution. , 1999, Physics in medicine and biology.

[3]  C. Byrne,et al.  Reconstruction of 2D PET data with Monte Carlo generated system matrix for generalized natural pixels , 2006, Physics in Medicine and Biology.

[4]  Steven G. Ross,et al.  Application and Evaluation of a Measured Spatially Variant System Model for PET Image Reconstruction , 2010, IEEE Transactions on Medical Imaging.

[5]  Long Zhang,et al.  Fast and memory-efficient Monte Carlo-based image reconstruction for whole-body PET. , 2010, Medical physics.

[6]  F. Fazio,et al.  Implementation and evaluation of a 3D one-step late reconstruction algorithm for 3D positron emission tomography brain studies using median root prior , 2002, European Journal of Nuclear Medicine and Molecular Imaging.

[7]  Craig S. Levin,et al.  Fast, Accurate and Shift-Varying Line Projections for Iterative Reconstruction Using the GPU , 2009, IEEE Transactions on Medical Imaging.

[8]  J. Ollinger Model-based scatter correction for fully 3D PET. , 1996, Physics in medicine and biology.

[9]  Vladimir Y. Panin,et al.  Fully 3-D PET reconstruction with system matrix derived from point source measurements , 2006, IEEE Transactions on Medical Imaging.

[10]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[11]  George Kontaxakis,et al.  3D-OSEM iterative image reconstruction for high-resolution PET using precalculated system matrix , 2006 .

[12]  Markus Alber,et al.  A parallelizable compression scheme for Monte Carlo scatter system matrices in PET image reconstruction , 2007, Physics in medicine and biology.

[13]  R Myers,et al.  The biological application of small animal PET imaging. , 2001, Nuclear medicine and biology.

[14]  Jose Luis Rubio Guivernau,et al.  Efficient methodologies for system matrix modelling in iterative image reconstruction for rotating high-resolution PET , 2010 .

[15]  D. Newport,et al.  A Single Scatter Simulation Technique for Scatter Correction in 3D PET , 1996 .

[16]  Ulla Ruotsalainen,et al.  Attenuation correction for PET using count-limited transmission images reconstructed with median root prior , 1999 .

[17]  D R Haynor,et al.  The use of importance sampling techniques to improve the efficiency of photon tracking in emission tomography simulations. , 1991, Medical physics.

[18]  Gunnar Brix,et al.  Fast implementation of the single scatter simulation algorithm and its use in iterative image reconstruction of PET data. , 2002, Physics in medicine and biology.

[19]  M. Rafecas,et al.  Use of a Monte Carlo-based probability matrix for 3-D iterative reconstruction of MADPET-II data , 2004, IEEE Transactions on Nuclear Science.

[20]  R M Lewitt,et al.  Multidimensional digital image representations using generalized Kaiser-Bessel window functions. , 1990, Journal of the Optical Society of America. A, Optics and image science.

[21]  D. McElroy,et al.  Use of a Monte-Carlo based probability matrix for 3D iterative reconstruction of MADPET-II data , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[22]  Andrew J. Reader,et al.  EM algorithm system modeling by image-space techniques for PET reconstruction , 2003 .

[23]  Jürgen Scheins,et al.  Analytical calculation of volumes-of-intersection for iterative, fully 3-D PET reconstruction , 2006, IEEE Transactions on Medical Imaging.

[24]  C.A. Bouman,et al.  Quantitative comparison of FBP, EM, and Bayesian reconstruction algorithms for the IndyPET scanner , 2003, IEEE Transactions on Medical Imaging.

[25]  R. Leahy,et al.  High-resolution 3D Bayesian image reconstruction using the microPET small-animal scanner. , 1998, Physics in medicine and biology.

[26]  Arman Rahmim,et al.  Statistical dynamic image reconstruction in state-of-the-art high-resolution PET , 2005, Physics in medicine and biology.

[27]  Habib Zaidi,et al.  Determination of the attenuation map in emission tomography. , 2003, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[28]  M. Desco,et al.  rPET detectors design and data processing , 2005, IEEE Nuclear Science Symposium Conference Record, 2005.

[29]  R. Fontaine,et al.  Fast, accurate and versatile Monte Carlo method for computing system matrix , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[30]  Qian Liao,et al.  An efficient anti-aliasing spectral continuous window shifting technique for PIV , 2005 .

[31]  R. Lewitt Alternatives to voxels for image representation in iterative reconstruction algorithms , 1992, Physics in medicine and biology.

[32]  Xin-Min Hua Monte Carlo simulation of Comptonization in inhomogeneous media , 1997 .

[33]  P.E. Kinahan,et al.  Modeling and incorporation of system response functions in 3D whole body PET , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

[34]  Thomas K. Lewellen,et al.  Modeling and incorporation of system response functions in 3-D whole body PET , 2006, IEEE Transactions on Medical Imaging.

[35]  S. Shokouhi,et al.  Statistical 3D image reconstruction for the RatCAP PET tomograph using a physically accurate, Monte Carlo based system matrix , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

[36]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[37]  Arion F. Chatziioannou,et al.  Molecular imaging of small animals with dedicated PET tomographs , 2001, European Journal of Nuclear Medicine and Molecular Imaging.

[38]  J L Rubio,et al.  Efficient methodologies for system matrix modelling in iterative image reconstruction for rotating high-resolution PET , 2010, Physics in medicine and biology.

[39]  Andrew J. Reader,et al.  Impact of Image-Space Resolution Modeling for Studies with the High-Resolution Research Tomograph , 2008, Journal of Nuclear Medicine.

[40]  Colas Schretter,et al.  A fast tube of response ray-tracer. , 2006, Medical physics.

[41]  C. Tsoumpas,et al.  Scatter simulation including double scatter , 2005, IEEE Nuclear Science Symposium Conference Record, 2005.

[42]  Andres Santos,et al.  Performance analysis of a low-cost small animal PET/SPECT scanner , 2007 .

[43]  A. Chatziioannou,et al.  Fully 3D system model estimation of OPET by Monte Carlo simulation , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

[44]  C. Knoess,et al.  PET-based molecular imaging in neuroscience , 2003, European Journal of Nuclear Medicine and Molecular Imaging.