2016-03-02 153 views
1

我有4列在我的数据帧lat1,long1 ... lat2,long2。我需要计算这些对之间的距离。我正在尝试使用Distm功能。计算2 lat long之间的距离

当我尝试使用distm (c(mydata2$lst_upd_longitude,mydata2$lst_upd_latitude), c(mydata2$long,mydata2$lat), fun = distHaversine)

ř抛出了一个错误“在.pointsToMatrix错误(X):错误长度为一个矢量,应为2” 现在我使用以下代码来计算距离为每一点。但我相信应该有更好的解决方案。此代码也消耗大量时间。

for(i in 1:nrow(mydata2)){ 
    mydata2$distance[i] <- distm (c(mydata2$lst_upd_longitude[i],mydata2$lst_upd_latitude[i]), 
           c(mydata2$long[i],mydata2$lat[i]), 
           fun = distHaversine)} 

回答

1

尝试

df <- read.table(sep=",", col.names=c("lat1", "lon1", "lat2", "lon2"), text=" 
52,4,52,13 
39,116,52,13") 
library(geosphere) 
distHaversine(df[, 2:1], df[, 4:3])/1000 # Haversine distance in km 
+0

我Lat Long网数据已经在数据框中。当我直接应用distm函数得到错误.pointToMatrix(x)中的错误:向量的长度错误应为2 –

+0

我的纬度长数据也在数据帧中。当申请distHaverine我得到没有错误,但对之间的距离。 – lukeA

+0

工作真棒!谢谢 –

0

请尝试以下脚本代码:

function getDistanceFromLatLonInKm(lat1,lon1,lat2,lon2) { 
var R = 6371; // Radius of the earth in km 
var dLat = deg2rad(lat2-lat1); // deg2rad below 
var dLon = deg2rad(lon2-lon1); 
var a = 
Math.sin(dLat/2) * Math.sin(dLat/2) + 
Math.cos(deg2rad(lat1)) * Math.cos(deg2rad(lat2)) * 
Math.sin(dLon/2) * Math.sin(dLon/2) 
; 
var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a)); 
var d = R * c; // Distance in km 
return d; 
} 

function deg2rad(deg) { 
return deg * (Math.PI/180) 
} 

有关详细信息请按照下面的链接:

Calculate distance between two latitude-longitude points? (Haversine formula)