2

我需要使用函数 extract() 使用相等大小的正方形的网格单元从栅格中进行加权平均提取。我的多边形网格位于 UTM21n 中，栅格位于 GCS WGS84 基准 D 中。在将栅格用于 extract() 之前是否必须重新投影？或者该函数会正确处理它吗？

4

2 回答 2

5

你可以在这里extract找到函数的源代码。代码以以下代码段开头：SpatialPolygons 

setMethod('extract', signature(x='Raster', y='SpatialPolygons'), 
function(x, y, fun=NULL, na.rm=FALSE, weights=FALSE, cellnumbers=FALSE, small=FALSE, df=FALSE, layer, nl, factors=FALSE, sp=FALSE, ...){ 

    px <- projection(x, asText=FALSE)
    comp <- .compareCRS(px, projection(y), unknown=TRUE)
    if (!comp) {
        .requireRgdal()
        warning('Transforming SpatialPolygons to the CRS of the Raster')
        y <- spTransform(y, px)
    }
...

这表明它extract实际上执行了投影本身（将 SpatialPolygon 的投影更改为栅格的投影），尽管它没有记录在帮助页面中。

于 2013-10-11T07:32:12.667 回答
0

该文档没有提到自动重投影。所以，我认为假设该功能不这样做是节省的。因此，您需要在调用extract.

于 2013-10-11T06:49:19.957 回答