Le portail francophone de la géomatique

Bonjour,

la fonction st_distance() fait le travail. Il reste à extraire la valeur minimale pour chaque parcelle.
Ci-dessous j'utilise la condition de jointure st_dwithin avec une distance de 5000 mètres  pour ne pas faire le produit cartésien des deux tables, en considérant que chaque parcelle est à moins de 5000 m d'une rivière (à adapter à votre cas).

SELECT parcelle.id, min(st_distance(parcelle.geom, riviere.geom))
FROM parcelle JOIN rivieres ON st_dwithin(parcelle.geom, riviere.geom,5000)
GROUP BY parcelle.id

La piste de la recherche du plus proche voisin est intéressante aussi.

Bonjour,

Comme le dit Mathieu, PostGIS a un opérateur distance dédié à ce genre de requetes: <->:
Cela permet de chercher les n voisins plus proches (KNN) des géométries.

L'avantage par rapport a min(distance) est double:
• énorme performance qq ms contre plusieurs secondes ou minutes pour min(distance)
• pas besoin de préfiltrer par une distance qui dépend des données: l'opérateur marche tjs qq soit la distance min entre géométries.

Pour utiliser cet opérateur, il faut faire une jointure latérale entre les tables qui permet de référencer les colonnes d'une table depuis une sous requête:

select p.id, t.riv_id, t.dist
from parcelle p cross join lateral (
        select r.id as riv_id, p.geom <-> r.geom as dist
        from riviere r
        order by p.geom <-> r.geom
        limit 1
    ) as t;

La différence de perf sur une table parcelle de 300 000 et rivieres de 400 geoms:
knn: 128 ms
min(distance): 7200 ms

Nicolas