### 批量地址转地理坐标（地理编码）的实现方法 在实际工作中，经常需要将一批地址信息转换为对应的地理坐标（即经纬度），这一过程称为**地理编码**。Python 提供了多种方式来实现这一需求，其中一种常见的方法是调用第三方地图服务提供商的 API，例如高德地图或百度地图的 API。 #### 使用高德地图 API 实现批量地址转坐标 高德地图提供了地理编码 API，可以通过 HTTP/HTTPS 协议访问远程服务接口，将结构化地址转换为经纬度坐标。以下是一个使用 Python 调用高德 API 进行批量地址转换的示例： ```python import requests import gcjtowgs # 用于将 GCJ-02 坐标转换为 WGS-84 坐标 my_key_list = ['XXX', 'XXX'] # 替换为你的高德 API 密钥 Lng_all = [] Lat_all = [] level_all = [] for i in range(0, len(data)): print(i) if i < 3000: my_key = my_key_list[0] else: my_key = my_key_list[1] try: address = data['注册地址'].iloc[i] url = f'https://restapi.amap.com/v3/geocode/geo?address={address}&output=json&key={my_key}' html = requests.get(url) html.close() info = html.json() # 获取高德返回的 GCJ-02 坐标 Lng_gcj = info['geocodes'][0]['location'].split(',')[0] Lat_gcj = info['geocodes'][0]['location'].split(',')[1] # 转换为 WGS-84 坐标 wgs = gcjtowgs.gcj02towgs84(float(Lng_gcj), float(Lat_gcj)) Lng_all.append(wgs[0]) Lat_all.append(wgs[1]) level_all.append(info['geocodes'][0]['level']) except: Lng_all.append('') Lat_all.append('') level_all.append('') # 将结果添加到原始数据中 data['Lng'] = Lng_all data['Lat'] = Lat_all data['level'] = level_all ``` 高德 API 返回的坐标为 **GCJ-02（火星坐标系）**，若需将其转换为国际通用的 **WGS-84 坐标系**，可以使用第三方库 `gcjtowgs` 完成坐标转换[^2]。 #### 多线程处理提升效率 对于大量地址数据的处理，可以使用 Python 的 `multiprocessing` 模块创建多线程处理机制，以加快请求速度并提高效率： ```python from multiprocessing import Pool if __name__ == '__main__': start = time.time() pool = Pool(20) urls = urlList results = pool.map(scrape, urls) csvfile = open('temp.csv', 'w') spamwriter = csv.writer(csvfile, delimiter=',', quotechar=',', quoting=csv.QUOTE_MINIMAL) print('over') print(len(results)) for index in range(len(results)): line = results[index] if index == 0: spamwriter.writerow(['', 'count_address', 'lng_wgs', 'lat_wgs']) else: if line: spamwriter.writerow([line['index'], line['address'].encode(encoding='utf-8'), line['lng_wgs'], line['lat_wgs']]) else: print('empty') print(time.time() - start) ``` 通过创建多线程池，可以并行处理多个地址请求，显著减少整体处理时间[^3]。 #### 使用百度地图 API 进行逆地理编码 若已有经纬度信息，需要获取对应的地址信息（即逆地理编码），可使用百度地图 API。以下是一个使用 Python 调用百度地图 API 的示例： ```python import pandas as pd import time import random # 读取经纬度数据 df = pd.read_excel('C:/VV/学习/读研/毕业论文/数据/中间稿/经纬度-congqing&sichuan.xlsx', skiprows=5000, nrows=5000) city_list = [] address_list = [] for index, row in df.iterrows(): lat = row['latitude'] lng = row['longitude'] city1, address1 = get_address(lat, lng, ak) print(city1, address1) city_list.append(city1) address_list.append(address1) sleep_time = random.uniform(0.04, 0.2) time.sleep(sleep_time) # 控制请求频率，避免超限 # 将结果写入原 Excel 文件 df['city'] = city_list df['address'] = address_list with pd.ExcelWriter('C:/VV/学习/读研/毕业论文/数据/中间稿/经纬度-congqing&sichuan.xlsx') as writer: df[['city', 'address']].to_excel(writer, startrow=5000, index=False, header=False) ``` 百度地图 API 对免费用户设置了每日请求限制（如 5000 次），因此在处理大量数据时需合理安排请求频率，避免超出限制[^4]。 ### 地理编码的应用场景 - **物流配送**：将客户地址转换为坐标，用于路径规划。 - **城市规划**：分析不同区域的地理分布情况。 - **商业分析**：将门店地址转换为坐标，用于可视化地图展示。 - **科研数据处理**：用于地理信息系统（GIS）中的空间分析。 ### 注意事项 - **API 限制**：多数地图服务提供商会对免费用户设置请求次数限制，需合理安排请求频率。 - **坐标系转换**：高德地图返回的坐标为 GCJ-02，若需与 GPS 数据兼容，应转换为 WGS-84 坐标系。 - **异常处理**：在调用 API 时，应加入异常处理机制，避免因单个请求失败导致整个程序中断。