恢复默认的 R sf map plot 默认图例的最简单方法是什么？

Question

我一直在尝试在等值线图上显示图例，但发现图例的 R 文档有点难以理解。已经读了几个小时了。

我看到下面有一个图例：

library(rnaturalearth)
world <- rnaturalearth::ne_countries(scale='medium',returnclass = 'sf')
class(world)
plot(world[,"adm0_dif"])

但以下不是：

plot(world[,"adm0_dif"], col=sf.colors(n = nrow(world), alpha=0.5))

在plot for sf object 的 R 帮助时代中，它声明：

“指定 col 禁止绘制图例键。”

所以，我想知道如果这就是我想要的，那么恢复默认图例的最简单方法是什么。

Answer 1

我使用 ggplot 比使用 base R plot 多，所以我不知道如何用 base R 添加图例。此外，我不熟悉 rnaturalearth 包和数据，所以我所做的只是复制你的情节使用 ggplot 并添加一个图例。我不确定这是否是您想要的图例，但我所做的只是用 ggplot 制作完全相同的图并添加了图例。如果它不是您想要的，您可以自定义图例。

library(ggplot2)

ggplot(data = world) +
  geom_sf(aes(fill = 1:nrow(world))) +
  scale_fill_viridis_c(name = "", alpha = 0.5, option = "plasma", trans = "sqrt") +
  theme_classic()

对于多个变量的并排图，我们可以使用gridExtra 包。

library(gridExtra)

plot1 <-  ggplot(data = world) +
  geom_sf(aes(fill = 1:nrow(world))) +
  scale_fill_viridis_c(name = "", alpha = 1, option = "plasma", trans = "sqrt") +
  ggtitle("adm0_dif")+
  theme(plot.title = element_text(hjust = 0.5),
        panel.background = element_rect(fill = "white", colour = "white"))

plot2 <-  ggplot(data = world) +
    geom_sf(aes(fill = labelrank)) +
    scale_fill_viridis_c(name = "", alpha = 1, option = "plasma", trans = "sqrt") +
    ggtitle("labelrank") +
    theme(plot.title = element_text(hjust = 0.5),
          panel.background = element_rect(fill = "white", colour = "white"))

grid.arrange(plot1, plot2, ncol=2)

对于未按数字编码的变量，如featurecla，请改用scale_fill_manual

plot_fea <-  ggplot(data = world) +
  geom_sf(aes(fill = featurecla), show.legend = F) +
  scale_fill_manual(name = "", values = c("mediumaquamarine"))+
  ggtitle("featurecla") +
  theme(plot.title = element_text(hjust = 0.5),
        panel.background = element_rect(fill = "white", colour = "white"))

plot_fea

注意我只为scale_fill_manual 中的值提供了值（“mediumaquamarine”），因为它只包含一个唯一变量。传递给 value 的颜色数量必须与列中唯一变量的数量相匹配。