{"id":10835,"date":"2023-06-22T09:00:00","date_gmt":"2023-06-22T07:00:00","guid":{"rendered":"https:\/\/www.sigterritoires.fr\/?p=10835"},"modified":"2023-06-18T07:44:10","modified_gmt":"2023-06-18T05:44:10","slug":"eliminacao-de-sobreposicoes-e-lacunas-entre-poligonos-numa-camada-com-qgis-e-geopackage","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/pt\/eliminacao-de-sobreposicoes-e-lacunas-entre-poligonos-numa-camada-com-qgis-e-geopackage\/","title":{"rendered":"Elimina\u00e7\u00e3o de sobreposi\u00e7\u00f5es e lacunas entre pol\u00edgonos numa camada (com QGis e Geopackage)"},"content":{"rendered":"\n

Aqui retomamos o exemplo abordado no artigo Elimina\u00e7\u00e3o de sobreposi\u00e7\u00f5es e espa\u00e7os entre pol\u00edgonos numa camada (com QGis e Postgis) <\/a><\/em><\/strong>para ver o mesmo tipo de solu\u00e7\u00e3o quando n\u00e3o se tem uma base de dados PostgreSQL\/Postgis. Aqui s\u00f3 vamos utilizar as possibilidades oferecidas pelo QGis ao trabalhar com um formato de dados Geopackage<\/a>.<\/p>\n\n\n\n\n\n\n\n

Se tiver a oportunidade de trabalhar com os dados cadastrais franceses correspondentes aos Localidades, verificar\u00e1 que a consist\u00eancia espacial (topologia) n\u00e3o \u00e9 respeitada:<\/p>\n\n\n\n

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

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

De facto, o limite de cada localidade \u00e9 o limite da estrada adjacente, o que deixa espa\u00e7os vazios onde se encontram as estradas. Embora em alguns casos isto n\u00e3o seja um problema, torna-se um problema quando se adiciona a camada \u00ab\u00a0Comunas\u00a0\u00bb (do mesmo cadastro) ao mapa:<\/p>\n\n\n\n

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

A esta escala n\u00e3o parece muito limpo, mas se fizermos um zoom \u00e9 ainda pior.<\/p>\n\n\n\n

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

E se olharmos para a camada de localidade propriamente dita, encontramos numerosas anomalias topol\u00f3gicas:<\/p>\n\n\n\n

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

Aqui vamos ver como resolver o primeiro problema: Como corrigir a topologia da camada de localidade, removendo pol\u00edgonos sobrepostos e espa\u00e7os vazios entre eles<\/strong>.<\/p>\n\n\n\n

Para o segundo problema: Como recriar uma camada comum consistente com a nova camada de localidades e sobreposta exatamente aos limites das localidades<\/strong>, consulte o artigo citado acima (Elimina\u00e7\u00e3o de sobreposi\u00e7\u00f5es e espa\u00e7os entre pol\u00edgonos numa camada (com QGis e Postgis)<\/em><\/strong>, uma vez que esta parte utiliza simplesmente o QGis e o Geopackage.<\/p>\n\n\n\n

O exemplo abaixo foi criado utilizando os dados do registo predial oficial do departamento 71, descarregado a partir do link https:\/\/cadastre.data.gouv.fr\/data\/etalab-cadastre\/2023-01-01\/shp\/departements\/71\/<\/a><\/p>\n\n\n\n

Como corrigir a topologia da camada de localidade<\/h2>\n\n\n\n

A solu\u00e7\u00e3o aqui proposta n\u00e3o requer uma base de dados PostgreSQL\/Postgis. No entanto, dadas as limita\u00e7\u00f5es do SQLite3 e um bug no DB Manager do QGis 3.30, haver\u00e1 algumas manipula\u00e7\u00f5es adicionais. Deve tamb\u00e9m contar com tempos de processamento longos, raz\u00e3o pela qual n\u00e3o \u00e9 adequado para o processamento de grandes camadas.<\/p>\n\n\n\n

1-Criar uma camada com um buffer \u00e0 volta das localidades<\/h3>\n\n\n\n

Para gerir os espa\u00e7os entre pol\u00edgonos, mesmo que certas ferramentas de encaixe possam ser utilizadas para lidar com certos problemas, a solu\u00e7\u00e3o mais simples e mais radical \u00e9 construir um buffer \u00e0 volta dos pol\u00edgonos de modo a transformar os espa\u00e7os vazios em sobreposi\u00e7\u00f5es. No caso das localidades, os espa\u00e7os vazios s\u00e3o da ordem dos 10m. Constru\u00edmos um buffer de 6m \u00e0 volta das localidades, o que dar\u00e1 uma sobreposi\u00e7\u00e3o de cerca de 2-3 metros.<\/p>\n\n\n\n

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

Escolha um ficheiro de sa\u00edda permanente no formato Geopackage, que lhe permitir\u00e1 seguir o procedimento no passo seguinte diretamente sem modifica\u00e7\u00f5es.<\/p>\n\n\n\n

O resultado da camada tamp\u00e3o \u00e9 o seguinte:<\/p>\n\n\n\n

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

J\u00e1 n\u00e3o temos de gerir os espa\u00e7os entre os pol\u00edgonos, apenas temos sobreposi\u00e7\u00f5es.<\/p>\n\n\n\n

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

Aparecer\u00e1 a janela do DB Manager<\/p>\n\n\n\n

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

Se o seu ficheiro geopackage n\u00e3o aparecer na lista de ficheiros, clique com o bot\u00e3o direito do rato em Geopackage->New connection<\/em><\/strong> e carregue-o.<\/p>\n\n\n\n

Clique na camada armazenada em buffer na etapa anterior e, em seguida, no menu Database (Base de dados) -> SQL Window (Janela SQL).<\/em><\/strong><\/p>\n\n\n\n

3- Execute as linhas do gui\u00e3o SQL<\/h3>\n\n\n\n

Todas as linhas do script seguinte correspondem \u00e0s seguintes tarefas<\/p>\n\n\n\n

    \n
  1. Criar uma vista com duas colec\u00e7\u00f5es de geometrias. Inicialmente, estas s\u00e3o os pol\u00edgonos das localidades em duplicado.<\/li>\n\n\n\n
  2. A primeira cole\u00e7\u00e3o cont\u00e9m as partes dos pol\u00edgonos que n\u00e3o se sobrep\u00f5em a outras<\/li>\n\n\n\n
  3. As partes que se sobrep\u00f5em s\u00e3o mantidas na segunda cole\u00e7\u00e3o<\/li>\n\n\n\n
  4. Por fim, adicionamos a parte sobreposta ao mais pequeno dos dois pol\u00edgonos afectados pela sobreposi\u00e7\u00e3o.<\/li>\n<\/ol>\n\n\n\n

    Ao contr\u00e1rio da janela SQL no pgAdmin e PostgreSQL\/Postgis, na janela SQL do DB Manager \u00e9 necess\u00e1rio executar o script consulta por consulta.<\/p>\n\n\n\n

    Para adaptar o script ao seu caso, \u00e9 necess\u00e1rio substituir os seguintes elementos utilizando um processador de texto:<\/p>\n\n\n\n

      \n
    1. todas as ocorr\u00eancias de tampon <\/em><\/strong>pelo nome da sua tabela que cont\u00e9m os localidades em buffer,<\/li>\n\n\n\n
    2. todas as ocorr\u00eancias de table_corrigee<\/em><\/strong> pelo nome da tabela de resultados que pretende,<\/li>\n\n\n\n
    3. e verifique se o identificador da sua tabela armazenada em buffer \u00e9 id<\/em><\/strong>. Por defeito, pode ter o nome fid<\/em><\/strong>. Neste caso, substitua todas as ocorr\u00eancias de id<\/em><\/strong> por fid<\/em><\/strong>.<\/li>\n<\/ol>\n\n\n\n

      Aqui est\u00e3o as linhas do script:<\/p>\n\n\n\n

      Primeira consulta (se j\u00e1 tiver executado o script)<\/h3>\n\n\n\n
      DROP TABLE IF EXISTS table_corrigee ; 
      <\/code><\/pre>\n\n\n\n
      Segunda consulta (se j\u00e1 tiver executado o script)<\/h3>\n\n\n\n
      DROP VIEW auto_jointure;<\/code><\/pre>\n\n\n\n
      Terceira 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-- recupera todos os pol\u00edgonos da tabela buffer em duas colec\u00e7\u00f5es: uma contendo as geometrias completas e outra contendo as partes sobrepostas. De momento, ambas as colec\u00e7\u00f5es cont\u00eam as geometrias 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      --verificar a validade do pol\u00edgono como precau\u00e7\u00e3o\n      ST_IsValid(tampon.geom) AND ST_IsValid(tampon_bis.geom)\n      --filtro para reter os pol\u00edgonos que se intersectam\n      AND ST_intersects(tampon.geom,tampon_bis.geom)   \n      --eliminar os pol\u00edgonos que se auto-intersectam\n      AND tampon.id <> tampon_bis.id          \n      --para uma intersec\u00e7\u00e3o de 2 pol\u00edgonos, manter apenas o mais pequeno\n      AND ST_Area(tampon_bis.geom) < ST_Area(tampon.geom)\n      --Como estamos a fazer \"uni\u00f5es\", temos de efetuar um agrupamento nos outros atributos\n    GROUP BY tampon.id,tampon.\"nom\" , tampon.commune ;<\/em><\/code><\/pre>\n\n\n\n
      Quarta consulta<\/h3>\n\n\n\n
      \/*\nVamos criar uma tabela que conter\u00e1 o resultado final: o primeiro passo \u00e9 remover todas as \u00e1reas sobrepostas. O segundo passo consiste em adicionar as \u00e1reas sobrepostas ao menor dos dois pol\u00edgonos. Os buracos criados pelo primeiro passo s\u00e3o preenchidos por este segundo passo. Os pol\u00edgonos finais est\u00e3o unidos, sem sobreposi\u00e7\u00f5es ou espa\u00e7os vazios.\n*\/\n\n    CREATE TABLE table_corrigee AS SELECT\n    id,\n    \"nom\",\n    commune,\n\/* As intersec\u00e7\u00f5es s\u00e3o subtra\u00eddas da camada que cont\u00e9m todos os pol\u00edgonos; as partes em conflito s\u00e3o assim eliminadas*\/\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\/*Devemos ent\u00e3o adicionar as intersec\u00e7\u00f5es que acab\u00e1mos de subtrair para preencher os espa\u00e7os criados.*\/\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);<\/em>\n\n<\/code><\/pre>\n\n\n\n
      Resolu\u00e7\u00e3o de problemas do DB Manager<\/h2>\n\n\n\n
      Depois de executar as consultas, poder\u00e1 ficar surpreendido ao ver que nem as vistas est\u00e1veis nem as vistas criadas durante os scripts aparecem na lista de tabelas do Geopackage..<\/p>\n\n\n\n
      \"\"<\/a><\/figure>\n\n\n\n
      No entanto, se escrever uma consulta que fa\u00e7a refer\u00eancia \u00e0s tabelas de script, elas funcionam. As tabelas est\u00e3o algures, mas o DB Manager n\u00e3o as v\u00ea\u2026 \u00c9 um erro, conhecido e reconhecido, e\u2026 aceite! Como os criadores do QGis est\u00e3o a preparar um substituto para o DB Manager, o Data Source Manager, n\u00e3o t\u00eam inten\u00e7\u00e3o de resolver este problema.<\/p>\n\n\n\n
      Em suma, se n\u00e3o podemos resolver o problema, podemos contorn\u00e1-lo.<\/p>\n\n\n\n
      Com o DB Manager, vamos criar uma tabela, com a op\u00e7\u00e3o padr\u00e3o<\/p>\n\n\n\n
      \"\"<\/a><\/figure>\n\n\n\n
      E criamos uma tabela com os mesmos campos que a camada buffer<\/p>\n\n\n\n
      \"\"<\/a><\/figure>\n\n\n\n
      Aparece na lista de tabelas do geopackage<\/p>\n\n\n\n
      \"\"<\/a><\/figure>\n\n\n\n
      Agora, tudo o que temos de fazer \u00e9 carregar o resultado do nosso script nesta tabela \u00ab\u00a0oficial\u00a0\u00bb utilizando a seguinte 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
      Carregamos esta nova camada no QGis:<\/p>\n\n\n
      \n
      \"\"<\/a>

      \u00c0 esquerda a camada corrigida, \u00e0 direita a camada de localidade original<\/figcaption><\/figure>\n<\/div>\n\n\n
      Resolvemos o problema da consist\u00eancia topol\u00f3gica das localidades.<\/p>\n\n\n\n
      Para os mais puristas, uma explica\u00e7\u00e3o final sobre o problema da visualiza\u00e7\u00e3o das tabelas no DB Manager.<\/p>\n\n\n\n
      O problema n\u00e3o vem do QGIS, mas do padr\u00e3o de geopackage.
      O OGR (e o qgis) apenas procuram tabelas na tabela gpkg_contents<\/em><\/strong>. Para que sejam vis\u00edveis, \u00e9 necess\u00e1rio acrescentar uma nova linha a essas tabelas:<\/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
      Neste caso, no entanto, \u00e9 de esperar que o QGIS DB Manager execute estas consultas automaticamente.<\/p>\n","protected":false},"excerpt":{"rendered":"
      Aqui retomamos o exemplo abordado no artigo Elimina\u00e7\u00e3o de sobreposi\u00e7\u00f5es e espa\u00e7os entre pol\u00edgonos numa camada (com QGis e Postgis) para ver o mesmo tipo de solu\u00e7\u00e3o quando n\u00e3o se tem uma base de dados PostgreSQL\/Postgis.…<\/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_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":"","jetpack_post_was_ever_published":false},"categories":[2164,2134],"tags":[2138,2140,2156,2162,2142,2148,2144,2146,2160,2150,2152,2154],"class_list":["post-10835","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-geopackage","category-qgis","tag-analise-espacial","tag-dados-espaciais","tag-geometria-pt","tag-geopackage-pt","tag-gestao-territorial","tag-investigacao-geografica","tag-modelacao-geoespacial","tag-poligonos","tag-qgis-pt","tag-sig-sistema-de-informacao-geografica","tag-territorios","tag-visualizacao-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-2OL","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/10835","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=10835"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/10835\/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=10835"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=10835"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=10835"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}