Форумы GIS-Lab.info

Можно ли в gdal конвертировать векторные данные в растровые?
Есть необходимость построить матрицу высот(mtw) исходя из заданных векторных данных (sxf).

gdal_grid, но mtw нельзя, gdal его не пишет

и это не "конвертация векторных данных в растровые", для этого gdal_rasterize, это "построение ЦМР, интерполяция"

Как отдельная утилита gdal_grid есть, а как метод какого-либо класса пока не реализован...

Это утверждение или вопрос?

Утверждение. Как я понял поддержка будет с версии 2.1. Пытаюсь пока с самой утилитой разобраться. пускаю так:
gdal_grid.exe -a invdist:power=2.0:smoothing=1.0 -of RMF -co "MTW=ON" -ot Float64 -l Not_Classified F39035.SXF F39035.mtw
Долго и усердно что-то считается, но в результате открываю файл, а он весь белый(

А что в ГИС Карта, почему не сделать?

[ Сообщение с мобильного устройства ]

Ну не хотят начальники использовать ГИС Карта...

Я так полагаю что у меня ничего не получается, потому что я использую не то zfield. Но не могу нигде не могу найти список возможных имен для данного параметра

Все имеется:
GDAL v.1.11
https://github.com/OSGeo/gdal/blob/1.11 ... d.cpp#L493
GDAL v.2.0
https://github.com/OSGeo/gdal/blob/trun ... n.cpp#L179

Вот еще можно почитать: http://www.gdal.org/grid_tutorial.html

Уже нашел. Но вопрос про zfield остается. Как получить список имен данного параметра из datasource?

Либо http://www.gdal.org/classOGRPoint.html# ... 98bb5f4075
либо любое поля типа real (на крайняк integer)

evgen_sobolev писал(а):gdal_grid.exe -a invdist:power=2.0:smoothing=1.0

а использовать нормальную программу для интерполяции (ту же SAGA), и результат конвертировать что не позволяет? боюсь, что результат invdist вашим начальникам не понравится, даже с эстетической точки зрения, не говоря о содержательной (но могу и ошибаться).

P.S. интерполяция требует некоторого понимания, что и зачем делается, а не просто нажатия "волшебной кнопки". Если нужна только красивая картинка, то подойдет кригинг с линейной вариограммой, MBA и tension splines.

я не про тип координаты я про то как узнать название поля в котором храниться координата Z. Вот к примеру тут: http://nextgis.ru/blog/gdal_grid_geometry_support/
дана строка:
gdal_grid -a invdist:power:1.5smoothing:1.0:radius1=500:radius2=500 -outsize 100 100 -of GTiff -ot Float64 -l isolines -zfield height c:\map.sxf c:\dem.tiff —config GDAL_NUM_THREADS ALL_CPUS
из которой видно, что поле в котором храниться координата Z называется heigh ( -zfield heigh)
Вот собственно и вопрос как мне узнать как у меня называется это поле?

evgen_sobolev писал(а): Вот собственно и вопрос как мне узнать как у меня называется это поле?

А почему вы думаете, что кто-то ответит, не видя ваших данных? Откройте данные в любой ГИС или используйте утилиту <ogrinfo>.

Спасибо за наводку orginfo, что-то не догадался. Все получается. Теперь буду возиться с настройкой параметров алгоритма. Такой еще вопрос, как в флаге слоев -l прописать все имеющиеся слои сразу, а не перечисляя?