{"id":10834,"date":"2023-06-22T09:00:00","date_gmt":"2023-06-22T07:00:00","guid":{"rendered":"https:\/\/www.sigterritoires.fr\/?p=10834"},"modified":"2023-06-18T07:44:39","modified_gmt":"2023-06-18T05:44:39","slug":"eliminacion-de-solapamientos-y-huecos-entre-poligonos-de-una-capa-con-qgis-y-geopackage","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/es\/eliminacion-de-solapamientos-y-huecos-entre-poligonos-de-una-capa-con-qgis-y-geopackage\/","title":{"rendered":"Eliminaci\u00f3n de solapamientos y huecos entre pol\u00edgonos de una capa (con QGis y Geopackage)"},"content":{"rendered":"\n

Aqu\u00ed retomamos el ejemplo tratado en el art\u00edculo Eliminaci\u00f3n de solapamientos y espacios entre pol\u00edgonos de una capa (con QGis y Postgis)<\/a><\/em><\/strong> para ver el mismo tipo de soluci\u00f3n cuando no se dispone de una base de datos PostgreSQL\/Postgis. Aqu\u00ed s\u00f3lo vamos a utilizar las posibilidades que ofrece QGis trabajando sobre un formato de datos Geopackage<\/a>.<\/p>\n\n\n\n\n\n\n\n

Si tiene la oportunidad de trabajar con los datos catastrales franceses correspondientes a Localidades, observar\u00e1 que no se respeta la coherencia espacial (topolog\u00eda):<\/p>\n\n\n\n

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

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

De hecho, el l\u00edmite de cada localidad es el borde de la carretera colindante, lo que deja espacios vac\u00edos donde est\u00e1n las carreteras. Aunque en algunos casos esto no es un problema, se convierte en uno cuando se a\u00f1ade la capa \u00ab\u00a0Comunas\u00a0\u00bb (del mismo catastro) al mapa:<\/p>\n\n\n\n

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

A esta escala no parece muy limpio, pero si se ampl\u00eda es a\u00fan peor.<\/p>\n\n\n\n

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

Y si nos fijamos en la propia capa de localidad, encontramos numerosas anomal\u00edas topol\u00f3gicas:<\/p>\n\n\n\n

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

Aqu\u00ed veremos c\u00f3mo resolver el primer problema: C\u00f3mo corregir la topolog\u00eda de la capa de localidad eliminando los pol\u00edgonos superpuestos y los espacios vac\u00edos entre ellos<\/strong>.<\/p>\n\n\n\n

Para el segundo problema: C\u00f3mo volver a crear una capa com\u00fan coherente con la nueva capa de localidad y superpuesta exactamente a los l\u00edmites de las localidade<\/strong>s, consulte el art\u00edculo citado anteriormente (Eliminaci\u00f3n de superposiciones y espacios entre pol\u00edgonos en una capa (con QGis y Postgis)<\/em><\/strong>), ya que en esta parte se utilizan simplemente QGis y Geopackage.<\/p>\n\n\n\n

El ejemplo siguiente se cre\u00f3 utilizando datos del catastro oficial del departamento 71, descargados del enlace https:\/\/cadastre.data.gouv.fr\/data\/etalab-cadastre\/2023-01-01\/shp\/departements\/7<\/a>1\/.<\/p>\n\n\n\n

C\u00f3mo corregir la topolog\u00eda de la capa de localidad<\/h2>\n\n\n\n

La soluci\u00f3n aqu\u00ed propuesta no requiere una base de datos PostgreSQL\/Postgis. Sin embargo, dadas las limitaciones de SQLite3 y un error en el DB Manager de QGis 3.30, habr\u00e1 algunas manipulaciones adicionales. Tambi\u00e9n debe esperar largos tiempos de procesamiento, por lo que no es adecuado para el procesamiento de grandes capas.<\/p>\n\n\n\n

1-Crear una capa con un buffer alrededor de las localidades<\/h3>\n\n\n\n

Para gestionar los espacios entre pol\u00edgonos, aunque se pueden utilizar algunas herramientas de ajuste para tratar ciertos problemas, la soluci\u00f3n m\u00e1s sencilla y radical es construir un b\u00fafer alrededor de los pol\u00edgonos para transformar los espacios vac\u00edos en superposiciones. En el caso de las localidades, los espacios vac\u00edos son del orden de 10 m. Construimos un b\u00fafer de 6 m alrededor de las localidades, lo que dar\u00e1 lugar a una superposici\u00f3n de unos 2-3 metros.<\/p>\n\n\n\n

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

Elija un archivo de salida permanente en formato Geopackage, que le permitir\u00e1 seguir directamente el procedimiento del paso siguiente sin modificaciones.<\/p>\n\n\n\n

El resultado de la capa intermedia es el siguiente:<\/p>\n\n\n\n

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

Ya no tenemos que gestionar los espacios entre pol\u00edgonos, s\u00f3lo tenemos superposiciones.<\/p>\n\n\n\n

2- Abrir el QGis DataBase Manager<\/h3>\n\n\n\n
\"\"<\/a><\/figure>\n\n\n\n

Aparecer\u00e1 la ventana DB Manager<\/p>\n\n\n\n

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

Si su archivo geopackage no aparece en la lista de archivos, haz clic con el bot\u00f3n derecho del rat\u00f3n en Geopackage->Nueva conexi\u00f3n<\/em><\/strong> y c\u00e1rguelo.<\/p>\n\n\n\n

Haz clic en la capa almacenada en el paso anterior, y luego en el Men\u00fa Base de Datos -> Ventana SQL<\/em><\/strong>.<\/p>\n\n\n\n

3- Ejecuta las l\u00edneas del script SQL<\/h3>\n\n\n\n

Todas las l\u00edneas del siguiente script corresponden a las siguientes tareas<\/p>\n\n\n\n

    \n
  1. Crear una vista con dos colecciones de geometr\u00edas. Inicialmente, se trata de los pol\u00edgonos de las localidades por duplicado.<\/li>\n\n\n\n
  2. La primera colecci\u00f3n contiene las partes de los pol\u00edgonos que no se solapan con otras<\/li>\n\n\n\n
  3. Las partes superpuestas se guardan en la segunda colecci\u00f3n<\/li>\n\n\n\n
  4. Por \u00faltimo, a\u00f1adimos la parte superpuesta al menor de los dos pol\u00edgonos afectados por la superposici\u00f3n.<\/li>\n<\/ol>\n\n\n\n

    A diferencia de la ventana SQL en pgAdmin y PostgreSQL\/Postgis, en la ventana SQL del DB Manager es necesario ejecutar el script consulta por consulta.<\/p>\n\n\n\n

    Para adaptar el script a tu caso, necesitas reemplazar los siguientes elementos utilizando un procesador de textos:<\/p>\n\n\n\n

      \n
    1. todas las apariciones de tampon <\/em><\/strong>por el nombre de la tabla que contiene las localidades con su buffer,<\/li>\n\n\n\n
    2. todas las apariciones de table_corrigee<\/em><\/strong> por el nombre de la tabla de resultados que desees,<\/li>\n\n\n\n
    3. y comprueba que el identificador de tu tabla de buffers es id<\/em><\/strong>. Por defecto puede llamarse fid<\/em><\/strong>. En este caso, sustituye todas las apariciones de id<\/em><\/strong> por fid<\/em><\/strong>.<\/li>\n<\/ol>\n\n\n\n

      Estas son las l\u00edneas del script:<\/p>\n\n\n\n

      Primera consulta (si ya ha ejecutado el script)<\/h3>\n\n\n\n
      DROP TABLE IF EXISTS table_corrigee ; 
      <\/code><\/pre>\n\n\n\n
      Segunda consulta (si ya ha ejecutado el script)<\/h3>\n\n\n\n
      DROP VIEW auto_jointure;<\/code><\/pre>\n\n\n\n
      Tercera consulta<\/h3>\n\n\n\n
      CREATE OR REPLACE VIEW auto_jointure AS SELECT\n\n    tampon.id as id,\n    tampon.nom as \"nom\",\n    tampon.commune as \"commune\",\n-- recuperar todos los pol\u00edgonos de la tabla buffer en dos colecciones: una que contiene las geometr\u00edas completas y otra que contiene las partes superpuestas. Por el momento, ambas colecciones contienen las geometr\u00edas completas.\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      --comprobar la validez del pol\u00edgono como precauci\u00f3n\n      ST_IsValid(tampon.geom) AND ST_IsValid(tampon_bis.geom)\n      --filtro para retener los pol\u00edgonos que se intersecan\n      AND ST_intersects(tampon.geom,tampon_bis.geom)   \n      --eliminar los pol\u00edgonos de intersecci\u00f3n con ellos mismos\n      AND tampon.id <> tampon_bis.id          \n      --para una intersecci\u00f3n de 2 pol\u00edgonos, conservar s\u00f3lo el m\u00e1s peque\u00f1o\n      AND ST_Area(tampon_bis.geom) < ST_Area(tampon.geom)\n      --ya que estamos haciendo 'uniones', necesitamos realizar una agrupaci\u00f3n en los otros atributos\n    GROUP BY tampon.id,tampon.\"nom\" , tampon.commune ;<\/em><\/code><\/pre>\n\n\n\n
      Cuarta consulta<\/h3>\n\n\n\n
      \/*\nVamos a crear una tabla que contendr\u00e1 el resultado final: el primer paso consiste en eliminar todas las \u00e1reas superpuestas. El segundo paso consiste en a\u00f1adir las \u00e1reas solapadas al menor de los dos pol\u00edgonos. Los huecos creados por el primer paso se rellenan con este segundo paso. Los pol\u00edgonos finales est\u00e1n unidos, sin solapamientos ni huecos.\n*\/\n\n    CREATE TABLE table_corrigee AS SELECT\n    id,\n    \"nom\",\n    commune,\n\/* Las intersecciones se restan de la capa que contiene todos los pol\u00edgonos; de este modo, se eliminan las partes conflictivas*\/\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\/*A continuaci\u00f3n, debemos a\u00f1adir las intersecciones que acabamos de restar para rellenar los espacios creados.*\/\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\n<\/code><\/pre>\n\n\n\n
      Soluci\u00f3n de los problemas de DB Managerr<\/h2>\n\n\n\n
      Despu\u00e9s de ejecutar las consultas, puede que se sorprenda al ver que ni el establo ni las vistas creadas durante los scripts aparecen en la lista de tablas de Geopackage.<\/p>\n\n\n\n
      \"\"<\/a><\/figure>\n\n\n\n
      Y sin embargo, si escribes una consulta referida a las tablas del script, funcionan. Las tablas est\u00e1n en alguna parte, pero DB Manager no las ve\u2026 Es un error, conocido y reconocido, y\u2026 \u00a1aceptado! Como los desarrolladores de QGis est\u00e1n preparando un sustituto para DB Manager, Data Source Manager, no tienen intenci\u00f3n de solucionar este problema.<\/p>\n\n\n\n
      En resumen, si no podemos resolver el problema, podemos sortearlo.<\/p>\n\n\n\n
      Con DB Manager vamos a crear una tabla, con la opci\u00f3n est\u00e1ndar<\/p>\n\n\n\n
      \"\"<\/a><\/figure>\n\n\n\n
      Y creamos una tabla con los mismos campos que la capa bufferizada<\/p>\n\n\n\n
      \"\"<\/a><\/figure>\n\n\n\n
      Que aparece en la lista de tablas del geopackagee<\/p>\n\n\n\n
      \"\"<\/a><\/figure>\n\n\n\n
      Ahora todo lo que tenemos que hacer es cargar el resultado de nuestro script en esta tabla \u00ab\u00a0oficial\u00a0\u00bb utilizando la siguiente consulta SQL:<\/p>\n\n\n\n
      insert into recup
      select fid,geom,nom,commune from lieux_dits_corr<\/code><\/pre>\n\n\n\n
      Cargamos esta nueva capa en QGis:<\/p>\n\n\n
      \n
      \"\"<\/a>

      A la izquierda la capa corregida, a la derecha la capa de localidades original<\/figcaption><\/figure>\n<\/div>\n\n\n
      Hemos resuelto el problema de la consistencia topol\u00f3gica de las localidades.<\/p>\n\n\n\n
      Para los puristas, una explicaci\u00f3n final del problema con la visualizaci\u00f3n de tablas en DB Manager.<\/p>\n\n\n\n
      El problema no viene de QGIS, sino del est\u00e1ndar de geopackage.
      OGR (y qgis) s\u00f3lo buscan tablas en la tabla gpkg_contents<\/em><\/strong>. Para que sean visibles, es necesario a\u00f1adir una nueva l\u00ednea a estas tablas:<\/p>\n\n\n\n
      INSERT INTO gpkg_contents (table_name, data_type) values (table_name,'attributes'). # attributes or features (if geometry); INSERT INTO gpkg_geometry (table_name, column_name, geometry_type_name, srs_id, z, m) (table_name, 'geom', 'Point', 4326, 1, 1);<\/code><\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n
      En este caso, sin embargo, cabr\u00eda esperar que QGIS DB Manager ejecutara estas consultas autom\u00e1ticamente.<\/p>\n","protected":false},"excerpt":{"rendered":"
      Aqu\u00ed retomamos el ejemplo tratado en el art\u00edculo Eliminaci\u00f3n de solapamientos y espacios entre pol\u00edgonos de una capa (con QGis y Postgis) para ver el mismo tipo de soluci\u00f3n cuando no se dispone de una base…<\/p>\n","protected":false},"author":1,"featured_media":10696,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"give_campaign_id":0,"_bbp_topic_count":0,"_bbp_reply_count":0,"_bbp_total_topic_count":0,"_bbp_total_reply_count":0,"_bbp_voice_count":0,"_bbp_anonymous_reply_count":0,"_bbp_topic_count_hidden":0,"_bbp_reply_count_hidden":0,"_bbp_forum_subforum_count":0,"sfsi_plus_gutenberg_text_before_share":"","sfsi_plus_gutenberg_show_text_before_share":"","sfsi_plus_gutenberg_icon_type":"","sfsi_plus_gutenberg_icon_alignemt":"","sfsi_plus_gutenburg_max_per_row":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[2164,1340,1357],"tags":[2166,2168,1359,1989,2170,2174,2172,2180,1369,2176,2182,2178],"class_list":["post-10834","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-geopackage","category-postgresql-es","category-qgis-es","tag-analisis-espacial","tag-datos-espaciales","tag-geometria","tag-geopackage-es","tag-gestion-territorial","tag-investigacion-geografica","tag-modelizacion-geoespacial","tag-poligonos-es","tag-qgis-es","tag-sig-sistema-de-informacion-geografica","tag-territorios-es","tag-visualizacion-cartografica"],"aioseo_notices":[],"jetpack_featured_media_url":"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2023\/06\/a2.png?fit=3013%2C1793&ssl=1","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-2OK","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/10834","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/comments?post=10834"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/10834\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media\/10696"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=10834"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=10834"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=10834"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}