Форумы GIS-Lab.info

У меня нет такого...

Версия 6.2.0?

Очень не хочется ставить сейчас gdal

QGIS не стоит?

В последнее время, чтобы получить channel network в SAGA я соблюдаю такую последовательность:

1. Sink drainage route detection (Terrain Analysis - Preprocessing)
2. Sink removal (с параметром depeen drainage routes) (Terrain Analysis - Preprocessing) (вот его результат используется дальше как исходный DEM)
3. Catchment areas (parallel), результат которого потом в channel network использую в качестве Initiation Grid (перед выполнением 4 пункта сам результат открываю и проверяю значения, определяю собственно порог значений (Initiation Threshold) - на сколько "маленькие" русла мне будут нужны, значение полностью зависят от ваших данных, адекватно использовать чье-то значение получится если только случайно). (Garden - Introducing Module Programming - 08...)
4. Channel Network с упомянутыми настройками

В общем-то, похожий подход описан у Дарьи Свидзинской в презентации, которую вам прислали.
У нас с вами, видимо, разные версии saga, но, надеюсь, получится попробовать.

С перепроецированием у SAGA действительно беда, проще воспользоваться QGIS.
На всех материалах получала адекватные результаты, но, к сожалению, сама не гидролог, может и есть косяки у такого подхода (в основном смущает sink removal).

Пока что обрабатываю перепроецированный растр (уже несколько часов кажется).
Переконвертировал с помощью GDAL в grd, единицы нормальные, сделал из него sgrd, если уж SAGA так их любит.
Я уж подумал, может nodata мешает? Правда я взял метод Planchon, Darboux, как в вебинаре. До этого использовал Wang&Liu
Вот что происходит:

Вы похоже с разрешением 10 метров или типа того наинтерполировали, это перебор. У SRTM порядка 50 метров нужно ставить ...

Хорошо, тогда скажите пожалуйста, где в GUI выставляется размер пиксела?
Или каков нормальный вид командной строки, потому что gdalwarp.exe -tr -t_srs "EPSG:32643" путькисх. путькрезульт
не работает.

Давайте по порядку, т.к. вы не ответили на несколько вопросов, которые вам задавали ранее:

1. Какая версия SAGA у вас установлена
2. Стоит ли QGIS и если да, то какой версии

Если хотите через командную строку, то изучите параметры gdalwarp http://www.gdal.org/gdalwarp.html
вам нужны s_srs, t_srs, r и tr как минимум

SAGA 2.3.2
QGIS 18.15

Я сейчас перепроецировал с разрешением 30 м, но уже через GDAL (в тот раз через QGIS).

но уже через GDAL (в тот раз через QGIS)

это одно и тоже

Не совсем. в QGIS нельзя указать размер ячейки.

Ну вот, опять пошли 88% бесконечные...
Кто-нибудь может сказать, сколько этих итераций обычно бывает?

QGIS использует gdalwarp для перепроецирования и размер ячйки там указать можно - нажимаете на карандашик и редактируете строку. Если вызывать перепроецирование через Processing, то там эти параметры вынесены отдельно.

Итераций может быть много - все зависит от размера растра. Попробуйте воспользоваться более новой версией SAGA - https://sourceforge.net/projects/saga-g ... A%20-%206/ уже 6 пошла.
Ну или выложите кусок вашего растра

Запустил с новой версией.
Растр здесь: https://yadi.sk/d/R1nqh_tJ3RVdgx

Ну вот, методом Wang&Liu получился нулевой растр!

Может быть формат нужен определённый? Я думал, какой читается тот и годится.

теперь у меня (из Украины) нет доступа к сервисам yandex'а - попробуйте dropbox или что-то еще
На всякий случай - лучше давать исходный растр GeoTIFF в WGS 84