聚合栅格中距给定点(即圆)的给定距离内的值
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)
我有一个栅格文件 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)