PBR反射率公式推导过程

反射率公式

$L_o(p,{\omega }_o)=\underset{\Omega }{\int }{f}_r(p,{\omega }_i,{\omega }_o)L_i(p,{\omega }_i)n\cdot {\omega }_{i}d{\omega }_i$

辐射率L解析

$L_i(p,{\omega }_i)$在指定方向上的单位立体角和垂直此方向的单位面积上的辐射通量。
能量
$Q=\frac{hc}{\lambda }$
一个波长为$\lambda$的光子携带的能量。
$h$：普朗克常数， $c$：光速
辐射通量
在单位时间通过表面的能量
$\Phi =\frac{\Delta Q}{\Delta t}$
辐照度
单位面积到达表面的通量
$E=\frac{\Delta \Phi }{\Delta A}$
A 面积
任意球面的极小面积
$\Delta A=(rsin\theta )(r\Delta \theta )=r^2(sin\theta \Delta \theta \Delta \varphi )$
极小立体角
$\Delta \omega =\frac{\Delta A}{r^2}=sin\theta \Delta \theta \Delta \varphi$
由前面的计算推出辐射率公式
$L(p,\omega )=li{m}_{\Delta \omega \to 0}\frac{\Delta E_{\omega }(p)}{\Delta \omega }=\frac{d\frac{d\varphi }{dA}}{d\omega }=\frac{d^2\Phi }{d\omega dA}$
这里的$dA$应该与立体角$\omega$垂直，所以应该乘以$cos\theta$。所以最终的辐射率公式：
$L(p,\omega )=\frac{d^2\Phi }{d\omega dAcos\theta }$

BRDF 双向反射分布函数

反射率公式中的$f_r$中是BRDF，反映了光线对平面最终反射出的光线所作出的贡献。
${\omega }_i$入射光，${\omega }_o$观察方向， $n$法线， $\alpha$粗糙度
$f_r=k_d{f}_{lambert}+k_s{f}_{cook-torrance}$
$k_d$是入射光线中被折射部分的能量所占的比率， $k_s$是被反射部分的比率。
BRDF的左侧表示漫反射部分
$f_{lambert}=\frac{c}{\pi }$
$c$表示颜色， 除以$\pi$是为了对漫反射光进行标准化。
** 镜面反射部分**
BRDF的右侧是镜面反射部分
$f_{cook-torrance}=\frac{DFG}{4({\omega }_o\cdot n)({\omega }_i\cdot n)}$
法线分布函数（D)(Normal Distribution Function）：估算在受到表面粗糙度的影响下，取向方向与中间向量一致的微平面的数量。这是用来估算微平面的主要函数。

几何函数（G）(Geometry Function)：描述了微平面自成阴影的属性。当一个平面相对比较粗糙的时候，平面表面上的微平面有可能挡住其他的微平面从而减少表面所反射的光线。

菲涅尔方程（F）(Fresnel Rquation)：菲涅尔方程描述的是在不同的表面角下表面所反射的光线所占的比率。

法线分布函数

$ND{F}_{GGXTR}(n,h,\alpha )=\frac{{\alpha }^2}{\pi ((n\cdot h{)}^2({\alpha }^2-1)+1{)}^2}$
$h$中间向量， $\alpha$粗糙度

几何函数
$G_{SchlickGGX}(n,v,k)=\frac{n\cdot v}{(n\cdot v)(1-k)+k}$
$k$是$\alpha$基于几何函数是针对直接光照还是针对IBL光照的重映射
$k_{direct}=\frac{(\alpha +1{)}^2}{8}$
$k_{IBL}=\frac{{\alpha }^2}{2}$
菲涅尔方程
$F_{Schlick}(h,v,F_0)=F_0+(1-F_0)(1-(h\cdot v){)}^5$
$F_0$表示平面的基础反射率， 利用折射指数计算所得。

反射率方程最终公式

$L_o(p,{\omega }_o)=\underset{\Omega }{\int }(k_d\frac{c}{\pi }+k_s\frac{DFG}{4({\omega }_o\cdot n)({\omega }_i\cdot n)})L_i(p,{\omega }_i)n\cdot {\omega }_{i}d{\omega }_i$