Как добавить таблицу цветов в TIF с одной полосой, используя GDAL?


11

У меня есть 1-полосный файл GeoTIFF. Интерпретация цвета для полосы - серый. Я хочу таблицу цветов, чтобы преобразовать мои Col.tiffв Col.jpegили Col.png.

Я попытался поместить таблицу цветов в этот tiff через .vrt:

gdal_translate -of VRT  Col.tiff Col.vrt

затем добавив в .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>

И переведите еще раз:

gdal_translate Col.vrt Col.tiff

Но он ответил:

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.

Как мне это сделать?

Для информации есть gdalinfo для группы:

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

Какие максимальные / минимальные значения имеет ваш тиф? Это Float64, так как это десятичное значение должно соответствовать целочисленным значениям для цветовой палитры?
AndreJ

Ошибка говорит о том, что выходные данные должны иметь тип Byte или UInt16. Возможно, вам следует добавить «-ot UInt16» в команду gdal_translate.
user30184

Спасибо, это работает сейчас, мне просто нужно найти хорошие цвета :)
Faquarl

@AndreJoost То, что я сделал, взял значения max / min и разделил их на 255 классов. Затем для каждого значения я ставлю цвет.
Faquarl

Наконец, я написал скрипт на Python, чтобы создать новый tiff с той же информацией, за исключением 4-полосных (rgba)
Faquarl

Ответы:


5

с помощью gdal вы можете раскрасить изображение на основе gdal_dem ( color_relief )

синтаксис файла конфигурации цвета получен из синтаксиса, поддерживаемого утилитой GRASS r.colors . Файлы таблицы цветов ESRI HDR (.clr) также соответствуют этому синтаксису. Альфа-компонент и поддержка табуляции и запятой в качестве разделителей являются специфическими расширениями 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

Смотрите также мой расширенный пример здесь: gis.stackexchange.com/questions/130199/…
AndreJ

прямая ссылка на ответ @ AndreJ: gis.stackexchange.com/a/132690/49134
Том Салеба,

да, спасибо, эти вопросы 2015 года должны быть помечены как дубликаты, но ответ AndreJ более подробный, чем мой.
Радуксу
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.