Efficient methodologies for system matrix modelling in iterative image reconstruction for rotating high-resolution PET

Ortuño Fisac, Juan Enrique and Kontaxakis Antoniadis, Georgios and Rubio Guivernau, José Luis and Guerra Gutiérrez, Pedro and Santos Lleo, Andres de (2010). Efficient methodologies for system matrix modelling in iterative image reconstruction for rotating high-resolution PET. "PHYSICS IN MEDICINE AND BIOLOGY", v. 55 (n. 7); pp. 1833-1861. ISSN 0031-9155. https://doi.org/10.1088/0031-9155/55/7/004.

Description

Title: Efficient methodologies for system matrix modelling in iterative image reconstruction for rotating high-resolution PET
Author/s:
  • Ortuño Fisac, Juan Enrique
  • Kontaxakis Antoniadis, Georgios
  • Rubio Guivernau, José Luis
  • Guerra Gutiérrez, Pedro
  • Santos Lleo, Andres de
Item Type: Article
Título de Revista/Publicación: PHYSICS IN MEDICINE AND BIOLOGY
Date: 2010
ISSN: 0031-9155
Volume: 55
Subjects:
Faculty: E.T.S.I. Telecomunicación (UPM)
Department: Ingeniería Electrónica
Creative Commons Licenses: Recognition - No derivative works - Non commercial

Full text

[img]
Preview
PDF - Requires a PDF viewer, such as GSview, Xpdf or Adobe Acrobat Reader
Download (1MB) | Preview

Abstract

A fully 3D iterative image reconstruction algorithm has been developed for high-resolution PET cameras composed of pixelated scintillator crystal arrays and rotating planar detectors, based on the ordered subsets approach. The associated system matrix is precalculated with Monte Carlo methods that incorporate physical effects not included in analytical models, such as positron range effects and interaction of the incident gammas with the scintillator material. Custom Monte Carlo methodologies have been developed and optimized for modelling of system matrices for fast iterative image reconstruction adapted to specific scanner geometries, without redundant calculations. According to the methodology proposed here, only one-eighth of the voxels within two central transaxial slices need to be modelled in detail. The rest of the system matrix elements can be obtained with the aid of axial symmetries and redundancies, as well as in-plane symmetries within transaxial slices. Sparse matrix techniques for the non-zero system matrix elements are employed, allowing for fast execution of the image reconstruction process. This 3D image reconstruction scheme has been compared in terms of image quality to a 2D fast implementation of the OSEM algorithm combined with Fourier rebinning approaches. This work confirms the superiority of fully 3D OSEM in terms of spatial resolution, contrast recovery and noise reduction as compared to conventional 2D approaches based on rebinning schemes. At the same time it demonstrates that fully 3D methodologies can be efficiently applied to the image reconstruction problem for high-resolution rotational PET cameras by applying accurate pre-calculated system models and taking advantage of the system's symmetries.

More information

Item ID: 13623
DC Identifier: http://oa.upm.es/13623/
OAI Identifier: oai:oa.upm.es:13623
DOI: 10.1088/0031-9155/55/7/004
Official URL: http://iopscience.iop.org/0031-9155/55/7/004/
Deposited by: Memoria Investigacion
Deposited on: 21 Nov 2012 11:38
Last Modified: 21 Apr 2016 12:57
  • Logo InvestigaM (UPM)
  • Logo GEOUP4
  • Logo Open Access
  • Open Access
  • Logo Sherpa/Romeo
    Check whether the anglo-saxon journal in which you have published an article allows you to also publish it under open access.
  • Logo Dulcinea
    Check whether the spanish journal in which you have published an article allows you to also publish it under open access.
  • Logo de Recolecta
  • Logo del Observatorio I+D+i UPM
  • Logo de OpenCourseWare UPM