A (μ + λ) evolutionary algorithm for reconstruction of SPECT data

Abstract

Algorithms used to reconstruct single photon emission computed tomography (SPECT) data are based on one of two principles: filtered back projection or iterative methods. In this paper, an evolution strategy (ES) was applied to reconstruct transaxial slices of SPECT data. Evolutionary algorithms are stochastic global search methods that have been used successfully for many kinds of optimization problems. The newly developed reconstruction algorithm consisting of /spl mu/ parents and /spl lambda/ children uses a random principle to readjust the voxel values, whereas other iterative reconstruction methods use the difference between measured and simulated projection data. The (/spl mu/ + /spl lambda/)-ES was validated against a test image, a heart, and a Jaszczak phantom. The resulting transaxial slices show an improvement in image quality, in comparison to both the filtered back projection method and a standard iterative reconstruction algorithm.