15

別の SpatialPolygon をカバーする SpatialPolygons オブジェクト (Hexagonal グリッド) を作成できるようにしたいと考えています。

すべての六角形の直径を 1 km にし (理想的にはこれを変更できます)、すべての六角形を合わせてオブジェクト全体をカバーしたいと考えています。以下の方法は、そのほんの一部をカバーしているようです...

sp以下は、パッケージを使用した私の試みです:

require(sp)
data(meuse.riv)
meuse.sr = SpatialPolygons(list(Polygons(list(Polygon(meuse.riv)), "x")))
plot(meuse.sr)

HexPts <-spsample(meuse.sr,type="hexagonal",cellsize=1000)
HexPols <- HexPoints2SpatialPolygons(HexPts)
plot(HexPols, add=TRUE)

いつものようにどんな助けも大歓迎です...

4

1 に答える 1