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

Можно ли установить тип данных (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' ...)

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

Источник

2 ответа

Решение

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

opts_str = '-ot Byte'
rast_opts = gdal.RasterizeOptions(options=opts_str)
Источник

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

Источник
Другие вопросы по тегам