Plot 填充区域基于 {osmdata} 使用 {ggplot2} [重复]

Question

这个问题在这里已经有了答案：

Sf package：关闭复杂形状的多边形（1个答案）

9 小时前关闭。

下面的表达式显示了我想如何通过 {osmdata} 和 {ggplot2} 创建一个 map，其中包含海洋/海洋。 我想用颜色填充陆地和/或海域。 然而，这样做似乎出乎意料地困难。 这篇博文甚至声称无法做到。

{osmplotr} 的这个小插图似乎必须解决： “因为 OpenStreetMap 将海岸线表示为线对象，所有海岸线数据都包含在$osm_lines object 中osm_line2poly() function 然后可以将这些线转换为多边形plot 填充区域。 "。 然而，就像在这个类似的 StackOverflow 问题中一样，function 会抛出一个错误，如在 reprex 的底部所示。 我还在这里发现 {tigris} package 可以提供必要的多边形数据 - 但仅适用于美国。

那么我怎样才能让它工作呢？

library(osmdata)
library(osmplotr)
library(sf)
library(tidyverse)

# define example bbox
bb <- tribble(
  ~xy,  ~min,  ~max,
  "x", 12.00, 12.18,
  "y", 54.08, 54.20
) %>% column_to_rownames("xy") %>% as.matrix()

# get "water"
water <- opq(bb) %>%
  add_osm_feature(key = "natural", value = "water") %>%
  osmdata_sf()

# get "coastline"
coast <- opq(bb) %>%
  add_osm_feature(key = "natural", value = "coastline") %>%
  osmdata_sf()

# ggplot
ggplot() +
  geom_sf(
    data = water$osm_multipolygons,
    fill = "navy",
    color = NA
  ) + 
  geom_sf(
    data = coast$osm_lines,
    fill = "navy",
    color = "blue"
  )

# trying osm_line2poly()
osmplotr::osm_line2poly(coast$osm_lines, bb)
#> Error in FUN(X[[i]], ...): unbenutztes Argument (V = c(3, 1, 6, 7, 2, NA, 5))

^{使用reprex v2.0.2创建于 2022-09-23}

Answer 1

感谢@JindraLacko ，我能够让我的 reprex 工作。 基本上，我们创建一个矩形/多边形，它是我们的 bbox 大小，然后通过海岸线将其分割。

library(lwgeom)
library(osmdata)
library(osmplotr)
library(sf)
library(tidyverse)

### define example bbox
lon_min <- 12.00 # xmin
lon_max <- 12.18 # xmax
lat_min <- 54.08 # ymin
lat_max <- 54.20 # ymax

bb <- get_bbox(c(lon_min, lat_min, lon_max, lat_max))

### get "water" that is not sea as polygons
water <- opq(bb) %>%
  add_osm_feature(key = "natural", value = "water") %>%
  osmdata_sf()

### get sea & land as polygons
# 1. get coastline (as line)
coast <- opq(bb) %>%
  add_osm_feature(key = "natural", value = "coastline") %>%
  osmdata_sf()

# 2. get overall rectangle for bbox
bb_rect <- data.frame(
  lat = c(lat_min, lat_max),
  lon = c(lon_min, lon_max)
) %>% 
  st_as_sf(coords = c("lon", "lat"), crs = 4326) %>% 
  st_bbox() %>% 
  st_as_sfc()

# 3. split overall rectangle for bbox via coastline
bb_rect_split <- bb_rect %>% 
  st_split(coast$osm_lines) %>% 
  st_collection_extract("POLYGON")

# 4. extract splitted rectangle parts
land <- bb_rect_split[1]
sea  <- bb_rect_split[2]

### ggplot
ggplot() +
  geom_sf(
    data = land,
    fill = "bisque",
    color = NA
  ) + 
  geom_sf(
    data = sea,
    fill = "navy",
    color = NA
  ) + 
  geom_sf(
    data = water$osm_multipolygons,
    fill = "navy",
    color = NA
  )

^{使用reprex v2.0.2创建于 2022-09-26}