A list-mode multi-energy window low-count SPECT reconstruction method for isotopes with multiple emission peaks.

SPECT provides a mechanism to perform absorbed-dose quantification tasks for $\alpha$-particle radiopharmaceutical therapies ($\alpha$-RPTs). However, quantitative SPECT for $\alpha$-RPT is challenging due to the low number of detected counts, the complex emission spectrum, and other image-degrading artifacts. Towards addressing these challenges, we propose a low-count quantitative SPECT reconstruction method for isotopes with multiple emission peaks. Given the low-count setting, it is important that the reconstruction method extract the maximal possible information from each detected photon. Processing data over multiple energy windows and in list-mode (LM) format provide mechanisms to achieve that objective. Towards this goal, we propose a list-mode multi-energy window (LM-MEW) OSEM-based SPECT reconstruction method that uses data from multiple energy windows in LM format, and includes the energy attribute of each detected photon. For computational efficiency, we developed a multi-GPU-based implementation of this method. The method was evaluated using 2-D SPECT simulation studies in a single-scatter setting conducted in the context of imaging [$^{223}$Ra]RaCl${_2}$. The proposed method yielded improved performance on the task of estimating activity uptake within known regions of interest in comparison to approaches that use a single energy window or use binned data. The improved performance was observed in terms of both accuracy and precision and for different sizes of the region of interest. Results of our studies show that the use of multiple energy windows and processing data in LM format with the proposed LM-MEW method led to improved quantification performance in low-count SPECT of isotopes with multiple emission peaks. These results motivate further development and validation of the LM-MEW method for such imaging applications, including for $\alpha$-RPT SPECT.