Le portail francophone de la géomatique

Bonjour,

Si les lignes se touchent, vous pouvez utiliser un groupement spatial:
group by st_touches(ligne1, ligne2)

si elles ne se touchent pas, le groupement se ferait plutot avec l'operateur distance, si tant est que les lignes a joindre entre elles sont a une distance minimale entre elles. Par exemple, si deux lignes sont plus proches de X metres, alors il faut les joindre; la condition serait:
group by st_distance(ligne1.startpoint, ligne2.startpoint) < X or st_distance(ligne1.endpoint, ligne2.startpoint) < X ...
etc... pour les combinaisons startpoint/endpoint des deux lignes.

(Pas sur d'etre tres clair )
Je tente de trouver un vrai exemple...

Nicolas

Ou alors grouper sur le critère d'identité entre coordonnées des extrémités. ST_UNION étant l'équivalent Postgis de SDO_AGGR_CONCAT_LINES chez Oracle, je pense moi aussi que c'est la bonne solution.

Pardon, j'ai ete un peu court dans mes explications:

plutot une requete du style:

select
case when st_dwithin(st_endpoint(t1.geometry), st_startpoint(t2.geometry), 10) then
st_linemerge(
    st_collect(
        st_addpoint(
            t1.geometry, 
            st_startpoint(t2.geometry)), 
        t2.geometry))
when st_dwithin(st_startpoint(t1.geometry), st_endpoint(t2.geometry), 10) then
st_linemerge(
    st_collect(
        st_addpoint(
            t2.geometry, 
            st_startpoint(t1.geometry)), 
        t2.geometry)) 
END
from toto t1, toto t2 
where t1.id <> t2.id 
and  (st_dwithin(st_startpoint(t1.geometry), st_startpoint(t2.geometry), 10) 
  or st_dwithin(st_startpoint(t1.geometry), st_endpoint(t2.geometry), 10)
  or st_dwithin(st_endpoint(t1.geometry), st_startpoint(t2.geometry), 10)
  or st_dwithin(st_endpoint(t1.geometry), st_endpoint(t2.geometry), 10));

Ici donc, je considere que les lignes a joindre ne se touchent pas.

Le group by permet de garder les lignes correspondant au critere de proximité.

La partie select regarde les conditions pour ajouter un point a une des deux lignes,histoire de les rendre jointives.
Le linemerge fabrique ensuite une seule ligne a partir des deux lignes rendues jointives.

A affiner pour:

• Supporter les multi
• merger le resultat final pour ne pas avoir des lignes qui se superposent.

Je joins une image: en haut, e reseau initial, en bas: le resultat de la requete.

Nico

PS. (merci, ca va me faire un bon sujet de TP pour la prochaine formation Postgis )

Hihi, oui effectivement. Je viens de me cramer, là

A propos de ces fonctions non existantes dans postgis: les dev postgis se concentrent surtout sur les fonctions "de base", style, Union, linemerge etc. bien définies dans un cadre précis.
Dès qu'il s'agit de faire des traitements plus SIG, avec prise en compte de bcp de conditions, c'est plus difficile de faire des fonctions génériques.

Mais je suis d'accord avec toi qu'un entrepot central ou ces procedures pourraient etre telechargées, ca serait bien. C'est l'esprit du Wiki Postgis.

(là, par ex, j'ai un peu triché dans la mesure ou j'ai retourné (st_reverse) les lignes qui n'etaient pas dans le bon sens et me donnaient des résultats faux. Ce retournement pourrait se faire dans la requete je pense).

Nicolas

Bonsoir,

j'ai un problème similaire au sujet évoqué ci dessus. A savoir:
  -une table qui comprend un ensemble de linestring
  - certaines se touchent, d'autres non

Je souhaiterai obtenir une table qui agrège en une linestring toutes les linestring  qui se touchent.
j'ai utilisé la requête proposée en remplaçant st_dwithin par st_touche. Mais je n'obtient pas le résultat attendu. Ainsi pour 3 linestring (1,2 et 3) qui se touchent j'obtiens
linestring 1
linestring 1+2
linestring 1+2+3
linestring 2 +3

alors que j’attends
linestring 1 +2 +3


Merci d'avance pour votre aide

Bonjour,

Hmm, solution elegante. Mais je n'avais jamais remarqué que st_union sur des lignes qui se touchent ne crée pas une nouvelle ligne soudant les  autres objets, mais plutot une MULTILINESTRING ou les elements constitutifs restent les "petites" lignes !

St_linemerge permet de "souder" ces lignes:

select (st_dump(st_linemerge(st_union(the_geom)))).geom from lines;

Nicolas

En fait lorsque cela bug lorsque les linestrings se croisent. De ce fait il faut travailler avec les extrémités et cela semble fonctionner

J'ai donc exécuté la requête suivante

select st_linemerge(st_collect(cs1.the_geom))
from (select the_geom, st_startpoint(cs_s.the_geom) as linepoint
          from [ma-table] cs_s
         UNION DISTINCT
         select the_geom, st_endpoint(cs_e.the_geom) as linepoint
          from [ma-table] cs_e
) as cs1
group by linepoint
having count(*) > 1

Bonjour,
je relance la discussion.
j'ai un cours d'eau avec des multilignes disjointes, forcément quand je fais un linemerge celui me renvoi la multiligne initiale.

select row_number () over(),cdeau_ppx1.topooh,st_linemerge(st_union(cdeau_ppx1.geom))as geom,cdeau_ppx1.code_ctxt
from pdpg.cdeau_ppx1
group by cdeau_ppx1.code_ctxt,cdeau_ppx1.topooh
order by geom desc

avez vous une solution?

Bonjour, merci de votre réponse
oui je comprend, l'idée si je comprend bien est d'utiliser st_dump afin d'exploser la multiligne et vérifier avec dwithin si elles sont connectées.
ci joints les infos

Si les lignes sont ordonnées comme il faut dans chaque multilinestring (ML), comme dans votre exemple, vous pouvez combler les gaps en dumpant les points de chaque ML et en reconstruisant une linestring avec: ca connecte alors tous les points:

Dump des ML en LN, puis dump des LN en points, puis construction d'un tableau des points ordonnés pour chaque LN à construire:

with tmp as (
    select id, (st_dump(geom)).path[1] as idln, (st_dump(geom)).geom
    from testln
), tmp1 as (
    select id, idln, (st_dumppoints(geom)).path[1] as idpt, (st_dumppoints(geom)).geom
    from tmp
) select id, st_makeLine(array_agg(geom order by idln, idpt)) as geom
from tmp1
group by id;

Si les LN ne sont pas ordonnées, il faut les ordonner avant de les joindre. Pour cela, il y a une solution générique avec WITH RECURSIVE pour traiter ligne par ligne (pas mal d'exemples de cette construction sur le forum)

Nicolas

Dernière modification par Nicolas Ribot (Thu 19 November 2020 11:06)

Bonsoir,

startpoint/endpoint similaires, c'est dans la même linestring (ligne fermée), ou entre deux lignes de la collection (lignes se touchant sans gap) ?

Pour identifier les doublons,  vous pouvez faire un group by geom (attention à votre dataset, car dans ce cas les bbox des géométries sont prises en compte: meme bbox => meme geom. Il faut donc etre sur que deux geom du dataset ne peuvent pas avoir la meme exacte bbox)

La requête que je vous ai donnée génère des doublons ?

Nicolas