Страница 1 из 3

Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 13:05
Jasen
Всем добрый день!

Поскольку жизнь заставила приобщиться теперь ещё и к SAGA, иду в этот раздел и открываю тему по гидрологии в SAGA. Особенно не рассчитываю, что вы сразу ответите, поэтому пойду искать помощи и на других ресурсах, а если что найду полезное - в этой теме напишу.
Задача с гидрологической точки зрения простая: построить водосбор и речную сеть для хорошо расчлененного рельефа (Таджикистан, р. Вахш, гидропост Комсомолабад). В QGIS или GRASS наверняка мне бы удалось это сделать, но, ходят слухи, что SAGA строит точнее, поэтому была выбрана SAGA.

Что я попробовал (работаю из командной строки, конечно, иначе немыслимо долго).
1. Импортировал hgt 30-метрового SRTM (Import USGS SRTM Grid).
2. Сшил их в единый ГРИД: Grid - Grid System - Mozaicking
3. Убрал бессточные ячейки: Terrain Analysis - Preprocessing - Fill Sinks (Wang$Liu). Почему именно этот? Он сразу строит водосборы и строит их отлично, но сразу на весь ГРИД.
4. Далее попытался построить речную сеть через Terrain analisys - Channels - Channel Network. В качестве Initiation Grid использовал ГРИД с бассейнами, полученными на шаге 3. В качестве ГРИД Flow Directions использовал ГРИД полученный на шаге 3. Initiation threshold 1000000, хотя потом пробовал и 10000 (медленно, принудительно вышел), и 100000 - результат тот же. Результат: пустой грид (значения 0, 0).
5. Попробовал инструмент Terrain Analysis - Channels - Channel Network and Drainage Basins. Получились три шейпа: с небольшими водотоками без основных русел, с точками соединений русел и с водосборчиками. Напоминаю, что хотел я получить один большой водосбор по створу Комсомолабад и умеренно прореженную речную сеть на нем.
6. Узнал, что есть UpslopeArea (Terrain Analysis - Hydrology - Flow Accumulation - Upslope area). После просчёта получился странный грид с очень маленькой upslope areа, однако координаты были заданы правильно, потому что эта площадка прилегала к заданной точке. Точку я нанес на грид с помощью инструмента Shapes - Conversion - Convert Table to Points (из таб файла тхт).

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

Ещё раз сформулирую задачу: на основе файлов ЦМР получить векторные площадь водосбора и речную сеть (первые 3 - 4 порядка) для конкретного створа, координаты которого известны.

Для тех, кто захочет тоже попробовать, выложил свой FilledDEM на Яндекс: https://yadi.sk/d/FKuIKd9v3RPcH3
В нем самом проблемы точно нет, ибо Accumulation дает идеальные водотоки, а водосборы через FillSinks также рисуются замечательно. Значит, проблема в моем незнании инструментов SAGA.

Заранее спасибо всем за помощь и участие.

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 14:38
darsvid
Вот презентация - там пошагово все расписано
https://www.slideshare.net/DariaSvidzin ... a-32605566

Из вашего описания - есть 2 ошибки:
1. в качестве initation grid должна быть flow accumulation/catchment area
2. flow directions - это flow directions а не гидрокорректированная ЦМР

Есть еще хорошая штучка TAS Whitebox - там тоже гидрологический анализ хорошо сделан
https://whiteboxgeospatial.wordpress.co ... tebox-gat/

Еще важный момент - из вашего описания непонятно - перепроецирован ли растр из географической системы координат в спроектированную. Если нет - то это нужно сделать

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 14:51
gamm
я бы посоветовал посмотреть на TAS / WhiteBox (http://www.uoguelph.ca/~hydrogeo/Whitebox/), и предварительно "прожечь" артефакты (у него в меню есть пункт).

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 15:38
Jasen
darsvid писал(а):
12 янв 2018, 14:38
flow directions - это flow directions а не гидрокорректированная ЦМР
Я понимаю, этот модуль генерирует не только filled sinks :)

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 15:40
Jasen
darsvid писал(а):
12 янв 2018, 14:38
Вот презентация - там пошагово все расписано
https://www.slideshare.net/DariaSvidzin ... a-32605566
Пожалуйста, дайте рабочую ссылку, тут пишут, что презентация перемещена.
Я смотрел англоязычные презентации, но они не помогли.

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 15:42
Jasen
darsvid писал(а):
12 янв 2018, 14:38
Еще важный момент - из вашего описания непонятно - перепроецирован ли растр из географической системы координат в спроектированную. Если нет - то это нужно сделать
Нет, все в исходных координатах SRTM, то есть WGS географический.
А точно надо? Почему тогда другие инструменты работали нормально? Вообще, не очень понимаю, почему, если Flow Accumulation нормально считается с таким растром, не работает Channel Network.

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 15:46
Jasen
darsvid писал(а):
12 янв 2018, 14:38
в качестве initation grid должна быть flow accumulation/catchment area
Я попробовал ещё Garden - Introducing Module Programming - 08: тоже не помогло в качестве Initiation Grid.
А из выпавшего у меня списка это какой инструмент, стало быть? (см. приложенную картинку)

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 15:54
darsvid
Странно, у меня ссылка на презентацию открывается без проблем, но на всякий случай еще раз презентация

Перепроецировать надо

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 16:15
Jasen
@darsvid, все понятно, на территории РФ сайт этот недоступен, как и сайт linkedin.

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 16:18
Jasen
Нашел презентацию в Ютубе:
https://www.youtube.com/watch?v=IDRs8d6Ma6k

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 16:32
darsvid
А, понятно - тогда прикрепляю файл к сообщению

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 16:59
Jasen
darsvid писал(а):
12 янв 2018, 15:54
Перепроецировать надо
А чем перепроэцировать? Proj4 я не нашел.

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 17:07
trir
на территории РФ сайт этот недоступен, как и сайт linkedin.
кому это мешало?
А чем перепроэцировать? Proj4 я не нашел.
GDAL

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 17:33
Jasen
trir писал(а):
12 янв 2018, 17:07
GDAL
а через SAGA это никак нельзя сделать? proj4 раздела больше нет. Очень не хочется ставить сейчас gdal

Re: Гидрологическая обработка рельфа

Добавлено: 12 янв 2018, 17:50
rhot
Geoprocessing > Projection > Coordinate transformation Grid :?: