Доступно с лицензией Spatial Analyst.
Локальное понижение – это ячейка или набор пространственно смежных ячеек, направлению стока из которых не может быть присвоено ни одного из восьми корректных значений направлений в растре направления стока. Такое происходит в том случае, если все соседние ячейки выше, чем обрабатываемая ячейка, либо когда две ячейки перетекают друг в друга, образуя петлю из двух ячеек.
Считается, что локальные понижения имеют неопределенное направление стока и им присваивается значение, которое представляет собой сумму их возможных направлений. Например, если самый крутой спуск и, следовательно, направление стока одинаковы как вправо (1), так и влево (16), значение 17 будет присвоено в качестве направления стока для этой ячейки.
Чтобы создать точное представление о направлении стока и, следовательно, суммарном стоке, лучше всего использовать набор данных, не содержащий локальных понижений. Цифровая модель рельефа (ЦМР), которая была обработана для удаления всех локальных понижений, называется ЦМР без понижений.
Локальные понижения в данных высот чаще всего возникают из-за ошибок в данных. Эти ошибки часто вызваны эффектами выборки и округлением высот до целых чисел. Естественные локальные понижения в данных высот с размером ячейки 10 метров и более встречаются редко (Mark 1988), за исключением ледниковых или карстовых областей, и обычно могут считаться ошибками. По мере увеличения размера ячейки количество локальных понижений в наборе данных также часто увеличивается.
Идентификация и удаление локальных понижений при создании ЦМР без понижений, представляет собой итеративный процесс. Когда локальное понижение заполнено, границы заполненной области могут создавать новые понижения, которые необходимо заполнить. Темы, касающиеся инструмента Заполнение зон, помогают объяснить, как создать ЦМР без понижений.
Иногда полезно знать глубину локального понижения или группы понижений. Эта информация может быть использована для определения подходящего значения ограничения по оси Z для инструмента Заполнение или для понимания типа ошибок, присутствующих в данных.
Пример
Процедура использования инструментов геообработки для заполнения локальных понижений следующая:
- Используйте Локальное понижение для создания растра понижений, кодированных глубиной.
Входной растр направления стока : flowdir
Выходной растр : sinks
- Используйте Водосборная область, чтобы создать растр участвующей области для каждого локального понижения.
Входной растр направления стока : flowdir
Входные растровые или векторные данные точек устьев : sinks
Выходной растр : sink_areas
- Используйте Зональную статистику с минимальной статистикой, чтобы создать растр минимальной высоты в водоразделе каждого локального понижения.
Входные растровые или векторные данные зон : sink_areas
Поле зоны : Value
Входной растр значений : elevation
Выходной растр : sink_min
Тип статистики : MINIMUM
- Используйте Заполнение зон, чтобы создать растр максимальной высоты в водоразделе каждого локального понижения.
Выходной растр зон : sink_areas
Входной растр весов : elevation
Выходной растр : sink_max
- Используйте Вычитание, чтобы вычесть минимальное значение из максимального значения и найти глубину.
Входной растр 1 : sink_max
Входной растр 2 : sink_min
Выходной растр : sink_depth
Это та же процедура с использованием Алгебры карт в 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)
База
Mark, D. M. 1988. "Network Models in Geomorphology." Modelling Geomorphological Systems, ed. M. G. Anderson. New York: John Wiley. 73–97.