Cómo funciona Sumidero
Un sumidero es una celda o un 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 vecinas son mayores que la celda de procesamiento o cuando dos celdas fluyen entre sí, 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 |
Ejemplo
A continuación, se presenta un procedimiento con herramientas de geoprocesamiento para rellenar sumideros:
- Utilice Sumidero para crear un ráster de sumideros codificado con la profundidad.
Ráster de entrada de dirección de flujos: flowdir
Ráster de salida: sumideros
- Utilice Cuenca hidrográfica para crear un ráster del área de contribución para cada sumidero.
Ráster de entrada de dirección de flujos: flowdir
Datos de corrientes de ráster o entidad de entrada: sumideros
Ráster de salida: sink_areas
- 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 entrada de valores: elevación
Ráster de salida: sink_min
Tipo de estadística: MINIMUM
- 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 entrada de zona: sink_areas
Ráster de entrada de pesos: elevación
Ráster de salida: sink_max
- 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.