关于r：具有分类数据的栅格地图的传说

Legend of a raster map with categorical data

我想绘制一个包含4个不同值(1)的栅格，并使用描述类别(例如2)但带有颜色框的分类文字图例：

我尝试使用诸如：

的图例

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

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

Initial

相关讨论

rasterVis软件包包括用于levelplot()的Raster方法，该方法绘制分类变量并生成适当的图例：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15

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="")

enter

默认情况下，栅格图中使用的颜色由rev(terrain.colors())生成(请参见?raster::plot)。您可以使用它来重新创建图例的4种颜色序列-或选择随机的颜色序列：

1 2	my_col = rev(terrain.colors(n = 4)) # my_col = c('beige','red','green','blue')

首先使用颜色顺序绘制地图。 legend = FALSE摆脱了标准颜色条：

1	plot(my_raster, legend = FALSE, col = my_col)

在左下方添加自定义图例。使用fill自变量生成彩色框：

1	legend(x='bottomleft', legend = c("land","ocean/lake","rivers","water bodies"), fill = my_col)