Enhanced 3D PET OSEM reconstruction using inter-update Metz filtering.

We present an enhancement of the OSEM (ordered set expectation maximization) algorithm for 3D PET reconstruction, which we call the inter-update Metz filtered OSEM (IMF-OSEM). The IMF-OSEM algorithm incorporates filtering action into the image updating process in order to improve the quality of the reconstruction. With this technique, the multiplicative correction image--ordinarily used to update image estimates in plain OSEM--is applied to a Metz-filtered version of the image estimate at certain intervals. In addition, we present a software implementation that employs several high-speed features to accelerate reconstruction. These features include, firstly, forward and back projection functions which make full use of symmetry as well as a fast incremental computation technique. Secondly, the software has the capability of running in parallel mode on several processors. The parallelization approach employed yields a significant speed-up, which is nearly independent of the amount of data. Together, these features lead to reasonable reconstruction times even when using large image arrays and non-axially compressed projection data. The performance of IMF-OSEM was tested on phantom data acquired on the GE Advance scanner. Our results demonstrate that an appropriate choice of Metz filter parameters can improve the contrast-noise balance of certain regions of interest relative to both plain and post-filtered OSEM, and to the GE commercial reprojection algorithm software.

[1]  T. Hebert,et al.  A generalized EM algorithm for 3-D Bayesian reconstruction from Poisson data using Gibbs priors. , 1989, IEEE transactions on medical imaging.

[2]  Gengsheng L. Zeng,et al.  A study of reconstruction artifacts in cone beam tomography using filtered backprojection and iterative EM algorithms , 1990 .

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

[4]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[5]  Anni Cai,et al.  Incremental backprojection algorithm: modification of the searching flow scheme and utilization of the relationship among projection views , 1993, IEEE Trans. Medical Imaging.

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

[7]  T. Hebert,et al.  Fast methods for including attenuation in the EM algorithm , 1990 .

[8]  Hakan Erdogan,et al.  Ordered subsets algorithms for transmission tomography. , 1999, Physics in medicine and biology.

[9]  Patrick Dupont,et al.  Evaluation of maximum-likelihood based attenuation correction in positron emission tomography , 1998 .

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

[11]  H. Murayama,et al.  Performance of a new 3D-only PET scanner-the EXACT3D , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[12]  C. Comtat,et al.  Fast reconstruction of 3-D PET data with accurate statistical modeling , 1997 .

[13]  L. Shepp,et al.  A Statistical Model for Positron Emission Tomography , 1985 .

[14]  J. Ollinger,et al.  Maximum likelihood reconstruction in fully 3D PET via the SAGE algorithm , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[15]  D. Visvikis,et al.  Intercomparison of four reconstruction techniques for positron volume imaging with rotating planar detectors. , 1998, Physics in medicine and biology.

[16]  Vincent J. Cunningham,et al.  Development of 3D dynamic acquisition in a neuro-PET scanner , 1993, 1993 IEEE Conference Record Nuclear Science Symposium and Medical Imaging Conference.

[17]  J. D. Wilson,et al.  A smoothed EM approach to indirect estimation problems, with particular reference to stereology and emission tomography , 1990 .

[18]  C. M. Chen,et al.  Incremental algorithm-a new fast backprojection scheme for parallel beam geometries. , 1990, IEEE transactions on medical imaging.

[19]  S Grootoonk,et al.  A rotating PET scanner using BGO block detectors: design, performance and applications. , 1993, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[20]  M L Egger,et al.  Incremental beamwise backprojection using geometrical symmetries for 3D PET reconstruction in a cylindrical scanner geometry. , 1998, Physics in medicine and biology.

[21]  H. Herzog,et al.  High resolution and better quantification by tube of response modelling in 3D PET reconstruction , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[22]  Kris Thielemans,et al.  On various approximations for the projectors in iterative reconstruction algorithms for 3D-PET , 1999 .

[23]  P K Marsden,et al.  Algorithms for calculating detector efficiency normalization coefficients for true coincidences in 3D PET. , 1998, Physics in medicine and biology.

[24]  Paul Kinahan,et al.  Analytic 3D image reconstruction using all detected events , 1989 .