2016-03-02 137 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)