Артефакты в tiff (рельеф) и тонкости работы с растрами

[evgeny.shirinyan]

судя по картинке, там исходно были изолинии, и проглядываются структурные элементы (бровка обрыва, тальвеги, и т.д.). Если все это есть в векторном виде, то именно для рельефа методы есть, например в Сканексе. Общий смысл такой - строится системе прямых линий (например, по линиям будущего грида), определяются места их пересечения с изолиниями (и, соответственно, высоты) и места пересечений со структурными элементами. Затем между пересечениями со структурными элементами вдоль линий строятся одномерные сплайны, которые используются для сгущения точек. Когда точек станет очень много, используем МВА (есть в R и SAGA).

Исходные были просто точки. Видимо, съемка делалась более детально в этих местах, чтобы передать характер.
В векторе линии есть. В принципе, думаю, можно идти через Civil 3D с его функционалом.

Конечно, был бы признателен, если разместите ссылку на какой-нить туториал по этой теме для чайников.

[pendduduk]

Грид если строил, то делал его по исходным точкам, а не по TIN. Вопрос про TIN. Я подразумевал под TIN алгоритм внутри QGIS-функции "Интерполяция", и в результате получал растр (или грид - это не одно и то же ли?). Как мне казалось, что TIN-формат хранения данных - все же штука векторная (во всяком случае в ArcView TIN и грид - разные вещи). Или я заблуждаюсь?

TIN это способ представления поверхности в виде треулольников, которые строятся обычно методом треангуляции Делоне.

Боюсь что подсказать что то конкретное вам будет сложно... не зная вашей конечной цели и не видя исходных данных.

Попробуйте разные варианты интерполяции из SAGA (с параметрами по умолчанию) и выберите оптимальный.

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

В QGIS получается растеризованный TIN, в векторном виде можно вытащить только саму сеть треугольников в виде триангуляции Делоне без высотных отметок (т.е. почти бесполезная штука для рельефа).

В GRASS есть несколько проверенных способов построения ЦМР, всё зависит, как обычно от исходных данных. Если это изолинии (либо точки по изолиниям, как в данном случае), то лучше всего использовать спец. модуль r.surf.contour. Также есть модули v.surf.rst (регуляризованный сплайн с натяжением) и v.surf.bspline (тоже сплайны, лучше для лидарных данных); есть ещё дополн. модуль r.surf.nnbathy (интерполяция методом "естественного соседства" и линейная по типу TIN).

[evgeny.shirinyan]

TIN это способ представления поверхности в виде треулольников, которые строятся обычно методом триангуляции Делоне.

С SAGA пробовал, там, кстати, TIN идет особняком - явно не растр. Вообще, я довольно много работаю с 3D, и Делоне там хватает).
Параметры по умолчанию в SAGA в том же Шепарде давали странные результаты, поэтому я и решил оставить пока что до лучших времен.

[evgeny.shirinyan]

В QGIS получается растеризованный TIN, в векторном виде можно вытащить только саму сеть треугольников в виде триангуляции Делоне без высотных отметок (т.е. почти бесполезная штука для рельефа).

В GRASS есть несколько проверенных способов построения ЦМР, всё зависит, как обычно от исходных данных. Если это изолинии (либо точки по изолиниям, как в данном случае), то лучше всего использовать спец. модуль r.surf.contour. Также есть модули v.surf.rst (регуляризованный сплайн с натяжением) и v.surf.bspline (тоже сплайны, лучше для лидарных данных); есть ещё дополн. модуль r.surf.nnbathy (интерполяция методом "естественного соседства" и линейная по типу TIN).

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

[gamm]

С SAGA пробовал, там, кстати, TIN идет особняком - явно не растр. Вообще, я довольно много работаю с 3D, и Делоне там хватает).
Параметры по умолчанию в SAGA в том же Шепарде давали странные результаты, поэтому я и решил оставить пока что до лучших времен.

у вас какя-то странная SAGA ... вам нужен пункт Tools: Grid -Spline Interpolation (у меня v2.1.2), там практически все годное:
- Multilevel B-spline Interpolation (это тот самый MBA, там и ссылочка есть - он на самом деле Approximation, а не Interpolation, но при заданном по умолчанию числе уровней иерархии 11 одно от другого уже не отличается)
- Thin Plate Spline (global) - должно потянуть, точек не очень много
- Tin Plate Spline (TIN) - не совсем Ренка, но если задать побольше соседей (Level2) и поиграть с регуляризацией, должен сгладить переходы на границах треугольников

[evgeny.shirinyan]

gamm, SAGA у меня такая же 2.1.2. Спасибо за пояснение.

Кстати, почему-то при команде "Растр-Преобразовать формат" не удается использовать команды в -r bilinear, cubic итд. вот скрин -

[evgeny.shirinyan]

И еще один вопрос (хотел бы уточнить) - правильно ли я понимаю, что горизонтали зачастую являются производным продуктом от точек и сильно зависят от алгоритма, который использовался при создании горизонталей? И таким образом, горизонтали - лишь гипотеза, в отличие от точки, локального измерения?

[gamm]

а зачем картинка в дропбоксе - ее не видно. Есть же нормальный tag [img] для картинок ...

[evgeny.shirinyan]

Поправил исходное сообщ. Почему-то дб не показывает картинку.

[evgeny.shirinyan]

И еще один вопрос (хотел бы уточнить) - правильно ли я понимаю, что горизонтали зачастую являются производным продуктом от точек и сильно зависят от алгоритма, который использовался при создании горизонталей? И таким образом, горизонтали - лишь гипотеза, в отличие от точки, локального измерения?

Пообщался с одним дорожником, и тот мне потвердил, что это так. Т.е. полагаемся на точки. Со структурными линиями уже я разбирался в Civil3d. Для такой детализации это оптимальный инструмент, как мне кажется.