Le portail francophone de la géomatique

Salut,

tout d'abord,
le message d'erreur semble indiquer que la colonne géométrique
de votre table lignes_zone_test est de type multi linestring.
A ce moment là, il faudrait les extraire avec st_dump.

Ensuite, le résultat que vous obtenez vient du fait que
votre requête exécute pour st_split pour chaque combinaison
des géométries de vos deux tables.
Cette fonction, utilisée avec une géométrie et une ligne de découpe
qui sont géométriquement disjointes, aura pour résultat la géométrie originale.
Pour supprimer ces résultats inutiles, il faudrait rajouter une condition de jointure
entre vos deux tables pour ne retenir que les éléments qui s'intersectent :

SELECT ST_CollectionExtract(ST_Split(rect.geom, line.the_geom),3)
FROM public.polygones_zone_test as rect 
JOIN public.lignes_zone_test as line
ON st_intersects(rect.geom, line.the_geom)

Si vous ne l'avez pas déjà fait, vous devriez ajouter un index géométrique
sur l'une de vos deux tables (typiquement celles qui contient le plus
grand nombre de ligne), cela accélérera probablement le temps de traitement.

Mais je ne suis pas sur que cela soit le résultat que vous attendiez,
supposez un polygone traversé par plusieurs lignes,
chaque découpage se déroulera indépendamment avec toujours le polygone orignal en entier.
Pour découper avec une multi ligne, il vous  faudrait passer à postgis 2.5

Dernière modification par tumasgiu (Mon 13 May 2019 13:55)

Les polygones ayant disparus sont certainement ceux n'intersectant avec
aucune ligne.

Je crois comprendre à quel résultat vous voulez aboutir, il faudrait changer
un peu la requête :

SELECT 
    st_collectionextract(
        st_split(rect.geom,
                coalesce(ST_CollectionHomogenize(st_collect(line.the_geom))), 
                    'LINESTRING EMPTY'::geometry('LINESTRING')),
        3) AS geom
FROM
    polygones_zone_test as rect
LEFT JOIN
    public.lignes_zone_test as line
ON 
    st_intersects(rect.geom, line.the_geom)
GROUP BY
    rect.geom

L'idée est de regrouper, pour chaque polygone,
toutes les LINESTRING qui le traverse en une seule MULTILINESTRING
et d'appeler st_split avec celle ci.

Les polygones qui ne sont traversés par aucune ligne
sont quand même renvoyées grâce au LEFT JOIN,
et on transforme la MULTILINESTRING NULL qui leur sont associée
par une LINESTRING vide, pour que le st_split fonctionne et
renvoie le polygone tel quel.

Dernière modification par tumasgiu (Thu 16 May 2019 18:04)

Bonjour,

Merci pour votre réponse. Je n'ai pas eu de notification Georezo donc je ne l'avais pas vu jusqu'à aujourd'hui... désolé...
J'ai donc testé votre code.

DROP TABLE if exists blabla;
CREATE TABLE blabla AS
SELECT 
    st_collectionextract(
        st_split(rect.geom,
                coalesce((ST_CollectionHomogenize(st_collect(line.geom))), 
                    'LINESTRING EMPTY'::geometry('LINESTRING'))),
        3) AS geom
FROM
    polygones_zone_test as rect
LEFT JOIN
    public.lignes_zone_test as line
ON 
    st_intersects(rect.geom, line.geom)
GROUP BY
    rect.geom

J'ai ajouté des parenthèses dans la condition COALESCE puisque que le st_split avait trois conditions de géométrie ce qui ne lui plaisait pas, je ne sais pas si cette correction est bien juste mais il me semble que oui.
Il ne pose plus de problème sur le st_split, mais il me dit que le SRID n'est pas le même sur les deux couches (ERROR:  ERREUR:  Operation on mixed SRID geometries) alors que j'ai fait un st_srid et j'ai bien 2154 sur les deux couches... est-ce qu'une fonction préalable modifie le SRID?

Merci du temps consacré à me répondre !

Rebonjour,
Merci j'avais solutionné ce problème en faisant un setSRID sur le coalesce mais votre solution marche aussi.
Du coup la fonction marche, je conserve bien tous les polygones, et ils sont, d'apparence, bien découpés aux emplacements des lignes.
En revanche, j'ai toujours le même nombre de polygones que sur la couche initiale. Il me crée des polygones multi-part.
Je peux solutionner ça en éclatant ensuite mes polygones en single part mais peut-être y a t'il une solution sur le code d'origine.

Un grand merci !

Bonjour,
J'ai encore un petit message d'erreur que je n'avais pas sur ma petite zone de test...

ERROR:  ERREUR:  GEOSUnion: TopologyException: found non-noded intersection

Il s'agit d'endroits ou j'ai des lignes qui se croisent, ce qui est rare mais qui peut arriver.
J'ai essayé d'utiliser st_makevalid, d'après ce que j'ai pu trouver comme info sur les différents forums, en le plaçant à divers endroits du code, mais il me renvoie toujours la même erreur.
De ce que je comprend, il suffirait de créer des nœuds aux croisements des lignes pour avoir une topologie valide?

Oui, sans problème

Merci !

J'ai donc rajouté un st_snaptogrid sur la partie lignes de la fonction st_split. Je l'ai mis un peu arbitrairement je l'avoue, ne sachant pas trop ce que ça donnerait. Et ça marche bien. Le paramètre de 0.0001 est parfaitement arbitraire également.

DROP TABLE if exists blabla;
CREATE TABLE blabla AS
SELECT 
    st_collectionextract(
        st_split(rect.geom,
                st_snaptogrid(coalesce((ST_CollectionHomogenize(st_collect(st_makevalid(line.geom)))), 
                    'SRID=2154;LINESTRING EMPTY'::geometry('LINESTRING')),0.0001)),
        3) AS geom
FROM
    public.polygones as rect
LEFT JOIN
    public.lignes as line
ON 
    st_intersects(rect.geom, line.geom)
GROUP BY
    rect.geom