Python中Geopandas多边形米制长度计算错误排查与解决

错误原因 #

1. 坐标顺序颠倒：gpd.points_from_xy()要求传入的顺序是经度(lon)在前，纬度(lat)在后，你写反成df['lat'], df['lon']，导致生成的Point坐标完全错位，后续计算自然错误。
2. 直接用地理坐标系计算长度：你的原始数据是WGS84（EPSG:4326）地理坐标系，单位是度，不是米或英里。不同纬度下，1度对应的实际地面距离差异很大（比如赤道上1经度≈111km，北纬30度时1经度≈96km），用固定系数0.000621371转换完全不适用。

正确解决步骤 #

1. 修正GeoDataFrame创建（坐标顺序+指定坐标系） #

首先确保点的坐标顺序正确，同时明确指定原始数据的坐标系为EPSG:4326：

import geopandas as gpd
from shapely.geometry import Point, Polygon

# 正确创建GeoDataFrame：lon在前，lat在后，指定crs为WGS84
gdf = gpd.GeoDataFrame(df, geometry=gpd.points_from_xy(df['lon'], df['lat']), crs="EPSG:4326")

2. 投影到平面坐标系 #

将地理坐标系转换为平面投影坐标系（比如UTM坐标系，单位为米），这样计算的长度才是真实地面距离：

# 自动匹配当前区域的UTM投影坐标系
utm_crs = gdf.estimate_utm_crs()
# 转换坐标系
gdf_proj = gdf.to_crs(utm_crs)

3. 创建多边形并计算周长（米制） #

从投影后的点创建Polygon，其length属性就是多边形的周长（自动包含首尾点的直线距离）：

# 提取投影后的点
points_proj = list(gdf_proj.geometry)
# 创建多边形
polygon_proj = Polygon([(p.x, p.y) for p in points_proj])
# 计算周长（单位：米）
perimeter_m = polygon_proj.length
print(f"多边形周长（米）：{perimeter_m}")

简化写法（无需手动创建Polygon） #

也可以直接用GeoPandas的LineString闭合后计算长度，结果一致：

from shapely.geometry import LineString

# 将点连成线并闭合（首尾相连）
line = LineString(list(gdf_proj.geometry))
closed_line = line.buffer(0).boundary
# 计算周长
perimeter_m = closed_line.length

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