Message posté par : Nicolas Ribot ---------------------------------------- non pas reçu ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=351797#p351797 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : Nicolas Ribot ---------------------------------------- Il y a surement plus simple en passant par st_snap(geom, geom) Il faut identifier les geom à snaper entre elles (les plus proches), puis faire l'union des geom snappées (car il y aura immanquablement des trous) et virer les trous (st_dumpRings). avec un with recursive ... peut etre. Pas super facile comme problématique ;) Nicolas ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=351805#p351805 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : Sylvain M. ---------------------------------------- J'y avais pensé, mais je vais perdre la topologie des contours avec les multipolygones initiaux (ce que j'aimerais éviter). Mais je pense qu'on peut corriger ça ensuite avec le st_snap(geom, geom) que je découvre : qu'en penses-tu ? ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=351828#p351828 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : Sylvain M. ---------------------------------------- Je propose ici une réflexion de méthode : 1] Créer un polygone convexe du polygone (ST_ConvexHull) 2] Combler les trous de ce polygone convexe via des Polygones (ST_DumpRings je crois, mais je ne maitrise pas) 3] Pour chaque polygone de comblement, déterminer via un indice basé sur la forme, si il s'agit d'une surface à exclure de la fusion (les fameux trous à conserver), ou d'une emprise de route : caractérisée par une forme "linéaire" (mais je ne sais pas comment calculer cet indice) 4] fusion du Multipolygone source avec les petit polygones identifiés comme emprise routière ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=351800#p351800 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : Nicolas Ribot ---------------------------------------- Une premiere approche en se basant sur les points des polygones proches d'un autres polygone. On cherche ces points, on construit le buffer a x m de ces points, on fait l'union de ces buffer, puis l'union de ce buffer avec tous les polygones: on obtient un polygone avec pas mal de trous. On vire tous les trous de ce résultat, et on ajoute les trous présents dans les polygones initiaux. Pb: • certaines parties sont comblées: elles ne devraient pas l'être: il faut appliquer un critere de surface sur les trous qu'on vire. • ca ne prend pas en compte les routes qui coupent en partie un polygone (mais la requete peut etre adaptée pour trouver ces points et les combler avec un buffer. En SQL: ----------------- Code : -- la table frtdmp contient le dump du multipolygone à nettoyer sous forme de polygones simples create table frtdmp as select (st_dump(geom)).path[1] as id, (st_dump(geom)).geom from frt; -- la table frtpts contient le dump des points des polygones create table frtpts as select id as idpg, (st_dumppoints(st_exteriorring(geom))).path as path, (st_dumppoints(st_exteriorring(geom))).geom from frtdmp; -- index et tout... with holes as ( -- trous des polygones select st_dumprings(geom) as dmp from frtdmp ), array_holes as ( -- on regroupe les trous dans un tableau select array_agg(st_exteriorring((dmp).geom)) as holes from holes where (dmp).path[1] > 0 ),tmp as ( -- on fait l'union des buffers des points proches des contours select st_union(st_buffer(geom2, 30)) as geom from frtpts f cross join lateral ( select f2.id, f2.idpg as idpg2, f2.geom <-> f.geom as dist, f2.geom as geom2 from frtpts f2 where f.id < f2.id and f.idpg <> f2.idpg order by f.geom <-> f2.geom limit 1 ) as t ), tmp1 as ( -- on fabrique la liste des polygones initiaux + le polygone issu des buffers des points select geom from tmp UNION ALL select geom from frtdmp ) select st_makePolygon(st_exteriorring(st_union(geom)), holes) -- on refait un polygone en ajoutant les trous initiaux from tmp1, array_holes group by holes; ----------------- Image du résultat jointe Nicolas ---------------------------------------- Ce message est accompagné de fichiers, pour les télécharger, suivre le lien ci-dessous. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=351808#p351808 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : Nicolas Ribot ---------------------------------------- Oui effectivement, mais la topo est perdu puisqu'il manque une info au niveau des découpages par les routes. La topo peut etre reconstruite s'il le faut (grass, postgis topology) Oui snap ca peut marcher. Nicolas ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=351831#p351831 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net