Fast and efficient fully 3D PET image reconstruction using sparse system matrix factorization with GPU acceleration

Physics in Medicine and Biology

Volumn 56, Issue 20, 2011, Pages 6739-6757

  Zhou, Jian ^a   Qi, Jinyi ^a  

^a University of California   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

COMPUTATION COSTS; DATA SPACE; EFFICIENT IMPLEMENTATION; EXPERIMENTAL STUDIES; FACTORIZATION METHODS; FULLY 3-D PET; GEOMETRICAL PROJECTIONS; GRAPHICS PROCESSOR UNITS; HIGH RESOLUTION; IMAGE BLURRING; IMAGE SPACE; ITERATIVE IMAGE RECONSTRUCTION; LIMITED MEMORY; MATRIX; MATRIX SIZE; ORIGINAL SYSTEMS; PROJECTION MATRIX; RECONSTRUCTED IMAGE; SIMPLE LINES; SINOGRAMS; SPARSE MATRICES; SPARSE MATRIX FACTORIZATION; SPARSE SYSTEMS; STORAGE REQUIREMENTS; SYSTEM MATRICES;

COST REDUCTION; FACTORIZATION; IMAGE RECONSTRUCTION; POSITRON EMISSION TOMOGRAPHY; THREE DIMENSIONAL;

MATRIX ALGEBRA;

ACCELERATION; ANIMAL; ARTICLE; COMPUTER; METHODOLOGY; MOUSE; POSITRON EMISSION TOMOGRAPHY; THREE DIMENSIONAL IMAGING; TIME;

ACCELERATION; ANIMALS; COMPUTERS; IMAGING, THREE-DIMENSIONAL; MICE; POSITRON-EMISSION TOMOGRAPHY; TIME FACTORS;

EID: 80053616587     PISSN: 00319155     EISSN: 13616560     Source Type: Journal    
DOI: 10.1088/0031-9155/56/20/015     Document Type: Article

References (34)

1. Alessio A, Kinahan P and Lewellen T 2006 Modeling and incorporation of system response functions in 3-D whole body PET IEEE Trans. Med. Imaging 25 828-37 (Pubitemid 44088368)
2. Alessio A M, Stearns C W, Tong S, Ross S G, Kohlmyer S, Ganin A and Kinahan P E 2010 Application and evaluation of a measured spatially variant system model for PET image reconstruction IEEE Trans. Med. Imaging 29 938-49
3. Bai B and Smith A 2006 Fast 3D iterative reconstruction of PET images using PC graphics hardware IEEE Nucl. Sci. Symp. Conf. Rec. 5 2787-90
4. Bell N and Garland M 2008 Efficient sparse matrix-vector multiplication on CUDA NVIDIA Technical Report NVR-2008-004 NVIDIA Corporation
5. Buluç A, Fineman J T, Frigo M, Gilbert J R and Leiserson C E 2009 Parallel sparse matrix-vector and matrix-transpose-vector multiplication using compressed sparse blocks SPAA'09: Proc. of the 21st Annu. Symp. on Parallelism in Algorithms and Architectures (ACM, NY, USA) pp 233-44
6. Cloquet C, Sureau F C, Defrise M, Simaeys G V, Trotta N and Goldman S 2010 Non-Gaussian space-variant resolution modelling for list-mode reconstruction Phys. Med. Biol. 55 5045
7. Herraiz J, Espaa S, Garcia S, Cabido R, Montemayor A, Desco M, Vaquero J and Udias J 2009 GPU acceleration of a fully 3D iterative reconstruction software for PET using CUDA IEEE Nucl. Sci. Symp. Conf. Rec. 4064-7
8. Hu J, Qi J, Huber J, Moses W and Huesman R 2005 MAP image reconstruction for arbitrary geometry PET systems with application to a prostate-specific scanner Int. Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine (Salt Lake City, UT) pp 416-20
9. Huesman R, Klein G, Moses W, Qi J, Reutter B and Virador P 2000 List-mode maximum-likelihood reconstruction applied to positron emission mammography (PEM) with irregular sampling IEEE Trans. Med. Imaging 19 532-7
10. Kadrmas D 2008 Rotate-and-slant projector for fast LOR-based fully-3-D iterative PET reconstruction IEEE Trans. Med. Imaging 27 1071-83
11. Karuta B and Lecomte R 1992 Effect of detector weighting functions on the point spread function of high-resolution PET tomographs: a simulation study IEEE Trans. Med. Imaging 11 379-85 (Pubitemid 23573819)
12. Kotasidis F, Matthews J, Angelis G, Noonan P, Jackson A, Price P, Lionheart W and Reader A 2011 Single scan parameterization of space-variant point spread functions in image space via a printed array: the impact for two PET/CT scanners Phys. Med. Biol. 56 2917-42
13. Lee D D and Seung H S 1999 Learning the parts of objects by non-negative matrix factorization Nature 401 788-91
14. Matej S and Lewitt R 1992 Image representation and tomographic reconstruction using spherically-symmetric volume elements IEEE Nucl. Sci. Symp. Conf. Rec. 2 1191-3
15. Moehrs S, Defrise M, Belcari N, Guerra A D, Bartoli A, Fabbri S and Zanetti G 2008 Multi-ray-based system matrix generation for 3D PET reconstruction Phys. Med. Biol. 53 6925
16. Mumcuoglu E, Leahy R, Cherry S and Hoffman E 1996 Accurate geometric and physical response modelling for statistical image reconstruction in high resolution PET IEEE Nucl. Sci. Symp. Conf. Rec. 3 1569-73
17. NVIDIA 2009a CUDA Data Parallel Primitives Library http://gpgpu.org/ developer/cudpp
18. NVIDIA 2009b CUDA Programming Guide 3.0 http://developer.download.nvidia. com/compute/cuda/3-0/toolkit/docs/NVIDIA-CUDA-ProgrammingGuide.pdf
19. NVIDIA 2009c CUDA sdk 2.3 http://developer.nvidia.com/object/cuda-2-3- downloads.html
20. NVIDIA 2009d CUDA toolkit 2.3 http://developer.nvidia.com/object/cuda-2- 3-downloads.html
21. Panin V, Kehren F, Michel C and Casey M 2006 Fully 3-D PET reconstruction with system matrix derived from point source measurements IEEE Trans. Med. Imaging 25 907-21 (Pubitemid 44088375)
22. Pratx G, Chinn G, Olcott P and Levin C 2009 Fast, accurate and shift-varying line projections for iterative reconstruction using the GPU IEEE Trans. Med. Imaging 28 435-45
23. Qi J, Leahy R, Cherry S, Chatziioannou A and Farquhar T 1998 High-resolution 3D Bayesian image reconstruction using the microPET small-animal scanner Phys. Med. Biol. 43 1001-13 (Pubitemid 28168428)
24. Rafecas M, Mosler B, Dietz M, Pogl M, Stamatakis A, McElroy D and Ziegler S 2004 Use of a Monte Carlo-based probability matrix for 3-D iterative reconstruction of MADPET II data IEEE Trans. Nucl. Sci. 51 2597-605
25. Rapisarda E, Bettinardi V, Thielemans K and Gilardi M C 2010 Image-based point spread function implementation in a fully 3D OSEM reconstruction algorithm for PET Phys. Med. Biol. 55 4131
26. Saad Y 1994 Sparskit: a basic tool kit for sparse matrix computations http://www-users.cs.umn.edu/saad/ software/SPARSKIT/sparskit.html
27. Scheins J, Herzog H and Shah N 2011 Fully-3D PET image reconstruction using scanner-independent, adaptive projection data and highly rotation-symmetric voxel assemblies IEEE Trans. Med. Imaging 30 879-92
28. Schmitt D, Karuta B, Carrier C and Lecomte R 1988 Fast point spread function computation from aperture functions in high-resolution positron emission tomography IEEE Trans. Med. Imaging 7 2-12 (Pubitemid 18638879)
29. Siddon R 1985 Fast calculation of the exact radiological path for a three-dimensional CT array Medi. Phys. 12 252-5 (Pubitemid 15089015)
30. Sureau F C, Reader A J, Comtat C, Leroy C, Ribeiro M-J, Buvat I and Trebossen R 2008 Impact of image-space resolution modeling for studies with the high-resolution research tomograph J. Nucl. Med. 49 1000-8 (Pubitemid 351793442)
31. Tai Y, Chatziioannou A, Yang Y, Silverman R, Meadors K, Siegel S, Newport D, Stickel J and Cherry S 2003 MicroPET II: design, development and initial performance of an improved microPET scanner for small-animal imaging Phys. Med. Biol. 48 1519-37 (Pubitemid 36713333)
32. Tohme M and Qi J 2009 Iterative image reconstruction for positron emission tomography based on a detector response function estimated from point source measurements Phys. Med. Biol. 54 3709-25
33. Zeng G and Gullberg G 2000 Unmatched projector/backprojector pairs in an iterative reconstruction algorithm IEEE Trans. Med. Imaging 19 548-55
34. Zhou J and Qi J 2009 Theoretical analysis and simulation study of a high-resolution zoom-in PET system Phys. Med. Biol. 54 5193