Metal 透视投影

引言

Hi，大家好，我是一牛。在上一篇博客中，我们学习了如何在Metal 中使用MVP变换实现正交投影。正交投影常用于工程制图，它可以保持物体的比例不变，但是如果我们想要模拟人眼中的现实世界（近大远小），我们需要使用透视投影。今天，让我们一起学习下透视投影是如何实现的。

视椎体

透视投影矩阵

和正交投影使用长方体不同的是，在透视投影中我们使用了视椎体，经过模型变换、视图变换后，为了得到透视投影矩阵，我们只需要将是椎体压缩成长方体，之后再进行一次正交投影，我们就能够得到透视矩阵。

那么我们是如何通过挤压视椎体得到正交投影需要要的长方体？

假设我们从+x 方向观察，看向-x

根据相似三角形我们可以得到

$y' = \frac{ny}{z}$

同理，我们可以得到

$x' = \frac{nx}{z}$

对于空间中的任何一点 $\begin{pmatrix} x \\ y \\ z \\ 1 \end{pmatrix}$

经过挤压得到 $\begin{pmatrix} \frac{nx}{z} \\ \frac{ny}{z} \\ unknown \\ 1 \end{pmatrix}$

由于是齐次坐标，我们可以将乘z 得到 $\begin{pmatrix} nx \\ ny \\ unknown \\ z \end{pmatrix}$

我们可以初步得到这个矩阵

$\begin{pmatrix} n & 0 & 0 & 0 \\ 0 & n & 0 & 0 \\ A & B & C & D \\ 0 & 0 & 1 & 0\end{pmatrix}$

接下来我们只需要求出待定系数 A,B,C,D

我们知道，近平面上的点经过挤压保持不变。也就是说对于近平面任意一点 $\begin{pmatrix} x \\ y \\ n \\ 1 \end{pmatrix}$

由于它是齐次坐标，我们可以乘上n，得到$\begin{pmatrix} nx \\ ny \\ n^2 \\ n \end{pmatrix}$

我们观察待求矩阵的第三方，可得 A=0,B=0,即

$\begin{bmatrix} 0 & 0 & A & B \end{bmatrix}$ $\begin{pmatrix} x \\ y \\ n \\ 1 \end{pmatrix}$ = $n^2$

-> $A * n + B = n^2$ (1)

又因为远平面的中心经过挤压保持不变，所以对应远平面中心点 $\begin{pmatrix} 0 \\ 0 \\ f \\ 1 \end{pmatrix}$

由于它是齐次坐标，我们可以乘上f，得到$\begin{pmatrix} 0 \\ 0 \\ f^2 \\ f \end{pmatrix}$

-> $\begin{bmatrix} 0 & 0 & A & B \end{bmatrix}$ $\begin{pmatrix} 0 \\ 0 \\ f \\ 1 \end{pmatrix}$ = $f^2$

-> $A * f + B = f^2$ (2)

联立（1）（2）解一元二次方程组

得 $A = n + f$ 、$B = -n*f$

带人A,B,C,D 系数，得到挤压矩阵 $\begin{pmatrix} n & 0 & 0 & 0 \\ 0 & n & 0 & 0 \\ 0 & 0 & n + f & -n*f \\ 0 & 0 & 1 & 0\end{pmatrix}$

将挤压矩阵左乘正交投影矩阵 $\begin{pmatrix} \frac{2}{r-l} & 0 & 0 & -\frac{r+l}{r-l} \\ 0 & \frac{2}{t-b} & 0 & -\frac{t+b}{t-b} \\ 0 & 0 & \frac{1}{f-n} & -\frac{n}{f-n} \\ 0 & 0 & 1 & 0\end{pmatrix}$

最终透视投影矩阵是

$\begin{pmatrix} \frac{2n}{r-l} & 0 & -\frac{r+l}{r-l} & 0 \\ 0 & \frac{2n}{t-b} & -\frac{t+b}{t-b} & 0 \\ 0 & 0 & \frac{f}{f-n} & -\frac{n*f}{f-n} \\ 0 & 0 & 1 & 0\end{pmatrix}$

在图形学中我们一般用视场角fovY,和近平面宽高比来定义透视投影矩阵。假定视椎体是对称物体，$b=-t$ 且 $l=-r$

$tan fovY / 2 = t / n$ 且 $aspectRatio = r / t$

我们可以得到透视投影矩阵

    func createPerspectiveMatrix(fov: Float, aspectRatio: Float, nearPlane: Float, farPlane: Float) -> simd_float4x4 {
        let tanHalfFov = tan(fov / 2.0);
        var matrix = simd_float4x4(0.0);
        matrix[0][0] = 1.0 / (aspectRatio * tanHalfFov);
        matrix[1][1] = 1.0 / (tanHalfFov);
        matrix[2][2] = farPlane / (farPlane - nearPlane);
        matrix[2][3] = 1.0;
        matrix[3][2] = -(farPlane * nearPlane) / (farPlane - nearPlane);
        return matrix;
    }

效果

仔细观察，与正交投影变换不同的是，图片的比例在转动的过程中会发生改变，这是因为转动过程中图片的深度值在发生改变，对应不同的深度值，透视投影会产生近大远小的效果！

结语

掌握好透视投影是我们打开3D世界的钥匙，而掌握好透视投影的关键是了解透视投影的矩阵推导过程。 谢谢大家，欢迎大家点赞、收藏！

本项目已开源