r - 在现有 geom_sf 图层下方插入 geom_sf 图层

Question

我有一张印度的基本 map ，其中包含州和边界、一些标签以及存储为 gg 对象的许多其他规范。我想生成一些带有区域图层的 map ，这些 map 将包含来自不同变量的数据。

为了防止地区 map 覆盖州和国家边界，它必须在所有之前的代码之前，我想避免重复。

我想我可以按照 this answer 为 gg 对象调用 $layers 来做到这一点.但是，它会引发错误。代表如下:

library(ggplot2)
library(sf)
library(raster)

# Download district and state data (should be less than 10 Mb in total)

distSF <- st_as_sf(getData("GADM",country="IND",level=2))

stateSF <- st_as_sf(getData("GADM",country="IND",level=1))

# Add border

countryborder <- st_union(stateSF)

# Basic plot

basicIndia <- ggplot() +
  geom_sf(data = stateSF, color = "white", fill = NA) +
  geom_sf(data = countryborder, color = "blue", fill = NA) +
  theme_dark()

basicIndia

# Data-bearing plot

districts <- ggplot() +
  geom_sf(data = distSF, fill = "gold")

basicIndia$layers <- c(geom_sf(data = distSF, fill = "gold"), basicIndia$layers)

basicIndia
#> Error in y$layer_data(plot$data): attempt to apply non-function

预期结果

任何帮助将不胜感激!

Answer 1

我仍然不确定我是否遗漏了您要查找的内容的详细信息，但 ggplot2 会按照您提供的顺序绘制图层。所以像

ggplot(data) +
  geom_col() +
  geom_point(...) +
  geom_line(...)

将绘制列，然后在这些列的顶部绘制点，然后在前一层的顶部绘制线。

sf 地 block 也是如此，这样可以很容易地制作像这样的多个地理级别的地 block 。

(我在 sf 对象上使用 rmapshaper::ms_simplify 只是为了简化它们并加快绘图速度。)

library(dplyr)
library(ggplot2)
library(sf)
library(raster)

distSF <- st_as_sf(getData("GADM",country="IND",level=2)) %>% rmapshaper::ms_simplify()
...

然后您可以按照需要显示的顺序添加图层来进行绘图。请记住，如果您需要对这些 sf 中的任何一个进行其他计算，您可以提前或在您的 geom_sf 中执行此操作。

ggplot() +
  geom_sf(data = distSF, fill = "gold", size = 0.1) +
  geom_sf(data = stateSF, color = "white", fill = NA) +
  geom_sf(data = countryborder, color = "blue", fill = NA)

关于尝试将一个绘图添加到另一个绘图:ggplot2 在图层中工作，因此您创建一个基础 ggplot 对象，然后在其上添加几何图形。因此，例如，您可以制作两个有效的图:

state_plot <- ggplot(stateSF) + 
  geom_sf(color = "white", fill = NA)
country_plot <- ggplot(countryborder) + 
  geom_sf(color = "blue", fill = NA)

但您不能添加它们，因为您将有 2 个基本 ggplot 对象。这应该是你提到的错误:

state_plot + 
  country_plot
#> Error: Don't know how to add country_plot to a plot

相反，如果您需要制作绘图，则在其上添加其他内容，制作基础 ggplot，然后添加几何图层，例如 geom_sf使用不同的数据集。

state_plot +
  geom_sf(data = countryborder, fill = NA, color = "blue")

^{由 reprex package 创建于 2018 年 10 月 29 日(v0.2.1)}