High-resolution 3D bayesian image reconstruction using the microPET small-animal scanner

Jinyi Qi, Richard M. Leahy, Simon R Cherry, Arion Chatziioannou, Thomas H. Farquhar

Research output: Contribution to journalArticle

498 Citations (Scopus)

Abstract

A Bayesian method is described for reconstruction of high-resolution 3D images from the microPET small-animal scanner. Resolution recovery is achieved by explicitly modelling the depth dependent geometric sensitivity for each voxel in combination with an accurate detector response model that includes factors due to photon pair non-collinearity and inter-crystal scatter and penetration. To reduce storage and computational costs we use a factored matrix in which the detector response is modelled using a sinogram blurring kernel. Maximum a posteriori (MAP) images are reconstructed using this model in combination with a Poisson likelihood function and a Gibbs prior on the image. Reconstructions obtained from point source data using the accurate system model demonstrate a potential for near-isotropic FWHM resolution of approximately 1.2 mm at the center of the field of view compared with approximately 2 rnm when using an analytic 3D reprojection (3DRP) method with a ramp filter. These results also show the ability of the accurate system model to compensate for resolution loss due to crystal penetration producing nearly constant radial FWHM resolution of 1 mm out to a 4 mm radius. Studies with a point source in a uniform cylinder indicate that as the resolution of the image is reduced to control noise propagation the resolution obtained using the accurate system model is superior to that obtained using 3DRP at matched background noise levels. Additional studies using pie phantoms with hot and cold cylinders of diameter 1-2.5 mm and FDG animal studies appear to confirm this observation.

Original languageEnglish (US)
Pages (from-to)1001-1013
Number of pages13
JournalPhysics in Medicine and Biology
Volume43
Issue number4
DOIs
StatePublished - Apr 1998

Fingerprint

Computer-Assisted Image Processing
image reconstruction
Image reconstruction
scanners
animals
Noise
Animals
Likelihood Functions
Architectural Accessibility
Bayes Theorem
Information Storage and Retrieval
high resolution
Photons
Full width at half maximum
Costs and Cost Analysis
point sources
Detectors
penetration
Crystals
Acoustic variables control

ASJC Scopus subject areas

  • Biomedical Engineering
  • Physics and Astronomy (miscellaneous)
  • Radiology Nuclear Medicine and imaging
  • Radiological and Ultrasound Technology

Cite this

High-resolution 3D bayesian image reconstruction using the microPET small-animal scanner. / Qi, Jinyi; Leahy, Richard M.; Cherry, Simon R; Chatziioannou, Arion; Farquhar, Thomas H.

In: Physics in Medicine and Biology, Vol. 43, No. 4, 04.1998, p. 1001-1013.

Research output: Contribution to journalArticle

Qi, Jinyi ; Leahy, Richard M. ; Cherry, Simon R ; Chatziioannou, Arion ; Farquhar, Thomas H. / High-resolution 3D bayesian image reconstruction using the microPET small-animal scanner. In: Physics in Medicine and Biology. 1998 ; Vol. 43, No. 4. pp. 1001-1013.
@article{4ef2114507f0464db7d3abfdb5c93222,
title = "High-resolution 3D bayesian image reconstruction using the microPET small-animal scanner",
abstract = "A Bayesian method is described for reconstruction of high-resolution 3D images from the microPET small-animal scanner. Resolution recovery is achieved by explicitly modelling the depth dependent geometric sensitivity for each voxel in combination with an accurate detector response model that includes factors due to photon pair non-collinearity and inter-crystal scatter and penetration. To reduce storage and computational costs we use a factored matrix in which the detector response is modelled using a sinogram blurring kernel. Maximum a posteriori (MAP) images are reconstructed using this model in combination with a Poisson likelihood function and a Gibbs prior on the image. Reconstructions obtained from point source data using the accurate system model demonstrate a potential for near-isotropic FWHM resolution of approximately 1.2 mm at the center of the field of view compared with approximately 2 rnm when using an analytic 3D reprojection (3DRP) method with a ramp filter. These results also show the ability of the accurate system model to compensate for resolution loss due to crystal penetration producing nearly constant radial FWHM resolution of 1 mm out to a 4 mm radius. Studies with a point source in a uniform cylinder indicate that as the resolution of the image is reduced to control noise propagation the resolution obtained using the accurate system model is superior to that obtained using 3DRP at matched background noise levels. Additional studies using pie phantoms with hot and cold cylinders of diameter 1-2.5 mm and FDG animal studies appear to confirm this observation.",
author = "Jinyi Qi and Leahy, {Richard M.} and Cherry, {Simon R} and Arion Chatziioannou and Farquhar, {Thomas H.}",
year = "1998",
month = "4",
doi = "10.1088/0031-9155/43/4/027",
language = "English (US)",
volume = "43",
pages = "1001--1013",
journal = "Physics in Medicine and Biology",
issn = "0031-9155",
publisher = "IOP Publishing Ltd.",
number = "4",

}

TY - JOUR

T1 - High-resolution 3D bayesian image reconstruction using the microPET small-animal scanner

AU - Qi, Jinyi

AU - Leahy, Richard M.

AU - Cherry, Simon R

AU - Chatziioannou, Arion

AU - Farquhar, Thomas H.

PY - 1998/4

Y1 - 1998/4

N2 - A Bayesian method is described for reconstruction of high-resolution 3D images from the microPET small-animal scanner. Resolution recovery is achieved by explicitly modelling the depth dependent geometric sensitivity for each voxel in combination with an accurate detector response model that includes factors due to photon pair non-collinearity and inter-crystal scatter and penetration. To reduce storage and computational costs we use a factored matrix in which the detector response is modelled using a sinogram blurring kernel. Maximum a posteriori (MAP) images are reconstructed using this model in combination with a Poisson likelihood function and a Gibbs prior on the image. Reconstructions obtained from point source data using the accurate system model demonstrate a potential for near-isotropic FWHM resolution of approximately 1.2 mm at the center of the field of view compared with approximately 2 rnm when using an analytic 3D reprojection (3DRP) method with a ramp filter. These results also show the ability of the accurate system model to compensate for resolution loss due to crystal penetration producing nearly constant radial FWHM resolution of 1 mm out to a 4 mm radius. Studies with a point source in a uniform cylinder indicate that as the resolution of the image is reduced to control noise propagation the resolution obtained using the accurate system model is superior to that obtained using 3DRP at matched background noise levels. Additional studies using pie phantoms with hot and cold cylinders of diameter 1-2.5 mm and FDG animal studies appear to confirm this observation.

AB - A Bayesian method is described for reconstruction of high-resolution 3D images from the microPET small-animal scanner. Resolution recovery is achieved by explicitly modelling the depth dependent geometric sensitivity for each voxel in combination with an accurate detector response model that includes factors due to photon pair non-collinearity and inter-crystal scatter and penetration. To reduce storage and computational costs we use a factored matrix in which the detector response is modelled using a sinogram blurring kernel. Maximum a posteriori (MAP) images are reconstructed using this model in combination with a Poisson likelihood function and a Gibbs prior on the image. Reconstructions obtained from point source data using the accurate system model demonstrate a potential for near-isotropic FWHM resolution of approximately 1.2 mm at the center of the field of view compared with approximately 2 rnm when using an analytic 3D reprojection (3DRP) method with a ramp filter. These results also show the ability of the accurate system model to compensate for resolution loss due to crystal penetration producing nearly constant radial FWHM resolution of 1 mm out to a 4 mm radius. Studies with a point source in a uniform cylinder indicate that as the resolution of the image is reduced to control noise propagation the resolution obtained using the accurate system model is superior to that obtained using 3DRP at matched background noise levels. Additional studies using pie phantoms with hot and cold cylinders of diameter 1-2.5 mm and FDG animal studies appear to confirm this observation.

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

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

U2 - 10.1088/0031-9155/43/4/027

DO - 10.1088/0031-9155/43/4/027

M3 - Article

C2 - 9572523

AN - SCOPUS:0031981253

VL - 43

SP - 1001

EP - 1013

JO - Physics in Medicine and Biology

JF - Physics in Medicine and Biology

SN - 0031-9155

IS - 4

ER -