在 ggplot 和 shapefile 中创建的重叠图

Question

我试图在使用 ggplot 创建的地图上重叠 shp 层，当我绘制两个独立的图时，你建议我如何重叠两个层？

data.shape<-readOGR(dsn="departamentos",layer="DEPARTAMENTOS")


ggplot()+
  geom_tile(data = tx_trend, aes(x = longitude, y = latitude, fill = slope))+
  scale_fill_gradientn(colors = rev(pals::linearlhot(100)), name = "ºC/10y", limits = c(0.1,0.5)) +
  #scale_fill_gradientn(colors = (pals::isol(100)), name = "ºC/10y", limits = c(0.1,0.45)) +
#  geom_point(data = filter(tx_trend, sign < 0.01),aes(x = longitude, y = latitude, color = "Sign. trend \n p-value <0.01"),
  geom_point(data = filter(tx_trend, sign < 0.01),aes(x = longitude, y = latitude, color = "Sign. trend \n p-value <0.01"),
                        size = 0.7, show.legend = T) +
  scale_color_manual(values = c("black"), name = "")+
  coord_fixed(1.3)+
  xlab("Longitude") + ylab("Latitude")+
  labs(title = "Decadal trend Summer", 
       subtitle = "(1981-2016)",
       caption = "")+
  theme_bw() +
  guides(fill = guide_colourbar(barwidth = 9, barheight = 0.5, title.position="right"))+
  theme(legend.position = "bottom")

Answer 1

如果没有可重复的示例很难，但我试一试！所以这是我的建议：

1. 使用install.packages("sf") 安装“sf”包并使用library(sf) 加载库

2. 使用st_read() 以sf 格式导入您的图层：

data.shape

1. 使用专用的geom_sf() 在您的图块和点图层上绘制您的形状。您只需要在 ggplot() 代码块中添加以下代码行：

geom_sf(data = st_geometry(data.shape), fill = NA, color = "red") +

coord_sf(default_crs = sf::st_crs(4326)) +

所以我建议：

ggplot()+
  geom_tile(data = tx_trend, aes(x = longitude, y = latitude, fill = slope))+
  scale_fill_gradientn(colors = rev(pals::linearlhot(100)), name = "ºC/10y", limits = c(0.1,0.5)) +
  #scale_fill_gradientn(colors = (pals::isol(100)), name = "ºC/10y", limits = c(0.1,0.45)) +
#  geom_point(data = filter(tx_trend, sign < 0.01),aes(x = longitude, y = latitude, color = "Sign. trend \n p-value <0.01"),
  geom_point(data = filter(tx_trend, sign < 0.01),aes(x = longitude, y = latitude, color = "Sign. trend \n p-value <0.01"),
                        size = 0.7, show.legend = T) +
  geom_sf(data = st_geometry(data.shape), fill = NA, color = "red") + # ADDED HERE
  coord_sf(default_crs = sf::st_crs(4326)) +       # ADDED HERE
  scale_color_manual(values = c("black"), name = "")+
  coord_fixed(1.3)+
  xlab("Longitude") + ylab("Latitude")+
  labs(title = "Decadal trend Summer", 
       subtitle = "(1981-2016)",
       caption = "")+
  theme_bw() +
  guides(fill = guide_colourbar(barwidth = 9, barheight = 0.5, title.position="right"))+
  theme(legend.position = "bottom")