GeoPandas/Shapely计算EPSG:4326多边形相交面积的单位疑问

首先直接回答两个核心问题：

• 你在EPSG:4326坐标系下直接调用.area得到的结果单位是平方度，这个单位不具备实际的面积参考价值：经纬度是角度单位而非长度单位，1经度对应的实际地面距离从赤道处的约111km随纬度升高逐渐缩小到极点处的0，因此平方度无法直接对应平方米、平方公里这类通用面积单位，这就是计算结果不符合预期的根本原因。
• 直接在EPSG:4326下计算相交面积的方法确实无法得到正确的实际面积，你需要根据使用场景选择以下两种正确计算方式：

方式一：转换为等积投影后计算（精度稳定，适合批量空间运算） #

根据你的研究区域范围选择匹配的等积投影坐标系，将两个GeoDataFrame统一转换到该投影下后再执行相交、面积计算操作，此时输出的面积单位为投影对应长度单位的平方：

• 研究区覆盖全球：可选择EPSG:6933（WGS 84 全球等积圆柱投影），输出面积单位为平方米
• 研究区为局部区域：选择对应范围的等积投影即可，例如中国区域可选择对应分带的高斯-克吕格等积投影，小范围场景下使用对应UTM投影的面积误差也在可接受范围内
  参考代码：

# 以转换为EPSG:6933全球等积投影为例
left_eqc = left.to_crs("EPSG:6933")
right_eqc = right.to_crs("EPSG:6933")
# 计算相交面积，结果单位为平方米
intersect_area = left_eqc.iloc[i].geometry.intersection(right_eqc.iloc[idx].geometry).area

方式二：测地线面积计算（无需转投影，适合跨区域快速计算） #

Geopandas内置了基于WGS84椭球的测地线计算能力，不需要手动转换投影，直接在EPSG:4326坐标系下就能得到符合实际的高精度面积结果：

# 先获取相交几何
intersect_geom = left.iloc[i].geometry.intersection(right.iloc[idx].geometry)
# 直接计算测地线面积，结果单位为平方米
intersect_area = intersect_geom.geodesic_area(ellps="WGS84")

# 如果是低于0.9版本的Geopandas没有内置geodesic_area方法，可通过pyproj实现
import pyproj
geod = pyproj.Geod(ellps="WGS84")
intersect_area = abs(geod.geometry_area_perimeter(intersect_geom)[0])

注意：如果你的数据覆盖高纬度区域、或者跨多个经度带，不要使用局部UTM投影计算面积，投影变形会带来明显误差，优先选择全球等积投影或者测地线计算方法。

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