# NOT RUN {
require(sp)
data(meuse)
coordinates(meuse) <- ~x+y
hex.polys <- hexagons(meuse, res=100)
plot(hex.polys)
plot(meuse,pch=20,add=TRUE)
# Points intersecting hexagons
hex.pts <- na.omit(over(meuse,hex.polys))
(hex.pts <- data.frame(PTID=rownames(hex.pts), hex.pts))
# }
Run the code above in your browser using DataLab