计算机动画5-7

光照模型

光源

任何发出辐射能量的对象——发光体和反射体及其混合
属性：位置、颜色、方向、形状

分类

• 点光源：所有光向四周平均发射 适用于比对象小得多的光源或离场景不太近的大光源
• 无穷远出光源（近似平行光源）：模拟类似太阳光线

辐射强度衰减

$I$ 随 $1/d_l^2$衰减
问题:

• 接近光源是强度变化过大
• $d_l$ 过大时,变化太小

改进:$f_{radatten}(d_1)=\{\begin{array}{l}1.0光源在无穷远处\\ \frac{1}{a_0+a_1{d}_l+a_2{d}_l^2}光源为局部光源\end{array}$

方向和投射效果

光线方向向量$V_{light}$和方向范围${\theta }_l$构成一个圆锥体的可照射区域
记光源位置道一个对象位置的方向向量为$V_{obj}$
则判断其是否可被照射到:$\alpha =arccos(V_{obj}\cdot V_{light})<{\theta }_l$
沿着光源的方向向量光强最强,家督增加,强度衰减
$$\begin{gathered} f_{angatten}(\alpha )=co{s}^{a_l}\alpha , \\ 其中,{\alpha }_l为衰减指数 \end{gathered}$$

基本光照模型(反射)

漫反射

粗糙表面，将光向各个防线发散出去

理想漫反射体

入射光在各个方向以相同强度发散而与观察者位置无关

反射光强度可以根据朗伯余弦定理

$$\begin{gathered} 强度=\frac{单位时间辐射能}{投影面积}\to \frac{cos{\varphi }_N}{dAcos{\varphi }_N}, \\ 其中dA是单位面积,{\varphi }_N是表面法向量与辐射能方向夹角 \end{gathered}$$
假设每一表面都按照理想漫反射体对待,则可以将模型简化为$I_d=I\ast k_d$,$k_d$是漫反射系数

• 环境光漫反射
  记环境光为$I_a$,则环境光对漫反射的贡献可表示为:$I_{a}mbdiff=k_a{I}_a$,但是只有环境光缺乏立体感
• 光源漫反射
  记光源强度为$I_l$,入射光与反射平面的夹角为$\theta$则漫反射光强度为$I_{l,diff}=k_d{I}_{l}cos\theta$
  上述两者综合可得漫反射公式如下:
  $$I_{diff}=\{\begin{array}{l}{k}_a{I}_a+k_{d}Il(N\cdot L)N\cdot L>0\\ k_a{I}_{a}else\end{array}其中L为反射光线单位向量,N为反射表面单位法向量$$

镜面反射

光滑表面上的高光, 是由接近镜面反射角的一个汇聚区域内，入射光的全部或绝大部分成为反射光所导致的

非理想反射体

反射方向分布在向量R周围的有限范围内

• 较光滑表面的镜面反射范围较小
• 粗糙的对象表面有较大的反射范围
  

Phong镜面反射模型

镜面反射光强度与$co{s}^{n_s}\varphi$成正比,即$I_{l,spec}\propto I_{l}co{s}^{n_s}\varphi$
其中$n_s$为镜面反射参数


$$\begin{gathered} I_{l,spec}=\{\begin{array}{l}{k}_s{I}_l(V\cdot R{)}^{n_s}V\cdot R>0和N\cdot L>0\\ 0else\end{array} \\ 其中,V为视点方向,R为反射R=(2N\cdot L)N-L \end{gathered}$$
简化:

• L与V之间的半角向量H来计算镜面反射范围:$H=\frac{L+V}{|L+V|}$
• 用$cos\alpha$代替$cos\varphi$
• 用$N\cdot H$代替$V\cdot R$

漫反射和镜面反射合并:
$$I=I_{ambdiff}+\sum _{l=1}^n(I_{l,diff}+I_{l,spec})=k_a{I}_a+\sum _{l=1}^n{I}_l[k_d(N\cdot L)+k_s(N\cdot H{)}^{n_s}]$$

透射模型

透明对象表面同时产生反射光和折射光

反射光

• 根据Snell公式计算折射角${\theta }_r$:$sin{\theta }_r=\frac{{\eta }_i}{et{a}_r}sin{\theta }_i$
• 折射方向向量T:$T=(\frac{{\eta }_i}{{\eta }_r}cos{\theta }_i-cos{\theta }_r)N-\frac{{\eta }_i}{{\eta }_r}L$

折射光

假设：

• 不考虑折射导致的路径平移
• 各对象间的折射率不变
• 折射角=入射角

简化光强:$I=k_t{I}_{trans}$,其中$k_t$为该点的透射系数

总光强:$I=(1-k_t)I_{refl}+k_t{I}_{trans}$

全局光照

光线跟踪算法

跟踪光线在场景中的反射和折射,并计算对光强的作用
优点:真实
缺点:计算量大
是一种高度视相关的方法,

基本算法

1. 为光线跟踪算法建立一个投影参考点在z轴、投影平面为xy平面的坐标系统
2. 从投影参考点出发，穿过每个像素中心进入场景
3. 计算与该光线相交的所有表面，找到最近的交点（可见面）
4. 计算该光线在此交点上反射和折射光强度
5. 构造光线跟踪树
   
   构造的停止条件:

• 该光线不和任意表面相交
• 该光线与一个光源相交且该光源不是一个反射面
• 该树达到最大允许深度(常用)

光线与对象表面计算求交

• 光线由初始位置$P_0$和单位向量$u$描述$P=P_0+su$,$s$为$P_0$到光线与$P$点的距离
• 由球面方程$|P-P_c{|}^2-r^2=0$联立可得$s=u\cdot \Delta P\pm \sqrt{(u\cdot \Delta P{)}^2-|\Delta P{|}^2+r^2}$

减少对象求交计算量的方法

• 包围盒
• 空间分割法
  • 将整个场景包含在一个立方体中
  • 将立方体逐次分割直至每个子立方体所包含的对象表面或表面数目小于等于一个预定的最大值
  • 跟踪穿过立方体的光线，仅需对包含表面的单元执行求交测试。光线所交的第一个对象表面记为可见面
  • 该过程继续直到找到一个相交的对象表面或光线射出包围场景的立方体

纹理与表面细节

纹理映射

2D映射、3D映射

方法：

• 建立像空间和纹理空间映射关系
• 纹理坐标0到1.0

二维映射

不同的Map Shape产生不同的效果，从纹理表中查找对应的颜色

• 平面
• 圆柱体
• 球体
• 盒子

三维映射

利用对象每个点$(x,y,z)$直接计算纹理颜色,不使用map shape

噪声

凹凸映射

纹理映射无法模拟粗糙表面物体
使用扰动函数并在光照模型中计算使用扰动法向量

• 令$P(u,v)$表示一个参数曲面上的点,$N=P_u\times P_v$,其中$P_u$和$P_v$是关于参数$u$和$v$的偏导数
• 增加一个扰动函数(凹凸函数)$b$,在表面法向量$n=\frac{N}{|N|}$方向上增加凹凸效果:$P^{\prime }(u,v)=P(u,v)+b(u,v)n$
• 表面法向量为$N^{\prime }=P_u^{\prime }\times P_v^{\prime }$
• $P$关于$u$的偏导数:$P_u^{\prime }=\frac{\partial }{\partial u}(P+bn)=P_u+b_{u}n+b{n}_u$
• 假定凹凸函数$b$很小，可以忽略最后一项得到：$P_u^{\prime }\approx P_u+b_{u}n$
• 同理$P_v^{\prime }\approx P_v+b_{v}n$
• 扰动的表面法向量为:$N^{\prime }=P_u\times P_v+b_v(P_u\times n)+b_u(n\times P_v)+b_u{b}_v(n\times n)=P_u\times P_v+b_v(P_u\times n)+b_u(n\times P_v)$