Two-dimensional statistical model for regularized backprojection in SPECT.

In SPECT, both the noise affecting the data and the discretization of the inverse Radon transform are responsible for the ill-posed nature of the reconstruction. To constrain the problem, we propose a regularized backprojection method (RBP) which takes advantage of the relationships existing between the continuity properties of the projections and those of the reconstructed object. The RBP method involves two stages: first, a statistical model (the fixed-effect model) is used to estimate the noise-free part of the projections. Then, the filtered projections are reconstructed using a backprojection algorithm (spline filtered backprojection) which ensures that the reconstructed object belongs to a space consistent with that containing the projections. The method is illustrated using analytical simulations, and the RBP approach is compared to the conventional filtered backprojection. The effect on the reconstructed slices of the parameters involved in RBP is studied in terms of spatial resolution, homogeneity in uniform regions and quantification. It is shown that appropriate combinations of these parameters yield a better compromise between homogeneity and spatial resolution than conventional FBP, with similar quantification performances.

[1]  Jean-Pierre V. Guédon,et al.  Bandlimited and Haar filtered back-projection reconstructions , 1994, IEEE Trans. Medical Imaging.

[2]  Alfred K. Louis,et al.  Approximation of the Radon Transform from Samples in Limited Range , 1981 .

[3]  I. J. Schoenberg,et al.  SPLINE FUNCTIONS AND THE PROBLEM OF GRADUATION , 1964 .

[4]  Calyampudi Radhakrishna Rao,et al.  Linear Statistical Inference and its Applications , 1967 .

[5]  S. Wold Cross-Validatory Estimation of the Number of Components in Factor and Principal Components Models , 1978 .

[6]  R. Huesman,et al.  Emission computed tomography , 1979 .

[7]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

[8]  Yves J. Bizais,et al.  Statistical model for tomographic reconstruction methods using spline functions , 1994, Optics & Photonics.

[9]  M. Hutchinson,et al.  Smoothing noisy data with spline functions , 1985 .

[10]  C. R. Rao,et al.  Linear Statistical Inference and its Applications , 1968 .

[11]  G. Wahba Spline models for observational data , 1990 .

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

[13]  E. Hoffman,et al.  Quantitation in Positron Emission Computed Tomography: 1. Effect of Object Size , 1979, Journal of computer assisted tomography.

[14]  P. Besse,et al.  Sur l'usage de la validation croisée en analyse en composantes principales , 1993 .

[15]  B. C. Penney,et al.  Two-dimensional filtering of SPECT images using the Metz and Wiener filters. , 1984, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[16]  C. Reinsch Smoothing by spline functions , 1967 .

[17]  Grant T. Gullberg,et al.  Users Manual: Donner Algorithms for Reconstruction Tomography , 1977 .

[18]  Jean-Pierre V. Guédon,et al.  Spline-Based Regularization for Discrete FBP Reconstruction , 1991, IPMI.

[19]  Peter Craven,et al.  Smoothing noisy data with spline functions , 1978 .