Crear un DEM sin depresión

Un modelo digital de elevación (DEM) libre de sumideros, un DEM sin depresión, es la entrada deseada al proceso de dirección de flujo. La presencia de sumideros puede ocasionar un ráster de dirección de flujo con error. En algunos casos, pueden haber algunos sumideros legítimos en los datos. Es importante comprender la morfología del área lo suficientemente bien para saber qué entidades pueden ser realmente sumideros en la superficie de la tierra y cuáles son solamente errores en los datos.

Identificar sumideros

Los sumideros se pueden ubicar utilizando la herramienta Sumidero. Esta herramienta necesita un ráster de dirección que crea la herramienta Dirección de flujo. El resultado es un ráster que identifica cualquier sumidero existente en los datos. Según los resultados, puede rellenar los sumideros, o puede utilizar la salida para ayudar a determinar el límite de relleno. Los sumideros se pueden rellenar con la herramienta Relleno. Para utilizar la salida del Sumidero para determinar el límite de relleno, consulte "Búsqueda de profundidad de sumidero" en este tema (debajo).

Vista de perfil de un sumidero
Vista de perfil de un sumidero

Rellenar sumideros

La herramienta Rellenar utiliza una variedad de herramientas de ArcGIS Spatial Analyst, incluidas algunas de las herramientas de análisis hidrológico que se analizaron antes, para crear un DEM sin depresión. Esta herramienta requiere una superficie de entrada, un límite de relleno y un ráster de salida. Cuando se llena un sumidero, se rellena a su punto de fluidez, la elevación mínima a lo largo del límite de la cuenca hidrográfica.

La identificación y remoción de sumideros cuando crea un DEM sin depresión es un proceso iterativo. Cuando se rellena un sumidero, los límites del área rellenada pueden crear sumideros nuevos, que después deben llenarse. Para un DEM grande o uno con muchos sumideros, esto puede tardar desde varios minutos hasta horas.

Búsqueda de profundidad de sumidero

Es útil conocer la profundidad de un sumidero o un grupo de sumideros. Esta información se puede utilizar para determinar un límite z adecuado para la herramienta Rellenar, para comprender el tipo de errores que están presentes en los datos, y para determinar si los sumideros son entidades morfológicamente legítimas. Los pasos siguientes señalan el proceso general para buscar la profundidad del sumidero:

  1. Cree un ráster de sumideros con valores que identifican la profundidad al ejecutar la herramienta Rellenar para ubicar sumideros en el ráster.
  2. Utilice la herramienta Cuenca hidrográfica para crear un ráster del área de contribución para cada sumidero mediante la dirección de flujo desde el ráster de elevación y la salida desde la herramienta Sumidero como entrada para los puntos de fluidez.
  3. Con Álgebra de mapas en Python, utilice la herramienta Estadísticas zonales con la opción Mínima para crear un ráster de la elevación mínima en la cuenca hidrográfica de cada sumidero:
     sink_min = ZonalStatistics(sink_areas, "Value", elev_ras, "Minimum")
    La entrada sink_areas es la salida de la herramienta Cuenca hidrográfica.
  4. Cree un ráster que contenga la elevación más baja junto al límite de cada cuenca hidrográfica con la herramienta Relleno zonal (esto corresponde a la elevación en donde el flujo deja la cuenca después de rellenar hasta el borde):
     sink_max = ZonalFill(sink_areas, elev_ras)
  5. Utilice la herramienta Menos para sustraer el valor mínimo del valor máximo con el fin de buscar la profundidad nuevamente:
     sink_depth = Minus(sink_max, sink_min)

La secuencia anterior también se puede completar mediante las herramientas de geoprocesamiento Sumidero, Cuenca hidrográfica, Estadísticas zonales, Relleno zonal y Menos.

Temas relacionados


7/11/2012