Iterative reconstruction of Fourier-rebinned PET data using sinogram blurring function estimated from point source scans

Michel S. Tohme, Jinyi Qi

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

Purpose: The accuracy of the system model that governs the transformation from the image space to the projection space in positron emission tomography (PET) greatly affects the quality of reconstructed images. For efficient computation in iterative reconstructions, the system model in PET can be factored into a product of geometric projection and sinogram blurring function. To further speed up reconstruction, fully 3D PET data can be rebinned into a stack of 2D sinograms and then be reconstructed using 2D iterative algorithms. The purpose of this work is to develop a method to estimate the sinogram blurring function to be used in reconstruction of Fourier-rebinned data. Methods: In a previous work, the authors developed an approach to estimating the sinogram blurring function of nonrebinned PET data from experimental scans of point sources. In this study, the authors extend this method to the estimation of sinogram blurring function for Fourier-rebinned PET data. A point source was scanned at a set of sampled positions in the microPET II scanner. The sinogram blurring function is considered to be separable between the transaxial and axial directions. A radially and angularly variant 2D blurring function is estimated from Fourier-rebinned point source scans to model the transaxial blurring with consideration of the detector block structure of the scanner; a space-variant 1D blurring kernel along the axial direction is estimated separately to model the correlation between neighboring planes due to detector intrinsic blurring and Fourier rebinning. The estimated sinogram blurring function is incorporated in a 2D maximum a posteriori (MAP) reconstruction algorithm for image reconstruction. Results: Physical phantom experiments were performed on the microPET II scanner to validate the proposed method. The authors compared the proposed method to 2D MAP reconstruction without sinogram blurring model and 2D MAP reconstruction with a Monte Carlo based blurring model. The results show that the proposed method produces images with improved contrast and spatial resolution. The reconstruction time is unaffected by the new method since the blurring component takes a relatively negligible part of the overall reconstruction time. Conclusions: The proposed method can estimate sinogram blurring matrix for Fourier-rebinned PET data and can be used to improve contrast and spatial resolution of reconstructed images. The method can be applied to other human and animal scanners.

Original languageEnglish (US)
Pages (from-to)5530-5540
Number of pages11
JournalMedical Physics
Volume37
Issue number10
DOIs
StatePublished - Oct 2010

Fingerprint

Positron-Emission Tomography
Computer-Assisted Image Processing

Keywords

  • axial blurring kernel
  • detector response
  • iterative image reconstruction
  • positron emission tomography
  • sinogram blurring kernel
  • system modeling

ASJC Scopus subject areas

  • Biophysics
  • Radiology Nuclear Medicine and imaging

Cite this

Iterative reconstruction of Fourier-rebinned PET data using sinogram blurring function estimated from point source scans. / Tohme, Michel S.; Qi, Jinyi.

In: Medical Physics, Vol. 37, No. 10, 10.2010, p. 5530-5540.

Research output: Contribution to journalArticle

@article{c6b33067af2b47eb8ada42942cf1eb53,
title = "Iterative reconstruction of Fourier-rebinned PET data using sinogram blurring function estimated from point source scans",
abstract = "Purpose: The accuracy of the system model that governs the transformation from the image space to the projection space in positron emission tomography (PET) greatly affects the quality of reconstructed images. For efficient computation in iterative reconstructions, the system model in PET can be factored into a product of geometric projection and sinogram blurring function. To further speed up reconstruction, fully 3D PET data can be rebinned into a stack of 2D sinograms and then be reconstructed using 2D iterative algorithms. The purpose of this work is to develop a method to estimate the sinogram blurring function to be used in reconstruction of Fourier-rebinned data. Methods: In a previous work, the authors developed an approach to estimating the sinogram blurring function of nonrebinned PET data from experimental scans of point sources. In this study, the authors extend this method to the estimation of sinogram blurring function for Fourier-rebinned PET data. A point source was scanned at a set of sampled positions in the microPET II scanner. The sinogram blurring function is considered to be separable between the transaxial and axial directions. A radially and angularly variant 2D blurring function is estimated from Fourier-rebinned point source scans to model the transaxial blurring with consideration of the detector block structure of the scanner; a space-variant 1D blurring kernel along the axial direction is estimated separately to model the correlation between neighboring planes due to detector intrinsic blurring and Fourier rebinning. The estimated sinogram blurring function is incorporated in a 2D maximum a posteriori (MAP) reconstruction algorithm for image reconstruction. Results: Physical phantom experiments were performed on the microPET II scanner to validate the proposed method. The authors compared the proposed method to 2D MAP reconstruction without sinogram blurring model and 2D MAP reconstruction with a Monte Carlo based blurring model. The results show that the proposed method produces images with improved contrast and spatial resolution. The reconstruction time is unaffected by the new method since the blurring component takes a relatively negligible part of the overall reconstruction time. Conclusions: The proposed method can estimate sinogram blurring matrix for Fourier-rebinned PET data and can be used to improve contrast and spatial resolution of reconstructed images. The method can be applied to other human and animal scanners.",
keywords = "axial blurring kernel, detector response, iterative image reconstruction, positron emission tomography, sinogram blurring kernel, system modeling",
author = "Tohme, {Michel S.} and Jinyi Qi",
year = "2010",
month = "10",
doi = "10.1118/1.3490711",
language = "English (US)",
volume = "37",
pages = "5530--5540",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "10",

}

TY - JOUR

T1 - Iterative reconstruction of Fourier-rebinned PET data using sinogram blurring function estimated from point source scans

AU - Tohme, Michel S.

AU - Qi, Jinyi

PY - 2010/10

Y1 - 2010/10

N2 - Purpose: The accuracy of the system model that governs the transformation from the image space to the projection space in positron emission tomography (PET) greatly affects the quality of reconstructed images. For efficient computation in iterative reconstructions, the system model in PET can be factored into a product of geometric projection and sinogram blurring function. To further speed up reconstruction, fully 3D PET data can be rebinned into a stack of 2D sinograms and then be reconstructed using 2D iterative algorithms. The purpose of this work is to develop a method to estimate the sinogram blurring function to be used in reconstruction of Fourier-rebinned data. Methods: In a previous work, the authors developed an approach to estimating the sinogram blurring function of nonrebinned PET data from experimental scans of point sources. In this study, the authors extend this method to the estimation of sinogram blurring function for Fourier-rebinned PET data. A point source was scanned at a set of sampled positions in the microPET II scanner. The sinogram blurring function is considered to be separable between the transaxial and axial directions. A radially and angularly variant 2D blurring function is estimated from Fourier-rebinned point source scans to model the transaxial blurring with consideration of the detector block structure of the scanner; a space-variant 1D blurring kernel along the axial direction is estimated separately to model the correlation between neighboring planes due to detector intrinsic blurring and Fourier rebinning. The estimated sinogram blurring function is incorporated in a 2D maximum a posteriori (MAP) reconstruction algorithm for image reconstruction. Results: Physical phantom experiments were performed on the microPET II scanner to validate the proposed method. The authors compared the proposed method to 2D MAP reconstruction without sinogram blurring model and 2D MAP reconstruction with a Monte Carlo based blurring model. The results show that the proposed method produces images with improved contrast and spatial resolution. The reconstruction time is unaffected by the new method since the blurring component takes a relatively negligible part of the overall reconstruction time. Conclusions: The proposed method can estimate sinogram blurring matrix for Fourier-rebinned PET data and can be used to improve contrast and spatial resolution of reconstructed images. The method can be applied to other human and animal scanners.

AB - Purpose: The accuracy of the system model that governs the transformation from the image space to the projection space in positron emission tomography (PET) greatly affects the quality of reconstructed images. For efficient computation in iterative reconstructions, the system model in PET can be factored into a product of geometric projection and sinogram blurring function. To further speed up reconstruction, fully 3D PET data can be rebinned into a stack of 2D sinograms and then be reconstructed using 2D iterative algorithms. The purpose of this work is to develop a method to estimate the sinogram blurring function to be used in reconstruction of Fourier-rebinned data. Methods: In a previous work, the authors developed an approach to estimating the sinogram blurring function of nonrebinned PET data from experimental scans of point sources. In this study, the authors extend this method to the estimation of sinogram blurring function for Fourier-rebinned PET data. A point source was scanned at a set of sampled positions in the microPET II scanner. The sinogram blurring function is considered to be separable between the transaxial and axial directions. A radially and angularly variant 2D blurring function is estimated from Fourier-rebinned point source scans to model the transaxial blurring with consideration of the detector block structure of the scanner; a space-variant 1D blurring kernel along the axial direction is estimated separately to model the correlation between neighboring planes due to detector intrinsic blurring and Fourier rebinning. The estimated sinogram blurring function is incorporated in a 2D maximum a posteriori (MAP) reconstruction algorithm for image reconstruction. Results: Physical phantom experiments were performed on the microPET II scanner to validate the proposed method. The authors compared the proposed method to 2D MAP reconstruction without sinogram blurring model and 2D MAP reconstruction with a Monte Carlo based blurring model. The results show that the proposed method produces images with improved contrast and spatial resolution. The reconstruction time is unaffected by the new method since the blurring component takes a relatively negligible part of the overall reconstruction time. Conclusions: The proposed method can estimate sinogram blurring matrix for Fourier-rebinned PET data and can be used to improve contrast and spatial resolution of reconstructed images. The method can be applied to other human and animal scanners.

KW - axial blurring kernel

KW - detector response

KW - iterative image reconstruction

KW - positron emission tomography

KW - sinogram blurring kernel

KW - system modeling

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

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

U2 - 10.1118/1.3490711

DO - 10.1118/1.3490711

M3 - Article

VL - 37

SP - 5530

EP - 5540

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 10

ER -