中国物理B ›› 2015, Vol. 24 ›› Issue (1): 18702-018702.doi: 10.1088/1674-1056/24/1/018702

• INTERDISCIPLINARY PHYSICS AND RELATED AREAS OF SCIENCE AND TECHNOLOGY • 上一篇    下一篇

PET image reconstruction with a system matrix containing point spread function derived from single photon incidence response

樊馨a b, 王海鹏a b, 贠明凯a c, 孙校丽a b, 曹学香a c, 刘双全a c, 柴培a c, 李道武a c, 刘宝东a c, 王璐a b, 魏龙a c   

  1. a Key Laboratory of Nuclear Analytical Techniques, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China;
    b University of Chinese Academy of Sciences, Beijing 100049, China;
    c Beijing Engineering Research Center of Radiographic Techniques and Equipment, Beijing 100049, China
  • 收稿日期:2014-06-18 修回日期:2014-09-11 出版日期:2015-01-05 发布日期:2015-01-05
  • 基金资助:

    Project supported by the National Natural Science Foundation of China (Grant Nos. Y4811H805C and 81101175).

PET image reconstruction with a system matrix containing point spread function derived from single photon incidence response

Fan Xin (樊馨)a b, Wang Hai-Peng (王海鹏)a b, Yun Ming-Kai (贠明凯)a c, Sun Xiao-Li (孙校丽)a b, Cao Xue-Xiang (曹学香)a c, Liu Shuang-Quan (刘双全)a c, Chai Pei (柴培)a c, Li Dao-Wu (李道武)a c, Liu Bao-Dong (刘宝东)a c, Wang Lu (王璐)a b, Wei Long (魏龙)a c   

  1. a Key Laboratory of Nuclear Analytical Techniques, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China;
    b University of Chinese Academy of Sciences, Beijing 100049, China;
    c Beijing Engineering Research Center of Radiographic Techniques and Equipment, Beijing 100049, China
  • Received:2014-06-18 Revised:2014-09-11 Online:2015-01-05 Published:2015-01-05
  • Contact: Wei Long E-mail:weil@ihep.ac.cn
  • Supported by:

    Project supported by the National Natural Science Foundation of China (Grant Nos. Y4811H805C and 81101175).

摘要:

A point spread function (PSF) for the blurring component in positron emission tomography (PET) is studied. The PSF matrix is derived from the single photon incidence response function. A statistical iterative reconstruction (IR) method based on the system matrix containing the PSF is developed. More specifically, the gamma photon incidence upon a crystal array is simulated by Monte Carlo (MC) simulation, and then the single photon incidence response functions are calculated. Subsequently, the single photon incidence response functions are used to compute the coincidence blurring factor according to the physical process of PET coincidence detection. Through weighting the ordinary system matrix response by the coincidence blurring factors, the IR system matrix containing the PSF is finally established. By using this system matrix, the image is reconstructed by an ordered subset expectation maximization (OSEM) algorithm. The experimental results show that the proposed system matrix can substantially improve the image radial resolution, contrast, and noise property. Furthermore, the simulated single gamma-ray incidence response function depends only on the crystal configuration, so the method could be extended to any PET scanners with the same detector crystal configuration.

关键词: point spread function, single photon incidence, system matrix, positron emission tomography

Abstract:

A point spread function (PSF) for the blurring component in positron emission tomography (PET) is studied. The PSF matrix is derived from the single photon incidence response function. A statistical iterative reconstruction (IR) method based on the system matrix containing the PSF is developed. More specifically, the gamma photon incidence upon a crystal array is simulated by Monte Carlo (MC) simulation, and then the single photon incidence response functions are calculated. Subsequently, the single photon incidence response functions are used to compute the coincidence blurring factor according to the physical process of PET coincidence detection. Through weighting the ordinary system matrix response by the coincidence blurring factors, the IR system matrix containing the PSF is finally established. By using this system matrix, the image is reconstructed by an ordered subset expectation maximization (OSEM) algorithm. The experimental results show that the proposed system matrix can substantially improve the image radial resolution, contrast, and noise property. Furthermore, the simulated single gamma-ray incidence response function depends only on the crystal configuration, so the method could be extended to any PET scanners with the same detector crystal configuration.

Key words: point spread function, single photon incidence, system matrix, positron emission tomography

中图分类号:  (Reconstruction)

  • 87.57.nf
87.57.uk (Positron emission tomography (PET)) 87.57.C- (Image quality)