Fully 3D bayesian image reconstruction for the ECAT EXACT HR+ l

Jinyi Qi, Richard M. Leahy, Chinghan Hsu, Thomas H. Farquhar, Simon R Cherry

Research output: Contribution to journalArticle

113 Citations (Scopus)

Abstract

A fully 3D Bayesian method is described for high resolution reconstruction of images from the Siemens/CTI ECAT EXACT HR+ whole body positron emission tomography (PET) scanner. To maximize resolution recovery from the system we model depth dependent geometric efficiency, intrinsic detector efficiency, photon pair non-colinearity, crystal penetration and inter-crystal scatter. We also explicitly model the effects of axial rebinning and angular mashing on the detection probability or system matrix. By fully exploiting sinogram symmetries and using a factored system matrix and automated indexing schemes, we are able to achieve substantial savings in both the storage size and time required to compute forward and backward projections. Reconstruction times are further reduced using multi-threaded programming on a four processor Unix server. Bayesian reconstructions are computed using a Huber prior and a shifted-Poisson likelihood model that accounts for the effects of randoms subtraction and scatter. Reconstructions of phantom data show that the 3D Bayesian method can achieve improved FWHM resolution and contrast recovery ratios at matched background noise levels compared to both the 3D reprojection method and an OSEM method based on the shifted-Poisson model.

Original languageEnglish (US)
Pages (from-to)1096-1103
Number of pages8
JournalIEEE Transactions on Nuclear Science
Volume45
Issue number3 PART 2
StatePublished - 1998

Fingerprint

image reconstruction
Image reconstruction
recovery
Recovery
Crystals
Positron emission tomography
hubs
background noise
matrices
programming
Full width at half maximum
Computer programming
subtraction
scanners
crystals
Program processors
central processing units
positrons
Servers
penetration

ASJC Scopus subject areas

  • Electrical and Electronic Engineering
  • Nuclear Energy and Engineering

Cite this

Qi, J., Leahy, R. M., Hsu, C., Farquhar, T. H., & Cherry, S. R. (1998). Fully 3D bayesian image reconstruction for the ECAT EXACT HR+ l. IEEE Transactions on Nuclear Science, 45(3 PART 2), 1096-1103.

Fully 3D bayesian image reconstruction for the ECAT EXACT HR+ l. / Qi, Jinyi; Leahy, Richard M.; Hsu, Chinghan; Farquhar, Thomas H.; Cherry, Simon R.

In: IEEE Transactions on Nuclear Science, Vol. 45, No. 3 PART 2, 1998, p. 1096-1103.

Research output: Contribution to journalArticle

Qi, J, Leahy, RM, Hsu, C, Farquhar, TH & Cherry, SR 1998, 'Fully 3D bayesian image reconstruction for the ECAT EXACT HR+ l', IEEE Transactions on Nuclear Science, vol. 45, no. 3 PART 2, pp. 1096-1103.
Qi, Jinyi ; Leahy, Richard M. ; Hsu, Chinghan ; Farquhar, Thomas H. ; Cherry, Simon R. / Fully 3D bayesian image reconstruction for the ECAT EXACT HR+ l. In: IEEE Transactions on Nuclear Science. 1998 ; Vol. 45, No. 3 PART 2. pp. 1096-1103.
@article{23b5ea2919824c1cb9c8807dac755ed9,
title = "Fully 3D bayesian image reconstruction for the ECAT EXACT HR+ l",
abstract = "A fully 3D Bayesian method is described for high resolution reconstruction of images from the Siemens/CTI ECAT EXACT HR+ whole body positron emission tomography (PET) scanner. To maximize resolution recovery from the system we model depth dependent geometric efficiency, intrinsic detector efficiency, photon pair non-colinearity, crystal penetration and inter-crystal scatter. We also explicitly model the effects of axial rebinning and angular mashing on the detection probability or system matrix. By fully exploiting sinogram symmetries and using a factored system matrix and automated indexing schemes, we are able to achieve substantial savings in both the storage size and time required to compute forward and backward projections. Reconstruction times are further reduced using multi-threaded programming on a four processor Unix server. Bayesian reconstructions are computed using a Huber prior and a shifted-Poisson likelihood model that accounts for the effects of randoms subtraction and scatter. Reconstructions of phantom data show that the 3D Bayesian method can achieve improved FWHM resolution and contrast recovery ratios at matched background noise levels compared to both the 3D reprojection method and an OSEM method based on the shifted-Poisson model.",
author = "Jinyi Qi and Leahy, {Richard M.} and Chinghan Hsu and Farquhar, {Thomas H.} and Cherry, {Simon R}",
year = "1998",
language = "English (US)",
volume = "45",
pages = "1096--1103",
journal = "IEEE Transactions on Nuclear Science",
issn = "0018-9499",
publisher = "Institute of Electrical and Electronics Engineers Inc.",
number = "3 PART 2",

}

TY - JOUR

T1 - Fully 3D bayesian image reconstruction for the ECAT EXACT HR+ l

AU - Qi, Jinyi

AU - Leahy, Richard M.

AU - Hsu, Chinghan

AU - Farquhar, Thomas H.

AU - Cherry, Simon R

PY - 1998

Y1 - 1998

N2 - A fully 3D Bayesian method is described for high resolution reconstruction of images from the Siemens/CTI ECAT EXACT HR+ whole body positron emission tomography (PET) scanner. To maximize resolution recovery from the system we model depth dependent geometric efficiency, intrinsic detector efficiency, photon pair non-colinearity, crystal penetration and inter-crystal scatter. We also explicitly model the effects of axial rebinning and angular mashing on the detection probability or system matrix. By fully exploiting sinogram symmetries and using a factored system matrix and automated indexing schemes, we are able to achieve substantial savings in both the storage size and time required to compute forward and backward projections. Reconstruction times are further reduced using multi-threaded programming on a four processor Unix server. Bayesian reconstructions are computed using a Huber prior and a shifted-Poisson likelihood model that accounts for the effects of randoms subtraction and scatter. Reconstructions of phantom data show that the 3D Bayesian method can achieve improved FWHM resolution and contrast recovery ratios at matched background noise levels compared to both the 3D reprojection method and an OSEM method based on the shifted-Poisson model.

AB - A fully 3D Bayesian method is described for high resolution reconstruction of images from the Siemens/CTI ECAT EXACT HR+ whole body positron emission tomography (PET) scanner. To maximize resolution recovery from the system we model depth dependent geometric efficiency, intrinsic detector efficiency, photon pair non-colinearity, crystal penetration and inter-crystal scatter. We also explicitly model the effects of axial rebinning and angular mashing on the detection probability or system matrix. By fully exploiting sinogram symmetries and using a factored system matrix and automated indexing schemes, we are able to achieve substantial savings in both the storage size and time required to compute forward and backward projections. Reconstruction times are further reduced using multi-threaded programming on a four processor Unix server. Bayesian reconstructions are computed using a Huber prior and a shifted-Poisson likelihood model that accounts for the effects of randoms subtraction and scatter. Reconstructions of phantom data show that the 3D Bayesian method can achieve improved FWHM resolution and contrast recovery ratios at matched background noise levels compared to both the 3D reprojection method and an OSEM method based on the shifted-Poisson model.

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

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

M3 - Article

VL - 45

SP - 1096

EP - 1103

JO - IEEE Transactions on Nuclear Science

JF - IEEE Transactions on Nuclear Science

SN - 0018-9499

IS - 3 PART 2

ER -