{"id":10693,"date":"2023-06-20T09:00:00","date_gmt":"2023-06-20T07:00:00","guid":{"rendered":"https:\/\/www.sigterritoires.fr\/?p=10693"},"modified":"2023-06-22T10:28:25","modified_gmt":"2023-06-22T08:28:25","slug":"eliminer-les-superpositions-et-les-espaces-entre-polygones-dans-une-coucheavec-qgis-et-postgis","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/eliminer-les-superpositions-et-les-espaces-entre-polygones-dans-une-coucheavec-qgis-et-postgis\/","title":{"rendered":"\u00c9liminer les superpositions et les espaces entre polygones dans une couche(avec QGis et Postgis)"},"content":{"rendered":"\n

Si vous avez l’occasion de travailler avec les donn\u00e9es cadastrales correspondantes aux Lieux-dits <\/strong><\/em>vous constaterez que la coh\u00e9rence spatiale (topologie<\/a>) n’est pas respect\u00e9e:<\/p>\n\n\n\n

<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

En effet, chaque lieu-dit a comme limite le bord du chemin attenant, ce qui laisse des espaces vides l\u00e0 o\u00f9 les chemins se trouvent. Si dans certains cas ceci n’est pas un probl\u00e8me, il le devient quand vous ajoutez la couche \u00ab\u00a0Communes<\/em><\/strong>\u00a0\u00bb (du m\u00eame cadastre) \u00e0 votre carte:<\/p>\n\n\n\n\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

A cette \u00e9chelle cela ne fait pas tr\u00e8s propre, mais si on zoome c’est encore pire<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

Et si on consid\u00e8re la couche lieux-dits elle-m\u00eame, on trouve des anomalies topologiques nombreuses:<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

Nous allons voir ici comment r\u00e9soudre les deux probl\u00e8mes:<\/p>\n\n\n\n

    \n
  1. Comment corriger la topologie de la couche lieux-dits en supprimant les superpositions de polygones et les espaces vides entre eux;<\/li>\n\n\n\n
  2. Comment recr\u00e9er une couche communes coh\u00e9rente avec la nouvelle couche lieux-dits et se superposant exactement aux limites des lieux-dits.<\/li>\n<\/ol>\n\n\n\n

    L’exemple ci-apr\u00e8s a \u00e9t\u00e9 r\u00e9alis\u00e9 avec les donn\u00e9es du cadastre officiel du D\u00e9partement 71, t\u00e9l\u00e9charg\u00e9es \u00e0 partir du lien https:\/\/cadastre.data.gouv.fr\/data\/etalab-cadastre\/2023-01-01\/shp\/departements\/71\/<\/a><\/p>\n\n\n\n

    Comment corriger la topologie de la couche lieux-dits<\/h2>\n\n\n\n

    La solution propos\u00e9e ici n\u00e9cessite une base de donn\u00e9es PostgreSQL avec Postgis<\/a>. Il est possible de l’adapter \u00e0 SQLite pour utiliser un fichier geopackage<\/a> et le DB Manager de QGis. Mais il faut s’attendre \u00e0 avoir des temps de traitement tr\u00e8s longs et \u00e0 contourner quelques carences et bugs. Nous verrons ici la meilleure des solutions, mais nous pr\u00e9sentons la solution avec Geopackage dans un autre article.<\/p>\n\n\n\n

    1-Cr\u00e9er une couche avec un buffer autour des lieux-dits<\/h3>\n\n\n\n

    Pour g\u00e9rer les espaces entre polygones, m\u00eame si certains outils d’accrochage permettent de traiter certains probl\u00e8mes, il n’en reste pas moins que la solution la plus simple et radicale est de construire un buffer autour des polygones de mani\u00e8re \u00e0 transformer les espaces vides en superpositions. Dans le cas des lieux-dits, les espaces vides sont de l’ordre de 10m. Nous construisons un buffer de 6m autour des lieux-dits ce qui donnera une superposition de l’ordre de 2-3 m\u00e8tres minimum.<\/p>\n\n\n\n

    \"\"<\/a><\/figure>\n\n\n\n
    \"\"<\/a><\/figure>\n\n\n\n

    Choisissez un fichier en sortie permanent au format Geopackage, cela vous permettra de suivre la proc\u00e9dure de l’\u00e9tape suivante sans modification.<\/p>\n\n\n\n

    Le r\u00e9sultat de la couche tampon est le suivant:<\/p>\n\n\n\n

    \"\"<\/a><\/figure>\n\n\n\n
    \"\"<\/a><\/figure>\n\n\n\n

    Nous n’avons plus \u00e0 g\u00e9rer des espaces entre les polygones, nous avons seulement des superpositions.<\/p>\n\n\n\n

    2- Chargez la couche lieux-dits dans PostgreSQL<\/h3>\n\n\n\n

    Vous pouvez suivre le d\u00e9tail de comment charger la couche dans l’article Comment charger un geopackage dans Postgis avec QGis<\/a> . Vous utilisez la m\u00eame proc\u00e9dure pour charger des shapefiles.<\/p>\n\n\n\n

    \"\"<\/a><\/figure>\n\n\n\n

    Et vous retrouvez votre table dans PostreSQL avec pgAdmin<\/p>\n\n\n\n

    \"\"<\/a><\/figure>\n\n\n\n

    3- Ex\u00e9cuter le script SQL<\/h3>\n\n\n\n

    Nous allons utiliser un seul script pour faire l’ensemble du travail:<\/p>\n\n\n\n

      \n
    1. On cr\u00e9e une vue avec deux collections de g\u00e9om\u00e9tries. Au d\u00e9part ce sont les polygones des lieux-dits en double exemplaire<\/li>\n\n\n\n
    2. On garde dans une premi\u00e8re collection les parties des polygones qui ne se superposent pas avec d’autres<\/li>\n\n\n\n
    3. On garde dans la deuxi\u00e8me collection les parties superpos\u00e9es<\/li>\n\n\n\n
    4. On ajoute au plus petit des deux polygones concern\u00e9s par la superposition la partie superpos\u00e9e<\/li>\n<\/ol>\n\n\n\n

      Voici le script:<\/p>\n\n\n\n

      DROP TABLE IF EXISTS table_corrigee ; \n\nCREATE OR REPLACE VIEW auto_jointure AS SELECT\n\n    tampon.id as id,\n    tampon.nom as \"nom\",\n    tampon.commune as \"commune\",\n-- on r\u00e9cup\u00e8re tous les polygones de la table tampon dans deux collections : une o\u00f9 on aura les g\u00e9om\u00e9tries compl\u00e8tes et une autre o\u00f9 on aura les parties superpos\u00e9es. Pour l'instant les deux collections contiennent les g\u00e9om\u00e9tries compl\u00e8tes.\n    ST_Union(tampon.geom) AS full_geom,\n    ST_Union(tampon_bis.geom) AS shared_geom\n    FROM  tampon ,tampon AS tampon_bis\n    WHERE\n      --on v\u00e9rifie la validit\u00e9 des polygones comme pr\u00e9caution\n      ST_IsValid(tampon.geom) AND ST_IsValid(tampon_bis.geom)\n      --on filtre pour conserver les polygones qui s'intersectent\n      AND ST_intersects(tampon.geom,tampon_bis.geom)   \n      --on \u00e9limine les polygones en intersection avec eux m\u00eames\n      AND tampon.id <> tampon_bis.id          \n      --pour une intersection de 2 polygones, on garde que le plus petit\n      AND ST_Area(tampon_bis.geom) < ST_Area(tampon.geom)\n      --puisque l'on fait des \"unions\u201d, il faut effectuer un regroupement sur les autres attributs\n    GROUP BY tampon.id,tampon.\"nom\" , tampon.commune ;\n\/*\nOn va cr\u00e9er une table qui va contenir le r\u00e9sultat final : la premi\u00e8re \u00e9tape consiste \u00e0 enlever toutes les zones de superposition. La deuxi\u00e8me \u00e9tape consiste \u00e0 ajouter les zones de superposition au polygone le plus petit des deux. Les trous cr\u00e9\u00e9s par la premi\u00e8re \u00e9tape sont bouch\u00e9s par cette deuxi\u00e8me \u00e9tape. Les polygones finaux soint jointifs, sans superposition ni espace.\n*\/\n\n    CREATE TABLE table_corrigee AS SELECT\n    id,\n    \"nom\",\n    commune,\n\/* On soustrait les intersections \u00e0 la couche contenant tous les polygones; on supprime ainsi les parties en conflits*\/\n    ST_Multi(ST_Difference(full_geom,shared_geom)) as geom,\n    ST_Area(full_geom) as area\n    FROM auto_jointure\n    WHERE ST_IsValid(full_geom) AND ST_IsValid(shared_geom)\n\/*On doit ensuite rajouter les intersections que l'on vient de soustraire pour remplir les espaces cr\u00e9es.*\/\n    UNION \n      SELECT\n      id,\n      \"nom\",\n      commune,\n      geom,\n      ST_Area(geom)\n      FROM tampon\n      WHERE id NOT IN (SELECT id FROM auto_jointure);\n\/*On red\u00e9finit ensuite une cl\u00e9 primaire pour que QGIS puisse charger la couche.*\/\n\nALTER TABLE table_corrigee ADD CONSTRAINT pk_table_corrigee PRIMARY KEY (id);\n<\/code><\/pre>\n\n\n\n
      Pour adapter le script \u00e0 votre cas, vous devez remplacer avec un traitement de texte les \u00e9l\u00e9ments suivants:<\/p>\n\n\n\n