Карта горизонтальной расчлененности рельефа

[Ремесленник]

Уважаемые Господа!
В недавнем времени возник вопрос о возможности построения карты горизонтальной расчлененности рельефа используя GRASS и QGIS. В связи с эти хотелось бы спросить можно ли такое сделать используя только ГРАСС и Квантум Гис? Немного подробнее о сути проблемы. В ряде источников такая процедура выполняется в АркМап последовательностью Spatial Analyst>FlowDirection>FlowAccamulation>StremOrder затем Зональная статистика в квадратных ячейках (подсчитывается количество водотоков того или иного уровня (по Стралеру) на ячейку (например размером 2*2 километра), затем суммарная длина всех тальвегов на ту же ячейку (таким образом горизонтальная расчлененность получается км/км2 )). Итак вопрос : как получить такую карту используя только Грасс и КвантумГис? К стати r.fill.dir так и не заработал. А если точнее работает только на маленьких кусках размером 8-10 мегабайт. 100-мегабайтный AsterGDEM v2 обработать не удалось.
Заранее Благодарен за советы, подсказки и рекомендации.

[rhot]

можно ли такое сделать используя только ГРАСС и Квантум Гис?

Можно, используйте v.strahler или r.stream.order.

Вот примерная пошаговая инструкция: http://www.surfaces.co.il/?p=1321

[Александр Мурый]

Тема требует подробного разбора, после конференции попробую разобраться, как посчитать в GRASS горизонтальную расчлененность рельефа.
Насчёт порядка водотоков по Штралеру почитайте эту тему.

[Ремесленник]

А как установить эти расширения ( v.strahler r.stream.order )? Их что надо компилировать?

[Александр Мурый]

А как установить эти расширения ( v.strahler r.stream.order )? Их что надо компилировать?

Попробуйте скачать отсюда скомпилированные версии этих модулей для Windows и поместить их в соответствующие каталоги.

[Ремесленник]

Установлю и будем экспериментировать.

[Александр Мурый]

Вот примерный порядок расчёта горизонтальной расчленённости рельефа в GRASS.
Возьмём упомянутый выше простой алгоритм, только (пока) без учёта порядков водотоков (считаем для всех).

1) Выделение речной сети с порядками по схеме Штралера.
Примерно так, как описано в этой теме. На картинке ниже сделано через r.stream.order.
2) Зональная статистика длин водотоков:
Можно пойти по схеме из этой темы.
3) Считаем коэфф. гориз. расчленённости через простой SQL-запрос по атрибутам квадратов.

В итоге имеем (сверху — речная сеть с порядками по Штралеру и сетка квадратов; снизу — "карта" коэффициентов гориз. расчленённости рельефа, до 0,35 км/км2).

relief_horiz_dissection.png (112.21 КБ) 10028 просмотров