3

だから私は次のように見える複数の種の範囲を持っています (たとえば、青色に着色) これはアフリカを東から西に走っています:

種の範囲の例

パッケージで使用gAreaすることで総面積を取得できます。rgeos私が知りたいのは、このファイルを構成する個々のポリゴンの数です。つまり、合計範囲に対していくつの異なる領域があるか (これは島、または単に分離された集団である可能性があります)、およびそれらのポリゴンの範囲は何ですか。私は次のコードを使用しています:

#Load example shapefile
shp <- readShapeSpatial("species1.shp")

#How many polygon slots are there?
length(shp@polygons)

>2

#How many polygons are in each slot
length(shp@polygons[[1]]@Polygons
length(shp@polygons[[2]]@Polygons

特定の面積を取得するには:

shp@polygons[[1]]@Polygons[[1]]@area

これは正しいです?範囲の真ん中にある湖がそれ自体でポリゴンを構成するのではないかと心配していますか? おおよそ次のようなリストになりたいと思います。

           Species A    Species B
Polygon 1      12          11
Polygon 2      13          10
Polygon 2      14          NA

上記のコードが正しければ、ポリゴンの数とその個々の範囲のすべての種のリストをコンパイルしたい場合、ループに渡すのは非常に簡単です。

ありがとう

4

2 に答える 2