在 R 中从 25 公里的海岸线创建缓冲区并计算近岸面积-解网

问：

赏金将在 5 天后到期。这个问题的答案有资格获得 +100 声望赏金。Alyssa C 正在寻找来自信誉良好的来源的答案：

我无法在 R 中找到解决方法。在 R 中欢迎 - 首选解决方案，但将在 GIS 中采用详细的解决方案，该解决方案可以作为 shapefile 上传到 R

我正在尝试创建一张包含多边形的地图，这些多边形表示从每个国家/地区的海岸线进入海洋的 25 公里缓冲区，以便我可以计算每个国家/地区在此缓冲区内的面积。我在使用上传的海岸线文件和 naturalearth 数据时都遇到了麻烦。这是我到目前为止尝试过的：

library(rgdal)
library(sf)
library(rnaturalearthdata)
library(rnaturalearth)
library(rgeos)
library(maps)
library(countrycode)

# World polygons from the maps package

world_shp <- sf::st_as_sf(maps::map("world", plot = T, fill = TRUE))

world_shp_df <- world_shp %>%
  mutate(Alpha.3.code = countrycode(ID, "country.name", "iso3c", warn = FALSE))

# Get coastlines using rnaturalearth
coastline <- ne_coastline(scale = "medium", returnclass = "sf")

# Buffer the coastlines by 25 km
buffer_distance <- 25000  # 25 km in meters
coastline_buffers <- st_buffer(coastline, dist = buffer_distance) %>%
  st_make_valid() 

ggplot() +
  geom_sf(data = coastline_buffers , fill = "lightblue")

这将导致一张带有水平线的地图，请参见下图：

我尝试过使用不同的 crs 简化几何形状，但我似乎无法弄清楚。任何帮助将不胜感激！

R 映射地理空间 rnaturalearth

# Create a band from latitudes -175° to 180° 
cleanup <- data.frame(lon=c(-175,-175,180,180,-175),lat=c(90,-90,-90,90,90)) %>% 
           st_as_sf(coords = c("lon", "lat"), 
           crs = st_crs(p[1,])) %>% 
           st_bbox() %>% 
           st_as_sfc()

# Remove points in this band
coastline_buffers_cleaned <- st_difference(coastline_buffers,cleanup)


ggplot() +
  geom_sf(data = coastline_buffers_cleaned , fill = "lightblue")

这种清理不是 100% 准确的，因为地图两端的一些沿海区域消失了，但条纹消失了，与总缓冲区相比，移除的区域可以忽略不计。

上一个：如何高效加载图像？

下一个：毛伊岛地图中图钉的延迟加载会在每次更新时闪烁/重新绘制所有图钉

在 R 中从 25 公里的海岸线创建缓冲区并计算近岸面积

Create buffer from coastline 25km in R and calculate nearshore area

评论