摘要:北斗坐标系是北斗卫星导航系统的大地基准。本文在扼要叙述背景情况之后,首先给出了北斗坐标系原点、尺度、定向的定义和参考椭球的定义和导出常数;其次介绍北斗坐标系的首次实现,这次实现由北斗系统的8个监测站在历元2010.0的ITRF2014框架下的坐标和速度体现,它们是通过处理一个包括8个监测站的全球GNSS网4期GPS数据、框架对准与监测站坐标序列的线性回归拟合得到的;最后就参考框架的更新周期、监测站境外布设和连续观测、精密星历的生成与发布,北斗坐标系的标准化问题进行了讨论,并提出一些建议。关键词:北斗坐标系;北斗卫星导航系统;参考框架;监测站一个全球卫星导航定位系统的基本任务是为全球用户提供可用性、完好性、连续性和精度符合规定要求的导航定位服务。卫星导航系统赖以导航定位的大地基准是其使用的坐标系。坐标系对导航定位系统的重要性不言而喻。坐标系在很大程度上决定了导航系统的性能,特别是导航定位的精度。另外,坐标系也是决定一个导航系统的国际兼容性和系统互操作性的重要因素。毫无疑问,一个全球卫星导航系统的大地基准应是一个地心的全球坐标系。根据2013年12月发布的北斗卫星导航系统空间信号接口控制文件,北斗导航系统的大地基准采用2000中国大地坐标系(CGCS2000)[1]。CGCS2000是我国的国家大地坐标系,2008年开始正式实施。在北斗系统建设的早期阶段,笔者利用2007—2009年的监测站GPS观测数据,计算了监测站的坐标和速度,建立了北斗系统的试用坐标系。实际上,这一试用坐标系,在严格意义上并不是CGCS2000,实际上更接近ITRF2000[2-3]。这一试用坐标系,在北斗系统的前期建设中曾经发挥了积极作用。但是,随着时间的推移,我们发现,将北斗系统的大地基准与国家大地坐标系捆绑一起是不适当的,不切实际的,对北斗系统的长远发展也是不利的。从长远观点来看,一个卫星导航系统,应该采用属于自己的独立大地基准[2]。GPS、GLONASS及Galileo均是如此[4-5]。北斗导航系统采用独立于国家大地坐标系,采用自己的专用坐标系,当属正确决策。于是,就提出了“北斗坐标系”。2017年12月,中国卫星导航系统管理办公室发布的《北斗卫星导航系统空间信号接口控制文件公开服务信号B2a(1.0版)》正式规定,北斗导航系统采用专用坐标系——北斗坐标系[6]。本文简要介绍北斗坐标系的定义和首次实现,以及坐标精度评估,并就某些问题提出一些建议。定义1名称与标识北斗坐标系是北斗卫星导航系统的大地基准。英文名称是BeiDou coordinate system,缩写BDCS。北斗坐标系通过参考历元的地面监测站坐标和速度实现,称为参考架。北斗坐标系将通过重新实现使参考框架最现时化和精度最佳化。坐标系的每次实现,对应产生一个新的参考框架。随着时间的推移,北斗坐标系将出现多个参考架。不同参考架的标识是BDCS(W×××),括号内符号W×××标示该参考架开始执行的北斗系统时(BDT)第×××周的0秒。例如,BDCS(W465)、BDCS(W1002)分别标示从BDT时第465周0秒开始执行的参考框架和从BDT时第1002周0秒开始执行的参考框架。2坐标系原点、尺度与定向北斗坐标系的原点、尺度与定向的定义如下。原点:包括海洋和大气的整个地球的质量中心。尺度:长度单位是m(SI)。这一尺度同地心局部框架的TCG时间坐标一致。定向:在1984.0时初始定向与BIH的定向一致。定向时间演变:定向随时间的演变使得整个地球的水平构造运动无整体旋转。上述定义与IERS规范[7]一致。北斗坐标系为一右手直角坐标系。原点为地球质量中心,Z轴指向IERS参考极方向,X轴为IERS参考子午面与通过原点且同Z轴正交的赤道面的交线,Y轴完成右手直角坐标系,如图 1。图1 BDCS参考椭球3参考椭球北斗坐标系采用CGCS2000参考椭球,即CGCS2000坐标系的参考椭球,该参考椭球的定义常数和导出常数见表 1—表 3。常数数值长半轴/ma=6 378 137.0扁率f=1:298.257 222 101地心引力常数(包括大气层)/m3s-2GM=3 986 004.418×108地球自转角速度/rad s-1ω=7 292 115.0×10-11表1 BDCS参考椭球的定义常数常数数值短半轴b/m6 356 752.314 1线偏心率E/m521 854.009 700 25极曲率半径c/m6 399 593.625 9第一偏心率平方e20.006 694 380 022 90第一偏心率e0.081 819 191 042 816第二偏心率平方e′20.006 739 496 775 48第二偏心率e′0.082 094 438 151 917轴比b/a0.996 647 189 319子午圈一象限弧长Q/m10 001 965.729 3椭球体积V/km31 083 207 319 783.546椭球表面积S/km2510 065 621.718平均半径R1/m6 371 008.771 4同面积之球的半径R2/m6 371 007.180 9同体积之球的半径R3/m6 371 000.790 0表2 BDCS参考椭球的导出几何常数值常数数值椭球面正常位U0/m2s-262 636 851.714 92阶带谐系数J20.108 262 983 225 8×10-24阶带谐系数J4-0.237 091 125 614 1×10-56阶带谐系数J60.608 346 525 889 2×10-88阶带谐系数J8-0.142 681 100 979 8×10-1010阶带谐系数J100.121 439 338 334 3×10-13赤道正常重力γe/ms-29.780 325 336 1极正常重力γp/ms-29.832 184 937 9平均正常重力γ/ms-29.797 643 222 4m=ω2a2b/(GM)0.003 449 786 506 78k=bγp/aγe-10.001 931 852 619 31地球质量M(包括大气)/kg5.973 331 96×1024表3 BDCS参考椭球的导出物理常数值CGCS2000参考椭球的几何中心与坐标系的原点重合,旋转轴与坐标系的Z轴一致。CGCS2000参考椭球为一等位旋转椭球,参考椭球面既是大地经纬度、高程的几何参考面,又是地球外部正常重力场的参考面。首次实现北斗坐标系通过地面监测站在参考历元的坐标和速度实现。因此如何获得监测站的坐标和速度,以及如何进行框架对准是北斗坐标系实现的关键。在本次实现中,监测站的坐标通过GNSS数据处理得到,参考框架对准通过最小约束法实现。1监测站GPS观测因早期北斗系统不具备全球观测能力,为了实现并维持其坐标系,截至2016年,在北斗系统地面监测站进行了4期GPS观测。第1期为2007—2009年,各站观测依次单独进行。第2期从2011年12月16日至12月31日,连续观测15 d,8个监测站同步联测。第3期从2014年4月24日至5月8日,连续观测15 d,8个监测站同步联测。第4期于2016年5月至11月单独观测完成。2北斗坐标系实现实现北斗坐标系的数据处理分3个步骤。第1步是联合处理由全球IGS站、国内陆态网基准站和北斗地面监测站组成的全球GNSS网数据得到监测站坐标的单日松弛解;第2步采用最小约束法将北斗监测站坐标单日解体现的框架,对准ITRF2014框架;第3步是将经框架对准的监测站坐标序列进行线性回归拟合,得到历元2010.0的监测站坐标和速度,最终实现北斗坐标系。2.1 站坐标松弛解在全球选取均匀分布且稳定可靠的62个IGS站,国内选取27个陆态网基准站,与8个北斗地面监测站组成全球GNSS网。全球网GPS数据处理用GAMIT10.5软件。软件采用绝对天线相位中心模型、FES2004海潮模型。为了更好地保持网形结构同时不影响松弛解解算结果,国际IGS站坐标给予1 m约束,国内基准站坐标给予10 m约束。除其他参数,输出参数包括站坐标单日松弛解及其协方差矩阵。2.2 框架对准监测站坐标单日松弛解是有先验松约束的站坐标解。在框架对准之前,需要将先验松约束去除。为此,执行去约束运算[8-11] (1)式中,∑unc代表去约束后的坐标协方差阵;∑est代表坐标松弛解的协方差阵;∑const代表加入的先验约束值。然后,采用如下方程加入最小约束[10-12],将去约束解体现的框架对准于由62个IGS站ITRF2014(http://itrf.ensg.ign.fr/ITRF_solutions/2014/ITRF2014.php)坐标定义的参考基准 (2)式中,N=(∑unc)-1;ΔX=XC-Xapr, Xapr是站坐标近似值(松弛解采用的概略坐标),XC是松弛解的站坐标;K是松弛解法方程自由项;B=(ATA)-1AT,A是由概略坐标构成的偏导数设计矩阵 (3)根据式(2),得到北斗监测站4期观测时段的ITRF2014坐标序列。2.3 参考历元坐标和速度计算用线性回归方程拟合监测站4期站坐标序列,未知参数取坐标和速度,协方差矩阵由坐标时间序列的方差确定,从而得到历元2010.0的监测站ITRF2014框架下的坐标和速度。表 4列出得到的监测站速度及其标准差。此外,为了验证结果的可靠性,将此次计算结果与前3期观测数据的计算结果进行比较,发现北斗监测站坐标差异不到1 mm。m/a站名VXVYVZσVXσVYσVZ北京-0.010 93-0.040 10-0.038 600.000 70.000 50.000 5成都-0.033 07-0.003 87-0.007 630.000 20.000 20.000 3哈尔滨-0.029 93-0.007 43-0.007 630.000 40.000 40.000 4喀什-0.032 20-0.004 530.011 570.000 30.000 20.000 2汕头-0.032 77-0.007 53-0.010 330.000 30.000 20.000 4三亚-0.030 37-0.001 43-0.009 070.000 50.000 30.000 6乌鲁木齐-0.031 730.000 170.004 630.000 30.000 30.000 3拉萨-0.046 00-0.008 600.012 100.000 20.000 10.000 2表 4 北斗监测站速度及标准差讨论与建议(1) 关于坐标系重新实现周期。一个坐标系的每次实现无论多么精确,随着时间的推移,监测站坐标和速度的误差积累将越来越大,以致坐标系重新实现成为不可避免。坐标系重新实现的结果就是框架的更新。重新实现就是利用历时更长的观测数据重新确定监测站坐标,并将框架对准于最新的ITRF框架。现在的问题是,参考框架更新周期多长比较合适?一般而论,总的原则是,以北斗框架与ITRF框架的失准度小于规定的限差为准,或者简单的,北斗框架更新与ITRF框架更新的周期相同,更新时间可以稍有落后。当然,也可以根据其他实际情况而定,例如由于监测站的数量变动或升级改造,框架需要更新。笔者建议,一般情况下,北斗框架的更新周期可以采用ITRF的更新周期但稍有滞后。(2) 关于监测站境外布设与连续观测。鉴于北斗系统的全球性质,监测站的境外扩展是必要的,也是必然的,这无疑有利于改善坐标系的实现精度。另外,假若监测站采用长期连续观测,对改善框架精度和现势性同样有益处。如果拥有长期连续观测数据,数据处理方法当然需要相应改变,或许可以采用监测站运动非线性假设,有可能得到更逼真的监测站速度。因此,笔者建议,监测站的境外布设以及长期连续观测,应是今后努力方向之一。(3) 关于卫星精密星历。如前所述导航系统的坐标系通过参考历元的监测站坐标和速度实现。这是实现坐标系的最基本形式,因为监测站坐标是卫星轨道确定必不可少的输入参数。作为卫星定轨的输出,卫星星历也是坐标系的重要体现(或实现)。卫星星历分为精密星历和广播星历。精密星历的精度为厘米级,即代表厘米级精度的参考框架。广播星历实际上是精密轨道的外推,精度一般在米级水平,甚至更低一些,代表低精度的参考框架。广播星历是用户最容易得到的坐标系产品。精密星历甚至也可以从网站下载得到。为了北斗系统应用的推广,对于精密定位(点定位和相对定位),有必要组织例行化的北斗卫星精密星历的生成,并通过有效形式向用户发布。(4) 关于北斗坐标系的标准化。卫星星历(轨道)既然是重要的坐标系产品,那么卫星定轨所用的卫星动力学模型、地球重力场模型、天球参考系与地球参考系之间的转换模型和地球定向参数,以及其他一些重要常数,当属坐标系的重要内容。在实际工作中,应将这些元素纳入坐标系一起考虑。为了使坐标系产品生产的规范化,与北斗坐标系工作的正规化,建议在条件成熟时,有必要制定北斗坐标系标准,以便实际工作有所遵循。总结(1) 北斗系统采用专用的地心大地基准。(2) 北斗系统的大地基准中文名称是北斗坐标系,英文名称是BeiDou coordinate system,英文缩写BDCS。(3) 北斗坐标系的不同框架的标识是“BDCS(W×××)”,其中W×××标示该框架自北斗系统时第×××周0秒开始执行。(4) 北斗坐标系的原点、尺度、定向与IERS规范规定一致,参考椭球采用CGCS2000椭球。(5) 北斗坐标系通过监测站坐标和速度与卫星星历实现。北斗坐标系还包括参考椭球、重力场模型、地球定向参数以及坐标转换参数等。(6) 北斗坐标系的首次实现,包括参考历元2010.0时8个监测站ITRF2014框架下坐标和速度,通过8个监测站4期GNSS数据与62个国际IGS站和27个国内陆态网络基准站数据的联合解算得到监测站坐标松弛解、框架对准与监测站坐标序列线性回归拟合得到。坐标标准差小于2 mm,速度标准差小于1 mm/a。关于北斗坐标系工作,提出如下建议:(1) 北斗坐标系的重新实现周期一般与ITRF相同。(2) 北斗系统采取全球跟踪与连续观测。(3) 组织北斗卫星精密星历的生成与发布。(4) 注意北斗坐标系的标准化工作。致谢: 感谢相关部门提供北斗观测数据。END引文格式:魏子卿, 吴富梅, 刘光明. 北斗坐标系[J]. 测绘学报,2019,48(7):805-809. DOI: 10.11947/j.AGCS.2019.20190131第一作者简介:魏子卿,中国工程院院士,研究方向为大地坐标系、大地边值问题及GNSS测量编辑:张永超初审:齐 阳审核:彭震中-----END-----社群交流 / 原创投稿 / 商务合作(请添加下方小家微信:jiace03)来源:中国测绘学会推荐阅读【注意】测绘间谍如何识别?发现后该怎么办?【技术】从地方坐标系到2000国家大地坐标系的快速转换方法【科普】实景三维模型和地形三维模型【公告】关于2021年中国测绘学会团体标准(第一批)立项的公告【行业】河南省测绘地理信息技术中心挂牌成立推荐关注温馨提示:近期,微信公众号信息流改版。每个用户可以设置 常读订阅号,这些订阅号将以大卡片的形式展示。因此,如果不想错过“测绘之家”的文章,你一定要进行以下操作:进入“测绘之家”公众号 → 点击右上角的 ··· 菜单 → 选择「设为星标」↓↓↓点击下方“阅读原文”查看更多精彩内容...
南方cass教你告别熬夜画图!
白天外业测绘,晚上还要熬夜加班画图,这是很多测绘小伙伴的工作常态。 如图1这种地形比较简单,外业采集了400+个碎步点,内业绘图需要多少时间?图1 新手需要2个小时以上,而高手一般50分钟即可完成,这...
08月31日评论
南方cass常见问题解决办法汇总(下)
南方cass如何设置多图显示? 南方cass粘贴到原坐标不能用,怎么办? 南方cass顶部菜单栏关闭了,怎么打开? 双击CAD图标,不想启动南方cass怎么办? 南方cass加载质检模块 提示系统错误...
05月24日评论
【图文】南方cass等高线绘制教程
在地形图中,等高线是表示地貌起伏的一种重要手段。常规的平板测图,等高线是由手工描绘的,等高线可以描绘得比较圆滑但精度稍低。在数字化自动成图系统中,等高线是由计算机自动勾绘,生成的等高线精度相当高。 ...
04月15日评论
【测量软件】横断面测量数据批量转换成断面图【CASS】
完成横断面外业测量,断面线上的高程点都是按顺序采集的(下图所示),想批量输出横断面图。如果一条条横断面手工连点绘制断面线,然后绘制断面图…… 横断面高程点截图 1、编辑横断面简码坐标文件 ...
02月23日评论
南方cass地形地籍成图软件,测绘地形图快速入门
本文根据数字化测图的特点介绍用地籍成图软件南方cass测绘地形图的方法,内容包括测区首级控制、图根控制、测区分幅、碎部测量、人员安排等,最后利用南方cass绘制一幅地形图。 01 控制测量和碎部测量原...
02月04日评论
南方cass地形地籍成图软件如何进行公路曲线设计?
南方cass基于AutoCAD和国产CAD平台,集数据采集、编辑、成图、质检等功能于一体的成图软件,广泛应用于地形成图、工程测量应用和数据空间建库等领域,并专注于满足国土、房产、测绘、规划、水利、交通...
01月25日评论
南方cass地形地籍成图软件如何进行区域土方量平衡?
土方平衡的功能常在场地平整时使用。当一个场地的土方平衡时,挖掉的土石方刚好等于填方量。以填挖方边界线为界,从较高处挖得的土石方直接填到区域内较低的地方,就可完成场地平整。这样可以大幅度减少运输费用。此...
12月28日评论
南方cass等高线法土方计算教程
如果将白纸图扫描矢量化后可以得到图形,但这样的图没有高程数据文件,所以无法用前面的几种方法计算土方量。 一般来说,这些图上都会有等高线,所以,南方cass开发了由等高线计算土方量的功能,专为这类用户设...
12月22日评论
南方cass软件如何用方格网法计算土方?
由方格网来计算土方量是根据实地测定的地面点坐标(X,Y,Z)和设计高程,通过生成方格网来计算每一个方格内的填挖方量,最后累计得到指定范围内填方和挖方的土方量,并绘出填挖方分界线。 系统首先将方格的四个...
12月14日评论
高手必备的CASS快捷键和命令大全
CASS高手是建立在一种左手键盘,右手鼠标的工作模式,掌握绘图的基本命令和快捷命令,日积月累,循序渐进的提高。现将最常用、最实用的命令汇总如下。 CASS9.1常用快捷命令 CASS9.1系统 Aut...
05月09日评论