Le portail francophone de la géomatique

vincent52 · Thu 21 September 2017 09:20

Bonjour à tous,

Ayant jusqu'à maintenant travaillé uniquement avec Spatialite pour effectuer des requêtes spatiales, je souhaite me mettre à Postgis pour pouvoir coupler SQL et R.
Bref, j'ai deux couches de points a et b, et je souhaite trouver le point b le plus proche de chaque point a.

Je veux un output du style :id_a, id_b, distance_mini.
Sur spatialite, je fais ce code :

Code:

select pa.id as id_a, pb.id as pb_id, min(st_distance(pa.geom,pb.geom)) as dist
from point_a as pa, point_b as pb
group by id_a

et ça fonctionne bien.

Mais lorsque je fais cette requête sur Postgis : "ERREUR: la colonne « b.id » doit apparaître dans la clause GROUP BY ou doit être utilisé dans une fonction d'agrégat"...

Quelqu'un a une idée ?

Merci d'avance et bonne journée

vincent52 · Thu 21 September 2017 11:42

Finalement j'ai trouvé la solution :

Code:

select distinct on (pa.id) pa.id as id_a, pb.id as id_b, st_distance(pa.geom,pb.geom) as dist
from point_a as pa, point_b as pb
where st_dwithin(pa.geom,pb.geom,8000)
order by pa.id,dist

Ici on calcule toutes les distances puis on les classe dans l'ordre croissant puis on selectionne la première lignes de chaque ID (avec distinct on). J'ai rajouté la clause where st_dwithin qui permet de diminuer le temps de calcul en donnant une distance limite, et ainsi éviter de calculer les distances pour la totalité des points B.

tumasgiu · Thu 21 September 2017 11:44

Salut,

bête et méchant :

Code:

SELECT
id_a, (array_agg(id_b ORDER by d))[1] id_b , MIN(d) d
FROM (
    SELECT pa.id id_a, pb.id id_b, st_distance(pa.geom, pb.geom) d
    FROM point_a pa, point_b pb
    ) foo
GROUP BY id_a

EDIT: je n'avais pas vu que vous vous étiez répondu.

Dernière modification par tumasgiu (Thu 21 September 2017 11:54)

Nicolas Ribot · Thu 21 September 2017 13:11

Bonjour

Postgis dispose d'un opérateur plus proche voision (KNN), tres performant:

http://postgis.net/docs/geometry_distance_knn.html

Nicolas

vincent52 · Thu 21 September 2017 13:38

Merci pour vos réponses

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Thu 21 September 2017 09:20

Requête Spatialite et Postgis

Code:

#2 Thu 21 September 2017 11:42

Re: Requête Spatialite et Postgis

Code:

#3 Thu 21 September 2017 11:44

Re: Requête Spatialite et Postgis

Code:

#4 Thu 21 September 2017 13:11

Re: Requête Spatialite et Postgis

#5 Thu 21 September 2017 13:38

Re: Requête Spatialite et Postgis

Pied de page des forums