具有两个lat列和两个long列的pandas数据帧中的距离计算

Question

我有一个带有这4列的pandas Dataframe df：

• pickup_latitude
• pickup_longitude
• dropoff_latitude
• dropoff_longitude

我想创建一个新的列，其中拾取点和下降点之间的距离。

我创建了这个函数：

def calcul_de_distance_vol_oiseau(data):

R = 6373.0

""" je change les variables en radians car python parle en radian et pas en degrées """
data['pickup_longitude'] = data['pickup_longitude'].apply(radians)
data['pickup_latitude'] = data['pickup_latitude'].apply(radians)
data['dropoff_longitude'] = data['dropoff_longitude'].apply(radians)
data['dropoff_latitude'] = data['dropoff_latitude'].apply(radians)

data['diff_lon'] = data['dropoff_longitude'] - data['pickup_longitude']
data['diff_lat'] = data['dropoff_latitude'] - data['dropoff_latitude']

data['calcul_intermediaire']= (data['diff_lat'] / 2).apply(sin)**2 + data['pickup_latitude'].apply(cos) * data['dropoff_latitude'].apply(cos) * (data['diff_lon'] / 2).apply(sin)**2

data['distance'] = R*np.arctan2((data['calcul_intermediaire']).apply(sqrt),(1 - data['calcul_intermediaire']).apply(sqrt))

return data

但是当我想验证我的点之间的距离时，我没有与本网站https://www.sunearthtools.com/fr/tools/distance.php中gps点之间的计算距离相同的结果。我认为我的功能有一个数学错误，但我找不到。

提前致谢

Answer 1

我认为你只需要：

data['pickup_longitude'] = data['pickup_longitude'].apply(radians)

和其他列的类似代码（使用lambda或定义函数）。

Answer 2

（迟到但是......）我建议你使用geopy.distance vincenty，定义一个函数来连接它与pandas并使用这个函数创建新的dateframe列：

from geopy.distance import vincenty

def pandasVincenty(row):
    '''calculate distance (m) between two lat&long points using the Vincenty formula '''

    return vincenty((row.pickup_latitude, row.pickup_longitude), (row.dropoff_latitude, row.dropoff_longitude)).meters 


data['distance'] =  data.apply(lambda r: pandasVincenty(r), axis=1)

这使用度数作为输入，并将米作为输出...