Le portail francophone de la géomatique

Renaud · Tue 27 November 2018 17:13

Bonjour,

Nous cherchons à déterminer le cours d'eau sur lequel se situent un ensemble de points de mesure.
Je m'attendais à trouver une fonction du genre ST_nearestneighbour, mais impossible de mettre la main dessus.

On s'en sort par un calcul de la matrice des distances point-rivière suivi d'un "row_number() + partition by", suivi d'un filtre.
Mais c'est un peu laborieux.

Est-ce que vous auriez une solution plus élégante ?

Merci

Dernière modification par Renaud (Tue 27 November 2018 17:19)

tumasgiu · Tue 27 November 2018 17:28

Utilisation d'une requête latérale ?

Code:

SELECT id_point, 
( SELECT id_riviere FROM
    (SELECT st_distance(point.geom, riviere.geom) d,  id_riviere 
     FROM riviere ORDER BY point.geom <-> riviere.geom LIMIT 10) foo
  ORDER BY d LIMIT 1
)
FROM point

Dernière modification par tumasgiu (Tue 27 November 2018 17:32)

Nicolas Ribot · Tue 27 November 2018 17:53

Bonsoir,

Et comme on me le faisait remarquer sur le forum Postgis, avec un postgis récent, pas besoin d'utiliser st_distance: geom <-> geom renvoie la vraie distance:

SELECT id_point,
( SELECT id_riviere FROM
(SELECT point.geom <-> riviere.geom d, id_riviere
FROM riviere ORDER BY point.geom <-> riviere.geom LIMIT 10) foo
ORDER BY d LIMIT 1
)
FROM point

Nicolas

tumasgiu · Tue 27 November 2018 18:09

Super, j'ai du le lire mais j'ai aussi du l'oublier !
A ce moment là, la sous-requête de l'expression latérale est superflue non ?

Dernière modification par tumasgiu (Tue 27 November 2018 18:10)

Renaud · Wed 28 November 2018 09:53

Bonjour,

Merci pour le retour,

Je teste la solution

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 Tue 27 November 2018 17:13

PostGIS - Plus proche voisin

#2 Tue 27 November 2018 17:28

Re: PostGIS - Plus proche voisin

Code:

#3 Tue 27 November 2018 17:53

Re: PostGIS - Plus proche voisin

#4 Tue 27 November 2018 18:09

Re: PostGIS - Plus proche voisin

#5 Wed 28 November 2018 09:53

Re: PostGIS - Plus proche voisin

Pied de page des forums