我有一个带有气候数据的栅格图层(称为sms_av),我试图将其限制为仅限于美国大陆的区域。为此,我使用美国地图(称为r)创建了另一个rasterlayer。 我试图覆盖它们并使用maske但我收到一条错误消息。这是我用来屏蔽的代码:
> x<-overlay(r, sms_av, fun='mask')
Error in (function (x, fun, filename = "", recycle = TRUE, ...) :
cannot use this formula, probably because it is not vectorized
我也尝试过它,但情节不会改变。我用这段代码相交: intersect_climate&lt; - intersect(sms_av,r)
最后我使用了gridmask: gridmask&lt; - mask(sms_av,r,maskvalue = TRUE) 图(gridmask) compareRaster(x,mask)出错:不同的数字或列
我仍然无法获取有关气候数据的美国地图,并感谢您给我的任何建议。谢谢!!!