r - 将地理多边形拆分为特定区域的多个多边形

Question

我正在尝试将 R 上的一个多边形拆分为多个面积相等的多边形。

我有多边形的边界，我需要使用 R 将其拆分为 1 公里乘 1 公里的方框。我想知道是否可以使用 R。

例如:

x <- extent(c(40.97453103, 41.06321504, -92.47427103, -92.36617044))
plot(x)

这将创建一个具有给定边界的框。我正在尝试在 1km x 1km 的范围内创建多个框，然后使用 ggmap 将其合并到谷歌地图上。

Answer 1

您可以使用 sf 包中的 st_make_grid 函数，但我们不知道您的坐标引用系统和使用的单位。这是一个具有任意大小的网格的示例:

library(sf)
#> Linking to GEOS 3.5.1, GDAL 2.1.3, proj.4 4.9.2
x <- cbind(c(40.97453103, 41.06321504, 41.06321504, 40.97453103, 40.97453103),
           c(-92.47427103, -92.47427103, -92.36617044, -92.36617044, -92.47427103))

x <- st_sf(st_sfc(st_polygon(list(x))))
grid <- st_make_grid(x, cellsize = c(0.01,0.01))

par(mar = c(1,1,1,1))
plot(x)
plot(grid, add = T)

由 reprex package 创建于 2018-02-25 (v0.2.0).