2012-11-01 3 views
-1

각각 12 개의 하위 영역이있는 두 개의 다른 영역이있는 모양 파일이 있습니다.이 모양 파일에서이 24 개의 하위 영역의 개별 모양 파일을 원합니다. 패키지 maptools 및 rgeos를 사용하여 시도했지만 it.Any 로그 대단히 감사하겠습니다. 감사.큰 모양 파일에서 많은 모양 파일을 만듭니다.

샤리프

+1

이미 시도한 항목의 R 코드를 제공하십시오. 사용자가 붙어있는 지점까지 재생성 할 수있는 것이 좋습니다. '모르겠다 '라는 말은 너무 애매하다. –

답변

1

당신이 관심있는 열에서 고유 한 값을 기준으로 루프에서 데이터를 분할하고, 부분 집합의 데이터를 기록 할 수 있습니다. 나는 maptools의 leu에서 rgdal을 사용하고 있지만, shapefile을 읽고 쓰는 데 maptools 함수를 사용하도록 코드를 쉽게 변경할 수 있습니다.

require(sp) 
require(rgdal) 

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

# CREATE VECTOR OF UNIQUE SUBREGION VALUES 
y <- unique([email protected]$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) 
    } 
+0

친애하는 에반스 감사합니다. 좋은 제안입니다. –