之前因为同事需要几千个小区的经纬度信息,所以就帮同事写了一段Python代码,通过调取高德地图的api实现地址查询经纬度这个功能。对于如何使用经纬度查询地址的方法,我之前分享过博文:【Python入门教程】获取图片可视化精准定位(逆地理编码),如果大家感兴趣可以自己去看下。这次主要是分享如何通过地址查询经纬度。
我这里就不演示如何注册高德Api了,之前那篇文章里有,自己去看!
一、导入需要的库
import requests
import xlrd
import xlwt
二、主函数
我这里就不过多的介绍了,看看官方的api就都知道了。这里需要注意的就是输入的文件格式,我这里已经写死程序了。只能输入.xls表格文件,然后表格中只有1列数据(地址),数据从第1行开始。最后就是输出的文件路径也必须是.xls,这是xlwt和xlrd两个库的弊端。
还有一件事,就是修改par = {'address': address, 'key': '*********'},改成自己高德的key,不懂的看之前那篇文章!!!文章来源:https://www.toymoban.com/news/detail-768492.html
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:
continue
wb_result.save(out_path)
三、总结
本次分享的内容没什么难度,略懂爬虫知识会一点表格操作即可完成,所以不过多跟大家介绍了。有一点很遗憾,就是这个查询到的坐标系是高德坐标系,也就是火星坐标系(GCJ-02)。它和实际的WGS84会有一点偏移,但是问题不大。毕竟国家加密坐标系肯定有自己的考虑,有点点偏移对于我来说够用了。文章来源地址https://www.toymoban.com/news/detail-768492.html
到了这里,关于【Python&GIS】基于高德Api实现批量地址查询经纬度的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!