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

Nous reprenons ici l’exemple trait\u00e9 dans l’article \u00c9liminer les superpositions et les espaces entre polygones dans une couche(avec QGis et Postgis)<\/em><\/strong> <\/a>pour voir le m\u00eame type de solution quand vous ne disposez pas d’une base de donn\u00e9es PostgreSQL\/Postgis. Nous allons utiliser ici seulement les possibilit\u00e9s offertes par QGis en travaillant sur un format de donn\u00e9es Geopackage<\/a>.<\/p>\n\n\n\n\n\n\n\n

Si vous avez l’occasion de travailler avec les donn\u00e9es cadastrales correspondantes aux Lieux-dits 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\u00a0\u00bb (du m\u00eame cadastre) \u00e0 votre carte:<\/p>\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 le premier probl\u00e8me : Comment corriger la topologie de la couche lieux-dits en supprimant les superpositions de polygones et les espaces vides entre eux<\/strong>.<\/p>\n\n\n\n

Pour le deuxi\u00e8me probl\u00e8me : Comment recr\u00e9er une couche communes coh\u00e9rente avec la nouvelle couche lieux-dits et se superposant exactement aux limites des lieux-dits<\/strong>, r\u00e9f\u00e9rez-vous \u00e0 l’article cit\u00e9 plus haut (\u00c9liminer les superpositions et les espaces entre polygones dans une couche(avec QGis et Postgis<\/a>)<\/em><\/strong>, car cette partie utilise simplement QGis et Geopackage.<\/p>\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 ne n\u00e9cessite pas d’une base de donn\u00e9es PostgreSQL\/Postgis. Par contre, \u00e9tant donn\u00e9 quelques limitations de SQLite3 et un bug du DB Manager de QGis 3.30,on aura quelques manipulations suppl\u00e9mentaires. Il faut aussi s’attendre \u00e0 avoir des temps de traitement longs, raison pour laquelle il n’est pas adapt\u00e9 au traitement de couches trop importantes. <\/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.<\/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 directement 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- Ouvrez le DataBase Manager de QGis<\/h3>\n\n\n\n
\"\"<\/a><\/figure>\n\n\n\n

Vous aurez la fen\u00eatre du DB Manager<\/p>\n\n\n\n

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

Si votre fichier geopackage n’appara\u00eet pas dans les fichiers list\u00e9s, cliquez-droit sur Geopackage->Nouvelle connexion<\/em><\/strong> et chargez -le.<\/p>\n\n\n\n

Cliquez sur la couche bufferis\u00e9e dans l’\u00e9tape pr\u00e9cedente, puis sur le Menu Base de donn\u00e9es ->Fen\u00eatre SQL<\/em><\/strong>.<\/p>\n\n\n\n

3- Ex\u00e9cuter les lignes du script SQL<\/h3>\n\n\n\n

L’ensemble des lignes du script suivant correspondent aux t\u00e2ches suivantes<\/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. Pour finir, on ajoute au plus petit des deux polygones concern\u00e9s par la superposition la partie superpos\u00e9e<\/li>\n<\/ol>\n\n\n\n

    Contrairement \u00e0 la fen\u00eatre SQL de pgAdmin et PostgreSQL\/Postgis, dans la fen\u00eatre SQL de DB Manager vous devez ex\u00e9cuter requ\u00eate par requ\u00eate. <\/p>\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