模型变换

错切变换

缩放变换

旋转变换

平移变换

绕任意向量的旋转矩阵公式

$R(\mathbf n, \alpha)=cos(\alpha) I + (1-cos(\alpha)) \mathbf n \cdot \mathbf n ^ T + sin(\alpha) \begin{bmatrix} 0 & -n_z & n_y \\ n_z & 0 & -n_x\\ -n_y & n_x & 0 \end{bmatrix}$

推导过程

相机变换

相机变换矩阵的推导

一个相机由以下几部分所构成

相机位置 $\vec e$
相机视线方向 $\vec g$
相机的正上方向 $\vec t$
在相机的参考系中，定义相机视线方向 $\vec g$ 为相机坐标系的z轴方向
定义相机的正上方向 $\vec t$ 为相机坐标系的y轴方向
定义 $\vec g \times \vec t$ 为相机坐标系的x轴方向
通常相机的视线方向为世界坐标系的 $-z$ 轴方向

首先将相机平移到原点，可得到矩阵

$T_{view}=\begin{bmatrix} 1 & 0 & 0 & -x_e\\ 0 & 1 & 0 & -y_e\\ 0 & 0 & 1 & -z_e\\ 0 & 0 & 0 & 1 \end{bmatrix}$

将 $\vec g$ 旋转到 $-z$ 方向，将 $\vec t$ 旋转到y轴方向，将 $\vec g \times \vec t$ 旋转到x轴方向，可得到矩阵 $R_{view}$ ，但是这个矩阵并不好求，但我们可以反过来求。
我们可以按照变换后的基底直接求出逆变换的变换矩阵，由于是正交基底构成的矩阵，所以是正交矩阵，所以 $R_{view}^T=R^{-1}_{view}$

$R^{-1}_{view}=\begin{bmatrix} x_{\vec g \times \vec t} & x_{\vec t} & x_{\vec {-g}} & 0\\ y_{\vec g \times \vec t} & y_{\vec t} & y_{\vec {-g}} & 0\\ z_{\vec g \times \vec t} & z_{\vec t} & z_{\vec {-g}} & 0\\ 0 & 0 & 0 & 1 \end{bmatrix}$

故最终相机变换矩阵得到

$M_{view}=R_{view}T_{view}$

投影变换

平行投影变换矩阵推导

假设平行投影摄像机的视线范围由长方体 $[l, r] \times [b,t] \times [f,n]$ 所确定，其中l,r为x轴方向的两平面位置，b,t为y轴方向的两平面位置，f,n为z轴方向的两平面位置。我们希望将其标准化到 $[-1,1]^3$
首先将其平移到原点，然后再缩放它。

$M_{ortho}= \begin{bmatrix} \frac{2}{r-l} & 0 & 0 & 0 \\ 0 & \frac{2}{t-b} & 0 & 0 \\ 0 & 0 & \frac{2}{n-f} & 0 \\ 0 & 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} 1 & 0 & 0 & -\frac{r+l}{2} \\ 0 & 1 & 0 & -\frac{t+b}{2} \\ 0 & 0 & 1 & -\frac{n+f}{2} \\ 0 & 0 & 0 & 1 \end{bmatrix}$

透视投影变换矩阵推导

透视投影摄像机的视线范围由一个平头四棱台体所确定，使得内部的物体可以呈现近大远小的效果。
我们可以想象以下直接通过“挤压”将平头四棱台挤压为一个长方体，变换前后n平面上的点不变，n~f平面上的点将被压缩，越接近f平面，点被压缩的越厉害，所以当变换到长方体后相同的三维几何体离f平面越近投影后的图像就越小，故能够呈现出近大远小的效果。这也是一种空间变换，但它在三维空间中属于一种非线性变换，如果我们希望使用矩阵描述，则需要借助齐次坐标将其放置在四维空间中完成这个变换。
假设这种“挤压”变换可以使用矩阵A描述，则最终的透视投影变换矩阵可表述如下：

$M_{persp}=M_{ortho} \cdot M_{persp \to ortho}$

矩阵 $M_{persp \to ortho}$ 的推导如下：
$x'=\frac{n}{z}x$
$y'=\frac{n}{z}y$

$\begin{bmatrix} x \\ y \\ z \\ 1 \\ \end{bmatrix} \Rightarrow \begin{bmatrix} \frac{n}{z}x \\ \frac{n}{z}y \\ ? \\ 1 \\ \end{bmatrix} = \begin{bmatrix} nx \\ ny \\ ? \\ z \\ \end{bmatrix}$

由于n平面上的任何点在变换前后位置都仍然不变，设未知量?为k可得出以下关系：

$\begin{bmatrix} x \\ y \\ n \\ 1 \\ \end{bmatrix} = \begin{bmatrix} nx \\ ny \\ k \\ n \\ \end{bmatrix}$

解得 $k=n^2$
故矩阵 $M_{persp \to ortho}$ 可初步写为如下形式

$M_{persp \to ortho}= \begin{bmatrix} n & 0 & 0 & 0\\ 0 & n & 0 & 0\\ 0 & 0 & A & B\\ 0 & 0 & 1 & 0 \end{bmatrix}$

由于f平面和z轴的交点在变换前后都不变，故得出以下关系：

$\begin{bmatrix} 0 \\ 0 \\ f \\ 1 \\ \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \\ f^2 \\ f \\ \end{bmatrix}$

由这两个特殊点可列出关于A，B未知数的二元一次方程组

$\left\{\begin{matrix} An+B=n^2 \\ Af+B=f^2 \end{matrix}\right.$

解得

$\left\{\begin{matrix} A=n+f \\ B=-nf \end{matrix}\right.$

故得到矩阵

$M_{persp \to ortho}= \begin{bmatrix} n & 0 & 0 & 0\\ 0 & n & 0 & 0\\ 0 & 0 & n+f & -nf\\ 0 & 0 & 1 & 0 \end{bmatrix}$

深入思考

当我们从透视投影的平头四棱台空间压缩变换到平行投影的长方体后，对于内部的任意一点 $(x,y,z)$ 坐标，变换后对应的坐标为 $(x',y',z')$ ，显然是可以得到 $z$ 越远的点， $(x',y')$ 坐标被压缩的越厉害，关系式可在上述推导中体现。但是内部任意一个点的z坐标在变换前后就没有那么直观，我们可以进一步做一系列推导。

$\begin{bmatrix} x' \\ y' \\ z' \\ 1 \end{bmatrix} = \begin{bmatrix} n & 0 & 0 & 0\\ 0 & n & 0 & 0\\ 0 & 0 & n+f & -nf\\ 0 & 0 & 1 & 0 \end{bmatrix} \begin{bmatrix} x \\ y \\ z \\ 1 \end{bmatrix} = \begin{bmatrix} nx \\ ny \\ (n+f)z-nf \\ z \end{bmatrix} = \begin{bmatrix} n\frac{x}{z} \\ n\frac{y}{z} \\ (n+f)-\frac{nf}{z} \\ 1 \end{bmatrix}$

可得到函数 $z'=(n+f)-\frac{nf}{z}$ 关于自变量 $z$ 的函数关系， $z$ 的定义域为 $[n,f]$
为了得知变换前后的z坐标是变远了还是变近了，可以做差得到函数 $f(z)$

$f(z)=z'-z=(n+f)-(\frac{nf}{z}+z), z \in [n,f] > 0$

令 $f(z)=0$ ，则得到方程 $z^2-(n+f)+nf=0$ ，即 $(z-n)(z-f)=0$ ，得到方程有两个零点

$x_1=n, x_2=f$

函数 $f(z)$ 求导得 $f'(z)=\frac{nf}{z^2}-1$
令 $f'(z)>0$ 得， $n<z<\sqrt{nf}$ ，令 $f'(z)<0$ 得， $\sqrt{nf}<z<f$
得出以下函数及导函数分布表

z	$n$	$(n,\sqrt{nf})$	$\sqrt{nf}$	$(\sqrt{nf},f)$	$f$
f’(z)	-	>0	0	<0	-
f(z)	0	增函数	-	减函数	0

故 $f(z)>=0$ 恒成立（当 $z=n \space or \space z=f$ 时，等号成立）
故对于透视投影空间内的某一点被压缩后， $z'>z$ ，z坐标变远。