Computationally Efficient Collimator-Detector Response Compensation in High Energy SPECT using 1D Convolutions and Rotations

Modeling of the collimator-detector response (CDR) in SPECT reconstruction enables improved resolution and accuracy, and is thus important for quantitative imaging applications such as dosimetry. The implementation of CDR modeling, however, can become a 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. This work proposes a 1D convolution and rotation-based CDR model that reduces reconstruction times but maintains consistency with models that employ 2D convolutions. To enable open-source development and use of these models in image reconstruction, we release a SPECTPSFToolbox repository for the PyTomography project on GitHub.