【发布时间】:2015-09-28 20:00:26
【问题描述】:
我有一个光栅文件nightlights,它可以根据卫星图像绘制夜间光度图(300MB 压缩数据可用here)。它需要经纬度输入进行提取,例如extract(nightlights, data.frame(long, lat) 用于长纬度输入。
给定point <- (long, lat),我希望能够汇总所有位于point 5 英里范围内的值,并计算平均亮度,忽略 NA。
raster 包带有执行此功能的an aggregate method,但仅适用于矩形而不是球形。
【问题讨论】:
-
基本上你想用一个循环内核来处理它,并查询得到的图像。也许像this 这样的东西? (我假设你会做很多这样的查询 - 但也许不是这样的......)