本文主要是介绍【PythonGIS】基于高德Api实现批量地址查询经纬度,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
之前因为同事需要几千个小区的经纬度信息,所以就帮同事写了一段Python代码,通过调取高德地图的api实现地址查询经纬度这个功能。对于如何使用经纬度查询地址的方法,我之前分享过博文:【Python入门教程】获取图片可视化精准定位(逆地理编码),如果大家感兴趣可以自己去看下。这次主要是分享如何通过地址查询经纬度。
我这里就不演示如何注册高德Api了,之前那篇文章里有,自己去看!
一、导入需要的库
import requests
import xlrd
import xlwt
二、主函数
我这里就不过多的介绍了,看看官方的api就都知道了。这里需要注意的就是输入的文件格式,我这里已经写死程序了。只能输入.xls表格文件,然后表格中只有1列数据(地址),数据从第1行开始。最后就是输出的文件路径也必须是.xls,这是xlwt和xlrd两个库的弊端。
还有一件事,就是修改par = {'address': address, 'key': '*********'},改成自己高德的key,不懂的看之前那篇文章!!!
def CoordTransform(excel_path, out_path):""":param excel_path: 输入excel表格,里面只有1列全是地址:param out_path: 输出的表格名称:return: None"""wb = xlrd.open_workbook(excel_path)# wb.n sheets # sheets数量# wb.sheet_names() # sheets名称ws = wb.sheet_by_index(0) # 通过表单索引获取# ws = wb.sheet_by_name(' ') # 通过表名获取wb_result = xlwt.Workbook(encoding='utf-8')ws_result = wb_result.add_sheet("经纬度")ws_result.write(0, 0, '编号')ws_result.write(0, 1, '地址')ws_result.write(0, 2, '经度')ws_result.write(0, 3, '纬度')for n in range(1, ws.nrows+1):try:address = ws.cell_value(n-1, 0)par = {'address': address, 'key': '*******************'}base = 'http://restapi.amap.com/v3/geocode/geo'response = requests.get(base, par)response.raise_for_status()answer = response.json()gps = answer['geocodes'][0]['location'].split(",")ws_result.write(n, 0, n)ws_result.write(n, 1, address)ws_result.write(n, 2, str(gps[0]))ws_result.write(n, 3, str(gps[1]))except requests.HTTPError or KeyError:continuewb_result.save(out_path)
三、总结
本次分享的内容没什么难度,略懂爬虫知识会一点表格操作即可完成,所以不过多跟大家介绍了。有一点很遗憾,就是这个查询到的坐标系是高德坐标系,也就是火星坐标系(GCJ-02)。它和实际的WGS84会有一点偏移,但是问题不大。毕竟国家加密坐标系肯定有自己的考虑,有点点偏移对于我来说够用了。
这篇关于【PythonGIS】基于高德Api实现批量地址查询经纬度的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!