Message posté par : Nicolas Ribot ---------------------------------------- Alors voila une méthode qui marche mieux: Le principe est de travailler avec les arcs des deux couches (arc au sens topologique: linestrings connectées entre elles aux extrémités) et de reconstruire les polygones à partir de ces arcs: • on veut tous les arcs des provinces (car contours de référence) • et les arcs des communes internes aux provinces (ne suivant pas le contour des provinces) Pour obtenir ces arcs, st_linemerge, st_union et st_boundary sont utilisées (les boundaries d'un pg sont des linestrings): ----------------- Code : -- prepa tables ajout clé primaire (id plus facile alter table com add column id serial primary key ; alter table province add column id serial primary key ; select * from spatial_ref_sys where srid = 25830; -- ajout clés primaires alter table dts_commune add column id serial primary key ; alter table dts_province add column id serial primary key ; -- on tente de travailler avec les segments (arcs) connectés des commune: -- union et dump. select st_numgeometries(geom) from dts_province; select st_numgeometries(geom) from dts_commune; -- ok pg simple, on utilise st_geometryN(1) pour extraire le pg. drop table dts_com_arcs; create table dts_com_arcs as select (st_dump(st_lineMerge(st_union(st_boundary(st_geometryN(geom, 1)))))).geom from dts_commune; alter table dts_com_arcs add column id serial primary key ; -- et des provinces drop table dts_prov_arcs; create table dts_prov_arcs as select (st_dump(st_lineMerge(st_union(st_boundary(st_geometryN(geom, 1)))))).geom from dts_province; alter table dts_prov_arcs add column id serial primary key ; ----------------- Ensuite on trouve une méthode pour garder les arcs des communes internes aux provinces: par ex avec une distance moyenne des points de l'arc aux contours de la province: les arcs internes auront une grande distance moyenne par rapport aux arcs suivant le contour des provinces: ----------------- Code : -- on veut faire la reconstruction des arcs des provinces avec les arcs internes des communes -- il faut identifier les arcs internes des communes: -- plusieurs méthodes, par ex distance moyenne des points de l'arc au arcs des communes create table dts_com_arcs_pt as select id as idarc, (st_dumppoints(geom)).geom from dts_com_arcs; -- pk alter table dts_com_arcs_pt add column id serial primary key ; -- spatial index: create index on dts_com_arcs_pt using gist(geom); -- et sur les arcs province create index on dts_prov_arcs using gist(geom); vacuum analyse dts_prov_arcs; vacuum analyse dts_com_arcs_pt; -- methode empirique ici: on determine quelle dist moyenne permet de choisir les arcs internes: -- disons 62m :) drop table dts_com_internal_arc; create table dts_com_internal_arc as with tmp as ( select p.idarc, avg(t.dist) as avg_dist from dts_com_arcs_pt p cross join lateral ( select a.geom <-> p.geom as dist from dts_prov_arcs a order by a.geom <-> p.geom limit 1 ) as t group by p.idarc ) select tmp.*, a.geom from tmp join dts_com_arcs a on tmp.idarc = a.id where avg_dist > 62; ----------------- On corrige un peu la topologie en forcant la connexion des arcs internes des communes avec les arcs des provinces, pour etre sur que notre réseau soit fermé: les fn de réf linéaires marchent bien pour ca: on ajoute alors des points aux arcs internes, ce qui les connecte aux arcs des provinces. Il faut penser à cette étape à garder tous les arcs internes des communes: ceux prolongés et ceux qui ne le sont pas Une fois qu'on a les arcs internes prolongés et les arcs des provinces, on passe par st_polygonize pour reconstruire les polygones. A cette étape, les erreurs topo des communes peuvent engendrer des petits pg supplémentaires: pas grave, on fait une jointure spatiale entre les communes initiales et les centroid des communes corrigées, on retrouve alors nos communes: un st_union finale group by commune permet de nettoyer les petits polygones en erreurs: ----------------- Code : -- on corrige les erreurs topo éventuelles des arcs en prolongeant les arcs internes sur les contours des provinces, -- sauf si l'arc interne est connecté à un autre arc interne: -- on cherche l'arc province le plus proche de chaque start/end de nos arcs internes -- on reprend notre distance de 100m: si dist <, on snape ce point sur l'arc: -- ce snap est fait avec les fn de référencement linéaire, qui sont rapides et précises -- idx code le startpoint/endpoint: 0 pour start, -1 pour end: ca permettra de faire -- un addpoint avec cet idx. -- une fois les arcs internes prolongés, on peut construire le réseau de lignes -- arcs internes + arcs province et reconstruire les pg communes avec les arcs province -- on peut faire en 1 seule requete, ou détailler par étape, suivant le dataset -- attention a l'etape tmp1: on filtre certains arcs => il faut etre sur qu'on travaille avec tous les arcs -- internes avant la reconstruction des pg: drop table dts_newcom; create table dts_newcom as with tmp as ( select idarc, 0 as idx, st_startpoint(geom) as geom from dts_com_internal_arc UNION ALL select idarc, -1, st_endpoint(geom) as geom from dts_com_internal_arc ), tmp1 as ( select t.idarc, t.idx, t.geom, st_lineinterpolatepoint(b.geom, st_linelocatepoint(b.geom, t.geom)) as newpt from tmp t cross join lateral ( select a.id as idarcprov, t.geom <-> a.geom as dist, geom from dts_prov_arcs a order by t.geom <-> a.geom limit 1 ) as b where dist < 100 ), tmp2 as ( select st_addPoint(a.geom, t.newpt, t.idx) as geom from tmp1 t join dts_com_arcs a on t.idarc = a.id UNION ALL select geom from dts_prov_arcs -- on ajoute les éventuels arcs internes filtrés à l'étape précédente: UNION ALL select geom from dts_com_internal_arc a where not exists ( select null from tmp1 where a.idarc = tmp1.idarc ) ), tmp3 as ( select st_union(geom) as geom from tmp2 ) select (st_dump(st_polygonize(geom))).geom from tmp3; -- plus de pg que de communes: -- les erreurs topo entrainent la création de petits pg: facile a merger avec leur voisin -- pour les faire disparaitre: par ex en refaisant le lien spatial avec les communes et en faisant l'union des pg par commune: -- on peut faire un lien spatial (st_pointOnSurface pour retrouver les communes initiales: -- table finale des communes: select cod_mun, mun_cod_prov, nombre_com, c.id, st_union(n.geom) as geom from dts_commune c join dts_newcom n on st_contains(c.geom, st_pointOnSurface(n.geom)) group by 1, 2, 3, 4; ----------------- Et en image avec la piece jointe. Nico ---------------------------------------- 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=340426#p340426 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 ---------------------------------------- Bonjour, Oui j'ai bien le polygone découpé. J'ai du jouer avec la distance moyenne pour garder les arcs internes: 62m, puis 100m pour choisir les arcs province sur lesquels snapper les arcs internes. Si l'arc est là, mais que le polygone ne se crée pas, c'est effectivement que l'arc ne doit pas bien etre connecté aux autres arcs: vous pouvez tester la distance des start/end points de cet arc avec les autres arcs. J'ai fait mes tests sur PG 13, postgis 3.1, avec la derniere version de geos. Ca peut expliquer les différences, car cette nouvelle version est plus robuste topologiquement (ca vaut le coup d'essayer avec cette version chez vous). Pour raffiner le process, vous pouvez tester si les arcs internes touchent les autres arcs apres snapping (st_touches). si ce n'est pas le cas, vous pouvez prolonger artificiellement ces arcs pour forcer l'intersection avec les autres arcs: st_polygonize pourra alors reconstruire les pg. Nicolas ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=340481#p340481 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 : GEOFFROY DETRY ---------------------------------------- Merci Nicolas pour cette analyse. Je vais mettre à jour ma version de PostGIS/GEOS sur le serveur PG13 car je suis admin serveur. Par contre, ne contrôlant pas le PG10_12 où se trouve la BD de production, je vais sans doute devoir implémenter la prolongation des axes. Merci encore pour la méthodologie d'axes et les vérifications qui ont montré les évolutions de GEOS 3.10. Geoffroy ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=340499#p340499 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 ---------------------------------------- Elle est un peu pourrie ma méthode je pense :) Ce qui marcherait bien c'est de travailler avec les arcs des communes et provinces, et reconstruire le set des arcs provinces pour l'exterieur et arc communes internes aux provinces, puis de reconstruire les pg avec st_polygonize. .. a suivre Nicolas ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=340422#p340422 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 : GEOFFROY DETRY ---------------------------------------- Bonsoir, J'ai reproduit la méthode d'arcs décrites par Nicolas. Très clair et avec la garantie de respecter la géométrie de la couche référence. L'enjeu est de bien sélectionner les arcs internes. Pour généraliser la méthodologie à l'ensemble du dataset, je partirai sur un critère de sélection du style "arc communal interne = arc touchant le contour communal en maximum deux points. En effet je crains que la méthode basée sur la distance moyenne d'un arc communal à un arc provincial ne sélectionne pas tous les tronçon du dataset global. Je vais tester cela. Par contre, j'ai été attiré par un autre élément. Lors de la reconstruction des polygones à partir de la multiligne (via ST_Polygonize), je m'aperçois qu'une limite n'est pas prise en compte et donc qu'un polygone plus important est créé, au lieu d'être séparé en deux (cf. PJ). Pourtant, ladite limite est bien dans la polyligne retournée en "tmp3" et topologiquement tout semble correct. Dans le screenshot proposé par Nicolas, on retrouve bien les deux polygones. J'ai examiné ce tronçon pour comprendre ce qu'il se passait mais sans vraiment trouver une explication, si ce n'est que la séquence des vertex ferait en sorte que la polyligne n'est pas considérée comme fermée. @Nicolas, si vous avez l'occasion, pouvez-vous confirmer que chez vous, le polygone que j'identifie en rouge dans le screenshot est bien coupé en deux? Cela me permettrait de voir si il s'agit d'une limitation Merci Geoffroy ---------------------------------------- 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=340479#p340479 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 ---------------------------------------- Sur cette version: ----------------- Code : POSTGIS="3.0.3 6660953" [EXTENSION] PGSQL="120" GEOS="3.8.1-CAPI-1.13.3" SFCGAL="1.3.9" PROJ="7.2.0" GDAL="GDAL 3.2.0, released 2020/10/26" LIBXML="2.9.4" LIBJSON="0.15" LIBPROTOBUF="1.3.3" WAGYU="0.4.3 (Internal)" TOPOLOGY RASTER ----------------- J'obtiens le meme résultat que vous (pg manquants) ca semble lié à la version de GEOS: soit vous mettez à jour, soit vous prolongez artificiellement les extrémités des arcs internes (ex ici: https://lists.osgeo.org/pipermail/postgis-users/2013-April/036831.html) Nicolas ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=340498#p340498 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