开发者

Problem writing shapefiles with writeOGR in R

I've been getting an error in a script that iterates over some text files in a directory, does some data manipulation and writes to a ESRI shapefile using writeOGR from the rgdal package in R. I keep getting an error after a certain number of files are written (a开发者_开发问答ppears to be about 43-44 individual files, or ~14 complete shapefiles, as each shapefile creates a number of files).

The error is usually something like

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

Additionally, after the script has crashed I get an error like

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

When I try to read or do anything with another file.

showConnections

does not show any open files at this point.

Here is some code that reproduces the problem:

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

coordinates(meuse) = ~x + y

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

I'm using R 2.12.2 on OSX. Error happens on both 64 bit and 32 bit versions.


Your code worked for me on Ubuntu 11.04 and R 2.13, might be worth also to try writeSpatialShape from the maptools package.

0

上一篇:

下一篇:

精彩评论

暂无评论...
验证码 换一张
取 消

最新问答

问答排行榜