r - 如何使用 sf 从点按因子构造/绘制多边形的凸包？

Question

我有一个物种出现的数据集，我试图通过制作凸包将其转换为出现区域。我可以手动执行此操作（即一次一个物种），但我真的很希望能够通过物种名称自动处理它。

可以在此处找到精简的示例数据集：https ://pastebin.com/dWxEvyUB

这是我目前手动操作的方式：

library(tidyverse)
library(sf)
library(rgeos)
library(maps)
library(mapview)
library(mapdata)
library(ggplot2)


fd <- read_csv("occurrence.csv")

spA.dist <- fd %>%
  filter(species == "sp.A") %>%
  dplyr::select(lon,lat) %>%
  as.matrix() %>%
  coords2Polygons(ID="distribution") %>%
  gConvexHull() %>%
  gBuffer()

spB.dist <- fd %>%
  filter(species == "sp.B") %>%
  dplyr::select(lon,lat) %>%
  as.matrix() %>%
  coords2Polygons(ID="distribution") %>%
  gConvexHull() %>%
  gBuffer() 

wrld2 = st_as_sf(map('world2', plot=F, fill=T))
ggplot() + 
  geom_sf(data=wrld2, fill='gray20',color="lightgrey",size=0.07) +
  geom_polygon(aes(x=long,y=lat,group=group),color="red",data=spA.dist,fill=NA) +
  geom_polygon(aes(x=long,y=lat,group=group),color="blue",data=spB.dist,fill=NA) + 
  coord_sf(xlim=c(100,300), ylim=c(-60,60))

这显示了一张地图，其中包含基于他们观察的凸包的两个物种出现区域。我意识到我在这里混合了不同的空间库，所以如果可能的话，最好在 sf 中完成这一切。在我的真实数据中，我有两个以上的物种，我可以复制并粘贴我为每个物种获得的代码，但似乎应该可以简化这一点，因此多边形（以及随后的凸包）是按因子级别构建的自动地。更像这样的东西：

polys <- st_as_sf(fd) %>%
  group_by(species) %>%
  magically_make_polygons(lon,lat) %>%
  st_convex_hull() %>%
  st_buffer()

我一直在寻找几天以及挖掘大量的文档。很多这种空间的东西对我来说是不直观的，所以我希望我缺少很多基本的理解。这可以做到吗？

Answer 1

这是一个可能的解决方案，使用tidyverse（实际上只有dplyr）和sf-package（以及mapview一些快速查看的包）。

你非常接近你自己的解决方案（kudo's）。诀窍是summarise分组数据，然后创建船体..

library( tidyverse )
library( sf )

#create simple feature
df.sf <- df %>%
  st_as_sf( coords = c( "lon", "lat" ), crs = 4326 )
#what are we working with? 
# perform fast visual check using mapview-package
mapview::mapview( df.sf )

#group and summarise by species, and draw hulls
hulls <- df.sf %>%
  group_by( species ) %>%
  summarise( geometry = st_combine( geometry ) ) %>%
  st_convex_hull()

#result
mapview::mapview( list( df.sf, hulls ) )