球諧光照在高度場渲染中的應(yīng)用

趙克勝+倪桂強(qiáng)+羅健欣

DOI：10.16644/j.cnki.cn33-1094/tp.2016.02.010

摘 ?要： 為提高高度場渲染的真實(shí)感，增加視覺觀察的準(zhǔn)確性，提出了使用球諧光照渲染高度場的方法。使用OpenGL圖形接口進(jìn)行實(shí)驗(yàn)程序開發(fā)，采用Puget Sound高度圖作為原始數(shù)據(jù)，利用對比實(shí)驗(yàn)說明球諧光照在高度場渲染中對渲染效果、渲染效率的提升。對比結(jié)果表明，將球諧光照用于高度場渲染可以有效增強(qiáng)高度場真實(shí)感。

關(guān)鍵詞： 高度場渲染; 球諧光照; 蒙特卡洛積分; 關(guān)聯(lián)勒讓德多項(xiàng)式

中圖分類號：TP391.41 ? ? ? ? ?文獻(xiàn)標(biāo)志碼：A ? ? 文章編號：1006-8228（2016）02-33-05

Terrains rendering with spherical harmonic lighting

Zhao Kesheng， Ni Guiqiang， Luo Jianxin

（Laboratory of Military Network Technology， PLA University of Science and Technology， Nanjing， Jiangsu 21007， China）

Abstract： In order to improve the sense of reality of terrains rendering and the accuracy of visual observation， a method of rendering terrains with spherical harmonic lighting is proposed. In the experiment， the OpenGL graphics interface is used as the development environment and the Puget Sound terrains as the raw data. The advantage of terrains rendering with spherical harmonic lighting in rendering quality and efficiency is described by the comparative experiment. The result shows that spherical harmonic lighting can improve the rendering quality of terrains rendering.

Key words： terrains rendering; spherical harmonic lighting; Monte-Carlo integration; associated Legendre polynomials

0 引言

近年來球諧（Spherical Harmonic）在光照、BRDF、形狀識別等計(jì)算機(jī)圖形學(xué)領(lǐng)域內(nèi)被廣泛使用[1]。球諧光照[2-5]是利用球諧計(jì)算3D模型上光照的技術(shù)，它可以實(shí)時(shí)地渲染出全局光照模式的圖像[6]。相對于光線追蹤[7]、輻射度算法[8]等光照技術(shù)，球諧光照具有渲染效果好、實(shí)時(shí)性強(qiáng)、代碼易編寫等優(yōu)點(diǎn)。

高度場渲染在模擬飛行、虛擬戰(zhàn)場三維游戲等領(lǐng)域都有著廣泛應(yīng)用[9]。通過對高度場場景模型加入光照效果將會(huì)大大增強(qiáng)結(jié)果圖像的真實(shí)感[10]，為交互可視化和視覺觀察提供很大的便利。由于高度場本質(zhì)上也是3D模型，所以本文研究了將球諧光照應(yīng)用于高度場渲染，并對其優(yōu)缺點(diǎn)進(jìn)行了分析。

1 球諧光照

1.1 光照計(jì)算

圖1所示為自然界中的某個(gè)光照場景，光線自左側(cè)傳播至x處，經(jīng)物體反射至x'處。此時(shí)觀察者于x'觀察到的x點(diǎn)的顏色由兩部分組成：一部分為x點(diǎn)本身發(fā)射出的沿xx'方向的光的顏色，另一部分為入射光線于x處經(jīng)物體反射至x'處的光的顏色。在計(jì)算機(jī)圖形學(xué)中，模擬這種物理現(xiàn)象的公式稱為渲染方程[11]：

式⑴

其中是觀察者在方向看到的x點(diǎn)的顏色。是物體在x點(diǎn)沿著方向自身發(fā)射的光的顏色。將沿的入射光顏色轉(zhuǎn)化為沿方向的反射光顏色。L（x'，）是從x'點(diǎn)發(fā)射沿方向的入射光，G（x，x'）是x與x'的幾何關(guān)系。V（x，x'）是x點(diǎn)對x'點(diǎn)的可見測試。

<E：＼方正創(chuàng)藝5.1＼Fit201601＼圖＼zks圖1.tif>

圖1 ?自然界中的光照場景

可以看出，使用圖形硬件來計(jì)算的難點(diǎn)在于，渲染公式中存在一個(gè)球上的積分，計(jì)算基于符號的積分在GPU中很難完成。為了能快速地得到結(jié)果，人們使用了一個(gè)近似的解法，這個(gè)方法就是球諧投影與重建。球諧投影與重建的原理是，使用球諧函數(shù)求原函數(shù)的近似解。球諧函數(shù)是一類擁有特殊性質(zhì)的基函數(shù)。

1.2 基函數(shù)

基函數(shù)是用來求解原函數(shù)近似解的一系列函數(shù)。利用基函數(shù)求原函數(shù)近似解的過程如下：

首先在原函數(shù)的定義域上積分原函數(shù)f（x）i與基函數(shù)Bi（x）的乘積，得到一系列的近似系數(shù)Ci，如圖2。

然后用系數(shù)Ci乘以各自對應(yīng)的基函數(shù)Bi（x），如圖3。

最后將乘積結(jié)果累加，得到原函數(shù)的近似結(jié)果，如圖4。

1.3 關(guān)聯(lián)勒讓德多項(xiàng)式

球諧函數(shù)的核心部分是關(guān)聯(lián)勒讓德多項(xiàng)式，關(guān)聯(lián)勒讓德多項(xiàng)式是標(biāo)準(zhǔn)正交多項(xiàng)式的一族。標(biāo)準(zhǔn)正交多項(xiàng)式是擁有以下的性質(zhì)的一組基函數(shù)：積分任何兩個(gè)正交基函數(shù)的積時(shí)，如果它們相同，結(jié)果為1，如果它們不同，結(jié)果為0即：

式⑵

關(guān)聯(lián)勒讓德多項(xiàng)式通常使用符號P表示，關(guān)聯(lián)勒讓德多項(xiàng)式有兩個(gè)參數(shù)：l和m。l的定義域?yàn)閺?開始的整數(shù)，m的定義域?yàn)?到l的整數(shù)。系數(shù)l和m把關(guān)聯(lián)勒讓德多項(xiàng)式分成了不同的階，l是階的索引，同一階內(nèi)的多項(xiàng)式正交于同一個(gè)常數(shù)，不同階內(nèi)的多項(xiàng)式正交于不同的常數(shù)。這種關(guān)系可以用如圖5中的三角網(wǎng)格來表示：

圖5 ?三角網(wǎng)格

我們使用遞歸定義來求解給定參數(shù)的關(guān)聯(lián)勒讓德多項(xiàng)式，這個(gè)過程需要用到下面三個(gè)公式。

式⑶

首先使用這個(gè)式子和給定的m求解出，注意x！！為雙階乘運(yùn)算，其結(jié)果為所有小于等于x并與x有相同奇偶性的整數(shù)的乘積，例如5！！=5*3*1=15。如果此時(shí)l=m，就終止程序，如果l≠m則使用公式⑶來求：

式⑷

上一步由公式⑵求得的代入此式，即可求得。此式只使用一次，若此時(shí)m+1=l，則終止程序，如果m+1≠1則使用公式4遞歸求解：

式⑸

1.4 球諧函數(shù)

在球面上，我們使用球諧函數(shù)來求原函數(shù)的近似解。與關(guān)聯(lián)勒讓德多項(xiàng)式相似，球諧函數(shù)也擁有兩個(gè)參數(shù)l和m。l將球諧函數(shù)分為了不同的階，它的取值范圍從0開始的整數(shù)，與關(guān)聯(lián)勒讓德多項(xiàng)式不同的是球諧函數(shù)的參數(shù)m的取值范圍為-l到l的整數(shù)。球諧函數(shù)通常用符號y來表示，它的定義如下：

式⑹

P為關(guān)聯(lián)勒讓德多項(xiàng)式，K是比例因子：

式⑺

定義順序yi如下：

式⑻

1.5 球諧投影

觀察渲染方程發(fā)現(xiàn)，G（x，x'）和V（x，x'）均為常數(shù)，則渲染方程中積分的實(shí)際形式為，它的求解過程如下。使用球諧函數(shù)來近似L（s）和t（s），方法與使用基函數(shù)求原函數(shù)近似的過程相同。

用和來代替L（s）和t（s）：

展開：

繼續(xù)展開：

由于球諧函數(shù)是正交多項(xiàng)式的子族，所以擁有以下性質(zhì)：

故有：

其中：

求解這兩個(gè)定積分的方法是蒙特卡洛積分。

1.6 蒙特卡洛積分

在實(shí)際中，許多需要計(jì)算多重積分的復(fù)雜問題，使用蒙特卡洛積分[12-14]都能有效地解決。它的推導(dǎo)過程如下。

已知函數(shù)f（x）的期望等于函數(shù)本身f（x）與它的概率密度函數(shù)p（x）的乘積在函數(shù)f（x）定義域上的積：

另一個(gè)計(jì)算函數(shù)期望的方法是，取大量的函數(shù)隨機(jī)采樣的平均值，若采樣數(shù)趨于無窮大，則最終結(jié)果收斂于函數(shù)期望：

所以，由此可推得：

將上式兩側(cè)都除以p（x），得到蒙特卡洛積分：

式⑼

通過這種轉(zhuǎn)換，一個(gè)函數(shù)的積分運(yùn)算就轉(zhuǎn)化成了四則運(yùn)算，在計(jì)算機(jī)中可以方便地求得結(jié)果。

運(yùn)用蒙特卡洛積分方法求解1.5節(jié)中的積分得：

由于渲染方程需要解決的是一個(gè)球面上的積分，所以利用蒙特卡洛積分求解渲染方程時(shí)需要在整個(gè)球面上隨機(jī)采樣離散分布的點(diǎn)。可以采用某種策略先得到2D的采樣點(diǎn)，然后將其映射到球面空間上。映射的公式如下：

式⑽

如圖6所示，類似于將平面的世界地圖貼到地球儀上，為了保證映射后的采樣點(diǎn)隨機(jī)且均勻，需要注意到采樣點(diǎn)在2D坐標(biāo)下的分布規(guī)律為沿兩極向赤道越來越密集。具體的采樣方法是創(chuàng)建一個(gè)1*1的采樣空間，劃分成N*N的網(wǎng)格，每個(gè)網(wǎng)格內(nèi)取一個(gè)坐標(biāo)，用公式⑽將其轉(zhuǎn)化為球坐標(biāo)，再利用公式⑸求得每個(gè)采樣點(diǎn)對應(yīng)的球諧函數(shù)值。

因?yàn)槭窃趩挝磺蛏暇鶆蚍植嫉牟蓸狱c(diǎn)，所以。因此有：

這樣利用球諧投影與重建求解3D模型上某一點(diǎn)經(jīng)過光照后顏色的整個(gè)過程就結(jié)束了。

2 利用球諧渲染高度場

2.1 實(shí)驗(yàn)準(zhǔn)備

實(shí)驗(yàn)硬件平臺為Intel（R） Core（TM） i5-4460 CPU @ 3.20GHz、8GB內(nèi)存，Radeon X1300/X1550 Seris顯卡的PC，采用VC++作為開發(fā)語言，VS2012作為開發(fā)平臺。實(shí)驗(yàn)數(shù)據(jù)來源于Puget Sound數(shù)據(jù)集，為1025*1025分辨率的原始高程數(shù)據(jù)。

2.2 同一四叉樹映射

高度場渲染常用四叉樹作為數(shù)據(jù)管理方式。現(xiàn)有的四叉樹結(jié)構(gòu)渲染節(jié)點(diǎn)的過程需要對節(jié)點(diǎn)關(guān)聯(lián)的網(wǎng)格進(jìn)行簡化并將簡化后的網(wǎng)格發(fā)送到GPU上。由于四叉樹節(jié)點(diǎn)關(guān)聯(lián)的網(wǎng)格三角形數(shù)目巨大，所以這一過程往往會(huì)限制渲染效率。

通過采用同一四叉樹映射[15]的方法可以有效提高渲染效率。同一四叉樹映射基于GPU支持頂點(diǎn)可編程性、強(qiáng)大的處理能力等特點(diǎn)。將統(tǒng)一的均勻網(wǎng)格駐留在GPU中，在渲染四叉樹節(jié)點(diǎn)時(shí)，通過簡單的頂點(diǎn)著色將網(wǎng)格頂點(diǎn)移動(dòng)到指定位置進(jìn)行渲染。這樣做雖然增加了三角形數(shù)量，但是因?yàn)橥耆珌G棄了網(wǎng)格簡化及將網(wǎng)格數(shù)據(jù)發(fā)送到GPU中的過程，所以反而提高了渲染效率。

2.3 渲染結(jié)果

圖7為無陰影的球諧光照渲染的高度場圖片，圖8為相同渲染場景的對比實(shí)驗(yàn)效果。對比畫框處，球諧光照的陰影效果更好，更易于分辨地形的起伏變化。

圖9為球諧光照渲染的高度場的平視圖，圖10為相同場景的對比實(shí)驗(yàn)效果。對比畫框處，球諧光照背光處的立體感很真實(shí)，視覺觀察的準(zhǔn)確性更高。

圖11為球諧光照渲染的草地場景，圖12為相同渲染場景的對比實(shí)驗(yàn)效果。可以看出，球諧光照把草地陰影的效果渲染得很好。

3 結(jié)束語

實(shí)驗(yàn)結(jié)果表明，將球諧光照用于高度場的渲染能有效提升高度場圖片的真實(shí)感，為使用者的視覺觀察與可視化交互增加準(zhǔn)確性。球諧光照用于高度場的缺點(diǎn)在于，球諧光照的預(yù)處理時(shí)間隨著網(wǎng)格節(jié)點(diǎn)的增長呈O（n2）的增長趨勢，所以，如何有效簡化網(wǎng)格成為下一步研究的重點(diǎn)。

參考文獻(xiàn)（References）：

[1] Lisle， Ian & Huang， Tracy（2007）.Algorithms for spherical

harmonic lighting[R]. GRAPHITE 2007： Proceedings of the 5th International Conference on Computer Graphics and Interactive Techniques in Australia and Southeast Asia，2007.12：1-4

[2] P.P.Sloan，J.Kautz，and J.Snyder，"Precomputed radiance

transfer for real-time rendering in dynamic，low-frequency lighting environments"[J].in SIGGRAPHPH' 02：Proceedings of the 29th annual conference on computer graphics and interactive techniques.New York，NY，USA：ACM Press，2002：527-536

[3] X.Liu， P.Sloan， H-Y.Shum， J.Synder. All-frequency

precomputed radiance transfer for glossy objects[J].Proc of Eurograhis Symposium on Rendering'04，2004：337-344

[4] P.Sloan， J.Hall， J.Hart， J.Synder. Clustered principal

components for precomputed radiance transfer[J].Proc.OF SIGGRAPH'03，2003：382-391

[5] Sloan P， Liu X， Shum H， et al. Bi-scale radiance transfer[R].

ACM-Transactions onGraphics， 2003.22（3）：370-375

[6] Robin Green，Spherical Harmonic Lighting： The Gritty

Details.

[7] Whitted J.T.，An Improved Illumination model for shaded

display[J]，CACM，1980.23：343-349

[8] Shirley P， A Ray Tracing Method For Illumination

Calculation In Diffuse-specular Scense[J]. In：Proc of CG International'90， New York，1990：469-476

[9] 劉思江，翁璐斌，秦樹鑫，楊一平，蔣永實(shí).保持邊界到的地形

網(wǎng)格簡化算法與組織格式[R].計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào)，2013.25（1）.

[10] Luo JX，Ni G Q，Cui P，et al. Quad-tree atlas ray casting：a

gpu based frameword for terrain visualization and its applications [M].Lectrue Notes in Computer Science. Heidelberg：Springer，2012，7145：74-85

[11] Kajiya，J.The Rending Equation[J]. SIGGRAPH' 86，

143-150

[12] L.Szirmay-Kalos， Monte-carlo global Illumination

Methods-State of the Art and New Developments[R].Budimerce， Slovakia：Comenius University，1999：3-21

[13] Wald I， Realtime Ray Tracing and Interactive Global

Illmination[D]. Phd thesis， Saarland University，2004.

[14] L.Szormay-Kalos. Monte-carlo global illumination

methods-state of the art and new developments[R]. In proceedings of the Fifteenth Spring Confernce on Computer Graphics. Budimerce， Slovakia：Comenius University，1999：3-21

[15] 羅健欣，胡谷雨，倪桂強(qiáng).平行流形空間光線投射高度場可

視化算法[J].計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào)，2013.25（3）.