如何通过邮政编码创建GeoDataFrame？无需手动查询经纬度

首先明确：GeoPandas本身没有内置邮政编码转经纬度的功能，但可以结合第三方地理编码库实现自动匹配，无需手动查询。以下是具体操作步骤：

1. 规整DataFrame结构 #

你的原始数据是每行Name1; Postcode1; Name2; Postcode2;...的格式，先将数据转换为每行对应一个名称+邮政编码的规整结构：

import pandas as pd

# 示例原始数据，替换为你的实际DataFrame
df = pd.DataFrame({
    'raw_data': ['Name1; Postcode1; Name2; Postcode2', 'Name3; Postcode3']
})

# 拆分每行数据并重塑为长格式
split_data = df['raw_data'].str.split(';', expand=True)
# 生成成对的列名
columns = []
for i in range(0, split_data.shape[1], 2):
    columns.extend([f'Name_{i//2+1}', f'Postcode_{i//2+1}'])
split_data.columns = columns

# 转长格式，得到单一的name和postcode列
melted_df = split_data.melt(var_name='type', value_name='value')
melted_df['id'] = melted_df['type'].str.extract('(\d+)')
melted_df['category'] = melted_df['type'].str.split('_').str[0]
clean_df = melted_df.pivot(index='id', columns='category', values='value').reset_index(drop=True)
clean_df = clean_df.dropna().rename(columns={'Name': 'name', 'Postcode': 'postcode'})

2. 自动获取经纬度（地理编码） #

推荐使用geopy库的Nominatim（基于OpenStreetMap），免费无需API密钥：

from geopy.geocoders import Nominatim
from geopy.extra.rate_limiter import RateLimiter

# 初始化地理编码器，建议填写你的项目名称/邮箱作为标识
geolocator = Nominatim(user_agent="your_project_identifier")
# 添加限速器，避免请求频率过高被封禁
geocode = RateLimiter(geolocator.geocode, min_delay_seconds=1)

# 对邮政编码进行地理编码
clean_df['location'] = clean_df['postcode'].apply(geocode)
# 提取经纬度信息
clean_df['lat'] = clean_df['location'].apply(lambda loc: loc.latitude if loc else None)
clean_df['lon'] = clean_df['location'].apply(lambda loc: loc.longitude if loc else None)

# 过滤匹配失败的行
clean_df = clean_df.dropna(subset=['lat', 'lon'])

3. 转换为GeoDataFrame #

用你熟悉的方式创建GeoDataFrame即可：

import geopandas as gpd

gdf = gpd.GeoDataFrame(
    clean_df,
    geometry=gpd.points_from_xy(clean_df['lon'], clean_df['lat']),
    crs='EPSG:4326'
)

# 可选：保存为地理数据文件
gdf.to_file('postcode_points.geojson', driver='GeoJSON')

注意事项 #

• 若处理大量邮政编码，建议使用付费地理编码服务（如Google Maps Geocoding API），速度更快、准确率更高，但需API密钥；
• Nominatim有严格的请求频率限制，必须添加RateLimiter；
• 部分地区的邮政编码可能无法匹配到精准位置，需后续人工核对补充。

内容的提问来源于stack exchange，提问作者p_vall_do