Le portail francophone de la géomatique

Bonjour Pascal,

Merci ! J'ai appliqué la formule de Nicolas Granier et ça fonctionne.

Mais pour info, ça ne pouvait pas fonctionner en complétant ce que j'avais déjà fait ? Avec le ST_Split ?

Bonne journée

Bonsoir,
Je dois pouvoir vous aider
Dans la requete que vous citez (il en manque un morceau dans votre citation), le terme "point" designe en realité un agglomérat de buffer carré autour de tous les points assemblés en une seule géométrie. Le terme point provient de la sous requete suivante passée dans le from

(select ST_Multi(ST_Union(st_expand(couche_de_point.the_geom, 0.05))) as point from couche_de_point ) as t1

Pourquoi un buffer ? car topologiquement parlant, pour utiliser vos termes si votre pont n'est pas sur la ligne celle ci ne pourra jamais être coupée. On fait donc un buffer dont la distance correspond à la réponse à la question " Dans mes données, jusqu'à quelle distance dois-je considérer qu'un pont enjambe un cours d'eau ? 1mm, 1m, 100m ?" Si je décide de faire un buffer de 1m pour faire un carre de 1m, les polygones vont pouvoir intersecter le cours d'eau distant de moins d'un mètre.

Pourquoi un buffer carré ? Là je dois dire que je ne m'en souviens plus mais il doit y avoir une histoire non négligeable de vitesse de calcul, sur un buffer carré vous n'avez toujours que 4 points qui constituent le polygone d'enveloppe, sur un rond vous en avez une multitude donc derrière en vitesse de calcul vous êtes cuits avec vos ronds.

Pourquoi un agglomérat de petits carrés ?

Ce qu'il faut garder en tete lors d'une requete SQL c'est que les opérations que vous lancez se font pour chaque enregistrement (chaque ligne de la table), et si vous faites une opération de type st_intersects entre deux tables les résultats en sortie s'en trouvent multipliés ce qui donne au final quelque chose comme ça :

ligne 1 table 1 intersects ligne 1 table 2
ligne 1 table 1 intersects ligne 2 table 2
ligne 1 table 1 intersects ligne 3 table 2
ligne 2 table 1 intersects ligne 1 table 2
ligne 2 table 1 intersects ligne 2 table 2

Tout ca pour dire que la requete fait apparaitre le terme st_multi(st_union ...) afin de n'avoir qu'une seule ligne de geométrie de point face à chaque ligne géométrique

si l'on décompose maintenant  la requete de découpage elle même :
st_difference (geométrie d'une ligne, géométrie de tous les carrés/buffers autour des points)
Admettons que sur un cours d'eau se trouvent deux ponts. Vous avez donc votre ligne qui se retrouve superposée avec deux petits carrés  distants mais géométriquement ne représentant qu'un seul objet géométrique. L'outil St_différence va effacer les morceaux de ligne qu'il y a sous les carrés. Vous avez donc à la fin de cette étape une ligne avec deux trous dedans mais stockée dans une seule géométrie multiligne.

st_dump(...).geom permet d'éclater la multiligne en 3 lignes (comme on a fait deux trous dedans, il y a donc 3 morceaux). Sans cette opération vous auriez autant de ligne dans la table d'arrivée que dans la table de départ.

enfin
ST_SNAP(mes morceaux de ligne, tous mes points, distance) permet de recoller les morceaux entre eux. Je m'explique. Pour l'instant vous avez bien des lignes mais elles ne sont pas jointives car vous avez fait un buffer sur les ponts et effacé la ligne sous le buffer. Si vous avez fait un buffer de 1 m, vous avez donc des trous dans vos cours d'eau d'environ 1m. L'outil St Snap permet de venir coller une géométrie sur une autre. On vient donc coller l'extrémité de nos cours d'eau sur nos ponts. Là encore on se retrouve avec le cas de 2 tables qui vont discuter ensemble, il faut donc faire un agglomérat des ponts pour n'avoir qu'une seule géométrie de ponts mais sans réaliser de buffer cette fois.

A+
Nicolas GRANIER

Dernière modification par Nicolas Granier (Tue 21 March 2017 23:29)