![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2023\/06\/a1.png?ssl=1\")
<\/a><\/figure>\n\n\n\nDe hecho, el l\u00edmite de cada localidad es el borde de la carretera o camino colindante, lo que deja espacios vac\u00edos donde est\u00e1n las carreteras o los caminos. 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 esta escala no parece muy limpio, pero si se ampl\u00eda es a\u00fan peor.<\/p>\n\n\n\n Y si nos fijamos en la propia capa de localidades, encontramos numerosas anomal\u00edas topol\u00f3gicas:<\/p>\n\n\n\n As\u00ed es, es el registro oficial de la propiedad\u2026<\/p>\n\n\n\n Aqu\u00ed veremos c\u00f3mo resolver ambos problemas:<\/p>\n\n\n\n El ejemplo siguiente se ha creado a partir de los datos del catastro oficial del departamento 71, descargados del enlace https:\/\/cadastre.data.gouv.fr\/data\/etalab-cadastre\/2023-01-01\/shp\/departements\/71\/<\/a><\/p>\n\n\n\n La soluci\u00f3n aqu\u00ed propuesta requiere una base de datos PostgreSQL con Postgis. Es posible adaptarla a SQLite para utilizar un archivo geopackage<\/a> y el QGis DB Manager. Pero hay que contar con tiempos de procesamiento muy largos y sortear algunas deficiencias y errores. Aqu\u00ed veremos la mejor soluci\u00f3n, pero presentaremos la soluci\u00f3n con Geopackage en otro art\u00edculo.<\/p>\n\n\n\n Para gestionar los espacios entre pol\u00edgonos, aunque algunas herramientas de encaje pueden tratar ciertos problemas, la soluci\u00f3n m\u00e1s sencilla y radical es construir un buffer 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 10m. Construimos un b\u00fafer de 6 m alrededor de las localidades, lo que dar\u00e1 una superposici\u00f3n de al menos 2-3 metros.<\/p>\n\n\n\n Elije un archivo de salida permanente en formato Geopackage, que te permitir\u00e1 seguir el procedimiento del paso siguiente sin modificaciones.<\/p>\n\n\n\n El resultado de la capa b\u00fafer es el siguiente:<\/p>\n\n\n\n Ya no tenemos que gestionar los espacios entre los pol\u00edgonos, s\u00f3lo los solapamientos.<\/p>\n\n\n\n Puedes seguir los detalles de c\u00f3mo cargar la capa en el art\u00edculo C\u00f3mo cargar un geopackage en Postgis con QGis<\/a> . Utiliza el mismo procedimiento para cargar shapefiles.<\/p>\n\n\n\n Y encuentras tu tabla en PostreSQL con pgAdmin<\/p>\n\n\n\n Vamos a utilizar un \u00fanico script para hacer todo el trabajo:<\/p>\n\n\n\n Este es el script:<\/p>\n\n\n\n Para adaptar el script a su caso, deber\u00e1 utilizar un procesador de textos para sustituir los siguientes elementos<\/p>\n\n\n\n Carga esta nueva capa en QGis:<\/p>\n\n\n\n Hemos resuelto el problema de la coherencia topol\u00f3gica de las localidades.<\/p>\n\n\n\n Si superponemos nuestra nueva capa con los l\u00edmites de los municipios del catastro, veremos que no son coherentes con nuestras nuevas localidades, al igual que no lo son con las localidades originales.<\/p>\n\n\n\n Para obtener unos l\u00edmites comunales coherentes, bastar\u00eda con agrupar las localidades por comunas. El problema es que no todos los municipios est\u00e1n divididos en localidades. Por tanto, tenemos \u00ab\u00a0lagunas\u00a0\u00bb en la capa de localidades, que son perfectamente normales.<\/p>\n\n\n\n Para construir municipios completos, primero hay que crear los pol\u00edgonos que faltan para las localidades, asignarles el c\u00f3digo de municipio correspondiente y, a continuaci\u00f3n, agrupar las localidades utilizando el atributo municipio.<\/p>\n\n\n\n En primer lugar, vamos a crear una nueva capa de pol\u00edgonos, emprise<\/em><\/strong>, con un \u00fanico pol\u00edgono que rodee completamente nuestra capa de localidades.<\/p>\n\n\n\n A continuaci\u00f3n, utilizamos el tratamiento \u00ab\u00a0Diferencia\u00a0\u00bb para conservar en nuestra tabla_corregida \u00fanicamente las zonas de emprise <\/em><\/strong>que no corresponden a las localidades existentes:<\/p>\n\n\n\n Puede dejar la salida de los tratamientos como capas temporales porque no tendr\u00e1 m\u00e1s uso para ellas despu\u00e9s del resultado final.<\/p>\n\n\n\n El resultado de la diferencia es:<\/p>\n\n\n\n El resultado es un \u00fanico multipol\u00edgono. Antes de continuar, necesitamos convertir esta capa multipoligonal en pol\u00edgonos simples, para poder eliminar el contorno del \u00e1rea.<\/p>\n\n\n\n Para ello utilizamos el procesamiento de pol\u00edgonos individuales<\/p>\n\n\n\n Una vez finalizado el procesamiento, cambie el resultado al modo de edici\u00f3n y seleccione el pol\u00edgono exterior.<\/p>\n\n\n\n Y se le suprime:<\/p>\n\n\n\n As\u00ed que ahora tenemos una capa de Geometr\u00edas Simples<\/em><\/strong> que contiene todos los lugares donde no hay localidades, y una capa de localidades corregidas<\/em><\/strong>. Fusionamos las dos capas en una \u00fanica capa sin huecos:<\/p>\n\n\n\n Queda por resolver el problema de los atributos de los pol\u00edgonos que no son localidades. Necesitan tener un atributo que contenga el c\u00f3digo del municipio, para que puedan agruparse por municipios.<\/p>\n\n\n\n Vamos a utilizar la capa de municipios del catastro y a realizar una uni\u00f3n por localidad con la nueva capa fusionada.<\/p>\n\n\n\n Ahora tenemos la tabla adjunta:<\/p>\n\n\n\n En esta imagen, tiene los cinco primeros registros correspondientes a zonas que no son localidades. El atributo municipio, que nos interesa para el siguiente paso, no est\u00e1 rellenado. Sin embargo, el atributo id2, que proviene de la uni\u00f3n por localidad con la capa municipios, s\u00ed est\u00e1 rellenado correctamente.<\/p>\n\n\n\n Abra la tabla en modo edici\u00f3n. Utilice una expresi\u00f3n para seleccionar los registros que no tienen rellenado el atributo commune <\/em><\/strong>y, a continuaci\u00f3n, abra la calculadora de campos y actualice el campo commune <\/em><\/strong>con el contenido de id2<\/em><\/strong>.<\/p>\n\n\n\n Ahora podemos agrupar las localidades por municipio para obtener nuestra capa municipios que es coherente con las localidades.<\/p>\n\n\n\n Utilizamos el tratamiento Grupo, con la capa resultante de la uni\u00f3n que acabamos de modificar.<\/p>\n\n\n\n As\u00ed que tenemos nuestras nuevas comunas. A\u00fan nos queda una \u00faltima operaci\u00f3n. Como en cualquier operaci\u00f3n de fusi\u00f3n, el resultado contendr\u00e1 una serie de peque\u00f1as escorias que habr\u00e1 que eliminar:<\/p>\n\n\n\n Para ello, utilizaremos el tratamiento Eliminar agujeros.<\/p>\n\n\n\n No olvides guardar el resultado en un archivo permanente: \u00a1ser\u00eda duro tener que rehacerlo todo!<\/span><\/p>\n\n\n\n Ahora podemos comparar el resultado de nuestra nueva capa com\u00fan con la capa original<\/p>\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n\n
C\u00f3mo corregir la topolog\u00eda de la capa de localidad<\/h2>\n\n\n\n
1-Crear una capa con un buffer alrededor de las 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- Cargar la capa localidades en PostgreSQL<\/h3>\n\n\n\n
<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n3- Ejecutar el 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 los pol\u00edgonos de la tabla buffer en dos colecciones: una donde tenemos las geometr\u00edas completas y otra donde tenemos las partes superpuestas. De momento, las dos 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 de los pol\u00edgonos de la siguiente manera\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\n AND tampon.id <> tampon_bis.id \n --para una intersecci\u00f3n de 2 pol\u00edgonos, mantener 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 ;\n\/*\nVamos a crear una tabla que contenga 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 son contiguos, 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; as\u00ed 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\/*Las intersecciones que se acaban de restar deben a\u00f1adirse a continuaci\u00f3n 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\/*A continuaci\u00f3n se redefine una clave primaria para que QGIS pueda cargar la capa.*\/\n\nALTER TABLE table_corrigee ADD CONSTRAINT pk_table_corrigee PRIMARY KEY (id);\n<\/code><\/pre>\n\n\n\n\n
<\/a><\/figure>\n\n\n<\/a>
A la izquierda la capa corregida, a la derecha la capa de localidad original<\/figcaption><\/figure>\n<\/div>\n\n\nC\u00f3mo recrear una capa de municipios coherente con la nueva capa de lieux-dits<\/h2>\n\n\n\n
<\/a><\/figure>\n\n\n\n<\/a><\/figure>\n\n\n\n1-Creaci\u00f3n de los pol\u00edgonos que faltan<\/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