在GeoPandas中绘制地面半径圆：结果不符合预期

核心原因 #

你用的平面投影（如EPSG:3395墨卡托）的buffer是基于平面几何计算的，而地球是球面，大距离下平面投影的距离会产生严重偏差——墨卡托投影在高纬度会拉伸距离，导致你生成的buffer转回到WGS84后，实际覆盖的球面距离远小于4800km，同时形状也会变形。

解决方案 #

直接在地理坐标系（EPSG:4326）下使用测地线缓冲区（geodesic buffer），GeoPandas的buffer方法支持geodesic=True参数，会基于球面大圆距离计算等范围区域，完全符合飞行距离的实际需求。

修改后的代码 #

import geopandas as gpd

lat = 33.975
lon = -118.405

world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))

# 构建点数据
data = [{"lat": lat, "lon": lon, "city": "SFA"}]
df = gpd.GeoDataFrame(data)

# 创建WGS84坐标系下的点几何
df["geometry"] = gpd.points_from_xy(df.lon, df.lat, crs="EPSG:4326")
df["center"] = df["geometry"]

# 直接生成测地线缓冲区，单位为米（4800*1000）
df["geometry"] = df.geometry.buffer(4800 * 1000, geodesic=True)

# 绘图
ax = df.plot(facecolor='none', edgecolor='red', figsize=(15,15))
df.center.plot(ax=ax, color="r")
world.plot(ax=ax, alpha=.5)

关键说明 #

• geodesic=True会让GeoPandas使用球面几何计算缓冲区，确保每个点到中心的大圆距离严格等于你设定的4800km，不会受投影变形影响。
• 无需来回转换投影，直接在WGS84下操作即可，避免了平面投影带来的距离偏差。
• 之前尝试的haversine公式如果没生效，大概率是没有基于球面构建多边形，而测地线buffer是GeoPandas原生支持的可靠方案。

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