利用SH函数绘制空间中动态物体

公式

面片上的光照情况：

$$L=Le+TL$$
$L$：最终要求的光照值；
$Le$ ：直接从光源照射过来的光；
$T$：光照传播因子；
$TL$：从其他面片传播过来的光照

球谐光照常用光照模型 1：

$$L(x,\vec{w_o}）=L_e(x,\vec{w_o})+\int_S f_r(x, \vec{w_i} \rightarrow \vec{w_o})L(x',\vec{w_i})G(x,x')V(x,x')dw_i[1]$$
$L(x,\vec{w_o})$：点 $x$处在$w_o$方向上的光密度
$L_e(x,\vec{w_o})$：物体自身在点 $x$处发出的光线
$f_r(x,\vec{w_i}\rightarrow \vec{w_o})$：表面上 $x$处的BRDF，他将从$w_i$ 方向上进入到该点的光反射并变换到 $w_o$ 方向上
$L(x',\vec{w_i})$：从其他物体上的 $x'$ 处沿 $w_i$ 方向到达的光线
$G(x,x')$：两 $x,x'$ （分别属于不同的物体）间的几何关系
$V(x,x')$：两点 $x,x'$ 间的可见性关系，取值0或1

蒙特卡洛积分

2D点 -> 3D球面点转换：

$$(2arcos(\sqrt{1-\xi_x}),2\pi \xi_y) \rightarrow(\theta,\varphi)$$
近似积分：
$$\int f(x)dx \approx \frac{1}{N}\sum_{i=1}^N f(x_i)w(x_i)=\frac{4\pi}{N}\sum_{i=1}^{N} f(x_i)$$
$N$：采样点数量；
$w(x_i)$权重函数，这里因为均匀随机分布，所以退化为一个常数，大小取决于单位球表面积。

勒让德多项式

勒让德多项式主要属性：

• $(l-m)P_l^m = x(2l-1)pP_{l-1}^m-(l+m-1)P_{l-2}^m$
• $P_m^m=(-1)^m(2m-1)!!(1-x^2)^{m/2}$
• $P_{m+1}^m=x(2m+1)P_m^m$

实数球谐函数

通过勒让德多项式得到的球谐函数实部表示如下：

$$\begin{aligned} y_l^m(\theta,\varphi) = \begin{cases} \sqrt{2}K_l^mcos(m\varphi)P_l^m(cos\theta), &&m>0\\ \sqrt{2}K_l^msin(-m\varphi)P_l^{-m}(cos\theta),&&m<0\\ K_l^0P_l^0(cos\theta),&&m=0 \end{cases} \end{aligned}$$
其中P为相应的伴随勒让德多项式，K为SH球谐函数的随访因子，如下：
$$K_l^m=\sqrt{\frac{(2l+1)(l-|m|)!}{4\pi(l+|m|)!}}$$
由上述参数化表示之后即可以求得蒙特卡洛采样后的球面空间上的球谐函数，其可视化形状一般表示为下图：

球谐投影与重建

对于一个分布于球面空间的函数$f(x)$, 其用球谐函数表示后为$\tilde f(x)$，则 2：

$$\tilde f(x)=\sum_{t=0}^n\sum_{m=-l}^l c_l^my_l^m(s), \qquad c_l^m=\int_Sf(s)y_l^m(s)\mathtt ds$$
对其系数应用蒙特卡洛离散化后得到 3：
$$c_l^m = \int_Sf(s)y_l^m(s)\mathtt ds \approx \frac{1}{N}\sum_{j=1}{N}f(s_j)w(s_j)y_l^m(s_j)\\ c_l^m=\frac{4\pi}{N}\sum_{j=1}^{N}f(s_j)y_l^m(s_j)$$
由此便可以对有一个球面空间上的函数进行球谐投影并还原，投影与重建情形如下：

特性：两个函数的乘积在球面空间上的积分值与它们的球谐系数向量组的点积相同。
从而有[^4]：

$$\int_Sf_1(s)f_2(s)\mathtt ds \approx \int_S\tilde f_1(s) \tilde f_2(s)\mathtt ds = \vec{c_i} \cdot \vec{c_2} = \sum_{i=0}^{N}c_{1_i}c_{2_i}$$

四面体剖分线性插值

假设原本点的位置为$\vec{r}$,对应的重心坐标为$\vec{\lambda}$，则有

$$\left[\begin{matrix} p_{0_x}-p_{3_x} & p_{1_x}-p_{3_x} & p_{2_x} - p_{3_x} \\ p_{0_y}-p_{3_y} & p_{1_y}-p_{3_y} & p_{2_y} - p_{3_y} \\ p_{0_z}-p_{3_z} & p_{1_z}-p_{3_z} & p_{2_z} - p_{3_z} \\ \end{matrix}\right] \left[\begin{matrix} \lambda_0 \\ \lambda_1 \\ \lambda_2 \\ \end{matrix}\right] = \left[\begin{matrix} r_x - p_{3_x} \\ r_y - p_{3_y} \\ r_z - p_{3_z} \\ \end{matrix}\right]$$

规则

固定采样点的不足：

1. 如果采样点比较稀疏则对场景的采样就粗略，影响绘制效果，设置得密集则每帧的采样工作计算量大，影响整体计算速度；
2. 对于动态场景较容易出现采样点失效的情况；
3. 每帧都必须对所有采样点建立立方体并计算该立方体的球谐系数，在此基础上插值计算所有的面片的光照．

改进方法：

1. 采样点随物体的运动而更新，大大降低了采样点失效的概率；
2. 每帧选用合适的采样点进行计算，既保证采样效率又减少了采样工作的计算量；
3. 利用相邻帧的变化较小、每帧可以只更新场景中的一部分物体的光照，将原本大量的计算分散到各帧，提高了效率．

算法

分帧计算：

初始化 ：首先对场景进行初始化计算，包括细分面片、建立面片列表和相互对应关系、建立立方体形状因子查找表、设置物体包围盒的8个顶点为人射参考点、设置出射参考点等
帧计算：

• Step1. 计算到场景边界物体的直接光照．直接光照是具有自我发光能力的物体向外发射出的直接照射到其他物体的光照，为了避免与后续计算重复，此处不计算光源照射到场景内部物体的直接光照，只计算光源到场景边界上的直接光照，并将其记录在每个边界面片的数据结构中．
• Step2. 选择一个场景内部的物体，计算它的光照．与以往的方法不同，每一帧都不是计算所有物体的光照，而是有所选择地同时对接收到的直接光照和间接光照进行计算，采用球谐方法并插值．
• Step3. 计算所选择的物体对边界的间接光照．场景内部的物体从周围环境吸收能量，再将部分能量反射到环境中去，这部分就是间接光照．将物体视为一个整体发射间接光能量，这个过程和直接光照非常类似．由于计算时每帧所选择的物体一般不同，所以维护每个边界面片接收到的间接光照能量时只更新当前物体对它的影响．
• Step4. 绘制场景面片．场景每个面片的光能量包括物体自身拥有的光能量和从其他物体处吸收来的光能量2个部分．将能量值转化为颜色值并绘制．

实现

球谐光照

通过以下变换转换$\theta, \varphi$到笛卡尔坐标系：

$$cos(\theta) = \frac{x}{r}$$
$$sin(\varphi) = \frac{y}{rsin(\theta)}=1$$
$$cos(\varphi) = \frac{x}{rsin(\theta)}=\frac{x}{y}$$
使用笛卡尔坐标系表示结果如下：

• $y_0^0(\theta, \varphi) = \frac{1}{2}\cdot\sqrt{\frac{1}{\pi}}$
• $y_1^{-1}(\theta,\varphi) = -\frac{1}{2}\cdot\sqrt{\frac{3}{\pi}}\cdot\frac{y}{r}$
• $y_1^0(\theta,\varphi) = \frac{1}{2}\cdot\sqrt{\frac{3}{\pi}}\cdot\frac{z}{r}$
• $y_1^1(\theta,\varphi)=-\frac{1}{2}\cdot\sqrt{\frac{3}{\pi}}\cdot\frac{x}{r}$
• $y_2^{-2}(\theta,\varphi)=-\frac{1}{2}\cdot\sqrt{\frac{15}{\pi}}\cdot\frac{yx}{r^2}$
• $y_2^{-1}(\theta,\varphi)=-\frac{1}{2}\cdot\sqrt{\frac{15}{\pi}}\cdot\frac{yz}{r^2}$
• $y_2^0(\theta,\varphi)=\frac{1}{4}\cdot\sqrt{\frac{5}{\pi}}\cdot\frac{2z^2-x^2-y^2}{r}$
• $y_2^1(\theta,\varphi)=-\frac{1}{2}\cdot\sqrt{\frac{15}{\pi}}\cdot\frac{zx}{r^2}$
• $y_2^2(\theta,\varphi)=\frac{1}{4}\cdot\sqrt{\frac{15}{\pi}}\cdot\frac{x^2-y^2}{r^2}$

对于光照方程[1],考虑一个diffuse表面时，各个方向行的反射光线均相同，用一个衡量代替。从而有：

$$L(x,\vec{w_0}) = k\int_SL(x',\vec{w_i})G(x,x')V(x,x')\mathtt dw_i$$
根据特性 [4]可以将 $L(x', \vec{w_i})$ 与其他部分分开处理，从而进一步简化模型 [5]：
$$L(x',\vec{w_0})=\int_S G(x,x')V(x,x')\mathtt dw_i$$
其中 $G(x,x')$ 也是简单的光线与定点法向量之间的夹角关系，从而计算夹角cos值进行储存即可。

脚注

---

1. 常用光照模型 ↩
2. 球面空间函数的球谐表示 ↩
3. 蒙特卡洛方法离散后的系数表示方法 ↩