我试图了解如何从R中的栅格图层中提取数据。
我已下载tavg 10m data set，我想提取某个地点的平均温度。
到目前为止，我的方法是阅读1月份的图层并提取法兰克福的平均温度：
layer.jan <- raster("data/tmean_10m_bil/tmean1.bil")
layer.jan <- layer.jan/10
# get value for Frankfurt
loc <- as.data.frame(c(50.114928, 8.702152))
loc <- coordinates(loc)
temp <- extract(layer.jan, matrix(loc, ncol = 2))
然而，这会返回 23.5 的值，遗憾的是不是1月份法兰克福的平均气温。
我已尝试了解this示例中的输入值，但我并未真正了解这些值如何映射到经度/纬度坐标。
请假设我对GIS数据知之甚少，这是我第一次尝试探索这个区域。