Le portail francophone de la géomatique

Salut,

j'ai du mal à concevoir le problème, vous pourriez poster
quelques copies d'écran ?

Merci.

Je pense que le problème réside dans la topologie de votre couche initiale.

Les scories ont tout l'air d'être des interstices entre polygones qui ne sont pas jointifs.
Plusieurs solutions sont possibles:

Utiliser st_snaptogrid sur votre couche de départ afin d'arrondir les coordonnées
des sommets de vos géométries sur une grille régulière.
C'est au prix d'une petite perte de précision.
Cela ne garantie pas que toutes les interstices seront comblées mais c'est peut être un
bon début.
Vous pouvez mettre à jour votre couche ou bien juste transformer votre couche
à la volée dans votre script, si vous voulez conserver la donnée originale.

Transformer votre couche géométrique en couche topologique en utilisant l'extension topologique
de Postgis avec la précision adaptée. La aussi perte de précision possible.

Utiliser le correcteur topologique de QGIS comme décrit dans le chapitre 5.3.1 du rapport dont vous
avez posté le lien (ou d'un autre logiciel..)

Sinon vous pouvez rester sur votre méthode, en conservant en plus
les anneaux intérieurs qui ne sont pas des scories.
Un moyen de distinguer erreur de topologie de réel anneau,
au vu de votre capture d'écran, serait de mesurer l'aire des anneaux.

SELECT
st_buildarea(
    st_collect(
         st_exteriorRing(g) 
        || 
        (array_agg(r) FILTER (WHERE st_area(st_buildArea(r)) > a_definir)
        )
     )
)
FROM (
    SELECT g, st_interiorRing(g, id) r
    FROM couche LATERAL JOIN st_numInteriorRing(g) id
    ) foo 
GROUP BY g
UNION
SELECT g FROM couche WHERE st_numInteriorRing(g) = 0

Dernière modification par tumasgiu (Wed 08 August 2018 22:51)

Bonjour Lucas,

Pour travailler sur des "scories", tu peux essayer d'extraire les trous de tes polygones, de les filtrer pour ensuite recréer les bons trous.


Voici une exemple du type de traitement qu'il faudrait faire qu'il faudrait généraliser pour pouvoir traiter une couche complète.

with
--Création d'une géométrie exemple
une_geom as  (
SELECT ST_PolyFromText('POLYGON((0 0,20 0,20 20,0 20,0 0),(10 10,10 11,11 11,11 10,10 10),(5 5,5 7,7 7,7 5,5 5))') AS geom
)

--Récupération du contour via le st_dumprings
-- La valeur path = 0 correspond au contour
ext as (
select row_number() over () as id,geom from (
SELECT (ST_DumpRings(a.geom)).geom,(ST_DumpRings(a.geom)).path[1] as p
          FROM une_geom
) as b where p=0
)

--Récuparation des trous via le st_dumprings
-- La valeur de path différente de 0 indique un trou
,trou as (
select st_union(geom) as geom from (
SELECT (ST_DumpRings(a.geom)).geom,(ST_DumpRings(a.geom)).path[1] as p
          FROM une geom
) as b where p!=0
--ici tu pourrais ajouter une condition pour filtrer sur la taille des trous
)

--Enfin, tu fais la différence entre le contour et les trous ce qui te donneras un polygon troué
select st_difference(ext.geom,trou.geom) from ext,trou


Alors attention, dans cet exemple simple, je me  base sur un polygone simple... Dans le cas où tu aurais à travailler avec des multipolygones, il faudrait passer par un st_dump (simple .. ) pour avoir chaque polygone pour ensuite en sortir les trous

Autre attention, des fois les "scories" proviennent d'autres problèmes de géométrie lié aux traitements. Aussi, il peut être intéressant  de faire une table  de travail contenant tes contours et tes trous, nettoyer les trous que tu ne veux pas conserver puis refaire les trous dans les polygones.

A+ et bonne journée.

Christophe BREDEL

Bonjour à vous,

Je réponds un peu en retard, de retour de congés. Merci pour vos multiples réponses, cela m'a donné plusieurs pistes à explorer.
Concernant les traitements des trous, j'ai essayé les commandes qui marchent plutôt bien avec des ajustements, mais il y a trop de risques de déperdition d'information. En effet, j'ai des anneaux de taille variables, qui peuvent être petits, et j'ai surtout beaucoup trop d'enregistrements pour pouvoir vérifier manuellement tous les traitements.

J'ai donc plutôt tenté de "d'arrondir les coordonnées des sommets de vos géométries sur une grille régulière" avec st_snaptogrid comme l'a suggéré Thomasgiu.

update matable
set geom = st_collectionExtract(st_makeValid(st_snapToGrid(geom,1)),3);

J'ai ensuite essayé de "topologiser" la couche comme détaillé dans les liens suivants :

https://goutergeo.hypotheses.org/173
https://makina-corpus.com/blog/metier/2 … voirie#id3

Malheureusement, j'ai une erreur de type "ERREUR:  Spatial exception - geometry intersects edge" et je ne peux aller au bout de la procédure (j'avoue ne pas avoir pu appliquer le script pl/pgsql recommandé, faute de pouvoir définir le chemin du répertoire d'extensions pour installer la bibliothèque)

Au final, j'ai fait le test des traitements juste après le snaptogrid et les scories disparaissent quasiment toutes, sans deperdition d'information en terme de nombres de polygones et de surfaces.

Après effectivement GRASS à l'air d'être quand même la solution opensource la plus adaptée pour la correction des géometries et topologies...
Merci à vous pour vos commentaires et votre aide, ça permet de mieux orienter les recherches quand on pense avoir fait le tour!

Dernière modification par Rlucas (Tue 21 August 2018 20:29)