小编Bla*_*ake的帖子

我将raster带有 4 个类的 a转换为 a SpatialPolygonsDataFrame，其中表示类的特征并不总是连接的。我使用disaggregate(package sp) 来获取单个功能。

我现在想摆脱所有这些小于 3000 平方米的功能。我已经使用poly.areas（包 GISTools）生成了该区域，并将信息作为属性添加到 Shapefile 中。

小特征应与其周围的特征合并。

有谁知道如何解决这个问题？unionSpatialPolygons仅适用于 ID，不适用于属性条件。

我有几个具有不同几何/轮廓的裁剪栅格.具体而言,空间产量图来自同一领域的几年,但程度各不相同 - 测量结果并不总是整个领域,但在某些年份只是其中的一部分.我想计算这些地图的平均值,并将它们组合成一个平均值栅格.然而,这确实意味着,不是每个像素都说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.

有谁知道如何解决这个问题？