Jak wygenerować sześciokątną siatkę w R

Question

Chciałbym móc utworzyć obiekt SpatialPolygons (który jest sześciokątną siatką), który obejmuje inny SpatialPolygon.

Chciałbym, aby wszystkie sześciokąty miały średnicę 1 km (najlepiej mogę to zmienić) i dla wszystkich sześciokątów razem, aby pokryć cały obiekt. Poniższa metoda wydaje się obejmować tylko niewielką ilość tego ...

Poniżej moja próba przy użyciu pakietu 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) 

Każda pomoc jak zawsze mile widziana ...

Answer 1

zastąpić meuse.sr z niektórymi buforowanymi wersjami, takimi jak rgeos::gBuffer(meuse.sr, width = 2000) w wywołaniu spsample. Oto pełna przykładów, które wybiera tylko przecinające sześciokątów:

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

library(rgeos) 
meuse.large = gBuffer(meuse.sr, width = 2000) 
HexPts <-spsample(meuse.large, type="hexagonal", cellsize=1000) 
HexPols <- HexPoints2SpatialPolygons(HexPts) 
plot(HexPols[meuse.sr,], add=TRUE)