1

もともと SpatialPolygons オブジェクトに含まれていた多数のポリゴンから座標を抽出しようとしています:

 Sr1 = Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))
 Sr2 = Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))
 Srs1 = Polygons(list(Sr1), "s1")
 Srs2 = Polygons(list(Sr2), "s2")
 SpP = SpatialPolygons(list(Srs1,Srs2), 1:2)

SpP オブジェクトから Sr1 と Sr2 の座標を抽出しようとしています。スタック交換の他の場所でこのコードを見ました:

Coords<-SpP@polygons[[2]]@Polygons[[1]]@coords

角括弧内のインデックスの不一致を理解できませんが、実行されます。しかし、出力は、Sr1 または Sr2 で指定した座標と一致しません。インデックスのすべての組み合わせを試しましたが、探している答えが得られません!

4

1 に答える 1

1

よろしいですか？それらは同じように見えます（コメントするには長すぎるため、「回答」として投稿するだけです）：

library(sp)

Sr1 <- Polygon(cbind(c(2, 4, 4, 1, 2), c(2, 3, 5, 4, 2)))
Sr2 <- Polygon(cbind(c(5, 4, 2, 5), c(2, 3, 2, 2)))
Srs1 <- Polygons(list(Sr1), "s1")
Srs2 <- Polygons(list(Sr2), "s2")
SpP <- SpatialPolygons(list(Srs1, Srs2), 1:2)

SpP@polygons[[1]]@Polygons[[1]]@coords

##      [,1] [,2]
## [1,]    2    2
## [2,]    1    4
## [3,]    4    5
## [4,]    4    3
## [5,]    2    2

Sr1@coords
##      [,1] [,2]
## [1,]    2    2
## [2,]    4    3
## [3,]    4    5
## [4,]    1    4
## [5,]    2    2


SpP@polygons[[2]]@Polygons[[1]]@coords
##      [,1] [,2]
## [1,]    5    2
## [2,]    2    2
## [3,]    4    3
## [4,]    5    2

Sr2@coords
##      [,1] [,2]
## [1,]    5    2
## [2,]    4    3
## [3,]    2    2
## [4,]    5    2
于 2015-06-08T02:14:36.877 に答える