概述

如果我們有一個(gè)上半球的余弦分布函數(shù)$(D_o)$，并對(duì)這個(gè)余弦分布函數(shù)進(jìn)行一個(gè)矩陣$(M)$變換，變換為一個(gè)新的分布函數(shù)$(D)$。

$D=M{D}_o$

這就是基本思想：用一個(gè)余弦分布函數(shù)來擬合出不同的分布函數(shù)，并且誤差很小。

LTC基本定義與性質(zhì)

我們定義$D_o$為我們的基函數(shù)，也就是我們的原始分布函數(shù)。
我們?cè)?span id="ozvdkddzhkzd" class="katex--inline">$D_o$上采樣一個(gè)方向向量${\omega }_o$，那么我們對(duì)這個(gè)方向向量做一個(gè)線性變換，也就是做一個(gè)$M$矩陣乘法變換為一個(gè)新的方向向量$\omega$。

$\omega =\frac{M{\omega }_o}{||M{\omega }_o||}$,

${\omega }_o=\frac{M^{?1}\omega }{||M^{?1}\omega ||}$,

那么一個(gè)目標(biāo)分布函數(shù)$D(\omega )$如何變換為$D_o({\omega }_o)$呢？

即就是一個(gè)變量替換：

$D(\omega )=D_o({\omega }_o)\frac{?{\omega }_o}{?\omega }=D_o(\frac{M^{?1}\omega }{||M^{?1}\omega ||})\frac{|M^{?1}|}{||M^{?1}\omega |{|}^3}$

$\frac{?{\omega }_o}{?\omega }=\frac{|M^{?1}|}{||M^{?1}\omega |{|}^3}$是Jacobian項(xiàng)，可以從幾何關(guān)系上進(jìn)行推導(dǎo)。

那就讓我們簡(jiǎn)單推導(dǎo)一下吧！(假如想要了解詳細(xì)，請(qǐng)參考原文哦！)

這要從立體角開始推導(dǎo)啦。我們定義一個(gè)原始分布上的立體角$?{\omega }_o$，將其進(jìn)行$M$變換，得到一個(gè)新的立體角$?\omega$。

$?\omega =?{\omega }_{o}A\frac{cos\theta }{r^2}$

這就是一個(gè)立體角的定義。讓我們來探討一下里面的面積與距離的表達(dá)吧！

偷懶一下，可以從圖中看到面積與距離的表達(dá)，那么讓我們來化簡(jiǎn)一下吧。

$M{\omega }_1\times M{\omega }_2=|M|M^{?T}({\omega }_1\times {\omega }_2)=|M|M^{?T}{\omega }_o$

$\frac{?\omega }{?{\omega }_o}=\frac{|M|}{||M{\omega }_o|{|}^3}$

$\frac{?{\omega }_o}{?\omega }=\frac{M^{?1}}{||M^{?1}\omega |{|}^3}$

如果$M$是縮放矩陣或者旋轉(zhuǎn)矩陣，這個(gè)Jacobian項(xiàng)等于一。即分布形狀是不變的。

讓我們假設(shè)矩陣是縮放矩陣($M=\lambda I$)，

$|M^{?1}=|M|=|\lambda I|={\lambda }^3$

$||M^{?1}\omega ||=||\lambda I\omega ||=\lambda$

$\frac{?{\omega }_o}{?\omega }=\frac{|M^{?1}|}{||M^{?1}\omega |{|}^3}=1$

讓我們假設(shè)矩陣是旋轉(zhuǎn)矩陣$M=R$。

$|M|=|M^{?1}|=1$

$||M^{?1}\omega ||=1$

$\frac{?{\omega }_o}{?\omega }=1$

那么被變換后的分布$D$有哪些性質(zhì)呢？

首先是范數(shù)：

${\int }_{\Omega }D(\omega )\mathrm{d}\omega ={\int }_{\Omega }{D}_o\frac{?{\omega }_o}{\omega }\mathrm{d}\omega ={\int }_{\Omega }{D}_o({\omega }_o)\mathrm{d}{\omega }_o$

其次對(duì)于多邊形積分相等：

${\int }_{P}D(\omega )\mathrm{d}\omega ={\int }_{P_o}{D}_o({\omega }_o)\mathrm{d}{\omega }_o$
$P_o=M^{?1}P$

使用LTC近似BRDF

我們?cè)趯?duì)LTC有了一些基本的認(rèn)識(shí)之后，讓我們用它來進(jìn)行擬合BRDF吧！ 這里作者使用了各項(xiàng)同性的GGX微表面分布模型。

對(duì)于一個(gè)原分布$D_o$，我們選擇一個(gè)上半球的余弦分布，當(dāng)然還可以選擇$cosin{e}^2$，均勻半球分布等：

$D_o({\omega }_o=(x,y,z))=\frac{1}{\pi }max(0,z)$

我們來定義一下目標(biāo)分布吧，也就是我們將要擬合的BRDF分布函數(shù)帶著余弦項(xiàng)的：

$D\approx \rho ({\omega }_v,{\omega }_l)cos{\theta }_l$

對(duì)于一個(gè)各項(xiàng)同性的材質(zhì)而言，BRDF函數(shù)分布只依賴于入射方向${\omega }_v=(sin{\theta }_v,0,cos{\theta }_v)$和粗糙度$\alpha$。

我們發(fā)現(xiàn)可以用一個(gè)矩陣$M$來很好的擬合$D$，矩陣被寫為：

$M=?\begin{array}{ccc}a& 0& b\\ 0& c& 0\\ d& 0& 1\end{array}?$

可以看見這個(gè)矩陣只有四個(gè)參數(shù)需要被確定。

可以從圖中看出擬合效果很不錯(cuò)！

使用LTC進(jìn)行多邊形光源的著色

我們?cè)诓煌亩噙呅喂庠?span id="ozvdkddzhkzd" class="katex--inline">$P$上計(jì)算直接光照的著色，${\omega }_v$代表視口方向(view direction)，$\rho$代表BRDF函數(shù)，$L$是多邊形光源出射的radiance。

$I={\int }_{P}L({\omega }_l)\rho ({\omega }_v,{\omega }_l)cos{\theta }_l\mathrm{d}{\omega }_l$

我們使用LTC對(duì)BRDF函數(shù)進(jìn)行擬合，可得：

$I={\int }_{P}L({\omega }_l)\rho ({\omega }_v,{\omega }_l)cos{\theta }_l\mathrm{d}{\omega }_l\approx {\int }_{P}L({\omega }_l)D({\omega }_l)\mathrm{d}{\omega }_l$

這里我們假設(shè)多邊形光源出射的radiance是一個(gè)常數(shù)，即$L({\omega }_l)=L$，那么我們可以進(jìn)一步簡(jiǎn)化：

$I={\int }_{P}L({\omega }_l)D({\omega }_l)\mathrm{d}{\omega }_l=L{\int }_{P}D({\omega }_l)\mathrm{d}{\omega }_l$

我們?cè)龠M(jìn)行積分變換：

${\int }_{P}D(\omega )\mathrm{d}\omega ={\int }_{P_o}{D}_o({\omega }_o)\mathrm{d}{\omega }_o=E(P_o)$

對(duì)于一個(gè)余弦分布的半球面函數(shù)上的多邊形積分是有表達(dá)式的，即：

$E(p_1,...,p_n)=\frac{1}{2\pi }{\sum }_{i=1}^{n}acos(<p_i,p_j>)?\frac{p_i\times p_j}{||p_i\times p_j||},?\begin{array}{c}0\\ 0\\ 1\end{array}??$

這里的$j=(i+1)mod(n)$，注意將多邊形光源進(jìn)行變換后，要將光源進(jìn)行上半球面的裁剪。

這里的證明就不再贅述了。

對(duì)于帶著紋理的光源來說，作者給出的方法很奇怪，就不再討論。

總結(jié)

對(duì)于這篇paper，將一種數(shù)學(xué)方法引入了圖形學(xué)，是很不錯(cuò)的。

優(yōu)點(diǎn)：

• 對(duì)于簡(jiǎn)單多邊形光源(矩形，三角形等)來說，對(duì)于場(chǎng)景的著色很快，效果很棒，誤差小。

缺點(diǎn)：

• 復(fù)雜多邊形光源(心形等)的著色速度慢，已經(jīng)很難達(dá)到實(shí)時(shí)。
• 只能針對(duì)二維平面上的多邊形光源。
• 只能進(jìn)行直接光照著色。

實(shí)現(xiàn)流程：

通過視口方向${\omega }_v，\alpha$進(jìn)行采樣LUT得到矩陣$M$。將面積光源進(jìn)行旋轉(zhuǎn)，即面積光源坐標(biāo)系轉(zhuǎn)變?yōu)橐韵袼攸c(diǎn)為原點(diǎn)的局部坐標(biāo)系。即構(gòu)造一個(gè)TBN矩陣。
mul(Minv, transpose(mat3(T, Ｂ, N)));

將世界坐標(biāo)系下的面積光源的頂點(diǎn)變換。

裁剪多邊形。

進(jìn)行線積分。

這里給出LTC的網(wǎng)站，想了解更詳細(xì)的內(nèi)容可以看看哦！

$對(duì)那充滿不安的日子說再見吧，$
$向著未知的夢(mèng)境之地前行，$
$我想總有一日我會(huì)變得更加堅(jiān)強(qiáng)！$

總結(jié)

以上是生活随笔為你收集整理的Polygonal-Light Shading with LTC的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。