16、渲染方程：定义、BRDF 及数值求解方法

d6e7f8

于 2025-09-17 11:56:30 发布

阅读量119

点赞数

CC 4.0 BY-SA版权

分类专栏：视觉效果的数学之美文章标签：渲染方程 BRDF 蒙特卡罗积分

本文链接：https://blog.youkuaiyun.com/d6e7f8/article/details/151980291

视觉效果的数学之美专栏收录该内容

21 篇文章 ¥499.90

订阅专栏¥69.90

会员秒杀 ¥9.9 重磅福利

超级会员免费看

渲染方程：定义、BRDF 及数值求解方法

1. 渲染方程的另一种定义

在某些情况下，原有的渲染方程可能并不适用，此时需要一个不同但等价的版本。新的渲染方程表达式如下：
[L(x, v) = \int_{A} f(x, u, v)L(x, u)V (x, x’)G(x, x’) dA_{x’} + L_e(x, v)]
其中：
- (x \in A) 是场景中接收来自 (x’ \in A) 点光线的点，(A) 表示场景中的所有表面。
- (V (x, x’))：若 (x) 和 (x’) 之间的路径无场景表面阻挡，其值为 1；否则为 0。
- (G(x, x’)) 是一个几何项，定义为：
[G(x, x’) = \frac{\langle u, n\rangle\langle u’, n’\rangle}{|x’ - x|^2}]
这里：
- (u \in R^3) 定义为 (u = \frac{x’ - x}{|x’ - x|})。
- 向量 (n \in S_x) 是点 (x) 处表面的单位法向量。
- (u’ \in R^3) 定义为 (u’ = \frac{x - x’}{|x - x’|})。
- (n’ \in S_{x’}) 是点 (x’) 处表面的单位向量。

下面用流程图展示光线传输中这些几何元素的关系：

graph LR
    A[x 点] -->|光线传输| B[x' 点]
    A -->|单位法向量| C[n]
    B -->|单位法向量| D[n']
    A -->|向量 u| E[u]
    B -->|向量 u'| F[u']
    E & F -->|计算| G[G(x, x')]
    A & B -->|判断| H[V (x, x')]
    G & H -->|参与积分| I[渲染方程]

2. 双向反射分布函数（BRDF）示例

我们主要考虑三种类型表面的 BRDF：

2.1 完美朗伯表面 BRDF

朗伯表面会将所有入射光均匀地散射到半球上，没有任何优先方向。这意味着 BRDF 在整个半球上是恒定的。因此，无论观察者从哪个方向观察，看到的结果都是相同的。朗伯表面的行为由光穿透表面边界后在表面上的散射来定义。

2.2 完美镜面 BRDF

完美镜面表面的 BRDF 定义为：
[f(x, u, v) = \frac{k_s\delta(r - v)}{\langle n, u\rangle}]
其中 (r \in S) 是该点半球上的完美镜面反射，反射方向 (r = u + 2n\langle n, u\rangle)，(k_s) 指定了反射光的比例。

2.3 修改后的 Blinn - Phong BRDF

基于修改后的 Blinn - Phong 模型的 BRDF 由两部分线性组合而成：
- 漫反射项：对应于完美朗伯分量。
- 镜面反射：可看作光遇到不同折射率介质边界时被表面反射的光，具有由表面法线定义的优先方向。

定义 (h(u, v) \in S) 为 (h(u, v) = \frac{(v + u)}{|(v + u)|})，其中 (v \in S) 是输出方向，(u \in S) 是光的入射方向。则修改后的 Blinn - Phong 模型定义为：
[f(x, u, v) = k_s\langle n, h(u, v)\rangle^p + k_d]
其中 (p \in N) 控制表面光的散射程度。增大 (p) 时，光在反射方向上会更加集中。

下面用表格总结三种 BRDF 的特点：
| BRDF 类型 | 特点 | 公式 |
| ---- | ---- | ---- |
| 完美朗伯表面 | 光均匀散射，BRDF 恒定 | 无特定公式，BRDF 为常数 |
| 完美镜面 | 遵循镜面反射规律 | (f(x, u, v) = \frac{k_s\delta(r - v)}{\langle n, u\rangle}) |
| 修改后的 Blinn - Phong | 漫反射和镜面反射组合 | (f(x, u, v) = k_s\langle n, h(u, v)\rangle^p + k_d) |

3. 数值近似

在实际应用中，渲染方程通常使用数值方法求解，常见的是蒙特卡罗积分法。蒙特卡罗积分的基本思想是从积分域中随机采样点，并使用这些点处被积函数的值来估计积分。

对于渲染方程 (L(x, v) = \int_{A} f(x, u, v)L(x, u)V (x, x’)G(x, x’) dA_{x’} + L_e(x, v))，使用蒙特卡罗积分估计积分的步骤如下：
1. 从集合 (A) 中随机采样 (N) 个点 (x’ 1, x’_2, …, x’_N)。
2. 计算这些采样点处被积函数的值 (f(x, u, v)L(x, u)V (x, x’_i)G(x, x’_i))。
3. 积分估计值为：
[I \approx \frac{1}{N} \sum {i = 1}^{N} f(x, u, v)L(x, u)V (x, x’_i)G(x, x’_i) \cdot |A|]
其中 (|A|) 是集合 (A) 的面积。

蒙特卡罗估计的准确性取决于采样点数 (N)。随着 (N) 的增加，估计值会收敛到积分的真实值，但计算成本也会随之增加。因此，需要在准确性和计算成本之间进行权衡。

下面用流程图展示蒙特卡罗积分估计渲染方程积分的过程：

graph LR
    A[渲染方程积分] --> B[随机采样 N 个点 x' 从 A]
    B --> C[计算采样点处被积函数值]
    C --> D[计算平均值并乘以 |A|]
    D --> E[得到积分估计值 I]

4. 蒙特卡罗积分方法的改进

蒙特卡罗积分方法可以通过重要性采样进一步改进。重要性采样是一种通过从更可能对积分有贡献的分布中采样点来减少蒙特卡罗估计方差的技术。

设 (p(x’)) 是集合 (A) 上的概率密度函数（PDF），则重要性采样的积分估计为：
[I \approx \frac{1}{N} \sum_{i = 1}^{N} \frac{f(x, u, v)L(x, u)V (x, x’_i)G(x, x’_i)}{p(x’_i)} \cdot |A|]
选择合适的 PDF (p(x’)) 对重要性采样的有效性至关重要。一个好的 (p(x’)) 应该与被积函数相似，这样采样点更有可能对积分有贡献。

在渲染方程的上下文中，常见的 (p(x’)) 选择是场景中光源的概率密度函数。因为光源是场景中辐射度的主要贡献者，从光源分布中采样点可以减少蒙特卡罗估计的方差，提高渲染的准确性。

5. 路径追踪

路径追踪是一种使用蒙特卡罗积分求解渲染方程的渲染算法。其基本思想是从相机追踪光线穿过场景，沿途与表面和光源进行交互。具体步骤如下：
1. 从相机发射光线进入场景。
2. 在每个交互点，路径追踪算法随机采样一个方向让光线继续传播。选择特定方向的概率由该交互点表面的 BRDF 决定。
3. 算法递归地在选定方向上追踪光线，直到光线到达光源或离开场景。
4. 相机处的辐射度通过对所有到达相机的光线的贡献求和来估计。每条光线的贡献由沿光线路径的辐射度和光线到达相机的概率的乘积决定。

路径追踪相对于其他渲染算法有几个优点：
- 能更准确地处理复杂的光照效果，如全局光照和软阴影。
- 与其他基于蒙特卡罗的渲染算法相比，能产生噪声更少的高质量图像。

然而，路径追踪的计算成本可能很高，特别是对于包含大量物体和光源的场景。

综上所述，渲染方程的求解涉及到多种方法和技术，从不同的方程定义到各种 BRDF 模型，再到数值求解方法和渲染算法，每种都有其特点和适用场景。在实际应用中，需要根据具体需求和场景复杂度选择合适的方法。

6. 各方法的对比与选择

为了更清晰地了解不同方法在求解渲染方程时的特点，我们将渲染方程的不同定义、BRDF 类型、数值求解方法以及渲染算法进行对比，以下是详细的对比表格：
| 类别 | 具体方法 | 优点 | 缺点 | 适用场景 |
| ---- | ---- | ---- | ---- | ---- |
| 渲染方程定义 | 传统定义 | - | 某些情况下不适用 | - |
| | 另一种定义 (L(x, v) = \int_{A} f(x, u, v)L(x, u)V (x, x’)G(x, x’) dA_{x’} + L_e(x, v)) | 适用于传统定义不适用的情况 | - | 传统定义不适用的场景 |
| BRDF 类型 | 完美朗伯表面 | 无论观察方向如何结果相同，计算相对简单 | 无法模拟镜面反射等复杂效果 | 简单的漫反射场景 |
| | 完美镜面 | 能准确模拟镜面反射 | 只能处理镜面反射，适用范围窄 | 有镜面反射的场景 |
| | 修改后的 Blinn - Phong | 结合漫反射和镜面反射，可控制光的散射程度 | 计算相对复杂 | 既有漫反射又有镜面反射的场景 |
| 数值求解方法 | 蒙特卡罗积分 | 基本思想简单，可用于求解复杂积分 | 估计准确性依赖采样点数，计算成本可能高 | 复杂积分的求解 |
| | 蒙特卡罗积分（重要性采样） | 减少估计方差，提高准确性 | 选择合适的 PDF 较难 | 场景中有明显主要贡献因素（如光源）的情况 |
| 渲染算法 | 路径追踪 | 能处理复杂光照效果，图像质量高、噪声少 | 计算成本高 | 对光照效果要求高、允许较高计算成本的场景 |

在实际应用中，我们可以根据以下步骤选择合适的方法：
1. 确定场景特点：分析场景中是否有镜面反射、漫反射、复杂光照效果等。
2. 考虑计算资源：评估可用于渲染的计算资源，包括时间和硬件性能。
3. 结合需求：根据对图像质量和渲染速度的要求进行权衡。

例如，如果场景是简单的漫反射场景，且对计算资源有限制，那么可以选择完美朗伯表面的 BRDF 和基本的蒙特卡罗积分方法；如果场景有复杂的光照效果且计算资源充足，那么路径追踪算法可能是更好的选择。

7. 技术细节深入分析

7.1 渲染方程中几何项 (G(x, x’)) 的作用

几何项 (G(x, x’) = \frac{\langle u, n\rangle\langle u’, n’\rangle}{|x’ - x|^2}) 考虑了光线在两点之间传输时的几何关系。其中 (\langle u, n\rangle) 和 (\langle u’, n’\rangle) 分别表示光线方向与表面法向量的夹角余弦，它们影响着光线与表面的交互强度。(|x’ - x|^2) 则体现了光线传播距离对辐射度的衰减。例如，当光线垂直照射表面时，(\langle u, n\rangle) 接近 1，光线与表面的交互较强；而当光线传播距离较远时，(|x’ - x|^2) 较大，辐射度会相应减小。

7.2 BRDF 中参数的影响

完美镜面 BRDF 中的 (k_s) ：(k_s) 指定了反射光的比例。当 (k_s = 1) 时，表示所有入射光都被反射，形成完美的镜面反射；当 (k_s < 1) 时，部分光被吸收，反射光强度减弱。
修改后的 Blinn - Phong BRDF 中的 (p) ：(p) 控制表面光的散射程度。当 (p) 较小时，光的散射范围较广，镜面反射效果较柔和；当 (p) 增大时，光在反射方向上更加集中，镜面反射效果更尖锐。

7.3 蒙特卡罗积分中采样点数的选择

蒙特卡罗积分的准确性与采样点数 (N) 密切相关。在实际操作中，可以通过以下步骤选择合适的 (N)：
1. 进行初步测试：选择一个较小的 (N) 值进行渲染，观察图像的噪声水平和估计的准确性。
2. 逐步增加 (N)：根据初步测试的结果，逐步增加 (N) 值，直到图像质量满足要求或达到计算资源的限制。
3. 权衡准确性和计算成本：在增加 (N) 提高准确性的同时，要考虑计算成本的增加。可以通过比较不同 (N) 值下的渲染时间和图像质量来找到最佳平衡点。