Please wait a minute...
Chin. Phys. B, 2015, Vol. 24(1): 018702    DOI: 10.1088/1674-1056/24/1/018702
INTERDISCIPLINARY PHYSICS AND RELATED AREAS OF SCIENCE AND TECHNOLOGY Prev   Next  

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
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
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.

Keywords:  point spread function      single photon incidence      system matrix      positron emission tomography  
Received:  18 June 2014      Revised:  11 September 2014      Accepted manuscript online: 
PACS:  87.57.nf (Reconstruction)  
  87.57.uk (Positron emission tomography (PET))  
  87.57.C- (Image quality)  
Fund: 

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

Corresponding Authors:  Wei Long     E-mail:  weil@ihep.ac.cn

Cite this article: 

Fan Xin (樊馨), Wang Hai-Peng (王海鹏), Yun Ming-Kai (贠明凯), Sun Xiao-Li (孙校丽), Cao Xue-Xiang (曹学香), Liu Shuang-Quan (刘双全), Chai Pei (柴培), Li Dao-Wu (李道武), Liu Bao-Dong (刘宝东), Wang Lu (王璐), Wei Long (魏龙) PET image reconstruction with a system matrix containing point spread function derived from single photon incidence response 2015 Chin. Phys. B 24 018702

[1] Gambhir S S, Czermin J, Schwimmer J, Silverman D H, Coleman R E and Phelps M E 2001 J. Nucl. Med 42 1S
[2] Natterer F and Wuebbeling F 2001 Mathematical Methods in Image Reconstruction (Vol. 5) (Philadelphia: SIAM) p. 81
[3] Hudson H and Larkin R 1994 IEEE Trans. Med. Imaging 13 601
[4] National Electrical Manufacturers Association 2001 NEMA Standards Publication NU 2-2001: Performance Measurements of Positron Emission Tomographs (Virginia: National Electrical Manufacturers Association)
[5] Mawlawi O and Townsend D W 2009 Eur J. Nucl. Med. Mol. Imaging 36 S15
[6] Nestle U, Weber W, Hentschel M and Grosu A L 2009 Phys. Med. Biol. 54 R1
[7] Pan T and Mawlawi O 2008 Med. Phys. 35 4955
[8] Wiant D B, Gersh J A, Bennett M C and Bourland J D 2009 Nuclear Science Symposium Conference Record, October 24-November 1, 2009 Orlando, USA, p. 3758
[9] Chien M K, Yun D, Qing G X and Chin T C 2008 IEEE Trans. Med. Imaging 27 1346
[10] Lecomte R, Schmitt D and Lamoureux G 1984 IEEE Trans. Nucl. Sci. NS-31 556
[11] Liang Z 1994 IEEE Trans. Med. Imaging 13 314
[12] Rahmim A, Tang J, Lodge M A, Lashkari S, Ay M R, Lautamaki R, Tsui B M W and Bengel F M 2008 Phys. Med. Biol. 53 5947
[13] Alessio A M, Kinahan P E and Lewellen T K 2006 IEEE Trans. Med. Imaging 25 828837
[14] Mumcuoglu E U, Leahy R M, Cherry S R and Hoffman E 1996 Nuclear Science Symposium, 1996, Conference Record, 1996 IEEE, November 2-9, 1996 Anaheim, USA, p. 1569
[15] Qi J, Leahy R M, Cherry S R, Chatziioannou A and Farquhar T H 1998 Phys. Med. Biol. 43 1001
[16] Alessio A M and Kinahan P E 2008 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro, May 14-17, 2008 Paris, France, p. 1315
[17] Bernardi E D, Zito F and Baselli G 2007 Engineering in Medicine and Biology Society, 2007. EMBS 2007. 29th Annual International Conference of the IEEE, August 22-26, 2007 Lyon, France, p. 6547
[18] Bernardi E D, Zito F, Michelutti L, Mainardi L, Gerundini P and Baselli G 2003 Engineering in Medicine and Biology Society, 2003. EMBS 2003.25th Annual International Conference of the IEEE, September 17-21, 2003 Cancun, Mexico, p. 975
[19] Fin L, Bailly P, Daouk J and Meye M E 2009 Med. Phys. 36 3072
[20] Panin V Y, Kehren F, Michel C and Casey M 2006 IEEE Trans. Med. Imaging 25 907
[21] Panin V Y, Kehren F, Rothfuss H, Hu D, Michel C and Casey M E 2006 IEEE Trans. Nucl. Sci. 53 152
[22] Tohme M S and Qi J 2009 Phys. Med. Biol. 4 3709
[23] Alessio A M, Stearns C W, Tong S, Ross S G, Ganin A and Kinahan P E 2010 IEEE Trans. Med. Imaging 29 938
[24] Qi J, Leahy R, Cherry S R, Chatziioannou A and Farquh-ar T 1998 Phys. Med. Bio. 43 1001
[25] Wang L, Wu L W, Wei L, Gao J, Sun C L Chai P and Li D W 2014 Chin. Phys. B 23 2027802
[26] Badawi R D, Lodge M and Marsden P K 1998 Phys. Med. Bio. 43 189
[27] Tohme M S 2011 Iterative Image Reconstruction for Positron Emission Tomography Based on Measured Detector Response Function Tomography Based on Measured Detector Response Function (Ph. D. dissertation) (California: University of California)
[28] Joseph P M 1982 IEEE Trans. Med. Imaging MI-1 192
[29] Johnson C, Yan Y, Carson R, Martino R and Daube-Witherspoon M 1995 IEEE Trans. Nucl. Sci. 42 1223
[30] Chen C, Lee S and Cho Z 1991 IEEE Trans. Med. Imaging 10 513
[31] http://www.opengatecollaboration.org/
[32] Leonard S 2005 Spatial resolution Study of PET Detector Modules Based on LSO Crystals and Avalanche Photodiode Arrays (Ph.D. dissertation) (Brussel: Vrije Universiteit Brussel)
[33] Defrise M, Kinahan P, Townsend D, Michel C, Sibomana M and Newport D 1997 IEEE Trans. Med. Imaging 16 145
[1] Reduced aliasing artifacts using shaking projection k-space sampling trajectory
Zhu Yan-Chun (朱艳春), Du Jiang (杜江), Yang Wen-Chao (杨文超), Duan Chai-Jie (段侪杰), Wang Hao-Yu (王浩宇), Gao Song (高嵩), Bao Shang-Lian (包尚联). Chin. Phys. B, 2014, 23(3): 038702.
[2] An effective scatter correction method based on single scatter simulation for a 3D whole-body PET scanner
Gao Fei(高飞), Yamada Ryoko, Watanabe Mitsuo, and Liu Hua-Feng(刘华锋). Chin. Phys. B, 2009, 18(7): 3066-3072.
No Suggested Reading articles found!