Abstract
The popular Radon transform approximation used in the modeling and reconstruction of positron emission tomography (PET) images fails to account for the non-trivial size of PET detectors. Currently, all reconstruction algorithms which account for detector width are based on the iterative EMML (expectation maximization maximum likelihood) method which uses computational approximations of the true point spread function. Starting with the angle-of-view model for the point spread function, a new integral equation model is obtained by applying a simple linear transformation to the emission sinogram which produces detections in arcs of various widths, instead of the usual detector tubes. In this paper, an exact mathematical representation for the new point spread function is obtained in terms of line integrals involving the Green's function and Poisson kernel for the unit disk. By applying appropriate mathematical transforms, this representation leads to a method for reconstructing two dimensional PET images in terms of an orthogonal basis consisting of tensor products of classical harmonics and Bessel functions. The paper contains details of the mathematical derivation of the representation for the point spread function, and the resulting reconstruction technique. It also contains numerical results which indicate that this new method produces images which are comparable (and sometimes superior) in quality to the EMML method, with a reconstruction speed which is similar that of the filtered back projection method.
Similar content being viewed by others
References
J.M.M. Anderson, B.A. Mair, and M. Rao, “Positron emission tomography: A review,” Jour. ofControl & Cybernetics, Vol. 25, pp. 1089–1111, 1996.
A. Bjorck, “Least squares methods,” inHandbook of Numerical Analysis, P. Ciarlet and J. Lions (Eds.), Vol. 1, Elsevier: Amsterdam, 1990.
F. Bowman, Introduction to Bessel Functions, Dover Publications: New York, 1958.
G. Brix, J. Doll, M.E. Bellemann, H. Trojan, U. Haberkorn, P. Schmidlin, and H. Ostertag, “Use of scanner characteristics in iterative image reconstruction for high resolution positron emission tomography studies of small animals,” European Journal of Nuclear Medicine, Vol. 24, pp. 779–786, 1997.
R.B. Carroll, “An orthogonal series approach to positron emissiontomography,” Ph.D. Thesis, University of Florida, Gainesville, FL, Oct. 1998.
S.R. Deans, The RadonTransform and Some of Its Applications, John Wiley & Sons: New York, 1983.
J.A. Fessler and A.O. Hero,“Penalized maximum-likelihood image reconstruction using space-alternating generalized EM algorithms,” IEEE Trans. Image Proc., Vol. 4, pp. 1417–1429, 1995.
J.B. Garnett, Bounded Analytic Functions, Academic Press: NewYork, 1981.
W. Gautschi, “Attenuation factors in practical Fourier analysis,” Numer. Math., Vol. 18, pp.373–400, 1972.
S. Geman and D.E. McClure, “Bayesian image analysis: An application to single photonemission tomography,” in Proc. Statist. Comput. Sect., Amer. Statist. Assn., 1985, pp. 12–18.
L.L. Helms,Introduction to Potential Theory, Wiley-Interscience: New York, 1969.
G.T. Herman, Image Reconstructionfrom Projections: The Fundamentals of Computerized Tomography, Academic Press: San Francisco, 1980.
P.E.Kinahan, S. Matej, J.S. Karp, G.T. Herman, and R.M. Lewitt, “A comparison of transform and iterative reconstruction techniques for a volume–imaging pet scanner with a large axial acceptance angle,” IEEE Trans. Nucl. Science, Vol. 42, pp. 2281–2287, 1995.
R. Leahy and X. Yan, “Incorporation of anatomical MR data for improved functional imaging withPET,” in International Conf. on Information Processing in Medical Imaging, 1991, pp. 105–120.
E. Levitanand G.T. Herman, “A maximum a posteriori probability expectation maximization algorithm for reconstruction in emission tomography,” IEEE Trans. on Medical Imaging, Vol. 6, pp. 185–192, 1987.
R.M. Lewitt,“Multidimensional digital image representations using generalized Kaiser-Bessel window functions,” J. Opt. Soc. Am. A, Vol. 7, pp. 1834–1846, 1990.
R.M. Lewitt, “Alternatives to voxels for image representation in iterativereconstruction algorithms,” Phys. Med. Biol., Vol. 37, pp. 705–716, 1992.
B.A. Mair, “A mathematicalmodel incorporating the effects of detector width in 2D PET,” Inverse Problems, Vol. 16, pp. 223–246, 2000.
S. Matej and R.M. Lewitt, “Practical considerations for 3-D image reconstruction using spherically symmetric volumeelements,” IEEE Trans. Med. Imag., Vol. 15, pp. 68–78, 1996.
E.U. Mumcuoğlu, R. Leahy, S.R. Cherry,and Z. Zhou, “Fast gradient-based methods for Bayesian reconstruction of transmission and emission PET images,” IEEE Trans. Med. Imag., Vol. 13, pp. 687–701, 1994.
F. Natterer, The Mathematics of ComputerizedTomography, Wiley: New York, 1986.
R. Plato, “Optimal algorithms for linear ill-posed problems yieldregularization methods,” Numer. Funct. Anal. Optimiz., Vol. 11, pp. 111–118, 1990.
A.J. Reader, K.Erlandsson, M.A. Flower, and R.J. Ott, “Fast accurate iterative reconstruction for low-statistics positron volume imaging,” Phys. Med. Biol., Vol. 43, pp. 835–846, 1998.
T.J. Rivlin, Chebyshev Polynomials: From ApproximationTheory to Algebra and Number Theory, 2nd edn., John Wiley & Sons: New York, 1990.
L.A. Shepp and Y.Vardi, “Maximum likelihood reconstruction for emission tomography,” IEEE Trans. on Medical Imaging, Vol. 1, pp. 113–122, 1982.
B.W. Silverman, M.C. Jones, D.W. Nycha, and J.D. Wilson, “A smoothed EM approach toindirect estimation problems, with particular reference to stereology and emission tomography,” J. R. Statist. Soc., Series B, Vol. 52, pp. 271–324, 1990.
A.V.D. Sluis and H.A.V.D. Vorst, “SIRT and CG methods for the iterative solutionsof sparse linear least squares problems,” Lin. Alg. Appl., Vol. 130, pp. 76–83, 1990.
D.L. Snyder and M.Miller, “The use of sieves to stabilize images produced with EM algorithm for emission tomography,” IEEE Trans. on Nucl. Sci., Vol. 32, pp. 3864–3872, 1985.
A. Terstegge, S. Weber, H. Herzog, H.W. Müller-Gärtner, and H.Halling, “Resolution and better quantification by tube of response modelling in 3D PET reconstruction,” in 1996 IEEE Nuclear Science Symposium, 1996.
P. Toft and J.J. Jensen, “A very fast implementation of 2D iterative reconstructionalgorithms,” in 1996 IEEE Nuclear Science Symposium, 1996, Vol. 3, pp. 1742–1746.
Y. Vardi, L.A.Shepp, and L. Kaufman, “A statistical model for positron emission tomography,” Journal of the American Statistical Association, Vol. 80, pp. 8–37, 1985.
C.R. Vogel, “Solving ill-conditioned linear systems using theconjugate gradient method,” Department of Mathematical Sciences, Montana State University, Tech. Report, 1987.
G.N. Watson, A Treatise on the Theory of Bessel Functions, Cambridge University Press: London, 1962.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Carroll, R., Mair, B. A New Model and Reconstruction Method for 2D PET Based on Transforming Detector Tube Data into Detector Arc Data. Journal of Mathematical Imaging and Vision 14, 165–185 (2001). https://doi.org/10.1023/A:1011263332663
Issue Date:
DOI: https://doi.org/10.1023/A:1011263332663