Traditional discrimination techniques for Euler deconvolution use only the color spectrums of structural indexes, without considering the spatial distribution characteristics and inherent relationships among the Euler solutions to separate adjacent causative sources. In the present study, a new approach was developed for discriminating uncorrelated Euler solutions from coherent solutions based on the focusing levels indicated by the probability density distributions generated using multivariate kernel density estimations (KDE). A novel multiple coverage technique was proposed by using a series of different sized moving windows over gridded gravity data, which formed tight clusters of Euler solutions for different sized causative sources. The results of the probability density distributions were obtained using a 3-D KDE method for the Euler solution subsets {x, y, z} of synthetical models, and real data from a survey conducted in British Columbia (Canada) which had successfully established more credible and meaningful geological models when compared with three other subsets.