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

R 将地理多边形拆分为特定区域的多个多边形,r,geospatial,polygon,R,Geospatial,Polygon,我正在尝试将R上的一个多边形拆分为多个面积相等的多边形 我有一个多边形的边界,我需要用R把它分割成1公里乘1公里的盒子。我想知道用R是否可行 例如: x <- extent(c(40.97453103, 41.06321504, -92.47427103, -92.36617044)) plot(x) x您可以使用sf软件包中的st_make_grid功能,但我们不知道您的坐标参考系统和使用的单位。 以下是具有任意大小网格的示例: 库(sf) #>链接至GEOS 3.5.1、GDAL 2

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

我有一个多边形的边界,我需要用R把它分割成1公里乘1公里的盒子。我想知道用R是否可行

例如:

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

x您可以使用
sf
软件包中的
st_make_grid
功能,但我们不知道您的坐标参考系统和使用的单位。 以下是具有任意大小网格的示例:

库(sf)
#>链接至GEOS 3.5.1、GDAL 2.1.3、项目4.9.2

一个有趣的问题。您将遇到的一个问题是,您的扩展框不完全是1km的单位(我在扩展框中测试的坐标彼此之间的距离为9.075km)。那么,1km箱子的起点在哪里?也就是说,你将有9个完整的1km正方形,但随后是一个0.075km范围的小溢流盒(对于我测试的一侧)。这是一个非常好的观察结果。引用将位于框的边缘。在这种情况下,我只能忽略0.075公里。下面的答案是正确的。为什么不将范围精确地修剪到9km x 9km,然后使用下面的答案。您可以使用geosphere package destPoint确定修剪坐标。我的另一个想法是在修剪后使用光栅化(raster())并使用res()的分辨率设置。是的,这就是我的想法。这是一个很好的想法,但要注意这一点,因为正如您所说的,这些框的大小是任意的,而不是1kmx1kmt。我想我可以用这个让我的小项目工作!!