在R中用writeOGR编写shapefile的问题

Question

我在脚本中遇到错误，该脚本遍历目录中的某些文本文件，执行一些数据操作并使用来自R中rgdal包的writeOGR写入ESRI shapefile。我在一定时间后仍然收到错误写入的文件数量（似乎是大约43-44个单独的文件，或~14个完整的shapefile，因为每个shapefile创建了许多文件）。

错误通常类似于

GDAL Error 4: Unable to open ./test_8.shp or ./test_8.SHP

此外，在脚本崩溃后，我收到类似

的错误

 cannot open file 'file_test.R': Too many open files

当我尝试阅读或对其他文件执行任何操作时。

showConnections

此时未显示任何打开的文件。

以下是一些可以重现问题的代码：

library(gstat)
library(rgdal)
data(meuse)

coordinates(meuse) = ~x + y

for (i in 1:50){
    writeOGR(meuse,'.',paste('test',i,sep="_"),driver="ESRI Shapefile")
}

我在OSX上使用R 2.12.2。 64位和32位版本都会出错。

Answer 1

您的代码在Ubuntu 11.04和R 2.13上为我工作，也值得尝试writeSpatialShape包中的maptools。