Расчёт бассейна Вычегды средствами свободного ПО

[rhot]

Вот что получилось.
Оранжевый полигон делал Александр по SRTM30, а зелёный полигон - это данные по бассейну С.Двины.

[gamm]

SRTM не подойдёт - исток и устье находятся севернее 60°.

вроде как вам по SRTM все уже сделали?

gamm, прошу прошения за мою невежественность, но что такое "ляптахи"?

а там видно, если вывести - поскольку рельеф строился по разному количеству снимков, то он не стыкуется в местах, где количество снимков менялось. Кстати, попробуйте по GTOTO30 построить, под ваш масштаб может и хватить. Только проекцию нужно правильно выбрать на такую большую область, может что-нибудь азимутальное, вдоль реки, чтобы искажений было поменьше.

дождаться немецкой ЦМР, сделанной по радарной съемке на весь "шарик"

Вот так новость,.. ссылку в студию!

попутал, у них модель геоида (гравиметрия). Вам не подойдет ...

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

Александр. А может технологию раскроете? Как в QGIS так сделать? (С помощью Sextante например или встроенным GRASS)

Хитростей тут особых нет. Технология, может, и не совсем корректная, зато работает. Надо статейку написать (постараюсь, как появится время).

Кратко:
-- импорт ЦМР в GRASS;
-- подготовка ЦМР (r.fill.dir);
-- получение растров направления стока и речной сети (r.watershed);
-- векторизуем речную сеть (r.to.vect);
-- смотрим на странице в Википедии примерные координаты устья;
-- выбираем на выделенной по ЦМР линии реки ближайшую к устью точку (вертекс);
-- используем растр направления стока и координаты выбранной точки для выделения бассейна реки (r.water.outlet);
-- векторизуем границы бассейна (r.to.vect).

[dime1]

gdb - закрытый формат Esri, если не ошибаюсь?.. Кто-нибудь может вытащить оттуда только водосбор Вычегды и перезалить в нормальный формат?

Не могу. "Your connection credentials are not trusted by our Firewall, please contact us directly"

TAS делаете как минимум breaching

gamm, Вам удалось понять что делает breaching?

Александр Мурый, нашел пару текстов про алгоритмы ArcHydro, добавил во вложения.
Особенно интересны и пригодны (в то же время просты) варианты дендритных (древовидных) бассейнов без бессточных областей (артефактов ЦМР - sinks, либо внутренних бассейнов endorheic basin) в вариантах:
а) когда неизвестно ничего и есть только DEM
б) когда есть синтетическая гидросеть (спорный вариант, ибо непонятно откуда взять более точную гидросеть)
в) когда известны расположения рек (?), как я понимаю, если вектор рек делается на заказ или оцифровуется вручную со спутниковых снимков

P.S. Снова "загорелся" вариантом почистить ASTER GDEM2 (просто подменив максимальные и минимальные значения из SRTM на ту же территорию (чтоб убрать дыры) и пройдясь "плавающим кругом") для того, что сравнить качество получаемой гидросети с SRTM4 (никто такого не пробовал делать?). Но снова почитав про "ляптахи", поугас.
P.S.S.: у SRTM30 plus разрешение пикселя 1 км., а у SRTM4 - 90 метров?

[_lenik]

Данные по бассейну Вычегды с сайта CCM (26 мб, шейп-файлы сегментов водотоков и элементарных водосборных участков):
http://freespace.by/download/ff4bfd2aba

описание данных здесь:
http://ccm.jrc.ec.europa.eu/documents/C ... 07_STD.pdf

[rhot]

_lenik,

Услуги FreeSpace доступны только для белорусских сетей. Ваша сеть не зарегистрирована в РУП "Белтелеком" как белорусская.

[_lenik]

упс, куда вам залить?

[_lenik]

Отсюда попробуйте
Скачать Вычегда.rar с WebFile.RU

[rhot]

Отсюда попробуйте
Скачать Вычегда.rar с WebFile.RU

Спасибо!

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

Спасибо, _lenik. Если что, теперь знаю, к кому обращаться по поводу данных CCM

Сделал границу бассейна Вычегды одним полигоном из шейп-файла бассейнов по данных CCM ("CATCHMENTS.shp").
P.S. Просто было интересно протестировать на этих данных скриптик для получения в GRASS внешней границы для серии полигонов.

[gamm]

gamm, Вам удалось понять что делает breaching?

да - "пропиливает" реки не имея векторной гидросети, и "форсирует" направление потока на площинах. Но добиться от автора описания алгоритма или ссылки на статью не удалось. А разбирать тысячи строк кода - стар я уже для такого, да и зачем ...

[dime1]

gamm, интересно. Своего рода подготовительный этап для получения гидрологически-корретной ЦМР.
А бреачинг нужно применять до Fill Sinks (присвоение локальным высотным депрессиям соседних значений) или после?

[rhot]

dime1 +1
Я б тоже не отказался услышать мнения об алгоритмах построения гидрологически-корретной ЦМР.

[gamm]

А бреачинг нужно применять до Fill Sinks (присвоение локальным высотным депрессиям соседних значений) или после?

вместо - никакие Fill Sinks не нужны, он их "пропилит". Мы много строили бассейнов (тысячи), и ни один метод даже в сравнение не идет.

[dime1]

Breaching depression

This tool can be used to pre-process a digital elevation model (DEM) prior to being used for hydrological analysis. It will remove topographic depressions and apparent dams in the DEM in order to enforce continuous flow. This is achieved by breaching through topographic barriers. As such, depression breaching provides an alternative to depression filling. Unlike depression filling, which has one unique solution for each depression, there may be several competing solutions for depression breaching. This tool uses a cost-distance criteria for deciding upon the breach target, i.e. the cell to which the tool will trench a connecting path, and for determining the breach path itself. As such, the trenched path will not necessarily follow a straight line, but rather, will follow a path of lowest cost, where the cost is determined by the amount of change made to the original DEM.

Параметры: maximum breach channel lenght (pixel) и maximum elevation decrement value.

Т.е. алгоритм может быть примерно следующим:
1) Breaching depression
2) Flow Direction grid (D8)
3) Flow Accumulation grid
4) Stream Definition grid (specify Flow Accumulation Grid threshhold)
5) Stream Segmentation grid (stream tributaries)
6) Drainage Line Processing (stream vector)
7) Catchment Grid Delineation
Catchment Polygon Processing
9) Adjoint Catcment Processing

После шага 5 можно перевести в вектор, "дорисовать"/"почистить" реки, конвертировать обратно в grid и продолжать.