基于GNSS基准站数据的控制点CGCS2000坐标计算.pdf

2023年11月引用格式：张吉君,刘彦杰，包暑光.基于GNSS基准站数据的控制点CGCS2000坐标计算J.水利水电快报,2 0 2 3,44（11）：7 3-78.水利水电快报EWRHI第 44 卷第 11 期基于 GNSS 基准站数据的控制点 CGCS2000坐标计算张吉君，刘彦杰，包暑光（长江三峡勘测研究院有限公司（武汉），湖北武汉430 0 7 4）摘要：针对中国部分区域CGCS2000控制点坐标较难获取的问题，提出了基于公开的GNSS数据计算控制点CGCS200O坐标的方法。选用GNSS基准站进行GNSS静态测量联测，利用GAMIT/GLOBK软件进行数据处理，获取未知点当前框架和当前历元下的空间三维坐标以进行框架历元转换，并计算未知控制点的CGCS2000坐标。结果表明：采用该方法计算控制点的平面精度总体保持在3cm以内，高程精度总体保持在5cm以内，可用于项目前期快速获取控制点 CGCS2000坐标和项目中控制点CCCS2000 坐标的校核。关键词：CGCS2000坐标；GAMIT/GLOBK；框架转换；历元转换；GNSS基准站中图法分类号：P226文章编号：10 0 6-0 0 8 1(2 0 2 3)11-0 0 7 3-0 6文献标志码：AD0I:10.15974/ki.slsdkb.2023.11.0120 引 言CGCS2000是一个严密的坐标系统,是ITRF97框架、2 0 0 0.0 历元下的瞬时坐标，也是中国目前正在推广并使用的新一代大地坐标系,所有测绘成果均需采用或者转换到CGCS2000坐标。但是，普通用户很难获得具有CGCS2000坐标成果的控制点，需要向测绘成果主管部门申请或者把原始观测数据送到有关部门解算才能获得；由于测绘成果保密的原因,这两种方式既耗时又繁2 。目前,对于获取控制点CGCS2000 坐标的研究,尹业彪2 通过GAMIT软件对IGS站观测数据进行解算,研究了通过现有资料和技术手段,实现获取CCCS2000坐标的可行性,提高小型测绘单位CGCS2000坐标的自给率,降低小型测绘活动的生产成本、提高生产效益。胡思华和曾庆春等在工程项目中联测IGS站获取CGCS2000坐标,并分析其精度与实用性3-4。黎鹏等5 将待测点与 ICS站组网解算得到高精度坐标，然后通过动态转换获取CGCS2000坐标,并探讨其在实际工程应用中的可行性。刘明波等6 采用联测IGS 站方法,为计算偏远地区控制点的CGCS2000坐标提供实用思路。收稿日期:2 0 2 3-0 3-14作者简：张吉君，男，工程师，硕士，主要从事工程勘测、大地测量与摄影测量方面的工作。Ema i l:z h a n g j i j u n c j w s j y.c o 通信作者：刘彦杰，男，高级工程师，主要从事水利水电勘测、地理信息系统方面的管理和研究工作。Ema i l:l i u y a n j i e .73但是,以上研究都只对某一特定区域GNSS数据进行解算和分析，具有一定的局限性。本文提出一种选用公开的GNSS基准站并对其坐标成果进行验证的方法，通过将分布于5个不同省份的控制点联测GNSS基准站，根据CGCS2000定义,结合坐标参考框架转换和历元归算，将对应ITRF框架当前历元下的坐标转换为ITRF97框架（参考历元为2000.0）下的坐标,即CGCS2000坐标。GAMIT/GLOBK是受到国内外学者一致认可的高精度 GNSS数据处理软件之一，处理基线和连续观测的静态数据时,解算时间快，能够有效地消除或者削弱GNSS测量误差带来的影响,获得高精度基线解算结果7-10 1。基线结果在平差时不受测区已知控制点的约束,采用GNSS基准站作为已知控制点，解算过程中使用卫星精密星历、对流层等改正模型,获取测站高精度的三维坐标信息，有效解决了常规GNSS测量依赖已知控制点、短基线解算的缺点 。因此，利用CAMIT/CLOBK软件进行数据处理和平差，可方便快捷地直接获取控制点CGCS2000坐标。经验证，该方法可直接获取大陆地区任意区域的控制点CGCS2000坐标,解决工程项目2023年11月呕需的已知点坐标问题,具有快捷性和高可靠性。随着中国深入推进“双碳”目标，加快能源结构调整，新疆、西藏等具备丰富新能源资源的偏远地区正在加紧推进抽水蓄能电站等大批水利水电项目。这些项目所在地通常远离城镇，没有通讯信号。因国家高等级控制点分布不均匀且项目前期申请控制点需要一定周期，对于急需完成各项报批流程的项目，如何尽快完成前期地形图测绘、地质钻孔放样是关注的重点。通过公开的 GNSS数据可以及时获取无网络信号区域控制点的CGCS2000系坐标，开展各项测绘工作，完成项目申报流程。1控制点计算方法公开的GNSS数据包括全球IGS站观测文件、香港 CORS（C o n t i n u o u s l y O p e r a t i n g Re f e r e n c e St a-tions)网观测文件以及其在各参考框架下的三维坐标和速度信息，还包括广播星历、精密星历、精密钟差等GNSS解算相关文件。利用公开的GNSS数据，将控制点联测公开的GNSS基准站，采用GAMIT软件进行基线处理；目前GAMIT求解的基线结果是ITRF2014框架、当前历元下的成果，从ITRF网站(https:/itrf.ign.fr/）获取 GNSS 基准站 ITRF2014（当前历元）下的精确坐标作为起算点，利用GLOBK软件进行空间三维平差，获取未知控制点ITRF2014（当前历元）下空间三维坐标，然后进行框架历元转换，由ITRF14框架（当前历元）下坐标转换到CGCS2000（IT RF97 框架、2 0 0 0.0 历元）坐标成果。图1为联测GNSS基准站计算控制点CGCS2000坐标的流程。框架转换参数可从ITRF网站上获取，其精度有保障；历元转换必须要获得观测控制点的速度值，其数值是利用控制点ITRF2014框架当前历元下的三维空间坐标,基于全国2 30 0 多个控制点多年实测速度拟合而成,或者利用 CPMC G C S2 0 0 0板块欧拉矢量计算控制点速度值12 2选用公开的GNSS基准站由于中国各省CORS站数据不公开,IGS站数据存在着分布数量偏少,且数据发布不连续的问题。因此，选用中国及其周边地区数据发布连续和稳定的IGS站点作为基准控制点，根据 CH/T2014-2016大地测量控制点坐标转换技术规范，基准控制点需要遵循连续性、稳定性、高精度、多种解、均衡性以及精度一致性原则。通过查阅InternationalGNSSService网站，核实IGS站点的数据信息和点位分布情况，如表1所示，最终选用41个IGS站点作.74水利水电快报EWRHI为常用的 GNSS 基准站。控制点GNSS观测数据ITRF2014框架（当前历元）空间三维坐标CGCS2000空间三维坐标高斯投影CGCS2000平面坐标与大地高程图1联测GNSS基准站计算控制点CGCS2000坐标的流程Fig.1 Process of calculating control point CGCS2000coordinate with GNSS reference station表 1常用的 GNSS 基准站Tab.1Common GNSS reference stations序号站名1AIRA2BJFS3CHAN4CKSV5CMUM6CPNM7DAEJ8DAE29GAMG10GUAM11HKCL12HKKT13HKNP14HKOH15HKQT16HKSC17HKSL18HKST19HKWS20HYDE21IISC为了检验所选IGS站点的稳定性，首先通过ITRF网站获取IGS站点对应框架和历元下的坐标，基准站ITRF2020框架（2 0 15.0 历元）和ITRF2014框架(2 0 10.0 历元）下的三维坐标及速度作为转换前输入成果，基准站ITRF2014框架（2 0 0 0.0 历元）和ITRF97框架（2 0 0 0.0 历元）下的坐标作为已知成果,然后通过以下3种方式检验IGS站点的稳定性和框架历元转换的正确性。转换结果见表2 3。（1）方式一：获取站点ITRF2020框架（2 0 15.0历元）下的三维坐标以及速度，将其框架历元转换为ITRF2014框架（2 0 0 0.0 历元）下的坐标。（2）方式二：获取站点ITRF2014框架（2 0 10.0历元）下的三维坐标以及速度，将其框架历元转换为ITRF2014框架（2 0 0 0.0 历元）下的坐标。第 44 卷第 11 期联测基准站观测数据精密星历等表文件GAMIT基线解算GLOBK平差计算框架历元转换序号2223242526272829303132333435363738394041站名JFNGKIT3KMNMLCK3LCK4LHAZPIMOPOL2SEJNSELESUWNTAEJTCMSTNMLTWTFULABURUMWUH2YAKTYONS张吉君等基于CNSS基准站数据的控制点CCCS2000坐标计算表2 方式一和方式二转换结果与已知成果比较差值Tab.2Difference between the conversion result of method 1 and method 2 and the known result方式一转换结果与已知成果比较差值序号站名1AIRA2BJFS3CHAN4GUAM5HYDE6IISC7KIT38LHAZ9PIMO10POL211SELE12SUWN13TAEJ14TCMS15TNML16ULAB17YAKT注：dX,dY,dz分别是框架历元转换后坐标与已知成果坐标在X,Y,Z方向的差值，（dP)2=(dX)+（d Y)+（d Z),下同。（3）方式三：获取站点ITRF2014框架（2 0 10.0历元)下的三维坐标以及速度，将其框架历元转换为ITRF97框架（2 0 0 0.0 历元）下的坐标。将框架历元转换后的坐标与已知成果进行比较,即可验证所选IGS站点的稳定性。从表2 和表3可知，通过框架历元转换获得的坐标精度为厘米级，方式一最大点位误差为3.7 5cm（站点BJFS），方式二比较结果最大点位误差为4.94 cm（站点PIMO）,方式三比较结果最大点位误差为2.7 6 cm（站点IISC），证明了此法获取CGCS2000坐标精度的可靠性。最后根据工程项目所在位置,在此41个IGS站点中选用合适的GNSS基准站作为联测点,用于计算未知控制点CGCS2000坐标。表3方式三转换结果与已知成果比较差值Tab.3Difference between the conversion result of method3 and the known result站名dXBJFS0.21GUAM1.54IISC0.91KIT30.25LHAS0.13POL20.66SELE0.04SHAO0.73SUWN0.46TAEJ0.593基线解算策略和精度评定方法3.1基基线解算策略利用GAMIT软件,采用IGS精密星历，使用批cm方式二转换结果与已知成果比较差值dXdY-2.360.18-2.08-2.490.370.171.050.34-1.310.18-1.171.050.950.111.320.582.500.630.16-0.660.000.630.291.510.06-0.02-0.161.110.25-0.33-0.630.301.16-0.38dYdz1.722.101.31-1.482.51-0.69-0.89-0.83-1.01-1.10-0.77-0.100.370.950.12-1.800.260.671.041.33dZ-1.23-1.881.29-1.06-1.730.130.280.191.091.020.070.620.170.210.071.05-2.31cmdP2.722.502.761.241.501.011.021.950.851.78dP2.673.751.361.532.171.581.001.452.801.230.641.660.181.140.421.262.61处理方式进行基线解算。在 Station.info文件中配置好各个控制点测站信息，包括控制点编号和点名、观测时间段、仪器类型、天线类型、天线高以及其量测方式等。在Process.defaults文件中设置好采样率、观测历元时长以及基线开始处理时间等。Site.defaults设置好rinex数据来源方式、参与平差的控制点和已知点。在Sestbl表中,观测值选项（Choiceof observable）设定为LC_AUTCLN；试验选项(（C h o i c e o f Ex p e r i me n t)设定为BASELINE。所有文件其余主要选项采用软件推荐的缺省值。3.2米精度评定方法标准化均方根误差NRMS值一般用来描述在单位时段内基线解算值与其加权平均值的偏离程度，是从历元的模糊度解算中得出的残差值，也是衡量解算质量的最主要指标之二7 。NRMS值越小,基线解算精度越高。一般认为NRMS值不大于0.3说明解算成功，若数值大于0.5说明解算过程中部分周跳可能未被探测修复，出现了某些参数数值偏差较大等异常,需检查原因,重新解算7 。4控制点CGCS2000坐标的计算与验证挑选5个区域的控制点进行静态测量，区域14控制点已知成果由权威机构解算获得，区域5已知成果为向测绘主管部门申请获取。5个区域的静态测量和GNSS基准站信息如表4所示。利用GAMIT对各个区域静态测量数据和联测GNSS基准站数据进行基线处理，采用NRMS值来.75.dX-0.49-0.010.150.88-0.76-0.740.00-1.07-0.190.000.29-0.370.00-0.110.31-0.011.11dY3.470.00-0.490.510.12-0.200.000.37-0.51-1.00-0.260.160.000.960.170.00-0.38dz0.430.00-0.41-1.44-1.200.010.010.124.910.05-0.250.40-0.01-0.40-0.350.00-2.22dP3.530.010.661.761.430.770.011.134.941.010.460.570.011.050.500.012.512023年11月Tab.4 Static measurement information and GNSS reference station information of five regions控制点观测时间区域编号1SYOI SY042CY01 CY052021-07-103LY01 LY052022-07-15 s4KL01 KL062021-10-28 5WZ01WZ032021-09-24衡量基线处理的质量好坏，如表5所示，GAMIT软件基线解算的NRMS值均为0.19 0.2 1,这说明各个区域GNSS网的整体外业观测质量较高,基线解的精度好。利用CLOBK软件对各个区域基线处理后的结果进行平差处理，平差后获得各个控制点ITRF2014框架（当前历元）下坐标，接着计算控制点速度值,最后进行框架历元转换，获取观测控制点CGCS2000空间三维坐标和高斯平面坐标。表5基线处理精度统计Tab.5Statistics of baseline processing accuracy区域年积日1106(2021 年)107(2021年)108(2021 年)111(2021 年)112(2021年)2191(2021年)3196(2022年)197(2022 年)4301(2021 年)302(2021年）5267(2021 年)将控制点计算的CGCS2000坐标与已知成果比较，统计各个控制点点位较差及其均值如图2、表6 所示。从比较成果可知：采用联测GNSS基准站方法求取控制点的CGCS2000坐标是可行的，平面精度效果较好，高程精度效果不佳；区域1 4空间三维点位较差dp基本都保持在5cm以内,点位较差均值最大的是区域2（4.16水利水电快报EWRHI表45个区域的静态测量信息和GNSS基准站信息观测时长/(年月日)h2021-04-16 72.02021 04 225.56.13.7 4.42022 07-168.1 8.32021-10 29约8LHAZ,PIMONRMS0.193050.201790.195350.196760.202780.193840.199200.200940.205410.205470.20474第 44 卷第 11 期平差固定地理位置联测基准站基准站BJFS,CPNM,HKSL,JFNG,LHAZ,BJFS,LHAZ,PIMO,POL2PIMO,POL2BJFS,CHAN,DAE2,GUAM,HKSL,CHAN,GUAM,JFNG,LCK3,PIMO,POL2,ULABLCK3,POL2AIRA,BJFS,CHAN,DAEJ,GUAM,BJFS,CHAN,JFNG,LHAZ,POL2DAEJ,JFNGBJFS,CHAN,CPNM,DAE2,HKSL,BJFS,JFNG,JFNG,LCK3,POL2,URUMLCK3,URUMBJFS,CMUM,HKSL,JFNG,KIT3,BJFS,CMUM,LHAZ,PIMOcm）,最小的是区域4（2.15cm）；区域1.4高斯平面点位较差ds基本保持在3cm以内，点位较差均值最大的是区域1（2.32 cm），最小的是区域3(0.47 cm）;高程精度均值dH保持在4cm以内;由于区域2 4的已知成果也是采用GAMIT/GLOBK计算所得，平面点位较差最大是1.88 cm,精度在2 cm以内；但是因为联测的GNSS基准站和GAMIT/GLOBK软件数据处理设置参数不一致，高程精度比较差，基本保持在5cm以内；区域1已知成果是使用其他数据处理软件计算出来的，平面坐标存在X方向差别较大的现象,平面点位精度和高程精度都在3cm以内；区域5已知成果是向测绘主管部门申请获得，由于已知控制点都在高层房屋楼顶，埋设年代久远，计算成果与已知成果比较存在一定偏差，点位精度不大于5.5cm。6dP均值日ds均值5F口dH绝对值均值4.1643.943.513.473.1832.322.16210区域1区域2图2 点位较差均值统计Fig.2 Mean value statistics of point calibration difference（区域中心位置）广东广州市和东莞市(1133242E,224945N)辽宁朝阳市(1200543E,413516N)山东临沂市(1180418E,353049N)甘肃康乐县(1034309E,350557N)广西梧州市(1110850E,232601N)5.234.232151.73.281200.47区域3区域区域4区域5.76.张吉君等基于GNSS基准站数据的控制点CGCS2000坐标计算表6 联测CNSS基准站计算CGCS2000坐标成果与已知成果的比较Tab.6Comparison of CGCS2000 coordinate results calculated by joint survey GNSS reference station with known results区域点名1SY01SY02SY03SY042CYO1CY02CY03CY04CY053LYO1LY02LY03LY04LY054KLO1KL02KL03KL04KL05KL065WZ01WZ02WZ03注：(ds)=(dx)+（d y）,dH表示大地高程之差。5结语联测公开的GNSS基准站可以方便获取已知点对应ITRF框架当前历元下的三维空间坐标，通过GAMIT/GLOBK软件基线处理和平差计算，可获得未知控制点对应ITRF框架当前历元下的三维空间坐标，再进行框架历元转换，即可获得未知点的CGCS2000坐标。比较不同省份的5个区域计算成果与已知成果，结果表明：平面精度效果较好,基本保持在3cm以内；高程精度效果略差，基本保持在5cm以内。该方法利用GAMIT/GLOBK软件计算出控制点的CGCS2000坐标具有快捷性和高可靠性，其精度满足项目申报前期需要快速获取控制点CGCS2000坐标的需求，可用于项目控制点CCCS2000坐标的校核,也可用于某些中小型工程或者测区周边没有已知高等级控制点的区域。参考文献：1 走赵忠海,蒋志楠,朱李忠.WGS-84（G 16 7 4）与 CGCS2000坐标转换研究J.测绘与空间地理信息,2 0 15,38(4):188-189,192.cm空间三维点位较差平面点位较差dXdY0.810.621.01-1.100.930.960.70-0.64-1.611.212.721.98-2.031.57-2.141.78-1.881.39-1.341.96-2.253.56-1.152.58-1.792.98-1.391.94-1.810.44-1.810.64-1.770.22-1.800.57-1.970.21-1.640.38一一一一dZ2.82-2.96-3.10-2.902.863.643.233.313.091.362.221.302.291.391.181.310.921.151.080.89一一一dP3.003.313.383.053.504.954.134.323.882.734.763.114.172.762.202.322.012.212.261.90一一2尹业彪.基于连续观测站的2 0 0 0 国家大地坐标获取方法研究D.西安：西安科技大学,2 0 14.3胡思华.控制点CGCS2000坐标的联测及解求J.矿山测量,2 0 2 1,49(1):30-34.4曾庆春，徐龙华.在工程项目中引入南昌2 0 0 0 坐标系统的探讨J.城市勘测,2 0 19（5）：146-148.5 黎鹏，何雯，罗晔.2 0 0 0 国家大地坐标动态转换方法探讨J.水利水电快报,2 0 2 0,41(7)：2 6-2 9.6刘明波，李祖锋，鹿恩锋.基于ICS站联测的控制点CGCS2000坐标求取J.工程勘察,2 0 2 1,49(3:47-49,67.7慕仁海，常春涛，党亚民，等.GAMIT10.71解算GNSS长基线精度分析J.全球定位系统，2 0 2 0,45（5）：14-19,83.8李建涛,朱兰艳,李永梅，等.基于GAMIT的不同参数对北斗长基线精度的影响分析J.全球定位系统，2018,43(5):23-28.9杨登科,安向东,汤.GAMIT数据处理中基线解算模式的对比分析J.测绘地理信息,2 0 16,41（2）：18-21.10 付强，王超,黎鹏.Trimble RTX-PP在长江中游河段静态数据解算中的精度分析J.水利水电快报，2022,43(7):62-65,77.77.dx2.23-2.16-2.352.330.940.700.840.720.830.290.620.640.12-0.210.390.400.310.360.460.67-3.113.144.75dy0.520.49-0.49-0.390.751.320.950.950.890.240.36-0.240.140.361.661.561.651.621.821.700.322.68-2.62ds2.292.212.402.361.201.501.271.191.220.370.710.680.180.421.711.611.681.661.881.823.134.135.43dH-1.912.44-2.36-1.923.294.713.924.153.652.714.713.054.162.731.371.621.041.461.180.53-5.35-5.17-5.172023年11月11杨德芳，丁开华，许伟.GAMIT/CLOBK在高精度长基线像控点解算中的应用J.西安科技大学学报，2019,39(5):849 858.Calculation of CGCS2000 coordinates of control pointsbased on GNSS reference station data水利水电快报EWRHI12 刘光明.CGCS2000坐标转换M.北京：测绘出版社,2 0 2 1.(编辑：江焘，高小雲）第44卷第11期ZHANG Jijun,LIU Yanjie,BAO Shuguang(Three Gorges Geotechnical Engineering Co.,Ltd.,(Wuhan),Wuhan 430074,China)Abstract:In view of the difficulty in obtaining the coordinates of CGCS2000 control points in some regions of Chi-na,a method of calculating the coordinates of CGCS2000 control points based on the published GNSS data was pro-posed.GNSS reference stations were selected for GNSS static measurement and joint survey;GAMIT/GLOBK soft-ware was used to process data and obtain the spatial three-dimensional coordinates of the current frame and currentepoch of the unknown point,so as to convert the frame epoch and calculate the CGCS2000 coordinates of the un-known control points.The results showed that the plane accuracy of the control points calculated by this method wasmostly kept within 3 cm and the elevation accuracy was mostly kept within 5 cm,which could be used to quickly ob-tain the CGCS2000 coordinates of the control point in the early stage of the project and check the CGCS2000 coordi-nates of the control points in the project.Key words:CGCS2000 coordinate;GAMIT/GLOBK;frame transformation;epoch transformation;GNSS refer-ence station欢迎关注“长江数字图书馆”微信公众号.78.