Keywords: arcpy env arcpy.env.workspace arcpy.env.overwriteOutput datetime time Kc Gravelius Dd Dc arcpy.GetCount_management() arcpy.Describe() shapeType arcpy.ListFields() matplotlib arcpy.Statistics_analysis arcpy.SearchCursor() getValue() arcpy.Dissolve_management() arcpy.JoinField_management arcpy.AddField_management() arcpy.CalculateGeometryAttributes_management() arcpy.Intersect_analysis() arcpy.Statistics_analysis() arcpy.JoinField_management() arcpy.TableToExcel_conversion() open() .write() close()
La zonificación hidrográfica de Colombia desde el punto de vista hidrológico, tiene sus inicios en el HIMAT mediante la Resolución 0337 del 1978, la cual establece que el país está conformado por cinco Áreas hidrográficas (1-Caribe, 2- Magdalena - Cauca, 3- Orinoco, 4- Amazonas y 5-Pacífico) que a su vez están divididas en Zonas Hidrográficas y subdivididas en Subzonas Hidrográficas. En ese entonces, el propósito de la zonificación fue de adoptar un sistema de codificación para estaciones Hidrometerológicas. Posteriormente, el IDEAM introduce esta zonificación para otros fines, tales como estudios y análisis hidrológicos relacionados con los informes ambientales, p.ej, el Índice de Aridez, el Escurrimiento y el Rendimiento Hídrico.1
La zonificación de cuencas hidrográficas corresponde a tres niveles de jerarquía: áreas, zonas y subzonas hidrográficas. Las áreas hidrográficas corresponden a las regiones hidrográficas o vertientes que, en sentido estricto, son las grandes cuencas que agrupan un conjunto de ríos con sus afluentes que desembocan en un mismo mar. Ahora bien, en Colombia se distinguen cuatro vertientes, dos de ellas asociadas a ríos de importancia continental (vertiente del Orinoco y vertiente del Amazonas) y las vertientes del Atlántico y del Pacífico. Se delimita adicionalmente como áea hidrográfica la cuenca Magdalena-Cauca, que aunque tributa y forma parte de la vertiente del Atlántico, tiene importancia socioeconómica por su alto poblamiento y aporte al producto interno bruto.2
| AH | Área Hidrográfica |
|---|---|
| 1 | Caribe |
| 2 | Magdalena-Cauca |
| 3 | Orinoco |
| 4 | Amazonas |
| 5 | Pacífico |
Las cuencas hidrográficas que entregan o desembocan sus aguas superficiales directamente de una área hidrográfica se denominaran zonas hidrográficas. Agrupan varias cuencas que se presentan como un subsistema hídrico con características de relieve y drenaje homogéneo y sus aguas tributan a través de un afluente principal hacia un área hidrográfica. Están integradas por cuencas de las partes altas, medias o bajas de una zona hidrográfica que captan agua y sedimentos de los tributarios de diferente orden tales como nacimientos de agua, arroyos, quebradas y ríos. Las cuencas que tributan sus aguas a su vez a las zonas hidrográficas se denomina subzonas hidrográficas. Ahora bien, respecto a la toponimia con que se identifican zonas y subzonas hidrográficas, a estas unidades se les asignó la toponimia de acuerdo con el nombre de la corriente más representativa o río principal o con el nombre heredado de la zonificación del HIMAT, que puede corresponder al espacio geográfico o región a la cual drenan las aguas superficiales.2
| AH | Área Hidrográfica | ZH | Zona Hidrográfica |
|---|---|---|---|
| 1 | Caribe | 11 | Atrato - Darién |
| 1 | Caribe | 12 | Caribe - Litoral |
| 1 | Caribe | 13 | Sinú |
| 1 | Caribe | 15 | Caribe - La Guajira |
| 1 | Caribe | 16 | Catatumbo |
| 1 | Caribe | 17 | Islas del Caribe |
| 2 | Magdalena - Cauca | 21 | Alto Magdalena |
| 2 | Magdalena - Cauca | 22 | Saldaña |
| 2 | Magdalena - Cauca | 23 | Medio Magdalena |
| 2 | Magdalena - Cauca | 24 | Sogamoso |
| 2 | Magdalena - Cauca | 25 | Bajo Magdalena - Cauca - San Jorge |
| 2 | Magdalena - Cauca | 26 | Cauca |
| 2 | Magdalena - Cauca | 27 | Nechí |
| 2 | Magdalena - Cauca | 28 | Cesar |
| 2 | Magdalena - Cauca | 29 | Bajo Magdalena |
| 3 | Orinoco | 31 | Inírida |
| 3 | Orinoco | 32 | Guaviare |
| 3 | Orinoco | 33 | Vichada |
| 3 | Orinoco | 34 | Tomo |
| 3 | Orinoco | 35 | Meta |
| 3 | Orinoco | 36 | Casanare |
| 3 | Orinoco | 37 | Arauca |
| 3 | Orinoco | 38 | Orinoco Directos |
| 3 | Orinoco | 39 | Apure |
| 4 | Amazonas | 41 | Guainía |
| 4 | Amazonas | 42 | Vaupés |
| 4 | Amazonas | 43 | Apaporis |
| 4 | Amazonas | 44 | Caquetá |
| 4 | Amazonas | 45 | Yarí |
| 4 | Amazonas | 46 | Caguán |
| 4 | Amazonas | 47 | Putumayo |
| 4 | Amazonas | 48 | Amazonas - Directos |
| 4 | Amazonas | 49 | Napo |
| 5 | Pacífico | 51 | Mira |
| 5 | Pacífico | 52 | Patía |
| 5 | Pacífico | 53 | Tapaje - Dagua - Directos |
| 5 | Pacífico | 54 | San Juan |
| 5 | Pacífico | 55 | Baudó - Directos Pacífico |
| 5 | Pacífico | 56 | Pacífico - Directos |
| 5 | Pacífico | 57 | Islas del Pacífico |
En el presente análisis no se han incluido resultados para la ZH - zona hidrográfica 57, correspondiente a las Islas del Pacífico, debido a que la capa geográfica SZH - subzonas hidrográficas no contiene el polígono de delimitación.
Estudiar la forma y densidad de las áreas, zonas y subzonas hidrográficas de Colombia a partir de la delimitación geográfica realizada por el Instituto de Hidrología, Meteorología y Estudios Ambientales - IDEAM de Colombia, adscrito al Ministerio de Medio Ambiente - Minambiente y la red de drenajes sencillos digitalizada a escala 1:100k por el Instituto Geográfico Agustín Codazzi - IGAC.
Análisis de forma y densidad
| Parámetro | Descripción | Fórmula |
|---|---|---|
| Kc | Índice de compacidad Gravelius | Kc = P / (2π ( A / π ) ^ 0.5 ) Kc = 0.28209479179826 * P / A ^ 0.5 |
| Dd | Densidad de drenajes, km/km² | Dd = Σ LCi / A |
| Dc | Densidad de corrientes, 1/km² | Dd = n / A |
Donde,
- P, perímetro de la zona hidrográfica en kilómetros.
- A, área de la zona hidrográfica hen km².
- LCi, longitud de cada cauce o corriente de agua.
- n, número de cauces o corrientes de agua.
Para la marcación complementaria del coeficiente de compacidad Kc, se utilizará la clasificación de la Autoridad Nacional de Licencias Ambientales de Colombia - ANLA, cuya clasificación en el modelo nacional de datos geográficos establecido en el dominio Dom_Forma_IndComp relacionado con la capa geográfica CuencaHidrografica del componente HIDROLOGIA, indica que se subclasifican en los siguientes 3 tipos: 3
| Código | Descripción | Valor de referencia 4 |
|---|---|---|
| 1001 | Casi redonda a oval redonda | 1 a 1.25 |
| 1002 | Oval-redonda a oval oblonga | 1.25 a 1.5 |
| 1003 | Oval-oblonga a rectangular-oblonga | > 1.5 |
- Identificar las SZH - subzonas hidrográficas de Colombia - Suramérica.
- A partir de las SZH - subzonas hidrográficas, crear las capas disueltas de AH - áreas hidrográficas y ZH - zonas hidrográficas.
- Mediante intersección espacial, determinar el número y longitud de drenajes para cada zonificación.
- Estimar los coeficientes de forma y densidad para cada zonificación.
- Ejecutar varias de las funciones de geoprocesamiento disponibles en ArcPy de ArcGIS.
- Python 2.7+ de ArcGIS for Desktop 10.6+
- ArcGIS Pro 2.9+
- PyCharm 2021.3+ Commmunity
- Sistema operativo Microsoft Windows
- matplotlib
Para la ejecución del script principal, no es necesaria la instalación de librerías adicionales en ArcGIS, sin embargo, para la versión Desktop es requerida como mínimo la versión 10.6 para el cálculo directo de los parámetros geométricos de las entidades a través de la función
CalculateGeometryAttributes_management().
- Creación automática de log de resultados en formato Markdown con impresión de versión en archivo de salida.
- Parametrización general con definición de rutas, nombres de campos de entrada y activación o desactivación de subprocesos intermedios principales. El usuario puede modificar el arreglo de datos para la sub-clasificación de las zonificaciones por rango de área y los clasificadores del índice de compacidad.
- Impresión de propiedades y entidades encontradas para las capas de entrada.
- Estadísticos generales de drenajes por subtipo.
- Evaluación a partir de todos los subtipos de drenajes disponibles (sin subtipo, permanentes, intermitentes) o solo para drenajes permanentes. Los archivos de registro de resultados contienen la etiqueta
AllyPerque identifica el tipo de drenajes utilizados. - Reproyección automática de la capa original SZH - subzonas hidrográficas al sistema de coordenadas
9377 MAGNA-SIRGAS / Origen-Nacional. La capa original de drenajes no se reproyecta debido a que se encuentra en este sistema. - Disolución espacial de SZH - subzonas hidrográficas para obtención de entidades multiparte.
- A partir de la capa de subzonas se generan las capas de orden superior ZH - Zona hidrográfica y AH - Área hidrográfica.
- Inclusión de campos y cálculo automático de propiedades geométricas de área, perímetro y longitud de cada drenaje en sistema internacional de unidades.
- Intersección espacial de SZH - subzonas hidrográficas con drenajes.
- Estadísticos de conteo y sumatoria de longitudes de drenaje para cada zonificación (AH, ZH, SZH).
- Unión de resultados estadísticos a capas geográficas.
- Análisis de distribución de áreas por rangos para niveles superiores de zonificación (AH, ZH).
- Análisis de forma y densidad para Kc - Coeficiente de compacidad, Dd - Densidad de drenajes y Dc - Densidad de corrientes.
- Etiqueta de marcado para Kc - Coeficiente de compacidad indicando la forma.
- Visualización detallada en pantalla de los resultados obtenidos.
- Exportación de resultados a formato de Microsoft Excel.
La capa de drenajes utiliza el sistema de proyección
9377 MAGNA-SIRGAS / Origen-Nacionalque contiene sistema proyectado con unidades en metros y para el cálculo de las longitudes requeridas por este script, no es necesaria su reproyección. La capa de subzonas hidrográficas requiere ser proyectada al origen nacional debido a que se encuentra referenciada en4686 MAGNAsolo en sistema geográfico.
Para el desarrollo de este microcontenido se recomienda que los scripts y demás archivos requeridos se alojen en D:\R.GISPython\HydroGeoZone\ utilizando la siguiente estructura de directorios.
| Directorio | Descripción |
|---|---|
| /Datos | Directorio de datos de entrada. |
| /Graph | Directorio con gráficas e ilustraciones. |
| /Map | Directorio con mapas para visualización de resultados en ArcGIS for Desktop y ArcGIS Pro. |
| /Output | Directorio de salida de capas generadas y archivos de resultados estadísticos. |
| /Reference | Referencias bibliográficas. |
| /Screenshot | Capturas de pantalla con resultados de ejecución. |
Este mapa representa las unidades de análisis para el ordenamiento ambiental de territorio definidas por el Ideam en convenio con el Instituto Geográfico Agustín Codazzi (IGAC), a escala 1:500.000 llamadas zonificación hidrográfica de Colombia. 5
Procedimiento:
-
Ingrese al portal http://www.ideam.gov.co/en/capas-geo y en el cuadro de búsqueda escriba Zonificación Hidrográfica, observará que a 2021.12.23 existen dos versiones de la capa de zonificación correspondientes al año 2010 y 2013.
-
Para las dos capas, realice la descarga del archivo de formas Shapefile y consulte sus metadatos y el catálogo de objetos disponible.
-
Descomprima solo los datos contenidos en las carpetas
/shapeen la carpeta/Datadentro de D:\R.GISPython\HydroGeoZone.
Catálogo de objetos en Subzonas 6
| Nombre | Alias | Definición | Tipo de dato |
|---|---|---|---|
| OBJECTID | OBJECTID | Identificador de objeto geográfico. | Texto |
| Shape | Shape | Tipo de geometría. | Geometría |
| COD_AH | Código Area | Código del Area hidrográfica a la que corresponde. | Entero |
| COD_ZH | Código Zona | Código de la Zona hidrográfica a la que corresponde. | Entero |
| COD_SZH | Código Subzona | Código de Subzona hidrográfica a la que corresponde. | Entero |
| NOM_AH | Nombre Área | Nombre del área hidrográfica a la que corresponde. Dominio Área Hidográfica. | Texto |
| NOM_ZH | Nombre Zona | Nombre de la zona hidrográfica a la que corresponde. | Texto |
| NOM_SZH | Nombre Subzona | Nombre de la Subzona hidrográfica a la que corresponde. | Texto |
| Shape_Length | Shape_Length | Perímetro en las unidades del sistema de referencia espacial. | Entero |
| Shape_Area | Shape_Area | Área en las unidades del sistema de referencia espacial. | Entero |
| RULEID | RULEID | Id único asignado por el sistema a la representación gráfica. | Entero |
| Override | Override | Representación gráfica. | Blob |
Flujo de agua superficial que depende de la precipitación pluvial y/o afloramiento de aguas subterráneas y va a desembocar en otra corriente, en una laguna o en el mar. Los drenajes dispersos son aquellos que no desembocan en otro cuerpo de agua, o desaparecen al ser no fotointerpretables, por ejemplo en corrientes subterráneas.7
Procedimiento:
- Ingrese al portal https://www.colombiaenmapas.gov.co/, en temática seleccione Cartografía Básica y busque
Base de datos vectorial básica. Colombia. Escala 1:100.000. Año 2019. En la parte inferior del Detalle del Servicio seleccione en Formato de descargaGeodatabasey de clic en Descargar, automáticamente iniciará la descarga a través de una orden de servicio. La GDB comprimida tiene un tamaño apoximado de 661 MB.
- Descomprima la base de datos geográfica en la carpeta de descargas, cree un mapa nuevo en blanco en ArcGIS for Desktop o en ArcGIS Pro, agregue el mapa base World Light Gray Canvas Base y desde el dataset
Superficies_Agua, agregue la capaDrenaje_Sencillo.
- Exporte la clase de entidad Drenaje_Sencillo a un archivo de formas en formato Shapefile dentro de la carpeta
/Dataen D:\R.GISPython\HydroGeoZone. La versión descargada contiene 426719 entidades.
Para el desarrollo del análisis no se ha utilizado la digitalización de la Base de datos vectorial básica - Colombia a Escala 1:25.000 del Año 2018 debido a que aún no se encuentran todas las planchas del país digitalizadas y almacenadas en la GDB disponible para descarga como se muestra en la siguiente imagen.
Información disponible a escala 1:25k
Catálogo de objetos en Drenaje_Sencillo para capa en formato Shapefile
| Nombre | Alias | Definición | Tipo de dato |
|---|---|---|---|
| OBJECTID | OBJECTID | Identificador único de objeto. | Numérico |
| Shape | Shape | Tipo de geometría. | Geometría |
| ESTADO_DRE | TEDD | Indica si el drenaje es permanente o intermitente. Subtipo. | Texto |
| SYMBOL | Symbol | En este campo se especifica la geometría como se representará el objeto (punto, línea o polígono). | Numérico y texto, 255 |
| PROYECTO | PROYECTO | Proyecto en el cual se desarrollaron los datos. | Numérico y texto, 255 |
| FECHA | FECHA | Fecha de realización de los datos. | Dato |
| DISPERSION | Dispersión | Indica si el drenaje es disperso no no. | Dato |
| NOMBRE_GEO | NMG | Nombre de la entidad geográfica. | Numérico y texto, 255 |
| PK_CUE | PK_CUE | Identificador único global de cada elemento. | Numérico |
| RULEID | RuleID | Identificador único de la representación. | Texto |
| GLOBALID | GLOBALID | Identificador global en la base de datos espacial. | Texto |
| SHAPE_Leng | SHAPE_Leng | Longitud del elemento. | Numérico |
Estado de drenajes - Subtipos
| Code | Description |
|---|---|
| 5101 | Permanente |
| 5102 | Intermitente |
- Para la ejecución correcta es necesario cerrar las aplicaciones de ArcGIS for Desktop antes de iniciar. En ArcGIS Pro a través de ejecución usando un Notebook desde la GUI, no es necesario cerrar la aplicación antes de iniciar la ejecución.
- Para la ejecución completa del análisis para drenajes permanentes, intermitentes y no clasificados, establecer las variables
intersectActive = Truepara volver a realizar la intersección espacial y calcular las longitudes de los drenajes intersecados ystatisticActive = Truepara volver a generar estadísticos en DBF y convertirlos a Excel. - Para analizar solo a partir de drenajes permanentes, establecer en
Truelas variables anteriores y establecer adicionalmente la variableonlyPermanentDrainActive = True. - Para realizar modificaciones en el Script e incluir nuevas funcionalidades y ejecutar pruebas de funcionamiento, se recomienda ejecutar todos los procesos incluyendo todos los drenajes (permanentes e intermitentes) y luego de obtener las capas principales del análisis espacial
DrenajeSencilloFiltro.shpyDrenajeSencilloIntersect.shp. Luego desactivar la ejecución de drenajes intersecados y recálculo de estadísticos detallados, de esta forma no tendrá que esperar (aproximadamente 10 minutos para 500k drenajes) a la creación completa de las capas principales resultantes del análisis.
Tenga en cuenta que la definición de las subzonas hidrográficas se realizó a escala 1:500k y los drenajes a escala 1:100k, por lo que espacialmente pueden existir fragmentos de tramos de drenaje que cruzan entre zonas. Los cálculos de densidad y forma se realizan a partir de la intersección espacial fraccionada de drenajes dentro de cada área teniendo en cuenta la consideración anterior, por tal motivo, el número de tramos drenajes de la capa original
Drenaje_Sencillo.shppuede ser diferente al número de tramos de la capa de intersección.
El detalle de los resultados obtenidos del análisis, ha sido generado automáticamente por el script en Python y puede ser consultado o descargado en formato Markdown a través de los enlaces presentados en la siguiente tabla.
| Reporte | Descripción |
|---|---|
| HydroGeoZoneDrainAll20211229.md | Análisis para SZH - Subzonas hidrográficas IDEAM 2013 y Drenajes Sencillos IGAC 2019 - Drenajes de todos los subtipos. |
| HydroGeoZoneDrainPer20211229.md | Análisis para SZH - Subzonas hidrográficas IDEAM 2013 y Drenajes Sencillos IGAC 2019 - Drenajes solo permanentes. |
A continuación se presentan los resultados obtenidos a partir de la intersección espacial entre las SZH - Subzonas hidrográficas IDEAM del año 2013 con los drenajes permanentes e intermitentes IGAC a 2019.
Debido a que existen algunas subzonas hidrográficas con el mismo código y como polígonos independientes, y debido a que el script principal realiza el proceso de unificación en entidades multiparte para garantizar que el conteo y estadísticos de análisis a partir de drenajes sea consistente en el proceso de intersección espacial, el número de polígonos de la capa original (316) corresponde a un total de 314 subzonas únicas como se muestra en los siguientes resultados.
Los resultados detallados de las sub-clasificaciones de área espacial para cada AH - Área hidrográfica se encuentran en los reportes detallados.
AH - Áreas hidrográficas año 2013 con drenajes permanentes e intermitentes a 2019. .zip.
AH - Áreas hidrográficas año 2013 solo drenajes permanentes a 2019. .zip
ZH - Zonas hidrográficas año 2013 con drenajes permanentes e intermitentes a 2019. .zip.
Mapa de áreas por ZH - Zona hidrográfica
Mapa de número de drenajes por ZH - Zona hidrográfica
Mapa de sumatoria de longitud de drenajes por ZH - Zona hidrográfica
Mapa de Kc - Coeficiente de Compacidad por ZH - Zona hidrográfica
Debido a que algunas ZH - Zonas Hidrográficas están compuestas por entidades multiparte, se pueden presentar valores altos en el cálculo del índice de compacidad, p.ej, la zona
38 - Orinoco Directosse compone de 5 entidades polígonos.
Mapa de Dd - Densidad de Drenajes por ZH - Zona hidrográfica
Mapa de Dc - Densidad de Corrientes por ZH - Zona hidrográfica
ZH - Zonas hidrográficas año 2013 con solo drenajes permanentes 2019. .zip
SZH - Subzonas hidrográficas año 2013 con drenajes permanentes e intermitentes a 2019. .zip.
Mapa de áreas por SZH - Zona hidrográfica
Mapa de número de drenajes por SZH - Subzona hidrográfica
Mapa de sumatoria de longitud de drenajes por SZH - Subzona hidrográfica
Mapa de Kc - Coeficiente de Compacidad por SZH - Subzona hidrográfica
Mapa de Dd - Densidad de Drenajes por SZH - Subzona hidrográfica
Mapa de Dc - Densidad de Corrientes por SZH - Subzona hidrográfica
Todos los subtipos de drenaje
Solo drenajes subtipo permanente
Script HydroGeoZone.py
# -*- coding: UTF-8 -*-
# Nombre: HydroGeoZone.py
# Descripción: Zonificación hidrográfica de Colombia - Análisis de forma y densidad usando Python
# Requerimiento: PyCharm 2021.3+, ArcGIS 10.6+, ArcGIS Pro 2.9.0
# Librerías
import arcpy
from arcpy import env
import sys
import matplotlib
import matplotlib.pyplot as plt
from datetime import datetime
from datetime import date
import time
# Variables y datos de entrada
absolutePath = r'D:/R.GISPython/HydroGeoZone' # Usar r'.' para retornar a ruta relativa
arcpy.env.workspace = absolutePath+'/Data/'
arcpy.env.overwriteOutput = True
dataPath = absolutePath+'/Data/'
outputPath = absolutePath+'/Output/'
outputPathGraph = absolutePath+'/Graph/'
urlGitHubFile = 'https://github.com/rcfdtools/R.GISPython/blob/main/HydroGeoZone/'
hydroSubZoneLayerIn = dataPath+'Zonificacion_hidrografica_2013.shp'
drainageLayerIn = dataPath+'Drenaje_Sencillo.shp'
drainageLayer = outputPath+'DrenajeSencilloFiltro.shp' # Capa drenajes filtro solo permanentes
hydroAreaLayer = outputPath+'AreaHidrografica.shp'
hydroZoneLayer = outputPath+'ZonaHidrografica.shp'
hydroSubZoneLayerProject = outputPath+'SubZonaHidrograficaProject.shp'
hydroSubZoneLayerCopy = outputPath+'SubZonaHidrografica.shp'
drainageLayerIntersect = outputPath+'DrenajeSencilloIntersect.shp'
statisticsTableDrainageDBF = outputPath+'DrenajeSencilloSubtipoEstadistica.dbf'
statisticsTableAHDBF = outputPath+'AreaHidrograficaEstadistica.dbf'
statisticsTableZHDBF = outputPath+'ZonaHidrograficaEstadistica.dbf'
statisticsTableSZHDBF = outputPath+'SubZonaHidrograficaEstadistica.dbf'
statisticsTableAHXLS = outputPath+'AreaHidrograficaEstadistica.xls'
statisticsTableZHXLS = outputPath+'ZonaHidrograficaEstadistica.xls'
statisticsTableSZHXLS = outputPath+'SubZonaHidrograficaEstadistica.xls'
fieldAHCode, fieldZHCode, fieldSZHCode = 'COD_AH', 'COD_ZH', 'COD_SZH'
fieldAHName, fieldZHName, fieldSZHName = 'NOM_AH', 'NOM_ZH', 'NOM_SZH'
drainageSubtype, drainageLen = 'ESTADO_DRE', 'SHAPE_Leng'
drainageSubtypePerm = 5101 # Código de drenajes subtipo permanente
outCoordinateSystem = "PROJCS['MAGNA-SIRGAS / Origen-Nacional',GEOGCS['GCS_MAGNA',DATUM['D_MAGNA',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',5000000.0],PARAMETER['False_Northing',2000000.0],PARAMETER['Central_Meridian',-73.0],PARAMETER['Scale_Factor',0.9992],PARAMETER['Latitude_Of_Origin',4.0],UNIT['Meter',1.0]] # GEOGCS['GCS_MAGNA',DATUM['D_MAGNA',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
evalValueKc = [[1.25, 'Casi redonda a oval redonda'], [1.5, 'Oval-redonda a oval oblonga'], [999999, 'Oval-oblonga a rectangular-oblonga']] # Rangos coeficiente de compacidad Kc según ANLA - Colombia en modelo de datos GDB Nacional
evalAreaSZH = [[300, 0, 0], [700, 0, 0], [900, 0, 0], [1100, 0, 0], [1300, 0, 0], [1500, 0, 0], [2000, 0, 0], [2500, 0, 0], [3500, 0, 0], [5000, 0, 0], [10000, 0, 0], [20000, 0, 0], [999999, 0, 0]] # Valores de corte para evaluar número de subzonas, posición 0 corresponde al valor de corte, posición 1 para conteo y posición 2 para acumulado.
decimalPos = 2 # Posiciones decimales para impresión de tablas en formato Markdown
consKc = 0.28209479179826
intersectActive = False # Volver a realizar la intersección espacial y calcular las longitudes de los drenajes intersecados.
statisticActive = False # Volver a generar estadísticos en DBF y convertir a Excel.
onlyPerDrainActive = False # Analizar solo para drenajes permanentes.
# Log file creation
currentDate = date.today()
currentDateTxt = str(currentDate.year)+str(currentDate.month)+str(currentDate.day)
if onlyPerDrainActive:
fileNameAux = 'DrainPer'
titleAuxTxt = 'Solo drenajes subtipo permanente'
else:
fileNameAux = 'DrainAll'
titleAuxTxt = 'Todos los subtipos de drenaje'
fileLogName = absolutePath+'/HydroGeoZone'+fileNameAux+currentDateTxt+'.md'
fileLog = open(fileLogName, 'w+') # w+ para crear el archivo si no existe
timeStart = time.time()
# Función para crear separador de filas cabecera en formato Markdown
def TableHeadMarkdown(n=2):
lineSep = '|---'
PrintLog(lineSep * n + '|', True)
# Función de impresión en pantalla y log de resultados
def PrintLog(txtPrint, onScreen=True):
if onScreen: print(txtPrint)
fileLog.write(txtPrint + '\n')
# Función para consultar los campos de atributos disponibles
def CapaPropiedades(capa):
cont = 1
totalEntidades = arcpy.GetCount_management(capa)
descGeometria = arcpy.Describe(capa)
tipoGeometria = descGeometria.shapeType
PrintLog('### Campos en ' + capa + ' (' + tipoGeometria + 's ' + str(totalEntidades) + ')\n', True)
PrintLog('| # | Campo | Tipo |', True)
TableHeadMarkdown(3)
campos = arcpy.ListFields(capa)
for campo in campos:
PrintLog('| ' + str(cont) + ' | ' + campo.name + ' | ' + campo.type + ' |', True) # Print field properties
cont += 1
# Cabecera
PrintLog('## Zonificación hidrográfica de Colombia - Análisis de forma y densidad usando Python - ' + titleAuxTxt, True)
PrintLog('\n')
PrintLog('\n* Archivo de resultados: ' + fileLogName +
'\n* Fecha y hora de inicio de ejecución: ' + str(datetime.now()) +
'\n* Script compatible con: ArcGIS for Desktop 10.6+ y ArcGIS Pro'
'\n* Python versión: ' + str(sys.version) +
'\n* Python rutas: ' + str(sys.path[0:5]) +
'\n* matplotlib versión: ' + str(matplotlib.__version__) +
'\n* Encuentra este script en https://github.com/rcfdtools/R.GISPython/tree/main/HydroGeoZone'
'\n* Cláusulas y condiciones de uso en https://github.com/rcfdtools/R.GISPython/wiki/License'
'\n* Créditos: r.cfdtools@gmail.com\n', True)
PrintLog('```\nSistema de coordenadas: ' + outCoordinateSystem + '\n```', False) # Mostrar como código en Markdown
print('\nAntes de iniciar cierre las aplicaciones de ArcGIS for Desktop...')
PrintLog('\n### Parámetros y datos de entrada ([/Data](https://github.com/rcfdtools/R.GISPython/tree/main/HydroGeoZone/Data))\n'
'\n* Ruta absoluta: ' + absolutePath +
'\n* Espacio de trabajo ArcPy: ' + arcpy.env.workspace +
'\n* Capa de entrada SZH - Subzonas hidrográficas: ' + hydroSubZoneLayerIn +
'\n* Capa de entrada Drenajes: ' + drainageLayerIn +
'\n* Intersección SZH & Drenajes: ' + str(intersectActive) +
'\n* Recálculo de estadísticos: ' + str(statisticActive) +
'\n* Analizar solo drenajes permanentes: ' + str(onlyPerDrainActive))
PrintLog('\n## Propiedades y entidades encontradas para las capas de entrada\n')
CapaPropiedades(hydroSubZoneLayerIn)
PrintLog('\n', True)
CapaPropiedades(drainageLayerIn)
PrintLog('\n### Evaluación de drenajes por subtipo\n', True)
PrintLog('> (0 - Sin asignación, 5101 - Permanente, 5102 - Intermitente)\n', True)
arcpy.Statistics_analysis(drainageLayerIn, statisticsTableDrainageDBF, [[drainageLen, 'SUM']], [drainageSubtype])
cursor = arcpy.SearchCursor(statisticsTableDrainageDBF)
PrintLog('| Código | # Drenajes |', True)
TableHeadMarkdown(2)
for fila in cursor:
PrintLog('| ' + str(fila.getValue(drainageSubtype)) + ' | ' + str(fila.getValue('FREQUENCY')) + ' |', True)
if onlyPerDrainActive:
print('Filtrado de drenajes permanentes activado.'
'\nFiltrando drenajes solo permanentes en ' + drainageLayer +
'\nEste proceso tardará varios minutos...')
whereFilter = '"'+drainageSubtype+'"='+str(drainageSubtypePerm)
arcpy.Select_analysis(drainageLayerIn, drainageLayer, whereFilter)
else:
print('Filtrado de drenajes solo permanentes desactivado...')
drainageLayer = drainageLayerIn
print('\n')
# Procesos y cálculos
print('### Reproyección de subzonas')
print('Reproyectando SZH a '+hydroSubZoneLayerProject+'...')
print('Sistema de coordenadas: '+outCoordinateSystem)
arcpy.Project_management(hydroSubZoneLayerIn, hydroSubZoneLayerProject, outCoordinateSystem)
print('\n')
print('### Disolución de subzona, zona y área hidrográfica')
print('La capa de subzonas hidrográficas contiene polígonos con el mismo código y en entidades separadas por lo que se requiere su conversión a multiparte.')
print('Disolviendo a multiparte SZH - subzonas hidrográfica '+hydroSubZoneLayerCopy+'...')
arcpy.Dissolve_management(hydroSubZoneLayerProject, hydroSubZoneLayerCopy, fieldSZHCode, '', 'MULTI_PART', '')
print('Incorporando atributos a capa disuelta SZH - subzona hidrográfica ' + hydroSubZoneLayerCopy+'...')
arcpy.JoinField_management(hydroSubZoneLayerCopy, 'COD_SZH', hydroSubZoneLayerProject, 'COD_SZH') # Siempre se ejecuta antes de las demás disoluciones
print('Disolviendo a multiparte AH - áreas hidrográficas '+hydroAreaLayer+'...')
arcpy.Dissolve_management(hydroSubZoneLayerCopy, hydroAreaLayer, fieldAHCode, '', 'MULTI_PART', '')
print('Disolviendo a multiparte ZH - zonas hidrográficas '+hydroZoneLayer+'...')
arcpy.Dissolve_management(hydroSubZoneLayerCopy, hydroZoneLayer, fieldZHCode, '', 'MULTI_PART', '')
print('\n')
print('### Inclusión de campos para cálculo de áreas, perímetros, longitudes, forma y densidad')
print('Agregando campo area (Area) en km²')
print('Agregando campo perimetro (Perm) en km')
print('\tAH - áreas hidrográficas...')
arcpy.AddField_management(hydroAreaLayer, 'Area', 'DOUBLE')
arcpy.AddField_management(hydroAreaLayer, 'Perm', 'DOUBLE')
print('\tZH - zonas hidrográficas...')
arcpy.AddField_management(hydroZoneLayer, 'Area', 'DOUBLE')
arcpy.AddField_management(hydroZoneLayer, 'Perm', 'DOUBLE')
print('\tSZH - subzonas hidrográficas...')
arcpy.AddField_management(hydroSubZoneLayerCopy, 'Area', 'DOUBLE')
arcpy.AddField_management(hydroSubZoneLayerCopy, 'Perm', 'DOUBLE')
print('Agregando Kc - Coeficiente de Compacidad')
print('Agregando KcTag - Coeficiente de Compacidad - Rótulo')
print('Agregando Dd - Densidad de Drenaje km/km²')
print('Agregando Dc - Densidad de corrientes 1/Km²')
print('\tAH - áreas hidrográficas...')
arcpy.AddField_management(hydroAreaLayer, 'Kc', 'DOUBLE')
arcpy.AddField_management(hydroAreaLayer, 'KcTag', 'TEXT')
arcpy.AddField_management(hydroAreaLayer, 'Dd', 'DOUBLE')
arcpy.AddField_management(hydroAreaLayer, 'Dc', 'DOUBLE')
print('\tZH - zonas hidrográficas...')
arcpy.AddField_management(hydroZoneLayer, 'Kc', 'DOUBLE')
arcpy.AddField_management(hydroZoneLayer, 'KcTag', 'TEXT')
arcpy.AddField_management(hydroZoneLayer, 'Dd', 'DOUBLE')
arcpy.AddField_management(hydroZoneLayer, 'Dc', 'DOUBLE')
print('\tSZH - subzonas hidrográficas...')
arcpy.AddField_management(hydroSubZoneLayerCopy, 'Kc', 'DOUBLE')
arcpy.AddField_management(hydroSubZoneLayerCopy, 'KcTag', 'TEXT')
arcpy.AddField_management(hydroSubZoneLayerCopy, 'Dd', 'DOUBLE')
arcpy.AddField_management(hydroSubZoneLayerCopy, 'Dc', 'DOUBLE')
print('\n')
print('### Cálculo de áreas y perímetros')
print('Compatible con ArcGIS for Desktop 10.6+ y ArcGIS Pro.')
print('Calculando area (Area) en km² y perimetro (Perm) en km')
print('\tAH - áreas hidrográficas...')
arcpy.CalculateGeometryAttributes_management(hydroAreaLayer, [['Area', 'AREA'], ['Perm', 'PERIMETER_LENGTH']], 'KILOMETERS', 'SQUARE_KILOMETERS')
print('\tZH - zonas hidrográficas...')
arcpy.CalculateGeometryAttributes_management(hydroZoneLayer, [['Area', 'AREA'], ['Perm', 'PERIMETER_LENGTH']], 'KILOMETERS', 'SQUARE_KILOMETERS')
print('\tSZH - subzonas hidrográficas...')
arcpy.CalculateGeometryAttributes_management(hydroSubZoneLayerCopy, [['Area', 'AREA'], ['Perm', 'PERIMETER_LENGTH']], 'KILOMETERS', 'SQUARE_KILOMETERS')
print('\n')
print('### Intersección de drenajes con subzonas hidrográficas y cálculo de longitud por segmento')
print('Tenga en cuenta que la definición de las subzonas hidrográficas se realizó a escala 1:500k y los drenajes a escala 1:100k, por lo que espacialmente pueden existir fragmentos de tramos de drenaje que cruzan entre zonas. Los cálculos de densidad y forma se realizan a partir de la intersección espacial fraccionada de drenajes dentro de cada área teniendo en cuenta la consideración anterior, por tal motivo, el número de tramos drenajes de la capa original `Drenaje_Sencillo.shp` puede ser diferente al número de tramos de la capa de intersección.')
if intersectActive:
print('Intersecando drenajes con subzonas ' + drainageLayerIntersect+'...')
print('Este proceso tardara varios minutos...')
arcpy.Intersect_analysis([hydroSubZoneLayerCopy, drainageLayer], drainageLayerIntersect, 'ALL', '', 'LINE')
print('Intersección completada...')
print('Agregando campo longitud (LDre) en km a drenajes intersecados...')
arcpy.AddField_management(drainageLayer, 'LDre', 'DOUBLE')
print('Calculando longitud (LDre) en km para cada segmento de drenaje intersecado...')
print('Este proceso tardara varios minutos...')
arcpy.CalculateGeometryAttributes_management(drainageLayerIntersect, [['LDre', 'LENGTH']], 'KILOMETERS')
print('Calculo de longitudes completado...')
else:
print('Actualización de intersección de drenajes con subzonas desactivada...')
print('\n')
print('### Estadísticos de análisis')
if statisticActive:
print('Generando estadísticos en DBF')
print('\tAH - área hidrográfica ' + statisticsTableAHDBF)
arcpy.Statistics_analysis(drainageLayerIntersect, statisticsTableAHDBF, [['LDre', 'SUM']], [fieldAHCode, fieldAHName])
print('\tZH - zona hidrográfica ' + statisticsTableZHDBF)
arcpy.Statistics_analysis(drainageLayerIntersect, statisticsTableZHDBF, [['LDre', 'SUM']], [fieldAHCode, fieldAHName, fieldZHCode, fieldZHName])
print('\tSZH - subzona hidrográfica ' + statisticsTableSZHDBF)
arcpy.Statistics_analysis(drainageLayerIntersect, statisticsTableSZHDBF, [['LDre', 'SUM']], [fieldAHCode, fieldAHName, fieldZHCode, fieldZHName, fieldSZHCode, fieldSZHName])
print('\n')
else:
print('Actualización de estadísticos detallados desactivada...\n')
print('Unión de capas geográficas y estadísticos')
print('Tenga en cuenta que en algunas subzonas hidrográficas pueden no existir drenajes restituidos.')
print('Uniendo capa con estadísticos ' + hydroAreaLayer)
print('\tAH - área hidrográfica ' + hydroAreaLayer)
arcpy.JoinField_management(hydroAreaLayer, 'COD_AH', statisticsTableAHDBF, 'COD_AH')
print('\tZH - zona hidrográfica ' + hydroZoneLayer)
arcpy.JoinField_management(hydroZoneLayer, 'COD_ZH', statisticsTableZHDBF, 'COD_ZH')
print('\tSZH - subzona hidrográfica ' + hydroSubZoneLayerCopy)
arcpy.JoinField_management(hydroSubZoneLayerCopy, 'COD_SZH', statisticsTableSZHDBF, 'COD_SZH')
print('\n')
PrintLog('\n## Distribución de áreas en km² de subzonas por área hidrográfica', True)
PrintLog('\n### Total nacional de SZH - subzonas hidrográficas por rango de área\n', True)
print('Imprimiendo en formato Markdown...')
PrintLog('| Rango km² | # Subzonas | Acumulado |', True)
TableHeadMarkdown(3)
cont = 0
for i in evalAreaSZH:
cursor = arcpy.SearchCursor(hydroSubZoneLayerCopy)
for fila in cursor:
if fila.getValue('Area') <= i[0]: i[2] += 1
if cont != 0:
evalAreaSZH[cont][1] = evalAreaSZH[cont][2] - evalAreaSZH[cont-1][2]
PrintLog('| ' + str(evalAreaSZH[cont - 1][0]) + '-' + str(evalAreaSZH[cont][0]) + ' | ' + str(i[1]) + ' | ' + str(i[2]) + ' |', True)
else:
evalAreaSZH[cont][1] = evalAreaSZH[cont][2]
PrintLog('| 0-' + str(evalAreaSZH[cont][0]) + ' | ' + str(i[1]) + ' | ' + str(i[2]) + ' |', True)
cont += 1
print('\n')
PrintLog('\n### Total SZH - subzonas hidrográficas por rango de área para cada AH - área hidrográfica\n', True)
print('Imprimiendo en formato Markdown...')
cursorZH = arcpy.SearchCursor(hydroAreaLayer)
for filaZH in cursorZH:
for j in evalAreaSZH: # Reinicializar acumuladores
j[1] = 0
j[2] = 0
PrintLog('\n#### AH - Área hidrográfica ' + str(filaZH.getValue(fieldAHCode)) + ' - ' + str(filaZH.getValue(fieldAHName))+'\n', True)
PrintLog('| Rango km² | # Subzonas | Acumulado |', True)
TableHeadMarkdown(3)
cont = 0
for i in evalAreaSZH:
cursorSZH = arcpy.SearchCursor(hydroSubZoneLayerCopy)
for fila in cursorSZH:
if fila.getValue('Area') <= i[0] and fila.getValue(fieldAHCode) == filaZH.getValue(fieldAHCode): i[2] += 1
if cont != 0:
evalAreaSZH[cont][1] = evalAreaSZH[cont][2] - evalAreaSZH[cont-1][2]
PrintLog('| ' + str(evalAreaSZH[cont - 1][0]) + '-' + str(evalAreaSZH[cont][0]) + ' | ' + str(i[1]) + ' | ' + str(i[2]) + ' |', True)
else:
evalAreaSZH[cont][1] = evalAreaSZH[cont][2]
PrintLog('| 0-' + str(evalAreaSZH[cont][0]) + ' | ' + str(i[1]) + ' | ' + str(i[2]) + ' |', True)
cont += 1
print('\n')
print('### Análisis de forma y densidad')
codeEvalKc = '''def getKcTag(Kc):
for i in evalValueKc[:]:
if Kc <= i[0]:
return i[1]'''
print('Calculando Kc - Coeficiente de Compacidad, Dd - Densidad de Drenaje km/km² y Dc - Densidad de corrientes 1/Km²')
print('\tAH - áreas hidrográficas...')
arcpy.CalculateField_management(hydroAreaLayer, 'Kc', 'consKc*!Perm!/(!Area!**0.5)', 'Python')
arcpy.CalculateField_management(hydroAreaLayer, 'KcTag', 'getKcTag(float(!Kc!))', 'Python', codeEvalKc)
arcpy.CalculateField_management(hydroAreaLayer, 'Dd', '!SUM_LDre!/!Area!', 'Python')
arcpy.CalculateField_management(hydroAreaLayer, 'Dc', '!FREQUENCY!/!Area!', 'Python')
print('\tZH - zonas hidrográficas...')
arcpy.CalculateField_management(hydroZoneLayer, 'Kc', 'consKc*!Perm!/(!Area!**0.5)', 'Python')
arcpy.CalculateField_management(hydroZoneLayer, 'KcTag', 'getKcTag(float(!Kc!))', 'Python', codeEvalKc)
arcpy.CalculateField_management(hydroZoneLayer, 'Dd', '!SUM_LDre!/!Area!', 'Python')
arcpy.CalculateField_management(hydroZoneLayer, 'Dc', '!FREQUENCY!/!Area!', 'Python')
print('\tSZH - subzonas hidrográficas...')
arcpy.CalculateField_management(hydroSubZoneLayerCopy, 'Kc', 'consKc*!Perm!/(!Area!**0.5)', 'Python')
arcpy.CalculateField_management(hydroSubZoneLayerCopy, 'KcTag', 'getKcTag(float(!Kc!))', 'Python', codeEvalKc)
arcpy.CalculateField_management(hydroSubZoneLayerCopy, 'Dd', '!SUM_LDre!/!Area!', 'Python')
arcpy.CalculateField_management(hydroSubZoneLayerCopy, 'Dc', '!FREQUENCY!/!Area!', 'Python')
print('\n')
print('### Convirtiendo estadísticos y resultados a XLS')
if statisticActive:
print('\tAH - área hidrográfica ' + statisticsTableAHXLS)
arcpy.TableToExcel_conversion(hydroAreaLayer, statisticsTableAHXLS)
print('\tZH - zona hidrográfica ' + statisticsTableZHXLS)
arcpy.TableToExcel_conversion(hydroZoneLayer, statisticsTableZHXLS)
print('\tSZH - subzona hidrográfica ' + statisticsTableSZHXLS)
arcpy.TableToExcel_conversion(hydroSubZoneLayerCopy, statisticsTableSZHXLS)
else:
print('Actualización de conversión a XLS desactivada...')
print('\n')
PrintLog('\n## Visualización de tablas resultados para análisis de forma y densidad', True)
PrintLog('\n### AH - área hidrográfica\n', True)
PrintLog(statisticsTableAHDBF, True)
PrintLog('\n| AH | Nombre AH | Área, km² | Perm, km | n Drenajes | Sum. LCi, km | Kc | Dd | Dc | Kc Tag |', True)
TableHeadMarkdown(10)
cursor = arcpy.SearchCursor(hydroAreaLayer)
for fila in cursor:
PrintLog('| ' + str(fila.getValue(fieldAHCode)) + ' | ' + fila.getValue(fieldAHName) + ' | ' + str(round(fila.getValue('Area'), decimalPos)) + ' | ' + str(round(fila.getValue('Perm'), decimalPos)) + ' | ' + str(fila.getValue('FREQUENCY')) + ' | ' + str(round(fila.getValue('SUM_LDre'), decimalPos)) + ' | ' + str(round(fila.getValue('Kc'), decimalPos)) + ' | ' + str(round(fila.getValue('Dd'), decimalPos)) + ' | ' + str(round(fila.getValue('Dc'), decimalPos)) + ' | ' + fila.getValue('KcTag') + ' |', True)
PrintLog('\n### ZH - zona hidrográfica\n', True)
PrintLog(statisticsTableZHDBF, True)
PrintLog('\n| AH | Nombre AH | ZH | Nombre ZH | Área, km² | Perm, km | n Drenajes | Sum. LCi, km | Kc | Dd | Dc | Kc Tag |', True)
TableHeadMarkdown(12)
cursor = arcpy.SearchCursor(hydroZoneLayer)
for fila in cursor:
PrintLog('| ' + str(fila.getValue(fieldAHCode)) + ' | ' + fila.getValue(fieldAHName) + ' | ' + str(fila.getValue(fieldZHCode)) + ' | ' + fila.getValue(fieldZHName) + ' | ' + str(round(fila.getValue('Area'), decimalPos)) + ' | ' + str(round(fila.getValue('Perm'), decimalPos)) + ' | ' + str(fila.getValue('FREQUENCY')) + ' | ' + str(round(fila.getValue('SUM_LDre'), decimalPos)) + ' | ' + str(round(fila.getValue('Kc'), decimalPos)) + ' | ' + str(round(fila.getValue('Dd'), decimalPos)) + ' | ' + str(round(fila.getValue('Dc'), decimalPos)) + ' | ' + fila.getValue('KcTag') + ' |', True)
PrintLog('\n### SZH - Subzona hidrográfica\n', True)
PrintLog(statisticsTableSZHDBF, True)
PrintLog('\n| AH | Nombre AH | ZH | Nombre ZH | SZH | Nombre SZH | Área, km² | Perm, km | n Drenajes | Sum. LCi, km | Kc | Dd | Dc | Kc Tag |', True)
TableHeadMarkdown(14)
cursor = arcpy.SearchCursor(hydroSubZoneLayerCopy)
for fila in cursor:
PrintLog('| ' + str(fila.getValue(fieldAHCode)) + ' | ' + fila.getValue(fieldAHName) + ' | ' + str(fila.getValue(fieldZHCode)) + ' | ' + fila.getValue(fieldZHName) + ' | ' + str(fila.getValue(fieldSZHCode)) + ' | ' + fila.getValue(fieldSZHName) + ' | ' + str(round(fila.getValue('Area'), decimalPos)) + ' | ' + str(round(fila.getValue('Perm'), decimalPos)) + ' | ' + str(fila.getValue('FREQUENCY')) + ' | ' + str(round(fila.getValue('SUM_LDre'), decimalPos)) + ' | ' + str(round(fila.getValue('Kc'), decimalPos)) + ' | ' + str(round(fila.getValue('Dd'), decimalPos)) + ' | ' + str(round(fila.getValue('Dc'), decimalPos)) + ' | ' + fila.getValue('KcTag') + ' |', True)
PrintLog('\n### Graficas generales')
scatterVarX, scatterVarXLabel = 'Area', 'Área, km²'
scatterVarY = ['FREQUENCY', 'SUM_LDre', 'Kc', 'Dd', 'Dc']
scatterVarYLabel = ['# Drenajes', 'Sum. Long. Drenajes, km', 'Kc - Índice de Compacidad', 'Dd - Densidad de drenajes', 'Dc - Densidad de corrientes']
iLabel = 0
for iY in scatterVarY[:]:
xPlotValues, yPlotValues = [], []
PrintLog('\nGrafica ' + scatterVarXLabel + ' vs. ' + scatterVarYLabel[iLabel] + ' - ' + titleAuxTxt)
PrintLog('\n')
cursor = arcpy.SearchCursor(hydroSubZoneLayerCopy)
for fila in cursor:
xPlotValues.append(fila.getValue(scatterVarX))
yPlotValues.append(fila.getValue(iY))
plt.grid(color='0.95')
plt.scatter(xPlotValues, yPlotValues, c='k', label=scatterVarYLabel[iLabel], s=12)
plt.xlabel(scatterVarXLabel)
plt.ylabel(scatterVarYLabel[iLabel])
plt.legend(loc='lower right')
plt.title(scatterVarX + ' vs. ' + scatterVarYLabel[iLabel] + ' - ' + titleAuxTxt)
plt.savefig(outputPathGraph+'Plot'+scatterVarX+'Vs'+iY+fileNameAux+'.png')
plt.show()
iLabel += 1
PrintLog('\n### Archivos de resultados ([/Output](https://github.com/rcfdtools/R.GISPython/tree/main/HydroGeoZone/Output) [/Graph](https://github.com/rcfdtools/R.GISPython/tree/main/HydroGeoZone/Graph))\n'
'\n* Capa AH - Área hidrográfica: ' + hydroAreaLayer +
'\n* Capa ZH - Zona hidrográfica: ' + hydroZoneLayer +
'\n* Capa SZH - Subzona hidrográfica: ' + hydroSubZoneLayerCopy +
'\n* Capa Drenajes filtro permanentes: ' + drainageLayer +
'\n* Capa intersección SZH & Drenajes: ' + drainageLayerIntersect +
'\n* Tabla resultados AH - Área hidrográfica (dBase): ' + statisticsTableAHDBF +
'\n* Tabla resultados ZH - Zona hidrográfica (dBase): ' + statisticsTableZHDBF +
'\n* Tabla resultados SZH - Subzona hidrográfica (dBase): ' + statisticsTableSZHDBF +
'\n* Tabla resultados AH - Área hidrográfica (Excel): ' + statisticsTableAHXLS +
'\n* Tabla resultados ZH - Zona hidrográfica (Excel): ' + statisticsTableZHXLS +
'\n* Tabla resultados SZH - Subzona hidrográfica (Excel): ' + statisticsTableSZHXLS +
'\n* Tabla resultados Drenajes por subtipo: ' + statisticsTableDrainageDBF)
PrintLog('\n> Los mapas desplegados en [/Graph](https://github.com/rcfdtools/R.GISPython/tree/main/HydroGeoZone/Map) han sido generados manualmente en ArcGIS a partir de los datos obtenidos utilizando los mapas de proyecto disponibles en [/Map](https://github.com/rcfdtools/R.GISPython/tree/main/HydroGeoZone/Map).')
PrintLog('\nFecha y hora de terminación de ejecución: '+str(datetime.now()), True)
timeEnd = time.time()
PrintLog('\nProceso completado (dt = ' + str(round(timeEnd - timeStart, 1)) + ' sec o ' + str(round((timeEnd - timeStart)/60, 1)) + ' min)', True)
fileLog.close()
- http://www.ideam.gov.co/capas-geo
- http://www.siac.gov.co/catalogo-de-mapas
- http://visor.ideam.gov.co/geovisor/#!/profiles/3
- https://desktop.arcgis.com/en/arcmap/10.6/tools/data-management-toolbox/calculate-geometry-attributes.htm
- Hidrografía Colombiana - IDEAM y SiGaia (versión no oficial de zonificación a 2018)
- ESRI - Versiones de Python según la versión de ArcGIS
- ESRI - Python, NumPy, and MatPlotlib en ArcGIS 10.6
- Compatible con ArcGIS for Desktop 10.6 o superior.
- Compatible con ArcGIS Pro.
| Versión | Descripción | Autor | Horas |
|---|---|---|---|
| 2021.12.31 | Documentación. | rcfdtools | 12 |
| 2021.12.30 | Documentación. | rcfdtools | 12 |
| 2021.12.29 | Optimización de script. Inclusión de rutas de datos de entrada y salida en archivo Markdown de resultados. Exportación automática de gráficas de matriz de dispersión con nombramiento independiente para análisis a partir de todos los subtipos de drenaje o solo para permanentes. | rcfdtools | 12 |
| 2021.12.28 | Versión inicial. | rcfdtools | 12 |
R.GISPython 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.
| Anterior | 🏠 Inicio | 🔰 Ayuda / Colabora | Siguiente |
|---|
Footnotes
-
http://www.ideam.gov.co/web/agua/zonificacion-hidrografica ↩
-
http://documentacion.ideam.gov.co/openbiblio/bvirtual/022655/MEMORIASMAPAZONIFICACIONHIDROGRAFICA.pdf ↩ ↩2
-
https://www.cvc.gov.co/sites/default/files/Planes_y_Programas/Planes_de_Ordenacion_y_Manejo_de_Cuencas_Hidrografica/La%20Vieja%20-%20POMCA%20en%20Ajuste/Fase%20Diagnostico/7_CapituloI_Diagnostico_Morfometria.pdf ↩
-
http://geoservicios.ideam.gov.co/geonetwork/srv/eng/catalog.search#/metadata/7696695f-ae9c-4780-a6d0-d3cd1808819a ↩
-
http://geoservicios.ideam.gov.co/CatalogoObjetos/queryByUUID?uuid=bcd645c9-0f11-4770-926e-1e1fdfbf5ce6 ↩
-
https://www.igac.gov.co/sites/igac.gov.co/files/anexo_1.1_catalogo_objetos_cartografiabasica_v1.0_.pdf ↩