在 R 中为经纬度数据创建半径内的质心标记

问题描述

我有经纬度数据（220000 个观测值 [包括重复项]），我想创建两个新列，定义 5 公里和 10 公里半径内的点数。

我的问题是向量大小，这意味着算法不会针对这个大小的数据运行。

这是我创建第一列 radius_5km 的代码。

long <- runif(n=220000,min=-153,max=174)
lat <- runif(n=220000,min=-42,max=67)
df <- data.frame(long,lat)

library(geosphere)
dfgesphe <- cbind(df,Radius_5km=rowSums(distm(df[,c("long","lat")],fun = disthaversine) / 1000 <= 5)) # number of points within distance 5 km

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）