使用geopy包提取城市经纬度并生成空间距离矩阵

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

代码来源：zhugby的python使用geopy包提取城市经纬度及计算距离矩阵

我自己跑的时候发现有一点点报错，也有人有一样的错误，所以分享一下，感谢zhugby的分享

1. #没有安装的话需要安装一下
   
2. pip install geopy
   
3. # 导入pandas模块
   
4. import pandas as pd
   
5. # 读入城市数据
   
6. df=pd.read_excel(r"你的文件")
   
7. df.head()   #打印前五行
   
8. 
   
9. # 使用geopy模块获取城市的经纬度坐标
   
10. from geopy.geocoders import Nominatim
    
11. import time
    
12. gps = Nominatim(user_agent='myuseragent')
    
13. 
    
14. loc_list_x=[]
    
15. loc_list_y=[]
    
16. coord=[]
    
17. for i in list(df['城市']):
    
18.     location = gps.geocode(i)
    
19.     loc_list_x.append(location.longitude)
    
20.     loc_list_y.append(location.latitude)
    
21.     coord.append((location.latitude,location.longitude))
    
22.     time.sleep(2)  # 暂停2秒钟
    
23. df['经度']=loc_list_x
    
24. df['维度']=loc_list_y
    
25. df['坐标']=coord
    
26. df.head()  #打印前五行
    
27. 
    
28. # 计算城市之间的空间距离
    
29. from geopy.distance import distance
    
30. import numpy as np
    
31. arr=np.zeros([len(coord),len(coord)])
    
32. for i in range(len(coord)):
    
33.     for j in range(len(coord)):
    
34.         arr[i,j]=distance(coord[i],coord[j]).km
    
35. dis_matrix=pd.DataFrame(arr,index=list(df['城市']),columns=list(df['城市']))
    
36. print(dis_matrix.head())  #打印前五行
    
37. 
    
38. # 将结果写入文件
    
39. df.to_excel('空间距离坐标.xlsx')
    
40. dis_matrix.to_excel('空间距离矩阵.xlsx')

复制代码

准备一个你需要的所有城市的excel，利用geopy模块获取每个城市的经纬度坐标，然后计算城市之间的空间距离，并将结果写入两个Excel文件中。
之前说的报错就是在获取城市经纬度的时候报错HTTPSConnectionPool(host='nominatim.openstreetmap.org', port=443)，搜了一下，看到说虽然nominatim是免费的，但是每秒最多一个请求，所以在每次循环中暂停了2秒（我试了1秒还是报错，所以增加成2秒，虽然时间久一点）
也可以去看下zhugby的原分享，有图片更直观

___________更新：有问题各位，好像相同读音的会串，比如佛山的坐标会识别成釜山，还不知道咋弄

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：距离矩阵 经纬度 GEO connection longitude 城市经纬度 空间距离矩阵

我又发现亳州也不对，他识别成博尔塔拉蒙古自治州，我不懂他这个怎么匹配的了，而且附件我好像删不掉