我知道要找到两个经纬度点之间的距离,我需要使用haversine函数:
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
我有一个数据框,其中一列是纬度,另一列是经度。我想知道这些点离设定点有多远,-56.7213600,37.2175900。如何从数据帧中获取值并将其放入函数中?
数据帧示例:
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
我无法确认计算是否正确,但以下方法有效:
在这样小的数据帧上,下面的代码实际上速度较慢,但我将其应用于100000行df:
与4.3秒的apply函数相比,它的速度快了近250倍,值得注意的是
如果我们将以上所有内容压缩为一行:
我们现在观察到进一步的加速大约是341倍。
相关问题 更多 >
编程相关推荐