GeoPandas中基于指定经纬度创建环形缓冲区及坐标异常排查

你的问题核心是直接在WGS84（EPSG:4326，经纬度坐标系）下用千米作为缓冲区半径导致的——这个坐标系的单位是「度」，不是米/千米，直接输入2km、2.1km会被当成2度、2.1度，换算成实际距离是几百公里，自然会出现坐标异常、地图无法正常显示的情况。

正确的操作步骤如下：

1. 初始化点数据并转换投影 #

先把经纬度点转换为米为单位的投影坐标系（比如UTM 19S，EPSG:32719，对应你提供的南纬33度、西经71.5度区域）：

import geopandas as gpd
from shapely.geometry import Point

# 原始经纬度点
center_point = Point(-71.5524818, -33.0679433)
# 创建GeoDataFrame，指定原始坐标系为WGS84
gdf = gpd.GeoDataFrame({'geometry': [center_point]}, crs="EPSG:4326")
# 转换为UTM投影（单位：米）
gdf_utm = gdf.to_crs("EPSG:32719")

2. 创建内外缓冲区 #

在投影坐标系下，直接用米作为半径创建缓冲区：

# 外缓冲区：2.1km = 2100米
buffer_outer = gdf_utm.buffer(2100)
# 内缓冲区：2km = 2000米
buffer_inner = gdf_utm.buffer(2000)

# 转换为GeoDataFrame，保留投影信息
buffer_outer_gdf = gpd.GeoDataFrame({'geometry': buffer_outer}, crs="EPSG:32719")
buffer_inner_gdf = gpd.GeoDataFrame({'geometry': buffer_inner}, crs="EPSG:32719")

3. 执行差值运算得到环形 #

# 外缓冲区减去内缓冲区，得到环形区域
ring_gdf = buffer_outer_gdf.overlay(buffer_inner_gdf, how='difference')

4. 可选：转换回WGS84坐标系 #

如果需要回到经纬度坐标系：

ring_gdf_wgs84 = ring_gdf.to_crs("EPSG:4326")

5. 绘制验证 #

import matplotlib.pyplot as plt

fig, ax = plt.subplots(figsize=(8,8))
ring_gdf.plot(ax=ax, color='lightblue', edgecolor='black')
gdf_utm.plot(ax=ax, color='red', markersize=50)
plt.show()

这样得到的环形缓冲区几何坐标会完全符合预期，地图也能正常显示。

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