我是 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 ()
使用曲线网格,有人能指出我正确的方向吗?非常感谢。
最佳答案
好吧,我想通了,这很简单。我能够通过以下简单代码使用 shapefile 对星星对象进行子集化:test[wrst]
。无需变形或重新采样。
关于r - R 程序星中曲线数据的空间运算,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/67894439/