从列表中SF更换几何

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) 

但是，这将引发一个错误：

"Error in vapply(lst, class, rep(NA_character_, 3)) : values must be length 3, but FUN(X[[2]]) result is length 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不是一个点而是一个合乎逻辑的，这（可能）不以任何方式“被迫”转化为一个sf对象。

你可以通过用mylyst的NA元素来预先清空POINTS来绕过这个问题。例如：

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。