Abstract: We propose a novel statistical method to extend Fermi-LAT catalogues of highlatitude -y-ray sources below their nominal threshold. To do so, we rely on the determination of the differential source -count distribution of sub -threshold sources which only provides the statistical flux distribution of faint sources. By simulating ensembles of synthetic skies, we assess quantitatively the likelihood for pixels in the sky with relatively low -test statistics to be due to sources, therefore complementing the source -count distribution with spatial information. Besides being useful to orient efforts towards multi -messenger and multi -wavelength identification of new -y-ray sources, we expect the results to be especially advantageous for statistical applications such as cross -correlation analyses.