基变换

理论部分

在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. Django项目文件

    Django项目文件 Pycharm创建django项目: 创建一个有一个app01应用的jangotest项目: 创建步骤: pycharm左上角文件file①---->新建项目new pro ...

  2. 迟早要知道的JS系列之常用数组方法

    常用数组方法 一.不会改变原始数组的方法: 即访问方法,下面的这些方法绝对不会改变调用它们的对象的值,只会返回一个新的数组或者返回一个其它的期望值. 1. concat() ** 语法:** Java ...

  3. python3.x 基础四:json与pickple

    每次打开一个文件,只dump1次 json.dump(dump的内容,文件句柄) json.load(文件句柄) json可以处理列表/字典/字符串等简单数据类型,但是不能处理复杂的数据类型,如函数的 ...

  4. POJ1436

    题目链接:https://vjudge.net/problem/POJ-1436 解题思路:基于y轴建立线段树. 如图是根据样例画出的图.下面都以题目样例为例. 但是,如果仅仅以给出的y1, y2为边 ...

  5. Spring MVC必须知道的执行流程

    Spring MVC的执行流程 一.名词解释 1.前端控制器(DispatcherServlet) 接收请求,响应结果,相当于转发器,中央处理器 2.处理器映射器(HandlerMapping) 根据 ...

  6. 实现.Net程序中OpenTracing采样和上报配置的自动更新

    前言 OpenTracing是一个链路跟踪的开放协议,已经有开源的.net实现:opentracing-csharp,同时支持.net framework和.net core,Github地址:htt ...

  7. 包装类的使用与Junit单元测试类

    包装类: 针对八种基本数据类型定义相应的引用类型,使之有了类的特点,就可以调用类的方法 基本数据类型 包装类 boolean Boolean byte Byte short Short int Int ...

  8. idea本地Maven仓库不能下载依赖jar包的解决方案

    1.确认maven是否正正常安装,是否配置了环境变量,可以通过命令 mvn -version 看是否显示maven的版本信息. 2.检查maven的setting.xml配置文件中本地仓库位置配置是否 ...

  9. vue 中引入使用jquery

    1.首先在配置文件中添加 cnpm install  下载jquery文件 2.在webpack配置文件中添加下面代码 3.接着在webpack.base.conf.js中module.exports ...

  10. Rocket - debug - Example: Read Memory

    https://mp.weixin.qq.com/s/ChXNTbx94WDC72GvmE9bGA 介绍riscv-debug的使用实例:使用三种方法读取内存. 1. Using System Bus ...