r - R 程序星中曲线数据的空间运算

Question

我是 R 包 stars 的新手，我正在尝试对曲线数据进行基本的空间操作。我正在使用 netCDF 气候数据。我能够将 netcdf 连同我想用来指定我想要进行分析的区域的 shapefile 读入 r。我尝试使用 st_crop() 直接裁剪文件，但收到以下错误:

Warning message:
In st_crop.stars(test, wrst) : crop only crops regular grids

然后我尝试使用如下代码扭曲星星对象:

warp <- test %>% st_set_crs(3338) %>% st_warp(st_as_stars(st_bbox(), dx = 2))

但是我得到这个错误:

Error in colrow_from_xy(pts, x, NA_outside = TRUE) : 
  colrow_from_xy not supported for curvilinear objects

为了在给定区域进行分析，我是否需要“展平”我的曲线网格？如果是这样，我该怎么做？或者，相反，如果我能够执行类似 st_crop() 的操作或等效于 raster 操作的 calc() 或 stackApply () 使用曲线网格，有人能指出我正确的方向吗？非常感谢。

Answer 1

好吧，我想通了，这很简单。我能够通过以下简单代码使用 shapefile 对星星对象进行子集化:test[wrst]。无需变形或重新采样。