Delimitar una cuenca hidrográfica Automáticamente con Modelbuilder

19 de mar. de 2011 5 min de lectura

En artículos anteriores se habló sobre la automatización de los procedimientos en ArcGIS mediante Modelbuilder y la delimitación de cuencas con hydrology de Spatial Analyst.

En este artículo combinaremos estas dos herramientas para realizar un modelo automatizado para delimitar esas cuencas a partir de un modelo digital de elevaciones, el cual se podrá seguir utilizando las veces que queramos.

..lo que necesitamos…

1. El insumo principal es el modelo digital del terreno que ya hemos dicho como se obtiene (descargarlo de un servido o interpolando datos puntuales).

2. Lo siguiente es agregar las sub-herramientas que vienen con Hydrology de la extensión Spatial Analyst, las cuales son las siguientes:

Herramienta
Descripción
Entradas necesarias
Archivo de salida
MDE
Insumo para el cálculo
 
MDE
Fill
Corrige el modelo digital de elevación
MDE
MDE corregido
Flow Direction
Calcula el raster de dirección de flujo
MDE corregido
Dirección de flujo
Flow Accumulation
Calcula el flujo acumulado en una celda
Dirección de flujo
Acumulación de flujo
Single Output Map Algebra (se agrega de Algebra de mapas)
Calcula el raster de corrientes a partir de la acumulación del flujo y la especificación de un número de celdas que queremos que conformen la corriente (esto fue explicado en más detalle en el artículo anterior).
Acumulación de flujo
Red de drenaje
Stream Link
Divide los drenajes creados, en segmentos continuos a partir del raster de corrientes y la dirección de flujo
Red de drenaje Dirección de flujo
Drenaje segmentado
Stream Order
Crea un raster del orden de las corrientes
Red de drenaje Dirección de flujo
Orden de los drenajes
Stream Feature
Crea un shape de drenajes a partir de la dirección de flujo y el raster de corrientes
Red de drenaje Dirección de flujo
Shape de la red de drenaje
Feature Vertice To Point (se agrega Data Management tools)
Permite determinar los puntos a partir del raster de corrientes
Shape de la red de drenaje
Vértices de los drenajes
Watershed
Permite la delineación de las cuencas
Dirección de flujo
Vértices de los drenajes
Cuencas

Paso 1.

Para empezar damos clic derecho sobre arcToolbox y seleccionar la opción New Toolbox, para agregar una nueva barra de herramientas  a la que le damos el nombre de ModelCuencas.

Damos clic derecho sobre modelCuencas y seleccionamos la opción New seguido de model, aparece la siguiente ventana…

Paso 2

Damos clic derecho sobre la ventana de ModelBuilder, elegimos la opción Create Variable

De la ventana que aparece seleccionamos y agregamos un Raster Dataset

El resultado es el siguiente…

Paso 3

En adelante, de ArcToolbox nos vamos a Spatial Analyst Tools para luego seleccionar Hydrology y arrastramos las herramientas: Fill, Flow Direction, Flow Accumulation, Stream Link y Watershed.

De Data Management tools nos vamos a Feature y arrastramos Feature Vertice To Point

De Spatial Analyst Tools nos vamos a Map Algebra y seleccionamos Single Output Map Algebra (Algebra de mapas)

El resultado será el siguiente.

Paso 4

Ahora, con la herramienta Add Conection conectamos las herramientas y salidas agregadas

Se verá de la siguiente forma…

Paso 5

Damos clic derecho sobre cada una de las entradas y las salidas (estoy hablando de las elipses); del menú que aparece seleccionamos las opciones Model Parameter y Add To Display, el resultado se muestra en la siguiente figura.

Como pueden ver, aparece una “P” junto a estos elementos que indica que son "parámetros" del modelo que se puede ajustar en tiempo de ejecución.

Paso 6

Ahora ingresamos modelo digital de elevaciones con que queremos trabajar y los parámetros en la herramienta Single Output Map Algebra. Procedemos de la siguiente forma…

Damos clic derecho sobre Raster layer y seleccionamos la opción Open.


Aparece una ventana en la que debemos buscar o asignar el Raster con el que queremos trabajar, en este caso asigné dtm_cuencas.

Una vez seleccionado el archivo raster que alimenta el modelo; algunas de las herramientas agregadas al modelo se activarán cambiando de color.

Paso 7

Para que las demás salidas y herramientas se activen, damos clic derecho sobre el cuadro de la herramienta llamada Single Output Map Algebra y seleccionamos la opción Open…, se mostrará una ventana en la que especificamos el tamaño mínimo de las celdas que abarca la cuenca, en este caso escribiremos Con([Output accumulation raster] > 20000, 1).

Nota: Output accumulation raster es la salida que se obtiene después de aplicar la herramienta Flow Accumulation, se debe tener cuidado con este dato, pues de lo contrario el modelo no funciona.

Damos clic en OK y obtenemos el resultado es el siguiente.

Paso 8

Damos clic derecho sobre el cuadro de la herramienta Feature Vertice To Points y seleccionamos la opción Open, en la ventana que aparece y en campo Point Type (optional) seleccionamos la opción END.

El modelo se verá de la siguiente forma.

Paso 9

Ahora le indicamos al modelo la ruta donde deben ser guardados los resultados, para ellos: Seleccionamos el cuadro de la herramienta Fill y damos clic derecho; seleccionamos la opción Properties…

Aparece la siguiente ventana, activamos la opción General Settings.

Después de dar clic en el botó Values…, se mostrará la siguiente ventana.

Seleccionamos la opción General settings, se desplegará una serie de opciones en las cuales debemos llenar Current Workpace y Scratch Workspace con la ubicación de la carpeta donde queremos guardar los resultados.

Paso 10

Finalmente ejecutamos el modelo…después de unos minutos… obtenemos el siguiente resultado.

Quizá parezca un poco laborioso construir el modelo, sin embargo la ventaja que nos proporciona es que lo podemos seguir utilizando cuantas veces queramos.

Al modelo se le podría cambiar el nombre de los archivos de salida, quedaría como se muestra.

Pedro Villegas

Ingeniero Civil, Master en Ingeniería con énfasis en hidrogeología