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

Spatial operations with curvilinear data in program stars for R

我是 Rstars 的新手,我正在尝试使用曲线数据进行基本空间操作。我正在使用 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

我是否需要 'flatten' 我的曲线网格才能在给定区域进行分析?如果是这样,我该怎么做?或者,相反,如果我能够使用曲线网格进行 st_crop()raster 等同于 calc()stackApply() 的操作,有人能指出我正确的方向吗方向?非常感谢。

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