Sophus降维、升维与欧拉角、旋转向量的爱恨情仇

news2025/1/19 3:36:35

0. 简介

在面对二维与三维之间的转换时,我们常常会困惑该如何去转换,在G2O中存在有理想的坐标转换工具,但是在Sophus中却缺乏这样的手段。之前在Sophus处简要的介绍了一下SE(2)与SE(3)的转换,最近发现之前的文章这部分需要拿出来详细的说一说。

1. 欧拉角与旋转向量

欧拉角、旋转向量、四元数和旋转矩阵是Sophus中常常提到的几个名词,欧拉角和旋转向量是类似的,SO(3)的旋转矩阵有9个量,但是只有3个自由度,并且是单位正交矩阵,具有冗余性,对其估计或优化问题的求解不方便。我们可以用一个旋转轴和一个旋转角描述任意旋转。一个方向与旋转轴一致,长度(模)等于旋转角的向量,我们称之为旋转向量(或轴角)

旋转向量到旋转矩阵:

R = c o s θ I + ( 1 − c o s θ ) n n T + s i n θ n ^ R = cos\theta I+(1-cos\theta)nn^T+sin\theta \hat{n} R=cosθI+(1cosθ)nnT+sinθn^

其中提到的 n ^ \hat{n} n^是向量的旋转矩阵,由于是反对称矩阵,所以只存在三个自由度,我们可以 X − Y − Z X-Y-Z XYZ轴来规定其反对称矩阵轴:
n ^ = [ 0 − n z n y n z 0 − n x − n y n x 0 ] \hat{n}=\begin{bmatrix} 0 && -{n_z} && {n_y} \\ {n_z} && 0 && -{n_x} \\ -{n_y} && {n_x} && 0 \end{bmatrix} n^= 0nznynz0nxnynx0

所以我们可以从se,so中得到旋转向量数值。而欧拉角对于轴角表示情况,转轴具有2个自由度,转角1个自由度。 根据三次基本转动选取的坐标轴的不同,欧拉角共有12种组合。如 果再考虑到可选取原始坐标系的坐标轴,也可选取“新”坐标系的坐标轴,则共有24种欧拉角表示。一般规定原始坐标系为静坐标系,每个基本转动后形成的新坐标系为动坐标系。

  • 24 种欧拉角表示列举如下:
  • 静轴(即转轴选静坐标系的坐标轴):
    s X Y Z , s X Z Y , s X Y X , sXYZ,sXZY,sXYX, sXYZ,sXZY,sXYX,
    s X Z X , s Y X Z , s Y Z X , sXZX,sYXZ,sYZX, sXZX,sYXZ,sYZX,
    s Y X Y , s Y Z Y , s Z X Y , sYXY,sYZY,sZXY, sYXY,sYZY,sZXY,
    s Z Y X , s Z X Z , s Z Y Z sZYX,sZXZ,sZYZ sZYX,sZXZ,sZYZ

动轴(即转轴选动坐标系的坐标轴):

r Z Y X , r Y Z X , r X Y X , rZYX,rYZX,rXYX, rZYX,rYZX,rXYX,
r X Z X , r Z X Y , r X Z Y , rXZX,rZXY,rXZY, rXZX,rZXY,rXZY,
r Y X Y , r Y Z Y , r Y X Z , rYXY,rYZY,rYXZ, rYXY,rYZY,rYXZ,
r X Y Z , r Z X Z , r Z Y Z rXYZ,rZXZ,rZYZ rXYZ,rZXZ,rZYZ

静轴欧拉角和动轴欧拉角有如下规律:
绕静轴 X Y Z XYZ XYZ 分别 转 α , β , γ α,β,γ α,β,γ 角度的转动与绕动轴 Z Y X ZYX ZYX分别转 γ , β , α γ,β,α γ,β,α 角度的转动等价,其他形式的欧拉角亦有此类似规律。

对于不同的坐标系定义,有不同的转换关系。我们只讨论常用的一种情况:如上图,右手系,Z轴朝上,X轴朝前,y轴朝左。绕Z轴作偏航(Yaw)运动,绕Y轴作俯仰(Pitch)运动,绕X轴作滚转(Roll)运动,运动正方向如上图所示。在这里插入图片描述
对于欧拉角计算公式我们可以得到 α , β , γ α,β,γ α,β,γ 三个角度计算得到的旋转矩阵
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
上式中 R R R与旋转次序有关,即当 θ z , θ y , θ x θz,θy,θx θzθyθx不都为小角时,对应于不同的旋转次序,空间坐标系b的最终位置时不同的,这就是有限转动的不可交换性。但是当 θ z , θ y , θ x θz,θy,θx θzθyθx都为小角时,忽略小角间的高阶小量,即: Δ x → 0 , c o s Δ x → 1 , s i n Δ x → Δ x Δx→0,cosΔx→1,sinΔx→Δx Δx0,cosΔx1sinΔxΔx
R ≈ n ^ = [ 0 − n z n y n z 0 − n x − n y n x 0 ] R \approx \hat{n}=\begin{bmatrix} 0 && -{n_z} && {n_y} \\ {n_z} && 0 && -{n_x} \\ -{n_y} && {n_x} && 0 \end{bmatrix} Rn^= 0nznynz0nxnynx0

其中 θ z , θ y , θ x θz,θy,θx θzθyθx角度为弧度,此时 θ z , θ y , θ x θz,θy,θx θzθyθx构成的列向量 [ θ z , θ y , θ x ] T [θz,θy,θx]^T [θzθyθx]T可视为三维空间的(旋转)矢量,此时旋转后的坐标系的最终角位置与旋转次序无关:无限转动与旋转次序无关

2. SE(2)与SE(3)的转换

SE(2),通常是作为二维向量的表示形式,基本的组成部分为 x , y , y a w x,y,yaw x,y,yaw三参数,其格式如下:
在这里插入图片描述
SE(3)则是在上式的基础上加入 z z z轴,其格式如下:

[ x w y w z w 1 ] = [ c 1 c 3 + s 1 s 2 s 3 c 3 s 1 s 2 − c 1 s 3 c 2 s 1 t w x c 2 s 3 c 2 c 3 − s 2 t w y c 1 s 2 s 3 − s 1 c 3 s 1 s 3 + c 1 c 3 s 2 c 1 c 2 t w z 0 0 0 1 ] × [ x y z 1 ] \begin{bmatrix} x_w \\ y_w \\ z_w \\ 1\end{bmatrix} =\begin{bmatrix} {c_1 c_3 + s_1 s_2 s_3} && {c_3 s_1 s_2 - c_1 s_3} && {c_2 s_1} && t_wx \\ {c_2 s_3} && {c_2 c_3} && -{s_2} && t_wy \\ {c_1 s_2 s_3 - s_1 c_3} && {s_1 s_3 + c_1 c_3 s_2} && {c_1 c_2} && t_wz \\ 0 && 0 && 0 && 1\end{bmatrix} \times \begin{bmatrix} x \\ y \\ z \\ 1\end{bmatrix} xwywzw1 = c1c3+s1s2s3c2s3c1s2s3s1c30c3s1s2c1s3c2c3s1s3+c1c3s20c2s1s2c1c20twxtwytwz1 × xyz1

在这里插入图片描述

李群SO(3)SE(3)
旋转矩阵构建Sophus::SO3d SO3(R)Sophus::SE3d SE3(R,t)
四元数构建Sophus::SO3d SO3(q)Sophus::SO3d SO3(q,t)
输出SO3.matrix()SE3.matrix()
对数映射Vector3d so3=SO3.log()Vecotr6d se3=SE3.log()
指数映射SO3d::exp(so3)SE3d::exp(se3)
向量到反对称矩阵SO3d::hat(so3)SE3d::hat(se3)
反对称矩阵到向量SO3d::vee(hat)SE3d::vee(hat)

3. 示例代码

…详情请参照古月居

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

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

相关文章

Leetcode:654. 最大二叉树(C++)

目录 问题描述: 实现代码与解析: 直接模拟(递归): 原理思路: 索引版本: 问题描述: 给定一个不重复的整数数组 nums 。 最大二叉树 可以用下面的算法从 nums 递归地构建: 创建一个根节点&a…

Thymeleaf从入门到清晰使用

文章目录什么是thymeleaf?第一个Thymeleaf程序Thymeleaf详解配置常用标签最后什么是thymeleaf? 模板引擎: 前端交给我们的页面,是html页面,如果是以前开发,我们需要把他们转成jsp页面,jsp的好处…

ABB AC500 系列 PLC 与上位机iFix 的通讯配置

ABB PLC IP 及 MODBUS TCP/IP 协议设置 通过 IP config tool 配置设备 IP 在 软件中,有 3 种方式可以进入 IP config tool 的配置界面  双击 IP_settings,点击 IP config tool,即可进入 IP config tool 界面 点击菜单栏的 Tool→IP confi…

k8s 微服务spring boot JVM 监控

目录 1.前言 1.1 JVM监控方案 1.2 接入操作步骤 2. 服务开启actuator prometheus监控 2.1 示例参考添加依赖 2.2 配置prometheus监控 3 配置 prometheus 监控 4 配置jvm grafana dashboard 1.前言 服务 部署环境 监控 spring-demo k8s v1.22 Prometheus Operator …

进程相关概念

1、什么是程序,什么是进程,有什么区别 程序是静态的概念。例如:gcc xxx.c -o pro,磁盘中生成pro文件,叫做程序(还未运行起来) 进程是程序的一次运行活动,也就是程序跑起来了&#xf…

【Linux】信号机制(非实时信号)

目录 前言 一.信号的概念以及产生 1.什么是信号 2.信号分为两类 3.查看信号的命令 4.信号如何产生 1).通过软件产生 2).通过硬件产生 3).通过键盘组合键产生 二.信号的发送以及保存 1.信号如何发送 2.信号如何保存 1).概念 2).底层实现结构&&内核中的实现…

WampServer服务器设置控件

WampServer服务器设置控件 WampServer是一个web开发环境,是一个用于Windows操作系统的类似服务器的软件,由Romain Bourdon构建。该工具允许您开始构建web应用程序,并在Windows上使用Apache服务器2、PHP编程语言和MySQL数据库在本地网络上启动…

java 实现 springboot项目 使用socket推送消息,前端实时进行接收后端推送的消息

1 后端 1.1 添加依赖 在我们的springboot项目里面&#xff0c;添加依赖&#xff1b; <dependency><groupId>org.springframework.boot</groupId><artifactId>spring-boot-starter-websocket</artifactId><version>2.7.0</version>…

企业微信在职继承功能如何理解?怎么使用?

企业微信的在职继承功能就是企业管理员将在职员工的客户分配给其他在职员工跟进的一种方式&#xff0c;很方便&#xff0c;可以在工作出现变更时使用。 前言 企业微信之前的离职继承功能受到很多企业的青睐&#xff0c;企业能够将离职员工的客户分配给其他员工接手&#xff0c…

云天励飞在科创板获准注册:计划募资30亿元,陈宁为实际控制人

2023年1月10日&#xff0c;证监会发布公告&#xff0c;同意深圳云天励飞技术股份有限公司&#xff08;下称“云天励飞”&#xff09;首次公开发行股票注册。据贝多财经了解&#xff0c;云天励飞于2020年12月8日在科创板递交上市申请&#xff0c;2021年8月6日过会。 本次冲刺上市…

蓝队攻击的四个阶段(一)

蓝队的攻击是一项系统的工作&#xff0c;整个攻击过程是有章可循、科学合理的涵盖了从前期准备、攻击实施到靶标控制的各个步骤和环节。按照任务进度划分&#xff0c;一般可以将蓝队的工作分为4个阶段:准备工作、目标网情搜集、外网络向突破和内网横向拓展。第一阶段准备工作&a…

App开发提效法宝之插件技术

近年来技术革新频率越来越高&#xff0c;最近工作中经常有小伙伴问到插件技术的相关内容&#xff0c;今天就来跟大家系统的说清楚什么是插件技术以及它的好处。欢迎评论区交流哦&#xff01; 什么是插件技术&#xff1f; 插件技术指的是一种应用程序&#xff0c;遵循程序接口…

高温热水解预处理对厌氧消化期间污泥腐殖化的调控机制

期刊&#xff1a;Water Research 影响因子&#xff1a;9.13 发表时间&#xff1a;2022样本类型&#xff1a;污泥客户单位&#xff1a;同济大学凌恩生物客户同济大学发表在《Water Research》上的文章“The neglected effects of polysaccharide transformation on sludge humif…

振弦采集模块参数配置工具VMTool 的使用

振弦采集模块参数配置工具VMTool 的使用 准备工作 &#xff08; 1&#xff09; 将 VMXXX 模块的 UART_TTL、 RS232&#xff08; 或 RS485&#xff09; 接口与计算机的 COM 端口连接&#xff1b; &#xff08; 2&#xff09; 连接振弦传感器及温度传感器到 VMXXX 的对应接口&…

三菱FX3U系列PLC运动控制_伺服回原点的3种方法示例

三菱FX3U系列PLC运动控制_伺服回原点的3种方法示例 方法1: 运动的方向为圆形、环形、电机往一个方向转动; 只有一个原点开关,没有极限开关 如下图所示, 原点回归的方式为:启动回原点后,电机开始寻找原点开关,在找到原点开关的瞬间,开始减速;在离开原点开关的瞬间,电…

23. 反爬案例:不登录不给,要数据请先登录我的站点

登录之后&#xff0c;可以查看数据&#xff0c;是部分站点常用规则&#xff0c;本篇博客将在爬虫训练场中实现该需求。 文章目录安装必备模块建立 models建立 login_form 表单文件flask_wtf 中 FlaskForm 类建立登录视图函数配置 login.html 页面安装必备模块 实现 Python Fla…

Qt基于CTK Plugin Framework搭建插件框架--插件通信【事件监听】

文章目录一、前言二、事件三、类通信3.1、新建接收插件3.2、新建发送插件3.3、启用插件四、信号槽通信4.1、新建接收插件4.2、新建发送插件4.3、启用插件五、类通信和信号槽通信的区别六、插件依赖七、获取元数据一、前言 CTK框架中的事件监听&#xff0c;其实就是观察者模式&…

分享微信答题抽奖小程序制作步骤_可以做答题后抽奖活动吗

知识答题小程序如何制作&#xff1f;现在越来越多的企业和组织逐步进行各种获奖知识问答小程序。那么&#xff0c;如何制作一个答题小程序呢&#xff1f;今天&#xff0c;我们一起来看看~需要的老板不要走开哦~既然点进来了&#xff0c;那就请各位老板耐心看到最后吧~怎么做一个…

Linux工具学习之【git】

✨个人主页&#xff1a; Yohifo &#x1f389;所属专栏&#xff1a; Linux学习之旅 &#x1f38a;每篇一句&#xff1a; 图片来源 &#x1f383;操作环境&#xff1a; CentOS 7.6 阿里云远程服务器 Whatever is worth doing is worth doing well. 任何值得去做的事情&#xff0…

实现高并发秒杀的方式

实现高并发秒杀的方式 引言商品秒杀-超卖解决商品超卖 方式一&#xff08;改进版加锁&#xff09; 方式二&#xff08;AOP版加锁&#xff09; 方式三&#xff08;悲观锁一&#xff09; 方式四&#xff08;悲观锁二&#xff09; 方式五&#xff08;乐观锁&#xff09; 方式六&a…