<abstract xml:lang="en" xmlns="http://www.w3.org/1999/xhtml">A method for system matrix calculation in the case of iterative reconstruction algorithms in SPECT was implemented and tested. Due to a complex mathematical description of the geometry of the detector set-up, we developed a method for system matrix computation that is based on direct measurements of the detector response. In this approach, the influence of the acquisition equipment on the image formation is measured directly. The objective was to obtain the best quality of reconstructed images with respect to specified measures. This is indispensable in order to be able to perform reliable quantitative analysis of SPECT images. It is also especially important in non-hybrid gamma cameras, where not all physical processes that disturb image acquisition can be easily corrected. Two experiments with an 131<italic>I</italic> point source placed at different distances from the detector plane were performed allowing the detector response to be acquired as a function of the point source distance. An analytical Gaussian function was fitted to the acquired data in both the one- and the two-dimensional case. A cylindrical phantom filled with a water solution of 131<italic>I</italic> containing a region of "cold" spheres as well as a uniform solution (without any spheres) was used to perform algorithm evaluation. The reconstructed images obtained by using four different of methods system matrix computation were compared with those achieved using reconstruction software implemented in the gamma camera. The contrast of the spheres and uniformity were compared for each reconstruction result and also with the ranges of those values formulated by the AAPM (American Association of Physicists in Medicine). The results show that the implementation of the OSEM (Ordered Subsets Expectation Maximization) algorithm with a one-dimensional fit to the Gaussian CDR (Collimator-Detector Response) function provides the best results in terms of adopted measures. However, the fit of the two-dimensional function also gives satisfactory results. Furthermore, the CDR function has the potential to be applied to a fully 3D OSEM implementation. The lack of the CDR in system matrix calculation results in a very noisy image that cannot be used for diagnostic purposes.</abstract>

A method for system matrix calculation in the case of iterative reconstruction algorithms in SPECT was implemented and tested. Due to a complex mathematical description of the geometry of the detector set-up, we developed a method for system matrix computation that is based on direct measurements of the detector response. In this approach, the influence of the acquisition equipment on the image formation is measured directly. The objective was to obtain the best quality of reconstructed images with respect to specified measures. This is indispensable in order to be able to perform reliable quantitative analysis of SPECT images. It is also especially important in non-hybrid gamma cameras, where not all physical processes that disturb image acquisition can be easily corrected. Two experiments with an 131I point source placed at different distances from the detector plane were performed allowing the detector response to be acquired as a function of the point source distance. An analytical Gaussian function was fitted to the acquired data in both the one- and the two-dimensional case. A cylindrical phantom filled with a water solution of 131I containing a region of "cold" spheres as well as a uniform solution (without any spheres) was used to perform algorithm evaluation. The reconstructed images obtained by using four different of methods system matrix computation were compared with those achieved using reconstruction software implemented in the gamma camera. The contrast of the spheres and uniformity were compared for each reconstruction result and also with the ranges of those values formulated by the AAPM (American Association of Physicists in Medicine). The results show that the implementation of the OSEM (Ordered Subsets Expectation Maximization) algorithm with a one-dimensional fit to the Gaussian CDR (Collimator-Detector Response) function provides the best results in terms of adopted measures. However, the fit of the two-dimensional function also gives satisfactory results. Furthermore, the CDR function has the potential to be applied to a fully 3D OSEM implementation. The lack of the CDR in system matrix calculation results in a very noisy image that cannot be used for diagnostic purposes.

System matrix computation for iterative reconstruction algorithms in SPECT based on direct measurements

Institute of Automatic Control, Silesian University of Technology, ul. Akademicka 16, 44-100 Gliwice, Poland

Nuclear Medicine and Endocrine Oncology Department, Comprehensive Cancer Centre, Maria Skłodowska-Curie Memorial Institute, ul. Wybrzeże Armii Krajowej 15, 44-101 Gliwice, Poland

International Journal of Applied Mathematics and Computer Science

A method for system matrix calculation in the case of iterative reconstruction algorithms in SPECT was implemented and tested. Due to a complex mathematical...

System matrix computation for iterative reconstruction algorithms in SPECT based on direct measurements

Published Online: Mar 28, 2011

Page range: 193 - 202

DOI: https://doi.org/10.2478/v10006-011-0014-1

Keywords
SPECT, iterative reconstruction, quantitative imaging

This content is open access.

System matrix computation for iterative reconstruction algorithms in SPECT based on direct measurements

Published Online: Mar 28, 2011

Page range: 193 - 202

DOI: https://doi.org/10.2478/v10006-011-0014-1

KeywordsSPECT, iterative reconstruction, quantitative imaging

This content is open access.

Keywords
SPECT, iterative reconstruction, quantitative imaging