r - R 程序星中曲线数据的空间操作

Question

我是包stars的新手，R并且正在尝试使用曲线数据进行基本的空间操作。我正在使用 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()使用曲线网格，有人可以指出我正确的方向吗？非常感谢。

标签： rsfr-rasterr-stars

Answer 1

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