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

Question

我正在尝试使用以下代码在 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)