so che per trovare la distanza tra due, i punti di longitudine latitudine ho bisogno di usare la funzione haversine:Vectorised Haversine formula con un dataframe panda
def haversine(lon1, lat1, lon2, lat2):
lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2])
dlon = lon2 - lon1
dlat = lat2 - lat1
a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
c = 2 * asin(sqrt(a))
km = 6367 * c
return km
Ho un dataframe dove una colonna è la latitudine e un'altra colonna è la longitudine. Voglio scoprire quanto sono lontani questi punti da un set point, -56.7213600, 37.2175900. Come prendo i valori dal DataFrame e li inserisco nella funzione?
esempio dataframe:
SEAZ LAT LON
1 296.40, 58.7312210, 28.3774110
2 274.72, 56.8148320, 31.2923240
3 192.25, 52.0649880, 35.8018640
4 34.34, 68.8188750, 67.1933670
5 271.05, 56.6699880, 31.6880620
6 131.88, 48.5546220, 49.7827730
7 350.71, 64.7742720, 31.3953780
8 214.44, 53.5192920, 33.8458560
9 1.46, 67.9433740, 38.4842520
10 273.55, 53.3437310, 4.4716664
Ho inviato una risposta che risolve il problema, ma penso che sia sub-ottimale – EdChum
Per un piccolo df la mia risposta sarebbe OK ma per un df molto più grande sarebbe meglio convertire i gradi in radianti e memorizzare questo e poi eseguire il calcolo sull'intero dataframe – EdChum
Ho 143 valori nel df e questo sembra aver funzionato bene. – user3755536