计算 2 个经纬度之间的距离

Calculate distance between 2 lat longs

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

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

R 抛出错误 "Error in .pointsToMatrix(x) : Wrong length for a vector, should be 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)}

请尝试以下脚本代码:

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

有关更多详细信息,请遵循以下 link:

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

尝试

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