从PPP到国家坐标系：基于速度场的希腊HGRS1987转换方法

精密单点定位（PPP）技术使GNSS接收机无需依赖本地参考站即可获得厘米级坐标，但PPP坐标并不自动等同于测量员、地籍机构或工程项目所使用的国家坐标系坐标。现代GNSS定位基于当代地球参考框架，而许多国家测绘系统建立于数十年前，与特定的参考框架实现和历元绑定。希腊的HGRS1987便是典型例子：它通过HTRS07与ETRS89关联，官方转换在2007.5历元固定，全国RMS约8.3厘米。 问题在于，PPP解算得到的是观测历元下ITRF框架的坐标，而该框架本身相对于旧的国家框架在运动。简单应用静态转换远远不够，必须先将坐标关联到正确的参考框架和历元，并顾及站点在不同历元间的运动。在变形区域，速度场不是次要改正，而是转换本身的一部分。希腊所在的愛琴海地区构造变形强烈且空间变化显著，因此需要能表示地壳三维运动的速度模型，并能提供任意点的速度。 本研究开发并测试了一套工作流：结合空间插值速度场、连续参考框架转换和最终投影坐标转换，而非提出新的七参数转换。速度模型基于Briole等人（2021）发布的GPS速度场（2000-2020年观测，希腊及邻国329个站点，其中希腊282个），速度相对于ITRF2014。为在任意国家三角网（NTN）标石处获取速度，研究采用双线性插值，选取302个站点构建模型，27个站点用于交叉验证，并注意地理分布以确保验证覆盖全国。 交叉验证结果令人鼓舞：预测速度与官方速度的RMS差异为东分量0.6毫米/年、北分量0.7毫米/年、垂直分量1.0毫米/年、水平分量0.9毫米/年。考虑到希腊非均匀的地球动力学行为和不可避免的插值误差，作者认为这些水平可接受。速度模型用于将坐标从观测历元传播到国家转换定义的历元，若存在系统偏差，会累积到转换坐标中；反之，足够真实的速度模型可使相隔多年的PPP坐标在同一局部基准下一致比较。

图1展示了模型基础的水平速度场，其空间变化明显：箭头并非单一刚性运动，而是大小和方向在希腊地区变化。这正是单一静态转换无法满足严格PPP到局部基准工作流的原因。 一旦估计出感兴趣点的速度，转换通过三个主要阶段进行。首先，PPP解算结果表达在与PPP处理相关的参考框架和观测历元下；若初始框架不是ITRF2014，则需将速度信息一致地转换到ITRF2014，因为采用的希腊速度场定义在该框架下。其次，将坐标转换到2007.5历元的ETRF2005——这是希腊国家转换关联的区域框架和历元，研究使用EUREF的ETRF/ITRF坐标转换工具完成此步骤。

最后，将ETRF2005坐标通过官方转换应用到HGRS1987投影坐标系。整个流程的关键在于速度场的准确性和框架转换的一致性。该方法不仅适用于希腊，也为其他动态变形区域的国家坐标转换提供了参考。欲了解完整技术细节、公式和验证结果，请点击原文链接查看。