R簡化shapefile

Question

我有一個amazonic大河流的shapefile。單獨的shapefile具有37.9 MB，與屬性表一起高達42.1 MB。我生成所有巴西亞馬遜的PNG圖像，每個像素爲1260x940像素，shapefile中的所有這些數據只會減慢每張地圖的繪製速度，所以我想簡化它。

gSimplify函數，在rgeos包中，似乎只是簡化每個多邊形，而不是擺脫較小的。我嘗試了0.1和1000的容差，並且總是得到長度（shp @ polygons）相同的值：27633.最終的繪圖幾乎需要同一時間繪製。我需要一個函數，我告訴最終的光柵將是1260x940像素，因此它可以刪除每個不必要的點。有沒有這樣的功能？

在此先感謝。

Answer 1

相當全面的解決方案在這裏：http://www.r-bloggers.com/simplifying-polygon-shapefiles-in-r/

總之，你需要讓你的多邊形區域：

area <- lapply(rivers@polygons, function(x) sapply(x@Polygons, function(y) y@area)) 

江河是您在R.

shape文件對象

然後你找出大的多邊形，並保留它們：

sizeth <- 0.001 #size threshold of polygons to be deleted 
    mainPolys <- lapply(area, function(x) which(x > sizeth)) 

    rivers@data <- rivers@data[-c(1:2),] 
    rivers@polygons <- rivers@polygons[-c(1:2)] 
    rivers@plotOrder <- 1:length(rivers@polygons) 
    mainPolys <- mainPolys[-c(1:2)] 

    for(i in 1:length(mainPolys)){ if(length(mainPolys[[i]]) >= 1 && 
    mainPolys[[i]][1] >= 1){ 
     rivers@polygons[[i]]@Polygons <- rivers@polygons[[i]]@Polygons[mainPolys[[i]]] 
     rivers@polygons[[i]]@plotOrder <- 1:length(rivers@polygons[[i]]@Polygons) } } 

這可能還不夠好，你可能不希望刪除任何多邊形，在這種情況下，從shapefiles包DP（）函數就可以了：

res <- 0.01 #the argument passed to dp() which determines extent of simplification. Increase or decrease as required to simplify more/less  
for(i in 1:length(rivers@polygons)){ 
     for(j in 1:length(rivers@polygons[[i]]@Polygons)){ 
     temp <- as.data.frame(rivers@polygons[[i]]@Polygons[[j]]@coords) 
     names(temp) <- c("x", "y") 
     temp2 <- dp(temp, res) 
     rivers@polygons[[i]]@Polygons[[j]]@coords <- as.matrix(cbind(temp2$x, temp2$y)) 
     } 
    }