我希望能够创建一个覆盖另一个 SpatialPolygon 的 SpatialPolygons 对象(这是一个六角形网格)。
我希望所有六边形的直径为 1 公里(理想情况下我可以改变这个)并且所有六边形一起覆盖整个物体。下面的方法似乎只涵盖了一小部分......
下面是我尝试使用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)
一如既往的任何帮助将不胜感激......
replace meuse.sr
有一些缓冲版本,比如rgeos::gBuffer(meuse.sr, width = 2000)
在通话中spsample
。这是一个仅选择相交六边形的完整示例:
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)
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)