r - 具有分类数据的栅格 map 图例

Question

我想绘制一个包含 4 个不同值( 1 )的栅格，并带有描述类别的分类文本图例，例如 2但有彩盒:

我试过使用图例，例如:

legend( 1,-20,legend = c("land","ocean/lake", "rivers","water bodies"))

但我不知道如何将一个值与显示的颜色相关联。有没有办法检索用“绘图”显示的颜色并在图例中使用它？

Answer 1

光栅可视 包裹包括一个 Raster levelplot() 的方法，绘制分类变量并生成适当的图例:

library(raster)
library(rasterVis)

## Example data
r <- raster(ncol=4, nrow=2)
r[] <- sample(1:4, size=ncell(r), replace=TRUE)
r <- as.factor(r)

## Add a landcover column to the Raster Attribute Table
rat <- levels(r)[[1]]
rat[["landcover"]] <- c("land","ocean/lake", "rivers","water bodies")
levels(r) <- rat

## Plot
levelplot(r, col.regions=rev(terrain.colors(4)), xlab="", ylab="")