2016-11-04 60 views
-1

我有一个200个县的形状文件。我应该如何抽样以从现有200的子集县?我曾尝试使用下面的R代码里面尝试:R空间多边形采样错误

library(maptools) 
TXcounties <- readShapePoly("C:/Users/Rvg296/Downloads/TXCountiesShapeFiles/TXCounties.shp") 

idx <- sample(1:250, 25, replace = FALSE) 
df.TXcounties <- as.data.frame(TXcounties) 
SpatialPolygonsDataFrame(idx, df.TXcounties). 

但这是抛出一个错误,如:在SpatialPolygonsDataFrame(IDX,df.TXcounties)

错误:试图从获得槽“多边形”没有时隙的基本类(“整数”)的对象

+0

请你能提供一个[最小的可重复的例子(http://stackoverflow.com/questions/5963269 /如何对做 - 一个伟大-R-重复性,例如)?而不是使用'maptools :: readShapePoly()'我建议'rgdal :: readOGR()',因为它读取投影。 – Phil

回答

0

的问题是,使用的是IDX,整数向量,作为用于SpatialPolygonsDataFrame()的第一个参数,但该功能需要一个空间的多边形对象作为它的第一个论点。在任何情况下,你应该能够做到整个事情很多更轻松地与这样的:

result <- TXcounties[idx,]