无人机搭载的RTK获得的经纬高坐标要转换为东北天坐标,才能用于局部的导航和定位。为了这个目的,查阅资料,越查越懵逼,竟然这么多的坐标系,略懂之后,将学到的信息记录如下,很多跟我的目的:”RTK的LBH坐标转换为东北天坐标“有点远,但是把周边知识搞的清楚了,还是有意义的。所以记录在下,最终的经纬高坐标系-ECEF坐标系-ENU坐标系的公式和原理在另一个博客里https://blog.csdn.net/j_____j/article/details/102636211有精简介绍。

RTK获得LBH坐标,LBH坐标系给出的坐标是(Ψ,λ,h)、(longitude,latitude,height),LBH坐标系(LBH坐标洗、经纬高坐标西、大地坐标系、地理坐标系),国际常用的是WGS-84坐标(坐标系的原点位于地球质心,z轴指向(国际时间局)BIH1984.0定义的协议地球极(CTP)方向,x轴指向BIH1984.0的零度子午面和CTP赤道的交点,y轴通过右手规则确定。),我国采用的是“2000国家大地坐标系”,经纬高坐标转换为ECEF坐标(ECEF坐标系、earth centered eath fixed、地心地固坐标系)(该坐标系以地球质心为原点,Z轴向北沿地球自转轴方向,X轴指向经纬度的(0,0)位置,右手系Y轴指向90度经线。该系与地球一同转动。因此地固坐标系是非惯性系),ECEF坐标系转换为站心地平直角坐标,站心地平直角坐标系(北东地坐标系、ENU坐标系)。

与ECEF坐标系比较相关的,有个ECI坐标系(ECI坐标系、earth centered inertial,地心惯性坐标系)(坐标原点取在地心,X轴指向春分点,Z轴指向北极,Y轴与前者构成右手系。该系不与地球一同转动,因此可以应用牛顿定律。是惯性系)

LBH坐标(WGS_84)→ ECEF坐标→ENU坐标系坐标转换公式,参考:https://wenku.baidu.com/view/e6e4b4c533d4b14e85246886.html,https://wenku.baidu.com/view/22239a0a763231126edb11ea.html

LBH坐标(WGS_84)→ ECEF坐标,

,其中,,因为扁率,所以也可以是

为经度(-180°-180°),为纬度(-90°-90°),N为寅卯圈半径,h为大地高,e为椭球偏心率(与地球长半径、短半径有关),a为椭球长半径,b为椭球短半径。在这个公式里,要提供给程序椭球长半径a和椭球短半径b。

WGS_84的参数:a=6378137.0,f=1/298.257223563。

ECEF坐标 → ENU坐标

其中(X0,Y0,Z0)为站心(第一个ECEF坐标点),LB为站心的经纬。

另外,ECEF坐标→ LBH坐标(WGS_84)

因为h的计算包含有待求解的,而的计算反过来又含有待求解的h,所以一般借助迭代法来逐次逼近。 迭代法一般计算过程如下:先假设的值等于0,分别计算出N,h和,然后再将刚得到的重新代入方程,再一次更新N,h和的值,如此循环。上述三式的迭代运算通常收敛的很快,一般经过3~4次的循环迭代就可结束运算。

LBH坐标(WGS84),RTK给的高度是大地高h(大地高)(地面点沿参考椭球面法线到参考椭球面的距离
),但是RTK设备可以计算输出海拔高。

此处要区分的是海拔高H(正高)(地面点沿重力线(铅垂线)到大地水准面的距离)(在中国,可以将似大地水准面的高程近似为海拔高)。

因此这里有个转换公式,RTK的高程数据(大地高)- 高程异常 = 海拔高(正高) 。

大地高和海拔高的区别在于,有椭球面和水准面,水准面是这样的(受地球自传、地球密度不均匀、引力等作用,所以是不规则的)

大地水准面 geoid
由静止海水面并向大陆延伸所形成的不规则的封闭曲面。它是重力等位面,即物体沿该面运动时,重力不做功(如水在这个面上是不会流动的)。大地水准面是指与全球平均海平面(或静止海水面)相重合的水准面。大地水准面是描述地球形状的一个重要物理参考面,也是海拔高程系统的起算面。大地水准面的确定是通过确定它与参考椭球面的间距-大地水准面差距(对于似大地水准面而言,则称为高程异常)来实现的。

区别于大地水准面,有个似大地水准面,从地面点沿正常重力线量取正常高所得端点构成的封闭曲面。似大地水准面严格说不是水准面,但接近于水准面,只是用于计算的辅助面。它与大地水准面不完全吻合,差值为正常高与正高之差。但在海洋面上时,似大地水准面与大地水准面重合。

关于水准面,我国规定采用的高程系统是正常高系统。如果不是进行科学研究,只是一般使用,正常高系 统结果在国内也可以称为海拔高度。正常高系统和正高系统是有区别的,主要是由于重力场的影响不同,重力线就会产生一些偏移。

美国人研发的全球定位系统GPS对空间坐标的描述,仍然采用WGS1984大地坐标系,他的高程为目标地物距离WGS84椭球体表面的法线距离。我国的地形图上的高程,标定的是目标地物距离大地水准面的铅锤距离

似大地水准面是前苏联地球物理学家、测量学家莫洛金斯基研究地球形状理论时,为避免大地水准面无法精确确定而引进的辅助面,为一与大地水准面十分接近、在海洋上两者完全重合、而在大陆上有2~4米的微小差异的曲面。由于正高与大地水准面的确定涉及到地球内部密度的假定,在理论上存在着不严密性,莫洛金斯理论作为现代大地测量里程碑,可以应用地面测量数据直接确定地球表面形状而不需要对地球密度作任何假设。似大地水准面只是通过一定的数学关系对应于地面的一个几何曲面,它既不是具有物理意义的水准面,也不是对于所有空间各点都为唯一的高程起算面。

我理解就是大地水准面是以海平面算出来的,似大地水准面是以陆地上的重力算出来的。

将水准面拟合成椭球面(将地球拟合成椭球体),拟合标准不同,则有了各种拟合的椭球面,比如克拉索夫斯基椭球体、1975年IUGG推荐的椭球体、美国WGS84椭球体,则有了不同标准下的大地高,但是因为各个拟合的椭球体有固定的旋转和偏移,因此是很容易就相互转换的。科研和国际常用的是WGS84定义的椭球体,美国规定的WGS84大地坐标系,力图让其WGS84椭球体在全球范围内整体拟合最优。中国规定的大地坐标系(如现行的2000国家大地坐标系),则尽量让椭球体对我国大陆地区地表拟合最优,因为我们国家陆地区域起伏太大。

我国现在用的是“2000国家大地坐标系”(其原点为包括海洋和大气的整个地球的质量中心。Z轴指向BIH1984.0定义的协议极地方向(BIH国际时间局),X轴指向BIH1984.0定义的零子午面与协议赤道的交点),我国的最佳拟合点,也称为大地原点,位于陕西省西安市泾阳县永乐镇

2000国家大地坐标系与WGS84坐标系都是大地坐标系

在定义上,CGCS2000和WGS84是一致的,坐标原点、尺度、定向及定向演变的定义都是相同的,两个坐标系使用的参考椭球也是非常相近的,在坐标系定义和实现上的比较,我们可以认为,WGS84和CGC2000是相容的,在坐标系的实现精度范围内,

WGS84和CGC2000的坐标是一致的。

北京54,西安80是参心坐标系,已经被淘汰,现在国内用2000国家大地坐标系,是地心坐标系。WGS84是地心坐标系。

大地坐标系是大地测量中以参考椭球面为基准面建立起来的坐标系。地面点的位置用大地经度、大地纬度和大地高度表示:(L, B, H)。

空间直角坐标系是以参考椭球中心为原点,以原点到0度经线与赤道交点的射线为x轴,原点到90度经线与赤道交点的射线为y轴,以地球旋转轴向北为z轴:(x, y, z)

共同点:显然,这两种坐标系都必须基于一个参考椭球。

不同点:大地坐标系以面为基准,所以还需要确定一个标准海平面。而空间直角坐标系则以一个点为基准,所以还需要确定一个中心点。

只要确定了椭球基本参数,则大地坐标系和空间直角坐标系就相对确定了,只是两种不同的表达而矣,这两个坐标系的点是一一对应的。

h54为1954年北京系大地高;h84为GPS WGS-84系大地高;H2为正常高;ε54为1954年北京系高程异常;ε84为WGS-84系高程异常;H1为正高;N54为大地水准面差异;N84为WGS-84系大地水准面差距。

由图,有以下关系式成立:

GPS测定的是大地高,要求解正常高必须先知道高程异常。在局部GPS网中巳知一些点 的高程异常(它由GPS水准算得), 考虑地球重力场模型,利用多面函数拟合法求定其它点的高程异常和正常高

即便我们忽略法线距离与铅锤距离的差异(即忽略垂线偏差角度),我们也会看到,这两个距离的起算基准可能不是同一个面。

GPS测出来的高程是大地高,地图上标的高程是正高或者正常高,这就是为什么同一个位置,我们用GPS测出来的高程与地形图上读出来的高程数值可能(通常)不一致。

但是我们一般忽略垂线偏差就角度,而将关系近似为如下:

大地水准面、似大地水准面、椭球面,示意图:

H大=H正+N

H大=H常+ξ

实践上,现在野外工作多采用GPS测定高度。GPS直接测定的是大地高,非海拔高。而一般的GPS似乎不能很理想地转化为海拔高,除非用控制点和已知参数做差分测量。那么,显然,用“高度”比用“海拔”更合适。当然,从效用的角度看,局部高度的变化,在大地高和海拔高的数值上,可能是相同的规律。

参考资料:

https://blog.csdn.net/weixin_33857230/article/details/86010321

https://blog.csdn.net/liuguobo/article/details/47093629

https://wenku.baidu.com/view/572e2a70b14e852459fb57a5.html

https://zhuanlan.zhihu.com/p/59743409

经纬高坐标系-ECEF坐标系-ENU坐标系相关推荐

  1. 经纬高、ECEF、ENU坐标系相互转换与实现

    ENU与ecef与WGS84相互转换(c++和python) - CodeAntenna 原理: 坐标系之间的转换关系(ECEF.LLA.ENU)_可即的博客-CSDN博客_lla坐标系

  2. LLA(经纬高)坐标转换成ENU(东北天)坐标的详细推导

    这是一篇经纬高(LLA)坐标转东北天坐标(ENU)的详细推导,并给出近似转换的过程和结果 参考资料: https://blog.csdn.net/qq_34213260/article/details ...

  3. 经纬高坐标转东北天坐标

    目录 1 问题描述 2 解决方案 2.1 经纬高转ECEF 2.2 ECEF转东北天 2.3 代入求解 2.3.1 东向化简 2.3.2 北向化简 2.3.3 天向化简 2.4 总结 1 问题描述 已 ...

  4. 经纬高坐标系转到东北天坐标系

    经纬高坐标系转到东北天坐标系 基本思路:首先把经纬高(大地坐标系.lla.llh)转到直角坐标系(地心地固直角坐标系(ECEF).xyz),然后再转为局部坐标系下(东北天坐标系.以第一点作为东北天坐标 ...

  5. 测站坐标系、地心非惯性系、经纬高互转

    目录 一.坐标系转换 1.1 测站坐标系rae(极坐标系)和空间直角坐标系xyz(东北天坐标系)互转 1.1.1.原理 1.1.2.公式 1.1.3.python代码 1.2 空间直角坐标系xyz(东 ...

  6. 高精地图生产中的坐标系

    漫话地图之高精地图生产中的坐标系_高德技术-程序员宝宝 - 程序员宝宝 (cxybb.com) 1.如何定义地球 1.1 地球椭球体 众所周知地球表面凸凹不平,是一个无法用数学公式表达的曲面,这样的曲 ...

  7. 漫话地图之高精地图生产中的坐标系

    导读 高精地图区别于普通电子地图,一般来说高精地图的精度要求在亚米级甚至厘米级,其所含有的道路交通元素也更丰富和细致.这就造成高精地图数据生产过程中会遇到很多普通地图没有的问题和困难,更需要有专业的测 ...

  8. 【cartographer】园区场景自动驾驶-利用GPS实现ENU坐标系下的建图与定位

    园区场景自动驾驶-利用GPS实现当地地理坐标系下的建图与定位 概述 主要实现如下的功能: 1.利用gps信息,修正题图为当地地理ENU坐标系: 2.在纯定位模式下,利用gps与imu(rtk)信息,实 ...

  9. ROS中ENU坐标系与无人机中NED坐标系的转换关系理解

    ROS中ENU坐标系与无人机中NED坐标系的转换关系理解 项目地址 无人机中NED坐标理解 ENU与NED转换 无人机中NED坐标理解 机体坐标系:机体坐标系固连飞机,其原点 取在多旋翼的重心位置上. ...

最新文章

  1. 阿里云服务器部署项目
  2. Response.Redirect:无法在发送 HTTP 标头之后进行重定向
  3. php 常用数组操作
  4. 最佳实践系列:前端代码标准和最佳实践
  5. CSS keylogger:攻击与防御
  6. Python的序列化与反序列化
  7. Android控件— — —ImageView
  8. [2018湖南省队集训] 6.28 T3 simulate
  9. linux的基础知识——进程组
  10. 关于JS特效的兼容问题。
  11. 吴恩达深度学习——浅层神经网络
  12. 日志记录总结要整理的地方
  13. 跨境电商个人物品清单申报开发代码
  14. [游泳] 游泳学习课程
  15. cad图纸怎么转换成图片?
  16. 2019最应该投资什么?是你明年的北大核心
  17. Mybatis-Plus eq、ne、gt、lt、ge、le分别代表含义
  18. 一文看懂语音技术商业化逻辑:对品牌意味着什么 | 综述
  19. PCB软件对比分析(AD、Pads、Allegro)
  20. C语言与C++初学“你好 世界”

热门文章

  1. GD32F4xx系列的看门狗配置
  2. python之items()的用法
  3. 浏览器急速模式跟兼容模式导致网站显示不同的问题
  4. Nginx 通过 certbot 为网站自动配置 SSL 证书并续期
  5. echarts v5去除内置的geoJSON
  6. docker 镜像的push 和 pull
  7. 基于SVM-支持向量机对鸢尾花数据进行分类
  8. 分割网络模型(FCN、Unet、Unet++、SegNet、RefineNet)
  9. oracle网络日志,oracle
  10. normalized与Normalize区别