如何在光栅图中重现原始配色方案？

Question

我正在尝试使用此代码在R中绘制this raster file：

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)