PBR物理渲染：Cook-Torrance与Disney模型

可视化演示

展示PBR能量守恒、Cook-Torrance微表面BRDF三函数、GGX粗糙度参数及Fresnel-Schlick近似效应，侧边指标实时切换深入理解物理渲染核心原理。

概念起源：从「凭手感调灯」到「科学布光」

老式摄影棚里，灯光师拍一瓶香水，靠的是几十年的手感：高光太亮就减一盏灯，反光太硬就加一张柔光纸。这套打法在固定场景里没问题，可一旦把香水搬到户外、换到黄昏，之前调好的参数全部作废，还得重新凭感觉来。传统的 Phong 等经验光照模型就是这样——它们是美术师手调的「临时配方」，看着差不多，但物理上经不起推敲：光可能被无中生有地放大，也可能在多个光源下叠出荒谬的亮度。

PBR 的核心理念，是给渲染世界装上一台「物理测光表」。它不再问「这个材质看起来应该多亮」，而是从真实世界的光学规律出发，建立一套自洽的数学描述。其中最重要的法则就是能量守恒：一束光打到表面，要么被反射出去，要么被吸收转化成热量，反射能量加上吸收能量永远等于入射能量，绝不允许光凭空变多。

反射方程（渲染方程）就是这台测光表的核心读数：

有了能量守恒和反射方程这两块基石，灯光师就不必再「盲调」了。只要把材质的 BRDF 描述准确，无论换到烈日还是烛光下，渲染引擎都能算出物理上可信的亮度。测光表已经就位，接下来要把「材质加工说明书」写清楚，这就是微表面理论登场的时候。

微表面理论：把表面放大一万倍看

摄影师想拍出金属拉丝和哑光塑料的区别，靠的是灯光技巧，但本质原因是两种材料表面的微观结构完全不同。请你想象：任何物体表面如果放大一万倍，都不再是光滑的平面，而是布满沟壑与凸起，像一地破碎的镜片——每一块小镜片都是一个「微面元」（微表面元素），它们是完美的微观镜面，各自朝向不同的方向。

微表面理论的核心假设正是如此：宏观表面由无数微小的完美镜面组成。一束平行光打上去，有的微面元朝左反射、有的朝右反射，宏观上我们看到的就是一团模糊的反射。如果所有微面元的朝向都差不多，反射就聚集成一个清晰明亮的高光点；如果朝向七零八落，反射就散成一片柔和的哑光。这里的关键词是「统计」——宏观观感是海量微面元方向分布的统计结果。

这个统计由法线分布函数 D 来描述。D 是一个概率密度，它在整个半球上必须满足归一化条件：

微表面理论让「粗糙度」从一句模糊的口头禅变成有精确数学含义的物理量：粗糙度越大，D 的分布越扁平，光被散射得越开；粗糙度越小，D 的分布越尖锐，高光如同一颗钉子一样锐利。现在，测光表和显微镜都准备好了，是时候请出 PBR 的核心人物——Cook-Torrance 模型。

Cook-Torrance 模型：灯光师的三件法宝

当摄影师布光时，桌上通常摆着三件法宝：柔光箱负责铺出均匀的底光，硬聚光灯负责打出锐利的轮廓，还得有反光板在暗部补一点细节。Cook-Torrance 模型也把材质对光的回应拆成了类似的三个部分，装进同一个框架——BRDF。它优雅地分开了「底光」和「高光」，用一套统一的物理公式描述从塑料到钢铁几乎所有材质。

BRDF 就是前面反射方程里的 f_r，它定量描述「一束光从 ω_i 方向打来，有多少能量被反射到 ω_o 方向」。Cook-Torrance 模型给出的 BRDF 分为两项：漫反射项 f_lambert 加上镜面反射项。后者是真正的重头戏，由 D、G、F 三个因子相乘，再除以一个几何校正分母：

这三件法宝各司其职，但又紧密咬合。缺少任何一件，高光要么亮得不真实，要么在掠射角产生违和的能量激增。接下来，我们把镜头推向 D、G、F 三位主演的特写，看看它们各自是如何工作的。

三项核心因子：GGX、Smith、Schlick 三兄弟

假如给 D、G、F 三兄弟分别贴一个标签，D 是「调度员」、G 是「保安」、F 是「变色龙」。调度员 D（法线分布函数）负责统计有多少微面元恰好朝向半向量 h；保安 G（几何遮蔽函数）负责检查这些微面元有没有被邻近凸起挡住，挡住的部分要把光扣下；变色龙 F（菲涅尔项）则随视角变化改变颜色和强度——你斜着看水面时，反光刺眼得睁不开眼，正是它在作怪。

D 的常见实现有 Beckmann 和 GGX 两种，其中 GGX（Trowbridge-Reitz）因为高光带有漂亮的「长尾」——从中心到边缘柔和衰减，而非生硬截断——成为游戏和影视业界的事实标准。GGX 的数学形式如下：

G 三弟的常见实现是 Smith 模型，它把遮蔽和阴影分成两个方向分别计算：视线方向的遮蔽 $G_1(\omega_o)$ 和光线方向的阴影 $G_1(\omega_i)$，两者相乘得到总的几何衰减。F 三弟则常用 Schlick 近似——一个成本极低却相当准确的公式：$F_{Schlick} = F_0 + (1-F_0)(1-n\cdot\omega_i)^5$，其中 F0 是垂直入射时的反射率，非金属约 0.04，金属则由 baseColor 直接决定。三兄弟协作，任何材质在不同光线、不同视角下，都能给出物理可信的高光表现。

Disney原则化：把物理公式变成「傻瓜式旋钮」

物理学家给出的 Cook-Torrance 公式虽然精确，但对美术师来说就像一本天书：D 用什么分布？G 用哪个模型？F0 取多少？每个参数背后都是一大串数学。Disney 公司在 2012 年提出了「原则化 BRDF」（Principled BRDF），目标是把物理公式封装成调音台上的旋钮——美术师不需要理解 GGX 和 Smith，只需要扭动几个直觉化的参数，就能得到物理正确的材质。

Disney 模型的精妙之处，在于用少量语义明确的参数驱动底层复杂的物理计算：baseColor 决定材质的基础色，metallic 在「非金属/金属」两个通道间平滑过渡，roughness 控制微表面粗糙度，其余还有 IOR（折射率）、clearcoat（清漆层）等附加旋钮。每个旋钮的物理意义明确，互不干扰。它的漫反射项长这样：

Disney 模型的推广，让「物理正确」不再是研究员的专利。美术师调整材质时，凭借直觉拧旋钮即可获得可信结果，团队协作效率大幅提升——这正是它成为行业标准的原因。旋钮已经就位，下一节我们去片场看看这套流程如何真正运转。

PBR实践流程：从摄影棚到工业流水线

还记得那位拿着测光表的灯光师吗？他一个人的专业能力，如今被做成了标准化的工业流程：剧组里每个美术师都用同一套旋钮、同一套贴图规范，拍出来的画面自然风格统一。在 PBR 工作流中，材质不再是一堆意义不明的数值，而是几张语义明确的贴图：Albedo（基础色）、Normal（法线）、Metallic（金属度）、Roughness（粗糙度），四图合一就能驱动任意 PBR 渲染器。

整个流程中，最容易出错的是工作空间。渲染必须在「线性空间」中计算——因为物理光照是线性叠加的，而显示器需要 sRGB 非线性输出。PBR 管线因此包含一个关键公式：金属度工作流中，材质的基础反射率 F0 由 baseColor 和 metallic 线性插值得到：

生产流程最后一步是光照环节：PBR 场景需要配合 HDR（高动态范围）环境贴图，提供远超 0~1 范围的真实光照能量，再经过曝光和色调映射输出到屏幕。至此，从贴图生产、线性空间计算到 HDR 光照，一条完整的 PBR 流水线就搭建完毕了——灯光师已经退居幕后，规则替他完成了所有辛苦的推断。

AI 应用场景

Python 演示

import math W, H = 200, 100CX, CY, R = 100, 50, 46  def normalize(x, y, z):    length = math.sqrt(x * x + y * y + z * z)    return x / length, y / length, z / length  def schlick_fresnel(f0, ndv):    return f0 + (1.0 - f0) * math.pow(1.0 - ndv, 5.0)  def ggx_ndf(ndh, roughness):    alpha = roughness * roughness    a2 = alpha * alpha    denom = ndh * ndh * (a2 - 1.0) + 1.0    return a2 / (math.pi * denom * denom)  def smith_g(ndl, ndv, roughness):    alpha = roughness * roughness    k = alpha * 0.5    g1v = ndv / (ndv * (1.0 - k) + k)    g1l = ndl / (ndl * (1.0 - k) + k)    return g1v * g1l  def cook_torrance(nx, ny, nz, vx, vy, vz, lx, ly, lz, base_color, roughness):    ndl = max(nx * lx + ny * ly + nz * lz, 0.0)    if ndl <= 0.0:        return 0.0, 0.0, 0.0    ndv = max(nx * vx + ny * vy + nz * vz, 0.0)    hx, hy, hz = normalize(lx + vx, ly + vy, lz + vz)    ndh = max(nx * hx + ny * hy + nz * hz, 0.0)    d = ggx_ndf(ndh, roughness)    g_smith = smith_g(ndl, ndv, roughness)    f = schlick_fresnel(0.04, ndv)    spec = d * g_smith * f / (4.0 * ndv * ndl + 1e-4)    diff = (1.0 - f) * ndl / math.pi    red = base_color[0] * diff + spec    grn = base_color[1] * diff + spec    blu = base_color[2] * diff + spec    return red, grn, blu  def main():    lx, ly, lz = normalize(0.7, -0.4, 0.9)    vx, vy, vz = 0.0, 0.0, 1.0    base = (0.62, 0.34, 0.72)    rough = 0.28    pixels = bytearray()    for y in range(H):        for x in range(W):            nx = (x - CX) / R            ny = (y - CY) / R            m2 = nx * nx + ny * ny            if m2 > 1.0:                pixels += bytes((28, 30, 44))                continue            nz = math.sqrt(1.0 - m2)            r, g, b = cook_torrance(nx, ny, nz, vx, vy, vz, lx, ly, lz, base, rough)            pixels += bytes((int(min(r, 1.0) * 255), int(min(g, 1.0) * 255), int(min(b, 1.0) * 255)))    with open("pbr_ball.ppm", "wb") as f:        f.write(b"P6\n%d %d\n255\n" % (W, H))        f.write(bytes(pixels))  if __name__ == "__main__":    main()

1. d = ggx_ndf(ndh, roughness)计算有多少微面元的法线恰好对齐半向量，这是高光形成的核心。
2. g_smith = smith_g(ndl, ndv, roughness)用 Smith 模型计算微表面凹凸互相遮光的衰减。
3. f = schlick_fresnel(0.04, ndv)用 Schlick 近似算出菲涅尔反射率，0.04 是非金属的垂直反射基线。
4. spec = d * g_smith * f / (4.0 * ndv * ndl + 1e-4)三项因子相乘再除以分母，得到物理可信的镜面高光。
5. diff = (1.0 - f) * ndl / math.pi漫反射要扣除被菲涅尔反射「截胡」的能量，保证总能量守恒。
6. nx = (x - CX) / R把每个像素坐标换算成球面局部坐标，用来推出球面法线方向。
7. nz = math.sqrt(1.0 - m2)球面方程解出法线的 z 分量，m2 大于 1 说明像素在球外。
8. normalize(0.7, -0.4, 0.9)把光源方向缩放到单位长度，后续点积计算才成立。
9. f.write(b"P6\n%d %d\n255\n" % (W, H))写出 PPM 二进制格式的文件头，声明图片尺寸和颜色深度。

常见误区

面试题