2017-09-15 17 views
0

R 사용 - 특정 닫힌 컨투어 레벨 (폴리곤을 형성하기 위해) 내에서 둘러싸인 영역을 계산하고자합니다. 이것은 다른 곳에서 다루어졌지만이 닫힌 등고선이 서로 분리되어 있다면 영역을 따로 계산할 수도 있습니다.개별 닫힌 컨투어 내 R 영역

예 :이 영역을 제공

library(splancs) 
x <- c(1,0,0,1,1,1.5,1.5,3,3,1) 
y <- c(0,0,1,1,0,0,2,2,0,0) 
m <- cbind(x, y) 

plot(m, type="b") 
areapl(m) 

= 4 .... 그러나 I는 각각의 별도의 영역에 대응하는 부분 (C)의리스트 (1,3)로서 출력하고 싶은 고립 된 다각형.

이것은 장난감의 예일뿐입니다. 내 실제 데이터에서 나는 먼저 폴리곤의 좌표를 알지 못하고 거기에 고립 된 폴리곤이 몇 개인 지 모릅니다.

다른 아이디어가 유용 할 것입니다. 다른 R 공간 패키지를 사용하고 있습니다.

+0

데이터가 'sp' 패키지에서 Spatial *으로 쉽게 가져올 수있는 형식으로 제공되고 각 격리 된 폴리곤에 고유 한 ID가 있으면'rgeos :: gArea (m, byID = T)'가 당신이 필요로하는 것을하십시오. – Jul

답변

0

나는 래스터로 변환 래스터 패키지에서 덩어리를 사용하여이 문제를 해결 할 수 있었다 :

r <- raster(ncols=12, nrows=12) 
set.seed(0) 
r[] <- round(runif(ncell(r))*0.7) 
rc <- clump(r) 
freq(rc) 
plot(r) 
plot(rc) 

이 예제는 래스터 패키지에서 직접 - 그것은 다른 사람에게 관련 될 수있는대로 내가 그것을 게시!