基于CORS的实时坐标进行高精度转换的探究

基于CORS的实时坐标进行高精度转换的探究

陈建玺

中煤航测遥感集团有限公司，陕西西安，710000

摘　要：

现今区域坐标转换将GPS测量成果转化为1980年国家大地坐标系或CGCS2000国家大地坐标系下的平面直角坐标，在测量中及其重要。CORS系统的坐标转换虽存在许多不足，但基于CORS能够实现双向数据通信。根据狄洛尼(Delaunay)三角形构建基准站网络的特点，提出了通过实时获取坐标转换参数而获得高精度平面坐标的方法，从数据分析和理论的角度进行了可行性探究。

关键词：CORS;狄洛尼三角形;坐标转换

一、引言

GPS测量作为由传统测绘向新型技术测绘革新的一种方式，它是由连续运行的卫星定位系统(Continuous　Operational　Reference　System，CORS)得到流动站的初始测量成果也就是WGS-84坐标，这就需要通过区域坐标转换将GPS测量成果转化为1980 年国家大地坐标系或CGCS2000国家大地坐标系下的平面直角坐标，以满足各种工程需要。CORS系统利用坐标转换参数进行坐标转换存在着诸多不足：对于一个面积较大的城市，由于只使用一套转换参数，造成有些区域的坐标转换精度很差[1]。如果能获取实时动态、最佳拟合测区的坐标转换参数，则能够获得高精度的转换坐标。

因此，实时获取流动站高精度平面直角坐标是一个急需解决的问题。本文提出一种实时更新高精度坐标转换参数的方法，就目前CORS理论、应用状况及软硬件条件进行可行性探究。

二、区域坐标转换

WGS-84坐标系与国家坐标系之间的三维坐标转换，主要是采用Bursa-Wolf模型[2](布尔莎七参数)或Molodensky(莫洛金斯基三参数)进行坐标转换。对于大区域坐标转换常用Bursa-Wolf模型，测区较小的采用Modolensky模型。

令某点的WGS-84坐标系下的坐标为(，，)T，其在1980年国家大地坐标系下的坐标为为()T。虽然我国2008年7月1日起正式启用CGCS2000国家大地坐标系，但是各地该坐标系下的成果较少，还未形成使用规模，这里仍然使用1980国家大地坐标系进行坐标转换，其布尔莎七参数坐标转换公式为：

(1)

式(1)中，()T为3个平移参数，(α，β，γ)T为3个旋转角参数，k为尺度参数。将WGS-84坐标系下的(，，)T转换为平面直角坐标的步骤为：①将WGS-84坐标系下的大地坐标()T转换为WGS-84坐标系下的空间直角坐标(，，)T；

②按公式(1)WGS-84坐标系下的空间直角坐标(，，)T；转换为1980国家大地坐标系下的空间直角坐标()T；

③依据1975国际椭球参数，将1980国家大地坐标系下的空间直角坐标()T转换为1980国家大地坐标系下的大地坐标()T；

④设置中央子午线经度，由Gauss投影正算公式(2)求得平面坐标(x，y)，与此同时CORS完成大地高到正常高的高程转换得到h。

三、坐标转换实验

3.1　蒲城七参数下的坐标转换

利用陕西CORS网对蒲城县奉先街道的5个二等控制点进行测量，点位分布如图1 所示，采样率为1s，平滑60次，并与已知坐标比较，分析其外符合精度。

图1　控制点分布图

表1　蒲城七参数下SDCORS测量结果与已知坐标残差(单位：m)

表1中可以看出，基于陕西省CORS系统利用蒲城七参数测得成果的外符合精度较差，北向残差为3cm～6cm，东向为7cm～13cm，高程方向上为2cm～6cm。测区观测环境较好，且实验操作规范严谨，考虑到因实验人员轻微晃动等引起的轻微误差，水平和垂直方向上仍存在较大差异，因而大区域的坐标转换无法满足精密工程建设的要求。

3.2　小区域七参数下的坐标转换

蒲城县奉先街道内控制点的点位分布如图1所示，经过分析发现由I002、I004和I005组成的三角形近似Delaunay三角形[3]且能覆盖整个测区，利用坐标转换软件反算出坐标转换参数，解算出另外2个点的坐标，与已知点坐标进行比较。所用数据为SDCORS测得5个点的WGS-84坐标与已知1980国家大地坐标，与表1中数据相同，保证了实验的严密性，表2是I001、I003的坐标残差。

表2　小区域坐标转换参数下的坐标成果与已知点坐标残差(单位：m)

由表2可知，坐标成果与已知点坐标的残差值较小，在水平和垂直方向上有较好的定位精度。

小区域坐标转换参数下，定位精度要好，但不易实现;而大区域的坐标转换参数，虽然易于实现，但转换精度相对较差。针对此情况，提出了基于CORS实时更新坐标转换参数获取高精度区域坐标转换的方法，并从理论上进行了可行性探究。

四、实时获取坐标转换参数实现高精度区域坐标转换的方法

4.1　CORS基准站网络组网形式与精度分布情况

我国大部分地区基准站间距离多为100km～200km，江苏省CORS基准站间平均距离为50km～80km，而广东省由于位于低纬度地区，电离层活动活跃，建立了高密度的CORS基准站网络，站间距离平均为30km～50km。当然，在相同的服务范围内CORS基准站距离越长，需要的基准站数量越少，建设成本也越低，与此同时，空间相关性差，内插精度低，初始化时间要长，影响可用性。

利用CORS进行动态测量，使用精度为10m的广播星历[4]，轨道误差对50km 基线最大影响为2.5cm，由于GPS系统的载波相位波长都在15cm以上(L1为19.03cm，L2为24.42cm)，对于距离小于50km基准站之间，这一误差不会影响模糊度的求解和初始化时间。而且，数据处理中心不对电离层延迟、对流层延迟等误差进行区分，统一集中所有基准站观测数据，选择、计算和播发用户的综合误差改正[5]信息，这就使得网内和网外的部分地区精度均匀分布。综上所述，考虑到CORS系统的精度及稳定性，基准站间最大距离应为30km～50km，其精度覆盖范围为网内及网外30km，且均匀分布。

4.2　CORS实时更新坐标转换参数的基本思路

CORS实时获取坐标转换参数需要2个前提条件：①数据处理中心与流动站实时双向数据通信;②CORS基准站有一定密度。双向数据通信是虚拟参考站(Virtual Reference Station，VRS)技术要求，所以，CORS系统满足第1个条件，而部分地区满足第2个条件。

基本思路如下，流动站将概略坐标发送给数据处理中心，数据处理中心根据其概略坐标判断出该点究竟位于由哪3个基准站组成的三角形内，并求出流动站至这3个基准站的距离，若流动站处在边缘地区，不在任何基准站三角形内，则选取距离流动站最近的3个基准站组成覆盖网络，并由基准站的WGS84坐标与1980年国家大地坐标或CGCS2000国家大地坐标反算出一套坐标转换参数，将其与虚拟参考站信息和差分改正数一同发送给流动站，用户只需套用坐标转换参数，完成相应历元的测量即可获得平面直角坐标;测量下一点时，数据处理中心判断流动站是否在三角网覆盖的一定范围(R覆盖)内，如果在覆盖范围内，不用发送新的坐标转换参数，否则重新计算一套新的坐标转换参数供用户使用。

五、结束语

本文结合CORS系统双向数据通信的技术要求和当前基准站网络的建设情况，研究了CORS基准站网络的组网形式和精度分布情况，对基于CORS实时获取高精度区域坐标转换进行了可行性探究。通过实时获取坐标转换参数实现区域坐标转换的方法能够消除服务缝隙，最佳拟合测区的覆盖范围，从而得到高精度的转换坐标，多快好省的服务与测绘行业。

参考文献

[1]　阳力，吕超.CORS系统卫星跟踪基准站建设的探讨[J].地理空间信息，2007，5(1)：30-31.

[2]黄俊华，陈文森.连续运行卫星定位综合服务系统建设与应用[M].北京：科学出版社，2009.

[3]梅胜强，宫煦利，姚宜斌，王泽民，黄石磊.基于Delaunay三角网的大规模CORS参考站组网技术研究[J].大地测量与地球动力学，2008，28(1)：131-135.

[4]　谷守周，等.基于CORS网的卫星钟差完备性监测[J].测绘科学，2013，38(3)：26-29.

[5]　高星伟.GPS/GLONASS网络RTK 算法研究及程序实现[D].武汉：武汉大学，2002.