我从这个 DataFrame 开始(因为我不想手动输入那么多,而且你没有为数据提供任何代码):
df
Index Number la ID
0 0 1 [43.3948, -23.9483] 1/1/90
1 1 2 [22.8483, -34.3948] 1/1/90
2 2 3 [44.9584, -14.4938] 1/1/90
3 3 4 [22.39458, -55.34924] 1/1/90
4 4 5 [33.9383, -23.4938] 1/1/90
5 5 6 [22.849, -34.397] 1/1/90
现在我介绍了一个人工列,它只是为了帮助我们获得距离的笛卡尔积
df['join'] = 1
df_c = pd.merge(df, df[['la', 'join','Index']], on='join')
下一步是通过.apply 应用vincenty 函数并将结果存储在额外的列中
df_c['distance'] = df_c.apply(lambda x: distance.vincenty(x.la_x, x.la_y).miles, 1)
现在我们有了原始矩阵的笛卡尔积,这意味着我们也有了每个城市与其自身的比较。但是我们将在下一步中通过执行-1 来考虑这一点。我们将Index_x 分组,并将所有小于 20 英里的距离相加。
df['num_close_cities'] = df_c.groupby('Index_x').apply(lambda x: sum((x.distance < 20))) -1
df.drop('join', 1)
Index Number la ID num_close_cities
0 0 1 [43.3948, -23.9483] 1/1/90 0
1 1 2 [22.8483, -34.3948] 1/1/90 1
2 2 3 [44.9584, -14.4938] 1/1/90 0
3 3 4 [22.39458, -55.34924] 1/1/90 0
4 4 5 [33.9383, -23.4938] 1/1/90 0
5 5 6 [22.849, -34.397] 1/1/90 1