光照模型和 BRDF

原创

本文参考了本篇博客，同时加入了 Phong 光照模型

一、预备知识

在介绍光照模型和 BRDF（双向反射分布函数）时，我们要先理解一些基本概念，这些概念与辐射度学有关

1. 能量

每个光子都具有一定的能量（Energy），能量大小与频率有关\(E = hv\)。用符号\(Q\)表示，单位是焦耳\(J\)

2. 功率

==功率、辐射通量、通量==

功率（Power），也被称为辐射通量（Radiant Flux）或者通量（Flux），指的是单位时间内通过表面或者空间区域的能量的总量。用符号\(\Phi\)表示，\(\Phi = \frac{dQ}{dt}\)，单位为瓦特\(W\)，或者焦耳/秒\(J/S\)

3. 辐照度和辐出度

==辐照度（辐射通量密度）==

==辐出度、辐射出射度、辐射度（辐射通量密度）==

• 辐照度（Irradiance）指单位时间内到达单位面积的辐射能量，或单位面积的辐射通量，即通量对于面积的密度。用符号\(E\)表示，\(E = \frac{d\Phi}{dA} = \frac{dQ}{dtdA}\)，单位\(W/m^{2}\)
• 辐出度（Radiant Existance），也称为辐射出射度、辐射度（Radiosity）。与辐照度不同的是，辐出度是衡量离开表面的通量密度，二者都可以被称为辐射通量密度（Radiant Flux Density）。用符号\(M\)表示

注意辐射通量和辐射通量密度之间的区别

如下图所示

一束平行光之间的间距为\(d\)，一个垂直入射，一个与表面法线成\(\theta_{i}\)的角度入射，观察表面的间距，其中左图间距为\(d\)，右图为\(\frac{d}{cos\theta_{i}}\)，相对于左图，右图的间距增大了。由于辐射通量是单位时间内通过表面或空间区域内的能量总和，所以辐射通量相同；而辐射通量密度是单位时间内通过单位面积的能量，右边面积大，所以右边辐射通量密度变小

让我们来看一下二者辐射通量密度的差异：假设右图中不垂直于光线方向的表面面积为\(A\)，将他投影到垂直于直线方向得到一个虚拟表面，面积为\(A^{\perp } = Acos\theta_{i}\)，通过这两个面积的通量相同，设为\(\Phi\)，则表面接收到的辐射通量密度\(E = \frac{\Phi}{A}\)，虚拟表面的辐射通量密度\(E_{\perp} = \frac{\Phi}{A^{\perp}} = \frac{\Phi}{Acos\theta_{i}}\)，得到\(E = E_{\perp}cos\theta_{i}\)

同理对于点光源，假想以点光源为中心不同半径的求包围着点光源，穿过这些球的辐射通量是相通的，均为\(\Phi\)，球的表面积为\(4\pi r^{2}\)，得到辐射通量密度\(E = \frac{\Phi}{4\pi r^{2}}\)，即通量密度与距离的平方成反比，即光的衰减与距离的平方成反比

4. 辐射强度

==辐射强度==

立体角：度量三维角度的量，用符号\(\omega\)表示，立体角计算为\(\omega = \frac{s}{r^{2}}\)，单位球的表面积为\(4\pi r^{2} = 4\pi\)，所以球面的立体角也为\(4\pi\)

辐射强度（Radiant Intensity），通过单位立体角的辐射能量，用符号\(I\)表示，\(I=\frac{d\Phi}{dw}\)，单位\(W/sr\)

引入辐射强度的原因是：有时候需要度量通过一个点的通量的密度，因为点的面积是 0，无法直接使用辐射通量密度，所以引入辐射强度。辐射强度不会随距离的增大而衰减，这是因为立体角不会随距离的变化而变化

5. 辐射率

==辐射率==

辐射率（Radiance），每单位面积每单位立体角的辐射通量密度，用符号\(L\)表示，\(L = \frac{d\Phi}{d\omega dA^{\perp}}\)，单位为\(W/m^{2}sr\)，其中\(dA^{\perp}\)是微分面积\(dA\)在垂直于光线方向的投影，如图

辐射率可以看成我们眼睛看到（相机拍到）的物体上一点的颜色，基于物理着色时，计算表面一点的颜色就是计算它的辐射率，辐射率不会随距离变化而衰减，这和我们日常感受一致，在没有雾霾的干扰时，我们看到的物体表面上一点的颜色并不会随距离变化而变化。为什么辐照度会随距离增大而衰减，但是我们看到的颜色却不会衰减呢？这是因为随着距离变大，我们看到的物体上的一块区域到达视网膜的通量密度会变小，同时这块区域在视网膜表面上的立体角也会变小，正好抵消了通量密度的变化

二、BRDF

我们看到一个表面，实际上是周围环境的光照射到表面上，然后表面将一部分光反射到我们眼睛里。 双向反射分布函数 BRDF（Bidirectional Reflectance Distribution Function） 就是描述表面入射光和反射光关系的

对于一个方向的入射光，表面会将光反射到表面上半球的各个方向，不同方向反射的比例是不同的，我们用 BRDF 来表示指定方向的反射光和入射光的比例关系，BRDF 定义为：

\[ f(l, v) = \frac{dL_{o}(v)}{dE(l)} \]

• \(f\)：BRDF，\(l\)是入射光方向，\(v\)是观察方向，即反射光方向
• \(dL_{o}(v)\)：表面反射到\(v\)方向的反射光的微分辐射率。表面反射到\(v\)方向的反射光的辐射率为\(L_{o}(v)\)，来自于表面上半球所有方向的入射光线的贡献，而微分辐射率\(dL_{o}(v)\)特指来自方向\(l\)的入射光贡献的反射辐射率
• \(dE(l)\)：表面上来自入射光方向\(l\)的微分辐照度。表面接收到的辐照度为\(E\)，来自上半球所有方向的入射光线的贡献，而微分辐照度\(dE(l)\)特指来自于方向\(l\)的入射光

表面对不同频率的光反射率可能不一样，因此BRDF和光的频率有关。在图形学中，将BRDF表示为RGB向量，三个分量各有自己的\(f\)函数

BRDF需要处理表面上半球的各个方向，如下图使用球坐标系定义方向更加方便。球坐标系使用两个角度来确定一个方向：

• 方向相对法线的角度\(\theta\)，称为极角（Polar Angle）或天顶角（Zenith Angle）
• 方向在平面上的投影相对于平面上一个坐标轴的角度\(\phi\)，称为方位角（Azimuthal Angle）

所以BRDF也可以表示成\(f(\theta_i, \phi_i, \theta_o, \phi_o)\)。对于各向同性材质，当\(l\)和\(v\)同时绕法线\(n\)旋转时，\(f\)值保持不变，此时可以用\(l\)和\(v\)在平面投影的夹角\(\phi\)来代替\(\phi _i\)和\(\phi _o: f(\theta_i, \theta_o, \phi)\)

为什么BRDF要定义成辐射率和辐照度的比值，而不是直接定义为辐射率和辐射率比值，有两种解释

==第一种解释参考BRDF为什么要定义为一个单位是sr-1的量？==

结合下面辐射通量密度\((A)\)和辐射率\((B)\)测量仪的示意图来看：辐照度(辐射通量密度)测量仪\((A)\)接受平面上半球的所有光线，可以测量一个较小面积来自于四面八方的所有光通量，光通量\(\Phi\)除以传感器面积\(A\)就可以得到辐照度(辐射通量密度)\(E\)。辐射度测量仪\((B)\)则有一个长筒控制光线只能从一个很小的立体角进入测量仪，光通量\(\Phi\)除以传感器面积\(A\)和立体角\(\omega\)就可以得到辐射率\(L\)

测平面上一点在某一个方向的出射辐射率很简单，只需要用仪器\((B)\)从该方向对准该点就可以了。而测平面一点入射的辐射率则没有那么简单，必须保证光源正好覆盖测量仪开口立体角，大了该点会接受到比测量值更多的光照，导致测量值比实际值小，小了则与仪器的设计立体角不一致，可在实际中是基本做不到光源大小正好覆盖测量仪开口立体角的。而测表面的辐照度则简单得多，只要保证光源很小，而且没有来自其他方向的光干扰，这时候测到的辐照度就是平面上来自光源方向的微分辐照度\(dE\)

==第二种解释从数学的角度出发==

对于现实世界中的非光学平面，一束光线射到表面上后，被表面反射到各个方向，其中一个出射方向的光通量只是整个反射光通量极小的一部分，当出射方向立体角趋于0时，\(\lim_{\omega_{o} \rightarrow 0}{\frac{dL_{o} }{Li} } = 0\)，所以在实际计算中使用辐射率和辐射率比值是没有意义的。而如果分母改成表面上接收到的来自光源方向的微分辐照度，我们知道\(dE = L_i(l) d\omega _{i} cos \theta _{i}\)，由于给入射辐射率乘了一个趋于零的微分立体角，\(dE\)的值会小很多，比值\(\frac{dL_o}{dE}\)是有意义的，而不是0

下面我们来看看怎么用BRDF来计算表面辐射率

我们考虑来自方向\(l\)的入射光辐射率\(L_i(l)\)，由辐射率和辐照度的定义：

\[ L_i(l) = \frac{d \Phi }{d\omega_i d A^{\bot } } = \frac{d \Phi }{d\omega_i dA cos \theta_i } = \frac{dE(l) }{d\omega_i cos \theta_i } \]

则照射到表面来自于方向l的入射光贡献的微分辐照度：

\[ dE(l) = L_i(l) d\omega_i cos \theta_i \]

表面反射到\(v\)方向的由来自于方向l的入射光贡献的微分辐射率：

\[ dL_o(v) = f(l, v) \otimes dE(l) = f(l, v) \otimes L_i(l) d\omega_i cos \theta_i \]
符号\(\otimes\)表示按向量的分量相乘，因为\(f\)和\(L_i\)都包含RGB三个分量

要计算表面反射到v方向的来自上半球所有方向入射光线贡献的辐射率，可以将上式对半球所有方向的光线积分：

\[ L_o(v) = \int_{\Omega }^{} f(l, v) \otimes L_i(l) cos \theta_i d\omega_i \]
上式称为反射方程（Reflectance Equation），用来计算表面反射辐射率

对于点光源、方向光等理想化的精准光源（Punctual Light），计算过程可以大大简化。我们考察单个精准光源照射表面，此时表面上的一点只会被来自一个方向的一条光线照射到（而面积光源照射表面时，表面上一点会被来自多个方向的多条光线照射到），则辐射率：

\[ L_o(v) = f(l, v) \otimes E_L cos \theta_i \]
对于多个精准光源，只需简单累加就可以了：

\[ L_o(v) = \sum_{k = 1}^{n}{f(l_k, v) \otimes E_{L_k} cos \theta_{i_k}} \]
这里使用光源的辐照度，对于阳光等全局方向光，可以认为整个场景的辐照度是一个常数，对于点光源，辐照度随距离的平方衰减，用公式\(E_{L} = \frac{\Phi }{4\pi r^2}\)就可以求出到达表面的辐照度，\(\Phi\)是光源的功率，比如100瓦的灯泡，r是表面离光源的距离

回头看看反射方程，是对表面上半球所有方向的入射光线积分，这里面包含了来自精准光源的光线，也包括周围环境反射的光线。处理来自周围环境的光线可以大幅提高光照的真实程度，在实时图形学中，这部分光照可以用基于图像的光照（Image Based Lighting）来模拟。我们将在下篇文章讨论基于图像的光照。

上面给出了BRDF的定义和使用BRDF计算表面反射辐射率的公式。但这个定义实际上是无法直接用于计算表面反射辐射率的，我们还要建立一个能模拟真实光照的模型，使得输入入射方向和出射方向，\(f(l, v)\)能输出表面反射微分辐射率和入射微分辐照度的比率。

1967年Torrance-Sparrow在Theory for Off-Specular Reflection From Roughened Surfaces中使用辐射度学和微表面理论建立了模拟真实光照的BRDF模型，1981年Cook-Torrance在A Reflectance Model for Computer Graphics中把这个模型引入到计算机图形学领域，现在这个模型已经成为基于物理着色的标准，被称为Cook-Torrance模型。下面我们来看看微表面理论和Cook-Torrance模型的推导过程。

三、Phong光照模型

光照到物体表面时，物体对光会发生反射、透射和吸收。简单光照模型只考虑物体对直接光照的反射作用，而物体间的光反射作用，只用环境光来表示。此时光源被假定为点光源，反射作用被细分为镜面反射和漫反射

==Phong光照模型==

根据Phong光照模型，物体被感知的亮度由环境光、漫反射光及镜面反射光组成，物体表面的反射光强度为
\[ I = K_{a}I_{a} + K_{d}I_{l}cos(\theta) + K_{s}I_{l}cos^{n}(\alpha) \]

• \(K_{a}~K_{d}~K_{s}\)：分别为环境光、漫反射光及镜面反射光系数
• \(I_{a}\)：入射的环境光强度
• \(I_{l}\)：光源的入射光强度
• \(\theta\)：入射光与物体表面法线之间的夹角
• \(\alpha\)：反射光与视线之间的夹角
• \(n\)：镜面反射参数

1. 漫反射

理想漫反射当光源来自一个方向时，漫反射光均匀向各方向传播，与视点无关，他是由表面的粗糙不平引起的，因而漫反射光的空间分布是均匀的。记入射光强为\(I_{p}\)，物体表面上点\(P\)的法向量为\(N\)，从点\(P\)指向光源的向量为\(L\)，当\(L, N\)为单位向量时，理想漫反射光强表示为：
\[ I_{d} = I_{p}K_{d}(L\cdot N) \]

\(K_{d}\)：与物体有关的漫反射系数，取值\((0, 1)\)。在RGB模型下，漫反射系数\(K_{d}\)有三个分量\(K_{dr}, K_{dg}, K_{db}\)分别表示RGB三原色的漫反射系数，反应物体颜色，通过调整三个分量可以改变物体的颜色，也可以吧入射光强\(I\)设为三个分量\(I_{r}, I_{g}, I_{b}\)通过调整这些分量来调整光源的颜色

2. 镜面反射

镜面反射光对于理想镜面，反射光集中在一个方向，并遵守反射定律。对一般光滑表面，反射光集中在一个范围内，且由反射定律决定的反射方向光强最大。因此对于同一点来说，从不同位置所观察到的镜面反射光强是不同的，镜面反射光强可以表示为：
\[ I_{s} = I_{p}K_{s}cos^{n}\alpha~~~\alpha \in (0, \frac{\pi}{2}) \]

• \(K_{s}\)：与物体有关的镜面反射系数
• \(\alpha\)：视线方向\(V\)与反射方向\(R\)的夹角
• \(n\)：反射指数，反映了物体表面的光泽程度，一般为1~2000，数目越大表示物体表面越光滑

镜面反射光将会在反射方向附近形成很亮的光斑，称为高光现象，将\(V\)和\(R\)看做单位向量，镜面反射光强可以重写为：
\[ I_{s} = I_{p}K_{s}(R\cdot V)^{n} \]

其中\(R = 2N(N\cdot L)-L\)，镜面反射光产生的高光区域只反映光源的颜色，例如在红光照射下，一个物体的高光域是红光，镜面反射系数\(K_{s}\)是一个与物体颜色无关的参数，在简单光照明模型中，只能通过改变物体的漫反射系数来控制物体的颜色

3. 环境光

环境光指光源间接对物体的影响，是在物体和环境之间多次反射，最终达到平衡时的一种光，可以近似地认为同一环境下的环境光，其光强分布是均匀的，他在任何一个方向上的分布都相同，在简单光照模型中，用一个常数来模拟环境光，用式子表示为：
\[ I_{e} = I_{a}K_{a} \]

本文由 Yonghui Wang 创作，采用 知识共享署名4.0 国际许可协议进行许可
本站文章除注明转载/出处外，均为本站原创或翻译，转载前请务必署名
最后编辑时间为: Apr 27, 2026 07:56 pm