Question

我有一些xy坐标作为SpatialPoints（points）对象，并使用它们从RasterLayer（栅格）中提取这些位置的温度值：

extract = extract(raster, points)

然而，有几个点落在栅格图层之外（即没有在陆地上绘图），我想使用缓冲参数将每个点周围的半径扩展10000米 -

extract2 = extract(raster, points, method="simple",buffer=10000, cellnumbers=TRUE)

This produces a "list" object i.e.- 

head(extract2)

[[1]]
  cell  value 
591332    165 

[[2]]
  cell  value 
475809     NA 

[[3]]
  cell  value 
534127     NA 

[[4]]
  cell  value 
534127     NA 

[[5]]
  cell  value 
534127     NA 

[[6]]
  cell  value 
534127     NA

但我想创建一个数据框，其中我在点位置（NA或温度值）和单元格编号具有栅格值，因此我可以访问栅格图层中感兴趣的单元格的原始xy坐标。我怎么能这样做？

Answer 1

当您提出问题时，请提供可重复的示例。这很容易，因为您可以在“提取”或其他帮助页面中构建示例。

# example data
r <- raster(ncol=36, nrow=18, crs='+proj=utm +zone=14 +datum=WGS84')
r[] <- 1:ncell(r)
xy <- cbind(x=-50, y=seq(-80, 80, by=20))

# extract
e <- extract(r, xy, buffer=10)
ee <- t(data.frame(e))
rownames(ee) <- NULL

data.frame(xy, ee)

上述情况在很多情况下都适用，但不是在边缘或lon / lat数据周围，因为可能存在不同数量的单元格。在这种情况下，您可以这样做：

 e <- extract(r, xy, buffer=10)
 m <- max(sapply(e, length))
 x <- rep(NA, m)
 ee <- t(sapply(e, function(y) {x[1:length(y)] <- y; x}))
 data.frame(xy, ee)

使用缓冲函数从栅格中提取xy坐标

1 个答案: