![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2023\/06\/a1.png?ssl=1\")
<\/a><\/figure>\n\n\n\nDe 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\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 E se olharmos para a camada de localidade propriamente dita, encontramos numerosas anomalias topol\u00f3gicas:<\/p>\n\n\n\n Aqui, veremos como resolver ambos os problemas:<\/p>\n\n\n\n O exemplo abaixo foi criado com base nos dados do registo cadastral oficial do departamento 71, descarregado a partir da seguinte liga\u00e7\u00e3o https:\/\/cadastre.data.gouv.fr\/data\/etalab-cadastre\/2023-01-01\/shp\/departements\/71\/<\/a><\/p>\n\n\n\n A solu\u00e7\u00e3o aqui proposta necessita de uma base de dados PostgreSQL\/Postgis. \u00c9 poss\u00edvel adapt\u00e1-la ao SQLite para utilizar um ficheiro geopackage<\/a> e o QGis DB Manager. Mas \u00e9 preciso contar com tempos de processamento muito longos e contornar algumas defici\u00eancias e bugs. Vamos ver a melhor solu\u00e7\u00e3o aqui, mas apresentaremos a solu\u00e7\u00e3o com o Geopackage noutro artigo.<\/p>\n\n\n\n Para gerir os espa\u00e7os entre pol\u00edgonos, ainda que algumas ferramentas de snapping possam resolver 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 pelo menos 2-3 metros.<\/p>\n\n\n\n Escolha um ficheiro de sa\u00edda permanente no formato Geopackage, que lhe permitir\u00e1 seguir o procedimento no passo seguinte sem modifica\u00e7\u00f5es.<\/p>\n\n\n\n O resultado da camada tamp\u00e3o \u00e9 o seguinte:<\/p>\n\n\n\n J\u00e1 n\u00e3o temos de gerir os espa\u00e7os entre os pol\u00edgonos, apenas as sobreposi\u00e7\u00f5es.<\/p>\n\n\n\n Pode seguir os pormenores de como carregar a camada no artigo Como carregar um geopackage no Postgis com o QGis<\/a> . Utiliza-se o mesmo procedimento para carregar shapefiles.<\/p>\n\n\n\n E encontra a sua tabela no PostreSQL com o pgAdmin<\/p>\n\n\n\n Vamos utilizar um \u00fanico script para fazer todo o trabalho:<\/p>\n\n\n\n O gui\u00e3o \u00e9 o seguinte:<\/p>\n\n\n\n Para adaptar o gui\u00e3o ao seu caso, \u00e9 necess\u00e1rio utilizar um processador de texto para substituir os seguintes elementos<\/p>\n\n\n\n Carregue esta nova camada no QGis:<\/p>\n\n\n\n Resolvemos o problema da coer\u00eancia topol\u00f3gica das localidades.<\/p>\n\n\n\n Se sobrepusermos a nossa nova camada com os limites das comunas do cadastro, podemos ver que n\u00e3o s\u00e3o coerentes com as nossas novas localidades, tal como n\u00e3o s\u00e3o coerentes com as localidades originais.<\/p>\n\n\n\n Para obter fronteiras comunais coerentes, bastaria agrupar as localidades por comuna. O problema \u00e9 que nem todas as \u00e1reas comunais est\u00e3o divididas em localidades. Por conseguinte, existem \u00ab\u00a0lacunas\u00a0\u00bb na camada de localidade, o que \u00e9 perfeitamente normal.<\/p>\n\n\n\n Para construir comunas completas, \u00e9 necess\u00e1rio primeiro criar os pol\u00edgonos em falta para as localidades, atribuir o c\u00f3digo de comuna correspondente e, em seguida, agrupar as localidades utilizando o atributo de comuna.<\/p>\n\n\n\n Em primeiro lugar, vamos criar uma nova camada de pol\u00edgonos, direito de passagem, com um \u00fanico pol\u00edgono que envolve completamente a nossa camada de localidade.<\/p>\n\n\n\n De seguida, utilizamos o tratamento \u00ab\u00a0Diferen\u00e7a\u00a0\u00bb para manter na nossa tabela corrigida apenas as zonas de direito de passagem que n\u00e3o correspondem \u00e0s localidades existentes:<\/p>\n\n\n\n Pode deixar o resultado dos tratamentos como camadas tempor\u00e1rias porque n\u00e3o ter\u00e1 mais nenhuma utiliza\u00e7\u00e3o para elas depois do resultado final.<\/p>\n\n\n\n O resultado da diferen\u00e7a \u00e9:<\/p>\n\n\n\n O resultado \u00e9 um \u00fanico multipol\u00edgono. Antes de continuar, temos de converter esta camada multipol\u00edgono em pol\u00edgonos simples, para podermos remover o contorno da \u00e1rea.<\/p>\n\n\n\n Utilizamos o processamento de pe\u00e7a \u00fanica e pe\u00e7a m\u00faltipla<\/p>\n\n\n\n Quando o processamento estiver conclu\u00eddo, passe o resultado para o modo de edi\u00e7\u00e3o e seleccione o pol\u00edgono exterior.<\/p>\n\n\n\n E estamos a acabar com ele:<\/p>\n\n\n\n Assim, temos agora uma camada de Geometrias Simples<\/em><\/strong> que cont\u00e9m todos os locais onde n\u00e3o existem localidades e uma camada de localidades corrigidas<\/em><\/strong>. Fundimos as duas camadas numa \u00fanica camada sem lacunas:<\/p>\n\n\n\n Ainda temos de resolver o problema dos atributos dos pol\u00edgonos que n\u00e3o s\u00e3o localidades. Estes precisam de ter um atributo que contenha o c\u00f3digo da comuna, para que possam ser agrupados por comuna.<\/p>\n\n\n\n Vamos utilizar a camada de comunas do registo predial e efetuar uma jun\u00e7\u00e3o por localidade com a nova camada fundida.<\/p>\n\n\n\n Temos agora a tabela em anexo:<\/p>\n\n\n\n Nesta imagem, temos os primeiros cinco registos correspondentes a \u00e1reas que n\u00e3o s\u00e3o localidades. O atributo comuna<\/em><\/strong>, que nos interessa para o passo seguinte, n\u00e3o est\u00e1 preenchido, mas o atributo id2<\/em><\/strong>, que prov\u00e9m da uni\u00e3o por localidade com a camada de comunas, est\u00e1 corretamente preenchido.<\/p>\n\n\n\n Abra a tabela no modo de edi\u00e7\u00e3o. Utilize uma express\u00e3o para selecionar os registos que n\u00e3o t\u00eam o atributo comuna <\/em><\/strong>preenchido e, em seguida, abra a calculadora de campos e actualize o campo comuna <\/em><\/strong>com o conte\u00fado de id2<\/em><\/strong>.<\/p>\n\n\n\n Agora podemos agrupar as localidades por comuna para obter a nossa camada de comunas que \u00e9 coerente com as localidades.<\/p>\n\n\n\n Utilizamos o tratamento Grupo<\/em><\/strong>, com a camada resultante da jun\u00e7\u00e3o que acab\u00e1mos de modificar.<\/p>\n\n\n\n Ent\u00e3o temos as nossas novas comunas. Resta apenas uma \u00faltima opera\u00e7\u00e3o a efetuar. Como em qualquer opera\u00e7\u00e3o de fus\u00e3o, o resultado conter\u00e1 uma s\u00e9rie de pequenas impurezas que precisam de ser eliminadas:<\/p>\n\n\n\n Para o fazer, vamos utilizar o tratamento Remover buracos.<\/p>\n\n\n\n N\u00e3o se esque\u00e7a de guardar o resultado num ficheiro permanente – seria dif\u00edcil ter de refazer tudo!<\/span><\/p>\n\n\n\n Podemos agora comparar o resultado da nossa nova camada comum com a camada original<\/p>\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n\n
Como corrigir a topologia da camada de localidade<\/h2>\n\n\n\n
1-Criar uma camada com um buffer \u00e0 volta das localidades<\/h3>\n\n\n\n
<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n2- Carregar o layer lieux-dits no PostgreSQL<\/h3>\n\n\n\n
<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n3- Executar o script SQL<\/h3>\n\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-- recuperamos todos os pol\u00edgonos da tabela de buffer em duas colec\u00e7\u00f5es: uma onde temos as geometrias completas e outra onde temos 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 dos pol\u00edgonos 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 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 ;\n\/*\nVamos criar uma tabela com 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 s\u00e3o cont\u00edguos, 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, removendo assim as partes em conflito*\/\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\/*As intersec\u00e7\u00f5es acabadas de subtrair devem ser adicionadas 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);\n\/*Uma chave prim\u00e1ria \u00e9 ent\u00e3o redefinida para que o QGIS possa carregar a camada*\/\n\nALTER TABLE table_corrigee ADD CONSTRAINT pk_table_corrigee PRIMARY KEY (id);<\/em>\n<\/code><\/pre>\n\n\n\n\n
<\/a><\/figure>\n\n\n<\/a>
\u00c0 esquerda a camada corrigida, \u00e0 direita a camada de localidade original<\/figcaption><\/figure>\n<\/div>\n\n\nComo recriar uma camada de comunas que seja coerente com a nova camada de lieux-dits<\/h2>\n\n\n\n
<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n1-Cria\u00e7\u00e3o dos pol\u00edgonos em falta<\/h3>\n\n\n\n
<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n