赞
踩
我们会遇到一些需要批量把地址数据转换成坐标的问题,一个个查的话通常会比较麻烦,那有没有比较方便的方法呢,我们可以利用高德的地理编码/逆地理编码 API功能, 通过 HTTP/HTTPS 协议访问远程服务的接口,提供结构化地址与经纬度之间的相互转化的能力;
这里我们先拿一些数据作为分析的基础,数据源来源:''上海市公共数据开放平台'',数据更新日期:2022年9月29日,网址:加油站-上海市经济和信息化委员会-上海市公共数据开放平台 (sh.gov.cn)
我们可以看到原始数据包含了油站名称、地址、电话等一些标签,但是没有坐标,我们要做的是把地址批量转化坐标。
直接利用批量地理坐标转经纬度工具:地理编码 (高德) 地址转经纬度 - 批量工具网 (piliang.tech)
直接把我们获取的位置放进去,城市选择数据所在的城市即可,这里是''上海'' ,这里需要申请一下高德是key,个人开发者地理编码api调用的限额是5000条/天,数量量大于这个的话,可以多天处理或者多个账号处理,这里放上开发文档,如何申请key看官方文档:创建应用和 Key-Web服务 API | 高德地图API (amap.com)
这里可以直接下载转换结果,坐标系GCJ-02,当然还有个别地址描述太模糊的,会查不出坐标,手动查一下坐标即可,大部分还是可以查到的,以上海加油站分布数据为例,830条数据仅有8条查不出来地址,可以手动调整一下地址,或者直接用高德坐标拾取器,拾取坐标:坐标拾取器 | 高德地图API (amap.com);
这里我们拿高德的坐标拾取器挑一个点验证一下,准确度还可以;
我们把它放到GIS里看一下,导入GIS里时要注意一下文件格式是.csv或者.xls,打开arcgis10.8——>添加数据——>右键显示xy数据——>选择x,y字段,arcgisPro的话是支持.xlsx的,流程一样;
结果如下图所示;
当然这里通过python脚本调用高德地图地理编码API也可以,只要把文件路径和 API 密钥改成自己的即可;
完整代码#运行环境 Python 3.11
- import pandas as pd
- import requests
-
- def geocode(address, key):
- """
- 使用高德地图API进行地理编码
- :param address: 地址字符串
- :param key: 高德地图API密钥
- :return: 返回地理编码结果的字典
- """
- base = 'https://restapi.amap.com/v3/geocode/geo'
- params = {
- 'address': address,
- 'key': key,
- 'output': 'json'
- }
- response = requests.get(base, params=params)
- answer = response.json()
- return answer
-
-
- def process_addresses_from_xlsx(file_path, key):
- """
- 从XLSX文件中读取地址,并进行地理编码。
- :param file_path: XLSX文件路径
- :param key: 高德地图API密钥
- :return: None
- """
- # 读取XLSX文件
- df = pd.read_excel(file_path)
-
- # 假设地址在名为'地址'的列中
- addresses = df['地址']
-
- # 创建一个空列表以保存地理编码的结果
- results = []
-
- # 遍历地址列中的每一项
- for address in addresses:
- # 调用地理编码API
- result = geocode(address, key)
-
- # 尝试获取结果
- try:
- location = result['geocodes'][0]['location']
- lng, lat = location.split(',')
- # 添加经度和纬度到结果列表
- results.append({
- 'Address': address,
- 'Longitude': float(lng),
- 'Latitude': float(lat)
- })
- except (KeyError, IndexError):
- # 如果地理编码失败,则跳过该地址
- continue
-
- # 将结果保存到一个新的DataFrame
- results_df = pd.DataFrame(results)
-
- # 保存结果到新的XLSX文件
- results_df.to_excel('加油站GCJ-02坐标.xlsx', index=False)
-
-
- if __name__ == '__main__':
- # 替换这里的key为您自己的API密钥
- api_key = 'your_api_key_here' # 替换为您的实际API密钥
- file_path = 'E:/data/加油站.xlsx' # 根据实际情况进行更改文件路径
-
- # 处理XLSX文件并进行地理编码
- process_addresses_from_xlsx(file_path, api_key)
因为当前坐标系是GCJ02,需要批量转成WGS84/BD09的话可以用免费这个网站:地图坐标系转换 - 在线工具 (try8.cn)
文章仅用于分享个人学习成果与个人存档之用,分享知识,如有侵权,请联系作者进行删除。所有信息均基于作者的个人理解和经验,不代表任何官方立场或权威解读。
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。