我有几个具有不同几何/轮廓的裁剪栅格.具体而言,空间产量图来自同一领域的几年,但程度各不相同 - 测量结果并不总是整个领域,但在某些年份只是其中的一部分.我想计算这些地图的平均值,并将它们组合成一个平均值栅格.然而,这确实意味着,不是每个像素都说5层/栅格有一个值.我可以接受这些缺失的值为NA,因此最终的平均值只能通过让我们说3个栅格来计算场的部分,其中地图不重叠.
我想用'extend {raster}'扩展栅格,用NA值填充非重叠部分:
y <- extend(y, shape, value=NA) #Shape是一个矩形形状,可以包含所有的yield map栅格
对于所有栅格,这都很好.但他们仍然没有相同的程度.即使我被调整的程度setExtent()或extent() <- extent()矩形shape文件或其他扩展光栅甚至达到一个程度,我仍然得到:
compareRaster(x)出错:不同的数字或列
..当我想堆叠它们并使用时calc(y, fun=mean,...).原始栅格范围与重新采样太不相同.但它们确实具有相同的分辨率和CRS.
有谁知道如何解决这个问题？