Question

假设我有以下sf数据框：

library(sf)

nrows <- 10
geometry = st_sfc(lapply(1:nrows, function(x) st_geometrycollection()))
df <- st_sf(id = 1:nrows, geometry = geometry)

我还有以下列表：

mylist = list('2'=st_point(c(-73,42)), '3'=NA)

我想用第二个观察中的几何替换列表中的点。我曾想过要做以下事情：

st_geometry(df[names(mylist),]) <- st_sfc(mylist)

但这会引发错误：

“vapply中的错误（lst，class，rep（NA_character_，3））：值必须长度为3，但FUN（X [[2]]）结果为长度1“

我首先通过消除NA值找到了以下解决方法：

condition <- mylist[!is.na(mylist)]
st_geometry(df[names(condition),]) <- st_sfc(condition)

有更好的方法吗？我可以强制mylist中的NA元素为空点吗？

Answer 1

因为在。

中无法正常工作

mylist = list('2'=st_point(c(-73,42)), '3'= NA)

3不是一个POINT，而是一个逻辑，它（可能）不能被强迫＆＃34;以任何方式进入sf对象。

您可以通过将mylyst的NA元素替换为预先清空的点数来避免这种情况。例如：

mylist[[which(is.na(mylist))]]  <- st_point()
st_geometry(df[names(mylist),]) <- st_sfc(mylist)

，给予：

> df
Simple feature collection with 10 features and 1 field (with 10 geometries empty)
geometry type:  GEOMETRY
dimension:      XY
bbox:           xmin: -73 ymin: 42 xmax: -73 ymax: 42
epsg (SRID):    NA
proj4string:    NA
   id                 geometry
1   1 GEOMETRYCOLLECTION EMPTY
2   2           POINT (-73 42)
3   3              POINT EMPTY
4   4 GEOMETRYCOLLECTION EMPTY
5   5 GEOMETRYCOLLECTION EMPTY
6   6 GEOMETRYCOLLECTION EMPTY
7   7 GEOMETRYCOLLECTION EMPTY
8   8 GEOMETRYCOLLECTION EMPTY
9   9 GEOMETRYCOLLECTION EMPTY
10 10 GEOMETRYCOLLECTION EMPTY

HTH。

从sf中的列表替换几何

1 个答案: