Conjugate-gradient preconditioning methods for shift-variant PET image reconstruction

IEEE Transactions on Image Processing

Volumn 8, Issue 5, 1999, Pages 688-699

  Fessler, Jeffrey A ^a,c   Booth, Scott D ^b  

^a UNIVERSITY OF MICHIGAN   (United States)

^b UNIVERSITY OF VIRGINIA MEDICAL CENTER   (United States)

^c IEEE   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

CONJUGATE GRADIENT PRECONDITIONING METHODS; HESSIAN MATRICES;

CONVERGENCE OF NUMERICAL METHODS; IMAGING TECHNIQUES; INVERSE PROBLEMS; ITERATIVE METHODS; LEAST SQUARES APPROXIMATIONS; MATRIX ALGEBRA; POISSON DISTRIBUTION; POSITRON EMISSION TOMOGRAPHY; SPURIOUS SIGNAL NOISE;

IMAGE RECONSTRUCTION;

EID: 0032655433     PISSN: 10577149     EISSN: None     Source Type: Journal    
DOI: 10.1109/83.760336     Document Type: Article

References (56)

1. E. Artzy, T. Elfving, and G. T. Herman, "Quadratic optimization for image reconstruction, II," Comput. Graph. Image Processing, vol. 11, pp. 242-261, 1979.
2. S. Kawata and O. Nalcioglu, "Constrained iterative reconstruction by the conjugate gradient method," IEEE Trans Med. Imag., vol. 4, pp. 65-71, June 1985.
3. A. R. Formiconi, A. Pupi, and A. Passeri, "Compensation of spatial system response in SPECT with conjugate gradient reconstruction technique," Phys. Med. Biol., vol. 34, pp. 69-84, 1989.
4. L. Kaufman, "Maximum likelihood, least squares, and penalized least squares for PET," IEEE Trans. Med. Imag., vol. 12, pp. 200-214, June 1993.
5. K. Sauer and C. Bouman, "A local update strategy for iterative reconstruction from projections," IEEE Trans. Signal Processing, vol. 41, pp. 534-548, Feb. 1993.
6. E. U. Mumcuoglu, R. Leahy, S. R. Cherry, and Z. Zhou, "Fast gradient-based methods for Bayesian reconstruction of transmission and emission PET images,"IEEE Trans. Med. Imag., vol. 13, pp. 687-701, Dec. 1994.
7. D. S. Lalush and B. M. W. Tsui, "A fast and stable maximum a posteriori conjugate gradient reconstruction algorithm," Med. Phys., vol. 22, pp. 1273-1284, Aug. 1995.
8. S. D. Booth and J. A. Fessler, "Combined diagonal/Fourier preconditioning methods for image reconstruction in emission tomography," in Proc. IEEE Int. Conf. Image Processing, 1995, vol. 2, pp. 441-444.
9. A. H. Delaney and Y. Bresler, "A fast and accurate Fourier algorithm for iterative parallel-beam tomography," IEEE Trans. Image Processing, vol. 5, pp. 740-753, May 1996.
10. J. A. Fessler and E. P. Ficaro, "Fully 3D PET image reconstruction using a Fourier preconditioned conjugate-gradient algorithm," in Proc. IEEE Nuclear Science Symp. Medical Imaging Conf., 1996, vol. 3, pp. 1599-1602.
11. J. A. Fessler, "Preconditioning methods for shift-variant image reconstruction," in Proc. IEEE Int. Conf. Image Processing, 1997, vol. 1, pp. 185-188.
12. O. Axelsson, Iterative Solution Methods, Cambridge, Cambridge, 1994.
13. L. Kaufman, "Implementing and accelerating the EM algorithm for positron emission tomography," IEEE Trans. Med. Imag., vol. 6, pp. 37-51, Mar. 1987.
14. D. Lalush and B. M. W. Tsui, "The importance of preconditioners in fast Poisson based iterative reconstruction algorithms for SPECT," in Proc. IEEE Nuclear Science Symp. Medical Imaging Conf., 1995, vol. 3, pp. 1326-1330.
15. R. H. Chan and M. K. Ng, "Conjugate gradient methods for Toeplitz systems," SIAM Rev., vol. 38, pp. 427-482, Sept. 1996.
16. T. F. Chan, "An optimal circulant preconditioner for Toeplitz systems," SIAM J. Sci. Stat. Comput., vol. 9, pp. 766-771, July 1988.
17. R. H. Chan, J. G. Nagy, and R. J. Plemmons, "FFT-based preconditioners for Toeplitz-block least-squares problems," SIAM J. Numer. Anal., vol. 30, pp. 1740-1768, Dec. 1993.
18. R. H. Chan and C. K. Wong, "Best-conditioned circulant preconditioners," Linear Algebra Applic., vol. 218, pp. 205-211, 1995.
19. R. H. Chan, M. K. Ng, and C. K. Wong, "Sine transform based preconditioners for symmetric Toeplitz systems," Linear Algebra Applicat., vol. 232, pp. 237-259, Jan. 1996.
20. N. H. Clinthorne et al., "Preconditioning methods for improved convergence rates in iterative reconstructions," IEEE Trans. Med. Imag., vol. 12, pp. 78-83, Mar. 1993.
21. R. H. Chan, M. K. Ng, and R. J. Plemmons, "Generalization of Strang's preconditioner with applications to Toeplitz least squares problems," Numer. Linear Algebra Applicat., vol. 3, pp. 45-64, 1996.
22. J. G. Nagy, R. J. Plemmons, and T. C. Torgersen, "Iterative image restoration using approximate inverse preconditioning," IEEE Trans. Image Processing, vol. 5, pp. 1151-1162, July 1996.
23. R. Chan, T. F. Chan, and C. Wong, "Cosine transform based preconditioners for total variation minimization problems in image processing," in Iterative Methods in Linear Algebra, II, V3, IMACS Series in Computational and Applied Mathematics, Proc. 2nd IMACS Int. Symp. Iterative Methods in Linear Algebra, S. Margenov and P. Vassilevski, Eds., IMACS, NJ, 1996, pp. 311-329.
24. C. R. Vogel and M. E. Oman, "Fast numerical methods for total variation minimization in image reconstruction," in Proc. SPIE 2563, Adv. Signal Process. Algorithms, 1995, pp. 359-367.
25. D. L. Snyder, A. M. Hammoud, and R. L. White, "Image recovery from data acquired with a charge-couple-device camera," J. Opt. Soc. Amer. A, vol. 10, pp. 1014-1023, May 1993.
26. D. L. Snyder et al., "Compensation for readout noise in CCD images," J. Opt. Soc. Amer. A, vol. 12, pp. 272-283, Feb. 1995.
27. J. A. Fessler, "Penalized weighted least-squares image reconstruction for positron emission tomography," IEEE Trans. Med. Imag., vol. 13, pp. 290-300, June 1994.
28. J. G. Nagy and D. P. O'Leary, "Restoring images degraded by spatially-variant blur," SIAM J. Sci. Comput., vol. 19, pp. 1063-1082, July 1998.
29. P. J. Green, "Iteratively reweighted least squares for maximum likelihood estimation, and some robust and resistant alternatives," J. R. Stat. Soc. B, vol. 46, pp. 149-192, 1984.
30. J. A. Fessler and W. L. Rogers, "Spatial resolution properties of penalized-likelihood image reconstruction methods: Space-invariant tomographs," IEEE Trans. Image Processing, vol. 5, pp. 1346-1358, Sept. 1996.
31. C. Bouman and K. Sauer, "A generalized Gaussian image model for edge-preserving MAP estimation," IEEE Trans. Image Processing, vol. 2, pp. 296-310, July 1993.
32. A. R. De Pierro, "A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography," IEEE Trans. Med. Imaging, vol. 14, pp. 132-137, Mar. 1995.
33. J. A. Fessler, "Conjugate-gradient preconditioning methods: Numerical results," Tech. Rep. 303, Commun. Signal Process. Lab., Dept. Elect. Eng., Comput. Sci., Univ. Michigan, Ann Arbor, MI, Jan. 1997, available from http://www.eecs.umich.edu/∼fessler.
34. E. U. Mumcuoglu, R. M. Leahy, Z. Zhou, and S. R. Cherry, "A phantom study of the quantitative behavior of Bayesian PET reconstruction methods," in Proc. IEEE Nuclear Science Symp. Medical Imaging Conf., 1995, vol. 3, pp. 1703-1707.
35. E. U. Mumcuoǧlu, R. M. Leahy, and S. R. Cherry, "Bayesian reconstruction of PET images: Methodology and performance analysis," Phys. Med. Biol., vol. 41, pp. 1777-1807, Sept. 1996.
36. E. U. Mumcuoǧlu and R. M. Leahy, "A gradient projection conjugate gradient algorithm for Bayesian PET reconstruction," in Proc. IEEE Nuclear Science Symp. Medical Imaging Conf., 1994, vol. 3, pp. 1212-1216.
37. J. J. Moré and G. Toraldo, "On the solution of large quadratic programming problems with bound constraints," SIAM J. Optim., vol. 1, pp. 93-113, Feb. 1991.
38. M. Bierlaire, P. L. Toint, and D. Tuyttens, "On iterative algorithms for linear least squares problems with bound constraints," Linear Algebra Applicat., vol. 143, pp. 111-143, 1991.
39. W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, Numerical Recipes in C. Cambridge, U.K.: Cambridge Univ. Press, 1988.
40. M. Al-Baali and R. Fletcher, "On the order of convergence of preconditioned nonlinear conjugate gradient methods," SIAM J. Sci. Comput., vol. 17, pp. 658-665, May 1996.
41. A. K. Jain, Fundamentals of Digital Image Processing. Englewood Cliffs, NJ: Prentice-Hall, 1989.
42. P. J. Huber, Robust Statistics. New York: Wiley, 1981.
43. P. Charbonnier, L. Blanc-Féraud, G. Aubert, and M. Barlaud, "Two deterministic half-quadratic regularization algorithms for computed imaging," in Proc. IEEE Int. Conf. Image Processing, 1994, vol. 2, pp. 168-171.
44. D. Geman and G. Reynolds, "Constrained restoration and the recovery of discontinuities," IEEE Trans. Pattern Anal. Machine Intell., vol. 14, pp. 367-383, Mar. 1992.
45. D. Geman and C. Yang, "Nonlinear image recovery with half-quadratic regularization," IEEE Trans. Image Processing, vol. 4, pp. 932-946, July 1995.
46. P. Charbonnier, L. Blanc-Féraud, G. Aubert, and M. Barlaud, "Deterministic edge-preserving regularization in computed imaging," IEEE Trans. Image Processing, vol. 6, pp. 298-311, Feb. 1997.
47. J. A. Fessler, "Grouped coordinate descent algorithms for robust edge-preserving image restoration," in Proc. SPIE 3071, Image Reconst. Restor. II, pp. 184-194, 1997.
48. M. Yavuz and J. A. Fessler, "New statistical models for randoms-precorrected PET scans," in Lecture Notes in Computer Science, Vol. 1230: Information Processing in Medical Imaging, J. Duncan and G. Gindi, Eds. Berlin, Germany: Springer-Verlag, 1997, pp. 190-203.
49. J. A. Fessler, "Resolution properties of regularized image reconstruction methods," Tech. Rep. 297, Commun. Signal Process. Lab., Dept. Elect. Eng. Comput. Sci., Univ. Michigan, Ann Arbor, MI, Aug. 1995.
50. J. A. Fessler, E. P. Ficaro, N. H. Clinthorne, and K. Lange, "Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction," IEEE Trans. Med. Imag., vol. 16, pp. 166-175, Apr. 1997.
51. J. A. Fessler, "Hybrid Poisson/polynomial objective functions for tomographic image reconstruction from transmission scans," IEEE Trans. Image Processing, vol. 4, pp. 1439-1450, Oct. 1995.
52. _, "Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography," IEEE Trans. Image Processing, vol. 5, pp. 493-506, Mar. 1996.
53. _, "ASPIRE 3.0 user's guide: A sparse iterative reconstruction library," Tech. Rep. 293, Commun. Signal Process. Lab., Dept. Elect. Eng. Comput. Sci., Univ. Michigan, Ann Arbor, July 1995.
54. K. Lange, "Convergence of EM image reconstruction algorithms with Gibbs smoothing," IEEE Trans. Med. Imag., vol. 9, pp. 439-446, Dec. 1990.
55. A. O. Hero, M. Usman, A. C. Sauve, and J. A. Fessler, "Recursive algorithms for computing the Cramer-Rao bound," IEEE Trans. Signal Processing, vol. 45, pp. 803-807, Mar. 1997.
56. W. Xia, R. M. Lewitt, and P. R. Edholm, "Fourier correction for spatially variant collimator blurring in SPECT," IEEE Trans. Med. Imag., vol. 14, pp. 100-115, Mar. 1995.