【Unity Shader入门精要 第4章】数学基础(二)

2024-05-04 17:04

本文主要是介绍【Unity Shader入门精要 第4章】数学基础(二),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

1. Unity中的坐标空间

1.1 五个坐标空间

模型空间

  • 模型自身的3D坐标系空间,左手坐标系
  • 是一个相对空间,坐标轴指向随模型旋转变化
  • 当物体有父节点时,Transform组件中各属性的值表示的即为该物体在其父物体的模型空间中的值
  • 当模型顶点传入顶点着色器时,其中的空间信息均为自身模型空间的信息

世界空间

  • 世界3D坐标系空间,左手坐标系
  • 是一个绝对空间,坐标轴指向不会发生变化
  • 当物体没有父节点时,Transform组件中各属性的值表示的即为该物体在世界空间中的值

观察空间

  • 摄像机的模型空间,也是一个3D坐标系空间
  • 是Unity中唯一的一个右手坐标系空间
  • 以摄像机位置为坐标原点
  • 摄像机面向的方向为Z轴负方向

齐次裁剪空间

  • 左手坐标系
  • 其本质是视锥体包围的空间,或摄像机渲染范围的空间
  • 齐次裁剪空间也被称作投影空间,但实际上这时候并没有进行真正的投影,因为还没有对矢量进行降维,该空间仍然是一个3D空间
  • 齐次空间范围内的部分可见,范围外的部分为不可见,需要从渲染中剔除
  • 在透视投影下:
    • 由近裁截面距离、远裁截面距离、FOV以及横纵比决定空间范围
    • 形状为平头四棱锥
    • FOV决定锥体竖直方向的张开角度
    • 根据近裁截面、远裁截面和FOV可以求出近裁截面和远裁截面的高度,根据横纵比可以求出近裁截面和远裁截面的宽度
  • 在正交投影下:
    • 由近裁截面距离、远裁截面距离、SIZE以及横纵比决定空间范围
    • 形状为立方体
    • SIZE决定立方体竖直方向上高度的一半
    • 根据近裁截面、远裁截面和SIZE可以求出近裁截面和远裁截面的高度,根据横纵比可以求出近裁截面和远裁截面的宽度
  • 齐次空间仍然为一个非归一化的三维空间

屏幕空间

  • 左手坐标系
  • 对应当前显示屏幕的2D像素空间
  • 空间范围为(0, 0 )到(ScreenMaxPixelX, ScreenMaxPixelY)

1.2 四个转换过程

模型变换

  • Model Transform —— M
  • 将顶点从模型空间转换到世界空间

观察变换

  • View Transform —— V
  • 将顶点从世界空间转换到观察空间

投影变换

  • Projection Tranform —— P
  • 将顶点从观察空间转换到投影空间(齐次裁剪空间)
  • 根据摄像机的设置不同,分为透视投影(Perspective Projection)和正交投影(Orthographic Projection)两种投影变换方式
  • 同样,这里所谓的投影变换并不是真正进行降维投影,仍然是3D空间到3D空间的变换
  • 所谓投影变换只是对矢量乘以一个变换矩阵(称为投影矩阵或裁剪矩阵),得到矢量在当前齐次空间中真正的齐次坐标
  • 变换矩阵的推导过程并不重要,变换矩阵长什么样也不重要,重要的是我们知道矢量经过与变换矩阵的计算后,其w分量发生了变化,不再是初始的1和0,而是变成了-Z(点)和1(向量)

屏幕映射

  • 以上MVP三种变换都是在顶点着色器中完成的,回忆第二章中渲染流水线的几何阶段,在顶点着色器和屏幕映射阶段之间还有一个必不可少的裁剪阶段,因此在做真正的屏幕映射之前,还需要先进行一步裁剪
  • 上面我们提到过,齐次裁剪空间是一个非归一化的空间,在不同的摄像机设置下,视锥体的范围是不同的,为了能够对不同的视锥体范围进行统一的裁剪处理,需要先将齐次坐标映射到一个统一的范围之内,这一步通过将齐次坐标的所有分量都除以w分量来完成,称为齐次除法(感觉齐次除法就是一个从四维投影到三维的过程)
  • 经过齐次除法后,我们得到了一个(x/w, y/w, z/w, 1)的坐标,称为归一化设备坐标(NDC),而原来的裁剪空间也被投影到了一个(-1,1)的正立方体空间范围(DX为0到1)
  • 此时就可以对所有的NDC坐标进行统一的裁剪处理,只需要考虑坐标是否超出(-1,1)的范围即可
  • 裁剪结束后,所有剩余的顶点都是在屏幕范围内的点,需要根据其X和Y的值映射到(0, 0)至(MaxScreenPixelX, MaxScreenPixelY)的范围内,可以看到,这时候才真正从3D空间投影到了2D空间
  • 剩余的Z值则表示该顶点到屏幕的距离,即深度值

1.3 空间转换图示

在这里插入图片描述

2. 坐标空间变换矩阵

对于坐标空间A和坐标空间B,当我们知道了:

  • B空间的原点 O 在A空间的表示为 OriginalB
  • B空间的 X 轴在A空间的表示为 AxisBx
  • B空间的 Y 轴在A空间的表示为 AxisBy
  • B空间的 Z 轴在A空间的表示为 AxisBz

则从坐标空间B到坐标空间A的转换矩阵MB→A(4*4)为:
( ∣ ∣ ∣ ∣ A x i s B x A x i s B y A x i s B z O r i g i n a l B ∣ ∣ ∣ ∣ 0 0 0 1 ) \left( \begin{matrix} | & | & | & |\\ AxisB~x~ & AxisB~y~ & AxisB~z~ & OriginalB\\ | & | & | & |\\ 0 &0&0&1 \end{matrix} \right) AxisB x 0AxisB y 0AxisB z 0OriginalB1
即将X轴、Y轴、Z轴、原点按列组合形成的4 * 4矩阵(推导过程懒得写了)

由于向量不存在空间位置属性,因此对于向量的变换矩阵可以不关心原点,故可以直接截取3 * 3的轴矩阵 MB→A(3*3):
( ∣ ∣ ∣ A x i s B x A x i s B y A x i s B z ∣ ∣ ∣ ) \left( \begin{matrix} | & | & | \\ AxisB~x~ & AxisB~y~ & AxisB~z~ \\ | & | & | \end{matrix} \right) AxisB x AxisB y AxisB z 

根据以上矩阵可以做如下推导:

  • 对于任意B空间内的向量VB,其在A空间的表示 VA = MB→A(3 * 3) * VB
  • 两边同时乘以 MB→A(3 * 3) 的逆矩阵 MB→A(3 * 3)-1
  • 得到等式 MB→A(3 * 3)-1 * VA = VB

即:MA→B(3 * 3) = MB→A(3 * 3)-1

同时,由于 MB→A(3 * 3)是由坐标空间B的三个坐标轴组合得到的矩阵,而坐标轴组合的矩阵一定是正交矩阵,因此 MB→A(3 * 3)的转置矩阵与逆矩阵相同

因此,此时从坐标空间A反向转换到坐标空间B的转换矩阵 MA→B(3 * 3) = MB→A(3 * 3)-1 = MB→A(3 * 3)T,故 MA→B(3 * 3) 为:
( — A x i s B x — — A x i s B y — — A x i s B z — ) \left( \begin{matrix} — & AxisB~x~ & — \\ — & AxisB~y~ & — \\ — & AxisB~z~ & — \end{matrix} \right) AxisB x AxisB y AxisB z 

3.法线变换

普通的转换矩阵适用于在不同的坐标空间之间对点和普通向量进行空间变换,但是如果直接用该矩阵对法线进行变换,就会出错

这是因为法线并不是空间内真实存在的向量,法线只是对空间内某一平面垂直的一个相对向量,或者说是对空间内真实存在的向量虚构出来的相对向量,具有一定的相对性

当进行空间变换后,这种相对性无法保证依然生效,因此可能导致变换后的法线与变换后的平面并不垂直

而切线是由空间内相邻两点确定的,也就是空间中真实存在的向量,因此可以通过普通的转换矩阵进行空间变换,并且我们可以利用切线来推导出使用于法线的转换矩阵

假设T为切线,N为法线,M为切线的转换矩阵,G为法线的转换矩阵
TB = M * TA
NB = G * NA

由切线与法线垂直可得:
TB · NB = (M * TA) · (G * NA) = 0
注意,上式表示的是两个垂直向量的点乘结果为0

如果将上式改成矩阵乘法来表示,就需要将TB 转置成列矩阵后再与行矩阵 NB 相乘
于是得到:
(M * TAT · (G * NA) = TAT * MT * G * NA = TAT * (MT * G) * NA = 0

由于TAT * NA = 0,因此要使上式成立,只要 MT * G = I,即G = (MT)-1 = (M-1)T

于是问题转化为求普通转换矩阵 M 的逆矩阵 M-1

  • 如果只存在旋转变换,由于旋转矩阵为正交矩阵,M-1 = MT,因此 G = (M-1)T = (MT)T = M
  • 如果只存在旋转变换和统一缩放变换,假设统一缩放系数为K,则G = 1/k * M(这里一直没明白,为什么不是1/k2 * M,有没有大佬给解释一下??)
  • 如果存在非统一缩放变换或平移变换,就只能求逆矩阵了

4. Unity内置变换矩阵

Unity内置了一些变量,可以轻松的获取对应的转换矩阵,这里提供的都是4*4的矩阵

变量作用
_Object2World模型空间到世界空间
_World2Object世界空间到模型空间
变量作用
UNITY_MATRIX_MVP模型空间到齐次裁剪空间
UNITY_MATRIX_V世界空间到观察空间
UNITY_MATRIX_P观察空间到齐次裁剪空间
UNITY_MATRIX_VP世界空间到齐次裁剪空间
变量作用
UNITY_MATRIX_MV模型空间到观察空间
UNITY_MATRIX_T_MV模型空间到观察空间的转置
UNITY_MATRIX_IT_MV模型空间到观察空间的逆转置

5. Unity内置摄像机参数

Unity内置了部分变量可以用来获取对应摄像机参数

变量类型作用
_WorldSpaceCameraPosfloat3摄像机世界坐标
变量类型作用
_ScreenParamsfloat4(width, height, 1 + 1/w, 1 + 1/h)
_ZBufferParamsfloat4(1 - F/N, F/N, x/F, y/F)
_ProjectinParamsfloat4(1/-1, N, F, 1 + 1/F),首位为-1时表示使用的是翻转的投影矩阵
Unity_OrthoParamsfloat4(w, h, _, 0/1) 最后位为0时为透视相机,为1时为正交相机
变量类型作用
Unity_CameraProjectionfloat4x4摄像机投影矩阵
Unity_CameraInvProjectionfloat4x4摄像机投影矩阵的逆矩阵
Unity_CameraWorldClipPlanes[6]float4左右上下远近

这篇关于【Unity Shader入门精要 第4章】数学基础(二)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/959681

相关文章

Spring Security 从入门到进阶系列教程

Spring Security 入门系列 《保护 Web 应用的安全》 《Spring-Security-入门(一):登录与退出》 《Spring-Security-入门(二):基于数据库验证》 《Spring-Security-入门(三):密码加密》 《Spring-Security-入门(四):自定义-Filter》 《Spring-Security-入门(五):在 Sprin

零基础学习Redis(10) -- zset类型命令使用

zset是有序集合,内部除了存储元素外,还会存储一个score,存储在zset中的元素会按照score的大小升序排列,不同元素的score可以重复,score相同的元素会按照元素的字典序排列。 1. zset常用命令 1.1 zadd  zadd key [NX | XX] [GT | LT]   [CH] [INCR] score member [score member ...]

uva 10014 Simple calculations(数学推导)

直接按照题意来推导最后的结果就行了。 开始的时候只做到了第一个推导,第二次没有继续下去。 代码: #include<stdio.h>int main(){int T, n, i;double a, aa, sum, temp, ans;scanf("%d", &T);while(T--){scanf("%d", &n);scanf("%lf", &first);scanf

uva 10025 The ? 1 ? 2 ? ... ? n = k problem(数学)

题意是    ?  1  ?  2  ?  ...  ?  n = k 式子中给k,? 处可以填 + 也可以填 - ,问最小满足条件的n。 e.g k = 12  - 1 + 2 + 3 + 4 + 5 + 6 - 7 = 12 with n = 7。 先给证明,令 S(n) = 1 + 2 + 3 + 4 + 5 + .... + n 暴搜n,搜出当 S(n) >=

uva 11044 Searching for Nessy(小学数学)

题意是给出一个n*m的格子,求出里面有多少个不重合的九宫格。 (rows / 3) * (columns / 3) K.o 代码: #include <stdio.h>int main(){int ncase;scanf("%d", &ncase);while (ncase--){int rows, columns;scanf("%d%d", &rows, &col

数论入门整理(updating)

一、gcd lcm 基础中的基础,一般用来处理计算第一步什么的,分数化简之类。 LL gcd(LL a, LL b) { return b ? gcd(b, a % b) : a; } <pre name="code" class="cpp">LL lcm(LL a, LL b){LL c = gcd(a, b);return a / c * b;} 例题:

【生成模型系列(初级)】嵌入(Embedding)方程——自然语言处理的数学灵魂【通俗理解】

【通俗理解】嵌入(Embedding)方程——自然语言处理的数学灵魂 关键词提炼 #嵌入方程 #自然语言处理 #词向量 #机器学习 #神经网络 #向量空间模型 #Siri #Google翻译 #AlexNet 第一节:嵌入方程的类比与核心概念【尽可能通俗】 嵌入方程可以被看作是自然语言处理中的“翻译机”,它将文本中的单词或短语转换成计算机能够理解的数学形式,即向量。 正如翻译机将一种语言

Java 创建图形用户界面(GUI)入门指南(Swing库 JFrame 类)概述

概述 基本概念 Java Swing 的架构 Java Swing 是一个为 Java 设计的 GUI 工具包,是 JAVA 基础类的一部分,基于 Java AWT 构建,提供了一系列轻量级、可定制的图形用户界面(GUI)组件。 与 AWT 相比,Swing 提供了许多比 AWT 更好的屏幕显示元素,更加灵活和可定制,具有更好的跨平台性能。 组件和容器 Java Swing 提供了许多

【IPV6从入门到起飞】5-1 IPV6+Home Assistant(搭建基本环境)

【IPV6从入门到起飞】5-1 IPV6+Home Assistant #搭建基本环境 1 背景2 docker下载 hass3 创建容器4 浏览器访问 hass5 手机APP远程访问hass6 更多玩法 1 背景 既然电脑可以IPV6入站,手机流量可以访问IPV6网络的服务,为什么不在电脑搭建Home Assistant(hass),来控制你的设备呢?@智能家居 @万物互联

poj 2104 and hdu 2665 划分树模板入门题

题意: 给一个数组n(1e5)个数,给一个范围(fr, to, k),求这个范围中第k大的数。 解析: 划分树入门。 bing神的模板。 坑爹的地方是把-l 看成了-1........ 一直re。 代码: poj 2104: #include <iostream>#include <cstdio>#include <cstdlib>#include <al