Keywords: rural-limit geopolitical-limit rural-zone dissolve geometry-calculator advanced-labeling
En esta actividad evaluaremos los límites veredales y sus diferencias a partir del análisis de diferentes fuentes de información. Luego, crearemos el límite territorial municipal global que será utilizado para el recorte de la información geo-espacial obtenida en actividades posteriores.
Tomado de: https://commons.wikimedia.org
- Evaluar los límites geopolíticos veredales e identificar y entender sus diferencias.
- Crear el límite municipal global.
- 🎓Actividad: Inventario de información geo-espacial recopilada del POT y diccionario de datos.
- 🧰Herramienta: Microsoft Excel 365.
- 🧰Herramienta: ESRI ArcGIS Pro 3.3.1 o superior.
- 🧰Herramienta: QGIS 3.38 o superior.
- 📂VeredaLimite.xlsx: libro para registro y análisis de áreas veredales
- 📂Veredas DANE 2020: Parte1, Parte2
- Abra el proyecto de ArcGIS Pro creado en la actividad anterior y desde el menú Insert cree un nuevo mapa New Map y renombre como CountyLimit. Podrá observar que ahora disponemos de dos mapas en el panel del catálogo.
- En la tabla de contenido o Contents, de clic derecho sobre el nombre del mapa CountyLimit, y seleccione la opción Propiedades o Properties. En la pestaña Coordinate Systems, busque, seleccione y asigne el sistema de proyección de coordenadas 9377, correspondiente al Orígen Único Nacional de Colombia.
- Desde el Panel de Catálogo (Catalog) localizado a la derecha, cargue al mapa desde la ruta
\R.SIGE\file\data\POT\Anexo_Acuerdo_012_2013\gdb\25899.gdb\Rural\, la clase de entidad denominada VEREDA. De clic en el ícono de simbología de la capa y establezca borde negro y sin relleno.
- Seleccione en Contents la capa VEREDA, en la parte superior podrá observar que se activan opciones específicas para esta capa, vaya al menú Labeling y establezca el rotulado de la capa utilizando a partir del nombre de la vereda.
- En el panel Contents, seleccione la capa, abra la tabla de atributos u oprima Ctrl+T y con la opción Add agregue los siguientes campos y guarde las modificaciones realizadas en la tabla cerrando la ventana de modificación o mediante la opción guardar o Save que aparece en el menú superior:
| Campo | Descripción | Tipo | Propiedad ArcGIS Pro |
|---|---|---|---|
| APha | Área planar en hectáreas | Double | Area |
| AGha | Área geodésica en hectáreas | Double | Area (geodesic) |
| PPm | Perímetro planar en metros | Double | Perimeter length |
| PGm | Perímetro geodésico en metros | Double | Perimeter length (geodesic) |
Tenga en cuenta que al crear los campos no se calculan automáticamente estas propiedades geométricas.
Los valores planares corresponden a los calculados a partir de la proyección de la capa sobre un plano horizontal.
Los valores geodésicos corresponden a los calculados teniendo en cuenta la curvatura terrestre.
En la tabla podrá observar que esta capa se compone de 14 entidades o veredas.
- En la tabla de atributos y dando clic sobre el campo
APha, seleccione la opción Calculate Geometry que le permitirá calcular las propiedades geométricas asociadas a las entidades.
Calcule las propiedades solicitadas indicando las unidades de cálculo y el sistema de proyección 9377, correspondiente a MAGNA Orígen Único Nacional de Colombia.
En actividades anteriores evidenciamos que el sistema de proyección de coordenadas de las capas recopiladas del POT es el 3116, correspondiente a MAGNA Sirgas Orígen Bogotá.
Tenga en cuenta que el cálculo de las propiedades geométricas puede variar entre sistemas de coordenadas, por lo cual, los valores calculádos automáticamente (por estar contenida la capa dentro de una GDB) en esta capa en los camposShape_AreayShape_length(a partir del CRS 3116), son diferentes a los obtenidos con el sistema 9377.
Incluya un campo adicional tipo entero largo con el nombre Nodos y desde el editor de geometría, calcule el número de nodos que componen cada entidad, utilizando para ello la propiedad geométrica número de vértices (Number of Vertices).
💡Recuerde que luego de incluir o modificar campos de atributos en tablas o capas, estos deben ser documentados en el diccionario de datos.
- Modifique el rótulo agregando el campo de área y perímetro geodésico, utilice una de las siguientes expresiones:
| Lenguaje | Sentencia |
|---|---|
| Arcade | $feature.nombre + textformatting.NewLine + "A (ha): " + Round($feature.AGha, 2) + textformatting.NewLine + "P (m):" + Round($feature.PGm, 2) |
| VBScript | [nombre] & VbNewLine & "A (ha): " & Round([AGha],2) & VbNewLine & "P (m):" & Round([PGm], 2) |
| Python | [nombre] + "\nA (ha): " + str(round(float([AGha]), 2)) + "\nP (m):" + str(round(float([PGm]), 2)) |
| JScript | [nombre] + "\nA (ha): " + parseFloat([AGha]).toFixed(2) + "\nP (m):" + parseFloat([PGm]).toFixed(2) |
En Arcade, el salto de línea también puede ser definido con la expresión "\n".
- En el panel lateral Contents, de clic derecho en VEREDA y seleccione la opción Create Chart / Bar Chart, cree una gráfica en orden descendente para los valores de área geográfica, podrá observar que la vereda Páramo de Guerrero es la de mayor extensión.
Para iniciar, en el panel lateral Contents, apague la capa VEREDA.
- Desde la misma ruta anterior
\R.SIGE\file\data\POT\Anexo_Acuerdo_012_2013\gdb\25899.gdb\Rural\, agregue la clase de entidad TERRENO_PREDIO_RURAL. Cambie la simbología a polígonos sin relleno.
- Abra la tabla de atributos e identifique el campo que contiene el código de identificación de vereda y cambie la simbología de representación de símbolo simple a valores únicos utilizando este campo. Podrá observar que los polígonos de predios de una misma vereda tienen el mismo color.
- En el menú superior, de clic en el grupo Analysis y de clic en Tools. Podrá observar que a la derecha aparece el panel Geoprocessing. En el campo de búsqueda, ingrese la palabra Dissolve, obtendrá diferentes resultados.
- En la ventana de resultados, seleccione la herramienta Dissolve asociada al grupo Data Management Tools e ingrese los parámetros definidos en la siguiente imagen. Nombre la capa como Vereda_DiagnosticoPOT_TerrenoPredioRural.shp
Para la ruta de salida de todas las capas en formato shapefile a generar en este curso, utilizaremos
\R.SIGE\file\shp.
Como puede observar, se han creado 14 veredas, sin embargo, contienen regiones internas vacías debidas a los corredores viales que no se incluyen en las capas de predios catastrales como unidades prediales, por lo que será necesario editar la geometría de la capa y eliminar manualmente los nodos internos para crear una región contínua. Para este ejercicio, no eliminaremos estos nodos internos.
- En la tabla de contenido, desplace el mapa VEREDA y ubíquelo sobre el mapa de veredas creadas a partir de predios, podrá observar que existe correspondencia exterior en casi todas las zonas.
Se recomienda recorrer todos los bordes perimetrales e identificar en que zonas existen diferencias importantes entre estas dos capas.
- En la tabla de atributos, agregue los campos requeridos (APha, AGha, PPm, PGm, Nodos) y realice los cálculos geométricos utilizando el CRS 9377.
Rotule la capa de veredas utilizando la expresión Arcade: $feature.vereda_id + textformatting.NewLine + "A (ha): " + Round($feature.AGha, 2) + textformatting.NewLine + "P (m):" + Round($feature.PGm, 2)
En la tabla de atributos podrá observar que a partir de la disolución hemos obtenido el conteo de predios en cada vereda, encontrando que la vereda Paso Ancho es la que tiene el mayor número de divisiones prediales con 2775 predios.
Como observa en la imagen, los valores de área de las veredas calculadas previamente son diferentes a los obtenidos a partir de predios. Para conocer el área correspondiente a vías + corredores de cauces + espacios públicos sin predio catastral, puede restar los valores de estas dos áreas.
En la tabla podrá observar que esta capa también se compone de 14 entidades o veredas.
Desde la ruta \R.SIGE\file\data\POT\Anexo_Acuerdo_012_2013\shp\, agregue la capa vectorial shapefile VEREDA.shp, cambie la simbología, agregue y calcule los campos de atributos geométricos indicados al inicio de esta actividad, y realice la rotulación. Siga el mismo procedimiento realizado en el numeral 1.1 de esta actividad.
Rótulo Arcade: $feature.nombre + textformatting.NewLine + "A (ha): " + Round($feature.AGha, 2) + textformatting.NewLine + "P (m):" + Round($feature.PGm, 2)
Para diferenciar las capas de veredas, utilice diferentes colores de borde y desde la tabla de contenido renombre las capas. Por ejemplo, la capa VEREDA de la formulación puede ser renombrada como Veredas (2013 formulación POT).
- Desde la plataforma de https://www.colombiaenmapas.gov.co/, descargue la capa de veredas en formato shapefile de Colombia creada por el DANE. Obtendrá un archivo comprimido, descomprima y guarde en la ruta
\R.SIGE\file\data\DANE\de su repositorio de datos.
Para la descarga se requiere de inicio de sesión de usuario utilizando, por ejemplo, una cuenta activa de Google.
- Desde la ruta de descompresión, cargue la capa de veredas de Colombia al mapa, cambie la simbología de representación y abra la tabla de atributos.
En la tabla de atributos podrá observar que existen diferentes campos asociados a los municipios de Colombia.
- Desde la tabla de contenido, de clic derecho en la capa y seleccione la opción de propiedades, de clic en la pestaña de creación de consultas o Definition Query y establezca un filtro para nombre de municipio igual a Zipaquirá.
Query filter: NOMB_MPIO = 'ZIPAQUIRÁ'
Podrá observar que se obtienen 14 resultados correspondientes a las veredas del municipio del caso de estudio.
- Desde la tabla de contenido, de clic derecho en la capa y seleccione la opción Data / Export Featues y exporte los resultados del filtro a una nueva capa, guarde en la ruta
\R.SIGE\file\shp\como VeredasDANE2020.shp.
- Cambie la simbología, agregue y calcule los campos de atributos geométricos indicados al inicio de esta actividad, y realice la rotulación. Siga el mismo procedimiento realizado en el numeral 1.1 de esta actividad.
Rótulo Arcade: $feature.nombre + textformatting.NewLine + "A (ha): " + Round($feature.AGha, 2) + textformatting.NewLine + "P (m):" + Round($feature.PGm, 2)
Si revisa con atención la tabla veredal, podrá observar que los códigos catastrales y los nombres de las veredas han sido modificados. En las versiones utilizadas para la elaboración del POT, las veredas iniciaban en el código catastral 02 y en esta versión en 01.
Con respecto a las áreas veredales del diagnóstico 2010 vs. formulación 2013 del POT, se puede observar en las tablas de resultados que el área aumentó en las veredas: Portachuelo, El Tunal, La Granja, El Cedro y Paso Ancho. No existen áreas con disminución debido a que en la modificación se ajustó el límite catastral veredal con incorporaciónes sobre el área catastral.
Con respecto a las áreas veredales del diagnóstico 2010 vs. DANE 2020, se puede observar en las tablas de resultados que el área aumentó y disminuyo, asjustándose los límites territoriales. También, los códigos catastrales han cambiado y se han renombrado algunas veredas.
Registre el resumen de cálculo de propiedades geométricas por vereda en el libro VeredaLimite.xlsx.
A partir de las capas de veredas evaluadas en los numerales anteriores y utilizando la capa del modelo de ordenamiento territorial - MOT disponible en la información recopilada del POT en la ruta \R.SIGE\file\data\POT\Anexo_Acuerdo_012_2013\shp\MOT.shp, cree los polígonos de delimitación del municipio. El proceso consiste en utilizar la herramienta Data Management Tools / Dissolve utilizada anteriormente en esta actividad.
- Disuelva las capas veredales, guarde en la ruta
\R.SIGE\file\shp\y nombre como:
| Capa veredal | Límite territorial disuelto |
|---|---|
| Veredas (2010 diagnóstico POT): \file\data\POT\Anexo_Acuerdo_012_2013\gdb\25899.gdb\RURAL\VEREDA | Mpio25899_DiagnosticoPOT2010.shp |
| Veredas (2013 formulación POT): \file\data\POT\Anexo_Acuerdo_012_2013\shp\VEREDA.shp | Mpio25899_FormulacionPOT2013.shp |
| Veredas (DANE 2020): \file\shp\VeredasDANE2020.shp | Mpio25899_DANE2020.shp |
| Modelo de ordenamiento territorial 2013 - MOT: \file\data\POT\Anexo_Acuerdo_012_2013\shp\MOT.shp | Mpio25899_MOT2013.shp |
Límite municipal a partir de veredas diagnóstico 2010 POT: Mpio25899_DiagnosticoPOT2010.shp
Límite municipal a partir de veredas formulación 2013 POT: Mpio25899_FormulacionPOT2013.shp
Límite municipal a partir de veredas DANE 2020: Mpio25899_DANE2020.shp
Límite municipal a partir de modelo de ordenamiento territorial 2013 MOT: Mpio25899_MOT2013.shp
Observará que a nivel interno aparece la región vacía correspondiente a los límites urbanos. A partir de diferencia entre el área del límite municipal y el límite veredal, podrá obtener el tamaño del área urbana.
En las ilustraciones, podrá observar que internamente existen otras subzonas o puntos correspondientes a errores de digitalización por falta de coalineación entre los bordes de las entidades a partir de las cuales se realizo la disolución.
- En la tabla de contenido, seleccione una de las capas de límites disueltos y desde el menú Edit, seleccione la herramienta Modify / Edit Vertices que le permitirá modificar los vértices de las entidades. Al seleccionar una entidad y activar el modo de edición de vértices, los nodos perimetrales son activados y también el despliegue de la barra inferior de edición de nodos.
La herramienta de edición de vértices solo se activará si primero selecciona una entidad con la herramienta Select.
Seleccione la herramienta de eliminación de nodos (puntero de edición con símbolo menos) y elimine los nodos internos. Tenga cuidado de eliminar nodos en la zona perimetral externa a conservar.
Para terminar la edición, en el panel superior del menú Edit, seleccione la opción de guardado o Save.
- Agregue los campos para el cálculo de las propiedades geométricas y obtenga estos valores. Rotule las entidades.
Rótulo Arcade: "A (ha): " + Round($feature.AGha, 2) + textformatting.NewLine + "P (m):" + Round($feature.PGm, 2)
Como observa en las ilustraciones, existen algunas diferencias en los valores de área calculados y los límites son diferentes en algunas zonas.
| Límite territorial disuelto | APha | AGha | PPm | PGm | Nodos | AGha veredal | AGha urbana | % cambio urbano |
|---|---|---|---|---|---|---|---|---|
| Mpio25899_DiagnosticoPOT2010.shp | 19497.37 | 19522.38 | 107553.74 | 107623.00 | 25951 | 18726.86 | 795.52 | 0.0% |
| Mpio25899_FormulacionPOT2013.shp | 19497.37 | 19522.38 | 107553.74 | 107623.00 | 25951 | 18582.34 | 940.04 | 18.2% |
| Mpio25899_MOT2013.shp | 19497.39 | 19522.39 | 107576.30 | 107645.57 | 22966 | 18586.90 | 935.49 | 17.6% |
| Mpio25899_DANE2020.shp | 19457.06 | 19482.02 | 107445.24 | 107514.37 | 13292 | 18457.49 | 1024.53 | 28.8% |
Como observa en la tabla, los valores del área municipal son similares o idénticas en la información utilizada para la elaboración del POT y el área DANE 2020 es inferior debido a actualización de límites catastrales perimetrales a partir de predios. Con respecto al tamaño del área urbana, el área veredal de la formulación aumento, sin embargo, el área final definida en el MOT es un poco inferior a la definida en la formulación. Con respecto al área DANE 2020, se puede evidenciar que catastralmente la zona urbana ha aumentado en 229 hectáreas (28.8%) con respecto al área veredal catastral 2010.
AGha veredal: obtenida de la sumatoria de valores obtenidos en áreas de cada vereda.
AGha urbana: obtenida a partir de la diferencia del área total municipal y área veredal.
Para los posteriores análisis realizados en este curso, utilizaremos como máscara de referencia, el polígono obtenido a partir de los límites del Modelo de Ordenamiento Territorial - MOT.
Para el desarrollo de las actividades desarrolladas en esta clase, se pueden utilizar en QGIS las siguientes herramientas o geo-procesos:
| Proceso | Procedimiento |
|---|---|
| Definición de CRS | En la parte inferior derecha de QGIS se encuentra la opción EPSG que permite cambiar el sistema de proyección del mapa. |
| Simbología | Modificable desde las propiedades de la capa en la pestaña Symbology. |
| Rotulado | Modificable desde las propiedades de la capa en la pestaña Labels. |
| Agregar campo | Modificable desde las propiedades de la capa en la pestaña Fields o desde la tabla de atributos. |
| Cálculos geométricos o de campo | Directamente desde la tabla de atributos mediante el botón Open Field Calculator o Ctr+I. La geometría de cálculo $area permite obtener el valor elipsoidal y area el valor proyectado. |
| Disolución | Se ejecuta desde el Processing Toolbox / Vector Geometry / Dissolve o desde el menú Vector / Geoprocessing Tools / Dissolve. |
| Edición geométrica para eliminación de zonas internas | Activar modo de edición en la capa o Toggle Editing, luego en la barra de edición seleccionar la herramienta Vertex Tool, oprimir y mantener la tecla Shift, seleccionar los nodos a eliminar, oprimir tecla Delete. |
Ejemplo rótulo en QGIS: 'A(ha): ' || round("AGha", 2) || '\n' || 'P (m): ' || round("PGm", 2)
Agregue a la tabla resúmen generada en la actividad Inventario de información geo-espacial recopilada del POT y diccionario de datos, las capas generadas en esta actividad que se encuentran listadas a continuación:
| Nombre | Descripción | Geometría | Registros |
|---|---|---|---|
| Vereda_DiagnosticoPOT_TerrenoPredioRural.shp | Veredas generadas a partir de la capa de predios catastro 2010 utilizados en el diagnóstico y formulación del POT (TERRENO_PREDIO_RURAL) | Polígono 2D | 14 |
| VeredasDANE2020.shp | Veredas DANE 2020 Municipio de Zipaquirá a partir de la capa de veredas de Colombia disponible en www.colombiaenmapas.gov.co | Polígono 2D | 14 |
| Mpio25899_DiagnosticoPOT2010.shp | Límite municipal obtenido a partir de la disolución del mapa veredal del diagnóstico 2010 del POT | Polígono 2D | 1 |
| Mpio25899_FormulacionPOT2013.shp | Límite municipal obtenido a partir de la disolución del mapa veredal de la formulación 2013 del POT | Polígono 2D | 1 |
| Mpio25899_DANE2020.shp | Límite municipal obtenido a partir del mapa veredal DANE 2020 | Polígono 2D | 1 |
| Mpio25899_MOT2013.shp | Límite municipal obtenido a partir del mapa del Modelo de Ordenamiento Territorial 2013 - MOT | Polígono 2D | 1 |
💡Para funcionarios que se encuentran ensamblando el SIG de su municipio, se recomienda incluir y documentar estas capas en el Diccionario de Datos.
En la siguiente tabla se listan las actividades que deben ser desarrolladas y documentadas por cada grupo de proyecto en un único archivo de Adobe Acrobat .pdf. El documento debe incluir portada (indicando el caso de estudio, número de avance, nombre del módulo, fecha de presentación, nombres completos de los integrantes), numeración de páginas, tabla de contenido, lista de tablas, lista de ilustraciones, introducción, objetivo general, capítulos por cada ítem solicitado, conclusiones y referencias bibliográficas.
| Actividad | Alcance |
|---|---|
| Avance P2 | Para su caso de estudio, evalúe las capas de veredas utilizadas en el diagnóstico y la formulación. |
| Avance P2 | A partir de la base predial utilizada en el diagnóstico y la formulación, cree la capa de veredas. |
| Avance P2 | Utilizando la base de datos predial más reciente, cree la capa de veredas. |
| Avance P2 | Descargue y extraiga las veredas a partir de la última versión DANE disponible. |
| Avance P2 | Cree las capas de límites territoriales, compare todos los límites veredales y territoriales obtenidos e identifique y evalúe en que zonas se han realizado mutaciones de suelo rural a urbano. Cree un libro de Excel como el mostrado en esta actividad. |
| Avance P2 | Cree diferentes perímetros urbanos a partir de la capas veredales y a partir también del mapa del modelo de ordenamiento territorial vigente, analice y explique sus diferencias. |
| Avance P2 | En una tabla y al final del informe de avance de esta entrega, indique el detalle de las sub-actividades realizadas por cada integrante de su grupo. Para actividades que no requieren del desarrollo de elementos de avance, indicar si realizo la lectura de la guía de clase y las lecturas indicadas al inicio en los requerimientos. Utilice las siguientes columnas: Nombre del integrante, Actividades realizadas, Tiempo dedicado en horas. |
No es necesario presentar un documento de avance independiente, todos los avances de proyecto de este módulo se integran en un único documento.
En el informe único, incluya un numeral para esta actividad y sub-numerales para el desarrollo de las diferentes sub-actividades, siguiendo en el mismo orden de desarrollo presentado en esta actividad.
- https://www.colombiaenmapas.gov.co/
- https://www.qgis.org/en/docs/index.html
- https://docs.qgis.org/3.34/en/docs/user_manual/processing_algs/qgis/vectorgeometry.html
- Geovisor del municipio de Zipaquirá
| Versión | Descripción | Autor | Horas |
|---|---|---|---|
| 2024.03.08 | Versión inicial con alcance de la actividad | rcfdtools | 4 |
| 2024.07.07 | Investigación y documentación para caso de estudio general | rcfdtools | 8 |
| 2024.07.08 | Investigación y documentación para caso de estudio general | rcfdtools | 8 |
R.SIGE es de uso libre para fines académicos, conoce nuestra licencia, cláusulas, condiciones de uso y como referenciar los contenidos publicados en este repositorio, dando clic aquí.
¡Encontraste útil este repositorio!, apoya su difusión marcando este repositorio con una ⭐ o síguenos dando clic en el botón Follow de rcfdtools en GitHub.
| 🏠 Inicio | 🔰 Ayuda / Colabora | Siguiente |
|---|