Accurate and efficient modeling of the detector response in small animal multi-head PET systems

In fully three-dimensional PET imaging, iterative image reconstruction techniques usually outperform analytical algorithms in terms of image quality provided that an appropriate system model is used. In this study we concentrate on the calculation of an accurate system model for the YAP-(S)PET II small animal scanner, with the aim to obtain fully resolution- and contrast-recovered images at low levels of image roughness. For this purpose we calculate the system model by decomposing it into a product of five matrices: (1) a detector response component obtained via Monte Carlo simulations, (2) a geometric component which describes the scanner geometry and which is calculated via a multi-ray method, (3) a detector normalization component derived from the acquisition of a planar source, (4) a photon attenuation component calculated from x-ray computed tomography data, and finally, (5) a positron range component is formally included. This system model factorization allows the optimization of each component in terms of computation time, storage requirements and accuracy. The main contribution of this work is a new, efficient way to calculate the detector response component for rotating, planar detectors, that consists of a GEANT4 based simulation of a subset of lines of flight (LOFs) for a single detector head whereas the missing LOFs are obtained by using intrinsic detector symmetries. Additionally, we introduce and analyze a probability threshold for matrix elements of the detector component to optimize the trade-off between the matrix size in terms of non-zero elements and the resulting quality of the reconstructed images. In order to evaluate our proposed system model we reconstructed various images of objects, acquired according to the NEMA NU 4-2008 standard, and we compared them to the images reconstructed with two other system models: a model that does not include any detector response component and a model that approximates analytically the depth of interaction as detector response component. The comparisons confirm previous research results, showing that the usage of an accurate system model with a realistic detector response leads to reconstructed images with better resolution and contrast recovery at low levels of image roughness.

[1]  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.

[2]  M. Rafecas,et al.  Effect of noise in the probability matrix used for statistical reconstruction of PET data , 2002, IEEE Transactions on Nuclear Science.

[3]  R. Laforest,et al.  Positron range modeling for statistical PET image reconstruction , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[4]  Magdalena Rafecas,et al.  Comparison of basis functions for 3D PET reconstruction using a Monte Carlo system matrix , 2012, Physics in medicine and biology.

[5]  F. Beekman,et al.  Monte Carlo-based statistical SPECT reconstruction: influence of number of photon tracks , 2005, IEEE Transactions on Nuclear Science.

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

[7]  S R Cherry,et al.  An improved analytical detector response function model for multilayer small-diameter PET scanners. , 2003, Physics in medicine and biology.

[8]  R Taschereau,et al.  NEMA NU-4 performance evaluation of PETbox4, a high sensitivity dedicated PET preclinical tomograph , 2013, Physics in medicine and biology.

[9]  Richard Laforest,et al.  Quantitative accuracy of MAP reconstruction for dynamic PET imaging in small animals. , 2012, Medical physics.

[10]  Jinyi Qi,et al.  Iterative image reconstruction for positron emission tomography based on a detector response function estimated from point source measurements , 2009, Physics in medicine and biology.

[11]  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.

[12]  Magnus Dahlbom,et al.  PET system calibrations and corrections for quantitative and spatially accurate images , 1989 .

[13]  Simon R. Cherry,et al.  Comparison of 3-D maximum a posteriori and filtered backprojection algorithms for high-resolution animal imaging with microPET , 2000, IEEE Transactions on Medical Imaging.

[14]  Jian Zhou,et al.  Fast and efficient fully 3D PET image reconstruction using sparse system matrix factorization with GPU acceleration. , 2011, Physics in medicine and biology.

[15]  A. Jackson,et al.  Single scan parameterization of space-variant point spread functions in image space via a printed array: the impact for two PET/CT scanners , 2011, Physics in medicine and biology.

[16]  K Thielemans,et al.  Image-based point spread function implementation in a fully 3D OSEM reconstruction algorithm for PET , 2010, Physics in medicine and biology.

[17]  I. Buvat,et al.  A review of partial volume correction techniques for emission tomography and their applications in neurology, cardiology and oncology , 2012, Physics in medicine and biology.

[18]  Thomas Beyer,et al.  The use of X-ray CT for attenuation correction of PET data , 1994, Proceedings of 1994 IEEE Nuclear Science Symposium - NSS'94.

[19]  Xiaoli Li,et al.  Performance Evaluation of Small Animal PET Scanners With Different System Designs , 2013, IEEE Transactions on Nuclear Science.

[20]  A. Martineau,et al.  A method for accurate modelling of the crystal response function at a crystal sub-level applied to PET reconstruction , 2011, Physics in medicine and biology.

[21]  R.L. Harrison,et al.  Measured spatially variant system response for PET image reconstruction , 2005, IEEE Nuclear Science Symposium Conference Record, 2005.

[22]  Dan J Kadrmas,et al.  LOR-OSEM: statistical PET reconstruction from raw line-of-response histograms , 2004, Physics in medicine and biology.

[23]  Bjorn De Sutter,et al.  A Fast Algorithm to Calculate the Exact Radiological Path through a Pixel or Voxel Space , 1998 .

[24]  Roger Lecomte,et al.  Detector response models for statistical iterative image reconstruction in high resolution PET , 1998 .

[25]  Pablo Aguiar,et al.  Geometrical and Monte Carlo projectors in 3D PET reconstruction. , 2010, Medical physics.

[26]  M Defrise,et al.  Non-Gaussian space-variant resolution modelling for list-mode reconstruction , 2010, Physics in medicine and biology.

[27]  A. Rahmim,et al.  Monte Carlo-based evaluation of inter-crystal scatter and penetration in the PET subsystem of three GE Discovery PET/CT scanners , 2011 .

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

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

[30]  Anthonin Reilhac,et al.  Analytical positron range modelling in heterogeneous media for PET Monte Carlo simulation , 2011, Physics in medicine and biology.

[31]  Paul Kinahan,et al.  Noise and signal properties in PSF-based fully 3D PET image reconstruction: an experimental evaluation , 2010, Physics in medicine and biology.

[32]  A. Rahmim,et al.  Space-variant and anisotropic resolution modeling in list-mode EM reconstruction , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[33]  Giovanni Corsini,et al.  Evaluation of the performance of the YAP-(S)PET scanner and its application in neuroscience , 2007 .

[34]  A. Del Guerra,et al.  Performance evaluation of the fully engineered YAP-(S)PET scanner for small animal imaging , 2006, IEEE Symposium Conference Record Nuclear Science 2004..

[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]  Andrew J. Reader,et al.  Impact of Image-Space Resolution Modeling for Studies with the High-Resolution Research Tomograph , 2008, Journal of Nuclear Medicine.

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

[38]  Ariela Sofer,et al.  Evaluation of 3D reconstruction algorithms for a small animal PET camera , 1996 .

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

[40]  R. Leahy,et al.  Accurate geometric and physical response modelling for statistical image reconstruction in high resolution PET , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[41]  I Buvat,et al.  Fully 3D Monte Carlo reconstruction in SPECT: a feasibility study , 2005, Physics in medicine and biology.

[42]  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.

[43]  R. Siddon.  Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[44]  A Geissbuhler,et al.  A normalization technique for 3D PET data. , 1991, Physics in medicine and biology.

[45]  Magdalena Rafecas,et al.  Comparison of different approaches based on Monte Carlo methods to calculate the system matrix for small animal PET , 2006 .

[46]  Gianluigi Zanetti,et al.  Multi-ray-based system matrix generation for 3D PET reconstruction , 2008, Physics in medicine and biology.

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

[48]  Jonathan A. Disselhorst,et al.  Characterization and optimization of image quality as a function of reconstruction algorithms and parameter settings in a Siemens Inveon small-animal PET scanner using the NEMA NU 4-2008 standards , 2011 .

[49]  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.