Как работает инструмент Направление стока (Flow Direction)

Для получения гидрологических характеристик поверхности является возможность определить направление стока из каждой ячейки растра. Это делается с помощью инструмента Направление стока (Flow Direction).

Этот инструмент использует поверхность в качестве входных данных и выдает растр, показывающий направление стока каждой ячейки. Если выбрана опция Выходной растр понижения (Output drop raster), создается выходной растр, показывающий отношение максимального изменения по высоте из каждой ячейки вдоль направления стока к расстоянию между центрами ячеек, выраженное в процентах. Если выбрана опция Сток из крайних ячеек направлен наружу (Force all edge cells to flow outward), сток из всех ячеек на ребре растра поверхности будет осуществляться наружу, за края растра.

Есть 8 корректных выходных направлений относительно 8 смежных ячеек, в которые может перейти сток. Этот подход обычно называется моделью восьминаправленного стока (D8), он следует принципам подхода в Jenson and Domingue (1988).

Присвоение кода направлению стока

Вычисление направления стока

Направление стока определяется направлением наиболее крутого спуска, или максимального понижения, из каждой ячейки. Оно вычисляется следующим образом:

 maximum_drop = change_in_z-value / distance * 100

Расстояние вычисляется между центрами ячеек. Следовательно, если размер ячейки принять за единицу, расстояние между двумя ортогональными ячейками будет равно 1, а расстояние между диагональными ячейками - 1,414 (квадратный корень из 2). Если максимальное понижение высоты до ближайших ячеек одинаково в нескольких направлениях, область соседства расширяется до тех пор, пока не будет найден самый крутой спуск.

Если найдено направление самого крутого понижения, выходной ячейке дается значение, представляющее это направление.

Если все соседние ячейки выше, чем обрабатываемая ячейка, такая ячейка будет рассматриваться как ошибка в данных; она должна быть заполнена до минимального значения высоты соседних ячеек. Сток будет осуществляться в эту ячейку. Однако в том случае, если локальное понижение размером в одну ячейку расположено на физическом краю растра, или в ее окрестностях есть хотя бы одна ячейка со значением NoData (нет данных), заполнения не происходит из-за недостаточной информации по соседним ячейкам. Чтобы ячейка могла рассматриваться как истинное локальное понижение размером в одну ячейку, для нее должна быть информация по всем соседним ячейкам.

Если сток из двух ячеек осуществляется друг в друга, они являются локальными понижениями с неопределенным направлением стока. Такой метод определения направления стока по цифровой модели рельефа рассмотрен в книге Дженсона и Доминика (Jenson and Domingue, 1988).

Ячейки, которые являются локальными понижениями, определяются с помощью инструмента Локальное понижение (Sink). Чтобы получить точное представление направления стока по поверхности, нужно заполнить локальные понижения до использования растра направления стока.

Ссылки

Greenlee, D. D. 1987. «Raster and Vector Processing for Scanned Linework». Photogrammetric Engineering and Remote Sensing 53 (10): 1383–1387.

Jenson, S. K., and J. O. Domingue. 1988. "Extracting Topographic Structure from Digital Elevation Data for Geographic Information System Analysis." Photogrammetric Engineering and Remote Sensing 54 (11): 1593–1600.