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

Jian Zhou, Jinyi Qi

Research output: Contribution to journalArticle

47 Citations (Scopus)

Abstract

Statistically based iterative image reconstruction has been widely used in positron emission tomography (PET) imaging. The quality of reconstructed images depends on the accuracy of the system matrix that defines the mapping from the image space to the data space. However, an accurate system matrix is often associated with high computation cost and huge storage requirement. In this paper, we present a method to address this problem using sparse matrix factorization and graphics processor unit (GPU) acceleration. We factor the accurate system matrix into three highly sparse matrices: a sinogram blurring matrix, a geometric projection matrix and an image blurring matrix. The geometrical projection matrix is precomputed based on a simple line integral model, while the sinogram and image blurring matrices are estimated from point-source measurements. The resulting factored system matrix has far less nonzero elements than the original system matrix, which substantially reduces the storage and computation cost. The smaller matrix size also allows an efficient implementation of the forward and backward projectors on a GPU, which often has a limited memory space. Our experimental studies show that the proposed method can dramatically reduce the computation cost of high-resolution iterative image reconstruction, while achieving better performance than existing factorization methods.

Original languageEnglish (US)
Pages (from-to)6739-6757
Number of pages19
JournalPhysics in Medicine and Biology
Volume56
Issue number20
DOIs
StatePublished - Oct 21 2011

Fingerprint

Computer-Assisted Image Processing
Positron-Emission Tomography
Costs and Cost Analysis

ASJC Scopus subject areas

  • Radiology Nuclear Medicine and imaging
  • Radiological and Ultrasound Technology

Cite this

Fast and efficient fully 3D PET image reconstruction using sparse system matrix factorization with GPU acceleration. / Zhou, Jian; Qi, Jinyi.

In: Physics in Medicine and Biology, Vol. 56, No. 20, 21.10.2011, p. 6739-6757.

Research output: Contribution to journalArticle

@article{de55341e43d249f5ab1dc93f493138d9,
title = "Fast and efficient fully 3D PET image reconstruction using sparse system matrix factorization with GPU acceleration",
abstract = "Statistically based iterative image reconstruction has been widely used in positron emission tomography (PET) imaging. The quality of reconstructed images depends on the accuracy of the system matrix that defines the mapping from the image space to the data space. However, an accurate system matrix is often associated with high computation cost and huge storage requirement. In this paper, we present a method to address this problem using sparse matrix factorization and graphics processor unit (GPU) acceleration. We factor the accurate system matrix into three highly sparse matrices: a sinogram blurring matrix, a geometric projection matrix and an image blurring matrix. The geometrical projection matrix is precomputed based on a simple line integral model, while the sinogram and image blurring matrices are estimated from point-source measurements. The resulting factored system matrix has far less nonzero elements than the original system matrix, which substantially reduces the storage and computation cost. The smaller matrix size also allows an efficient implementation of the forward and backward projectors on a GPU, which often has a limited memory space. Our experimental studies show that the proposed method can dramatically reduce the computation cost of high-resolution iterative image reconstruction, while achieving better performance than existing factorization methods.",
author = "Jian Zhou and Jinyi Qi",
year = "2011",
month = "10",
day = "21",
doi = "10.1088/0031-9155/56/20/015",
language = "English (US)",
volume = "56",
pages = "6739--6757",
journal = "Physics in Medicine and Biology",
issn = "0031-9155",
publisher = "IOP Publishing Ltd.",
number = "20",

}

TY - JOUR

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

AU - Zhou, Jian

AU - Qi, Jinyi

PY - 2011/10/21

Y1 - 2011/10/21

N2 - Statistically based iterative image reconstruction has been widely used in positron emission tomography (PET) imaging. The quality of reconstructed images depends on the accuracy of the system matrix that defines the mapping from the image space to the data space. However, an accurate system matrix is often associated with high computation cost and huge storage requirement. In this paper, we present a method to address this problem using sparse matrix factorization and graphics processor unit (GPU) acceleration. We factor the accurate system matrix into three highly sparse matrices: a sinogram blurring matrix, a geometric projection matrix and an image blurring matrix. The geometrical projection matrix is precomputed based on a simple line integral model, while the sinogram and image blurring matrices are estimated from point-source measurements. The resulting factored system matrix has far less nonzero elements than the original system matrix, which substantially reduces the storage and computation cost. The smaller matrix size also allows an efficient implementation of the forward and backward projectors on a GPU, which often has a limited memory space. Our experimental studies show that the proposed method can dramatically reduce the computation cost of high-resolution iterative image reconstruction, while achieving better performance than existing factorization methods.

AB - Statistically based iterative image reconstruction has been widely used in positron emission tomography (PET) imaging. The quality of reconstructed images depends on the accuracy of the system matrix that defines the mapping from the image space to the data space. However, an accurate system matrix is often associated with high computation cost and huge storage requirement. In this paper, we present a method to address this problem using sparse matrix factorization and graphics processor unit (GPU) acceleration. We factor the accurate system matrix into three highly sparse matrices: a sinogram blurring matrix, a geometric projection matrix and an image blurring matrix. The geometrical projection matrix is precomputed based on a simple line integral model, while the sinogram and image blurring matrices are estimated from point-source measurements. The resulting factored system matrix has far less nonzero elements than the original system matrix, which substantially reduces the storage and computation cost. The smaller matrix size also allows an efficient implementation of the forward and backward projectors on a GPU, which often has a limited memory space. Our experimental studies show that the proposed method can dramatically reduce the computation cost of high-resolution iterative image reconstruction, while achieving better performance than existing factorization methods.

UR - http://www.scopus.com/inward/record.url?scp=80053616587&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=80053616587&partnerID=8YFLogxK

U2 - 10.1088/0031-9155/56/20/015

DO - 10.1088/0031-9155/56/20/015

M3 - Article

VL - 56

SP - 6739

EP - 6757

JO - Physics in Medicine and Biology

JF - Physics in Medicine and Biology

SN - 0031-9155

IS - 20

ER -