Гидрологическая обработка рельфа
- Jasen
- Активный участник
- Сообщения: 100
- Зарегистрирован: 27 янв 2006, 18:33
- Репутация: 4
- Ваше звание: Специалист
- Откуда: Москва
- Контактная информация:
Гидрологическая обработка рельфа
Всем добрый день!
Поскольку жизнь заставила приобщиться теперь ещё и к 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.
Заранее спасибо всем за помощь и участие.
Поскольку жизнь заставила приобщиться теперь ещё и к 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.
Заранее спасибо всем за помощь и участие.
Что написал - то написал!
- darsvid
- Активный участник
- Сообщения: 163
- Зарегистрирован: 29 июн 2012, 12:40
- Репутация: 115
- Откуда: Kyїv, Ukraine
- Контактная информация:
Re: Гидрологическая обработка рельфа
Вот презентация - там пошагово все расписано
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/
Еще важный момент - из вашего описания непонятно - перепроецирован ли растр из географической системы координат в спроектированную. Если нет - то это нужно сделать
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/
Еще важный момент - из вашего описания непонятно - перепроецирован ли растр из географической системы координат в спроектированную. Если нет - то это нужно сделать
-
- Гуру
- Сообщения: 4067
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1062
- Ваше звание: программист
- Откуда: Казань
Re: Гидрологическая обработка рельфа
я бы посоветовал посмотреть на TAS / WhiteBox (http://www.uoguelph.ca/~hydrogeo/Whitebox/), и предварительно "прожечь" артефакты (у него в меню есть пункт).
- Jasen
- Активный участник
- Сообщения: 100
- Зарегистрирован: 27 янв 2006, 18:33
- Репутация: 4
- Ваше звание: Специалист
- Откуда: Москва
- Контактная информация:
Re: Гидрологическая обработка рельфа
Я понимаю, этот модуль генерирует не только filled sinks
Что написал - то написал!
- Jasen
- Активный участник
- Сообщения: 100
- Зарегистрирован: 27 янв 2006, 18:33
- Репутация: 4
- Ваше звание: Специалист
- Откуда: Москва
- Контактная информация:
Re: Гидрологическая обработка рельфа
Пожалуйста, дайте рабочую ссылку, тут пишут, что презентация перемещена.darsvid писал(а): ↑12 янв 2018, 14:38Вот презентация - там пошагово все расписано
https://www.slideshare.net/DariaSvidzin ... a-32605566
Я смотрел англоязычные презентации, но они не помогли.
Что написал - то написал!
- Jasen
- Активный участник
- Сообщения: 100
- Зарегистрирован: 27 янв 2006, 18:33
- Репутация: 4
- Ваше звание: Специалист
- Откуда: Москва
- Контактная информация:
Re: Гидрологическая обработка рельфа
Нет, все в исходных координатах SRTM, то есть WGS географический.
А точно надо? Почему тогда другие инструменты работали нормально? Вообще, не очень понимаю, почему, если Flow Accumulation нормально считается с таким растром, не работает Channel Network.
Что написал - то написал!
- Jasen
- Активный участник
- Сообщения: 100
- Зарегистрирован: 27 янв 2006, 18:33
- Репутация: 4
- Ваше звание: Специалист
- Откуда: Москва
- Контактная информация:
Re: Гидрологическая обработка рельфа
Я попробовал ещё Garden - Introducing Module Programming - 08: тоже не помогло в качестве Initiation Grid.
А из выпавшего у меня списка это какой инструмент, стало быть? (см. приложенную картинку)
- Вложения
-
- saga 1.png (70.25 КБ) 15479 просмотров
Что написал - то написал!
- darsvid
- Активный участник
- Сообщения: 163
- Зарегистрирован: 29 июн 2012, 12:40
- Репутация: 115
- Откуда: Kyїv, Ukraine
- Контактная информация:
Re: Гидрологическая обработка рельфа
Странно, у меня ссылка на презентацию открывается без проблем, но на всякий случай еще раз презентация
Перепроецировать надо
Перепроецировать надо
- Jasen
- Активный участник
- Сообщения: 100
- Зарегистрирован: 27 янв 2006, 18:33
- Репутация: 4
- Ваше звание: Специалист
- Откуда: Москва
- Контактная информация:
Re: Гидрологическая обработка рельфа
@darsvid, все понятно, на территории РФ сайт этот недоступен, как и сайт linkedin.
Что написал - то написал!
- Jasen
- Активный участник
- Сообщения: 100
- Зарегистрирован: 27 янв 2006, 18:33
- Репутация: 4
- Ваше звание: Специалист
- Откуда: Москва
- Контактная информация:
Re: Гидрологическая обработка рельфа
Нашел презентацию в Ютубе:
https://www.youtube.com/watch?v=IDRs8d6Ma6k
https://www.youtube.com/watch?v=IDRs8d6Ma6k
Что написал - то написал!
- darsvid
- Активный участник
- Сообщения: 163
- Зарегистрирован: 29 июн 2012, 12:40
- Репутация: 115
- Откуда: Kyїv, Ukraine
- Контактная информация:
Re: Гидрологическая обработка рельфа
А, понятно - тогда прикрепляю файл к сообщению
- Вложения
-
- Гидрологический анализ ЦМР в Открытой ГИС SAGA.pdf
- (1.84 МБ) 1038 скачиваний
- Jasen
- Активный участник
- Сообщения: 100
- Зарегистрирован: 27 янв 2006, 18:33
- Репутация: 4
- Ваше звание: Специалист
- Откуда: Москва
- Контактная информация:
Re: Гидрологическая обработка рельфа
А чем перепроэцировать? Proj4 я не нашел.
Что написал - то написал!
-
- Гуру
- Сообщения: 5298
- Зарегистрирован: 09 апр 2010, 19:30
- Репутация: 1015
- Ваше звание: просто мимо прохожу
- Откуда: Ё-бург
Re: Гидрологическая обработка рельфа
кому это мешало?на территории РФ сайт этот недоступен, как и сайт linkedin.
GDALА чем перепроэцировать? Proj4 я не нашел.
- Jasen
- Активный участник
- Сообщения: 100
- Зарегистрирован: 27 янв 2006, 18:33
- Репутация: 4
- Ваше звание: Специалист
- Откуда: Москва
- Контактная информация:
Re: Гидрологическая обработка рельфа
а через SAGA это никак нельзя сделать? proj4 раздела больше нет. Очень не хочется ставить сейчас gdal
Что написал - то написал!
- rhot
- Гуру
- Сообщения: 1727
- Зарегистрирован: 25 янв 2011, 17:50
- Репутация: 194
- Ваше звание: доктор
- Откуда: Архангельск
Re: Гидрологическая обработка рельфа
Geoprocessing > Projection > Coordinate transformation Grid
Последний раз редактировалось rhot 12 янв 2018, 17:52, всего редактировалось 1 раз.
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 0 гостей