SpatialPointsDataFrame在SpatialLinesDataFrame上

Question

我有几个数据集(~25和所有相同的列)，总共有400,000次观测。我正在运行一个循环来提取每个数据测量点的道路类。当我在一个小数据集上运行我的脚本时(大约4,000次)，一切都很好。但是，当我使用相同的脚本运行一个循环时，在前47,000个obs之后有一个错误(停止循环)。我搜索了那个错误，没有发现任何结果。有什么建议吗？

谢谢马丁

.Call中的错误(“rgeos_intersects_prepared”、.RGEOS_HANDLE、spgeom1、spgeom2 )：不允许负长度向量

下面是我用于分析的代码的一部分

proj4string(Road) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
Road_pc <- spTransform(Road, CRS = CRS("+init=epsg:2957"))

output <- NULL 

for( j in 1:length(import.list)){

dfb <- as.data.frame(import.list[j])

coordinates(dfb) <- ~ coords.x1 + coords.x2 
proj4string(dfb) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
pc <- spTransform(dfb, CRS = CRS( "+init=epsg:2957" ))
pc100m <- gBuffer(pc, width=100, byid=TRUE ) 

zzc <- over(pc100m, Road_pc)
dfc <- cbind(dfb, zzc)
dfc$DAdates <- j
output<- rbind(dfc, output)
}

Answer 1

我不确定原因，但是当我在每个脚本上单独运行脚本时。它起作用了。它需要更多的时间，但它是有效的。m。