Fully 3D PET image reconstruction using a fourier preconditioned conjugate-gradient algorithm
- Univ. of Michigan, Ann Arbor, MI (United States)
Since the data sizes in fully 3D PET imaging are very large, iterative image reconstruction algorithms must converge in very few iterations to be useful. One can improve the convergence rate of the conjugate-gradient (CG) algorithm by incorporating preconditioning operators that approximate the inverse of the Hessian of the objective function. If the 3D cylindrical PET geometry were not truncated at the ends, then the Hessian of the penalized least-squares objective function would be approximately shift-invariant, i.e. G{prime}G would be nearly block-circulant, where G is the system matrix. We propose a Fourier preconditioner based on this shift-invariant approximation to the Hessian. Results show that this preconditioner significantly accelerates the convergence of the CG algorithm with only a small increase in computation.
- OSTI ID:
- 513259
- Report Number(s):
- CONF-961123--; CNN: Grant CA-60711; Grant CA-54362
- Country of Publication:
- United States
- Language:
- English
Similar Records
Maximum likelihood reconstruction in fully 3D PET via the SAGE algorithm
Parallel preconditioning techniques for sparse CG solvers