使用 terra 函数 rast 创建栅格时的 nrow 限制

Question

我正在尝试使用以下数据框（没有 NA）创建栅格，其中的纬度、经度和值是从 SeaDAS 生成的 HDF4 文件中提取的：

> dim(df)
[1] 10538622        3
> head(df)
          x         y         z
1 -61.38239 -30.91730 -0.015534
2 -61.38239 -30.91730 -0.015534
3 -61.38239 -30.91730 -0.015534
4 -61.38239 -30.91730 -0.015534
5 -61.37815 -30.91839 -0.015534
6 -61.37815 -30.91839 -0.015534

但我收到以下错误：

>   ras = rast(df, type = "xyz", crs = prj, digits = 6)
Error in matrix(NA, nrow = ncell(r), ncol = nlyr(r)) : 
  invalid 'nrow' value (too large or NA)
In addition: Warning message:
In matrix(NA, nrow = ncell(r), ncol = nlyr(r)) :
  NAs introduced by coercion to integer range

我已经尝试使用边界框减小数据框的大小，但即使只有 23530 个点，我仍然会遇到同样的错误。有几个点彼此非常接近（因此 df 中的相同点 1、2、3 和 4）并且网格可能不规则。这可能是什么问题？

谢谢，马修斯

Answer 1

如果值不在规则的网格上，您或许应该使用 rasterize。我猜测检测到的是规则网格，但空间分辨率极高。 terra 的开发版本现在应该在错误前显示警告。

使用 terra 函数 rast 创建栅格时的 nrow 限制

nrow limit when creating a raster using terra function rast

r

raster

terra