Le portail francophone de la géomatique

Bonjour, C'est bien la fonction st_union qu'il faut utilisé. Mais le problème c'est de trouver un moyen qui détecte d'abord toutes les lignes qui se trouvent entre deux point et ensuite faire le st_union. Donc je me suis bloqué à l'étape: trouver toutes les lignes qui se trouvent entre deux points.

Cordialement

Bonjour,

En effet, la fonction st_relate permet de voir les relations de géométrie, elle peut contribuer à solutionner mon problème. Le seul problème c'est au moment où je fais le parcours (avec la méthode with recursive)  que j'arrive pas à fusionner les lignes sans point à leur extrémité. J'ai mis en pièce jointe le problème et la solution souhaitée.

Bonjour,
On ne sait pas vraiment le but métier de l'exercice mais ca n'enlève pas son intérêt.
Tu peux peut-être t'en sortir en accrochant les lignes "non accrochées" entre elles avec St_Snap si tu es capable de définir un seuil de tolérance.
Sinon, pour attribuer à chaque ligne, le couple de point auquel il correspond, tu pourrais créer un segment qui matérialise le couple de point (st_makeline) et faire un petit buffer et faire une relation entre ce buffer et tes lignes.
Difficile à expliquer! Le challenge m'intéresse mais j'ai peu de temps.

Bonjour Nicolas,

En fait j'ai une couche ligne et une couche point. Par exemple j'ai des lignes (par exemple ligne1,ligne2,ligne3,ligne4,ligne5) qui se trouvent entre deux points (point a et point b), en ce moment je cherche à fusionner ces cinq lignes pour que j'aie une ligne entre ces deux points. Je rappelle que le nombre de ligne entre deux points peut varié (minimum deux lignes). Dans mon exemple, j'ai montrer les lignes à fusionner (je les ai espace pour uniquement montrer que sont à fusionner, donc pas de problème de snap) en une seule ligne.

Bonjour ppluvinet,

Je tiens juste à rappeler que j'ai pas des problèmes de snap, c'est juste à titre d'information afin de montrer les lignes à fusionner Par contre ton idée de faire un st_makeline avec un buffer et regarder les lignes dans le buffer me parait une solution. Je vais le tester et voir.
J'ai testé ta méthode mais le souci est, avec la fonction st_makeline,il me crée beaucoup de ligne,y a t-il un moyen pour paramétrer ça afin qu’il me crée que des segments souhaiter.

Dernière modification par fassa222 (Tue 27 June 2017 14:15)

Bonjour,

Je relance ce sujet afin de savoir s'il y a eu une solution apportée ?

Je pense avoir besoin de cette astuce pour mon jeu de données, que j'ai déjà commencé à traiter (et discuté avec vous) ici : https://georezo.net/forum/viewtopic.php?id=104853

J'ai donc une série de lignes (des portions de cours d'eau), auxquelles je "colle" des points (ptO). Pour l'instant, je découpe ensuite ces lignes avec ces points.
Le résultat est parfait, sauf que j'aurais besoin de n'obtenir qu'un tracé de cours d'eau entre chaque ptO; c'est à dire fusionner les lignes qui étaient déjà coupées avant la découpe par les ptO, et qui sont donc en deux (ou plus) morceaux entre mes ptO.
Sachant qu'entre certains ptO, la ligne est bien unique.


J'avais donc pensé à tout ce qui a été évoqué ici :

--->> avant de découper les lignes avec les ptO: fusionner toutes mes lignes de cours d'eau en une seule (st_union) pour ensuite découper uniquement au niveau des points, mais je perds des informations utiles sur les lignes ET je ne peux pas appliquer la découpe car j'obtiens une multilinestring au lieu d'une line

DROP TABLE IF EXISTS unionligne CASCADE;
CREATE TABLE unionligne AS 
(SELECT ST_Union(geom) as geom 
FROM lignes_dumped -- couche sur laquelle un st_dump(geom).geom a été réalisé
);
ALTER TABLE unionligne ADD COLUMN id SERIAL PRIMARY KEY;

--- découpage -- Merci N.RIBOT sur Georezo :)
        -- premiere table des positions des points par rapport a "leurs" lignes
DROP TABLE IF EXISTS locus CASCADE;
CREATE TEMP TABLE locus as (
            -- on genere la position du debut de la ligne (index 0)  pour le premier segment à découper
    select l.id, 0 as l
    from unionligne l

    UNION ALL

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

    UNION ALL

            -- calcul des positions des points par rapport aux lignes
    select l.id, (st_linelocatepoint(l.geom, p.geom)) as  l
    from unionligne l, ptO p
    where st_dwithin(l.geom, p.geom, 0.1) 
AND st_linelocatepoint(l.geom, p.geom) NOT IN (0,1)
GROUP BY l.id, l
    order by  id, l
  
);

-- deuxieme table qui ajoute un index croissant par position croissante, pour permettre de faire la jointure plus tard sur elle-meme:
DROP  TABLE IF EXISTS loc_with_idx CASCADE;
CREATE TEMP TABLE
loc_with_idx as (
    select id, l, rank() over (partition by id 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--

DROP TABLE IF EXISTS splitang CASCADE;
CREATE TABLE splitang AS (select l.id, l.gid, l.__gid ,  l.eu_cd  ,  l.name  ,  l.ssbv ,  l.nom ,  l.localisati ,   l.code_hydro ,  l.cd ,   st_linesubstring(l.geom, loc1.l, loc2.l) as geom 
from loc_with_idx loc1, loc_with_idx loc2,  unionligne l
where l.id = loc1.id
and loc1.id = loc2.id
and loc2.idx = loc1.idx+1
GROUP BY l.id, l.gid, l.__gid ,  l.eu_cd  ,  l.name  ,  l.ssbv ,  l.nom ,  l.localisati ,  l.code_hydro ,  l.cd ,   l.geom, loc1.l, loc2.l );

ALTER TABLE splitang ADD COLUMN idsplit SERIAL PRIMARY KEY;

Qui donne donc le message d'erreur suivant :

ERREUR:  line_locate_point: 1st arg isn't a line
********** Erreur **********

ERREUR: line_locate_point: 1st arg isn't a line
État SQL :XX000

--->> au lieu de découper les lignes avec les pt : créer une couche de points à partir des lignes de cours d'eau avec ST_DumpPoints() ; utiliser st_makeline() entre cette nouvelle couche de points et la couche pt (initialement de découpe, qui servirait ici de points d'extrémités) mais je n'arrive pas à identifier de champ de liaison /group by dans mes données...
Sachant que je peux récupérer toutes les infos de la ligne coupée par le pt sur la couche pt.
Mais ensuite, comment faire la jointure st_makeline(geom1, geom2)?


--->> après la découpe des lignes (splitang) avec les pt : fusionner les lignes situées entre les pt ...
J'ai normalement réussi à identifier les lignes "orpheline" de pt à l'une ou aux deux extrémités (= celles à fusionner entre elles si une de leur extrémité sans pt est proche d'une autre) avec les requêtes suivantes :

DROP TABLE IF EXISTS extremsplita CASCADE;
CREATE TEMP TABLE extremsplita AS (
SELECT 'a'||(row_number() over()) as idextrem, idsplit,  __gid , ST_LineInterpolatePoint(s.geom, 0) as geom 
FROM splitang s);
DROP TABLE IF EXISTS extremsplitb CASCADE;
CREATE TEMP TABLE extremsplitb AS (
SELECT 'b'||row_number() over() as idextrem, idsplit,  __gid , ST_LineInterpolatePoint(s.geom, 1) as geom 
FROM splitang s
);


SELECT  a.idsplit , a.idextrem as a,oa.snapcd, b.idextrem as b, ob.snapcd, l.geom
FROM splitang l
LEFT OUTER JOIN extremsplita a
ON l.idsplit = a.idsplit  
LEFT OUTER JOIN ptO oa
ON ST_Dwithin(oa.geom, a.geom, 0.1)
LEFT OUTER JOIN extremsplitb b
ON l.idsplit = b.idsplit
LEFT OUTER JOIN ptO ob
ON ST_Dwithin(ob.geom, b.geom, 0.1)
WHERE (oa.snapcdu IS NULL OR ob.snapcd IS NULL)
GROUP BY a.idsplit , a.idextrem,oa.snapcd, b.idextrem, ob.snapcd, l.geom

Maintenant je ne sais pas comment "fusionner" ces lignes entre elles en une entité, seulement si elles partagent un point d'extrémité en commun (st_dwithin(..... 0.1) ?) ET qu'elles ont au PLUS un pt à leur extrémité, afin d'obtenir plusieurs entité de lignes groupées...

Pensez-vous qu'une ou l'autre méthode serait la plus adaptée ?
Avez-vous une ou plusieurs idées pour la /les développer ?

Merci d'avance pour vos retours !

Dernière modification par meonais (Mon 31 July 2017 16:42)