基变换

理论部分

在n维的线性空间中,任意n个线性无关的向量都可以作为线性空间的基,即空间基不唯一。对于不同的基,同一个向量的坐标一般是不同的。因为在计算机图形学中,主要研究三维的空间,所以可以简化问题倒三维空间中的基变换,也就是坐标变换。

假设有两组单位正交基(图形学中选取的坐标系一般是正交的,为了方便)\(x,y,z\)和\(u,v,w\),他们的关系是

\[\begin{cases}
u = a_{11}x+a_{21}y+a_{31}z\\
v = a_{12}x+a_{22}y+a_{32}z\\
w = a_{13}x+a_{23}y+a_{33}z\\
\end{cases}
\]

因为都是基,所以可以用另一组基来表示。我们可以将其写作矩阵的形式

\[[u,v,w] = [x,y,z]
\begin{bmatrix}
a_{11} & a_{21} & a_{31}\\
a_{12} & a_{22} & a_{32}\\
a_{13} & a_{23} & a_{33}\\
\end{bmatrix}
\]

\[A = \begin{bmatrix}
a_{11} & a_{21} & a_{31}\\
a_{12} & a_{22} & a_{32}\\
a_{13} & a_{23} & a_{33}\\
\end{bmatrix}
\]

一般教科书里把A称作过渡矩阵。设一个向量\(V\),在基\((x,y,z)\)下的坐标为\((x_1,y_1,z_1)^{T}\),在基\((u,v,w)\)下的坐标为\((u_1,v_1,w_1)^{T}\),则有

\[\begin{aligned}
(x,y,z)(x_1,y_1,z_1)^{T} & =(u,v,w)(u_1,v_1,w_1)^{T} \\
&=(x,y,z)A(u_1,v_1,w_1)^{T}\\
\end{aligned}
\]

\[(x_1,y_1,z_1)^{T} = A(u_1,v_1,w_1)^{T} ......(1)
\]

因为\((x,y,z)\)为正交矩阵,有逆矩阵,所以可以消除。所以得到变换公式。

矩阵A的构成

上面说了一些理论,那么A矩阵是什么,怎么用呢?坐标值如果脱离基的话,没有任何意义。我们平时说的坐标,一般都有一个默认的条件的,就是位于原点的直角坐标系,笛卡尔坐标系。当我们在描述某个坐标系的时候,我们一般用笛卡尔坐标系作为参考来描述。也就是用笛卡尔坐标系里的坐标值来描述另一个坐标系的基。

在理论部分,如果\([x,y,z]\)代表的是笛卡尔坐标,那么\((a_{11},a_{21},a_{31})\)恰好就是\(u\)在笛卡尔坐标系的坐标表示,\(v,w\)同理,所以A矩阵相当于是由基\((u,v,w)\)在笛卡尔坐标系中的坐标组成的。记\((x_u,y_u,z_u)\)为\(u\)的坐标,\(v,w\)同理,那么矩阵A可以写作

\[A = \begin{bmatrix}
x_u & x_v & x_w\\
y_u & y_v & y_w\\
z_u & z_v & z_w\\
\end{bmatrix}
\]

由于A是正交矩阵,那么A的逆矩阵等于其转置矩阵,则由

\[A^{-1}(x_1,y_1,z_1)^{T} = (u_1,v_1,w_1)^{T} ......(2)
\]

\[A^{-1} = \begin{bmatrix}
x_u & y_u & z_u\\
x_v & y_v & z_v\\
x_w & y_w & z_w\\
\end{bmatrix}

\]

如果令\(u_0 = (x_u,y_u,z_u)^{T},v_0 = (x_v,y_v,z_v)^{T},w_0 = (x_w,y_w,z_w)^{T}\),那么(1)公式可以写作

\[(x_1,y_1,z_1)^{T}=
\begin{bmatrix}
u_0 &v_0&w_0
\end{bmatrix}
(u_1,v_1,w_1)^{T}
=u_1u_0+v_1v_0+w_1w_0
\]

就相当于是用\((u,v,w)\)基来表示了向量\(V\),而\((u,v,w)\)基使用笛卡尔坐标系来表示,所以\(V\)向量自然就可以转换到笛卡尔坐标系了。

记\(V_x = (x_1,y_1,z_1)^{T}\),(2)公式可以写作

\[(u_1,v_1,w_1)^{T}=
\begin{bmatrix}
u_0^{T} \\v_0^{T}\\w_0^{T}
\end{bmatrix}
V_x
=
\begin{bmatrix}
u_0^{T}V_x \\v_0^{T}V_x\\w_0^{T}V_X
\end{bmatrix}
\]

相当于是向基\((u,v,w)\)做投影,得到各个轴上的投影值,而投影值本质上就是这个坐标系内的坐标值。

综上所述,利用三维空间的坐标,阐述了坐标系的变换过程,可能有些不严谨,但是对于计算机图形学中的坐标系转换应该够了。

实际应用,世界坐标系->相机坐标系

图中\(O,x,y,z\)坐标系一般是世界坐标系,而\((e,u,v,w)\)是相机坐标系,模拟人眼看到的东西,物体的描述现在是在世界坐标系下,所以我们现在要将其描述转换到相机坐标系下,所有的物体都没有动,只是表示形式的变换。

定义相机坐标系

首先我们需要定义三个量

  • 眼睛的位置 \(e\),就是相机放在什么地方。
  • 视线朝向的向量 \(g\),相机看向什么地方,类似于眼睛往哪里看。
  • 朝向上方的向量 \(t\),相机顶部的朝向,类似于人头顶的朝向,一般都朝向天空。

以上的信息足够让我们定义出一个相机坐标系来

\[\begin{aligned}
w &= -\frac{g}{||g||}\\
u &= \frac{t\times w}{||t\times w||}\\
v &= w \times u
\end{aligned}
\]

  • 首先\(g\)是相机看过去的方向,而一般相机坐标系都是朝着\(w\)轴的负方向看去,因此对\(g\)取反并归一化,得到单位向量\(w\)。
  • 接着\(t\)类似于\(v\)的朝向,然后\(t\)和\(w\)做叉乘得到\(u\),因为叉乘得到的是\(t\)和\(w\)构成平面的法向量,然后归一化。
  • 同理\(w\)和\(u\)叉乘得到\(v\)

构建变换矩阵

现在\(e,u,v,w\)都有了,言外之意就是他们在世界坐标系中可以用坐标值来表示了。

现在要把物体在世界坐标系内的坐标转到相机坐标系中。首先假设,我们的相机坐标系位于原点,那么利用之前推导的公式,可以利用公式(2)进行转换

\[(u_1,v_1,w_1)^{T}=
\begin{bmatrix}
u_0^{T} \\v_0^{T}\\w_0^{T}
\end{bmatrix}
V_x
=
\begin{bmatrix}
u_0^{T}V_x \\v_0^{T}V_x\\w_0^{T}V_X
\end{bmatrix}
\]

如果相机坐标系位于原点,那么问题就结束了。但是相机坐标系位于\(e\)点,所以还存在一个位移量。那么是加还是减呢?是减,为什么?因为是相对于新坐标系。例如5相对于0的坐标是5,而距离3是2个单位,而5和3都是相对0来说的。因此用5-3得到5相对于3的坐标是2。因此我们应该用物体在世界坐标系下的坐标减去相机坐标系的原点坐标吗,得到物体相对于相机坐标系原点的坐标。

因此还得有一个位移矩阵

\[B=
\begin{bmatrix}
1&0&0&-x_e\\
0&1&0&-y_e\\
0&0&1&-z_e\\
0&0&0&1\\
\end{bmatrix}
\]

这里开始使用齐次坐标了,通常使用的都是齐次坐标,可以将位移可其他变换统一起来。

那么\(A^{-1}\)矩阵为

\[A^{-1} = \begin{bmatrix}
x_u & y_u & z_u & 0\\
x_v & y_v & z_v & 0\\
x_w & y_w & z_w & 0\\
0 & 0 & 0 & 1
\end{bmatrix}
\]

归纳总结

归纳总结一下

  1. 定义 \(e,g,t\),通过计算得到\(e,u,v,w\),用世界坐标系的坐标值表示
  2. 利用矩阵B进行位移
  3. 利用矩阵\(A^{-1}\)投影到相机坐标系的坐标轴上

\[u_{相机坐标系坐标值}=\begin{bmatrix}
x_u & y_u & z_u & 0\\
x_v & y_v & z_v & 0\\
x_w & y_w & z_w & 0\\
0 & 0 & 0 & 1
\end{bmatrix}
\begin{bmatrix}
1&0&0&-x_e\\
0&1&0&-y_e\\
0&0&1&-z_e\\
0&0&0&1\\
\end{bmatrix}
x_{世界坐标系坐标值}
\]

[computer graphics]世界坐标系->相机坐标系详细推导的更多相关文章

  1. Computer Graphics Research Software

    Computer Graphics Research Software Helping you avoid re-inventing the wheel since 2009! Last update ...

  2. 水题 HDOJ 4716 A Computer Graphics Problem

    题目传送门 /* 水题:看见x是十的倍数就简单了 */ #include <cstdio> #include <iostream> #include <algorithm ...

  3. Mathematics for Computer Graphics数学在计算机图形学中的应用 [转]

    最近严重感觉到数学知识的不足! http://bbs.gameres.com/showthread.asp?threadid=10509 [译]Mathematics for Computer Gra ...

  4. Computer Graphics Thinking–texture tiling

    Here is one question: how to tile texture? One thing worth to notice: The DirectX and OpenGL stipula ...

  5. HDU 4716 A Computer Graphics Problem

    A Computer Graphics Problem Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (J ...

  6. Mathematics for Computer Graphics

    Mathematics for Computer Graphics 最近严重感觉到数学知识的不足! http://bbs.gameres.com/showthread.asp?threadid=105 ...

  7. Vector Math for 3D Computer Graphics (Bradley Kjell 著)

    https://chortle.ccsu.edu/VectorLessons/index.html Chapter0 Points and Lines (已看) Chapter1 Vectors, P ...

  8. Fundamentals of Computer Graphics 中文版(第二版) (Peter Shirley 著)

    1 引言 2 数学知识 3 光栅算法 4 信号处理 5 线性代数 6 矩阵变换 7 观察 8 隐藏面消除 9 表面明暗处理 10 光线追踪 11 纹理映射 12 完整的图形流水线 13 图形学的数据结 ...

  9. HDU 4716 A Computer Graphics Problem (水题)

    A Computer Graphics Problem Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (J ...

随机推荐

  1. 发现用System.Net.Mail发邮件(代码附后),附件稍微大一点就会造成程序假死. 有没有什么简单的解决办法呢? 多谢!!

    附件大,上传,发送一定会慢.程序卡,应该是主线程正在发送,邮件造成的.创建其他线程在后台去发.这样就不影响主线程做其他工作了   using System; using System.Collecti ...

  2. Java并发编程volatile关键字

    volatile理解 Java语言是支持多线程的,为了解决线程并发的问题,在语言内部引入了 同步块 和volatile 关键字机制.volatile具有synchronized关键字的“可见性”,vo ...

  3. layui 数据表格使用

    //第一个实例 table.render({ elem: '#demo' ,height: 400 ,url: '/My_NewsTest/LookServlet' //数据接口 ,page: tru ...

  4. NO.1 MSP-ESP432P4111开箱

    本人准备2020TI杯模拟电子邀请赛,预计参赛可能会使用TI平台,故从某宝购置一块MSP-ESP432P4111 LaunchPad为参赛做准备.TI官网40美刀,但我只能找国内二道贩子买,有点小贵& ...

  5. DFS与DP算法

    名词解释: DFS(Dynamic Plan):动态规划 DFS(Depth First Search):深度优先搜索 DFS与DP的关系 很多情况下,dfs和dp两种解题方法的思路都是很相似的,这两 ...

  6. PAT 1032 Sharing (25分) 从自信到自闭

    题目 To store English words, one method is to use linked lists and store a word letter by letter. To s ...

  7. JavaScript几种继承方式的总结

    1.原型链继承 直接将子类型的原型指向父类型的实例,即"子类型.prototype = new 父类型();",实现方法如下: //父类构造函数 function father(n ...

  8. spring框架中三层架构相关的注解

    做了这么多年的C++,再去学Java,确实发现,语言都是相通的,即使是Java的那么多生态,理解起来也并不费劲 Spring 框架目前还在学习中,处于 Tourist 阶段,目前只求会做,不求原理,等 ...

  9. 01 . 分布式存储之FastDFS简介及部署

    分布式存储简介 现代的互联网已经进入大数据时代,每天都有数以万计的数据产生,这些数据的规模轻轻松松地可以达到几P的级别,传统的的单机存储早已捉襟见肘,根本无法满足大数据对存储系统的要求.这时,各种分布 ...

  10. Spring boot Sample 003之spring-boot-configuration-properties

    一.环境 1.1.Idea 2020.1 1.2.JDK 1.8 二.目的 通过properties文件配置spring boot 属性文件 三.步骤 3.1.点击File -> New Pro ...