我最近写了一个小函数来做这件事。

library(tidyverse)
library(sf)
#> Linking to GEOS 3.10.2, GDAL 3.4.1, PROJ 7.2.1; sf_use_s2() is TRUE

# create data.frame
df <-
  tribble(
    ~ lat,     ~ long,
    27.81674, -82.05954,
    27.81499, -82.05964,
    27.81498, -82.05872,
    27.81484, -82.05872,
    27.81492, -82.05889
  )

# write function
fn_convert_to_polygon <- function(df){
  df %>% 
    as.matrix(ncol = 2) %>% 
    as.numeric() %>% 
    matrix(ncol = 2) %>% 
    st_linestring() %>% 
    st_cast('POLYGON') %>% 
    st_sfc()
}

# apply function
fn_convert_to_polygon(df) %>% 
  st_as_sf()
#> Simple feature collection with 1 feature and 0 fields
#> Geometry type: POLYGON
#> Dimension:     XY
#> Bounding box:  xmin: 27.81484 ymin: -82.05964 xmax: 27.81674 ymax: -82.05872
#> CRS:           NA
#>                                x
#> 1 POLYGON ((27.81674 -82.0595...

创建于2023-10-10使用reprex v2.0.2