Reconstruction of 2D PET data with Monte Carlo generated natural pixels

In discrete detector PET, natural pixels are image basis functions calculated from responses of detector pairs. By using reconstruction with natural pixels the discretization of the object into a predefined grid can be avoided. Instead of trying to find the object x from system matrix A and measured data p by solving Ax = p, one solves Mq = AA/sup T/q = p. The backprojection of the final q results in the estimated image of the object x. In previous work natural pixels were calculated by intersections of strip functions or combinations of strip functions. However, in PET systems the detector response is not correctly described by a strip function. In PET the detector response is better described by a Gaussian tube. It is quite difficult to calculate the intersections of all possible lines of response (LOR) in the scanner. This paper proposes an easy and efficient way to generate the matrix M directly by Monte Carlo simulation. Using the natural pixel matrix M has several advantages over using the conventional system matrix A. The discretisation of the object into pixels, which requires an angular dependent projector and backprojector, is avoided. Due to rotational symmetry in the PET scanner the matrix M is block circulant and only the first block row needs to be stored. The block circulant property of the matrix allows use of fast direct methods to calculate the solution. Iterative methods to solve the system are preferred because of numerical stability.

[1]  Grant T. Gullberg,et al.  An SVD reconstruction algorithm using a natural pixel representation of the attenuated Radon transform , 1996 .

[2]  R.M. Leahy,et al.  Fully 3D Bayesian image reconstruction for the ECAT EXACT HR+ , 1997, 1997 IEEE Nuclear Science Symposium Conference Record.

[3]  A. Macovski,et al.  A Natural Pixel Decomposition for Two-Dimensional Image Reconstruction , 1981, IEEE Transactions on Biomedical Engineering.

[4]  Gengsheng Lawrence Zeng,et al.  Projection space image reconstruction using strip functions to calculate pixels more "natural" for modeling the geometric response of the SPECT collimator , 1998, IEEE Transactions on Medical Imaging.

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

[6]  Robert M. Lewitt,et al.  Overview of methods for image reconstruction from projections in emission computed tomography , 2003, Proc. IEEE.

[7]  S. Matej,et al.  System and Gram matrices of 3D planogram data , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[8]  Michel Koole,et al.  A 3D model incorporating spatial detection uncertainty in continuous detector PET: Monte Carlo evaluation of a theoretical study , 2003 .

[9]  S. Matej,et al.  System and gram matrices of 3-D planogram data , 2004, IEEE Transactions on Nuclear Science.

[10]  Simon R. Cherry,et al.  Fully 3D Bayesian image reconstruction for the ECAT EXACT HR , 1997 .