r - 函数 "extract"如何处理不同的投影？

Question

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

Answer 1

您可以找到SpatialPolygons的函数extract的源代码 here 。代码以以下片段开头:

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 的投影更改为栅格的投影)，尽管它没有记录在帮助页面中。