gdal 2.1栅格化

Question

在gdal 2.1中使用gdal.Rasterize时是否可以设置数据类型（Byte，Float32）？目前，我使用 gdal.Translate 转换为 Byte，但效率很低。

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

但我希望尽可能避免这种情况。有什么想法吗？

Answer 1

这已在 GDAL 2.1.3 中修复。以前的数据类型是 hard-coded 并且没有在 Python 绑定中公开。参见 https://trac.osgeo.org/gdal/ticket/6710。

Answer 2

如果您使用的是 gdal 2.1.2，请像在命令行中一样使用 gdal.RasterizeOptions()：

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

gdal 2.1栅格化

gdal 2.1 Rasterize

python

gdal

rasterize