Si Li1, Jiahan Zhang2, Andrzej Krol3, C Ross Schmidtlein4, David Feiglin5, Yuesheng Xu6. 1. School of Data and Computer Science, Guangdong Province Key Lab of Computational Science, Sun Yat-sen University, Guangzhou 510275, China. Electronic address: lisi23@mail.sysu.edu.cn. 2. Department of Physics, Syracuse University, Syracuse, NY 13244, USA. Electronic address: jiahancheung@gmail.com. 3. Department of Radiology, SUNY Upstate Medical University, Syracuse, NY 13210, USA; Department of Pharmacology, SUNY Upstate Medical University, Syracuse, NY 13210, USA. Electronic address: krola@upstate.edu. 4. Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY 10065, USA. Electronic address: schmidtr@mskcc.org. 5. Department of Radiology, SUNY Upstate Medical University, Syracuse, NY 13210, USA. Electronic address: feiglind@upstate.edu. 6. School of Data and Computer Science, Guangdong Province Key Lab of Computational Science, Sun Yat-sen University, Guangzhou 510275, China; Emeritus Professor, Department of Mathematics, Syracuse University, Syracuse, NY 13244, USA. Electronic address: xuyuesh@mail.sysu.edu.cn.
Abstract
PURPOSE: The authors recently developed a preconditioned alternating projection algorithm (PAPA) for solving the penalized-likelihood SPECT reconstruction problem. The proposed algorithm can solve a wide variety of non-differentiable optimization models. This work is dedicated to comparing the performance of PAPA with total variation (TV) regularization (TV-PAPA) and a novel forward-backward algorithm with nested expectation maximization (EM)-TV iteration scheme (FB-EM-TV). METHODS: Monte Carlo technique was used to simulate multiple noise realizations of the fan-beam collimated SPECT data for a piecewise constant phantom with warm background, and hot and cold spheres with uniform activities at two noise levels. They were reconstructed using the aforementioned algorithms with attenuation, scatter, distance-dependent collimator blurring and sensitivity corrections. Noise suppressing performance, lesion detectability, lesion contrast, contrast recovery coefficient, convergence speed and selection of optimal parameters were evaluated. The conventional EM algorithms with TV post-filter (TVPF-EM) and Gaussian post-filter (GPF-EM) were used as benchmarks. RESULTS: The TV-PAPA and FB-EM-TV demonstrated similar performance in all investigated categories. Both algorithms outperformed TVPF-EM in terms of image noise suppression, lesion detectability, lesion contrast and convergence speed. We established that the optimal parameters versus information density approximately followed power laws, which offers a guidance in parameter selection for reconstruction methods. CONCLUSIONS: For the simulated SPECT data, TV-PAPA and FB-EM-TV produced qualitatively and quantitatively similar images. They performed better than the benchmark TVPF-EM and GPF-EM, with only limited loss of lesion contrast.
PURPOSE: The authors recently developed a preconditioned alternating projection algorithm (PAPA) for solving the penalized-likelihood SPECT reconstruction problem. The proposed algorithm can solve a wide variety of non-differentiable optimization models. This work is dedicated to comparing the performance of PAPA with total variation (TV) regularization (TV-PAPA) and a novel forward-backward algorithm with nested expectation maximization (EM)-TV iteration scheme (FB-EM-TV). METHODS: Monte Carlo technique was used to simulate multiple noise realizations of the fan-beam collimated SPECT data for a piecewise constant phantom with warm background, and hot and cold spheres with uniform activities at two noise levels. They were reconstructed using the aforementioned algorithms with attenuation, scatter, distance-dependent collimator blurring and sensitivity corrections. Noise suppressing performance, lesion detectability, lesion contrast, contrast recovery coefficient, convergence speed and selection of optimal parameters were evaluated. The conventional EM algorithms with TV post-filter (TVPF-EM) and Gaussian post-filter (GPF-EM) were used as benchmarks. RESULTS: The TV-PAPA and FB-EM-TV demonstrated similar performance in all investigated categories. Both algorithms outperformed TVPF-EM in terms of image noise suppression, lesion detectability, lesion contrast and convergence speed. We established that the optimal parameters versus information density approximately followed power laws, which offers a guidance in parameter selection for reconstruction methods. CONCLUSIONS: For the simulated SPECT data, TV-PAPA and FB-EM-TV produced qualitatively and quantitatively similar images. They performed better than the benchmark TVPF-EM and GPF-EM, with only limited loss of lesion contrast.