PostGIS 3.3与2.2版本ST_Transform坐标转换结果不一致求助

原因分析 #

• PROJ库版本迭代：PostGIS的坐标转换依赖PROJ库，从2.2.1到3.3.3，配套的PROJ库版本跨度极大（2.2.1对应PROJ4.x系列，3.3.3对应PROJ8.x及以上）。PROJ新版本对Web墨卡托（900913/EPSG:3857）与WGS84（4326）的转换算法做了精度优化，修正了旧版本中近似转换的误差，这是结果差异的核心原因。
• SRID映射规范调整：900913是EPSG:3857的非官方别名，旧版本PostGIS对该SRID的参数映射存在简化处理，而新版本严格遵循EPSG官方定义的投影参数，转换逻辑更严谨，输出结果更精准。

解决办法 #

如果需要完全匹配旧版本的转换结果，可通过以下方式实现：

• 自定义旧版转换参数：使用PROJ旧格式参数字符串手动指定转换逻辑，模拟PostGIS2.2.1的行为，示例查询：

SELECT ST_AsText(
  ST_Transform(
    '010100002031BF0D00CDCCCC1CD50E5F41EC51B81E6583F2C0'::geometry(Geometry,900913),
    '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs'::text
  )
);

• 接受新版本精准结果：如果业务对精度要求更高，新版本的结果更符合官方投影标准，建议逐步迁移适配新结果——旧版本的转换存在系统误差，新版本的输出才是更准确的坐标值。

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