[ x y ] [ a c b d ] \begin{bmatrix}
x & y
\end{bmatrix}
\begin{bmatrix}
a & c \\
b & d
\end{bmatrix} [ x y ] [ a b c d ]
矩阵乘法
首先复习一下线性代数,矩阵相乘,结果如下 ↓
[ a x + b y c x + d y ] \begin{bmatrix}
ax+by & cx+dy
\end{bmatrix} [ a x + b y c x + d y ]
用文字大概可以描述为 A 的第 n
行和 B 的第 m
列运算(相乘再相加)得出结果矩阵的 (n,m)
的值。
感觉也可以描述为一个值与自己同一行的值有一定关系,最后将这些关系映射成新的值。看不懂的话还是先看教科书吧,我也不懂描述数学 🤣,顺便推一个知乎链接:矩阵乘法的本质是什么?
缩放矩阵
在明白矩阵乘法的基础上,可以尝试从结果反推出矩阵。
假设( x , y ) (x,y) ( x , y ) 是一个点的坐标,要以( 0 , 0 ) (0,0) ( 0 , 0 ) 为中心缩放这个点,自然是( x s x , y s y ) (x_sx,y_sy) ( x s x , y s y ) ,要得出这个结果,我们改一下上面的例子就能得出:
[ x y ] [ x s 0 0 y s ] \begin{bmatrix}
x & y
\end{bmatrix}
\begin{bmatrix}
x_s & 0 \\
0 & y_s
\end{bmatrix} [ x y ] [ x s 0 0 y s ]
平移矩阵
要平移这个点自然是( x + x t , y + y t ) (x+x_t,y+y_t) ( x + x t , y + y t ) ,但是只以上面为基础是改不出来的,因为相加的另一项总与另一个轴有关,并不是一个常数,怎么办呢……
想不出来不如先看看下面这个式子:
[ x y o ] [ a d g b e h c f i ] \begin{bmatrix}
x & y & o
\end{bmatrix}
\begin{bmatrix}
a & d & g \\
b & e & h \\
c & f & i \\
\end{bmatrix} [ x y o ] ⎣ ⎡ a b c d e f g h i ⎦ ⎤
结果是( a x + b y + o c , d x + e y + f o , g o + h o + i o ) (ax+by+oc,dx+ey+fo,go+ho+io) ( a x + b y + oc , d x + ey + f o , g o + h o + i o ) ,当 o 是 1 的时候,就是( a x + b y + c , d x + e y + f , g + h + i ) (ax+by+c,dx+ey+f,g+h+i) ( a x + b y + c , d x + ey + f , g + h + i ) 不就每一项都能加个常量了。所以加常量的矩阵长这样:
[ x y 1 ] [ 1 0 0 0 1 0 x t y t 1 ] \begin{bmatrix}
x & y & 1
\end{bmatrix}
\begin{bmatrix}
1 & 0 & 0 \\
0 & 1 & 0 \\
x_t & y_t & 1 \\
\end{bmatrix} [ x y 1 ] ⎣ ⎡ 1 0 x t 0 1 y t 0 0 1 ⎦ ⎤
为了处理更方便可以把缩放的矩阵也扩展为 3×3 矩阵:
[ x y 1 ] [ x s 0 0 0 y s 0 0 0 1 ] \begin{bmatrix}
x & y & 1
\end{bmatrix}
\begin{bmatrix}
x_s & 0 & 0 \\
0 & y_s & 0 \\
0 & 0 & 1
\end{bmatrix} [ x y 1 ] ⎣ ⎡ x s 0 0 0 y s 0 0 0 1 ⎦ ⎤
旋转矩阵
旋转比上面两种转换要难一点,首先复习一下两角和 公式:
s i n ( α + β ) = s i n α c o s β + c o s α s i n β c o s ( α + β ) = c o s α c o s β − s i n α s i n β sin(α+β)=sinαcosβ+cosαsinβ \\
cos(α+β)=cosαcosβ-sinαsinβ s in ( α + β ) = s in α cos β + cos α s in β cos ( α + β ) = cos α cos β − s in α s in β
再借极坐标 转换一下:
x = r c o s α y = r s i n α x r = r c o s ( α + β ) = r c o s α c o s β − r s i n α s i n β = x c o s β − y s i n β y r = r s i n ( α + β ) = r s i n α c o s β + r c o s α s i n β = x s i n β + y c o s β x=rcosα \\
y=rsinα \\
x_r=rcos(α+β)=rcosαcosβ-rsinαsinβ=xcosβ-ysinβ \\
y_r=rsin(α+β)=rsinαcosβ+rcosαsinβ=xsinβ+ycosβ x = rcos α y = rs in α x r = rcos ( α + β ) = rcos α cos β − rs in α s in β = x cos β − ys in β y r = rs in ( α + β ) = rs in α cos β + rcos α s in β = x s in β + ycos β
这不就是上面熟悉的形式,立即造个矩阵处理旋转:
[ x y 1 ] [ c o s β s i n β 0 − s i n β c o s β 0 0 0 1 ] \begin{bmatrix}
x & y & 1
\end{bmatrix}
\begin{bmatrix}
cosβ & sinβ & 0 \\
-sinβ & cosβ & 0 \\
0 & 0 & 1 \\
\end{bmatrix} [ x y 1 ] ⎣ ⎡ cos β − s in β 0 s in β cos β 0 0 0 1 ⎦ ⎤
约定
[ 1 0 x t 0 1 y t 0 0 1 ] [ x s 0 0 0 y s 0 0 0 1 ] [ c o s β − s i n β 0 s i n β c o s β 0 0 0 1 ] [ x y 1 ] \begin{bmatrix}
1 & 0 & x_t \\
0 & 1 & y_t \\
0 & 0 & 1 \\
\end{bmatrix}
\begin{bmatrix}
x_s & 0 & 0 \\
0 & y_s & 0 \\
0 & 0 & 1
\end{bmatrix}
\begin{bmatrix}
cosβ & -sinβ & 0 \\
sinβ & cosβ & 0 \\
0 & 0 & 1 \\
\end{bmatrix}
\begin{bmatrix}
x \\
y \\
1
\end{bmatrix} ⎣ ⎡ 1 0 0 0 1 0 x t y t 1 ⎦ ⎤ ⎣ ⎡ x s 0 0 0 y s 0 0 0 1 ⎦ ⎤ ⎣ ⎡ cos β s in β 0 − s in β cos β 0 0 0 1 ⎦ ⎤ ⎣ ⎡ x y 1 ⎦ ⎤
很多文章大家都把 xyz 写成 3 行而不是 3 列,而且是最后乘 xyz 的,只是我一开始觉得横着写比较好理解就那么整了,有点误导人,所以这里补一下约定的写法,不然大家看到这里的矩阵跟其他地方看到的不一样可能会感到疑惑。
实战
虽然上面只是说明了某一个点的变换,但是一个图形或图像就是点的集合,对这些点都应用相同的矩阵就能对图片进行变换。
最简单的矩阵变换实战就藏在 CSS 里,就是那个 transform
的 matrix
,是个处理 2D 图像的矩阵。
2D
matrix
的取值有一点点奇葩,只有 6 个参数,我们需要把 matrix
的 6 个参数按下面位置填入矩阵,不过也确实,最后一行没有必要修改。
m a t r i x ( a , b , c , d , e , f ) [ a c e b d f 0 0 1 ] matrix(a,b,c,d,e,f) \\
\begin{bmatrix}
a & c & e \\
b & d & f \\
0 & 0 & 1
\end{bmatrix} ma t r i x ( a , b , c , d , e , f ) ⎣ ⎡ a b 0 c d 0 e f 1 ⎦ ⎤
例如默认的矩阵就是 a=d=1
,也就是 transform: matrix(1, 0, 0, 1, 0, 0);
。
对于多次平移、缩放、旋转的组合操作,可以将这些操作矩阵都相乘再与 xy 求值。
例如“放大 3 倍”、“下移 45”、“右移 15”、“旋转 30°”就是:
[ 3 0 0 0 3 0 0 0 1 ] [ 1 0 15 0 1 45 0 0 1 ] [ c o s ( π 6 ) − s i n ( π 6 ) 0 s i n ( π 6 ) c o s ( π 6 ) 0 0 0 1 ] \begin{bmatrix}
3 & 0 & 0 \\
0 & 3 & 0 \\
0 & 0 & 1
\end{bmatrix}
\begin{bmatrix}
1 & 0 & 15 \\
0 & 1 & 45 \\
0 & 0 & 1 \\
\end{bmatrix}
\begin{bmatrix}
cos(\frac{π}{6}) & -sin(\frac{π}{6}) & 0 \\
sin(\frac{π}{6}) & cos(\frac{π}{6}) & 0 \\
0 & 0 & 1 \\
\end{bmatrix} ⎣ ⎡ 3 0 0 0 3 0 0 0 1 ⎦ ⎤ ⎣ ⎡ 1 0 0 0 1 0 15 45 1 ⎦ ⎤ ⎣ ⎡ cos ( 6 π ) s in ( 6 π ) 0 − s in ( 6 π ) cos ( 6 π ) 0 0 0 1 ⎦ ⎤
算出来大概是这样(注意矩阵乘法不符合交换律,所以顺序改变得出的结果是不同的):
( 2.598075 − 1.5 45 1.5 2.598075 135 0 0 1 ) \begin{pmatrix}2.598075&-1.5&45\\ 1.5&2.598075&135\\ 0&0&1\end{pmatrix} ⎝ ⎛ 2.598075 1.5 0 − 1.5 2.598075 0 45 135 1 ⎠ ⎞
也就是 matrix(2.6,1.5,-1.5,2.6,45,135)
,可以看出还是有一定规律的(除非用了 skew)。
3D
虽然上面并没有提到 3D 的变换矩阵,但是在 2D 基础上理解并不难,尤其是平移和缩放:
[ 1 0 0 x t 0 1 0 y t 0 0 1 z t 0 0 0 1 ] [ x s 0 0 0 0 y s 0 0 0 0 z s 0 0 0 0 1 ] \begin{bmatrix}
1 & 0 & 0 & x_t \\
0 & 1 & 0 & y_t \\
0 & 0 & 1 & z_t \\
0 & 0 & 0 & 1 \\
\end{bmatrix}
\begin{bmatrix}
x_s & 0 & 0 & 0 \\
0 & y_s & 0 & 0 \\
0 & 0 & z_s & 0 \\
0 & 0 & 0 & 1
\end{bmatrix} ⎣ ⎡ 1 0 0 0 0 1 0 0 0 0 1 0 x t y t z t 1 ⎦ ⎤ ⎣ ⎡ x s 0 0 0 0 y s 0 0 0 0 z s 0 0 0 0 1 ⎦ ⎤
旋转是依然复杂,绕 3 个轴的旋转需要使用 3 个矩阵处理(下面的 s 和 c 代表 sin 和 cos):
[ c − s 0 0 s c 0 0 0 0 1 0 0 0 0 1 ] [ c 0 s 0 0 1 0 0 − s 0 c 0 0 0 0 1 ] [ 1 0 0 0 0 c − s 0 0 s c 0 0 0 0 1 ] \begin{bmatrix}
c & -s & 0 & 0 \\
s & c & 0 & 0 \\
0 & 0 & 1 & 0 \\
0 & 0 & 0 & 1
\end{bmatrix}
\begin{bmatrix}
c & 0 & s & 0 \\
0 & 1 & 0 & 0 \\
-s & 0 & c & 0 \\
0 & 0 & 0 & 1
\end{bmatrix}
\begin{bmatrix}
1 & 0 & 0 & 0 \\
0 & c & -s & 0 \\
0 & s & c & 0 \\
0 & 0 & 0 & 1
\end{bmatrix} ⎣ ⎡ c s 0 0 − s c 0 0 0 0 1 0 0 0 0 1 ⎦ ⎤ ⎣ ⎡ c 0 − s 0 0 1 0 0 s 0 c 0 0 0 0 1 ⎦ ⎤ ⎣ ⎡ 1 0 0 0 0 c s 0 0 − s c 0 0 0 0 1 ⎦ ⎤
第一个矩阵十分熟悉,其实就是 2D 多一个维度,可以说 2D 的旋转其实就是绕 Z 轴旋转,第三个 X 相关的数据是不变的,是绕 X 轴旋转。第二个 Y 相关的数据不变,就是绕 Y 轴旋转,不过需要注意绕 Y 旋转的 sin 正负有点不同。
回到 CSS 的 matrix3d
,它的参数是完整的 4×4 矩阵,计算方式跟 2D 是一样的,就不再赘述了:
m a t r i x 3 d ( a 1 , b 1 , c 1 , d 1 , a 2 , b 2 , c 2 , d 2 , a 3 , b 3 , c 3 , d 3 , a 4 , b 4 , c 4 , d 4 ) [ a 1 a 2 a 3 a 4 b 1 b 2 b 3 b 4 c 1 c 2 c 3 c 4 d 1 d 2 d 3 d 4 ] matrix3d(a1, b1, c1, d1, a2, b2, c2, d2, a3, b3, c3, d3, a4, b4, c4, d4) \\
\begin{bmatrix}
a1&a2&a3&a4 \\
b1&b2&b3&b4 \\
c1&c2&c3&c4 \\
d1&d2&d3&d4
\end{bmatrix} ma t r i x 3 d ( a 1 , b 1 , c 1 , d 1 , a 2 , b 2 , c 2 , d 2 , a 3 , b 3 , c 3 , d 3 , a 4 , b 4 , c 4 , d 4 ) ⎣ ⎡ a 1 b 1 c 1 d 1 a 2 b 2 c 2 d 2 a 3 b 3 c 3 d 3 a 4 b 4 c 4 d 4 ⎦ ⎤
都拓展到 3D 了,那就不得不说 WebGL 了!
WebGL
其实很久之前就知道 transform 的 matrix 是这个原理,但是现在要研究 WebGL 才算是顺道把这个原理弄得比较明白。
相对于 CSS 的 transform,转换的是元素的“图像”,WebGL 就真的是纯正的对一大堆点 进行转换。
attribute vec4 a_position;
attribute vec4 a_color;
uniform mat4 u_matrix;
varying vec4 v_color;
void main() {
// Multiply the position by the matrix.
gl_Position = u_matrix * a_position;
// Pass the color to the fragment shader.
v_color = a_color;
}
在一个 vertex shader 里,所有玄机都可以藏在 u_matrix
里,包括模型移动、转向,甚至摄像机参数,都可以通过矩阵处理,例如这里 的透视投影矩阵:
[ f a s p e c t 0 0 0 0 f 0 0 0 0 ( n e a r + f a r ) × r a n g e I n v − 1 0 0 n e a r × f a r × r a n g e I n v × 2 0 ] \begin{bmatrix}
\frac{f}{aspect} & 0 & 0 & 0 \\
0 & f & 0 & 0 \\
0 & 0 & (near + far)\times rangeInv & -1 \\
0 & 0 & near\times far\times rangeInv\times 2 & 0
\end{bmatrix} ⎣ ⎡ a s p ec t f 0 0 0 0 f 0 0 0 0 ( n e a r + f a r ) × r an g e I n v n e a r × f a r × r an g e I n v × 2 0 0 − 1 0 ⎦ ⎤
这里 还有把移动相机的逆矩阵应用到物体造成移动相机效果的操作,不过当然我现在还有一部分不能理解,以后看懂了再展开说吧 🤣
最近看计算机图形相关的资料又找到了学生时代觉得数学十分神奇的感觉,希望以后能慢慢看懂更多算法,先挖个坑,以后一定会写的(鸽子)!