Reconstruction From Uniformly Attenuated SPECT Projection Data Using the DBH Method

An algorithm was developed for the 2-D reconstruction of truncated and nontruncated uniformly attenuated data acquired from single photon emission computed tomography (SPECT). The algorithm is able to reconstruct data from half-scan (180deg) and short-scan (180deg+fan angle) acquisitions for parallel- and fan-beam geometries, respectively, as well as data from full-scan (360deg) acquisitions. The algorithm is a derivative, backprojection, and Hilbert transform (DBH) method, which involves the backprojection of differentiated projection data followed by an inversion of the finite weighted Hilbert transform. The kernel of the inverse weighted Hilbert transform is solved numerically using matrix inversion. Numerical simulations confirm that the DBH method provides accurate reconstructions from half-scan and short-scan data, even when there is truncation. However, as the attenuation increases, finer data sampling is required.

[1]  Michael A. King,et al.  Quantitative myocardial perfusion SPECT , 1998, Journal of nuclear cardiology : official publication of the American Society of Nuclear Cardiology.

[2]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .

[3]  D P Harrington,et al.  The inversion of the exponential Radon transform for quantitative brain SPECT. , 1996, Physics in medicine and biology.

[4]  R. Clackdoyle,et al.  Quantitative reconstruction from truncated projections in classical tomography , 2004, IEEE Transactions on Nuclear Science.

[5]  K Kose,et al.  Image reconstruction algorithm for single-photon-emission computed tomography with uniform attenuation. , 1989, Physics in medicine and biology.

[6]  H H Barrett,et al.  Objective assessment of image quality: effects of quantum noise and object variability. , 1990, Journal of the Optical Society of America. A, Optics and image science.

[7]  V. Hutson.  Integral Equations , 1967, Nature.

[8]  Donald C. Solmon,et al.  Inversion of the exponential X‐ray transform. I: Analysis , 1988 .

[9]  Hiroyuki Kudo,et al.  Truncated Hilbert transform and image reconstruction from limited tomographic data , 2006 .

[10]  Zhengrong Liang,et al.  A unified reconstruction framework for both parallel-beam and variable focal-length fan-beam collimators by a Cormack-type inversion of exponential Radon transform , 1999, IEEE Transactions on Medical Imaging.

[11]  Xiaochuan Pan,et al.  A unified analysis of exact methods of inverting the 2-D exponential radon transform, with implications for noise control in SPECT , 1995, IEEE Trans. Medical Imaging.

[12]  M. Defrise,et al.  Image reconstruction from truncated data in SPECT with uniform attenuation , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[13]  Alexander Katsevich,et al.  An improved exact filtered backprojection algorithm for spiral computed tomography , 2004, Adv. Appl. Math..

[14]  Frédéric Noo,et al.  Image reconstruction in 2D SPECT with 180° acquisition , 2001 .

[15]  J. H. Zar,et al.  Biostatistical Analysis (5th Edition) , 1984 .

[16]  R. Novikov.  An inversion formula for the attenuated X-ray transformation , 2002 .

[17]  T. Turkington,et al.  Simultaneous compensation for attenuation, scatter and detector response for SPECT reconstruction in three dimensions. , 1992, Physics in medicine and biology.

[18]  I. S. Gradshteyn,et al.  Table of Integrals, Series, and Products , 1976 .

[19]  B.M.W. Tsui,et al.  Three-dimensional Iterative Reconstruction Algorithms With Attenuation And Geometric Point Response Correction , 1990, 1990 IEEE Nuclear Science Symposium Conference Record.

[20]  C. Metz,et al.  The exponential Radon transform , 1980 .

[21]  William H. Press,et al.  Numerical recipes in C , 2002 .

[22]  D. Gilland,et al.  Implementation of simultaneous attenuation and detector response correction in SPECT , 1988 .

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

[24]  C E Floyd,et al.  Nonisotropic attenuation in SPECT: phantom tests of quantitative effects and compensation techniques. , 1987, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[25]  Jinghan Ye,et al.  Transmission scan truncation with small-field-of-view dedicated cardiac SPECT systems: Impact and automated quality control , 2005, Journal of nuclear cardiology : official publication of the American Society of Nuclear Cardiology.

[26]  J R Perry,et al.  Correction of nonuniform attenuation in cardiac SPECT imaging. , 1989, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[27]  F. Noo,et al.  A two-step Hilbert transform method for 2D image reconstruction. , 2004, Physics in medicine and biology.

[28]  T. J. Breen,et al.  Biostatistical Analysis (2nd ed.). , 1986 .

[29]  Donald C. Solmon,et al.  Filtered-backprojection and the exponential Radon transform , 1989 .

[30]  Zheng Chang,et al.  Problems created in attenuation-corrected SPECT images by artifacts in attenuation maps: a simulation study. , 2005, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[31]  L. Delves,et al.  Computational methods for integral equations: Frontmatter , 1985 .

[32]  Joseph Lipka,et al.  A Table of Integrals , 2010 .

[33]  D P Harrington,et al.  An analytical approach to quantitative reconstruction of non-uniform attenuated brain SPECT. , 1994, Physics in medicine and biology.

[34]  Zhengrong Liang,et al.  Quantitative Brain SPECT in three Dimensions , 1996 .

[35]  Theory for image reconstruction from divergent-beam projections in SPECT , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[36]  Jiangsheng You,et al.  Explicit finite inverse Hilbert transforms , 2006 .

[37]  A. Katsevich.  Analysis of an exact inversion algorithm for spiral cone-beam CT. , 2002, Physics in medicine and biology.

[38]  S. Bellini,et al.  Compensation of tissue absorption in emission tomography , 1979 .

[39]  Gengsheng Lawrence Zeng,et al.  Hilbert Transform Based FBP Algorithm for Fan-Beam CT Full and Partial Scans , 2007, IEEE Transactions on Medical Imaging.

[40]  Thomas F. Budinger,et al.  The Use of Filtering Methods to Compensate for Constant Attenuation in Single-Photon Emission Computed Tomography , 1981, IEEE Transactions on Biomedical Engineering.

[41]  Xiaochuan Pan,et al.  A family of π-scheme exponential Radon transforms and the uniqueness of their inverses , 2002 .

[42]  Andrew Markoe,et al.  Fourier Inversion of the Attenuated X-Ray Transform , 1984 .

[43]  An inversion method for an attenuated x-ray transform , 1996 .

[44]  H. Rullg̊ard.  An explicit inversion formula for the exponential Radon transform using data from 180° , 2004 .

[45]  Xiaochuan Pan,et al.  Exact image reconstruction on PI-lines from minimum data in helical cone-beam CT. , 2004, Physics in medicine and biology.

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

[47]  Gengsheng Lawrence Zeng,et al.  An Analytical Inversion of the 180degree Exponential Radon Transform with a Numerically Generated Kernel , 2007, Int. J. Image Graph..

[48]  N. Hassani,et al.  Principles of computerized tomography. , 1976, Journal of the National Medical Association.

[49]  N. Pelc,et al.  An attenuated projector-backprojector for iterative SPECT reconstruction. , 1985, Physics in medicine and biology.

[50]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[51]  Grant T. Gullberg,et al.  Two finite inverse Hilbert transform formulae for region-of-interest tomography , 2007 .

[52]  K. Lange,et al.  EM reconstruction algorithms for emission and transmission tomography. , 1984, Journal of computer assisted tomography.

[53]  K. Atkinson,et al.  A survey of numerical methods for the solution of Fredholm integral equations of the second kind , 1977 .

[54]  B. Tsui,et al.  Noise properties of filtered-backprojection and ML-EM reconstructed emission tomographic images , 1992, IEEE Conference on Nuclear Science Symposium and Medical Imaging.

[55]  A. Puro.  Cormack-type inversion of exponential Radon transform , 2001 .