地址转坐标:利用高德API进行批量地理编码

我们会遇到一些需要批量把地址数据转换成坐标的问题,一个个查的话通常会比较麻烦,那有没有比较方便的方法呢,我们可以利用高德的地理编码/逆地理编码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)

文章仅用于分享个人学习成果与个人存档之用,分享知识,如有侵权,请联系作者进行删除。所有信息均基于作者的个人理解和经验,不代表任何官方立场或权威解读。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值