Delimitación de cuencas hidrograficas con hydrology de Spatial Analyst y Modelbuilder

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.