论文推荐｜马下平：ITRF中GNSS/SLR并址站归心基线的“一步解”

《测绘学报》

构建与学术的桥梁拉近与权威的距离

马下平

西安科技大学测绘科学与技术学院,陕西西安710054

收稿日期：2017-02-17；修回日期：2017-11-17

基金项目：国家自然科学基金面上项目（46134017）；西安科技大学博士启动金（2016QDJ005）；西安科技大学培育基金（201714）

第一作者简介：马下平(1984-),男,博士生,讲师,研究方向为大地测量数据处理。E-mail:celiang0321@163.com

摘要：提出将SLR望远镜的参考点和两轴偏差作为未知参数，在ITRF中联合并址站归心测量中GNSS基线网和地面网观测量（水平方向、垂直角和边长），建立SLR站观测设备的参考点与观测标志、观测标志之间、参考点和两轴偏差与其他未知参数之间的多种约束条件来求解归心基线的“一步解”。利用“一步解”解算出“陆态网络”中北京、昆明和西安3个GNSS/SLR并址站在ITRF2014中的归心基线及其协方差阵。结果显示：归心基线的中误差优于2mm，与已有分步解相比，差值不超过2mm；水平轴和垂直轴之间的偏差分别为3.8、0.7和3.6mm，中误差分别为1.3、1.2和1.3mm。

OneStepSolutiontotheLocalTieVectorsofGNSS/SLRinITRF

MAXiaping

Abstract:Thispaperproposesaone-stepsolutionoflocaltievectors,thesolutiontakesRP(ReferencePoint)andaxesoffsetsofSLRtelescopeasunknownparameters,GNSSbaselinesnetworkandconventionalterrestrialobservation(horizontal,verticalangleanddistance)arecombinedinITRF(InternationalTerrestrialReferenceFrame),andmultiplyconstraintconditsinBeijing,Kunming,Xi'anofCMONOC(CrustalMovementObservationNetworkofChina)(rootmeansquareerror),,andtheoffset,0.7mm,3.6,1.2mm,1.3mm,respectively.

Keywords:localtievectorsonestepsolutionGNSS/SLRconditions

ITRF(internationalterrestrialreferenceframe)的解算是综合应用全球导航卫星系统GNSS、甚长基线干涉测量VLBI(verylongbaselineinterferometry)、卫星激光测距SLR(satellitelaserranging)和多普勒定轨与无线电定位技术DORIS(dopplerorbitdeterminationandradiopositioningintegratedonsatellite)4种空间大地测量技术的观测数据和并址站归心基线[1-4]。并址站归心基线作为联系ITRF中各观测技术的主要约束条件，已成为限制ITRF解算精度的关键因素之一。即使各观测技术的数据质量有了很大的提高，如果归心基线的精度没有相应地提高，ITRF解的精度也会受到一定程度的影响[5-6]。

并址站归心测量的实质就是获取ITRF中各观测设备参考点RP(referencepoint)的坐标及其协方差阵[7-9]。SLR望远镜的RP为绕垂直轴p的旋转中心与水平轴s和垂直轴p之间的最短距离的交点，如图1所示。GNSS的RP为其接收机天线的几何中心(如图2所示)，可以直接测定。而SLR望远镜的RP为其望远镜的旋转中心，并非一个可以直接观测到的实体点，因此必须设法以间接的方式得到。

图1SLR望远镜的参考点RP[9]

图选项

图2GNSS天线的参考点ARP[7]

图选项

对于SLR望远镜参考点的测定，通常是先在SLR望远镜周围布设控制点，采用常规三维控制测量精确测定各控制点坐标，通过极坐标法或前方交会测定固定在SLR望远镜上标志点的坐标[10-13]，也有在并址站设备上安装GNSS天线，通过静态或动态技术获取标志点的坐标[14-15]，以标志点坐标为虚拟观测量解算SLR望远镜参考点的坐标及其方差阵，目前解算方法主要有两种：分步解算方式和构建坐标转换模型的方式。分步解法主要有：①将标志点的三维空间直角坐标转换为站心坐标，通过曲线拟合求出绕垂直轴和绕水平轴旋转中心的站心坐标，再将旋转中心的站心坐标逆变换得到三维空间直角坐标[10-12,16-18]；②以标志点三维空间直角坐标作为虚拟观测量，建立两类约束条件获取参考的位置，具体可参考文献[13]；③为了提高解算精度，可引入额外的约束条件，如距离相等条件、旋转平面一致性条件、共线性条件、半径相等条件、旋转中心距离相等条件、旋转轴正交性条件、轴偏差相等[7,19-20]；④利用参考点与绕垂直轴和水平轴旋转中心之间的空间向量关系，获取参考点的位置，间接得到两轴之间的偏差[21]。

除分步解法外，目前最新的解算方法是通过构建包括参考点、轴偏差等未知参数的并址站设备的独立坐标系统与地面坐标系统之间的转换模型间接得到并址站参考点位置[15,20,22,24]，但这种方法同样是在独立坐标系中进行分步解算，而且解算时需要知道一些参数的值(绕垂直轴和水平轴的旋转角)等。

总结以上文献，仍存在着以下方面的特点与不足：①将观测量和确定参考点和两轴偏差分开计算，存在站心系和空间直角坐标系多次换算的问题，造成最终结果的不准确和最终精度的损失；②忽略或是未考虑观测标志与参考点或几何设备所应满足的几何约束条件，如忽略两轴偏差的影响，将会造成结果0.5mm的精度损失[19]；③虽然可以建立坐标转换模型获取参考点的位置，但同样需要进行分步解算，而且事先需要提供一些参数的值。

不同的文献有着不同的解算方法，不同的方法所得出结果各有不同，因此本文力图在一个坐标框架(ITRF)中，将SLR望远镜的参考点、两轴偏差等参数视为未知参数，联合GNSS基线和地面网观测数据，引入满足条件的各种几何约束条件方程，提出并址站归心基线的“一步解”。“一步解”经平差后可直接得到并址站归心基线及其方差协方差阵。

1“一步解”的数学模型1.1地面观测量和GNSS基线网的观测方程

1.1.1地面观测量的观测方程

地面观测量一般是指经全站仪通过极坐标法、前方交会等方法测量的水平方向、垂直角和边长。通过建立站心坐标可以实现与空间直角坐标之间的转换，可以将地面观测量归算至空间直角坐标系中。

根据法线站心坐标和空间直角之间的转换关系，测点j在以测站i为原点的站心坐标系中的坐标向量xj和空间直角坐标向量xi和xj关系为

(1)

式中，xj=(Nj,Ej,Uj)，Nj、Ej、Uj为站心坐标分量；同样，xi=(Xi,Yi,Zi),xj=(Xj,Yj,Zj)，Xi,Yi,Zi和Xj,Yj,Zj为测站i和测点j的空间直角坐标分量，本文后面的坐标分量统一以坐标向量的形式表示。其中，R为旋转矩阵

(2)

式中，(Bi,Li)为测站i的纬度和经度。

由极坐标和站心坐标之间的关系可以得到测站i至测点j的水平方向αij、垂直角βij和边长Sij在站心系中的观测方程为

(3)

(4)

(5)

式中，(εαij,εβij,εSij)为αij、βij和Sij的观测误差；gi为测站i的定向角；Ka为测站的大气折光；Ra为地球曲率。

对式(3)—式(5)进行线性化后可以得到αij、βij和Sij在空间直角坐标系中的观测方程见式(6)—式(8)

(6)

式中，变量上方带0均表示近似值；δj,δYj,δZj)为测点j的空间直角坐标向量改正数，δXj,δYj,δZj为各分量的改正数；δgi为测站i定向角改正数。其中，

特别说明的是，旋转矩阵R中3列向量对应的是δ的3个元素。

(7)

式中，δKa为大气折光参数的改正数。其中

(8)

式中，

1.1.2GNSS基线网向量的观测方程

经GNSS软件解算后的基线网向量视为观测量，则GNSS基线ij的观测方程为

(9)

式中，Δxij和εΔxij为基线ij的坐标向量及其误差，Δxij=(ΔXij,ΔYij,ΔZij)，εΔx=(εXij,εYij,εZij)，(ΔXij,ΔYij,ΔZij)和(εXij,εYij,εZij)为基线ij在X,Y,Z方向的3个分量及其观测误差；xi0=(Xi0,Yi0,Zi0)，xj0=(Xj0,Yj0,Zj0)。

式(6)—式(9)以向量形式表示的观测方程为

(10)

式中，x1为4类观测方程中所包括的未知参数向量改正数；A为4类观测方程系数阵组成的系数向量，l为4类观测方程常数组成的向量。

1.2约束条件

1.2.1观测标志所满足的条件

观测标志绕某一旋转轴旋转一周时，可构成一个球面和平面方程。当观测标志k绕垂直轴p旋转构成的平面法向量和球面半径分别为np和rp，同样绕水平轴s旋转构成的平面法向量和球面半径分别为ns和rs，则得到的两组平面和球面方程分别为

(11)

(12)

式中，xk为观测标志k的坐标向量；xp和xs为垂直轴p和水平轴s旋转中心的坐标向量。

1.2.2距离和半径相等条件

在旋转过程中，安装在SLR望远镜上的两个观测标志之间的距离是保持不变的。如：当观测标志1和2绕旋转轴第k1和k2次旋转时分别构成的距离相等

(13)

同样，同一观测标志绕同一旋转轴多次旋转时具有相同的中心和半径，如：观测标志1和2绕垂直轴p第k1次和k2次旋转时，构成的球面半径相等。则有

(14)

同样，绕水平轴s进行旋转时，得到

(15)

1.2.3法向量一致性条件

不同的观测标志绕同一旋转轴旋转时构成的平面法向量相等。如标志1和2绕垂直轴p构成的法向量np1、np2和水平轴s构成的法向量ns1、ns2分别相等，且为单位向量

(16)

(17)

同样，同一标志绕同一旋转轴多次旋转构成的法向量也相等。如：绕p轴和s轴进行第k1次和k2次旋转构成的法向量分别相等。

(18)

(19)

1.2.4参考点、轴偏差和旋转中心之间的约束条件

如图3所示，x0为参考点RP的坐标向量，xv为RP至s轴的最短距离与s轴交点的坐标向量，d为两轴的偏差。当绕垂直轴p旋转时可构成线性条件为

(20)

图3参考点与两轴偏差之间的关系

图选项

同样，绕水平轴s旋转时也可构成线性条件为

(21)

式中，λp和λs为垂直轴p和水平轴s的比例参数。

垂直轴p和水平轴s之间的最短距离必正交于两轴，因此，两轴之间的法向量n同时垂直于p和s，得出

(22)

由于x0和xv具有同一个法向量，得出线性关系

(23)

同样，两轴的偏差向量d为参考点x0和投影点的xv的关系为

(24)

式(13)—(21)和式(23)—(24)经过线性化后，这些约束条件可以表示为

(25)

式中，2为约束条件中引入的未知参数向量改正数，C2和w为其系数阵和常数项。

在观测方程式(10)中引入约束方程式(25)，经过最小二乘平差后可得到参考点的坐标及其方差阵。联立式(10)和式(25)，得

(26)

式(26)中涉及GNSS基线观测量、水平方向和垂直角、边长共3类观测量，对于不同类观测量采用经验公式定权是不够准确的。为了提高这3类观测量的方差协方差Σh估计的精度，本文采用赫尔默特方差分量验后估计定权法。具体的计算公式和步骤可参考广义测量平差的内容。

根据附有约束条件的间接平差可得法方程为

(27)

式中，N=ATΣh-1A。则未知参数向量2及其方差协方差阵

(28)

(29)

式中，

假设式(26)中有u个未知参数，其中1有t个独立未知参数(t＜u)；约束条件共有f个；误差方程的改正数为V，则的计算公式为

(30)

由于每一个SLR望远镜周围附近都有一个IGS站，因此可利用IGS站的ITRF2014成果作为基准约束，可将式(27)的平差成果计算至ITRF2014中。

2成果计算与分析

本文采用“一步解”处理了2011年笔者实际观测的“陆态网络”北京(BJFS-SLR)、昆明(KUNM-SLR)和西安(XIAA-SLR)3个并址站归心测量数据，获取了3个并址站的归心基线及其协方差阵，SLR并址站外业观测方案具体可参考文献[12-13]。GNSS基线网和地面观测量的定权采用方差分量估计法，具体可参考文献[25]。计算步骤如下：

首先用GAMIT软件解算出SLR望远镜周围GNSS强制观测墩构成控制网的基线向量；其次列出GNSS网基线向量、地面三维观测量的观测方程，再引入本文总结的多类约束条件方程，利用式(28)和式(29)解算出“陆态网络”3个SLR并址站的归心基线及其协方差矩阵。解算结果如下：

表1给出3个并址站的归心基线及其中误差；表2、表3和表4为3个并址站归心基线的方差协方差阵；表5给出与文献[12]的比较结果；表6给出与ITRF2014的比较结果；表7和表8分别给出两个并址站的两轴偏差及其中误差。

表13个并址站的归心基线及其中误差

站点ΔX/mΔY/mΔZ/mMX/mmMY/mmMZ/mmSLR-BJFS-16.5167118.3176-146.28350.51.00.9SLR-KUNM-20.2144-18.854345.77992.11.41.9SLR-XIAA-14.864714.6906-28.08021.00.30.3注：BJFS和KUNM为SLR站并址的IGS站。XIAA为“陆态网络”站点。

表选项

表2BJFS站归心基线的方差协方差阵

m2站点σΔX2/σΔYΔX/σΔZΔXσΔXΔY/σΔY2/σΔZΔYσΔXΔZ/σΔYΔZ/σΔZ2SLR-BJFS2.556e-7-1.03776172344734e-77.00495601096625e-8-1.03776172344735e-79.27569196990275e-7-6.07599483204138e-77.00495601096627e-8-6.07599483204136e-77.40643113442457e-7注：σΔX2、σΔY2、σΔZ2分别为归心基线向量各分量的方差；σΔXΔY、σΔXΔZ、σΔYΔZ、σΔYΔX、σΔZΔX、σΔZΔY分别为归心基线分量各分量间的协方差。

表选项

表3KUNM站归心基线的方差协方差阵

m2站点σΔX2/σΔYΔX/σΔZΔXσΔXΔY/σΔY2/σΔZΔYσΔXΔZ/σΔYΔZ/σΔZ2SLR-KUNM4.56612652089061e-64.82316732063208e-72.96876737075803e-74.82316732063209e-71.936e-6-1.199e-62.96876737075803e-7-1.199e-63.79381263579369e-6

表选项

表4XIAA站归心基线的方差协方差阵

m2站点σΔX2/σΔYΔX/σΔZΔXσΔXΔY/σΔY2/σΔZΔYσΔXΔZ/σΔYΔZ/σΔZ2SLR-XIAA9.88298843603568e-7-2.40588580169321e-72.43266817097115e-7-2.40588580169324e-79.33790564321807e-8-5.64534063719309e-82.43266817097115e-7-5.64534063719304e-89.27854101664978e-8

表选项

表5“一步解”和文献[12]结果的差值(2013)

mm站点差值ΔXΔYΔZSLR-BJFS-0.10.20.0SLR-KUNM1.61.71.5SLR-XIAA0.9-1.2-1.2

表选项

表63个并址站归心基线与ITRF2014的比较

mm站点差值ΔXΔYΔZSLR-BJFS-3.59.8-0.3SLR-KUNM19.6-32.3-19.1SLR-XIAA———

表选项

表73个SLR并址站两轴之间的偏差

mm并址站dXdYdZd北京SLR3.50.81.03.8昆明SLR0.10.3-0.60.7西安SLR1.52.2-2.43.6

表选项

表83个SLR并址站两轴之间偏差的中误差

mm并址站dXdYdZd北京SLR0.51.00.61.3昆明SLR0.80.50.91.2西安SLR1.00.80.31.3

表选项

由以上的结果可知：

(1)3个并址站归心基线中误差基本位于2mm左右，其中误差最大为2.1mm。

(2)本文“一步解”得出的归心基线与Gong的结果相比，基线向量的差值分量小于2mm，差值最大为1.7mm。出现差值的主要原因是Gong采用分步解法，而且未考虑两轴偏差、观测标志与参考点满足的各种几何条件，因此从理论上本文的算法更加严密，结果更加可靠。

(3)与ITRF2014参考值相比，差值普遍较大，其中差值最大为32.3mm。主要原因在于IFRF2014解算时未考虑GNSS、VLBI、SLR、DORIS这4种技术之间的系统偏差，导致本文测定的并址站归心基线与ITRF2014采用值的差异较大。

(4)3个SLR望远镜两轴之间的偏差最大为3.5mm，其中中误差最大不超过2.0mm。

3结论

针对已有的并址站归心基线解算中需要分步解算、涉及多次坐标转换和未充分考虑多种约束条件等问题，本文提出了并址站归心基线的“一步解”。“一步解”是将SLR望远镜参考点和轴偏差引入至约束方程中，按照附有约束条件的间接平差直接获取归心基线及其协方差阵，理论严密、数学模型简单、易于程序设计，经与已有算法的结果相比，差值基本在2mm左右；利用北京、昆明和西安SLR望远镜的实测数据，计算出两轴的偏差，偏差分别为3.8、0.7、3.6mm，本文的解算方法和结果对ITRF解算具有一定的参考价值。

【引文格式】马下平。ITRF中GNSS/SLR并址站归心基线的“一步解”[J].测绘学报，2018，47(1)：64-70.DOI:10.11947/

往期精彩回顾

招聘|武汉大学宇航科学与技术研究院专职科研岗

关于稿件“时间”安排那些事儿~

重磅！新增博士、硕士学位授权点名单出炉，有你的母校吗

组建“自然资源部”的来龙去脉

权威|专业|学术|前沿

微信投稿邮箱|song_qi_fan@163.com

进群请备注：姓名+单位+稿件编号