2. threejs源码阅读——math/matrix

Matrix接口

Vector接口一样,threejs也在Matrix3.d.ts文件中提供了一个Matrix接口,并且在Matrix3(3x3矩阵类)和Maxtrix4(4x4矩阵类)中实现了该接口。

Matrix接口

我们可以看到,矩阵中的元素是用数组的方式储存的。由于不同维度的矩阵内所用的元素个数不同,所以Matrix3和Matrix4的set方法是分开声明在每个类中的。这里我们以Matrix3.set为例:

Matrix3.set

可以看到,在set时是以行优先的顺序进行传值的,也就是第一行第一个,第一行第二个,第一行第三个,第二行第一个,第二行第二个...

M_{3\times3} = \begin{pmatrix}  n11 & n12 & n13 \\  n21 & n22 & n23 \\  n31 & n32 & n33 \end{pmatrix}

但是在set方法的实现中,可以看到elements数组内是以列优先的顺序进行存储的:

Matrix3.set方法实现

M_{3\times3} = \begin{pmatrix}  n11 \rightarrow e0 & n12 \rightarrow e3 & n13 \rightarrow e6  \\  n21 \rightarrow e1 & n22\rightarrow e4 & n23\rightarrow e7 \\  n31\rightarrow e2  & n32\rightarrow e5 & n33\rightarrow e8 \end{pmatrix}

简单点说,这里只是存储方式的不同。由于大多数人都习惯以行优先的方式考虑矩阵,所以threejs中所有的文档都是以行优先的方式表示的。但如果我们想要阅读源码,就要注意到[n11, n12, n13, n21 ... n33]并不是按顺序储存在elements中。事实上如果我们按顺序读取elements中的元素并以行优先的方式组成矩阵,它将是原本矩阵的转置矩阵。

https://en.wikipedia.org/wiki/File:Row_and_column_major_order.svg

在阅读后续的矩阵计算相关的方法中,我们一定要时刻注意threejs的矩阵元素是以列优先的方式存储的!

如同Vector的类型,Matrix类型中所有返回值为自身类对象的方法也是原地修改该对象。例如Transpose转置

Matrix3.prototype.transpose

所有的Matrix也提供了clone和copy方法,用于深拷贝不同Matrix对象。

那么threejs中的向量与矩阵进行计算是以什么形式进行的呢?这里我们举一个例子。


行向量还是列向量?

在Vector3类中提供了appyMatrix3这个方法,这个方法以一个3x3的矩阵作为参数,相乘后原地修改向量的值。

Vector3.applyMatrix3

如果我们将Vector3看做三行一列的矩阵(3x1)的话,那么就需要左乘3x3矩阵,最后得到列向量:

V =\begin{pmatrix}  e0 & e3 & e6 \\  e1 & e4 & e7 \\  e2 & e5 & e8\end{pmatrix}\begin{pmatrix}x\\y\\z\end{pmatrix}=\begin{pmatrix}e0*x + e3 *y+e6*z \\e1*x+e4*y+e7*z\\e2*x+e5*y+e8*z\end{pmatrix}


如果我们将Vecot3看做一行三列的矩阵(1x3)的话,那么就需要右乘3x3矩阵,最后得到行向量:

V =\begin{pmatrix}x&y&z\end{pmatrix}\begin{pmatrix}  e0 & e3 & e6 \\  e1 & e4 & e7 \\  e2 & e5 & e8\end{pmatrix}=\begin{pmatrix}x*e0+y*e1+z*e2 &x*e3+y*e4+z*e5 &x*e6+y*e7+z*e8\end{pmatrix}

在这里我们发现:threejs中的向量和矩阵的乘法是以列向量左乘矩阵的方式实现的。

接下来,我们再假设一个向量

V_{xyz} = (1,2,3)

加入齐次坐标w=1,来表示三维空间中的一个点(1,2,3)

P = (1,2,3,1)

然后我们构建一个仿射变换矩阵做缩放+平移运算

T = \begin{pmatrix}2 & 0 & 0 & 5\\0 & 2 & 0 & 5\\0 & 0 & 2 & 5\\0 & 0 & 0 & 1\end{pmatrix}

即沿着xyz轴均缩放2倍,然后沿着xyz移动5个单位距离。最终结果是:

P^ \prime = \begin{pmatrix}2 & 0 & 0 & 5\\0 & 2 & 0 & 5\\0 & 0 & 2 & 5\\0 & 0 & 0 & 1\end{pmatrix}\begin{pmatrix}1\\2\\3\\1\\\end{pmatrix}= \begin{pmatrix}2 *1+ 0*2+0*3 + 5 * 1 = & 7\\0 *1+ 2*2+0*3 + 5 * 1 = &9\\0 *1+ 0*2+2*3 + 5 * 1 = &11\\0 *1+ 0*2+0*3 + 1* 1 = &1\\\end{pmatrix}

这个计算我们会使用Vector3.applyMatrix4

Vector3.applyMatrix4

接下来验证一下:

初始化向量和矩阵并调用applyMatrix4

控制台打印结果:

控制台打印

结论

threejs中的矩阵以行优先的形式初始化,以列优先的形式储存。

threejs中的向量采用列向量,与矩阵做左乘运算

threejs中一个三维向量与4x4的矩阵进行运算时,会先给三维向量补充第四个值w,也就是齐次坐标,w值与矩阵的第四行有关系:

const w = 1 / ( e[ 3 ] * x + e[ 7 ] * y + e[ 11 ] * z + e[ 15 ] );

在上面的例子中,我们使用了一个缩放+平移的仿射变换矩阵作为例子,此时w为1。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 211,743评论 6 492
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,296评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 157,285评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,485评论 1 283
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,581评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 49,821评论 1 290
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,960评论 3 408
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,719评论 0 266
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,186评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,516评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,650评论 1 340
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,329评论 4 330
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,936评论 3 313
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,757评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,991评论 1 266
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,370评论 2 360
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,527评论 2 349