Как воспроизвести исходную цветовую схему на растровом графике?

Question

Я пытаюсь построить this raster file в R, используя этот код:

df = raster("/path.file.tif") 
plot(df, interpolate=T) 

Участок должен выглядеть следующим образом:

Но почему выход выглядит, как показано ниже:

Answer 1

Попробуйте plotRGB и brick, по умолчанию растровый только считывает один диапазон из многодиапазонного набора данных. С raster() вы, вероятно, получаете только первую полосу, которая обычно будет «красной» и, возможно, целым числом без знака (значения от 0 до 255). Затем plot(df) дает вам палитру по умолчанию terrain.color, растянутую по этим значениям.

Это проверяется с помощью файла:

require(raster) 

df = brick("C:/temp/ASTGTM2_N22E086_dem_shade_color.tif") 
plotRGB(df, interpolate = TRUE) 

Это меньший пример, воспроизводимые из ?plotRGB:

require(raster) 
b <- brick(system.file("external/rlogo.grd", package="raster")) 
plotRGB(b, interpolate = TRUE)