2017-09-15 4 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 ....しかし、私は、それぞれの別々の領域に対応する領域Cのリスト(1,3)として出力したいと思います分離されたポリゴン。

これは単なるおもちゃの一例です。私の実際のデータでは、最初にポリゴンの座標を知ることはできず、ポリゴンがいくつあるか分かりません。

任意のアイデアは素晴らしいことだ - 私は他のR空間のパッケージを使用することに開いています。

+0

のデータが 'sp'パッケージから容易*空間にインポートすることができるフォーマットで来る場合、各単離されたポリゴンが '、rgeosを自身のIDを持っている場合:: 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) 

この例では、ラスタパッケージから直接である - それは他の人に関連するかもしれないと私はそれを掲載しています!