Как работает инструмент Локальное понижение

Доступно с лицензией Spatial Analyst.

Локальное понижение – это ячейка или набор пространственно смежных ячеек, направлению стока из которых не может быть присвоено ни одного из восьми корректных значений направлений в растре направления стока. Такое происходит в том случае, если все соседние ячейки выше, чем обрабатываемая ячейка, либо когда две ячейки перетекают друг в друга, образуя петлю из двух ячеек.

Считается, что локальные понижения имеют неопределенное направление стока и им присваивается значение, которое представляет собой сумму их возможных направлений. Например, если самый крутой спуск и, следовательно, направление стока одинаковы как вправо (1), так и влево (16), значение 17 будет присвоено в качестве направления стока для этой ячейки.

Чтобы создать точное представление о направлении стока и, следовательно, суммарном стоке, лучше всего использовать набор данных, не содержащий локальных понижений. Цифровая модель рельефа (ЦМР), которая была обработана для удаления всех локальных понижений, называется ЦМР без понижений.

Локальные понижения в данных высот чаще всего возникают из-за ошибок в данных. Эти ошибки часто вызваны эффектами выборки и округлением высот до целых чисел. Естественные локальные понижения в данных высот с размером ячейки 10 метров и более встречаются редко (Mark 1988), за исключением ледниковых или карстовых областей, и обычно могут считаться ошибками. По мере увеличения размера ячейки количество локальных понижений в наборе данных также часто увеличивается.

Идентификация и удаление локальных понижений при создании ЦМР без понижений, представляет собой итеративный процесс. Когда локальное понижение заполнено, границы заполненной области могут создавать новые понижения, которые необходимо заполнить. Темы, касающиеся инструмента Заполнение зон, помогают объяснить, как создать ЦМР без понижений.

Иногда полезно знать глубину локального понижения или группы понижений. Эта информация может быть использована для определения подходящего значения ограничения по оси Z для инструмента Заполнение или для понимания типа ошибок, присутствующих в данных.

Вид профиля локального понижения
Вид профиля локального понижения

Пример

Процедура использования инструментов геообработки для заполнения локальных понижений следующая:

  1. Используйте Локальное понижение для создания растра понижений, кодированных глубиной.

    Входной растр направления стока : flowdir

    Выходной растр : sinks

  2. Используйте Водосборная область, чтобы создать растр участвующей области для каждого локального понижения.

    Входной растр направления стока : flowdir

    Входные растровые или векторные данные точек устьев : sinks

    Выходной растр : sink_areas

  3. Используйте Зональную статистику с минимальной статистикой, чтобы создать растр минимальной высоты в водоразделе каждого локального понижения.

    Входные растровые или векторные данные зон : sink_areas

    Поле зоны : Value

    Входной растр значений : elevation

    Выходной растр : sink_min

    Тип статистики : MINIMUM

  4. Используйте Заполнение зон, чтобы создать растр максимальной высоты в водоразделе каждого локального понижения.

    Выходной растр зон : sink_areas

    Входной растр весов : elevation

    Выходной растр : sink_max

  5. Используйте Вычитание, чтобы вычесть минимальное значение из максимального значения и найти глубину.

    Входной растр 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.

Связанные разделы


В этом разделе
  1. Пример
  2. База