PROJ坐标转换中的垂直基准面处理问题分析

问题背景

在PROJ坐标转换工具中，当使用EPSG:9464（GDA2020地理3D坐标系）向EPSG:9463（GDA94地理3D坐标系）转换时，发现了一个关于垂直高度处理的潜在问题。该转换过程应用了7参数赫尔默特变换，但未能正确保留AHD（澳大利亚高程基准）高度值。

技术细节分析

当前转换行为

当前PROJ生成的转换管道如下：

+proj=pipeline
  +step +proj=axisswap +order=2,1
  +step +proj=unitconvert +xy_in=deg +xy_out=rad
  +step +proj=cart +ellps=GRS80
  +step +proj=helmert +x=0.06155 +y=-0.01087 +z=-0.04019 +rx=-0.0394924
        +ry=-0.0327221 +rz=-0.0328979 +s=-0.009994 +convention=coordinate_frame
  +step +inv +proj=cart +ellps=GRS80
  +step +proj=unitconvert +xy_in=rad +xy_out=deg
  +step +proj=axisswap +order=2,1

这个转换管道完整地应用了7参数变换，包括对Z坐标（高度）的变换。然而，当坐标系包含基于大地水准面的高程（如AHD）时，这种处理方式是不正确的。

正确转换行为

当使用复合坐标系（如EPSG:28356+EPSG:5711）进行转换时，PROJ会生成不同的转换管道：

+proj=pipeline
  +step +inv +proj=utm +zone=56 +south +ellps=GRS80
  +step +proj=push +v_3
  +step +proj=cart +ellps=GRS80
  +step +proj=helmert +x=0.06155 +y=-0.01087 +z=-0.04019 +rx=-0.0394924
        +ry=-0.0327221 +rz=-0.0328979 +s=-0.009994 +convention=coordinate_frame
  +step +inv +proj=cart +ellps=GRS80
  +step +proj=pop +v_3
  +step +proj=utm +zone=56 +south +ellps=GRS80

这里的关键区别在于使用了push +v_3和pop +v_3操作，这些操作将垂直分量（Z值）从转换过程中排除，从而保留了原始的大地水准面高度。

问题本质

这个问题的核心在于坐标转换中如何处理基于大地水准面的高程值：

1. 对于纯水平基准转换（如GDA94到GDA2020），垂直分量应该保持不变
2. 大地水准面模型（如AHD）已经考虑了基准差异，不应再应用额外的垂直变换
3. 当前实现未能区分这种情况，导致垂直分量被错误地转换

解决方案建议

PROJ应该在处理包含大地水准面高程的坐标转换时：

1. 自动识别垂直基准类型
2. 对于基于大地水准面的高程，在转换管道中插入push/pop操作
3. 仅对椭球高（ellipsoidal height）应用完整的3D变换

这种处理方式将确保：

• 水平坐标正确转换
• 大地水准面高程保持不变
• 转换结果符合测量规范要求

实际影响

这个问题会影响所有使用PROJ进行澳大利亚坐标系转换的应用，特别是那些需要精确高程数据的应用场景，如：

• 测绘工程
• 水文测量
• 地质调查
• 基础设施建设

正确的处理方式对于确保高程数据的连续性和一致性至关重要。