Searching for as yet undetected γ -ray sources is one of the main stated goals of the Fermi Large Area Telescope Collaboration. In this paper, we explore the capability of a filtering method based on a finite mixture of von Mises–Fisher distributions. The proposed procedure is specifically designed to handle data with support on the unit sphere. The assumption of a parametric model for each high energy emitting source allows us to derive an explicit expression for both the direction of the sources and their angular resolutions. The corresponding measures are based on the directional mean and the quantiles of the single mixture components. Sound criteria of model selection can provide an automatic way to determine the number of detected sources. Additionally, a likelihood-ratio test is developed to evaluate their significance. The procedure is tested on simulated data sets of photon emissions from high energy sources within the energy range 10 − 1 , 000  GeV . A real data example consisting of a sample of the Fermi LAT data collected over a period of about 7.2 years within the energy range 10 − 1 , 000  GeV , in a subregion of the γ -ray sky, is furthermore provided.