Application of Krylov subspaces to SPECT imaging

The application of the conjugate gradient (CG) algorithm to the problem of data reconstruction in SPECT imaging indicates that most of the useful information is already contained in Krylov subspaces of small dimension, ranging from 9 (two‐dimensional case) to 15 (three‐dimensional case). On this basis, a new, proposed approach can be basically summarized as follows: construction of a basis spanning a Krylov subspace of suitable dimension and projection of the projector–backprojector matrix (a 106 × 106 matrix in the three‐dimensional case) onto such a subspace. In this way, one is led to a problem of low dimensionality, for which regularized solutions can be easily and quickly obtained. The required SPECT activity map is expanded as a linear combination of the basis elements spanning the Krylov subspace and the regularization acts by modifying the coefficients of such an expansion. By means of a suitable graphical interface, the tuning of the regularization parameter(s) can be performed interactively on the basis of the visual inspection of one or some slices cut from a reconstruction. © 2003 Wiley Periodicals, Inc. Int J Imaging Syst Technol 12, 217–228, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/ima.10026