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' ...)

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


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

Ответы (2)


arrow_upward
3
arrow_downward

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

person Benjamin    schedule 07.02.2017

arrow_upward
3
arrow_downward

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

opts_str = '-ot Byte'
rast_opts = gdal.RasterizeOptions(options=opts_str)
person Johan    schedule 08.02.2017