Форумы GIS-Lab.info

Hi all,

When I used Extract raster values at vector points location with Python/GDAL, I get error:



And when I opened point layer, new field created is zero




Can you help me?

Thanks

Поскольку спустя хз знает сколько лет, ни одна ГИС быстро и качественно так и не научилась экстрагировать данные из растров по точечным профилям, то расширил и углубил скрипт.

Теперь скрипт:
- умеет фильтровать растры в указанной папке по расширению
- умеет альтернативно использовать gdalallocinfo (полезно когда огромные растры не помещаются в память)

хочется придумать что делать с названиями растров которые не влезают в ограничение DBF, варианты:
- сохранять отдельно перекодировочную таблицу с вектором имен и кодов
- вообще забить на чтение/сохранение в shape, читать/сохранять в CSV/VRT

В гроб эти шейпы, CSV вполне бы хватило. Собственно, не зная о существовании этого скрипта я как-то такое писал под Windows и именно с CSV (cmd.exe в циклах for с ним только и умеет работать).

Максим Дубинин писал(а):ни одна ГИС быстро и качественно так и не научилась экстрагировать данные из растров

А как же GRASS? Он всегда это умел. Конечно растр придется импортировать в GRASS, а вот вывод можно организовать по всякому и в шейп и в CSV.

Семплировать растр по точкам иди координатам: v.what.rast или r.what (последний поддерживает семплирование сразу по множеству растров).

у CSV есть ряд существенных недостатков, но по хорошему нужно бы дописать блок работы с ним в скрипт тоже

GRASS хорош, но требует много подготовительных действий (одна загрузка огромных растров в область сожрет кучу времени), а хочется просто указать точки и растры и получить результат, что этот скрипт и делает.

Добавил ключ -f для экспорта названий полей в виде: Теперь в принципе легко результат втащить в R и изящно восстановить names в оригинальном виде, ну или где угодно еще это сделать: И поправил код для многоканальных растров, раза в три побыстрее стало работать.

Теперь можно собирать данные со списка растров
и со списка папок растров

ссылка на скрипт в статье выдает ошибку

gornak писал(а):ссылка на скрипт в статье выдает ошибку

Исправил.

архиватор 7z пишет- недопустимый метод сжатия , да и объем подозрительный 5 Кб , нет тестовых вложений(которые по -видимому подразумеваются), в итоге - распаковывается в 0 Кб

gornak писал(а):архиватор 7z пишет- недопустимый метод сжатия , да и объем подозрительный 5 Кб , нет тестовых вложений(которые по -видимому подразумеваются), в итоге - распаковывается в 0 Кб

Проверил — у меня всё нормально разархивируется разными архиваторами. 7z показывает, то архив целый.

тоже проверил, все норм, внутри 1 файл - скрипт на питоне

да прошу прощения!-давно 7z архиватор не обновлял

чем может быть вызвана ошибка?
установлен python27 и gdal к нему, простейшие скрипты работают

у вас GDAL 1.X, а я писал под двойку
удалите в скрипте ", gdal.OF_VECTOR" в 395 строке по идее должно заработать

самый простой способ получить GDAL 2.0 - поставить NGQ и запускаться из тамошнего шелла.