请求栅格网格点的值，知道它在 R 中的坐标

Question

我有一张在 R 中打开的栅格地图（数字高程模型），我需要知道一个特定网格点的值（即高度），知道其地理坐标。 我试过使用 extract() 函数，但它似乎不起作用，我试过将其转换为矩阵，但我只知道坐标，而不知道单元格的列号和行号。 这就是我的栅格数据的样子我认为这真的很简单，但我被困在这里。

知道如何得到它吗？

Answer 1

对于任何栅格r ，例如：

> r = raster(matrix(runif(100),10,10),xmn=3,xmx=13,ymn=10,ymx=20)

使用您的坐标创建一个包含两列的矩阵：

> m = matrix(c(12,16),ncol=2)
> m
     [,1] [,2]
[1,]   12   16

并提取：

> extract(r,m)
[1] 0.9597013

如果您有更多点，请向矩阵添加额外的行。

对于栅格之外的位置，您将返回 NA：

> m = matrix(c(12,21),ncol=2)
> extract(r,m)
[1] NA