Bezdepbonusu7Программирование: статьи и ответы
Меню

Вопрос

Gdal 2.1 Растеризация

Автор вопроса: JohanОтветов: 2

Можно ли установить тип данных (Byte, Float32) при использовании gdal.Rasterize в gdal 2.1? В настоящее время я использую gdal.Translate для преобразования в байт, но это неэффективно.

tif = my/target.tif
shp = my/source.shp
tiftemp = my/temp/solution.tif
rasterizeOptions = gdal.RasterizeOptions(xRes=20, yRes=20, allTouched=True etc.)
gdal.Rasterize(tiftemp, shp, options=rasterizeOptions)
#translate to Byte data type (not supported by Rasterize?)
gdal.Translate(tif, tiftemp, outputType=gdal.GDT_Byte,
                       creationOptions=['COMPRESS=PACKBITS')

Я знаю, что можно использовать

subprocess.check_call('gdal_rasterize', '-ot', 'byte' ...)

но я бы предпочел избежать этого, если это возможно. Любые идеи?

  • Пока tifftemp - набор данных, в который вы растеризуете, имеет правильный тип, все должно быть в порядке. Benjamin 19.07.2016
  • gdal.RasterizeOptions(outputType=gdal.GDT_Byte) вызывает TypeError: RasterizeOptions() получил неожиданный ключевой аргумент 'outputType' Johan 22.07.2016

pythongdalrasterize

Ответы

Ответ 1 · рейтинг 3 · Benjamin

Это было исправлено в GDAL 2.1.3. Ранее тип данных был жестко запрограммирован и не отображался в привязках Python. См. https://trac.osgeo.org/gdal/ticket/6710.

Ответ 2 · рейтинг 3 · Johan

Если вы используете gdal 2.1.2, используйте gdal.RasterizeOptions(), как если бы вы использовали его в командной строке:

opts_str = '-ot Byte'
rast_opts = gdal.RasterizeOptions(options=opts_str)

Перевод вопроса и ответов сообщества Stack Overflow, лицензия CC BY-SA. Оригинал на Stack Overflow.

Похожие материалы