我有一个光栅文件nightlights，用于映射卫星图像的夜间亮度（可用的300MB压缩数据here）。它需要用于提取的经度 - 纬度输入，例如， { - 1}}用于长期投入。
鉴于extract(nightlights, data.frame(long, lat)，我希望能够汇总位于point <- (long, lat) 5英里内的所有值，并计算平均亮度，忽略NAs。
光栅包附带an aggregate method执行此功能，但仅适用于矩形而非球形圆。
答案 0 :(得分：3)
我认为您可以使用缓冲区参数：
e <- extract(nightlights, cbind(long, lat), buffer=10)