2

我有一个只包含 1 和 NA 的大型光栅文件 (5GB)。我想将其转换为具有 1 的区域的多面体，相邻单元格溶解为一个多边形。

我已使用将文件导入 R

r = raster::raster(my_filename)
r
class      : RasterLayer 
dimensions : 17452, 45000, 785340000  (nrow, ncol, ncell)
resolution : 0.008, 0.008  (x, y)
extent     : -180, 180, -55.9875, 83.6285  (xmin, xmax, ymin, ymax)
crs        : NA 
source     : C://...binary_X01_januarysnow.asc 
names      : binary_X01_januarysnow

我尝试了几种方法来创建多边形：

  • rasterToPolygons 来自具有溶解==TRUE 选项的栅格（R 崩溃）
  • isoband 包中的 isoband（R 崩溃），

当我在我的栅格覆盖 appr 的一个子集上尝试它们时，这两种方法都按预期工作。西班牙地区，所以我认为问题只是数据的大小而不是我的代码。

  • 然后我尝试使用 read_stars 读取我的栅格，并使用stars::st_as_sf(st, as_points = FALSE, merge = TRUE, connect8 = TRUE). 这返回了一个空多边形，可能是因为该文件被读取为星星代理对象，但我不确定，我在网上找不到任何有关该文件的信息。然后我通过使用强制将栅格读取为星星而不是星星代理，read_stars(my_filename, proxy=FALSE)并尝试使用上述st_as_sf命令，但收到消息“错误：无法分配大小为 2.9 Gb 的向量”

我知道在最坏的情况下，我可能会降低光栅分辨率并因此降低尺寸，并且能够创建我想要的多边形（但分辨率较低），但我想知道是否有人有其他建议我可以尝试？1 和 NA 都位于大的连续区域中，因此如果有帮助的话，在边缘具有高分辨率就足够了。

PS 这是我在 StackOverflow 上的第一个问题，所以如果我的问题没有清楚地描述，我深表歉意。我不知道如何提供大型数据集的可重现示例。

4

1 回答 1

2

您正在寻找的是as.polygons()来自terra包，raster包的继任者。terra比处理大数据集更好raster。

于 2020-10-07T15:37:28.657 回答