计算2个纬度之间的距离

Question

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

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

时

R抛出错误＆＃34; .pointsToMatrix（x）出错：向量的长度错误，应为2＆＃34; 现在我使用下面的代码来计算每个点的距离。但我相信应该有更好的解决方案。此代码也耗费了大量时间。

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)}

Answer 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

Answer 2

请尝试以下脚本代码：

 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)
}

有关详细信息，请访问以下链接：