Missing head and color banding in low-count SPECT reconstructions

TY - JOUR

T1 - Missing head and color banding in low-count SPECT reconstructions

AU - de Nijs, Robin Vedel

AU - Jensen, Björn Neumann

AU - Mortensen, Jann

PY - 2014/12/1

Y1 - 2014/12/1

N2 - Due to low counts in an 111In single photon emission computed tomography (SPECT) scan, a large part of the head was missing in the reconstructed images on Philips Extended Brilliance Workspace (EBW) and IntelliSpace Portal (ISP) workstations. This problem occurred for the ordered subsets expectation maximization (OSEM) algorithm with and without resolution recovery (Astonish), but not for filtered backprojection (FBP) or maximum likelihood expectation maximization (MLEM). There were also underflow problems because the images are stored as integers resulting in a loss of intensity resolution and color banding. Philips EBW2.0 and ISP5.02 workstations upscale low-count images, but the result is not always optimal, for example, in the case of low counts in one part and more counts in another part of an image. On these workstations, the missing head artefact problem could be resolved by applying a Hann pre-filter (with a cutoff at the Nyquist frequency, which only influences the filtering) in the reconstruction process. Upscaling of the projection data prior to reconstruction did not recover the head in the images, neither did limiting the reconstructed volume to the low-count part of interest. Underflow problems were partially solved by the new version 2.0 of the Philips EBW and ISP stations, although situations could arise where underflow still poses a problem. A solution for the underflow problems is to upscale the raw projection data before reconstruction. While this results in a pure upscaling of the FBP reconstruction, the effect in iterative statistical reconstruction is not only upscaling of the intensities because the assumption of Poisson statistics of the data is violated. However, the influence of this last matter seems limited. Reconstruction of studies with low counts in relevant areas should be performed with care. Reconstruction artefacts and scaling issues can easily arise.

AB - Due to low counts in an 111In single photon emission computed tomography (SPECT) scan, a large part of the head was missing in the reconstructed images on Philips Extended Brilliance Workspace (EBW) and IntelliSpace Portal (ISP) workstations. This problem occurred for the ordered subsets expectation maximization (OSEM) algorithm with and without resolution recovery (Astonish), but not for filtered backprojection (FBP) or maximum likelihood expectation maximization (MLEM). There were also underflow problems because the images are stored as integers resulting in a loss of intensity resolution and color banding. Philips EBW2.0 and ISP5.02 workstations upscale low-count images, but the result is not always optimal, for example, in the case of low counts in one part and more counts in another part of an image. On these workstations, the missing head artefact problem could be resolved by applying a Hann pre-filter (with a cutoff at the Nyquist frequency, which only influences the filtering) in the reconstruction process. Upscaling of the projection data prior to reconstruction did not recover the head in the images, neither did limiting the reconstructed volume to the low-count part of interest. Underflow problems were partially solved by the new version 2.0 of the Philips EBW and ISP stations, although situations could arise where underflow still poses a problem. A solution for the underflow problems is to upscale the raw projection data before reconstruction. While this results in a pure upscaling of the FBP reconstruction, the effect in iterative statistical reconstruction is not only upscaling of the intensities because the assumption of Poisson statistics of the data is violated. However, the influence of this last matter seems limited. Reconstruction of studies with low counts in relevant areas should be performed with care. Reconstruction artefacts and scaling issues can easily arise.

U2 - 10.1186/2197-7364-1-10

DO - 10.1186/2197-7364-1-10

M3 - Journal article

C2 - 26501452

SN - 2197-7364

VL - 1

SP - 1

EP - 5

JO - E J N M M I Physics

JF - E J N M M I Physics

IS - 1

M1 - 10

ER -