我有五个经度和纬度，形成这样的形状。
df <- c(order=1:5,
lon=c(119.4,119.4,119.4,119.5,119.5),
lat=c(-5.192,-5.192,-5.187,-5.187,-5.191))
我怎样才能使用sf这样的包轻松地将它们转换为 sf 多边形数据框？
## Simple feature collection with 1 feature and 0 fields
## geometry type: POLYGON
## dimension: XY
## bbox: xmin: 119.4 ymin: -5.192 xmax: 119.5 ymax: -5.187
## epsg (SRID): 4326
## proj4string: +proj=longlat +datum=WGS84 +no_defs
## geometry
## 1 POLYGON ((119.4 ...
我看到这个问题出现在搜索结果中，所以我想我会提供一种更灵活的方法来sf从一系列lat和lon坐标中创建多边形。
st_as_sf有一个参数coords，它将把给定的点作为数据框中的坐标列，并将这些列转换为sf POINT几何。然后，因为sf与 配合得很好dplyr，我们可以st_combine将点转换为 aMULTIPOINT并st_cast转换为POLYGON。与使用 的“手动”构造相比st_polygon，它的优点是我们不必仔细考虑关闭环或传递给构造函数的嵌套列表的正确级别，并且如果我们有多个多边形一组坐标我们可以group_by用来一次创建所有的多边形。
NB 从技术上讲，您可以使用do_union=FALSEinside of来执行此操作summarise，但我认为这种语法更清晰一些，更类似于正常summarise的 .
df <- data.frame(
lon = c(119.4, 119.4, 119.4, 119.5, 119.5),
lat = c(-5.192, -5.192, -5.187, -5.187, -5.191)
)
library(tidyverse)
library(sf)
#> Linking to GEOS 3.6.1, GDAL 2.2.3, proj.4 4.9.3
polygon <- df %>%
st_as_sf(coords = c("lon", "lat"), crs = 4326) %>%
summarise(geometry = st_combine(geometry)) %>%
st_cast("POLYGON")
polygon
#> Simple feature collection with 1 feature and 0 fields
#> geometry type: POLYGON
#> dimension: XY
#> bbox: xmin: 119.4 ymin: -5.192 xmax: 119.5 ymax: -5.187
#> epsg (SRID): 4326
#> proj4string: +proj=longlat +datum=WGS84 +no_defs
#> geometry
#> 1 POLYGON ((119.4 -5.192, 119...
plot(polygon)
由reprex 包（v0.2.0）于 2018 年 10 月 5 日创建。
相当于@Yo B.的答案，但有sf
library(sf)
df <- data.frame(lon=c(119.4,119.4,119.4,119.5,119.5),
lat=c(-5.192,-5.192,-5.187,-5.187,-5.191))
# You need first to close your polygon
# (first and last points must be identical)
df <- rbind(df, df[1,])
poly <- st_sf(st_sfc(st_polygon(list(as.matrix(df)))), crs = 4326)
poly
## Simple feature collection with 1 feature and 0 fields
## geometry type: POLYGON
## dimension: XY
## bbox: xmin: 119.4 ymin: -5.192 xmax: 119.5 ymax: -5.187
## epsg (SRID): 4326
## proj4string: +proj=longlat +datum=WGS84 +no_defs
## st_sfc.st_polygon.list.as.matrix.df....
## 1 POLYGON ((119.4 -5.192, 119...
编辑以回答评论中的问题
请参阅主要 sf 小插图以获得清晰而详细的解释sf，sfc并将sfg对象总结为：
用于表示简单特征的三个类是：
- sf，具有特征属性和特征几何的表（data.frame），其中包含
- sfc，包含每个特征（记录）的几何图形的列表列，由
- sfg，单个简单特征的特征几何。
该st_sfc函数仅构建几何列（这是一个多边形列表 - 这里只有一个多边形）。中的“c”sfc代表“列”。该函数st_sf构建一个完整的sf对象（它也有一个data.frame类），它是一个带有几何列的数据框。在给定的示例中，多边形没有附加数据（无属性）。您可以通过构建 data.frame 来附加数据：
poly <- st_sf(data.frame(landuse = "Forest",
size = 23 ,
st_sfc(st_polygon(list(as.matrix(df))))),
crs = 4326)
poly
## ## Simple feature collection with 1 feature and 2 fields
## geometry type: POLYGON
## dimension: XYZ
## bbox: xmin: 1 ymin: 119.4 xmax: 5 ymax: 119.5
## epsg (SRID): 4326
## proj4string: +proj=longlat +datum=WGS84 +no_defs
## landuse size geometry
## 1 Forest 23 POLYGON Z ((1 119.4 -5.192,...
然后，您可以从空间对象中提取这些元素中的每一个并检查它们的类：
完整的 sf 对象：带有 sfc 几何列的 data.frame
class(poly)
## "sf" "data.frame"
提取为列表的第三列：sfc 对象
class(poly[[3]])
## "sfc_POLYGON" "sfc"
几何列的第一个元素：sfg 多边形对象
class(poly[[3]][[1]])
## "XY" "POLYGON" "sfg"
library(sfheaders)从 20191004 在 CRAN 上可以获取一个 data.frame 并将其转换为sf对象
library(sf)
library(sfheaders)
df <- data.frame(
lon = c(119.4, 119.4, 119.4, 119.5, 119.5),
lat = c(-5.192, -5.192, -5.187, -5.187, -5.191)
)
sfheaders::sf_polygon(
obj = df
)
## given only two columns of data are in df there's no need to specify lon & lat arguments
# Simple feature collection with 1 feature and 1 field
# geometry type: POLYGON
# dimension: XY
# bbox: xmin: 119.4 ymin: -5.192 xmax: 119.5 ymax: -5.187
# epsg (SRID): NA
# proj4string:
# id geometry
# 1 1 POLYGON ((119.4 -5.192, 119...
我不知道“sf”，但以防万一你的意思是“sp”，这里是 SPDF 的完整结构
df <- data.frame(lon=c(119.4,119.4,119.4,119.5,119.5),
lat=c(-5.192,-5.192,-5.187,-5.187,-5.191))
require(sp)
spdf <- SpatialPolygonsDataFrame(
SpatialPolygons(
Srl=list(
Polygons(srl=list(
Polygon(coords=df)
), ID=1)
)
),
data=data.frame(a=1)
)
plot(spdf)