如何将EPSG:4326的Shapely几何对象投影至EPSG:3857？

Shapely几何对象本身不携带坐标系信息，要完成投影转换并计算面积，你可以借助pyproj库（GeoPandas底层依赖的坐标转换库）来实现，具体操作如下：

步骤1：导入所需模块 #

确保已安装pyproj（通常安装GeoPandas时会自动配套安装），导入相关工具：

from pyproj import Transformer, CRS

步骤2：转换坐标系并计算面积 #

针对你手中的EPSG:4326格式Shapely几何，执行以下代码完成转换和面积计算：

# 定义源坐标系（4326）和目标坐标系（沿用你之前使用的3857）
src_crs = CRS.from_epsg(4326)
target_crs = CRS.from_epsg(3857)

# 创建转换器并转换几何对象，always_xy=True确保经纬度顺序正确映射
transformer = Transformer.from_crs(src_crs, target_crs, always_xy=True)
projected_geom = transformer.transform(feat.geometry)

# 计算转换后几何的面积
surface = projected_geom.area

或者你可以用更简洁的transform_geom函数，效果完全一致：

from pyproj import transform_geom, CRS

src_crs = CRS.from_epsg(4326)
target_crs = CRS.from_epsg(3857)
projected_geom = transform_geom(src_crs, target_crs, feat.geometry, always_xy=True)
surface = projected_geom.area

额外提示 #

• always_xy=True是必须的：EPSG:4326的坐标顺序是经度在前，纬度在后，而多数投影坐标系遵循XY对应横向/纵向的逻辑，这个参数能避免坐标顺序混乱导致的计算错误。
• 农田面积计算建议使用等面积投影：EPSG:3857（Web墨卡托）在高纬度区域面积失真严重，如果追求精度，可根据地块所在位置选择对应的UTM分区坐标系（比如EPSG:32633对应北半球33区），或使用当地官方的等面积投影。

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