聚合栅格中距给定点(即圆)的给定距离内的值

Aggregate the values in a raster that lie within a given distance from a given point (i.e. a circle)

我有一个栅格文件 nightlights,它映射来自卫星图像的夜间光度(300MB 压缩数据可用 here)。它需要经纬度输入进行提取,例如extract(nightlights, data.frame(long, lat) 用于长纬度输入。

给定 point <- (long, lat),我希望能够汇总位于 point 5 英里范围内的所有值,并计算平均光度,忽略 NA。

栅格包附带 an aggregate method 执行此功能,但仅适用于矩形而不是球形圆。

我认为你可以为此使用缓冲区参数:

e <- extract(nightlights, cbind(long, lat), buffer=10)