坐标转换(三)—— ENU坐标系与机体坐标系之间的转换

2024-03-01 23:50

本文主要是介绍坐标转换(三)—— ENU坐标系与机体坐标系之间的转换,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

今天我们要讲的是应用于卫星通信领域、雷达目标跟踪领域、无源目标定位领域的一个坐标转换过程——机体坐标系与ENU坐标系的转换。小D将以《机载卫星通信天线角度的计算推导》为基础,来分析和推导坐标转换公式。

开头部分,我们先讲一下ENU坐标系和机体坐标系的定义。

1、ENU坐标系(导航坐标系)

ENU坐标系中的E表示East,N表示North,U表示Up,翻译过来就是东北天坐标系,小D第一次看到东北天这个名字时,还以为这个坐标系跟咱们祖国大东北有啥关系,后来发现没有关系。

所以,为了不让别人产生误解,小D就叫的北天东,其他文献中也有叫北东天的。另外还有一个坐标系是END,D表示Down。

ENU坐标系是以某个点为原点,三个轴表示出正东方向,正北方向和指向天上的方向,其中的“指向天上”是指垂直于水平面或海平面向上。

在大多数的文献中,ENU坐标系Z轴一般表示指向天,如果X轴指向正东方向,Y轴指向正北方向,那么这个坐标系就是ENU右手坐标系;如果Z轴指向天,X轴指向正北,Y轴指向正东,那么这个ENU坐标系就是左手坐标系。

我们今天参考的文献《机载卫星通信天线角度的计算推导》中的NEU坐标系是一个左手坐标系,并且ENU坐标系在这篇文献中也被称为导航坐标系,记作n系。ENU左手坐标系示意图如下图所示:

图1 ENU左手坐标系(导航坐标系n系)

2、机体坐标系(载体坐标系)

机体坐标系也叫载体坐标系,记作b系。机体坐标系也有很多种定义方式,有的地方定义的是左手坐标系,有的地方是右手坐标系;有的机体Z轴垂直于水平面向上,有的垂直于水平面向下。

我们今天讲的机体坐标系是这样定义的:X轴正方向指向飞机前进的方向,与机身纵轴重合,Y轴指向机身右翼,Z轴垂直于机身平面向上,构成左手坐标系。

航向角、俯仰角、横滚角

机体坐标系中有个重要的概念,就是机体姿态角,姿态角包含三个角度:航向角(yaw),俯仰角(pitch),横滚角(roll),这三个角的定义如下:

航向角α机头方向在水平面的投影与正北之间的夹角,北偏东为正,范围为0-360°;也可定义右偏为正,左偏为负,范围为-180-180°。

俯仰角β:机头方向与水平面的夹角,抬头为正,范围为-90-90°。

横滚角γ:机体绕机体纵轴旋转过的角度,右滚为正,范围为-180-180°。

其实上面各个角度都是人为定义的,你也可以有不同的定义方式,重要的是在解算角度的时候,实际的偏离角度要和定义的角度一致。

机体坐标系的航向角、俯仰角和横滚角示意图如下所示:

图2 飞机航向角、俯仰角和横滚角图解

上面的图是小D用visio画的,应该可以很好地说明航向角、俯仰角和横滚角。 

3、ENU左手坐标系转换到机体左手坐标系

ENU左手坐标系(导航坐标系)转到机体左手坐标系(载体坐标系)是一个有点抽象的过程,小D第一次看这个转换过程的时候,觉得头都大了,想了好久才想明白是怎么旋转的。为了讲清楚转换原理,小D根据《机载卫星通信天线角度的计算推导》里的图,重新画了一个三维图,希望可以让大家更好地理解这个过程。

下图中,α是飞机的航向角,β是俯仰角,γ是横滚角,大家可以结合小D参考的这篇文献中的解释来看这几个角,小D把截图放在文末了。

紫色的XnYnZn坐标系代表NEU坐标系,灰色的XbYbZb坐标系代表机体坐标系,虚线坐标系是中间过程坐标系。为了不让示意图看起来很复杂,文献以及小D都把ENU坐标系XnYnZn的原点放到了机尾位置,它本来应该是和机体坐标系XbYbZb共用一个原点的。

图3 ENU坐标系转换到机体坐标系的过程示意图

现在我们来看看ENU坐标系XnYnZn转换到机体坐标系XbYbZb的整个旋转过程,一共有三个旋转变换:

①XnYnZn坐标系先绕Zn轴顺时针旋转α角,形成X´nY´nZ´n坐标系;

②X´nY´nZ´n坐标系绕Y´n轴逆时针旋转β角,形成X´´nY´´nZ´´n坐标系;

③X´´nY´´nZ´´n坐标系绕X´´n轴逆时针旋转γ角,与机体坐标系XbYbZb重合。

图中没有画出第③个过程,大家可以自己想像一下,或者用模型比划一下,模型小D是用便签纸和胶棒粘的一个三维的空间XYZ坐标系,用来理解空间中的变换过程,大家也可以像小D一样做一个模型。

图片

图4 自制模型

由于上面的过程是左手坐标系的旋转过程,所以顺时针旋转时,计算带入的角度是正的,逆时针旋转计算带入的角度是负的,根据小D上次发的公众号文章中的左手坐标系的旋转矩阵,我们可以计算ENU坐标系XnYnZn转换到机体坐标系XbYbZb的旋转矩阵:

图片

这就是最终的旋转矩阵了。

4、机体左手坐标系转换到ENU左手坐标系

机体坐标系转到ENU坐标系的过程是一个逆过程:

XbYbZb坐标系先绕Xb轴顺时针旋转γ角,形成X´bY´bZ´b坐标系;

②X´bY´bZ´b坐标系绕Y´b轴顺时针旋转β角,形成X´´bY´´bZ´´b坐标系;

X´´bY´´bZ´´b坐标系绕Z´´b轴逆时针旋转α角,与ENU坐标系XnYnZn重合。

所以,机体左手坐标系转换到ENU左手坐标系旋转矩阵为:

图片

5、文献中的公式推导

《机载卫星通信天线角度的计算推导》中的转换过程示意图如下图所示,其中U轴是水平面正上方,小D绘制的跟文献中的图有一些不同的地方,就是Zb轴小D觉得应该画在U轴的右边。

图片

图5 文献中的旋转过程示意图

除此之外,小D觉得文献中给出的绕X,Y,Z轴逆时针旋转的旋转矩阵是有点问题的,正确的旋转矩阵小D在上篇公众号文章中给出了,大家可以对比一下,如果你觉得小D的观点有问题,也可以告诉小D。

图片

图片

文献中关于ENU坐标系(导航坐标系)与机体坐标系(载体坐标系)的推导过程是这样的:

图片

图片

需要注意的是,上面的“逆时针旋转2π-α”与“顺时针旋转α”是等效的。该篇文献的作者讲解的旋转过程是导航坐标系转到载体坐标系的过程,但是计算的旋转矩阵是载体坐标系转到导航坐标系的公式。

还有一点就是,作者可能是笔误,三个矩阵相乘的结果也有一点问题,小D用matlab计算了一下作者给出的三个矩阵相乘的公式,计算过程如下:

图片

图6 旋转矩阵相乘的验证代码

图片

图7 代码计算结果 

代码计算结果如上图所示,和作者计算的不一样,字体有点小,大家可以自己在matlab里验证旋转矩阵计算结果。

虽然文献有一些小问题,但是文献中的示意图画得很形象,小D就是通过这篇文献理解了整个旋转过程。

如果大家觉得小D讲的坐标转换有问题,欢迎在评论区讨论。

参考文献

[1] 李川.《机载卫星通信天线角度的计算推导》. 天津市电子工业协会2022年年会.

[2] 党亚民.《大地坐标系及其应用》, 测绘出版社, 2010年10月1日第1版.

[3] 黄珹、刘林. 《参考坐标系及航天应用》, 电子工业出版社, 2015年5月1日.

6、期刊论文及书籍PDF下载链接

①《机载卫星通信天线角度的计算推导》:

https://pan.baidu.com/s/1dD94KeA8OFCJwxOcp3ozRA?pwd=nvl9 
提取码:nvl9 

②《大地坐标系及其应用》:

链接:https://pan.baidu.com/s/1Min-PZcR6_hl2uyQq7HTdQ?pwd=8r3h

提取码:8r3h

③《参考坐标系及航天应用》:

链接:https://pan.baidu.com/s/1ocTwQptGg0UY1qsBwfHpwg?pwd=ozmb 
提取码:ozmb 

这篇关于坐标转换(三)—— ENU坐标系与机体坐标系之间的转换的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/764243

相关文章

C#实现将Excel表格转换为图片(JPG/ PNG)

《C#实现将Excel表格转换为图片(JPG/PNG)》Excel表格可能会因为不同设备或字体缺失等问题,导致格式错乱或数据显示异常,转换为图片后,能确保数据的排版等保持一致,下面我们看看如何使用C... 目录通过C# 转换Excel工作表到图片通过C# 转换指定单元格区域到图片知识扩展C# 将 Excel

C++使用printf语句实现进制转换的示例代码

《C++使用printf语句实现进制转换的示例代码》在C语言中,printf函数可以直接实现部分进制转换功能,通过格式说明符(formatspecifier)快速输出不同进制的数值,下面给大家分享C+... 目录一、printf 原生支持的进制转换1. 十进制、八进制、十六进制转换2. 显示进制前缀3. 指

使用Python开发一个带EPUB转换功能的Markdown编辑器

《使用Python开发一个带EPUB转换功能的Markdown编辑器》Markdown因其简单易用和强大的格式支持,成为了写作者、开发者及内容创作者的首选格式,本文将通过Python开发一个Markd... 目录应用概览代码结构与核心组件1. 初始化与布局 (__init__)2. 工具栏 (setup_t

Java中Date、LocalDate、LocalDateTime、LocalTime、时间戳之间的相互转换代码

《Java中Date、LocalDate、LocalDateTime、LocalTime、时间戳之间的相互转换代码》:本文主要介绍Java中日期时间转换的多种方法,包括将Date转换为LocalD... 目录一、Date转LocalDateTime二、Date转LocalDate三、LocalDateTim

golang获取当前时间、时间戳和时间字符串及它们之间的相互转换方法

《golang获取当前时间、时间戳和时间字符串及它们之间的相互转换方法》:本文主要介绍golang获取当前时间、时间戳和时间字符串及它们之间的相互转换,本文通过实例代码给大家介绍的非常详细,感兴趣... 目录1、获取当前时间2、获取当前时间戳3、获取当前时间的字符串格式4、它们之间的相互转化上篇文章给大家介

Python实现AVIF图片与其他图片格式间的批量转换

《Python实现AVIF图片与其他图片格式间的批量转换》这篇文章主要为大家详细介绍了如何使用Pillow库实现AVIF与其他格式的相互转换,即将AVIF转换为常见的格式,比如JPG或PNG,需要的小... 目录环境配置1.将单个 AVIF 图片转换为 JPG 和 PNG2.批量转换目录下所有 AVIF 图

详解如何通过Python批量转换图片为PDF

《详解如何通过Python批量转换图片为PDF》:本文主要介绍如何基于Python+Tkinter开发的图片批量转PDF工具,可以支持批量添加图片,拖拽等操作,感兴趣的小伙伴可以参考一下... 目录1. 概述2. 功能亮点2.1 主要功能2.2 界面设计3. 使用指南3.1 运行环境3.2 使用步骤4. 核

Java实现时间与字符串互相转换详解

《Java实现时间与字符串互相转换详解》这篇文章主要为大家详细介绍了Java中实现时间与字符串互相转换的相关方法,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录一、日期格式化为字符串(一)使用预定义格式(二)自定义格式二、字符串解析为日期(一)解析ISO格式字符串(二)解析自定义

在java中如何将inputStream对象转换为File对象(不生成本地文件)

《在java中如何将inputStream对象转换为File对象(不生成本地文件)》:本文主要介绍在java中如何将inputStream对象转换为File对象(不生成本地文件),具有很好的参考价... 目录需求说明问题解决总结需求说明在后端中通过POI生成Excel文件流,将输出流(outputStre

python+opencv处理颜色之将目标颜色转换实例代码

《python+opencv处理颜色之将目标颜色转换实例代码》OpenCV是一个的跨平台计算机视觉库,可以运行在Linux、Windows和MacOS操作系统上,:本文主要介绍python+ope... 目录下面是代码+ 效果 + 解释转HSV: 关于颜色总是要转HSV的掩膜再标注总结 目标:将红色的部分滤