点到多边形最近边距离计算时PostGIS与Geopandas结果不一致求助

问题根因 #

两者结果差异的核心原因是你在Geopandas中使用的EPSG:25832坐标系多边形WKT本身存在错误，和PostGIS中参与计算的多边形不是同一个几何对象：
EPSG:25832为UTM 32N投影坐标系，对应东经6°到12°、北纬0°到84°的区域，该坐标系下东经6.7°、北纬51°附近的投影坐标，东向（第一数值）为34万左右的6位数值，北向（第二数值）为566万左右的7位数值。你提供的fieldCoordinatesAsTextInWKTInEPSG25832内的多边形坐标，东向为562万左右的7位数值、北向为100万左右的6位数值，和点坐标POINT(341849.5 5661622.5)完全不在同一个空间区域，明显是投影转换出错的无效WKT。

结果验证 #

• PostGIS返回的4.32797817574802是正确结果：对比EPSG:4326坐标系下的点和多边形坐标，点位置为(6.741952,51.08428)，多边形的最近边就在附近，实际距离换算为米后和PostGIS输出完全匹配。
• Geopandas返回的3.8954865274727614e-05米是错误结果：因为输入的多边形WKT本身位置错误，计算自然无法得到正确值。

正确计算步骤 #

1. 直接通过代码完成坐标系转换，避免手动转换WKT出错，示例代码如下：

from shapely import wkt
from shapely.ops import transform
from pyproj import Transformer

# 加载EPSG:4326坐标系的原始多边形和点
poly_4326 = wkt.loads("POLYGON((6.741879696309871 51.08423775429969,6.742907378503366 51.08158745820981,6.746964018740842 51.08233499299334,6.746152690693346 51.08440763989611,6.741879696309871 51.08423775429969))")
pt_4326 = wkt.loads("POINT(6.7419520458647835 51.08427961641239)")

# 定义4326转25832的坐标转换器，always_xy=True确保坐标顺序为经度/东向在前
transformer = Transformer.from_crs("EPSG:4326", "EPSG:25832", always_xy=True)
# 完成几何对象的坐标转换
poly_25832 = transform(transformer.transform, poly_4326)
pt_25832 = transform(transformer.transform, pt_4326)

# 计算距离，单位为米，结果与PostGIS输出一致
print(poly_25832.distance(pt_25832))

2. 若必须手动使用WKT参与计算，转换完成后先比对坐标范围是否和点坐标处于同一区间，确认WKT有效后再计算。

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