gdal 2.1 Rasteriser

Question

Est-il possible de définir le type de données (Byte, Float32) lors de l'utilisation de gdal.Rasterize dans gdal 2.1? Actuellement, j'utilise gdal.Translate pour convertir en Byte, mais c'est inefficace.

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

Je suis conscient qu'il est possible d'utiliser

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

mais je préférerais éviter cela si possible. Des idées?

Answer 1

Si vous êtes GDAL 2.1.2, utilisez gdal.RasterizeOptions() comme si vous utilisiez sur la ligne de commande:

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

Answer 2

Ceci a été corrigé dans GDAL 2.1.3. Auparavant, le type de données était codé en dur et non exposé dans les liaisons Python. Voir https://trac.osgeo.org/gdal/ticket/6710.