旋转变换(一)旋转矩阵

1. 简介

计算机图形学中的应用非常广泛的变换是一种称为仿射变换的特殊变换,在仿射变换中的基本变换包括平移、旋转、缩放、剪切这几种。本文以及接下来的几篇文章重点介绍一下关于旋转的变换,包括二维旋转变换、三维旋转变换以及它的一些表达方式(旋转矩阵、四元数、欧拉角等)。

2. 绕原点二维旋转

首先要明确旋转在二维中是绕着某一个点进行旋转,三维中是绕着某一个轴进行旋转。二维旋转中最简单的场景是绕着坐标原点进行的旋转,如下图所示:

如图所示点v 绕 原点旋转θ 角,得到点v’,假设 v点的坐标是(x, y) ,那么可以推导得到 v’点的坐标(x’, y’)(设原点到v的距离是r,原点到v点的向量与x轴的夹角是ϕ ) 
x=rcosϕ    y=rsinϕ 
x′=rcos(θ+ϕ)   y′=rsin(θ+ϕ
通过三角函数展开得到 
x′=rcosθcosϕrsinθsinϕ 
y′=rsinθcosϕ+rcosθsinϕ 
带入x和y表达式得到 
x′=xcosθysinθ 
y′=xsinθ+ycosθ 
写成矩阵的形式是:

尽管图示中仅仅表示的是旋转一个锐角θ的情形,但是我们推导中使用的是三角函数的基本定义来计算坐标的,因此当旋转的角度是任意角度(例如大于180度,导致v’点进入到第四象限)结论仍然是成立的。

3. 绕任意点的二维旋转

绕原点的旋转是二维旋转最基本的情况,当我们需要进行绕任意点旋转时,我们可以把这种情况转换到绕原点的旋转,思路如下: 
1. 首先将旋转点移动到原点处 
2. 执行如2所描述的绕原点的旋转 
3. 再将旋转点移回到原来的位置

也就是说在处理绕任意点旋转的情况下需要执行两次平移的操作。假设平移的矩阵是T(x,y),也就是说我们需要得到的坐标 v’=T(x,y)*R*T(-x,-y)(我们使用的是列坐标描述点的坐标,因此是左乘,首先执行T(-x,-y))

在计算机图形学中,为了统一将平移、旋转、缩放等用矩阵表示,需要引入齐次坐标。(假设使用2x2的矩阵,是没有办法描述平移操作的,只有引入3x3矩阵形式,才能统一描述二维中的平移、旋转、缩放操作。同理必须使用4x4的矩阵才能统一描述三维的变换)。

对于二维平移,如下图所示,P点经过x和y方向的平移到P’点,可以得到:

x′=x+tx

y′=y+ty 
由于引入了齐次坐标,在描述二维坐标的时候,使用(x,y,w)的方式(一般w=1),于是可以写成下面矩阵的形式

按矩阵乘法展开,正好得到上面的表达式。也就是说平移矩阵是

如果平移值是(-tx,-ty)那么很明显平移矩阵式

我们可以把2中描述的旋转矩阵也扩展到3x3的方式,变为:

从平移和旋转的矩阵可以看出,3x3矩阵的前2x2部分是和旋转相关的,第三列与平移相关。有了上面的基础之后,我们很容易得出二维中绕任意点旋转的旋转矩阵了,只需要把三个矩阵乘起来即可:

4. 三维基本旋转

我们可以把一个旋转转换为绕基本坐标轴的旋转,因此有必要讨论一下绕三个坐标值x、y、z的旋转。 
本文在讨论过程中使用的是类似于OpenGL中定义的右手坐标系,同时旋转角度的正负也遵循右手坐标系的约定。如下图所示

4.1 绕X轴的旋转

在三维场景中,当一个点P(x,y,z)绕x轴旋转θ角得到点P’(x’,y’,z’)。由于是绕x轴进行的旋转,因此x坐标保持不变,y和z组成的yoz(o是坐标原点)平面上进行的是一个二维的旋转,可以参考上图(y轴类似于二维旋转中的x轴,z轴类似于二维旋转中的y轴),于是有: 
x′=x 
y′=ycosθzsinθ 
z′=ysinθ+zcosθ 
写成(4x4)矩阵的形式

4.2 绕Y轴旋转

绕Y轴的旋转和绕X轴的旋转类似,Y坐标保持不变,除Y轴之外,ZOX组成的平面进行一次二维的旋转(Z轴类似于二维旋转的X轴,X轴类似于二维旋转中的Y轴,注意这里是ZOX,而不是XOZ,观察上图中右手系的图片可以很容易了解到这一点),同样有: 
x′=zsinθ+xcosθ 
y′=y 
z′=zcosθxsinθ 
写成(4x4)矩阵的形式

4.3 绕Z轴旋转

与上面类似,绕Z轴旋转,Z坐标保持不变,xoy组成的平面内正好进行一次二维旋转(和上面讨论二维旋转的情况完全一样)

4.4 小结

上面描述了三维变换中绕单一轴旋转的矩阵表达形式,绕三个轴旋转的矩阵很类似,其中绕y轴旋转的矩阵与绕x和z轴旋转的矩阵略有点不同(主要是三个轴向顺序和书写矩阵的方式不一致导致的,绕三个不同坐标旋转轴以及其他二个坐标轴组成平面的顺序是: XYZ(绕x轴) YZX(绕y轴) ZXY(绕z轴),其中绕y轴旋转,其他两个轴是ZX,这和我们书写矩阵按

的方式不一致,而导致看起来绕Y轴旋转的矩阵似乎是和其他两个矩阵不一致。如果我们颠倒写法,将公式写成

的方式,那么这三个旋转矩阵看起来在形式上就统一了,都是

这种表现形式了(左上角都是−sinθ

5. 绕任意轴的三维旋转

绕任意轴的三维旋转可以使用类似于绕任意点的二维旋转一样,将旋转分解为一些列基本的旋转。绕任意轴旋转如下图所示:

P点绕向量u旋转θ角,得到点Q,已知P点的坐标和向量u,如何求Q点的坐标。 
我们可以把向量u进行一些旋转,让它与z轴重合,之后旋转P到Q就作了一次绕Z轴的三维基本旋转,之后我们再执行反向的旋转,将向量u变回到它原来的方向,也就是说需要进行的操作如下: 
1. 将旋转轴u绕x轴旋转至xoz平面 
2. 将旋转轴u绕y轴旋转至于z轴重合 
3. 绕z轴旋转θ角 
4. 执行步骤2的逆过程 
5. 执行步骤1的逆过程 
原始的旋转轴u如下图所示:

第1、2、3步骤如下图所示:

步骤1将向量u旋转至xoz平面的操作是一个绕x轴的旋转操作,步骤2将向量u旋转到与z轴重合,第1、2步骤的示意图如下:

作点P在yoz平面的投影点q,q的坐标是(0, b, c),原点o与q点的连线oq和z轴的夹角就是u绕x轴旋转的角度。通过这次旋转使得u向量旋转到xoz平面(图中的or向量)【步骤1】 
过r点作z轴的垂线,or与z轴的夹角为β, 这个角度就是绕Y轴旋转的角度,通过这次旋转使得u向量旋转到与z轴重合【步骤2】

步骤1中绕x轴旋转的是一次基本的绕x轴的三维旋转,按照之前的讨论,旋转矩阵是:

这里的θ就是图中所示的α角 (注意α角度是绕x旋转的正的角度) 
从图中我们还可以得到:

于是旋转矩阵(记作 Rx(α))为:

在完成步骤1之后,向量u被变换到了r的位置,我们继续步骤2的操作,绕y轴旋转负的β角(注意:这里的β是负的),经过这次变换之后向量u与z轴完全重合,由于这一步也是执行的一次绕Y轴的基本旋转,旋转矩阵(记作 Ry(−β))为:

使用−β替换表达式中的θ,此外根据图中描述,我们可以计算得到:

带入上面的表达式,于是旋转矩阵(记作 Ry(−β))为:

在完成前面两个步骤之后,u方向和z轴完全重合,因此执行旋转θ角,执行的是一次绕z轴的基本三维旋转(记作 R(θ),根据之前的讨论,我们可以得到:

最后两步骤是前面1和2的逆操作,也就是绕Y轴旋转β 和绕X轴旋转−α,这两个矩阵分别记作 Ry(β) 和 Rx(−α),得到它们的方式很简单,只需要将上面步骤1和步骤2中的角度修改成相反数即可,也就是:

最终得到 绕任意轴u旋转的旋转矩阵是【因为使用的列向量,因此执行的是左乘(从右往左)】:

MR=Rx(−α)Ry(β)Rz(θ)Ry(−β)Rx(α)=

(注:式中的(u,v,w)对应上文中向量(a,b,c),公式我自己笔算过,为了减少编辑公式的时间(使用LaTex编辑太繁琐,因此找了一张公式的图片贴在此处)

如果向量是经过单位化的(单位向量),那么有a2+b2+c2=1,可以简化上述的公式,得到:

旋转(一):旋转矩阵推导相关推荐

  1. 点绕坐标系旋转旋转矩阵推导、理解

    点绕坐标系旋转旋转矩阵推导.理解 开会过程中,翻笔记本时突然发现了之前推导的点绕坐标系旋转的旋转矩阵求解过程.在直觉上,一个点绕3个坐标轴旋转的旋转矩阵在形式上看起来是一致的,观感上比较优雅.整齐,但 ...

  2. 运动学解析~旋转矩阵推导

    旋转矩阵的推导要注意最关键的一点,就是我们需要的是坐标系的旋转变换,而不是点的旋转变换.大部分解释是第二条推导,所以感觉上就不能说服人. 本文转自这个链接:旋转矩阵怎么推导_如何推导旋转矩阵_独孤小白 ...

  3. 3维旋转矩阵推导与助记

    旋转矩阵的应用范围比较广,是姿态变换,坐标变换等的基础.本篇先介绍旋转矩阵的推导过程与助记方法. 旋转矩阵的旋转其实包含两种意思,一是在同一个坐标系下,向量的旋转:二是坐标系的旋转,使得同一向量在不同 ...

  4. 坐标系旋转矩阵推导过程

    一.先来个平面旋转的分析: 两角和(差)公式 推导 旋转变换一般是按照某个圆心点,以一定半径 r 旋转一定的角度α,为了简单起见我们给出下面的情景 假定点A(x,y)想经过旋转变换到达B(x',y') ...

  5. 坐标系旋转与向量旋转(旋转矩阵)

    二维坐标系旋转与向量(坐标)旋转 坐标系的旋转和向量的旋转在工程应用过程中经常会遇到,在这里对二维坐标系的旋转和向量旋转做一个简单的推导,方便大家的理解. 坐标系旋转 一个平面坐标系逆时针旋转一个角度 ...

  6. 刚体在三维空间的旋转(关于旋转矩阵、DCM、旋转向量、四元数、欧拉角)

    最近学习了一些关于三维空间旋转相关的知识,借此梳理一下备忘. 三维空间的旋转(3D Rotation)是一个很神奇的东东:如果对某个刚体在三维空间进行任意次的旋转,只要旋转中心保持不变,无论多少次的旋 ...

  7. 【OpenGL-矩阵】旋转矩阵推导

    参考博文 https://blog.csdn.net/lady_killer9/article/details/89306082 绕x轴逆时针旋转的矩阵推导过程 这是绕x轴逆时针旋转的矩阵形式: 在O ...

  8. 三维坐标旋转矩阵推导过程(包看懂)

    推导前提: 坐标使用右手坐标系,角度逆时针旋转为正. 绕X轴旋转角度为 俯仰角 即Pitch 绕Y轴旋转角度为 偏航角 即Yaw(Head) 绕Z轴旋转角度为 翻滚角 即Roll 一.平面二维坐标点的 ...

  9. 详解坐标变换矩阵 - 绕 x 轴旋转的旋转矩阵

    在高级驾驶辅助系统(ADAS)领域,存在多种常用的坐标系:LiDAR 坐标系.车辆坐标系.相机坐标系.图像坐标系等. 在高级驾驶辅助系统(ADAS)领域,存在多种常用的坐标系:LiDAR 坐标系.车辆 ...

最新文章

  1. Flask框架(flask中的邮件发送Flask-Mail(邮件扩展))
  2. 基于倒排表的电力调度知识问答系统构建 (100行代码带你实现nlp问答系统)
  3. JQUERY的parent()
  4. springMVC 的工作原理和机制
  5. HTML select控件
  6. 来了!汇总数据库行业最新招聘、热门资讯、精选文章...
  7. python中config命令_python的logging.config使用详解
  8. 硬核!如何全面系统地自学 Java ?
  9. 深入理解浏览器的缓存机制
  10. 服务器的固态硬盘使用raid非ssd,畅谈固态硬盘搭建RAID的方法
  11. sqldbx mysql 乱码_使用PLSQL Developer和DbVisualizer、SQLDBx查询oracle数据库时出现乱码...
  12. Golang之自旋锁
  13. python 异步爬取必应搜索结果
  14. No7.装饰器的进阶
  15. java项目——防止羊毛党“薅羊毛”
  16. 牛客大回馈,华为云优惠购买与使用心得
  17. 2021第十二届蓝桥杯省赛一等奖Java B组总结
  18. 2022年QIMA LVMH Ethical Audit Preparation Document
  19. JVM--类加载器详解
  20. 基于STM32设计的指针式电子钟与万年历

热门文章

  1. 文本识别论文之Aster解读
  2. [整理]Breakpoint on arbitrary selector
  3. 剪辑视频,如何提取视频封面
  4. LaTeX一行排两个图片示例
  5. ubuntu 安装启动 openssh-server
  6. 体验提升1:提升用户体验,应该从哪里开始?
  7. c语言枚举类型enum用法参数,C语言枚举类型(Enum)
  8. 除视频水印的软件 md5修改
  9. y7000p内存是一个16还是8+8_又来一颗新8核16线程处理器?5000左右有啥触控笔记本可以关注?...
  10. 互联网凛冬,看大厂HR怎么说~