Disponible con una licencia de Spatial Analyst.
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 de mayor pendiente y, por lo tanto, la dirección del flujo es la misma tanto a la derecha (1) como a la izquierda (16), se asignará el valor 17 como la dirección del flujo para esa celda.
Para crear una representación precisa de la dirección del flujo y, por lo tanto, flujo acumulado, es mejor 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.
La presencia de sumideros en los datos de elevación suele deberse a errores en los datos. Estos errores suelen ser causados por efectos de muestreo y por el redondeo de 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 mayor son poco frecuentes (Mark 1988), excepto en áreas glaciales o karst, y se pueden considerar errores por lo general. A medida que aumenta el tamaño de celda, la cantidad de sumideros de un dataset también aumenta frecuentemente.
La identificación y remoción de sumideros cuando crea un DEM sin depresiones sin depresión es un proceso iterativo. Cuando se rellena un sumidero, los límites del área rellenada pueden crear nuevos sumideros que se deben rellenar. Los temas relacionados con la herramienta Relleno zonal explican cómo crear un DEM sin depresiones.
En ocasiones, resulta útil conocer la profundidad de un sumidero o grupo de sumideros. Esta información se puede utilizar para determinar un valor apropiado para establecer el límite Z para la herramienta Rellenar o para comprender el tipo de errores presentes en los datos.
Ejemplo
A continuación, se presenta un procedimiento que utiliza las herramientas de geoprocesamiento para rellenar los sumideros:
- Utilice Sumidero para crear un ráster de sumideros codificado con profundidad.
Ráster de dirección del flujo de entrada: flowdir
Ráster de salida: sinks
- Utilice Cuenca hidrográfica para crear un ráster del área de cuenca contribuyente para cada sumidero.
Ráster de dirección del flujo de entrada: flowdir
Datos de puntos de fluidez de ráster o entidad de entrada: sinks
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: value
Ráster de valores de entrada: elevation
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 zona de entrada: sink_areas
Ráster de peso de entrada: elevation
Ráster de salida: sink_max
- Utilice Resta para sustraer el valor mínimo del valor máximo para encontrar 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 utiliza Á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.