美文网首页
利用R根据经纬度计算地理距离

利用R根据经纬度计算地理距离

作者: 赵会成 | 来源:发表于2019-02-17 11:24 被阅读0次

计算两两地点之间的空间距离

经纬度表格示例

Pop Lon Lat

shenyang 123.3087 41.6467

yingkou 122.4622 40.7922

#!/usr/bin/env Rscript

# 根据经纬度计算距离

require(geosphere);

# 读取数据

lonlat = read.table('Book2.txt', header=TRUE);

# matrix

muer.lonlat = cbind(lonlat$Lon, lonlat$Lat);

# 精确计算,椭圆

muer.dists = distm(muer.lonlat, fun=distVincentyEllipsoid);

rownames(muer.dists) = lonlat$Pop;

colnames(muer.dists) = lonlat$Pop;

write.table(muer.dists, 'muer_geodists.txt', quote=F, sep="\t") # 写出去

相关文章

网友评论

      本文标题:利用R根据经纬度计算地理距离

      本文链接:https://www.haomeiwen.com/subject/ytcgeqtx.html