從大型文件中創建許多形狀文件。

Question

我有一個形狀文件，每個區域有兩個不同的區域，每個區域有12個子區域。我想要從這個形狀文件中分離出這24個子區域的獨立形狀文件。我也嘗試過使用包maptools和rgeos，但是無法計算出它。任何對數將非常感激。謝謝。

謝里夫

Answer 1

可以基於在感興趣的列中的唯一值的循環分割你的數據和寫出來的數據子集。我在maptools的leu中使用rgdal，但是您可以輕鬆更改代碼以使用maptools函數讀取/寫入shapefile。

require(sp) 
require(rgdal) 

# READ SHAPEFILE 
dat <- readOGR("C:/DATA", "dat") 

# CREATE VECTOR OF UNIQUE SUBREGION VALUES 
y <- unique(dat@data$SUBREGIONS) 

# CREATE SHAPEFILE FOR EACH SUBREGION AND WRITE OUT 
    for (i in 1:length(y)) { 
    temp <- dat[dat$SUBREGIONS == y[i], ] 
    writeOGR(temp, dsn=getwd(), y[i], driver="ESRI Shapefile", 
      overwrite_layer=TRUE) 
    }