Fast and Accurate Collimator-Detector Response Compensation in High-Energy SPECT Imaging with 1D Convolutions and Rotations

Abstract

Modeling of the collimator-detector response (CDR) in SPECT reconstruction enables improved resolution and more accurate quantitation, especially for higher energy imaging (e.g.Lu-177 and Ac-225). Such modeling, however, can pose a significant computational bottleneck when there are substantial components of septal penetration and scatter in the acquired data, since a direct convolution-based approach requires large 2D kernels. The present work presents an alternative method for fast and accurate CDR compensation using a linear operator built from 1D convolutions and rotations (1D-R). To enable open-source development and use of these models in image reconstruction, we release a SPECTPSFToolbox repository for the PyTomography project on GitHub.