9.5 电介质BSDF - Physically Based Rendering From Theory To Implementation(PBRT)

Physically Based Rendering From Theory To Implementation

翻译：https://pbr-book.org/4ed/Reflection_Models/BSDF_Representation

总目录： Physically Based Rendering From Theory To Implementation - 基于物理的渲染从理论到实践第四版原书翻译

9.5 Dielectric BSDF 电介质BSDF

在介质情况下，相对折射率是实值，除了反射外还必须考虑镜面透射。DielectricBxDF处理平滑和粗糙接口的这种情况。

图9.14显示了使用BxDF对玻璃材料建模的抽象模型图像。

图9.14:当镜面反射的BRDF和镜面透射的BTDF用 Fresnel 电介质公式调制时，反射和透射量的实际角度相关变化给出了玻璃的视觉精确表示。(场景由Beeple提供)

<<DielectricBxDF Definition>>= 
class DielectricBxDF {
  public:
    <<DielectricBxDF Public Methods>> 
  private:
    <<DielectricBxDF Private Members>> 
};

构造函数接受一个float值的eta参数和一个microfacet分布mfDistrib。光谱变化的IOR将光分散到不同的方向，通过随机采样单个波长来处理，然后实例化相应的DielectricBxDF。第10.5.1节更详细地讨论了这个主题。

<<DielectricBxDF Public Methods>>= 
DielectricBxDF(Float eta, TrowbridgeReitzDistribution mfDistrib)
    : eta(eta), mfDistrib(mfDistrib) {}

<<DielectricBxDF Private Members>>= 
Float eta;
TrowbridgeReitzDistribution mfDistrib;

Flags()方法处理三种不同的情况。第一种是当介电界面是折射率匹配的——也就是说，两边的IOR相等(在这种情况下)——光只被折射。否则，在另外两种情况下，BSDF同时具有反射和折射部分。在这两种情况下，TrowbridgeReitzDistribution的 EffectivelySmooth() 方法区分了镜面散射（specular）和光滑散射（glossy）。

<<DielectricBxDF Public Methods>>+=  
BxDFFlags Flags() const {
    BxDFFlags flags = (eta == 1) ? BxDFFlags::Transmission :
                      (BxDFFlags::Reflection | BxDFFlags::Transmission);
    return flags | (mfDistrib.EffectivelySmooth() ? BxDFFlags::Specular
                                                  : BxDFFlags::Glossy);
}

Sample_f()方法必须在采样完美镜面反射或透射之间进行选择。像以前一样，我们推迟处理粗糙的表面，现在只讨论完美的镜面情况。

<<DielectricBxDF Method Definitions>>= 
pstd::optional<BSDFSample>
DielectricBxDF::Sample_f(Vector3f wo, Float uc, Point2f u,
        TransportMode mode, BxDFReflTransFlags sampleFlags) const {
    if (eta == 1 || mfDistrib.EffectivelySmooth()) {
        <<Sample perfect specular dielectric BSDF>> 
    } else {
        <<Sample rough dielectric BSDF>> 
    }
}

图9.15:使用DielectricBxDF渲染的玻璃物体。(a)在每次散射事件中以等概率在镜面反射和透射之间进行选择。(b)根据菲涅耳方程的值进行概率选择，如Sample_f()方法所实现。选择与贡献概率成正比的散射模式，通过选择较少的低贡献路径，可以显著降低误差。

由于介质具有反射和透射两种特性，采样方案必须在这两种成分之间随机选择，从而影响密度函数。虽然理论上可以接受任何离散分布，但从蒙特卡洛方差的角度出发，一种有效的方法是根据这两个分量的贡献进行采样——换句话说，与菲涅耳反射率R和互补透射率1-R成比例。图9 - 15展示了这种采样方式与均匀分配反射和透射的好处。

<<Sample perfect specular dielectric BSDF>>= 
Float R = FrDielectric(CosTheta(wo), eta), T = 1 - R;
<<Compute probabilities pr and pt for sampling reflection and transmission>> 
if (uc < pr / (pr + pt)) {
    <<Sample perfect specular dielectric BRDF>> 
} else {
    <<Sample perfect specular dielectric BTDF>> 
}

由于可以通过sampleFlags参数选择性地启用或禁用BSDF组件，考虑到这一点，组件选择是基于调整后的概率pr和pt。

<<Compute probabilities pr and pt for sampling reflection and transmission>>= 
Float pr = R, pt = T;
if (!(sampleFlags & BxDFReflTransFlags::Reflection)) pr = 0;
if (!(sampleFlags & BxDFReflTransFlags::Transmission)) pt = 0;
if (pr == 0 && pt == 0)
    return {};

在最常见的情况下，同时对反射和透射进行采样，BSDF值和采样概率包含公共因子R或T，当取它们的比值时，它们会抵消。因此，所有采样的光线最终都做出了相同的贡献，菲涅耳因子体现在反射和透射光线的相对比例上。

把所有这些放在一起，与类似的片段<<Sample perfect specular conductor BRDF>>相比，以下代码的唯一变化是样本的离散概率的合并。

读者补充： 反射光谱为反射率除以|Cos(wi)|，Fresnel反射率需要除以|Cos(wi)|才能得到正确的f值,因为我们在Fresnel提前计算了光传输方程的cos项。

<<Sample perfect specular dielectric BRDF>>= 
Vector3f wi(-wo.x, -wo.y, wo.z);
SampledSpectrum fr(R / AbsCosTheta(wi));
return BSDFSample(fr, wi, pr / (pr + pt), BxDFFlags::SpecularReflection);

镜面透射是沿着类似的线处理，虽然使用折射光线的方向为wi。 相应的BTDF方程类似于完美镜面反射的情况，方程(9.9)，虽然有一个额外的技术细节：根据IOR $\eta$，折射在角域中压缩或扩展辐射，实现必须缩放ft来考虑这一点。这个修正并不改变场景中的辐射功率的数量，而是模拟相同的功率如何包含在不同的立体角中。这一步的细节因双向渲染算法中传播的方向而异，因此我们将相应的片段<<Account for non-symmetry with transmission to different medium>>推迟到9.5.2节。

<<Sample perfect specular dielectric BTDF>>= 
<<Compute ray direction for specular transmission>> 
SampledSpectrum ft(T / AbsCosTheta(wi));
<<Account for non-symmetry with transmission to different medium>> 
return BSDFSample(ft, wi, pt / (pr + pt),
                  BxDFFlags::SpecularTransmission, etap);

函数Refract()通过斯涅尔定律计算折射方向wi，在全内反射的情况下，该定律失效。原则上，这种情况不应该发生:传输情况的采样概率为T，在全内反射情况下为零。但是，由于浮点舍入误差，这里偶尔会出现不一致的情况。我们通过退回无效样品来处理这个极端情况。

<<Compute ray direction for specular transmission>>= 
Vector3f wi;
Float etap;
bool valid = Refract(wo, Normal3f(0, 0, 1), eta, &etap, &wi);
if (!valid) return {};

与ConductorBxDF一样，如果界面是光滑的，并且所有散射都是完美的镜面反射，则f()方法返回零。

<<DielectricBxDF Method Definitions>>+=  
SampledSpectrum DielectricBxDF::f(Vector3f wo, Vector3f wi,
                                  TransportMode mode) const {
    if (eta == 1 || mfDistrib.EffectivelySmooth())
        return SampledSpectrum(0.f);
    <<Evaluate rough dielectric BSDF>> 
}

此外，如果 BSDF 使用 delta 分布表示，则返回 PDF 值为零。

<<DielectricBxDF Method Definitions>>+= 
Float DielectricBxDF::PDF(Vector3f wo, Vector3f wi, TransportMode mode,
          BxDFReflTransFlags sampleFlags) const {
    if (eta == 1 || mfDistrib.EffectivelySmooth())
        return 0;
    <<Evaluate sampling PDF of rough dielectric BSDF>> 
}

缺失的三个片段<<Sample rough dielectric BSDF>>、<<Evaluate rough dielectric BSDF>>和<<Evaluate sampling PDF of rough dielectric BSDF>>将在9.7节中介绍。

9.5.1 Thin Dielectric BSDF 薄电介质BSDF

电介质界面很少孤立地出现：一种特别常见的配置包括两个附近的折射率变化，它们是平滑的、局部平行的、互逆的——也就是说，具有相对IOR $\eta$ 和具有$1/\eta$逆的相应界面。例如，窗户上的平板和丙烯酸玻璃或用于密封和保存食物的塑料箔。

这种重要的特殊情况被称为薄电介质，因为与较大的横向尺寸相比，薄电介质两个界面的空间接近。当入射光分裂成具有两个独立界面的反射和透射组件时，它在一个递归过程中散射，该过程将一些光困在两个界面内(尽管随着反射次数的增加，该量逐渐衰减)。

虽然一般电介质内部的散射过程可能令人望而生畏，但简单的解析解可以完全描述在这样薄的电介质内部发生的事情-即满足上述简化条件的界面对。pbrt提供了一个名为ThinDielectricBxDF的专用BSDF，该BSDF利用这些见解来有效地表示无限数量的内部交互。它进一步允许这样的表面用单个界面表示，节省了在两个表面之间跟踪光线路径的射线相交费用。

<<ThinDielectricBxDF Definition>>= 
class ThinDielectricBxDF {
  public:
    <<ThinDielectricBxDF Public Methods>> 
  private:
    Float eta;
};

这个BxDF的唯一参数是内部介质的相对折射率。

<<ThinDielectricBxDF Public Methods>>= 
ThinDielectricBxDF(Float eta) : eta(eta) {}

由于这个BxDF只模拟完美的镜面散射，它的f()和PDF()方法都只返回零，因此不包括在这里。

图9.16:薄平面平行介质中的光路。入射光线(红色)在玻璃板(蓝色)内产生无限的内部散射过程。在每次散射事件中，都有一部分光设法逃逸，导致捕获的能量连续衰减。

薄电介质BSDF的理论可以追溯到Stokes(1860)的开创性工作，他研究了光在玻璃板堆中的传播。图9.16说明了仅涉及单个玻璃板的最常见情况：入射光线(红色)按比例R和T=1-R反射和透射。当发射的光线遇到底部界面时，互易性使其以相同的比例反射和发射。这个过程不断重复。

当然，第一次相互作用之外的光线相对于入口点会发生位移。由于假设电介质很薄，这种空间位移可以忽略不计;只有反射或透射光的总量才是重要的。通过这种简化，就有可能把无数散射事件的影响汇总成一个简单的反射率和透射率因子的修正。

考虑从顶层反射出来的路径。它们的总反射率 $R^{\prime }$ 由几何级数给出，可以转换成显式形式:

$$R^{\prime }=R+TRT+TRRRT+...=R+\frac{T^{2}R}{1-R^2}$$

一个类似的数列给出了光的透射量，但由于能量守恒，它可以很容易地计算出来$T^{\prime }=1-R^{\prime }$。图9.17绘制$R^{\prime }$和$R$随入射角$\theta$的函数图。与单一 Fresnel 相互作用相比，第二个界面具有增加反射总量的效果。

图9.17:Fresnel 界面和薄电介质的反射率。该图比较了由 Fresnel 方程(9.6)确定的单个介质界面的反射率与根据式(9.10)匹配的薄介质界面的反射率。

Sample_f()方法计算$R^{\prime }$和$T^{\prime }$系数，然后计算采样反射和传输的概率，就像DielectricBxDF所做的那样，重用相应的代码片段。

<<ThinDielectricBxDF Public Methods>>+= 
pstd::optional<BSDFSample>
Sample_f(Vector3f wo, Float uc, Point2f u, TransportMode mode,
         BxDFReflTransFlags sampleFlags) const {
    Float R = FrDielectric(AbsCosTheta(wo), eta), T = 1 - R;
    <<Compute R and T accounting for scattering between interfaces>> 
    <<Compute probabilities pr and pt for sampling reflection and transmission>> 
    if (uc < pr / (pr + pt)) {
        <<Sample perfect specular dielectric BRDF>> 
    } else {
        <<Sample perfect specular transmission at thin dielectric interface>> 
    }
}

更新后的反射系数和透射系数很容易使用式(9.10)计算，但必须注意避免在$R=1$的情况下被零除。

<<Compute R and T accounting for scattering between interfaces>>= 
if (R < 1) {
    R += Sqr(T) * R / (1 - Sqr(R));
    T = 1 - R;
}

在该方法的实现中也重用了对完美镜面反射进行采样的DielectricBxDF片段，继承了计算得到的R值。传输情况稍微偏离DielectricBxDF，因为传输方向是简单的否定。

<<Sample perfect specular transmission at thin dielectric interface>>= 
// 采样在薄介质界面上具有完美的镜面透射
// 入射光wi 由物体表面指向外侧，因此折射光为-wi，于实际入射光方向相同
Vector3f wi = -wo;
SampledSpectrum ft(T / AbsCosTheta(wi));
return BSDFSample(ft, wi, pt / (pr + pt), BxDFFlags::SpecularTransmission);

9.5.2 Non-Symmetric Scattering and Refraction 非对称散射和折射 ❗❗❗

图9.18:Snell的窗口。如果在游泳池里向上看，只有通过圆形窗户才能看到天空，因为光线在临界角度之外没有折射。窗外，只能看到池底的倒影。

所有基于物理的BRDFs都是对称的:入射方向和出射方向可以在不改变函数值的情况下互换。但对于BTDFs，情况通常并非如此。与BTDFs的非对称性是由于当光折射成比入射介质折射率更高的材料时，能量被压缩成更小的角度集合(反之亦然，当方向相反时)。这种效果很容易看到你自己，例如，从一个安静的室外游泳池的水下看天空。因为光线的折射不能超过临界角(对于水来说约48.6°)，所以入射的半球被压缩成一个相当小的半球子集，这个子集包含了剩下的角度(图9 - 18)。因此，沿着折射的射线的辐射度必须增加，这样当光通过界面时，能量才能保存下来。

图9.19：不同折射率介质之间的边界处透射光量用两种折射率的平方比进行缩放。直观上，这可以理解为辐射的微分立体角被压缩或扩大作为传输的结果。

更正式地，考虑入射辐射到达两个介质的边界，折射率为${\eta }_i$和${\eta }_o$(图9.19)。现在假设所有入射光都被透射，则透射的微分通量的量为：

$$d^2{\Phi }_o=d^2{\Phi }_i$$

如果我们使用辐照度的定义，公式(4.3)，我们等效有：

$$L_{o}cos{\theta }_{o}dAd{\omega }_o=L_{i}cos{\theta }_{i}dAd{\omega }_i$$
注：原文这里${\omega }_i$应该为$d{\omega }_i$，应该是笔误

将立体角扩展为球面角得到

$$L_{o}cos{\theta }_{o}dAsin{\theta }_{o}d{\theta }_{o}d{\varphi }_o=L_{i}cos{\theta }_{i}dAsin{\theta }_{i}d{\theta }_{i}d{\varphi }_i，(9.11)$$

微分Snell定律，方程(9.2)，关于$\theta$给出有用的关系式

$${\eta }_{o}cos{\theta }_{o}d{\theta }_o={\eta }_{i}cos{\theta }_{i}d{\theta }_i$$

或

$$\frac{cos{\theta }_{o}d{\theta }_o}{cos{\theta }_{i}d{\theta }_i}=\frac{{\eta }_i}{{\eta }_o}$$

将Snell定律和这个关系式同时代入式(9.11)，然后化简，得到

$$L_o{\eta }_i^{2}d{\varphi }_o=L_i{\eta }_o^{2}d{\varphi }_i$$

最后，$d{\varphi }_i=d{\varphi }_o$，给出了入射和透射辐射之间的最终关系:

$$L_o=L_i\frac{{\eta }_o^2}{{\eta }_i^2}，（9.12）$$

因此，一个BTDF所满足的对称关系为

$${\eta }_o^2{f}_t(p,{\omega }_o,{\omega }_i)={\eta }_i^2{f}_t(p,{\omega }_i,{\omega }_o)，(9.13)$$

非对称散射对于从相机和灯光开始采样光路的双向光传输算法来说尤其成问题。如果不考虑非对称性，那么这种算法可能会产生错误的结果，因为这种算法的设计基本上是基于对称性原则的。

我们会说，从光源开始采样的光路携带importance，而从相机开始采样的光路携带radiance。这些术语对应于在路径起点处记录的量。对于importance传输，除非采取特别预防措施，否则bsdf的入射和出射方向参数将(不正确地)反转。

因此，我们定义了伴随BSDF $f^{\ast }$，它的唯一作用是用交换的参数对原始BSDF求值:

$$f^{\ast }(p,{\omega }_o,{\omega }_i)=f(p,{\omega }_i,{\omega }_o)$$

所有基于重要性传输的采样步骤都使用BSDF的伴随形式，而不是其原始版本。pbrt中的大多数bsdf都是对称的，因此 $f$ 和 $f^{\ast }$ 之间没有实际的区别。然而，非对称的情况需要额外的注意。

TransportMode枚举用于通知这些非对称bsdf关于传输数量的信息，以便它们可以正确地在伴随形式和非伴随形式之间切换。

<<TransportMode Definition>>= 
enum class TransportMode { Radiance, Importance };

伴随的BTDF是

$$f_{\mathrm{t}}^{\ast }\left(\mathrm{p},{\omega }_{\mathrm{o}},{\omega }_{\mathrm{i}}\right)=f_{\mathrm{t}}\left(\mathrm{p},{\omega }_{\mathrm{i}},{\omega }_{\mathrm{o}}\right)=\frac{{\eta }_{\mathrm{i}}^2}{{\eta }_{\mathrm{o}}^2}{f}_{\mathrm{t}}\left(\mathrm{p},{\omega }_{\mathrm{o}},{\omega }_{\mathrm{i}}\right),$$

有效地抵消了式(9.12)中的比例因子。

有了这些方程，我们现在可以在DielectricBxDF评估和采样方法的实现中定义剩下的缺失部分。每当辐射在折射边界上传输时，我们应用公式(9.12)中的比例因子。对于重要性传输，我们使用伴随的BTDF，由于公式(9.12)和式(9.13)的组合，它缺少比例因子。

<<Account for non-symmetry with transmission to different medium>>= 
if (mode == TransportMode::Radiance)
    ft /= Sqr(etap);

---

一般情况下，我们都是从相机开始向光源采样，因此mode = TransportMode::Radiance，因为实际光线传播是从光源指向相机，这与我们的追踪方式（从相机指向光源）方向相反，因此我们需要使用$f^{\ast }$。