如何从R中的栅格中提取单像素值

Question

我有一个栅格和气象站的位置。所以我想从我的栅格中提取对应于气象站位置的单像素值。我尝试了不同的方法，但当我与ArcMap交叉检查时，我无法获得确切的像素值。有没有人能帮我找出准确的像素值？谢谢

Answer 1

在提出R问题时，请包含一些数据和代码，以显示您尝试过的内容和正在使用的软件包。

示例数据

library(raster)
f <- system.file("external/test.grd", package="raster")
r <- raster(f)

如果要获取特定位置的值，可以执行以下操作

station <- cbind(179735, 331230)
extract(r, station)
# 259.9123 

或者首先计算单元格编号，并将其用作索引

i <- cellFromXY(r, station)
i
#[1] 5554
r[i]
# 259.9123 

至于你的后续问题(你应该问一个新的问题)“我如何从那个位置提取周围的9个像素的平均值？”

获取相邻单元格并提取

j <- adjacent(r, i, 8, pairs=F, include=T)
j
#[1] 5554 5473 5553 5633 5475 5555 5635 5474 5634
r[j]
#[1] 259.9123 273.5190 267.6167 218.8983 273.3555 252.1958 217.1910 274.5584
#[9] 223.2197

如果你有很多地方

stations <- rbind(cbind(179735, 331230), cbind(179061, 330212))
cells <- cellFromXY(r, stations)
# get adjacent values but set "pairs" to TRUE
adj <- adjacent(r, cells, 8, pairs=TRUE, include=TRUE, sorted=TRUE)
vadj <- r[ adj[,2] ]

现在使用tapply或aggregate来获得平均值

v <- tapply(vadj, adj[,1], mean)
v
#    5554     7537 
#251.1630 412.0748