根据经纬度计算两个地点的球面距离

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

地理空间距离计算及优化(根据两个点经纬度计算距离)

1. # coding: utf-8
   
2. 
   
3. # 1.球面距离公式
   
4. # 2.经纬度数据处理
   
5. 
   
6. # In[153]:
   
7. 
   
8. 
   
9. import pandas as pd
   
10. table = pd.read_excel('./data/行政区划数据库_with 经纬度-省市区-邮编-区号-拼音-简称.xls')
    
11. table_sf = table[table['LevelType']==1][['Name','lng','Lat']]
    
12. sf_a=list();sf_b=list()
    
13. for i in table_sf['Name']:
    
14. for j in table_sf['Name']:
    
15. sf_a.append(i),sf_b.append(j)
    
16. table_a_b=pd.DataFrame()
    
17. table_a_b['sf_a']=sf_a
    
18. table_a_b['sf_b']=sf_b
    
19. table_a_b_lng_Lat = table_a_b.merge(table_sf,how='inner',left_on='sf_b',right_on='Name').merge(table_sf,how='inner',left_on='sf_a',right_on='Name')[['sf_a','sf_b','lng_y','Lat_y','lng_x','Lat_x']]
    
20. 
    
21. 
    
22. # 3.计算距离
    
23. 
    
24. # In[147]:
    
25. 
    
26. 
    
27. df = table_a_b_lng_Lat[['sf_a','sf_b','lng_y','Lat_y','lng_x','Lat_x']]
    
28. df['distance']=pd.Series(df[['lng_y','Lat_y','lng_x','Lat_x']].apply(lambda row: getDistance(row['Lat_y'],row['lng_y'],row['Lat_x'],row['lng_x']), axis=1))
    
29. df=df[df['distance']!=0]
    
30. df.columns=['起点省份','终点省份','起点经度','起点纬度','终点经度','终点纬度','距离']
    
31. df=df[['起点省份','起点经度','起点纬度','终点省份','终点经度','终点纬度','距离']]
    
32. 
    
33. 
    
34. # In[156]:
    
35. 
    
36. 
    
37. df.to_excel("./result/根据经纬度计算两个地点的距离.xlsx")

复制代码

结果示例



省份编码-经纬度-计算.zip (408.23 KB, 需要: RMB 5 元)

2018-12-11 18:26:49 上传

需要: RMB 5 元  [购买]

附件文件夹
当前路径下进入cmd输入：

1. python 省份编码-经纬度.py

复制代码

即可完成！！！

扫码加我 拉你入群

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

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

分享0 收藏4 回帖

请问爬虫可以实现所有上市公司与7个中心的球面距离计算嘛，已知经纬度，谢谢！

面包味的锅包肉 发表于 2019-10-14 16:42
请问爬虫可以实现所有上市公司与7个中心的球面距离计算嘛，已知经纬度，谢谢！

7个中心是指？