我发现readWKT
封装rgeos
可能是解决此问题的最佳方法。但是,我在此过程中遇到了几个问题。
我似乎无法解决的问题是我需要接收多个多面体,因此在将空间数据添加回我的主数据帧时the_geom_sp
,我可以在命令中将其作为向量调用。但是,在我有一个预先存在的多面体向量之前SpatialPolygonsDataFrame
,我无法制作the_geom_sp
成 a 。所以我将它格式化为一个列表 - 但是当我在命令中调用它时,它无法识别列表中的任何多面体对象。运行循环时,我还收到一条错误消息,表明我需要被归类为 a才能无损地接收来自循环的条目。我也受到 的约束,因为它一次只能应用于一个观察。SpatialPolygonsDataFrame
SpatialPolygonsDataFrame
" In `[<-`(`*tmp*`, n, value = the_geom_sp_entry) : implicit list embedding of S4 objects is deprecated
the_sp_geom
SpatialPolygonsDataFrame
readWKT
有谁知道我可以如何解决这个问题?我已经继续并为项目的这一部分切换到 Python,但希望你们能提供任何帮助。
the_geom_sp_entry <- SpatialPolygons(list())
the_geom_sp <- list()
for (n in 1:length(ld_2016$PROP_ID)) {
the_geom_sp[n] <- the_geom_sp_entry
the_geom_sp[n] <- readWKT(ld_2016$the_geom[n])
}
sp_2016 <- SpatialPolygonsDataFrame(the_geom_sp, ld_2016)