別の 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)
いつものようにどんな助けも大歓迎です...