学术交流丨魏子卿院士:北斗坐标系

来源:微信公众号中国测绘学会

引文格式:魏子卿, 吴富梅, 刘光明. 北斗坐标系[J]. 测绘学报,2019,48(7):805-809. DOI: 10.11947/j.AGCS.2019.20190131

第一作者:魏子卿(中国工程院院士,研究方向为大地坐标系、大地边值问题及GNSS测量)

编辑:张永超,

初审:齐阳,

审核:彭震中

摘要:

北斗坐标系是北斗卫星导航系统的大地基准。本文在扼要叙述背景情况之后,首先给出了北斗坐标系原点、尺度、定向的定义和参考椭球的定义和导出常数;其次介绍北斗坐标系的首次实现,这次实现由北斗系统的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。

表1 BDCS参考椭球的定义常数

表2 BDCS参考椭球的导出几何常数值

表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。

表 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) 注意北斗坐标系的标准化工作。

未来智能实验室的主要工作包括:建立AI智能系统智商评测体系,开展世界人工智能智商评测;开展互联网(城市)云脑研究计划,构建互联网(城市)云脑技术和企业图谱,为提升企业,行业与城市的智能水平服务。

  如果您对实验室的研究感兴趣,欢迎加入未来智能实验室线上平台。扫描以下二维码或点击本文左下角“阅读原文”

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/484860.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【原创】论码农的财富修养

《喜剧之王》里看到了演员的自我修养,然后想聊下码农的财富修养,这里不是聊理财,更多是保险Insurance;保险是什么,是一种通过支付少量的保险费用来应对未知风险的机制,比如车险、财险、健康险等&#xff1b…

centos 安装MySQL报错You could try using --skip-broken to work around the problem You could try running:

centos 安装MySQL报错You could try using --skip-broken to work around the problem You could try running: rpm -Va --nofiles --nodigest wget http://repo.mysql.com/mysql-community-release-el7-5.noarch.rpm rpm -ivh mysql-community-release-el7-5.noarch.rpm yum …

Java中序列化的好处及意义

1、序列化是干什么的?简单说就是为了保存在内存中的各种对象的状态,并且可以把保存的对象状态再读出来。虽然你可以用你自己的各种各样的方法来保存Object States, 但是Java给你提供一种应该比你自己好的保存对象状态的机制,那就是序列化。2…

​马斯克的全球上网计划:月费646元,网速可达200Mbps

SpaceX创始人兼首席执行官埃隆马斯克来源:腾讯科技审校作者:金鹿腾讯科技讯 4月16日消息,美国太空探索技术公司SpaceX在六个月前向公众推出了Starlink Early Access计划。在最初的几个月里,这项卫星互联网服务的用户增长到了1万多…

剑指 Offer 68 - II. (二叉树)二叉树的最近公共祖先

题目描述 给定一个二叉树, 找到该树中两个指定节点的最近公共祖先。 百度百科中最近公共祖先的定义为:“对于有根树 T 的两个结点 p、q,最近公共祖先表示为一个结点 x,满足 x 是 p、q 的祖先且 x 的深度尽可能大(一个节点也可以…

从狼群智能到无人机群协同决策

来源:飞思实验室文:小翼生物群集行为是自然界存在的一种普遍现象 ,先前的研究已表明人们可以通过模仿鸟群、鱼群、蜂群等行为, 研究其群体性决策机制, 进而用于构建集群系统的决策系统框架. 狼群在认知、分工, 以及复杂环境适应性等方面具有突出的优势, …

解决问题 com.alibaba.fastjson.JSONObject cannot be cast to xxx

{"msg": "success","code": "success","data": [{"create_time": "2020-11-09 11:54:08","exec_status": "正常","host_id": 1527414,"id": 3136268,"ip…

CSS常用样式

边框 http://www.w3school.com.cn/cssref/pr_border.asp p {border:5px solid red; } 可以按顺序设置如下属性: border-width border-style border-colorobject.style.border"3px solid blue"边框阴影(CSS3 ) http://www.w3school.com.cn/cssref/pr_box-s…

华为自动驾驶实车实路测试视频曝光!

来源:量子位、物联网智库等自动驾驶的风又刮起来了,这次处在风暴中心的是华为。4月15日,消息称配备华为自动驾驶技术的北汽新能源极狐阿尔法S的HI版车型在上海进行了公开试乘,这也是华为自动驾驶技术公开试乘的全球首秀。就试乘效…

Leetcode--450. 删除二叉搜索树中的节点

给定一个二叉搜索树的根节点 root 和一个值 key,删除二叉搜索树中的 key 对应的节点,并保证二叉搜索树的性质不变。返回二叉搜索树(有可能被更新)的根节点的引用。 一般来说,删除节点可分为两个步骤: 首先…

未来智能:人有人用,机有机用

来源: 苇草智酷作者:王飞跃当今最好的人工智能程序到底有多智能?它们是如何工作的?能做些什么?我们有必要担心机器比人类聪明且将很快夺取我们的工作吗?以上这些问题都颇受人们的关注,梅拉妮米歇…

Leetcode--329. 矩阵中的最长递增路径

给定一个整数矩阵,找出最长递增路径的长度。 对于每个单元格,你可以往上,下,左,右四个方向移动。 你不能在对角线方向上移动或移动到边界外(即不允许环绕)。 示例 1: 输入: nums [ [9,9,…

排序及查找----[(冒泡,快速)(拉格朗日,二分)]

代码展示: 1 using System;2 using System.Collections.Generic;3 using System.Linq;4 using System.Text;5 using System.Threading.Tasks;6 7 namespace paixu8 {9 class Program10 {11 static void Main()12 {13 int[] Array { 1, 45…

因果作用推断、因果网络学习及其他

来源:久章智能 未来智能实验室的主要工作包括:建立AI智能系统智商评测体系,开展世界人工智能智商评测;开展互联网(城市)云脑研究计划,构建互联网(城市)云脑技术和企业图谱…

学习《图解HTTP》有感

HTTP协议的出现主要是为了解决文本传输的难题(大家互相共享文件没有一个规范),最初的设计理念是:借助多文档之间的关联形成超文本,成为相互参阅的WWW(万维网) 计算机和网络设备要通信&#xff…

mysql解决select * from 表名 (where + 约束条件为空)

mysql解决select * from 表名 (where 约束条件为空),示例如下: SELECT * from tableName WHERE name is NULL; 从 tableName 表中查询 name 字段为空的记录。转载于:https://www.cnblogs.com/zrmw/p/9887655.html

通信人工智能的下一个十年(上)

来源:电信科学欧阳晔1, 王立磊1, 杨爱东1, 马利克萨哈2, 大卫贝兰格3,4, 高同庆5, 韦乐平6, 张亚勤71 亚信科技(中国)有限公司,北京 1001932 美国威瑞森电信公司,美国 纽约 100363 美国电话电报公司,美国 得…

[Swift]LeetCode19. 删除链表的倒数第N个节点 | Remove Nth Node From End of List

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★➤微信公众号:山青咏芝(shanqingyongzhi)➤博客园地址:山青咏芝(https://www.cnblogs.com/strengthen/)➤GitHub地址&a…

CAS

cas:Compare and Swap,即比较再交换。 对CAS的理解,CAS是一种无锁算法,CAS有3个操作数,内存值V,旧的预期值A,要修改的新值B。当且仅当预期值A和内存值V相同时,将内存值V修改为B&…

大脑使用交叉存储,来区分现在和过去

来源:混沌巡洋舰在每一个清醒的时刻,人类和其他动物必须在我们对过去的记忆和对现在的意识两者中保持平衡。我们必须吸收关于周围世界的新的感官信息,同时保持对早期观察或事件的短期记忆。我们理解周围环境、学习、行动和思考的能力都依赖于…