
在构建光线追踪器时,使用视图-投影-viewport矩阵的逆矩阵将屏幕坐标映射回世界空间,会发现其逆矩阵第4行前两个元素(e41、e42)始终为0——这是由透视投影的数学结构决定的固有性质,而非偶然,且在任意相机变换下均严格成立。
在构建光线追踪器时,使用视图-投影-viewport矩阵的逆矩阵将屏幕坐标映射回世界空间,会发现其逆矩阵第4行前两个元素(e41、e42)始终为0——这是由透视投影的数学结构决定的固有性质,而非偶然,且在任意相机变换下均严格成立。
透视投影的本质是将三维世界点通过齐次坐标变换,映射到规范化的设备坐标(NDC)空间,并引入深度信息用于裁剪与Z缓冲。你实现的 perspectiveFovLH 方法生成的是一个典型的左手系、FOV驱动的OpenGL风格透视矩阵,其结构具有明确的稀疏性:
// 简化示意:标准左手系透视矩阵 P(列主序)
P = [
xscale, 0, 0, 0
0, yscale, 0, 0
0, 0, Q, -Q * znear
0, 0, 1, 0
]
注意:该矩阵第4行仅在第3列(e43)和第4列(e44)非零,即 e41 = e42 = 0,且 e43 = 1, e44 = 0。这一结构直接源于透视除法的设计逻辑——w 分量仅依赖于输入点的 z 坐标和固定常量(如 znear),与 x、y 完全无关。
当该透视矩阵与任意仿射变换矩阵 M(如视图变换、平移、旋转、缩放)复合得到总变换 T = P * M 时,由于 M 是仿射的,其最后一行恒为 [0, 0, 0, 1],因此 T 的第4行仅由 M 的第3行线性组合而来(因 P 第4行仅作用于 M 的第3行和第4行)。关键在于:P 的第4行不含对 M 第1、2行(即输入 x、y)的贡献,故 T 的第4行中,x 和 y 的系数(即 T₄₁ 和 T₄₂)必然为 0。
而逆矩阵 R = T⁻¹ 的第4行,本质上解的是“给定齐次输出 s = [sx, sy, sz, sw]ᵀ,如何恢复原始齐次输入 v = [vx, vy, vz, vw]ᵀ”。由上述结构可知,sw 仅由 vz 和 vw 决定;sz 同样仅由 vz 和 vw 线性组合而成。因此,在求解 vw 的表达式中,不涉及 sx 或 sy ——这意味着 R₄₁ = R₄₂ = 0 是恒成立的代数结果,与相机位置、朝向、焦距、近远平面等参数完全无关。
✅ 实践意义与安全假设:
- 在光线投射中,你将屏幕像素 (u,v) 映射为 NDC 点(如 [-1,1]² × [0,1]),再扩展为齐次坐标 s = [x,y,0 or 1,1]。此时计算 w = sₓ·R₄₁ + s_y·R₄₂ + s_z·R₄₃ + R₄₄ 可简化为 w = s_z·R₄₃ + R₄₄,跳过两次乘加运算。
- 该优化绝对安全:只要你的投影矩阵遵循标准透视构造(即 e41=e42=0),且视图变换为仿射(所有实际相机变换均满足),则 R₄₁=R₄₂=0 永真。
⚠️ 注意事项:
- 此性质不适用于正交投影(其 e43=0, e44=1,导致 R₄₁,R₄₂ 通常非零);
- 若引入非仿射变换(如畸变校正、自定义扭曲矩阵),则 M 不再满足 [0,0,0,1] 结构,该结论失效;
- 确保 perspectiveFovLH 中未意外覆盖 e41/e42(你的代码未赋值,默认为0,正确);
- 验证逆矩阵时,建议用数值方法(如 LU 分解)双重校验,避免手工推导误差。
总之,e41=e42=0 不是巧合,而是透视几何与齐次坐标的必然体现——它既是理论可证的不变量,也是光线追踪中值得信赖的性能优化支点。











