Cómo funciona Sumidero

Un sumidero es una celda o conjunto de celdas conectadas espacialmente a cuya dirección de flujo no se le puede asignar uno de los ocho valores válidos en un ráster de dirección de flujo. Esto puede ocurrir cuando todas las celdas próximas son mayores que la celda que se está procesando o cuando dos celdas fluyen una dentro de la otra, creando un bucle de dos celdas.

Se considera que los sumideros tienen direcciones de flujo indefinidas y se les asigna un valor que es la suma de sus posibles direcciones. Por ejemplo, si la caída más empinada y, por lo tanto, la dirección de flujo son las mismas para la derecha (1) y la izquierda (16), se asignará el valor 17 como la dirección de flujo para esa celda.

Para crear una representación precisa de la dirección de flujo y, por consiguiente, flujo acumulado, se recomienda utilizar un dataset que esté libre de sumideros. Un modelo digital de elevación (DEM) que se procesó para quitar todos los sumideros se denomina DEM sin depresiones.

Los sumideros en los datos de elevación ocurren generalmente debido a errores en los datos. Comúnmente, estos errores se producen debido a los efectos de muestreo y al redondeo de las elevaciones a números enteros. Los sumideros que ocurren naturalmente en los datos de elevación con un tamaño de celda de 10 metros o más grande son poco frecuentes (Mark 1988), excepto en áreas glaciales o de karst, y generalmente se pueden considerar como errores. A medida que el tamaño de celda aumenta, la cantidad de sumideros en un dataset también suele aumentar.

La identificación y remoción de sumideros durante la creación de un DEM sin depresiones es un proceso iterativo. Cuando se rellena un sumidero, los límites del área que se rellenó pueden crear nuevos sumideros que deben rellenarse. Los temas acerca de la herramienta Relleno zonal explican cómo se crea un DEM sin depresiones.

A veces resulta útil conocer la profundidad de un sumidero o un grupo de sumideros. Esta información se puede utilizar para determinar un valor adecuado a fin de establecer el Límite Z para la herramienta Relleno o para comprender el tipo de errores presentes en los datos.

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

Ejemplo

A continuación, se presenta un procedimiento con herramientas de geoprocesamiento para rellenar sumideros:

  1. Utilice Sumidero para crear un ráster de sumideros codificado con la profundidad.

    Ráster de dirección de flujos de entrada: flowdir

    Ráster de salida: sinks

  2. Utilice Cuenca hidrográfica para crear un ráster del área de contribución para cada sumidero.

    Ráster de dirección de flujos de entrada: flowdir

    Datos de corrientes de ráster o entidad de entrada: sinks

    Ráster de salida: sink_areas

  3. Utilice Estadísticas zonales con la estadística mínima para crear un ráster de la elevación mínima en la cuenca hidrográfica de cada sumidero.

    Datos de zonas de entidad o ráster de entrada: sink_areas

    Campo de zona: Valor

    Ráster de valores de entrada: elevación

    Ráster de salida: sink_min

    Tipo de estadística: MÍNIMO

  4. Utilice Relleno zonal para crear un ráster de la elevación máxima en la cuenca hidrográfica de cada sumidero.

    Ráster de zona de entrada: sink_areas

    Ráster de pesos de entrada: elevación

    Ráster de salida: sink_max

  5. Utilice Resta para sustraer el valor mínimo del valor máximo con el fin de buscar la profundidad.

    Ráster de entrada 1: sink_max

    Ráster de entrada 2: sink_min

    Ráster de salida: sink_depth

Este es el mismo procedimiento que se utiliza con Álgebra de mapas en ArcPy:

sinks = Sink(flowdir sink_areas = Watershed(flowdir, sinks) sink_min = ZonalStatistics(sink_areas, "Value", elevation, "Minimum") sink_max = ZonalFill(sink_areas, elevation) sink_depth = Minus(sink_max, sink_min)

Referencia

Mark, D. M. 1988. "Network Models in Geomorphology." Modelling Geomorphological Systems, ed. M. G. Anderson. New York: John Wiley. 73–97.

Temas relacionados


7/11/2012