Como adicionar uma tabela de cores a um tiff de uma banda usando GDAL?

11

Eu tenho o arquivo GeoTIFF de 1 banda. A interpretação das cores da banda é cinza. Quero uma tabela de cores para transformar meu Col.tiffem Col.jpegou Col.png.

Eu tentei colocar uma tabela de cores neste tiff através de um .vrt:

gdal_translate -of VRT  Col.tiff Col.vrt

adicionando no .vrt:

<ColorInterp>Palette</ColorInterp>
<ColorTable>
   <Entry c1="0" c2="0" c3="255" c4="240"/>
   <Entry c1="0" c2="0" c3="254" c4="240"/>
   <Entry c1="0" c2="0" c3="253" c4="240"/>
   <Entry c1="0" c2="0" c3="252" c4="240"/>
</ColorTable>

E traduza novamente:

gdal_translate Col.vrt Col.tiff

Mas respondeu:

Warning 1: 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.

Como eu posso fazer isso?

Para mais informações, há gdalinfo para a banda:

Band 1 Block=55x18 Type=Float64, ColorInterp=Gray
  Description = Something
  Metadata:
    GRIB_COMMENT=Something[p]
    GRIB_ELEMENT=TMP
    GRIB_FORECAST_SECONDS=0 sec
    GRIB_REF_TIME=1401926400 sec UTC
    GRIB_SHORT_NAME=SMTHG
    GRIB_UNIT=[p]
    GRIB_VALID_TIME=1401926400 sec UTC
Faquarl
fonte
Quais são os valores máx / min do seu tif? É Float64, então como esse valor decimal deve corresponder aos valores inteiros de uma paleta de cores?
precisa saber é
O erro diz que a saída deve ser do tipo Byte ou UInt16. Talvez você deva adicionar "-ot UInt16" ao comando gdal_translate.
user30184
Obrigado, ele funciona agora, eu só precisa encontrar as boas cores :)
Faquarl
@AndreJoost O que eu fiz foi pegar valores max / min e dividi-los entre 255 classes. Então, para cada valor, eu coloco uma cor.
Faquarl
Finalmente eu fiz um script Python para fazer uma nova tiff com as mesmas informações, exceto há 4bands (rgba)
Faquarl

Respostas:

5

com gdal, você pode colorir uma imagem com base em gdal_dem ( color_relief )

a sintaxe do arquivo de configuração de cores é derivada daquela suportada pelo utilitário GRASS r.colors . Os arquivos de tabela de cores ESRI HDR (.clr) também correspondem a essa sintaxe. O componente alfa e o suporte de tabulação e vírgula como separadores são extensões específicas da GDAL

aspect: aspect oriented grey colors
aspectcolr: aspect oriented rainbow colors
bcyr: blue through cyan through yellow to red
bgyr: blue through green through yellow to red
byg: blue through yellow to green
byr: blue through yellow to red
celsius: blue to red for degree Celsius temperature
corine: EU Corine land cover colors
curvature: for terrain curvatures (from v.surf.rst and r.slope.aspect)
differences: differences oriented colors
elevation: maps relative ranges of raster values to elevation color ramp
etopo2: colors for ETOPO2 worldwide bathymetry/topography
evi: enhanced vegetative index colors
gdd: accumulated growing degree days
grey: grey scale
grey.eq: histogram-equalized grey scale
grey.log: histogram logarithmic transformed grey scale
grey1.0: grey scale for raster values between 0.0-1.0
grey255: grey scale for raster values between 0-255
gyr: green through yellow to red
haxby: relative colors for bathymetry or topography
ndvi: Normalized Difference Vegetation Index colors
population: color table covering human population classification breaks
population_dens: color table covering human population density classification breaks
precipitation: precipitation color table (0..2000mm)
precipitation_monthly: precipitation color table (0..1000mm)
rainbow: rainbow color table
ramp: color ramp
random: random color table
rstcurv: terrain curvature (from r.resamp.rst)
rules: create new color table based on user-specified rules read from stdin
ryb: red through yellow to blue
ryg: red through yellow to green
sepia: yellowish-brown through to white
slope: r.slope.aspect-type slope colors for raster values 0-90
srtm: color palette for Shuttle Radar Topography Mission elevation
terrain: global elevation color table covering -11000 to +8850m
wave: color wave
radouxju
fonte
Veja também meu exemplo estendido aqui: gis.stackexchange.com/questions/130199/…
AndreJ 3/15/15
link direto para a resposta de @ AndreJ: gis.stackexchange.com/a/132690/49134
Tom Saleeba
Sim, obrigado. Essas perguntas de 2015 devem ser marcadas como duplicadas, mas a resposta de AndreJ é mais detalhada que a minha.
radouxju