-->

19 marzo 2011

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


 
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.















 

24 comentarios :

  1. es muy bueno amigo pero queria preguntarte una cosa tenes tutoriales para clasificar areas boscosas con imagens lansat digamos trabajando con las imagenes clasificar todas sus areas cuantas hectarias hay de cultivo de pasto de arboles bajos, medios, altos para clasificar las enfermedades de los cultivos y etc. si tienes tutoriales de eso te lo agradeseria mucho, cuate. y gracias este es mi correo toro_elpropio@hotmail.com

    ResponderEliminar
  2. Excelente post, super util y super claro!! Saludos!

    ResponderEliminar
  3. Excelente, realmente es muy bueno tu Blog, felicitaciones y gracias por tus aportes.
    Un afectuoso saludo.

    ResponderEliminar
  4. Gracias amigo, como siempre excelente aporte

    ResponderEliminar
  5. muy buen aporte pero como simpre ami me pasa qeu ejecuto todos los pasoss y al final algo no me sale en este caso lo hice todo como decia, pero en la parte donde debo darle fill y propiedades sigo las instruciones pero donde me dice curren workpace le doy en la carpeta para guardar y le doy un nombre y me dice que the object cuold not be fould y le cambio el nombre y siempre me sale ese mensaje y no lo pude guardar, que habra pasado, por favor ayudame, gracias.

    ResponderEliminar
    Respuestas
    1. Hola Pattry, es confuso lo que sucede, pero fijate bien si en la ventana te pregunta que señales una carpeta no que guardes un archivo. solo pones el cursor sobre la carpta y le das add o guardar...por favor me indicas si te sirvió la respuesta

      Eliminar
  6. Hola aplicare tu consejo y te escribo. Y muchas gracias realmente tus aportes son muy importantes y oportunos estes pendientes de estas publicaciones pues te estare molestando.

    ResponderEliminar
  7. Hola ya creo que se que me paso en el paso 8 cuando le doy clic Feature Vertice To Points y le agrego end se deben ver como con una especie de sombra las elipses verdes y en mi caso no pasa nada, yo creo que desde ahi viene el problema bueno creia parece que ese es otro inconveniente, pues continue asi para ver si me guardaba y le di como me dijistes y si se deja guardar pero al dar ok no se ejecuta el modelo, como te dije el inconveniente qeu tengo en el paso 8 yo creo que por eso no se ejeccuta, ahora si no se que hacer lo repeti varias veces y me da igual, y ahora quien podraa salvarme, tuuu.

    ResponderEliminar
  8. Hola nuevamente yo segui intentandolo y me paso lo mismo el el paso 8 aun asi continue con el proceso lo guarde le di ok y nada entonces le di la fecha de run y me corrio el modelo, entonces los rastes me quedaron guardados en una carpeta que yo determine, no se si ese procedimiento que ejecute estuvo bien, ahora me surgio una duda, mira como ese modelo ModelBuilder se le pueden cambiar los parametros cada vez que uno lo requiera, pero yo abri artoolbox y alli no aparece, entonces como se le van a cambiar los parametros, pues en la carpeta como te dije si quedaron guardados los rarters, pero el ModelBuilder como tal no lo veo, seria por el procedimiento que le hice explicame, hasta pronto y muchas gracias.

    ResponderEliminar
  9. hola Peter puedes por favor subir los datos para poder practicarlo mejor?

    ResponderEliminar
  10. Muy pero muy bueno el post! Una explicación muy sencilla para realizar algo que en principio parecería muy complicado. Excelente!

    ResponderEliminar
  11. se puede hacer este poceso para delimitar todas las cuencas de un pais

    ResponderEliminar
  12. tengo un problema al delimitar, solo me sale una raya desde mi punto de salida. el area es llano

    ResponderEliminar
  13. Este comentario ha sido eliminado por el autor.

    ResponderEliminar
  14. Muy buen proceso de automatizacion de cuanca, te felicito y te agradesco es lo que estaba buscando creo que seria util que en l aparte final se agrega una ultima herramienta para la conversion de raster a poligono.

    gracias y buen blog. felicidades

    ResponderEliminar
  15. Muy buen proceso de automatizacion de cuanca, te felicito y te agradesco es lo que estaba buscando creo que seria util que en l aparte final se agrega una ultima herramienta para la conversion de raster a poligono.

    gracias y buen blog. felicidades

    ResponderEliminar
  16. Saludos
    En argis 10.1 en map algebra solo esta la opcion de raster calculator y no puedo conectarlo en el model, cual es el procedimiento en este caso?
    Gracias

    ResponderEliminar
  17. hola tengo el mismo problema
    En argis 10..2.2 en map algebra solo esta la opcion de raster calculator y no puedo conectarlo en el model, cual es el procedimiento en este caso?
    Gracias

    ResponderEliminar
  18. hola ! muy buen post !!! me esta siendo de mucha ayuda. Yo llegue al punto donde tengo que darle la condición a fill y voy a propiedades pero no me aparece la opción de "general Setting " por que puede ser ?

    ResponderEliminar

Tu comentario es muy importante...¿Cómo te pareció el tema?