Come generare una griglia esagonale in R

Question

Vorrei poter creare un oggetto SpatialPolygons (che è una griglia esagonale) che copre un altro SpatialPolygon.

Vorrei che tutti gli esagoni avessero un diametro di 1 km (idealmente posso variare) e per tutti gli esagoni insieme per coprire l'intero oggetto. Il metodo seguito solo sembra coprire una piccola quantità di esso ...

Qui di seguito è il mio tentativo utilizzando il pacchetto 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) 

Qualsiasi aiuto come sempre è molto apprezzato ...

Answer 1

sostituire meuse.sr con una versione con buffer, come rgeos::gBuffer(meuse.sr, width = 2000) nella chiamata a spsample. Ecco un esempio completo che seleziona solo gli esagoni che si intersecano:

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)