Annonce
Pour sa 21ème année, l’association GeoRezo a toujours besoin de vous !
10€ = 1 mois de frais bancaires ; 15€ = 12 mois de nom de domaine ; 30€ = 1 semaine de location des serveurs …
Retrouver nos membres bienfaiteurs
#1 Wed 15 November 2023 07:36
- nico-29
- Participant assidu
- Lieu: Brest
- Date d'inscription: 2 Jun 2006
- Messages: 185
GDAL Python et ColorTable
Bonjour,
Avec GDAL sous Python, je crée un raster geotiff à partir de plusieurs dalles raster en passant par les raster virtuel (ff contient la liste des dalles initiales, geotiff le chemin du raster final).
Code:
vrt_options = gdal.BuildVRTOptions(resampleAlg='cubic', addAlpha=True, outputSRS="EPSG:2154") vrt = gdal.BuildVRT('my.vrt', ff, options=vrt_options) gdal.Translate(geotiff, vrt)
ça ça marche, mais je voudrais ajouter une colortable prédéfinie dans le raster final. j'ai essayé ce code :
Code:
stats = vrt.GetRasterBand(1).GetStatistics(True, True) min = round(stats[0]) max = round(stats[1]) colors = gdal.ColorTable() colors.CreateColorRamp(min, (94, 79, 162), max, (50, 136, 189)) vrt.GetRasterBand(1).SetRasterColorTable(colors) vrt.GetRasterBand(1).SetRasterColorInterpretation(gdal.GCI_PaletteIndex) gdal.Translate(geotiff, vrt)
Mais il me renvoie un message d'erreur:
Code:
ERROR 1: ./tmp/monfichier.tif: Unable to export color table to GeoTIFF file. Color tables can only be written to 1 band or 2 bands Byte or UInt16 GeoTIFF files.
Or mes données sont en Float32, ce sont des altitudes
J'ai bien essayé d'ajouter des options, mais ça ne change rien:
Code:
options_string = " ".join( [ "-b 1", "-ot Byte", "-expand rgb", ] ) gdal.Translate(geotiff, vrt, options=options_string)
Comment faire pour gérer des données en Float32?
Merci,
Nico
Hors ligne