计算机图形学 变换矩阵
为了完成我的 计算机图形学作业,有一道变换矩阵的问题,翻了好久,发现这篇博客讲的很清楚,感谢大佬的分享。
以下内容转自:http://m.blog.csdn.net/csxiaoshui/article/details/65446125
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+txy′=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旋转的旋转矩阵是【因为使用的列向量,因此执行的是左乘(从右往左)】:
(注:式中的(u,v,w)对应上文中向量(a,b,c),公式我自己笔算过,为了减少编辑公式的时间(使用LaTex编辑太繁琐,因此找了一张公式的图片贴在此处)
如果向量是经过单位化的(单位向量),那么有a2+b2+c2=1,可以简化上述的公式,得到:
参考文献:
- Wiki Rotation (mathematics)
- Euler’s rotation theorem
- Maths - Rotation Matrices
- 绕任意轴旋转
- Rotation About an Arbitrary Axis in 3 Dimensions
- Rotation about an Arbitrary Axis (Line)
计算机图形学 变换矩阵相关推荐
- 第一个使用计算机图形学术语,计算机图形学考题
下面哪个函数是OpenGl的工具函数( ) A.glBegin B.glEnd C. glClear D. gluLookAt Breseham算法直线算法中,当误差项d 0.5时,直线与列垂直网格的 ...
- java 矩阵_(JAVA)计算机图形学与矩阵
(JAVA)计算机图形学与矩阵 平移 旋转 变比 一)矩阵 1)矩阵的乘法运算 矩阵是计算机图形学的重要工具. 计算机图形学要常常用到 矩阵乘法(矩阵与矩阵相乘). 矩阵(Aij)m*n 与矩阵(Bi ...
- 3、计算机图形学——模型视图变换、投影变换与视口变换
一.模型视图变换 模型视图变换主要是为了让摄像机回归到世界坐标的原点并且和拍摄物体一起进行变换,便于计算 模型视图变换的根据就是物体和相机的相对位置不变,那么,投影得到的图片也是不变的 首先规定相机拍 ...
- 1、计算机图形学——2D变换与齐次坐标
一.2D变换 1.缩放变换 缩放变换就是将图形上的点的各个维度同时放大或者缩小,比如,讲一个图片的长宽都缩小0.5倍 则缩小之后,新的图片上的点(x', y')可以表示为 用写成矩阵的形式就是 如果长 ...
- 计算机图形学三大应用领域,计算机图形学作业
计算机图形学作业题 1. 计算机中由图形的形状参数(方程或分析表达式的系数,线段的端点坐标等)加属性参数(颜色.线型等)来表示图形称图形的参数表示:枚举出图形中所有的点称图形的点阵表示,简称为图像(数 ...
- 用C#实现计算机图形学算法
多数情况下计算机图形学算法都用C++实现,下面鄙人用C#实现一部分算法.并附上运行截图. 一 图案 1 金刚石 金刚石图案是每一个顶点都与其他顶点相连的正n边形.金刚石图案有时被用作计算机图形设备的测 ...
- 对计算机图形学知识总结第二版,计算机图形学入门总结!
图形学实验在百度网盘. 1. Why study transformation 1.1 Modeling translation (平移) rotation(旋转) scaling (缩放) proj ...
- 计算机图形学基础考试题,计算机图形学基础复习题
<计算机图形学基础复习题>由会员分享,可在线阅读,更多相关<计算机图形学基础复习题(8页珍藏版)>请在人人文库网上搜索. 1.计算机图形学基础复习题 一.判断题 1. PNG( ...
- 计算机有哪些应用领域 请举例说明,计算机图形学作业求助
计算机图形学作业求助 计算机图形学作业题 计算机中由图形的形状参数(方程或分析表达式的系数,线段的端点坐标等)加属性参数(颜色.线型等)来表示图形称图形的参数表示:枚举出图形中所有的点称图形的点阵表示 ...
最新文章
- 简历包装要避开哪些坑,资深面试官告诉你!
- android焦距跟实际距离,手机等效焦距一算便知_nubia Z5S Mini_手机Android频道-中关村在线...
- begin again
- php xml常用函数的集合及四种方法
- visio中公式太小_visio绘图中的数据计算
- 【Jmeter篇】导出Jmeter脚本及查看接口响应时间
- 分布式文件系统-HDFS( HDFS全称是Hadoop Distributed System)
- CCNP-EIGRP路由认证(附实验PDF下载)
- Java 的强引用、弱引用、软引用、虚引用
- win7 下 qwt安装教程
- logo设计好坏的12个标准
- 网易编程题——牛牛的闹钟
- x550网卡linux驱动,Intel英特尔X520/X540/X550/82599系列万兆网卡驱动5.10.2版For Linux(2021年2月1日发布)...
- [从零开始学习FPGA编程-22]:进阶篇 - 架构 - FPGA内部硬件电路的设计与建模
- 一种串口扩展电路应用
- 怎么拍照识别植物?手把手教你识别
- python小项目——2048小游戏(详解)
- ValueError: Unknown activation function: ReLU
- 【密码学|算法设计】拓展的欧几里得算法及理论证明 (Extended Euclidean Algorithm)
- c语言设计篮球计分器软件,会玩篮球更会玩篮球赛计时计分器,完整设计资料开源...
热门文章
- NVIDIA显卡虚拟化vGPU终于支持KVM了
- 史蒂芬·乔布斯amp;斯坦福大学
- 设置mysql中自增列的初始值和增加步长
- sql语句往某个字段指定位置追加或者插入值
- hfss和python_利用Python与HFSS联合仿真设计一个微带天线!
- $forceUpdate的使用详解
- Kaggle淋巴结病理切片有无癌细胞鉴别建模:Logistic+SVM+RandomForest+CNN
- 主流Webrtc流媒体服务器之Kurento Media Server
- 解决Electron Fiddle无法下载安装electron不同版本支持的问题 - Fiddle一直都是checking status
- JAVA基础(16) Servlet/Filter/Listener/Interceptor区别与联系