Keywords: land-index land-occupation
Utilizando las tablas del registro 1 de catastro, calcule el índice de ocupación y el índice de construcción municipal en cada vereda, en la zona urbana y en cada manzana. Una los resultados obtenidos con la capa de veredas DANE y manzanas IGAC y cree un mapa que simbolice o represente los valores obtenidos.
- Estimar el índice de ocupación y construcción municipal.
- Representar geográficamente los índices obtenidos.
- Identificar diferentes estrategias para el cálculo de índices y conocer sus limitaciones.
- 🎓Actividad: Inventario de información geo-espacial recopilada del POT y diccionario de datos.
- 🎓Actividad: Análisis de destinaciones económicas IGAC (creación de dominios).
- 🎓Actividad: Análisis veredal y límite territorial.
- 🧰Herramienta: ESRI ArcGIS Pro 3.3.1 o superior.
- 🧰Herramienta: QGIS 3.38 o superior.
| Índice | Descripción1 | Interpretación |
|---|---|---|
| Indice de ocupación | Es la proporción del área de suelo que puede ser ocupada por edificación en primer piso bajo cubierta, y se expresa por el cociente que resulta de dividir el área que puede ser ocupada por edificación en primer piso bajo cubierta por el área total del predio. | Se interpreta como el desarrollo horizontal de un inmueble. |
| Índice de construcción | Es el número máximo de veces que la superficie de un terreno puede convertirse por definición normativa en área construida, y se expresa por el cociente que resulta de dividir el área permitida de construcción por el área total de un predio. | Se interpreta como el desarrollo vertical de un inmueble. |
Para el análisis de estos índices, se pueden utilizar diferentes estrategias:
| Estrategia | Descripción |
|---|---|
| A partir de registros catastrales | Las tablas del registro catastral 1 y 2 del IGAC, contiene información detallada de las características generales y físicas de los inmuebles. Se actualizan con mayor frecuencia que las bases prediales vectoriales, sin embargo, las áreas de predios y construcciones registradas pueden no corresponder con el área geodésica o planar de las capas de predios y construcciones. |
| A partir de información cartográfica | A partir de la cartografía digital de predios y construcciones y sus áreas planares o geodésicas, se puede realizar el análisis de ocupación, sin embargo, no es posible calcular indice de construcción. Generalmente, las bases vectoriales prediales, no contienen el mismo número de predios que contienen las tablas de registros catastrales. |
| Por foto-interpretación | Utilizando imágenes satelitales históricas, se puede evaluar el desarrollo progresivo de las construcciones y calcular las ocupaciones, sin embargo, la resolución espacial de captura es una limitante en construcciones pequeñas y el tiempo de procesamiento puede ser extenso si no se implementan algoritmos de aprendizaje. |
En esta actividad, desarrollaremos el análisis a partir de registros catastrales combinados con información geográfica.
- Abra el proyecto de ArcGIS Pro, creado previamente y desde el menú Insert cree un nuevo mapa New Map, renombre como LandIndex y establezca el CRS 9377. Agregue al mapa la capa de predios creada en la actividad Análisis de destinaciones económicas IGAC, nombrada como
\file\gdb\SIGE.gdb\SIGE\TerrenoPredio_2013, ajuste la simbología a valores únicos representando el campo de atributosvereda_id.
Previamente, a los predios de la zona urbana, se les asignó el código
vereda_id= 2589901, con el propósito de poder crear análisis estadísticos y poblar completamente este campo.
- Utilizando la herramienta de geo-procesamiento Data Managenet Tools / Dissolve, disuelva los predios a partir del campo
vereda_id, nombre como\file\gdb\SIGE.gdb\SIGE\Vereda_TerrenoPredio_2013. Abra la tabla de atributos, podrá observar que se obtienen 15 polígonos y que no se encuentra el nombre de las veredas y del área urbana.
La disolución de predios para obtener los límites geopolíticos internos del municipio, no incluyen las áreas correspondientes a los corredores viales. Para el análisis de índices de ocupación y construcción, utilizaremos únicamente el área correspondiente a las unidades inmobiliarias.
- Para incluir los nombres de las veredas en la capa disuelta, cree un campo de texto de 100 caracteres de longitud e identifíquelo como
ZonaGeo. Agregue la capa de veredas utilizadas en el diagnóstico del POT desde la ruta\file\gdb\SIGE.gdb\POT2013Formulacion\VEREDA1y realice un Join o unión de registros en la capa disuelta.
- Utilizando el calculador de campo sobre
ZonaGeo, asigne el nombre de vereda contenido en el camponombrede la capa VEREDA1.
- Remueva la unión de campo e ingrese manualmente el nombre del área urbana. Rotule a partir del campo
ZonaGeo.
Recuerde que en el análisis veredal realizado en la actividad Análisis veredal y límite territorial, evidenciamos que los códigos veredales y nombres de veredas DANE, fueron actualizados en el año 2020. Para este ejemplo, utilizaremos la codificación utilizada en la formulación del POT del año 2013.
- Agregue al mapa la tabla del registro 1 de catastro del IGAC, disponible en la ruta
\file\gdb\SIGE.gdb\IGAC2009Registro1y cree los siguientes atributos:
| Campo | Descripción | Tipo | Propiedad ArcGIS Pro |
|---|---|---|---|
| Zona | Zona urbana o rural | text (7) | n/a |
| vereda_id | Código veredal | text (30) | n/a |
- Utilizando el calculador de campo sobre
Zona, obtenga los 7 primeros dígitos del campoPRE_CODde cada uno de los 46305 registros de esta tabla (Zona=!PRE_COD![:7]).
- Seleccione por atributos, los registros rurales cuyo código de zona es 2589900 (15063 registros catastrales) y utilizando el calculador de campo, asigne en
vereda_id, los 13 primeros dígitos del campoPRE_COD(vereda_id=!PRE_COD![:13]).
- Seleccione por atributos, los registros urbanos cuyo código de zona es 2589901 (31242 registros catastrales) y utilizando el calculador de campo, asigne en
vereda_id, los 7 primeros dígitos del campoPRE_COD(vereda_id=!PRE_COD![:7]).
De esta forma, ha obtenido en cada registro catastral, los mismos códigos de vereda de los polígonos generados previamente en la disolución de predios.
- En la tabla de atributos IGAC2009Registro1, y utilizando la herramienta de resúmen estadístico o Summarize sobre el campo de atributos
vereda_id, genere una tabla que consolide el total de metros construídos en cada zona geográfica definida. Nombre la tabla como\file\gdb\SIGE.gdb\IGAC2009Registro1_IndConstGeneral.
- Abra la tabla de resumen generada, podrá observar que contienen el conteo de registros utilizados en cada zona, la sumatoria de área catastral construída y la sumatoria del área de las construcciones.
Tenga en cuenta que en predios que contienen mejoras y propiedades horizontales, el área total de terreno obtenida, incluye la sumatoria múltiple de las áreas del mismo predio. Por ejemplo, si un predio de 100 m² posee internamente 4 mejoras, la sumatoria de terreno obtenida en el resumen estadístico de este predio será de 400 m². Al calcular los índices de construcción a partir de esta sumatoria, obtendrá valores errados y menores al índice de construcción real del predio. Es por ello que el cálculo del índice general se debe realizar a partir del área planar o geodésica de los polígonos disueltos a partir de las unidades prediales.
- En la capa Vereda_TerrenoPredio_2013, agregue 3 campos numéricos dobles con los nombres
AGm2,AreaConsm2,IndConst. Luego cree un Join o unión de tablas con los registros obtenidos del resúmen estadístico IGAC2009Registro1_IndConstGeneral
- Utilizando el calculador de geometría, calcule el área geodésica
AGm2de cada zona geográfica y asigne enAreaConsm2el valor obtenido en el campoSUM_area_construida.
- Remueva la unión y calcule el índice de construcción dividiendo el área total construída entre el área total de la zona (
IndConst = !AreaConsm2!/!AGm2!). Modifique el rótulo incluyendo el índice obtenido. Podrá observar que el área urbana tiene el mayor índice correspondiente a 0.65 o 65% de construcción, y que las veredas por ser rurales, presentan índices comparativamente menores.
Rótulo Arcade: $feature.ZonaGeo + '\nIndConst: ' + round($feature.IndConst, 4) + " (" + round($feature.IndConst*100, 4) + "%)"
- Cree una gráfica de barras representando los índices obtenidos.
Para la estimación del índice de ocupación, es necesario excluir de la tabla de registros catastrales, los propietarios secudarios y los identificadores de condición de propiedad y número de construcción repetidos; para ello, incluiremos en el análisis solo los valores con los códigos indicados en la siguiente tabla:
| Condición de propiedad |
Descripción2 | Código a incluir |
|---|---|---|
| 0 | Predio no reglamentado en propiedad horizontal - PH | 000 |
| 9 | Predio en propiedad horizontal - PH | 901 |
| 8 | Predio en condominio | 801 |
| 7 | Parques cementerios | 701 |
| 6 | Mejoras por edificaciones en terreno ajeno en propiedad horizontal - PH | 601 |
| 5 | Mejoras por edificaciones en terreno ajeno de propiedades no reglamentadas en propiedad horizontal - PH | 501 |
| 4 | Vías | 401 |
| 3 | Bienes de uso público diferentes a las vías | 301 |
| 2 | No Ley 14 de 1983 | 201 |
El código de inclusión dependerá de las mutaciones catastrales contenidas en el municipio, en municipios con formación catastral antigua no actualizada, algunos de los códigos primarios (201, 301, 401...) no fueron conservados dinámicamente.
- En la tabla IGAC2009Registro1, cree un campo de texto de 3 caracteres con el nombre CPropCons, utilizando el calculador de campo, obtenga los 3 dígitos finales del campo
CODPROPIEDAD.
- Desde el Definition Query de la tabla de registros catastrales, filtre las filas correspondientes a números de orden 001 y condición de propiedad y construcción primarios. Obtendrá 27108 registros que no coinciden exactamente con el número de unidades prediales presentes en la capa TerrenoPredio_2013, correspondientes a 26304 (804 registros más).
La diferencia principal se debe a que el registro 1 contienen más registros de propiedad primaria que las unidades prediales incluídas en la capa de predios. Generalmente, los procesos de actualización predial no son incorporados en la capa de vectores inmediatamente de la expedición de resolución de subdivisión.
Definition Query para este ejemplo: num_orden = '001' And CPropCons IN ('000', '201', '801', '901', '329')
- En la tabla de atributos IGAC2009Registro1, y utilizando la herramienta de resúmen estadístico o Summarize sobre el campo de atributos
vereda_id, genere una tabla que consolide el total de metros construídos en cada zona geográfica definida. Nombre la tabla como\file\gdb\SIGE.gdb\IGAC2009Registro1_IndOcupGeneral.
💡Para obtener el área ocupada en primer planta, cree una columna de atributos numérica doble y divida los metros construídos de cada registro entre el número de pisos, luego consolide el total de metros ocupados. De esta forma obtendrá un valor más representativo de la ocupación en cada vereda y en el área urbana.
- Abra la tabla de obtenida y verifique su contenido. Verifique que la sumatoria de las frecuencias sea igual al total de registros filtrados, correspondientes a 27108, para ello, en la cabecera de columna del campo
FRECUENCYde clic derecho y seleccione la opción Visualize Statistics.
- En la capa Vereda_TerrenoPredio_2013, agregue 2 campos numéricos dobles con los nombres
AreaOcupm2eIndOcup. Luego cree un Join o unión de tablas con los registros obtenidos del resúmen estadístico IGAC2009Registro1_IndOcupGeneral
- Utilizando el calculador de campo, asigne en
AreaOcupm2el valor obtenido en el campoSUM_area_construida. Como observa, el área ocupada es menor que el área construída.
- Remueva la unión y calcule el índice de ocupación dividiendo el área total ocupada entre el área total de la zona (
IndOcup = !AreaOcupm2!/!AGm2!). Modifique el rótulo incluyendo el índice obtenido. Podrá observar que el área urbana tiene el mayor índice correspondiente a 0.65 o 65% de construcción, y que las veredas por ser rurales, presentan índices comparativamente menores.
Rótulo Arcade: $feature.ZonaGeo + '\nIndConst: ' + round($feature.IndConst, 4) + " (" + round($feature.IndConst*100, 4) + "%)\nIndOcup: " + round($feature.IndOcup, 4) + " (" + round($feature.IndOcup*100, 4) + "%)"
- Cree una gráfica de barras representando los índices obtenidos.
- Desde la ruta
\file\gdb\SIGE.gdb\IGAC2013Urbano\. agregue al mapa las capas MANZANA, EDIFICACION y CONSTRUCCION_ANEXA. Simbolice utilizando colores que permitan diferencias las construcciones y rotule las manzanas a partir de campocodigo.
- Utilizando la herramienta de geo-procesamiento Data Management Tools / Merge, combine los polígonos de las construcciones en una única capa, nombre como
\file\gdb\SIGE.gdb\SIGE\Construccion_2013.
- En la capa combinada, cree un campo de numérico doble con el nombre
AGm2y calcule el área geodésica de cada construcción, también cree un campo de texto de 30 caracteres de longitud con el nombreManzanaIdy asigne los primeros 13 caracteres del campoterreno_predio_id.
- Utilizando la herramienta Summarize sobre el campo
ManzanaIdde la capa de construcciones integradas, genere un resúmen estadístico totalizando el área ocupada en cada manzana. Nombre como\file\gdb\SIGE.gdb\Construccion_2013_Manzana. Como observa, en cada manzana ha obtenido el número de construcciones que la ocupan y el área tota de ocupación.
- En la capa de manzanas urbanas, cree dos campos numéricos dobles con los nombres
AreaOcupm2eIndOcup. Luego realice un Join o unión con la tabla de resumen estadístico de construcciones por manzana.
- Utilizando el calculador de campo sobre
AreaOcupm2de la capa de manzanas, asigne el valorSUM_AGm2de la tabla de resúmen estadístico.
- Remueva la unión y en el campo
IndOcup, calcule el índice de ocupación dividiendo el total del área ocupada entre el área total de la manzana.
Tenga en cuenta que para las manzanas sin construcciones, el indicador no será calculado y Python devolverá un mensaje a advertencia indicando que no se calcularon todos los valores solicitados.
En algunos casos, el índice de ocupación puede ser mayor que 1.000, esto debido a que pueden existir construcciones que cubran toda la manzana y que tengan voladizos o cubiertas fuera de ella. Modifique el rótulo de manzana incluyendo el valor del índice calculado.
Rótulo Arcade: $feature.codigo + "\nIndOcup: " + round($feature.IndOcup, 4) + " (" + round($feature.IndOcup*100, 4) + "%)"
- Desactive las construcciones y modifique la simbología de las manzanas utilizando una representación por Colores Graduados por quantiles a partir del campo del índice de ocupación. Incluya también una gráfica de estadística visual a partir del campo del índice, y una gráfica de dispersión comparando las áreas de las manzanas vs. el área total construída. Podrá observar que las ocupaciones altas se encuentran en las manzanas pequeñas, que los índices no se encuentran distribuídos al rededor de una tendencia centrada y que no hay una relación directa dominante entre el tamaño de las manzanas y sus áreas ocupadas.
- Ajuste a simbología por colores graduados, el análisis realizado a nivel general en las veredas y el área urbana. Podrá observar que las zonas contiguas al área urbana se encuentran más ocupadas, y que las zonas de protección que previamente analizamos en otras actividades del curso, tienen una ocupación baja.
Para el desarrollo de las actividades desarrolladas en esta clase, se pueden utilizar en QGIS las siguientes herramientas o geo-procesos:
| Proceso | Procedimiento |
|---|---|
| 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. |
| Disolución (Dissolve) | Se ejecuta desde el Processing Toolbox / Vector Geometry / Dissolve o desde el menú Vector / Geoprocessing Tools / Dissolve. |
| Combinación de capas (Merge) | Herramienta disponible en el Processing Toolbox / Vector General / Merge vector layers. |
| Unión de tablas | En las propiedades de capa geográfica, seleccionar la pestaña Join y realizar la unión. Solo se mostrarán los registros correspondientes a la primer coincidencia encontrada. Para visualizar o representar otros valores, es necesario filtrar previamente la tabla a unir. |
Ejemplo rótulo en QGIS: 'A(ha): ' || round("AGha", 2) || '\n' || 'P (m): ' || round("PGm", 2)
📓QGIS training manual
📓Herramientas comúnmente utilizadas en QGIS
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_TerrenoPredio_2013 | Disolución predial por vereda y zona urbana a partir de la capa TerrenoPredio_2013. | Polígono 2D | 15 |
| Construccion_2013 | Combinación de edificaciones y construcciones anexas a partir de las capas EDIFICACION y CONSTRUCCION_ANEXA. | Polígono 2D | 41877 |
| IGAC2009Registro1_IndConstGeneral | Tabla de resumen estadístico de áreas construídas por vereda y zona urbana general, a partir de la tabla IGAC2009Registro1. | n/a (tabla) | 15 |
| IGAC2009Registro1_IndOcupGeneral | Tabla de resumen estadístico de áreas ocupadas por vereda y zona urbana general, a partir de la tabla IGAC2009Registro1. | n/a (tabla) | 15 |
| Construccion_2013_Manzana | Tabla de resumen estadístico de áreas ocupadas por manzana urbana, a partir de la tabla IGAC2009Registro1. | n/a (tabla) | 598 |
💡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 P4 | Realice el análisis de índices presentado en esta actividad. Complementariamente, realice el cálculo del índice de construcción en cada manzana urbana a partir de los valores reportados en la tabla del Registro 1 del IGAC. |
| Avance P4 | Descargue la base predial y registros de la última actualización catastral de su caso de estudio, realice un análisis de índices actualizado y compare con los valores obtenidos del estudio y adopción del POT. |
| Avance P4 | 🧭Mapa digital impreso P4-4: Indices generales de ocupación y construcción por vereda y en el área urbana. Incluir vectores, tablas y gráficos de análisis. Embebido dentro del informe final como una imágen y referenciados como anexo. |
| Avance P4 | 🧭Mapa digital impreso P4-5: Índices de ocupación y construcción por manzana urbana. Incluir vectores, tablas y gráficos de análisis. Embebido dentro del informe final como una imágen y referenciados como anexo. |
| Avance P4 | 🧭Mapa digital impreso P4-6: Comparación índices adopción POT vs. ultima actualización catastral. Incluir vectores, tablas y gráficos de análisis. Embebido dentro del informe final como una imágen y referenciados como anexo. |
| Avance P4 | 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://geoportal.igac.gov.co/contenido/datos-abiertos-catastro
- https://datos.icde.gov.co/search?collection=dataset&q=CATASTRO
- Reglamentación técnica de la formación catastral, la actualización de la formación catastral y la conservación catastral, IGAC.
| Versión | Descripción | Autor | Horas |
|---|---|---|---|
| 2024.03.21 | Versión inicial con alcance de la actividad | rcfdtools | 4 |
| 2024.07.26 | 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 |
|---|