Energy and position dependence of scatter functions in high resolution PET

The authors have shown in previous works that distinct non-stationary analytical scatter kernels can be extracted from line source measurements and used to independently subtract object scatter and subtract or restore detector scatter in high resolution PET. In this work, the dependence of the scatter components on energy threshold and source position was investigated. Line source measurements were acquired in multispectral mode using the Sherbrooke PET simulator. Scatter parameters were extracted from data cumulated in energy windows with a lower threshold varying from 129 keV to 516 keV in steps of 42 keV, and a fixed upper threshold of 644 keV. Decreasing the lower threshold from 344 keV to 129 keV increases the trues by only 25%, but increases object scatter by 136% and almost triples detector scatter. A gain in efficiency by a factor of 2 or more would result from recovering the latter by restoration in the broad window. The intensity and shape of the scatter functions for both object and detector are shown to have a significant dependence on energy and position. This dependence needs to be taken into account in the design of kernels for accurate scatter correction over a broad energy range