Jak dodać tabelę kolorów do tiff jednopasmowego za pomocą GDAL?

11

Mam 1-pasmowy plik GeoTIFF. Interpretacja kolorów dla pasma jest szara. Chcę tabelę kolorów, aby przekształcić moje Col.tiffdo Col.jpeglub Col.png.

Próbowałem umieścić tabelę kolorów w tym tiffie przez .vrt:

gdal_translate -of VRT  Col.tiff Col.vrt

następnie dodając .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>

I przetłumacz ponownie:

gdal_translate Col.vrt Col.tiff

Ale odpowiedział:

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.

Jak mogę to zrobić?

Aby uzyskać informacje, gdalinfo dla zespołu:

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
źródło
Jakie wartości maksymalne / minimalne ma twój tif? Jest to Float64, więc w jaki sposób ta wartość dziesiętna powinna odpowiadać wartościom całkowitym dla palety kolorów?
AndreJ
Błąd mówi, że dane wyjściowe muszą być typu Byte lub UInt16. Być może powinieneś dodać „-ot UInt16” do komendy gdal_translate.
user30184
Dzięki, teraz działa, muszę tylko znaleźć dobre kolory :)
Faquarl,
@AndreJoost To, co zrobiłem, to pobranie maksymalnych / minimalnych wartości i podzielenie ich między 255 klas. Następnie dla każdej wartości umieszczam kolor.
Faquarl,
W końcu stworzyłem skrypt w języku Python, aby utworzyć nowy tiff z tymi samymi informacjami, z wyjątkiem tego, że istnieją 4 pasma (rgba)
Faquarl

Odpowiedzi:

5

za pomocą gdal możesz pokolorować obraz na podstawie gdal_dem ( color_relief )

składnia pliku konfiguracji kolorów pochodzi od tego obsługiwanego przez narzędzie GRASS r.colors . Pliki tabeli kolorów ESRI HDR (.clr) również są zgodne z tą składnią. Składnik alfa oraz obsługa tabulatorów i przecinków jako separatorów są rozszerzeniami specyficznymi dla 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
źródło
Zobacz także mój rozszerzony przykład tutaj: gis.stackexchange.com/questions/130199/…
AndreJ
link bezpośrednio do odpowiedzi @ AndreJ: gis.stackexchange.com/a/132690/49134
Tom Saleeba
tak, dziękuję, te pytania z 2015 r. należy oznaczyć jako duplikaty, ale odpowiedź AndreJ jest bardziej szczegółowa niż moja.
radouxju,