Le portail francophone de la géomatique

Bonjour,

J'ai mis au moins 2 jours à faire ma requête de découpage alors si je peux vous en faire bénéficier c'est avec plaisir.

Tout d'abord voici les conclusions que j'ai pu tirer de mes tests :
-un point n'est jamais positionné exactement sur un tronçon de ligne, tout n'est qu'une question de tolérance.
- il faut effectivement regrouper les points en une couche sinon ca duplique les résultats en sortie.

voici la requete :

select couche_de_ligne.champ1,couche_de_ligne.champ2,
ST_SNAP((ST_DUMP(st_difference(couche_de_ligne.the_geom,point))).geom,all_point,0.1) as the_geom 
from (select ST_Multi(ST_Union(st_expand(couche_de_point.the_geom, 0.05))) as point from couche_de_point ) as t1, 
couche_de_ligne,(select ST_MULTI(ST_COLLECT(couche_de_point.the_geom))as all_point from couche_de_point) as t2

et voici l'explication :
Afin de palier au problème de tolérance, je convertis mes points en buffer (très petit 5cm). J'ai préalablement assemblé ces points en une seule couche de multipoint (avec ST_UNION). Je découpe avec ST_Différence puis ST_Dump comme dans votre solution 1, et je rajoute un ST_SNAP à 10cm, pour recoller les extrémités de mes lignes découpées sur les points de découpage eux-même, mais que j'ai assemblé avec ST_COLLECT ET non ST_UNION afin de leur conserver des géométries différenciées.

Nicolas GRANIER

Mais de rien. J'ai fait moi même pas mal de veille sans trouver un truc qui fonctionne à 100%, alors j'ai assemblé plusieurs morceaux de sql pour parvenir à cette solution.
A bientôt.

Nicolas GRANIER

Re

Je me suis demandé s'il etait possible de realiser la meme chose avec les fonctions de referencement linéaire, qui permettent de couper des lignes et de projeter des points sur des lignes proches, realisant un snapping tres precis et tres rapide.

Le principe est de trouver pour chaque ligne, la position de chaque point de decoupage par rapport a la ligne: 0 represente le debut de la ligne, 1 la fin de la ligne. Il s'agit de la localisation du point sur la ligne.

Puis, pour pouvoir realiser le decoupage des lignes initiales en se servant des differentes localisations des points sur la ligne, une table est construite en creant un champ index, incremental ( partie "rank() over (partition...)" ) . Cet index va permettre de fabriquer un record avec les champs: geom, localisation, localisation + 1

Enfin, les lignes initiales sont découpés par les localisations des points (st_line_substring), en realisant les jointures qui vont bien:

-- premiere table des positions des points par rapport a "leurs" lignes
with locus as (
        -- on genere la position du debut de la ligne (index 0)  pour le premier segment à découper
    select l.gid, 0 as l
    from lines l

    UNION ALL

        -- on genere la position de la fin de la ligne (index 1)  pour le dernier segment à découper
    select l.gid, 1 as l
    from lines l

    UNION ALL

        -- calcul des positions des points par rapport aux lignes
    select l.gid, st_line_locate_point(l.geom, p.geom) as  l
    from lines l, points p
    where st_dwithin(l.geom, p.geom, 0.01)
    
    order by gid, l
), 
-- deuxieme table qui ajoute un index croissant par position croissante, pour permettre de faire la jointure plus tard sur elle-meme:
loc_with_idx as (
    select gid, l, rank() over (partition by gid order by l) as idx
    from locus
) 
-- decoupage final de chaque ligne avec les positions des points sur la ligne, grace a st_line_substring
select l.gid, st_line_substring(l.geom, loc1.l, loc2.l) as geom 
from loc_with_idx loc1, loc_with_idx loc2, lines l
where l.gid = loc1.gid
and loc1.gid = loc2.gid
and loc2.idx = loc1.idx+1;

L'avantage de cette methode est que le snapping entre les points et les lignes est super rapide et precis, quelque soit la precision.

En terme de vitesse, c'est tres performant ! : 90ms avec le meme jeu de donnée que dans le post précédent.

Bonjour,

Oui 6h de traitement pour cette volumetrie, c'est un peu trop.
la methode la plus rapide (referencement lineaire) prend environ 20 s pour une volumetrie de 109133 lignes des communes croisées avec 800 000 points sur ces lignes.
La methode WITHIN prend environ 120% de plus.

Nicolas