相当于rbindlist？

Question

我需要循环查找SpatialPolygonDataFrames的负载。 这些是通过缓冲道路文件创建的。 但是，由于我的SpatialLinesDataFrame类的“ roads”文件长为15000行，因此rbind函数变得越来越慢。 如果它们是数据帧，我将使用rblindlist加快速度。 但这不适用于sp对象（对吗？）。 有人有个好主意吗？

for (i in 1:nrow(roads)) {

  temp <- gDifference(gBuffer(roads[i,], byid = T, width =  15, capStyle = 'ROUND'), 
                  gBuffer(roads[i,], byid = T, width =  10, capStyle = 'ROUND'))

  slot(slot(temp, "polygons")[[1]], "ID") <- as.character(roads[i,]$oid)

  if (i == 1) {difference <- temp}

  if (i > 1) {difference <- rbind(difference, temp)}

  rm(temp)
  print(i)

}

谢谢

詹姆士

Answer 1

我不确定您是否需要循环播放。 您可以使用包sf 。 这是您想要达到的目标的一个示例：

library(mapview) # for the example lines data
library(sf)
trails$diff = st_geometry(st_buffer(trails, dist = 1500)) / st_geometry(st_buffer(trails, 1000))
mapview(trails$diff)

这是您想要的结果吗？ 如果没有，您可能要提供一个可复制的示例。