R упростить шейп-файл

Question

У меня есть шейп-файл с амазонкими большими реками. Только шейп-файл содержит 37,9 МБ, а таблица атрибутов - до 42,1 МБ. Я генерирую PNG-изображения всех бразильских Amazon, размером 1260x940 пикселей, и все эти данные в шейп-файле только замедляют рисование каждой карты, поэтому я хочу упростить ее.

Функция 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.

Затем выяснить большие полигоны и хранить их:

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 пакета будет делать трюк:

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)) 
     } 
    }