PROJ坐标转换中GDA94到GDA2020水平偏移问题的技术解析

问题背景

在PROJ地理坐标转换工具中，用户发现从GDA94坐标系转换到GDA2020坐标系时，当使用EGM2008或EGM96垂直参考系统时，水平偏移量未能正确应用。这一现象特别出现在使用--3d选项进行三维坐标转换时。

现象描述

正常情况下，GDA94到GDA2020的转换应包含约1米左右的水平偏移。例如：

400500 6500000 60 → 400500.998 6500001.495 59.922

但当目标坐标系包含EGM2008(EPSG:3855)或EGM96(EPSG:5773)垂直参考时：

400500 6500000 60 → 400500.000 6500000.000 92.044

水平偏移量丢失，仅垂直分量被转换。

技术原因分析

这一现象的根本原因在于坐标转换的参考框架选择：

1. WGS84的中间角色：EGM2008和EGM96垂直参考系统都是基于WGS84椭球定义的。PROJ在进行这类转换时，会自动引入WGS84作为中间参考框架。

2. GDA与WGS84的关系：在EPSG定义中，GDA94和GDA2020都被视为与WGS84等效（尽管实际存在约3米的差异）。因此转换路径变为：

   • GDA94 → WGS84（视为等效，无偏移）
   • WGS84 → GDA2020（视为等效，无偏移）
   • 最后应用垂直转换

3. 转换路径选择：PROJ的转换路径选择算法会优先考虑官方定义的转换关系，当存在多个等效精度的转换路径时，可能不会自动选择包含水平偏移的路径。

解决方案

对于需要同时保持水平偏移和EGM垂直转换的场景，可以采用分步转换方法：

1. 首先进行GDA94到GDA2020的三维转换
2. 然后进行WGS84到EGM垂直系统的转换
3. 最后组合结果

示例命令：

echo "400500 6500000 60" | cs2cs --3d EPSG:28350 +to GDA2020 | cs2cs --3d "WGS 84" "WGS 84 + EGM2008 height" | cs2cs "GDA2020 + EGM2008 height" "GDA2020 / MGA zone 50 + EGM2008 height"

技术建议

1. 了解转换路径：使用projinfo命令查看所有可能的转换路径，选择最适合的路径。

2. 精度权衡：在需要高精度水平转换的场景下，可能需要牺牲部分垂直精度，或采用自定义转换路径。

3. 未来改进：PROJ团队可以考虑增加转换优先级设置选项，让用户能够指定优先保持水平或垂直精度。

总结

这一问题反映了地理坐标转换中参考框架选择的复杂性。在实际应用中，理解各种参考系统之间的关系和转换路径选择机制至关重要。对于高精度要求的应用场景，建议仔细验证转换结果，必要时采用分步转换方法确保各分量转换的正确性。