R 모양 파일 단순화

Question

나는 아마존 강이있는 형태 파일을 가지고 있습니다. shapefile 단독으로 37.9 MB, 속성 테이블과 함께 42.1 MB까지 올라갑니다. 저는 각 1260x940 픽셀의 모든 브라질 아마존의 PNG 이미지를 생성하고 있으며, shapefile의 모든 데이터는 각지도의 드로잉 속도를 늦추기 때문에 단순화하고 싶습니다.

gSimplify 함수는 rgeos 패키지에서 각 폴리곤을 단순화하는 것처럼 보이며 더 작은 폴리곤을 제거하지 않는 것처럼 보입니다. 나는 그것을 0.1과 1000의 허용 오차로 시도했다. 그리고 항상 길이 (shp @ polygons)에 같은 값인 27633을 얻는다. 그리고 마지막 플롯은 그릴 때 거의 같은 시간이 걸린다. 최종 래스터가 1260x940 픽셀이 될 것이라고 말하는 함수가 필요하므로 모든 불필요한 점을 제거 할 수 있습니다. 거기에 할 수있는 기능이 있습니까?

미리 감사드립니다. 여기

Answer 1

꽤 포괄적 인 솔루션 :

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

을 다음

알아낼 강 R.

에서 Shape 파일 객체는 당신입니다 : http://www.r-bloggers.com/simplifying-polygon-shapefiles-in-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 패키지에서 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)) 
     } 
    }