Грид уклонов из инструмента Flow Direction

[_lenik]

Инструментом Spatial Analyst «Направление стока» (Flow Direction) помимо самого грида направления стока (а по сути этот грид отображает направление максимального уклона из ячейки) можно создать еще и грид того самого максимального уклона {out_drop_raster}. Уклон рассчитывается по простой формуле - разница высот ячеек/расстояние между центрами ячеек*100.

При использовании грида рельефа из данных SRTM грид максимальных уклонов выглядит как положено: максимальные уклоны на крутых склонах, минимальные - на равнинных участках, значения плавно меняются. Для проверки я в некоторых местах вручную пересчитывал – все сходится, как заявлено в справке.


Но при использовании грида рельефа интерполированного по точкам наблюдается следующая картина:
Грид рельефа (в данном случае сплайн, но пробовал и ОВР и топогрид – результат тот же)перепад высот 25 м, размер ячейки 1x1 м, проекция Гаусса-Крюгера:

relief_spline.png (57.42 КБ) 5981 просмотр

Грид направления стока:

flow_direction.png (282.55 КБ) 5981 просмотр

Грид максимальных уклонов по идее должен выглядеть примерно так:

а выглядит так:

drop.png (137.92 КБ) 5981 просмотр

Т.е. имеет всего 3 значения: 0, 33,3333 и 50. А если сделать размер ячейки 2 м, то 16,666 и 25, а если 10 м — 3,3333, 5, 6,6666 и 10 и т.д.
Если кто-нибудь может прояснить суть явления, буду бесконечно признателен.

[swdartvader]

А Вы не пробовали триангрулировать существующие измерения и посмотреть уклоны на триангулированной поверхности?
Я могу предложить Вам использовать инструмент Slope из набора Spatial Analyst, предназначеный для вычисления уклона, как в градусах, так и в % падения. Результаты получаются аналогичные изображенным при триангулировании.

А вычислять максимальные из них уже можно растровой алгеброй


А по-существу я плохо понимаю, как явление максимального уклона может быть непрерывна, он же не повсеместно максимален. (это к одному из Ваших изображений-приложений)

[_lenik]

А Вы не пробовали триангрулировать существующие измерения и посмотреть уклоны на триангулированной поверхности?
Я могу предложить Вам использовать инструмент Slope из набора Spatial Analyst, предназначеный для вычисления уклона, как в градусах, так и в % падения. Результаты получаются аналогичные изображенным при триангулировании.

Ну собственно 4-й рисунок получен инструментом Slope, но я хочу разобраться почему дополнительный грид уклонов из инструмента Flow Direction получается таким странным при использовании интерполированного грида

[swdartvader]

Выложите кусочек исходников, мб так понятнее станет мозгу улья)

[_lenik]

Вот например

[swdartvader]

Из растра spline1m я извлек мультиточки и построил по ним TIN, посмотрел уклоны, все получается хорошо (рис 1). Также хорошо получаются и уклоны инструментом Slope (рис 2)

1.jpg (350.32 КБ) 5877 просмотров

2.jpg (331.04 КБ) 5877 просмотров

В инструменте Flow Direction растр снижения показывает отношение максимального изменения по высоте из каждой ячейки вдоль направления стока к расстоянию между центрами ячеек, выраженное в процентах. Что-то мне подсказывает, что ключевую роль здесь играет фраза "вдоль направления стока", т.е., это не совсем то же самое что и классический уклон, однако четкую формулировку разницы я и не могу привести, буду еще думать

[_lenik]

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

Если вам несложно - можете по моему растру spline1m построить растры направления стока и снижения и показать что получилось? Может что с моим ArcGISом.

А еще лучше, можете по этим точкам интерполировать грид - может на этой стадии проблема.