Implementation of a fully 3D system model for brain SPECT with fan- beam-collimator OSEM reconstruction with 3D total variation regularization

In order to improve tomographically reconstructed image quality, we have implemented a fully 3D reconstruction, using an ordered subsets expectation maximization (OSEM) algorithm for fan-beam collimator (FBC) SPECT, along with a volumetric system model-fan-volume system model (FVSM), a modified attenuation compensation, a 3D depth- and angle-dependent resolution and sensitivity correction, and a 3D total variation (TV) regularization. SPECT data were acquired in a 128x64 matrix, in 120 views with a circular orbit. The numerical Zubal brain phantom was used to simulate a FBC HMPAO Tc-99m brain SPECT scan, and a low noise and scatter-free projection dataset was obtained using the SimSET Monte Carlo package. A SPECT scan for a mini-Defrise phantom and brain HMPAO SPECT scans for five patients were acquired with a triple-head gamma camera (Triad 88) equipped with a low-energy high-resolution (LEHR) FBC. The reconstructed images, obtained using clinical filtered back projection (FBP), OSEM with a line-length system model (LLSM) and 3D TV regularization, and OSEM with FVSM and 3D TV regularization were quantitatively studied. Overall improvement in the image quality has been observed, including better axial and transaxial resolution, better integral uniformity, higher contrast-to-noise ration between the gray matter and the white matter, and better accuracy and lower bias in OSEM-FVSM, compared with OSEM-LLSM and clinical FBP.

[1]  M. Persson,et al.  Total variation norm for three-dimensional iterative reconstruction in limited view angle tomography , 2001, Physics in medicine and biology.

[2]  P B Hoffer,et al.  Computerized three-dimensional segmented human anatomy. , 1994, Medical physics.

[3]  Fan beam image reconstruction using generalized natural pixel bases [SPECT] , 1995, 1995 IEEE Nuclear Science Symposium and Medical Imaging Conference Record.

[4]  S. B. Lo,et al.  Strip and line path integrals with a square pixel matrix: a unified theory for computational CT projections. , 1988, IEEE transactions on medical imaging.

[5]  Grant T. Gullberg,et al.  Total variation regulated EM algorithm , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

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

[7]  L Johansson,et al.  A method for attenuation and scatter correction of brain SPECT based on computed tomography images , 2003, Nuclear medicine communications.

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

[9]  B C Lentle,et al.  Interactive three-dimensional region of interest analysis of HMPAO SPECT brain studies. , 1990, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.