Le portail francophone de la géomatique

• Si vous contrôlez les données OpenStreetMap et que vous travaillez sur une petite zone, il sera plus interessant et plus précis de travailler dans une projection locale (en France, la projection Lambert 93 est la projection officielle, métrique, qui couvre le territoire métropolitain).

• Si votre zone geographique est plus vaste que la france (et/ou si vous ne controlez pas le fond carto mais qu'il est issu d'un service externe) , il va etre difficile de trouver une projection qui soit précise sur un grand territoire. La plupart des services carto actuels utilisent cette projection monde EPSG:3857, Pseudo Mercator, qui déforme effectivement les surfaces, notamment près des pôles. Pour nos latitudes, les erreurs ne sont pas énormes. Le calcul de distance fait par ces services utilise souvent une formule précise pour calculer la distance, qui tient compte de l'ellipsoïde de la terre (cf. st_distanceSphere et st_distanceSpheroid, de PostGis).
Vous pourrez calculer des distances précises dans cette projection, si vous avez les bonnes coordonnées des points.

Concernant votre recherche de points, je vous invite plutot à utiliser la fonction st_dwithin (http://postgis.net/docs/manual-2.0/ST_DWithin.html).
C'est beaucoup plus performant que de passer par un buffer et st_intersects.

Elle supporte le type Geography (SRID=4326), que vous devriez utiliser pour le stockage et le traitement de vos données, vu que vous travaillez sur la terre entière.
Les calculs faits avec ce type d'objets sont précis car la forme de la terre est prise en compte. (http://postgis.net/docs/manual-2.0/usin … _Geography)

Concernant la précision de vos données, vous pouvez tester quelques points remarquables dont on trouve les coordonnées précises en WGS84 et regarder la différence avec vos données.