物理渲染学习笔记（二）——光的传播

一、假设
首先我们假设介质又许多非常小并且随机分布的微粒组成，微粒之间的距离相比与微粒本身大小可视为无限大，因此每个微粒的光照行为是独立事件，但我们不描述单个微粒的光照过程，我们只描述光穿过介质时整体的概率分布。

---

二、光子的交互

当一个光子穿越一堆微粒时，可能一个也碰不到直接通过，也有可能撞到微粒，撞到的概率与介质的系数 $\sigma_t$ 相关，撞到的时候，有两种可能，吸收或者散射，与其相关的介质系数分别为 $\sigma_a$ 和 $\sigma_s$, 两种都会导致该方向的 radiance 变化，最后，我们可以得到 $\sigma_t = \sigma_a + \sigma_s$ 。

---

2.1 Absorption
考虑一束极小的光束沿 $\vec \omega$ 穿过介质，则 $x$ 点射进去的光子数量与入射radiance $L(x\to \vec \omega)$ 成比例，光束上的每一部分 $\Delta t$， 一部分光子与介质交互并被吸收，若该步长内的吸收系数为 $\sigma_a(x+t\vec \omega)$， 则 $\sigma_a(x+t\vec \omega)\Delta t$部分的光子被吸收，因此，剩余的光子可被表示为：

$$L((x+t\vec \omega)\to \vec \omega) = L(x\to \vec \omega)(1-\sigma_a(x+t\vec \omega)\Delta t)$$ 整理可得： $$\frac{L((x+t\vec \omega)\to \vec \omega) -L(x\to \vec \omega)}{\Delta t} = -\sigma_a(x)L(x\to\vec \omega)$$
当步长无限小时，可得：
$$\lim_{\Delta t \to 0}\frac{L((x+t\vec \omega)\to \vec \omega) -L(x\to \vec \omega)}{\Delta t} = (\vec \omega\cdot \nabla_a)L(x\to \vec \omega)=-\sigma_a(x)L(x\to \vec \omega)$$
其中 $\nabla_a$ 表示光子被吸收的梯度， $(\vec \omega\cdot \nabla_a)$ 则表示方向导数

---

2.2 Out-scattering
同理我们可以得到out-scattering造成的radiance变化：

$$(\vec \omega\cdot \nabla_o)L(x\to \vec \omega)=-\sigma_s(x)L(x\to \vec \omega)$$

同时考虑吸收和Out-scattering，可以得到：

$$(\vec \omega\cdot \nabla_t)L(x\to \vec \omega)=-(\sigma_a+\sigma_t)L(x\to \vec \omega)=-\sigma_t(x)L(x\to \vec \omega)$$

上面的公式，其实可以表述为透射率 $T_r$， 透射率表示为能不与微粒进行碰撞的光子比例，对上式进行积分可得：

$$T_r(x' \leftrightarrow x) = e^{-\tau(x' \leftrightarrow x)}$$
参数 $\tau$ 被称为光学厚度，定义如下： $$\tau(x' \leftrightarrow x) = \int_0^d \sigma_t(x + t\vec \omega)dt$$
其中 $t \in [0,d]$ ,在均匀介质中， $\tau$ 可被计算为： $$\tau(x' \leftrightarrow x) = d\sigma_t$$

回到透射率，透射率通常在0到1之间，当只考虑吸收和out-scattering时，透射率可用来计算减少的radiance：

$$L(x\to \vec \omega) = T_r(x' \leftrightarrow x)L(x'\to -\vec \omega)$$
此外，一个点到自身的透射率总为1，直线上的三个点透射率满足下的等式： $$T_r(x \leftrightarrow x'') = T_r(x \leftrightarrow x')T_r(x' \leftrightarrow x'')$$

---

2.3 Emission和In-scattering

In-scattering：

$$(\vec \omega \cdot \nabla_i)L(x\to \vec \omega) = \sigma_s(x)L_i(x\to \vec \omega)$$
Emiision： $$(\vec \omega \cdot \nabla_e)L(x\to \vec \omega) = \sigma_a(x)L_e(x\to \vec \omega)$$

2.4 相函数

就像BRDF，相函数满足可交换性：

$$p(x,\vec \omega' \leftrightarrow \vec \omega)$$
归一性： $$\int_{\Omega_{4\pi}}p(x,\vec \omega' \leftrightarrow \vec \omega)d\vec \omega' = 1, \forall \vec \omega$$

各向同性相函数：

漫反射就是完美的各向同性散射，其相函数为常数:

$$p_I(x,\vec \omega' \leftrightarrow \vec \omega) = \frac{1}{4\pi}$$

Henyey-Greenstein相函数：
使用最为广泛的相函数，只取决于 $\theta$， 定义如下：

$$p_{HG}(x,\theta) = \frac{1-g^2}{4\pi(1+g^2-2gcos\theta)^{1.5}}$$ 其中 $g\in [-1,1]$， 决定了前向和后向散射的强度，在物理上， $g$ 表示散射方向的平均余弦值，大于零时是前向散射，小于零时是后向散射，等于零时则表示各向同性。同时，更加复杂的散射函数可用Henyey-Greenstein用如下计算做近似： $$g(x) = \int_{\Omega_{4\pi}}p(x,\vec \omega' \leftrightarrow \vec \omega)cos\theta d\vec \omega'$$

---

三、体渲染方程
综合上述的absorption、out-scattering、in-scattering、emission，我们可得得到光在介质中传播的完整函数：

$$(\vec \omega \cdot \nabla)L(x\to \vec \omega) = \underbrace{\underbrace{\sigma_a(x)L(x\to \vec \omega)}_{absorption} - \underbrace{\sigma_s(x)L(x\to \vec \omega)}_{out-scattering}}_{extinction} + \underbrace{\sigma_a(x)L_e(x\to \vec \omega)}_{emission} + \underbrace{\sigma_s(x)L_i(x\to \vec \omega)}_{in-scattering}$$

如图所示，对上式两边在全方向进行积分，并代入透射率 $T_r$ 可以得到：

$$L(x\to \vec \omega) = \underbrace{T_r(x\leftrightarrow x_s)L(x_s\to -\vec \omega)}_{reduced\space surface\space radiance} + \underbrace{\int_0^sT_r(x\leftrightarrow x_t)\sigma_a(x)L_e(x_s\to -\vec \omega)dt}_{accumulated \space emitted \space radiance} + \underbrace{\int_0^sT_r(x\leftrightarrow x_t)\sigma_s(x_t)L_i(x_s\to -\vec \omega)dt}_{accumulated \space in-scattering\space radiance}$$