Analytical algorithm for system matrix based on Monte Carlo simulation

Songyun Ling,Wentian Cao,Jun Li,Baiyin Cai,Shanglian Bao
DOI: https://doi.org/10.19650/j.cnki.cjsi.2009.09.001
2009-01-01
Abstract:System matrix is the key element in the iterative reconstruction algorithm for SPECT (single photon emission computed tomography). Monte Carlo (MC) simulation is the most accurate and common method to calculate system matrix. However, the major disadvantage of the standard MC algorithm is the long computing time it may take. Therefore, a new method is proposed here, which adopts the projection data attained from MC simulation to fit a point spread function (PSF) of the collimator. Then this function is taken as the analytical expression to calculate the system matrix. This new method not only shortens the computing time, but also contains the MC information, which makes the system matrix calculated using this analytical method comparable with the result of MC simulation. Thus the method can correct the influence of splay angle of the collimator on the calculation of system matrix. The system matrix calculated using the presented method was applied to reconstruct the MC projection of a flood source phantom and a Jaszczak phantom. The reconstructed images recover the primary images very well.
What problem does this paper attempt to address?