System matrix for OSEM SPECT with attenuation compensation in mesh domain

Levon Vogelsang, Andrzej Krol, David H. Feiglin, Edward Lipson

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

The purpose of this study was to develop and implement an accurate and computationally efficient method for determination of the mesh-domain system matrix including attenuation compensation for Ordered Subsets Expectation Maximization (OSEM) Single Photon Emission Computed Tomography (SPECT). The mesh-domain system matrix elements were estimated by first partitioning the object domain into strips parallel to detector face and with width not exceeding the size of a detector unit. This was followed by approximating the integration over the strip/mesh-element union. This approximation is product of: (i) strip width, (ii) intersection length of a ray central to strip with a mesh element, and (iii) the response and expansion function evaluated at midpoint of the intersection length. Reconstruction was performed using OSEM without regularization and with exact knowledge of the attenuation map. The method was evaluated using synthetic SPECT data generated using SIMIND Monte Carlo simulation software. Comparative quantitative and qualitative analysis included: bias, variance, standard deviation and line-profiles within three different regions of interest. We found that no more than two divisions per detector bin were needed for good quality reconstructed images when using a high resolution mesh.

Original languageEnglish (US)
Title of host publicationProgress in Biomedical Optics and Imaging - Proceedings of SPIE
Volume7622
EditionPART 3
DOIs
StatePublished - 2010
EventMedical Imaging 2010: Physics of Medical Imaging - San Diego, CA, United States
Duration: Feb 15 2010Feb 18 2010

Other

OtherMedical Imaging 2010: Physics of Medical Imaging
CountryUnited States
CitySan Diego, CA
Period2/15/102/18/10

Fingerprint

Single photon emission computed tomography
Single-Photon Emission-Computed Tomography
Set theory
set theory
mesh
tomography
attenuation
Detectors
strip
photons
Response Elements
matrices
Software
intersections
Bins
detectors
Image quality
unions
qualitative analysis
quantitative analysis

Keywords

  • Attenuation compensation
  • Mesh domain
  • OSEM
  • Reconstruction
  • SPECT
  • System matrix

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics
  • Electronic, Optical and Magnetic Materials
  • Biomaterials
  • Radiology Nuclear Medicine and imaging

Cite this

Vogelsang, L., Krol, A., Feiglin, D. H., & Lipson, E. (2010). System matrix for OSEM SPECT with attenuation compensation in mesh domain. In Progress in Biomedical Optics and Imaging - Proceedings of SPIE (PART 3 ed., Vol. 7622). [76225B] https://doi.org/10.1117/12.844427

System matrix for OSEM SPECT with attenuation compensation in mesh domain. / Vogelsang, Levon; Krol, Andrzej; Feiglin, David H.; Lipson, Edward.

Progress in Biomedical Optics and Imaging - Proceedings of SPIE. Vol. 7622 PART 3. ed. 2010. 76225B.

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Vogelsang, L, Krol, A, Feiglin, DH & Lipson, E 2010, System matrix for OSEM SPECT with attenuation compensation in mesh domain. in Progress in Biomedical Optics and Imaging - Proceedings of SPIE. PART 3 edn, vol. 7622, 76225B, Medical Imaging 2010: Physics of Medical Imaging, San Diego, CA, United States, 2/15/10. https://doi.org/10.1117/12.844427
Vogelsang L, Krol A, Feiglin DH, Lipson E. System matrix for OSEM SPECT with attenuation compensation in mesh domain. In Progress in Biomedical Optics and Imaging - Proceedings of SPIE. PART 3 ed. Vol. 7622. 2010. 76225B https://doi.org/10.1117/12.844427
Vogelsang, Levon ; Krol, Andrzej ; Feiglin, David H. ; Lipson, Edward. / System matrix for OSEM SPECT with attenuation compensation in mesh domain. Progress in Biomedical Optics and Imaging - Proceedings of SPIE. Vol. 7622 PART 3. ed. 2010.
@inproceedings{267b9e28559e466290d9740b100dd8bf,
title = "System matrix for OSEM SPECT with attenuation compensation in mesh domain",
abstract = "The purpose of this study was to develop and implement an accurate and computationally efficient method for determination of the mesh-domain system matrix including attenuation compensation for Ordered Subsets Expectation Maximization (OSEM) Single Photon Emission Computed Tomography (SPECT). The mesh-domain system matrix elements were estimated by first partitioning the object domain into strips parallel to detector face and with width not exceeding the size of a detector unit. This was followed by approximating the integration over the strip/mesh-element union. This approximation is product of: (i) strip width, (ii) intersection length of a ray central to strip with a mesh element, and (iii) the response and expansion function evaluated at midpoint of the intersection length. Reconstruction was performed using OSEM without regularization and with exact knowledge of the attenuation map. The method was evaluated using synthetic SPECT data generated using SIMIND Monte Carlo simulation software. Comparative quantitative and qualitative analysis included: bias, variance, standard deviation and line-profiles within three different regions of interest. We found that no more than two divisions per detector bin were needed for good quality reconstructed images when using a high resolution mesh.",
keywords = "Attenuation compensation, Mesh domain, OSEM, Reconstruction, SPECT, System matrix",
author = "Levon Vogelsang and Andrzej Krol and Feiglin, {David H.} and Edward Lipson",
year = "2010",
doi = "10.1117/12.844427",
language = "English (US)",
isbn = "9780819480231",
volume = "7622",
booktitle = "Progress in Biomedical Optics and Imaging - Proceedings of SPIE",
edition = "PART 3",

}

TY - GEN

T1 - System matrix for OSEM SPECT with attenuation compensation in mesh domain

AU - Vogelsang, Levon

AU - Krol, Andrzej

AU - Feiglin, David H.

AU - Lipson, Edward

PY - 2010

Y1 - 2010

N2 - The purpose of this study was to develop and implement an accurate and computationally efficient method for determination of the mesh-domain system matrix including attenuation compensation for Ordered Subsets Expectation Maximization (OSEM) Single Photon Emission Computed Tomography (SPECT). The mesh-domain system matrix elements were estimated by first partitioning the object domain into strips parallel to detector face and with width not exceeding the size of a detector unit. This was followed by approximating the integration over the strip/mesh-element union. This approximation is product of: (i) strip width, (ii) intersection length of a ray central to strip with a mesh element, and (iii) the response and expansion function evaluated at midpoint of the intersection length. Reconstruction was performed using OSEM without regularization and with exact knowledge of the attenuation map. The method was evaluated using synthetic SPECT data generated using SIMIND Monte Carlo simulation software. Comparative quantitative and qualitative analysis included: bias, variance, standard deviation and line-profiles within three different regions of interest. We found that no more than two divisions per detector bin were needed for good quality reconstructed images when using a high resolution mesh.

AB - The purpose of this study was to develop and implement an accurate and computationally efficient method for determination of the mesh-domain system matrix including attenuation compensation for Ordered Subsets Expectation Maximization (OSEM) Single Photon Emission Computed Tomography (SPECT). The mesh-domain system matrix elements were estimated by first partitioning the object domain into strips parallel to detector face and with width not exceeding the size of a detector unit. This was followed by approximating the integration over the strip/mesh-element union. This approximation is product of: (i) strip width, (ii) intersection length of a ray central to strip with a mesh element, and (iii) the response and expansion function evaluated at midpoint of the intersection length. Reconstruction was performed using OSEM without regularization and with exact knowledge of the attenuation map. The method was evaluated using synthetic SPECT data generated using SIMIND Monte Carlo simulation software. Comparative quantitative and qualitative analysis included: bias, variance, standard deviation and line-profiles within three different regions of interest. We found that no more than two divisions per detector bin were needed for good quality reconstructed images when using a high resolution mesh.

KW - Attenuation compensation

KW - Mesh domain

KW - OSEM

KW - Reconstruction

KW - SPECT

KW - System matrix

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

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

U2 - 10.1117/12.844427

DO - 10.1117/12.844427

M3 - Conference contribution

AN - SCOPUS:84873895568

SN - 9780819480231

VL - 7622

BT - Progress in Biomedical Optics and Imaging - Proceedings of SPIE

ER -