• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    多光源繪制方法綜述

    2020-01-09 03:39:46劉逸凡
    計算機研究與發(fā)展 2020年1期
    關(guān)鍵詞:著色光源繪制

    劉逸凡 徐 昆

    (清華大學(xué)計算機科學(xué)與技術(shù)系 北京 100084)

    全局光照效果(global illumination)包含直接光照效果(direct illumination)和間接光照效果(indirect illumination),可以表現(xiàn)反射、折射、焦散(caustics)、互反射等效果,繪制結(jié)果具有很高的真實感.近年來,隨著游戲、影視、動畫等圖形學(xué)相關(guān)產(chǎn)業(yè)的發(fā)展,全局光照效果成為普遍而重要的需求[1-2].

    計算某點的全局光照,既要考慮光源對該點的直接光照,又要考慮來自光源的光線與其他物體表面發(fā)生1次或多次反射折射后對該點的間接光照.在光線跟蹤框架[3]中,可以通過從視點處反向追蹤光路來計算到達成像平面的全局光照.如圖1所示,對于成像平面上的每一個像素,從視點發(fā)射若干條射線.對于每條射線,嘗試與場景中的物體求交:如相交,則可能在交點處發(fā)生反射或折射,產(chǎn)生新的反射折射射線,然后繼續(xù)與場景中的其他物體求交,如此遞歸下去.在圖1中,點x1的出射光亮度L1,由直接光照和間接光照2部分組成:通過依次計算場景中每個光源對點x1的貢獻并求和,得到直接光照;通過計算點x2的出射光亮度L2對點x1的貢獻,得到間接光照.

    Fig. 1 Global illumination in path tracing圖1 光線跟蹤中的全局光照

    一方面,通過引入虛擬點光源(virtual point lights, VPLs),全局光照的計算(包括間接光照分量)可以轉(zhuǎn)化為多光源下的直接光照計算;另一方面,隨著工業(yè)界對于圖像真實感訴求的提高,越來越多的繪制場景本身就含有大量光源:夜景的霓虹燈、火星四濺的粒子特效、復(fù)雜網(wǎng)格光源的模擬等.這使得多光源繪制成為計算機圖形學(xué)和繪制領(lǐng)域的一個重要研究內(nèi)容,涌現(xiàn)出了一大批重要方法.

    理論上,多光源繪制計算的時間復(fù)雜度與場景中光源的數(shù)量成正比.如此一來,隨著光源數(shù)量的增多,每個點的光照計算量也會成線性增加,造成總體繪制時間的線性增長.例如,圖2是動畫電影《尋夢環(huán)游記》中的一個場景,它包含了約800萬個光源,如果不采取更好的優(yōu)化算法,如此多光源所帶來的繪制時間的增加是難以處理的.

    Fig. 2 Complex illumination in Coco(? 2017 Disney·Pixar)[2]圖2 《尋夢環(huán)游記》中的復(fù)雜光照[2]

    為了提高多光源場景的繪制效率,研究人員一般從3個方面入手:1)使用蒙特卡洛方法對直接光照進行重要性采樣,以縮短繪制時間;2)從占據(jù)計算量主體的可見性測試入手,提高單個光源的光照計算效率;3)通過光源聚類等方法實現(xiàn)繪制時間對光源數(shù)量的可伸縮性(scalable),使總體繪制的時間復(fù)雜度由O(n)變?yōu)镺(logn),其中n為光源數(shù)量.

    根據(jù)單幀圖像的繪制時長,繪制可以分為實時繪制和離線繪制.實時繪制與離線繪制的實現(xiàn)方法有著很大差異,對于多光源問題的處理方式也有所不同.這里,本文著重分析離線繪制的優(yōu)化算法,簡要介紹實時繪制的相關(guān)方法.

    本文將首先闡述直接光照計算和重要性采樣的基本概念;然后從加速可見性計算、通過光源聚類加速光照計算、通過光源聚類優(yōu)化重要性采樣等方面,對現(xiàn)有的離線方法進行總結(jié)和分析;最后對實時繪制中的光源剔除方法進行介紹.

    1 直接光照與重要性采樣

    我們將物體表面需要進行繪制的點稱作著色點(shading point).如圖3所示,僅考慮直接光照,著色點x在反射方向ω的反射光亮度,可通過對所有光源表面A進行積分得到[3-4]

    (1)

    其中,積分項F(·)定義為

    Le(y→x)B(y→x→ω)V(y?x)G(y?x),

    (2)

    其中,y為光源表面一點.Le(y→x)為y朝x發(fā)射的光亮度.B(y→x→ω)為雙向反射分布函數(shù)(bidirectional reflectance distribution function, BRDF),描述x處的材質(zhì)反射性質(zhì).V(y?x)是二值可見性函數(shù):x,y之間可見為1,不可見為0.幾何項

    (3)

    其中d(y,x)為x,y之間的歐幾里得距離;θx為x,y連線與x法向之間的夾角;θy為x,y連線與y法向之間的夾角.

    由于式(1)沒有解析解,一般使用蒙特卡洛方法將其近似為N次采樣的積分項的加權(quán)平均:

    (4)

    其中,p(y|x,ω)表示給定著色點x和方向ω,采樣光源上一點y的概率密度函數(shù)(probability density function, PDF).概率密度函數(shù)與積分項越相似,蒙特卡洛估計的方差就越小.當(dāng)概率密度函數(shù)正比于積分項時,方差為0,此時,式(4)可以看作等式.這種采樣方法也稱作重要性采樣(importance sampling).

    Fig. 3 Illustration of direct illumination圖3 直接光照示意圖

    一般地,尤其是對于非環(huán)境光源,采樣光源上一點y,通常分解為2步[4]:先采樣一個光源l,再在光源上采樣一個點y.則式(4)變換為

    (5)

    關(guān)于在給定幾何形狀的光源上進行點的采樣這一問題(即得到p(y|l,x,ω)),許多技術(shù)已經(jīng)提供了近似最優(yōu)的解決方法[4-5],因此不再是研究重點.而如何合理對光源本身進行采樣(即得到p(l|x,ω)),是直接光照計算中一個開放性的研究問題.

    在對光源進行重要性采樣時,采樣的概率分布與光源對場景的貢獻越相近,則生成的圖片噪點越少.為了提高繪制質(zhì)量,往往需要為每個著色點分別計算光源采樣的概率密度函數(shù):估計每一個光源的重要性程度,從而形成離散的概率分布.關(guān)于光源重要性的估計,理想情況是考慮到式(2)中的所有項:發(fā)射光亮度、BRDF、可見性和幾何項.但為了保證計算效率,往往只能考慮其中的一部分.應(yīng)該考慮哪幾項?又如何進行快速估計?這些問題是重要性采樣研究中的主要研究目標(biāo).

    下面,本文將從提高單個光源的光照估計效率和實現(xiàn)光照估計的亞線性復(fù)雜度這2個角度,介紹離線繪制對于多光源場景的解決辦法.

    2 加速可見性計算

    2點之間可見性的判斷占據(jù)整個繪制計算量的主體.許多工作都是從可見性入手,提高繪制效率.

    由于暗光源對場景整體貢獻較小,1994年Ward[6]提出在直接光照計算中對暗光源跳過可見性測試以加速繪制.1996年Shirley等人[5]將該思想拓展到重要性采樣中,只對亮度大的光源進行采樣.由于不考慮暗光源,這些方法造成了繪制結(jié)果的有偏性.

    1995年Jensen等人[7]的方法和2000年Keller等人[8]的方法通過擴展光子圖(photon map)方法來減少繪制時的可見性計算.然而對于大規(guī)模場景,產(chǎn)生的光子常常不處于視點重要區(qū)域,難以取得理想效果.

    1999年Hart等人[9]利用圖像空間的連貫性提前記錄每個“像素-光源”對之間的潛在遮擋物體,后續(xù)計算著色點可見性時只需考慮這些記錄的潛在遮擋物體.2002年Fernandez等人[10]將該思想擴展到場景空間.他們對場景分塊,預(yù)先將區(qū)域塊與光源之間的可見性關(guān)系分為3種進行記錄:完全可見、完全不可見、不確定.在繪制時,只對那些相對著色點所在區(qū)域塊可見性關(guān)系為不確定的光源進行可見性測試.2003年Wald等人[11]將類似方法用到了計算采樣權(quán)重上.

    2009年Dong等人[12]提出了VPLs框架下的可見性加速計算方法.他們用K-Means方法將虛擬點光源聚類為虛擬面光源,采用虛擬面光源計算可見性的同時,仍使用虛擬點光源進行其他項的計算.

    Fig. 4 Algorithm overview of ref [13]圖4 文獻[13]的算法總覽

    相比于緩存1對1的可見性,Ulbrich等人[14]通過使用陰影圖記錄1對多的可見性結(jié)果.他們在物體表面離散地分布一些記錄點,每個記錄點存儲1張陰影圖.估計任意2點間的可見性時,尋找最近的陰影圖,然后將2點間的距離與陰影圖中的相應(yīng)值進行比較.

    這類緩存記錄的方法,通過犧牲一定的存儲,避免了一部分可見性測試的重復(fù)計算,從而提高了計算效率.然而對于復(fù)雜的場景,由于物體表面的相似性不高,緩存的方法往往難以得到很好的效果,并且會占用大量存儲.

    同年,Billen等人[15]提出一種基于隨機選取的可見性估計算法:計算2點之間可見性時,利用蒙特卡洛方法隨機選取一部分場景物體進行測試.理論上,2點之間的可見性值,是分別只考慮單個潛在遮擋物體時可見性值的累乘.該算法的核心思想是將累乘分解為累加,然后進行蒙特卡洛采樣.例如,當(dāng)只存在2個潛在的遮擋物體A和B時,考察點x和點y之間的可見性V(x,y),它等于VA(x,y)·VB(x,y).通過分解得到:

    (6)

    (7)

    假設(shè)采樣概率p1,p2,p3均為13,相比原來需要測試2個物體,現(xiàn)在平均只需要測試 (1+2+3)3=1.33個物體,減少了可見性測試的次數(shù).

    2014年Billen等人[16]進一步提出利用蒙特卡洛估計來降低可見性測試中潛在遮擋物體的平均網(wǎng)格復(fù)雜度.在測試2點之間的可見性時,隨機地選擇遮擋物體的簡化模型或原始網(wǎng)格模型.與傳統(tǒng)方法相比,該方法并不具有太多優(yōu)勢.主要因為,很多情況下簡化網(wǎng)格和原始網(wǎng)格的可見性測試的時間消耗相差不大.

    上面2種算法雖然在效果提升上不是十分顯著,但巧妙地將蒙特卡洛估計運用到可見性計算中,非常具有啟發(fā)意義.

    表1從準(zhǔn)確性、存儲占用、應(yīng)用范圍3個方面對上述加速可見性計算的工作進行了分析總結(jié).本節(jié)介紹的算法,通過加速可見性測試降低了單個光源的計算時間消耗,對于小規(guī)模光源場景可以取得較好的效果.然而,繪制時間與光源數(shù)量成線性的復(fù)雜度并沒有降低,隨著場景中的光源數(shù)量繼續(xù)增多,僅僅加速可見性測試難以顯著提高計算效率.這時,常見的做法是將光源聚類,相似的一組光源用一個單獨的大光源替代,從而使得繪制時間相比光源數(shù)量的復(fù)雜度從線性降為對數(shù).下面,我們將首先介紹光源聚類的基本方法,然后介紹基于光源聚類的重要性采樣方法.

    Table 1 Analysis of Visibility Test Acceleration Methods表1 可見性測試加速方法比較

    3 使用光源聚類加速直接光照計算

    場景中每個光源對繪制結(jié)果的貢獻往往是不均勻的.部分光源有著較低的重要性,例如,當(dāng)光源被遮擋時,或光源距離著色區(qū)域很遠;同時,部分光源是非常重要的,比如對于繪制結(jié)果中的高光有貢獻的光源,我們需要精確計算它們的光能貢獻.通過利用場景中光源重要程度的不均勻性和光源之間的相似性,對光源進行自適應(yīng)的聚類,并以類為單位進行直接光照計算,是光源聚類加速方法的主要思想.

    使用光源聚類來加速直接光照計算的相關(guān)研究主要是基于VPLs算法框架.VPLs算法源自Keller的即時輻射度方法[17],其核心思想是將復(fù)雜的全局光照計算問題近似為大量虛擬點光源(VPLs)下的直接光照計算問題.為了確保圖像的真實性,VPLs算法往往需要生成數(shù)量巨大的虛擬點光源.此時,通過依次遍歷光源來計算直接光照的蠻力方法不再適用,需要有更優(yōu)化的算法來降低線性時間復(fù)雜度,實現(xiàn)可伸縮性.近10多年,出現(xiàn)了許多研究成果,其中2014年的1篇綜述[18]詳細介紹了VPLs算法的原理以及相關(guān)算法,感興趣的讀者可以進行查閱.本文將著重介紹2014年之后的相關(guān)工作.

    需要注意的是,基于VPLs框架的光源聚類不一定能直接泛化到一般情形下的直接光照計算中.原因主要有2點:1)VPLs框架中的虛擬光源往往是點光源,基于此的一些光源聚類方法可能不適用于球光源、面光源等其他光源類型;2)VPLs算法的應(yīng)用場合決定了它往往不需要過多的像素采樣數(shù),基于此的一些光源聚類算法會預(yù)先計算和存儲所有著色點,用于后續(xù)的繪制過程.當(dāng)像素采樣數(shù)增多時,這種算法會占用大量的存儲空間,不適合泛化到一般情境中.

    光源聚類的方法大致可分為2類:基于樹的方法和基于矩陣的方法.不同聚類方法的主要區(qū)別在于,估計光源重要性的標(biāo)準(zhǔn)以及聚類的形式存在差異.但總體來說,光源聚類的核心思路都是將相似的光源聚成一類,形成一個光源簇(light cluster).每個光源簇的總體光能貢獻,可以通過將該光源簇看作一個更大、更亮的光源來進行近似估計,如圖5所示:

    Fig. 5 Lights clustering when evaluating VPLs lighting [18]圖5 估計VPLs光照時的光源聚類[18]

    3.1 基于樹的光源聚類

    Fig. 6 A light tree and three example lightcuts[20]圖6 1個光源樹和3個光源割[20]

    Fig. 7 Product graph [24]圖7 乘積圖 [24]

    Lightcuts算法需要為每個著色點選擇一個合適的光源割.當(dāng)光源數(shù)量很多且分布復(fù)雜時,從光源樹中搜索合適的光源割也需要花費不少的計算時間.針對這一問題,一些工作提出了改進方法.

    2006年Walter等人[24]將Lightcuts算法擴展,提出多維Lightcuts算法(multidimensional lightcuts),以便能夠處理高維空間的繪制積分,例如景深、運動模糊、空間抗鋸齒等.對于每個像素,他們將高維積分離散成一系列著色點的光照求和.由此產(chǎn)生的著色點數(shù)量非常多,對每個著色點都重新計算一個光源割將帶來巨大的時間成本.他們分別將同一個像素內(nèi)的著色點用二叉樹組織起來,并通過使用圖7所示的乘積圖(product graph)來隱式地表達“著色點-光源”對的層次結(jié)構(gòu),從而提高了效率.

    2013年王光偉等人[25]提出對所有的著色點進行空間聚類,同一類中的著色點使用同一個光源割作為初始光源割,然后對同一類中的每個著色點進一步細化光源割.這樣的做法避免了每次都從光源樹的根節(jié)點開始搜索光源割,一定程度上減少了計算量.

    2015年Bus等人[26]提出的IlluminationCut沿用多維Lightcuts[24]的思路,但不再只考慮同一像素內(nèi)的著色點,而是將所有著色點都用一個二叉樹組織起來,并用乘積圖(圖7所示)隱性地對“著色點-光源”對進行聚類.這樣一來,只需遍歷一次乘積圖(等價于同時遍歷一次著色點樹和光源樹),即可得到每個著色點的光源割,平攤了計算時間.

    Fig. 9 The light transport matrices of two complex scenes (subsampled from the original) [31]圖9 2個樣例場景的光能傳輸矩陣(欠采樣之后的結(jié)果)[31]

    2016年Rehfeld等人[27]提出一種Lightcuts插值算法.他們利用光源割的局部連貫性,先選取一部分著色點按照原始的Lightcuts算法計算光源割,然后通過插值來得到其余著色點的光源割.如此一來,減少了光源割的計算次數(shù),提高了計算效率.

    上面介紹的算法都是基于Lightcuts的聚類思路,即先建立光源二叉樹,從光源樹根節(jié)點自頂向下通過誤差上界標(biāo)準(zhǔn)來得到光源割.除此之外,一些工作提出了其他的基于樹的聚類方法.

    2015年Bus等人[28]基于VPLs框架提出了一種不依賴視角信息的點光源聚類算法.它的核心思路是:1)在預(yù)處理過程中,針對每個點光源,將其余點光源聚類;2)繪制過程中,對著色點計算直接光照時,使用最近鄰的點光源的聚類作為該著色點的初始聚類,然后進行細微調(diào)整.這樣一來,由于預(yù)處理過程中的聚類操作只依賴場景物體信息和光源信息,視點的變換不需要重新進行聚類操作,節(jié)省了大量的計算時間.在具體的聚類方法上,他們以文獻[29]作為理論依據(jù),使得聚類后的光源簇與對應(yīng)點滿足完全分離(well-separated)標(biāo)準(zhǔn).如圖8所示,這保證了對應(yīng)點到光源簇內(nèi)任意一點的距離和朝向都是相似的.然而,由于這種聚類標(biāo)準(zhǔn)沒有考慮到著色點的材質(zhì)信息,該算法更適用于漫反射表面.另外需要注意的是,VPLs框架中,虛擬點光源的分布與場景幾何分布非常接近,因此對著色點的光源聚類使用最近鄰點光源對應(yīng)的光源聚類來代替在直觀上是合理的.但這一思路不易推廣到光源分布與場景幾何分布不一致的一般情形.

    Fig. 8 Definition of well-separated[28]圖8 完全分離的定義[28]

    2018年Maria等人[30]基于Bus等人[28]算法做了進一步的改進,在聚類時考慮了可見性.他們提出在構(gòu)造完全分離對分解(well-separated pair decom-position, WSPD)結(jié)構(gòu)時直接查詢點對之間的可見性,這樣雖然增加了WSPD結(jié)構(gòu)的構(gòu)造時間,但總體上節(jié)約了45%左右的繪制時間,效果提升明顯.

    3.2 基于矩陣的光源聚類

    這一類方法,將多光源場景的繪制問題用光能傳輸矩陣(如圖9所示)來形式化地表示:矩陣的行索引指向著色點,列索引指向光源,第i行第j列的元素表示第j個光源對第i個著色點的光能貢獻.第i行所有元素之和代表了所有光源對第i個著色點的總貢獻.樸素的蠻力算法需要計算矩陣中每個元素的值,然后將它們沿著列方向進行累加以得到每個著色點的光照.當(dāng)存在m個著色點、n個光源時,蠻力算法所需的時間復(fù)雜度為O(mn).為了減少時間成本,這類方法通常的做法是僅計算矩陣的部分元素再估計矩陣的其余部分,而不是直接計算全部元素.需要注意的是,此類方法常常需要占據(jù)較多的存儲空間,當(dāng)像素數(shù)量非常大時,往往不太適用.

    2007年的MRCS(matrix row-column sampling)算法[32]是首個基于矩陣的多光源繪制方法.他們觀察到,由于著色點以及光源的空間連貫性,光能傳輸矩陣往往是一個低秩矩陣(例如圖9中大部分區(qū)域都是黑色),它的信息可以通過計算低維矩陣來近似.算法的主要思路如圖10所示:首先對傳輸矩陣的行做降采樣,把這些行合并為一個簡化矩陣(reduced matrix),然后對這個簡化矩陣的列索引(即對應(yīng)的光源)做聚類.根據(jù)傳輸矩陣的低秩假設(shè),簡化矩陣的列聚類結(jié)果可以近似為傳輸矩陣的列聚類結(jié)果.接下來在每一列簇中隨機選擇一個代表列進行繪制(繪制結(jié)果需要乘以列簇大小)來近似整個列簇的值.最后通過累加這些代表列即可得到最終繪制結(jié)果.該算法的優(yōu)點在于,部分計算步驟可以用陰影圖方法實現(xiàn),通過GPU運算提高了計算效率.

    Fig. 10 Overview of MRCS algorithm [32]圖10 MRCS算法總覽[32]

    Fig. 11 Algorithm overview of ref [33]圖11 文獻[33]的算法總覽

    將MRCS與Lightcuts算法進行比較: 前者讓所有著色點共享相同的光源聚類,平攤了聚類計算量,但容易忽視只對局部區(qū)域有貢獻的光源,難以捕捉光照細節(jié);后者對每個著色點都分別進行光源聚類,雖然能捕捉局部光照細節(jié),但存在不必要的重復(fù)計算,時間成本高.2011年提出的LightSlice算法[31]結(jié)合了兩者的思想,先對所有著色點劃分一個粗糙的光源聚類,然后以此為起始點,對每一部分著色點繼續(xù)細分光源聚類.

    2015年Huo等人[33]通過求解一個矩陣恢復(fù)問題來加速聚類之后的計算效率.他們觀察到,即使采用之前的LightSlice等方法進行聚類,最后要繪制的矩陣仍然很大,需要不少的計算時間.由于可見性測試是整個計算過程中最耗時的部分,他們決定從這部分入手.他們提出只對一小部分的矩陣元素的可見性進行精確計算,其余部分采用矩陣分離技術(shù)[34]恢復(fù)得到.具體的算法流程如圖11所示: 1)將矩陣按行聚類得到一系列切片(slice),對每個切片按列聚類,在每個列簇中選擇一個代表列合并成簡化矩陣.計算簡化矩陣中每個元素不包含可見性的光照值,由于不包含可見性,這部分的計算很快.2)離散采樣一些元素進行可見性測試,并利用可見性的局部性特點大致估計其余未采樣元素的可見性.3)得到的矩陣可看作是低秩目標(biāo)矩陣(希望恢復(fù)的矩陣)與稀疏誤差矩陣(估計值與實際值之差)的和,通過矩陣分離算法,恢復(fù)得到目標(biāo)矩陣.

    在VPLs框架中,為了生成云、霧等體繪制效果,一些工作將虛擬點光源替換為虛擬線光源(virtual ray lights, VRLs)[35]或虛擬束光源(virtual beam lights, VBLs)[36].學(xué)者們提出了LightSlice的改進算法[37]和自適應(yīng)的矩陣恢復(fù)方法[38]以應(yīng)對這類光源.

    表2從存儲占用、GPU友好性、對一般情形的可泛化性這3個方面對基于樹的光源聚類方法和基于矩陣的光源聚類方法進行對比.相比基于樹的算法,基于矩陣的聚類算法由于需要構(gòu)建光能傳輸矩陣,內(nèi)存占用量大.但由于傳輸矩陣的一整列能通過陰影圖的方法進行繪制,基于矩陣的方法GPU實現(xiàn)方便.基于矩陣的方法往往需要提前知曉并存儲所有著色點,可泛化程度較低.

    Table 2 Analysis of Light Clustering Methods表2 光源聚類方法比較

    總的來說,本節(jié)所介紹的方法均沒有使用重要性采樣,這樣雖然減少了噪聲的引入、降低了時間成本,但造成了圖像結(jié)果的不準(zhǔn)確性和有偏性.光源聚類也可以用來加速重要性采樣的概率分布的計算,從而提高計算效率,這部分工作將在第4節(jié)進行介紹.

    4 基于光源聚類的重要性采樣

    當(dāng)光源數(shù)量比較多時,如果對每個著色點都遍歷一遍所有光源來計算采樣概率分布,會花費巨大的時間成本.對此,許多工作將光源聚類之后再進行重要性采樣,使得直接光照采樣過程由2步變?yōu)?步:

    1) 采樣一個光源簇C;

    2) 在光源簇C的內(nèi)部采樣一個光源l;

    3) 在光源l上采樣一個點y.

    由于光源簇的數(shù)量遠少于光源數(shù)量,建立光源簇的采樣概率分布只需要消耗亞線性的時間復(fù)雜度;位于同一光源簇內(nèi)部的光源往往很相似,因此在光源簇內(nèi)部采樣光源可以使用較簡單的采樣策略,例如均勻采樣或只根據(jù)光源發(fā)射光亮度采樣.如此一來,通過聚類的方法,對每個著色點進行光源采樣的時間復(fù)雜度與光源數(shù)量呈亞線性,實現(xiàn)了伸縮性.

    為了在直接光照采樣中同時考慮光源光亮度(這里代指光源發(fā)射光亮度與幾何項的乘積)和BRDF的信息,常常使用1997年由Veach[39]提出的多重重要性采樣 (multiple importance sampling, MIS)方法:分別按照光源光亮度和著色點BRDF的分布進行直接光照的采樣計算,再將結(jié)果進行加權(quán)平均.這種方法雖然實現(xiàn)簡單方便,但會造成采樣的浪費.

    2007年Akerlund等人[40]在計算光源簇的重要性時直接考慮了BRDF信息.他們假設(shè)著色點的BRDF值在同一個光源簇內(nèi)的變化非常平滑且可以忽略,從而通過在光源簇內(nèi)隨機選擇一個光源計算BRDF值來作為權(quán)重,乘以該光源簇的光亮度估計,作為采樣的重要性.這種隨機選取光源簇內(nèi)部光源進行BRDF信息估計的方法計算量雖小,但精確度較低.

    Fig. 12 Detect the number of BRDF sample intersects for each cluster [41]圖12 檢查每個光源簇與BRDF采樣光線發(fā)生相交的數(shù)量 [41]

    BRDF常常具有高頻特性,在光源簇內(nèi)平滑變化的假設(shè)往往難以成立.2009年Wang等人[41]提出了更加實際的結(jié)合BRDF信息的光源簇采樣方法.如圖12所示,計算光源簇的重要性時,在著色點位置按照BRDF的分布采樣若干條射線,記錄與每個光源簇所在包圍盒發(fā)生相交的射線數(shù)量,以此近似BRDF在光源簇內(nèi)的積分,將其再除以包圍盒相對于著色點的立體角,得出BRDF在光源簇內(nèi)的均值,最后乘以光源簇的光亮度估計,作為該光源簇的重要性.這種BRDF的估計方法,相比Akerlund等人[40]的方法更加精確,但它依賴于BRDF采樣射線的命中率,并且難以拓展到在光源簇內(nèi)部采樣光源的過程中.

    2013年提出的VisibilityCluster方法[42]使用基于矩陣的聚類思想,在光源簇采樣中考慮可見性因素.他們按圖13所示的步驟來估計可見性:首先,根據(jù)幾何相似性分別對光源和著色點進行聚類;然后,隨機進行可見性測試,估計每個小矩陣塊的平均可見性;再通過分離光源簇,進一步細分矩陣塊;在估計得到可見性信息之后,結(jié)合前面上述采樣方法[41],即可在光源簇采樣中同時考慮光源光亮度、BRDF和可見性3項.然而,基于矩陣的聚類思想需要占用大量的存儲,難以適用于大規(guī)模場景的繪制.

    Fig. 13 Construction of VisibilityClusters[42]圖13 可見性簇的構(gòu)建[42]

    2006年Donikian等人[43]提出在繪制過程中迭代漸近地學(xué)習(xí)光源簇的采樣概率分布.他們首先用Lightcuts算法[20]對光源聚類,并將圖像空間劃分為若干塊.在繪制過程中,通過極大似然估計迭代地學(xué)習(xí)光源簇對每個塊以及每個像素的重要性,形成相對塊的采樣概率分布和相對像素的采樣概率分布.每一次迭代后,按照一定的經(jīng)驗性權(quán)重,將相對塊的采樣分布、相對像素的采樣分布以及均勻采樣分布進行加權(quán)平均,作為下一次迭代中光源簇的采樣概率密度函數(shù).對于光源簇內(nèi)部,他們根據(jù)光源發(fā)射亮度進行采樣.對于常見的光源,其發(fā)射亮度與著色點無關(guān),只需計算1次,不影響算法的可伸縮性.該算法的主要缺點在于,3個概率分布的結(jié)合權(quán)重是經(jīng)驗性的設(shè)定,不夠魯棒和可靠.

    2016年Vévoda等人[44]將上面的思想從圖像空間擴展到場景空間:他們將整個繪制場景劃分為均勻的網(wǎng)格,對每個網(wǎng)格塊單獨進行自適應(yīng)的光源聚類,并以此進行重要性采樣.為了避免不必要的計算,他們延遲到繪制過程中才進行聚類:對每個著色點,判斷其所屬網(wǎng)格塊是否已經(jīng)進行了光源聚類,如果沒有則立即進行聚類操作.當(dāng)區(qū)域塊包含聚類結(jié)果后,簡單估計每個光源簇到著色點的光能貢獻,以此生成離散的概率分布來采樣光源簇.

    2016年P(guān)harr等人[4]也給出了一種針對多光源采樣問題的實現(xiàn).他們也將繪制場景劃分為均勻的網(wǎng)格,與前面方法不同的是,他們沒有選擇聚類,而是在繪制過程中自適應(yīng)地構(gòu)建對所有光源的采樣概率分布:在著色點所屬區(qū)域塊中隨機采樣一些點,簡單計算每個光源到這些采樣點的光能貢獻,從而形成對該區(qū)域塊的光源采樣分布.該算法實現(xiàn)起來很簡單,但由于沒有進行光源聚類,當(dāng)光源數(shù)量增多時,每個區(qū)域快的光源采樣分布的構(gòu)造都需要花費大量時間,且采樣沒有考慮可見性以及BRDF的信息.

    2018年Vévoda等人[45]繼續(xù)研究改進他們在2016年提出的方法[44],提出使用貝葉斯回歸模型在線學(xué)習(xí)每個光源簇的采樣重要性.與之前相同,他們首先使用類似Lightcuts[7]的算法將光源用二叉樹組織起來,然后對場景區(qū)域均勻劃分,在繪制過程中對每個場景塊自適應(yīng)地進行光源聚類.不同的是,類似Donikian等人的方法[43],他們也在繪制過程中漸近地學(xué)習(xí)光源簇的采樣概率分布.

    1) 通過理論推導(dǎo),他們發(fā)現(xiàn)光源簇的采樣重要性既與光源簇內(nèi)部光源采樣點的光能貢獻的均值有關(guān),又與光能貢獻的方差有關(guān).

    2) 他們觀察到,不考慮BRDF信息和采樣點與著色點之間的距離時:若可見性為1,光源簇內(nèi)采樣點的光能貢獻近似平滑,可以用高斯函數(shù)來擬合;若可見性為0,光能貢獻為常值0,可以用狄拉克δ函數(shù)來擬合.因此,光源簇對著色點的光能貢獻可以建模為高斯函數(shù)與狄拉克δ函數(shù)的混合模型.

    3) 對每個著色點和每個光源簇:參考Lightcuts算法對光源簇簡單估計的重要性作為先驗,用當(dāng)前著色點所屬場景塊內(nèi)記錄的其他已繪制的著色點的光照信息作為后驗,通過最大后驗估計推斷模型得到均值與方差,進而計算出光源簇的采樣重要性.

    該方法的主要優(yōu)勢有2點:一是估計光源簇的重要性時,不僅考慮了光能貢獻的均值,還考慮到光能貢獻的方差,即內(nèi)部光源差異性越大的光源簇應(yīng)該更高概率地被采樣到;二是通過漸進式學(xué)習(xí),光源簇的重要性中包含了可見性信息,在遮擋關(guān)系復(fù)雜的場景中也能夠取得很好的繪制效果.

    同年,Estevez等人[46]從另一個角度對多光源的重要性采樣問題進行研究.算法的基本步驟非常簡潔:對光源建樹,然后對每個著色點采樣一個光源.他們的貢獻主要有3點:

    1) 對光源進行層次結(jié)構(gòu)的組織時,不再采用以往方法通常使用的Lightcuts[7]建樹方法,而是提出了表面積朝向準(zhǔn)則(surface area orientation heuristic, SAOH)來引導(dǎo)光源層次結(jié)構(gòu)的構(gòu)建,使得相似空間位置和朝向的光源位于同一節(jié)點之下,生成的光源樹質(zhì)量更好,并且能夠很好地適用于網(wǎng)格光源(mesh lights).

    2) 相比計算所有光源的重要性,構(gòu)建概率密度函數(shù)進行采樣,他們提出通過遍歷光源樹來采樣光源,即從根節(jié)點出發(fā)估計左右子節(jié)點的重要性,隨機采樣一個子節(jié)點,如此遞歸直到抵達葉節(jié)點.這樣一來,光源采樣消耗的時間與光源數(shù)量從線性關(guān)系降低為對數(shù)關(guān)系,實現(xiàn)了可伸縮性.

    3) 對于深度較低的光源樹節(jié)點,其包含的光源差異可能較大,由此估計的重要性程度會存在不準(zhǔn)確性;另外,當(dāng)著色點位于某個節(jié)點包圍盒內(nèi)部時,也無法較好地估計該節(jié)點的重要性程度.對于這些問題,他們提出在遍歷光源樹時,先不采樣子節(jié)點,而是進行分離(split)操作,同時對左右子節(jié)點進行遍歷.只有當(dāng)節(jié)點的光能貢獻誤差小于一定閾值,且著色點位于節(jié)點包圍盒之外時,才停止分離操作開始子節(jié)點的采樣,如圖14所示.其實從另一個角度來看,他們采取的分離操作等價于尋找一個光源割并對光源進行了聚類.不同的是,他們不是隨機選擇一個光源簇再選擇一個光源,而是對每個光源簇都分別采樣一個光源,計算光能貢獻并相加.從這樣的角度來看,該方法提出的遍歷樹采樣方法是一種光源簇內(nèi)采樣光源的策略.

    2019年“Ray Tracing Gems”中[47]介紹了該算法的GPU實現(xiàn),并比較分析了多種光源樹構(gòu)建方法.

    2019年Liu等人[48]基于Estevez等人的工作,提出在光源采樣中考慮BRDF信息.他們用現(xiàn)有的面光源解析方法[49-50]估計BRDF在光源樹節(jié)點對應(yīng)包圍盒的立體角下的積分,以此作為該光源樹節(jié)點的BRDF重要性權(quán)重.他們的方法可以有效解決在同時包含多光源和光澤材質(zhì)的場景中傳統(tǒng)BRDF采樣方法失效的問題.然而,由于在光源上進行點的采樣時無法考慮BRDF信息,他們的方法仍然需要與BRDF采樣方法進行結(jié)合.

    Fig. 14 Importance sampling with adaptive tree splitting[46]圖14 帶有自適應(yīng)樹分離的重要性采樣[46]

    基于光源聚類的重要性采樣方法,近年來已經(jīng)成為離線繪制中處理多光源場景的主流手段.本文之前介紹的加快可見性測試、光源聚類進行直接光照計算等工作,都可對光源聚類的重要性采樣算法起到貢獻作用.例如,加快可見性測試方面的工作可用于加快光源(簇)可見性項的重要性估計;光源聚類進行直接光照計算的相關(guān)工作可以用于構(gòu)建更好的光源樹和光源割,提高重要性采樣的效率.

    表3對本節(jié)介紹的光源采樣方法進行了比較.可以發(fā)現(xiàn),在采樣時考慮可見性項的方法往往是通過預(yù)處理或者數(shù)據(jù)驅(qū)動地在線學(xué)習(xí).而考慮BRDF項的方法仍然具有很大的改進空間.如何在采樣過程中又快又好地同時考慮光源的可見性和著色點的BRDF,仍然是一個值得探討的話題.

    Table 3 Analysis of Lights Sampling Methods表3 光源采樣方法比較

    5 實時繪制中的光源剔除(light culling)

    實時繪制往往將光源的影響范圍看作是有限的,這里將簡要介紹通過光源剔除來加速繪制速率的方法,感興趣的讀者可以查閱2015年的SIGGRAPH課程[51]進行深入了解.

    早先為了處理多光源場景,許多工作[52-54]提出使用延遲著色法(deferred shading)來替代正向著色法(forward shading),其分為2步操作:在幾何處理階段用G緩存(G-Buffer)存儲可見片元的幾何信息;在繪制階段,遍歷每個光源,對其影響范圍內(nèi)的每個可見片元進行繪制.由于每次繪制一個片元在一個光源下的光照時都需要讀寫G-Buffer,這類方法會占用非常高的帶寬.

    為了解決高帶寬占用的問題,一些工作[55-57]提出切片式延遲著色法(tiled deferred shading或deferred lighting),將原本繪制階段的內(nèi)外循環(huán)顛倒:對G-Buffer中的每個片元,查找可能受影響的光源,分別進行繪制.為了給片元分配受影響的光源,他們選擇將屏幕空間用網(wǎng)格劃分,每個網(wǎng)格共享相同的光源分配結(jié)果,從而攤銷計算時間.光源分配的方法為:分別計算每個光源投影到屏幕空間的包圍盒,對每個包圍盒內(nèi)所有網(wǎng)格區(qū)域記錄對應(yīng)的光源索引.

    為了避免延遲著色法本身所具有的高存儲要求,切片式的思想被運用到正向著色法,出現(xiàn)了切片式正向著色法(tiled forward shading, forward+)[58-59].

    2012年Harada[60]提出在使用切片式著色法時對每個切片保守地估計最小最大深度,裁去影響區(qū)域為空白的光源.

    2012年Olsson等人[61]將切片式著色法擴展為聚類著色法(clustered shading).他們發(fā)現(xiàn),原先的切片式方法是在屏幕空間進行光源分配,由于不同視角下屏幕空間的光源密度不同,繪制時間充滿不確定性.對此,他們提出為原本屏幕空間的切片(tile)引入深度,將視錐體用網(wǎng)格劃分(屏幕坐標(biāo)方向均勻劃分,深度方向?qū)?shù)劃分),然后對每個非空的網(wǎng)格單元進行光源分配,如圖15所示.之后,一些工作對具體的實現(xiàn)方法進行了探索和改進[62-63].

    Fig. 15 View frustum in tiled shading and clustered shading[63]圖15 切片式著色和聚類著色中的視錐體[63]

    2017年O’Donnell等人[64]提出通過沿深度方向建立一維的完全二叉樹來存儲每個屏幕切片或者視錐體網(wǎng)格塊中的光源,并且給出了GPU友好的構(gòu)造和遍歷算法.該算法與聚類著色法結(jié)合,可以很好地彌補視錐體后方網(wǎng)格塊深度范圍過大、光源數(shù)量過多的問題.

    上面的光源剔除算法只考慮有限的光源影響范圍,這會造成圖像結(jié)果的有偏性.2016年Tokuyoshi等人[65]提出隨機光源采樣方法,使得繪制結(jié)果保持無偏.他們分別為每個光源設(shè)定一個固定的隨機數(shù),通過俄羅斯輪盤賭(Russian roulette)決定每個光源的影響范圍,然后使用現(xiàn)有的光源剔除方法進行繪制.

    表4給出了本節(jié)介紹的5種光源剔除方法的特點.可以看出,傳統(tǒng)的延遲著色法會占用較多帶寬.相比切片式著色,聚類式著色不太依賴視角,不同視角的效果相近.另外,所有的延遲著色類算法都需要花費較高的存儲.

    Table 4 Analysis of Lights Culling Methods表4 光源剔除方法比較

    6 總結(jié)與展望

    本文對近年來離線繪制和實時繪制中的多光源繪制方法進行了綜述.

    在離線繪制方面,本文首先介紹了通過加速可見性測試以提高單個光源繪制效率的一類方法.這類方法通常可以運用在直接光照采樣的重要性估計中;然后介紹了基于VPLs框架的光源聚類方法,可以實現(xiàn)計算時間相對光源數(shù)量的可伸縮性;最后介紹了使用光源聚類思想的重要性采樣方法,這方面的工作往往受益于前2個領(lǐng)域的工作成果.如何在采樣中高效地同時考慮光源的可見性和著色點的BRDF仍然是一個開放性的問題.

    在實時繪制方面,本文介紹了光源剔除方面的工作.對于這一部分,如何結(jié)合GPU硬件更好地優(yōu)化具體的實現(xiàn)算法,往往起到非常關(guān)鍵的作用.

    針對多光源繪制問題,未來的研究重點可能是3個方向:

    1) 如何建立更優(yōu)的光源樹.提出更好的建樹啟發(fā)準(zhǔn)測,嘗試自底向上建立光源樹.

    2) 如何制定更優(yōu)的重要性采樣策略,盡可能地包含BRDF項和可見性項.一方面繼續(xù)優(yōu)化光源簇的采樣方法,另一方面也要考慮如何高效地在光源簇內(nèi)部采樣光源.

    3) 如何利用空間連貫性,讓相似的著色點共享重要性采樣的概率分布,攤銷計算量.

    猜你喜歡
    著色光源繪制
    Art on coffee cups
    蔬菜著色不良 這樣預(yù)防最好
    蘋果膨大著色期 管理細致別大意
    《光源與照明》征訂啟事
    光源與照明(2019年3期)2019-06-15 09:21:04
    10位畫家為美術(shù)片著色
    電影(2018年10期)2018-10-26 01:55:48
    放學(xué)后
    童話世界(2018年17期)2018-07-30 01:52:02
    綠色光源
    科技連載:LED 引領(lǐng)第三次照明革命——與傳統(tǒng)照明的對比(一)
    在轉(zhuǎn)變中繪制新藍圖
    Thomassen與曲面嵌入圖的著色
    99久久无色码亚洲精品果冻| 亚洲五月天丁香| 99精品久久久久人妻精品| 国产精品亚洲av一区麻豆| 精品第一国产精品| 两个人的视频大全免费| 亚洲av美国av| 麻豆av在线久日| 欧美黄色淫秽网站| 国产亚洲精品一区二区www| 亚洲精品国产一区二区精华液| 日本一本二区三区精品| 色老头精品视频在线观看| 久久欧美精品欧美久久欧美| 一二三四社区在线视频社区8| 国产亚洲精品久久久久5区| 日韩有码中文字幕| 最近最新中文字幕大全电影3| 99国产精品一区二区三区| 国产成人欧美在线观看| 亚洲一区二区三区不卡视频| 神马国产精品三级电影在线观看 | 亚洲欧美日韩高清在线视频| 亚洲欧美精品综合一区二区三区| 国产又色又爽无遮挡免费看| 欧美日韩乱码在线| 久久久久国产一级毛片高清牌| 看黄色毛片网站| 精品国产亚洲在线| 91麻豆精品激情在线观看国产| 制服丝袜大香蕉在线| 99久久国产精品久久久| 色综合婷婷激情| 一进一出抽搐gif免费好疼| 欧美黑人精品巨大| av超薄肉色丝袜交足视频| 国产成人欧美在线观看| 又黄又爽又免费观看的视频| 99久久精品热视频| 日韩高清综合在线| 日韩大尺度精品在线看网址| 熟女电影av网| 婷婷六月久久综合丁香| 18美女黄网站色大片免费观看| 一卡2卡三卡四卡精品乱码亚洲| 无限看片的www在线观看| 神马国产精品三级电影在线观看 | 男女之事视频高清在线观看| 首页视频小说图片口味搜索| 欧美日本视频| 久久人人精品亚洲av| 欧美性猛交╳xxx乱大交人| 一本一本综合久久| 麻豆一二三区av精品| 一区二区三区高清视频在线| 亚洲精品色激情综合| 禁无遮挡网站| 一级毛片高清免费大全| 久久国产乱子伦精品免费另类| 国产av不卡久久| 欧美色视频一区免费| 黄色a级毛片大全视频| 亚洲精品在线美女| 国产精品一区二区三区四区久久| 香蕉国产在线看| 国产成人av激情在线播放| 国产99白浆流出| 亚洲熟妇中文字幕五十中出| 午夜精品久久久久久毛片777| 欧美极品一区二区三区四区| 国产精品一区二区三区四区免费观看 | 免费在线观看亚洲国产| 在线十欧美十亚洲十日本专区| 亚洲欧美日韩无卡精品| 免费看十八禁软件| 午夜久久久久精精品| 午夜两性在线视频| 国产精品九九99| 黄色片一级片一级黄色片| 一级毛片女人18水好多| 亚洲真实伦在线观看| 久久中文字幕一级| 校园春色视频在线观看| 亚洲国产高清在线一区二区三| 男女床上黄色一级片免费看| 非洲黑人性xxxx精品又粗又长| 久久久久久久久中文| 好男人电影高清在线观看| or卡值多少钱| 精品久久久久久成人av| 精品少妇一区二区三区视频日本电影| 亚洲成av人片免费观看| 亚洲国产精品久久男人天堂| 国产精品久久久久久久电影 | 欧美3d第一页| 精品福利观看| 婷婷六月久久综合丁香| 免费搜索国产男女视频| 性欧美人与动物交配| 欧美性长视频在线观看| 丁香欧美五月| 在线观看免费日韩欧美大片| 久久午夜综合久久蜜桃| 18禁观看日本| 欧美又色又爽又黄视频| av中文乱码字幕在线| 一个人观看的视频www高清免费观看 | 大型黄色视频在线免费观看| ponron亚洲| 91麻豆精品激情在线观看国产| 久久精品影院6| 久久中文看片网| 久久亚洲精品不卡| 国内揄拍国产精品人妻在线| av欧美777| 亚洲无线在线观看| 亚洲成人久久爱视频| 制服人妻中文乱码| 在线十欧美十亚洲十日本专区| 国产精品综合久久久久久久免费| АⅤ资源中文在线天堂| 亚洲欧美日韩高清专用| 国产人伦9x9x在线观看| 亚洲成人免费电影在线观看| 少妇被粗大的猛进出69影院| 亚洲国产欧美人成| 国产69精品久久久久777片 | 日韩国内少妇激情av| 我的老师免费观看完整版| 国产精品一及| a级毛片在线看网站| 亚洲第一欧美日韩一区二区三区| 成人高潮视频无遮挡免费网站| 日本五十路高清| 亚洲精品美女久久av网站| 中亚洲国语对白在线视频| 狂野欧美激情性xxxx| 亚洲午夜理论影院| 日韩欧美一区二区三区在线观看| 欧美在线一区亚洲| 亚洲一区高清亚洲精品| 无限看片的www在线观看| 国产男靠女视频免费网站| 国产精品日韩av在线免费观看| 久久久国产成人免费| 国产av一区在线观看免费| 国内精品一区二区在线观看| 97人妻精品一区二区三区麻豆| 嫩草影院精品99| 亚洲成人久久爱视频| 看免费av毛片| ponron亚洲| 亚洲熟妇熟女久久| 黄色片一级片一级黄色片| xxx96com| 一级a爱片免费观看的视频| 久久中文字幕人妻熟女| av片东京热男人的天堂| 全区人妻精品视频| 国产精品一及| 久久久久久亚洲精品国产蜜桃av| 757午夜福利合集在线观看| 亚洲成人国产一区在线观看| 丰满人妻熟妇乱又伦精品不卡| 午夜激情av网站| 久久99热这里只有精品18| 亚洲欧美日韩高清专用| 久久精品亚洲精品国产色婷小说| 久久久久久久久久黄片| av福利片在线| 精品福利观看| 亚洲熟妇熟女久久| 成年女人毛片免费观看观看9| 熟妇人妻久久中文字幕3abv| 久久精品影院6| 日韩有码中文字幕| 国产视频一区二区在线看| 亚洲18禁久久av| 三级毛片av免费| 好看av亚洲va欧美ⅴa在| 午夜福利18| a级毛片在线看网站| 午夜成年电影在线免费观看| 国产精品98久久久久久宅男小说| av在线天堂中文字幕| 91麻豆av在线| 国产精品免费视频内射| 久久中文字幕人妻熟女| 国产片内射在线| 久久精品人妻少妇| 久久中文字幕人妻熟女| 99热这里只有精品一区 | 欧美不卡视频在线免费观看 | 国产精品一区二区免费欧美| 999久久久国产精品视频| а√天堂www在线а√下载| 亚洲欧美精品综合久久99| 免费电影在线观看免费观看| 在线观看免费午夜福利视频| 在线观看66精品国产| 黄片小视频在线播放| 香蕉国产在线看| 免费观看精品视频网站| 18禁黄网站禁片免费观看直播| 一个人观看的视频www高清免费观看 | 亚洲色图 男人天堂 中文字幕| 淫妇啪啪啪对白视频| 精品国产亚洲在线| 两个人视频免费观看高清| 高清在线国产一区| 三级国产精品欧美在线观看 | 精品久久蜜臀av无| 在线a可以看的网站| 丁香欧美五月| 午夜精品在线福利| 国产激情偷乱视频一区二区| 丰满的人妻完整版| 美女 人体艺术 gogo| 欧美日韩精品网址| 亚洲熟妇中文字幕五十中出| 欧洲精品卡2卡3卡4卡5卡区| 老汉色∧v一级毛片| 亚洲精品中文字幕一二三四区| 亚洲男人的天堂狠狠| 亚洲成人免费电影在线观看| 日日干狠狠操夜夜爽| 色尼玛亚洲综合影院| 又粗又爽又猛毛片免费看| 叶爱在线成人免费视频播放| 亚洲avbb在线观看| 97碰自拍视频| 1024视频免费在线观看| 国产午夜精品论理片| 一级a爱片免费观看的视频| 亚洲黑人精品在线| 搡老熟女国产l中国老女人| 日韩精品免费视频一区二区三区| 操出白浆在线播放| 午夜激情av网站| 久久久精品欧美日韩精品| 精品午夜福利视频在线观看一区| 亚洲精品在线美女| 日韩av在线大香蕉| 天堂av国产一区二区熟女人妻 | 久9热在线精品视频| 成人精品一区二区免费| 午夜a级毛片| 热99re8久久精品国产| 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产清高在天天线| 国产成人av教育| 少妇熟女aⅴ在线视频| 午夜免费观看网址| 欧美黄色淫秽网站| 一级毛片精品| 人人妻人人看人人澡| 麻豆国产av国片精品| 1024手机看黄色片| 亚洲国产欧洲综合997久久,| 在线观看免费日韩欧美大片| 99热这里只有精品一区 | 亚洲激情在线av| 久久久国产精品麻豆| 男插女下体视频免费在线播放| 色尼玛亚洲综合影院| 国产精品一区二区三区四区免费观看 | 欧美中文综合在线视频| 欧美日韩精品网址| 国产精品久久久人人做人人爽| 亚洲av第一区精品v没综合| 亚洲色图 男人天堂 中文字幕| 成人永久免费在线观看视频| 女人高潮潮喷娇喘18禁视频| 欧美日韩中文字幕国产精品一区二区三区| 午夜老司机福利片| 岛国在线观看网站| 午夜免费成人在线视频| 天堂动漫精品| 全区人妻精品视频| 一进一出好大好爽视频| 黄片大片在线免费观看| 国产成人精品久久二区二区免费| 国产男靠女视频免费网站| 日本免费a在线| 成人av在线播放网站| 久久 成人 亚洲| 美女免费视频网站| 一个人免费在线观看的高清视频| 精品久久久久久久毛片微露脸| 成年免费大片在线观看| 老司机深夜福利视频在线观看| 国产精品久久久久久亚洲av鲁大| 琪琪午夜伦伦电影理论片6080| 别揉我奶头~嗯~啊~动态视频| 午夜日韩欧美国产| 成人av一区二区三区在线看| 国产蜜桃级精品一区二区三区| 嫁个100分男人电影在线观看| 99国产精品99久久久久| 男人舔奶头视频| 亚洲精品久久成人aⅴ小说| 国产精品av久久久久免费| 国产精品美女特级片免费视频播放器 | 国产精品一及| 久久久久国内视频| 黄色丝袜av网址大全| 精品不卡国产一区二区三区| 黄色女人牲交| 1024手机看黄色片| 50天的宝宝边吃奶边哭怎么回事| 国产精品免费视频内射| 日韩欧美一区二区三区在线观看| 国产精品 国内视频| 级片在线观看| 亚洲av成人不卡在线观看播放网| 视频区欧美日本亚洲| 97超级碰碰碰精品色视频在线观看| 日日夜夜操网爽| 88av欧美| 校园春色视频在线观看| 人人妻,人人澡人人爽秒播| 欧美人与性动交α欧美精品济南到| 特级一级黄色大片| 露出奶头的视频| 国产高清videossex| 日韩欧美 国产精品| 午夜免费激情av| 欧美日韩国产亚洲二区| 日韩欧美三级三区| 无人区码免费观看不卡| 一边摸一边做爽爽视频免费| 久久久精品欧美日韩精品| 国产精品电影一区二区三区| 国产精品,欧美在线| 又爽又黄无遮挡网站| 人人妻人人澡欧美一区二区| 一级作爱视频免费观看| 精品不卡国产一区二区三区| 国产成人影院久久av| 亚洲一区二区三区不卡视频| 性色av乱码一区二区三区2| 午夜福利成人在线免费观看| 99国产精品一区二区三区| 成人国产一区最新在线观看| www日本黄色视频网| 女生性感内裤真人,穿戴方法视频| 无人区码免费观看不卡| 天天一区二区日本电影三级| 美女大奶头视频| 久久久久久久午夜电影| 国产精品免费一区二区三区在线| 男女视频在线观看网站免费 | 搞女人的毛片| 亚洲国产看品久久| 久久久久国内视频| 九色国产91popny在线| 黄色 视频免费看| videosex国产| 亚洲人成伊人成综合网2020| 国产亚洲精品av在线| 91字幕亚洲| 老汉色av国产亚洲站长工具| 久9热在线精品视频| 国产成人av激情在线播放| 全区人妻精品视频| 欧美黑人精品巨大| 久久久精品欧美日韩精品| 国产精品 国内视频| 成年免费大片在线观看| 又黄又粗又硬又大视频| 国产野战对白在线观看| avwww免费| www.www免费av| 成人亚洲精品av一区二区| 国产午夜精品论理片| 天天躁狠狠躁夜夜躁狠狠躁| 99国产综合亚洲精品| 国产欧美日韩一区二区精品| 可以在线观看毛片的网站| 人人妻,人人澡人人爽秒播| 俄罗斯特黄特色一大片| 国产成人欧美在线观看| 午夜福利成人在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 两性午夜刺激爽爽歪歪视频在线观看 | www日本黄色视频网| 麻豆国产av国片精品| 一级黄色大片毛片| 岛国在线免费视频观看| 搡老熟女国产l中国老女人| 亚洲性夜色夜夜综合| xxx96com| 久久精品国产清高在天天线| 久久伊人香网站| 午夜免费成人在线视频| 国产精品亚洲av一区麻豆| 国产高清视频在线播放一区| 亚洲五月天丁香| 久久精品亚洲精品国产色婷小说| 岛国视频午夜一区免费看| 熟女少妇亚洲综合色aaa.| 午夜福利在线在线| 青草久久国产| 白带黄色成豆腐渣| 久久人妻福利社区极品人妻图片| 成人一区二区视频在线观看| 91麻豆精品激情在线观看国产| 中出人妻视频一区二区| 男女午夜视频在线观看| 后天国语完整版免费观看| 悠悠久久av| 这个男人来自地球电影免费观看| 日韩欧美 国产精品| 久久久久久久精品吃奶| 国产熟女xx| 国产高清视频在线观看网站| 亚洲第一电影网av| 国产一区在线观看成人免费| 精品一区二区三区四区五区乱码| 天天躁夜夜躁狠狠躁躁| 欧美丝袜亚洲另类 | 日韩欧美精品v在线| 国产单亲对白刺激| 国产精品久久久人人做人人爽| 最新美女视频免费是黄的| 亚洲国产精品999在线| 午夜激情福利司机影院| 国产爱豆传媒在线观看 | 三级男女做爰猛烈吃奶摸视频| 俺也久久电影网| 国产精品国产高清国产av| 极品教师在线免费播放| 色尼玛亚洲综合影院| 国产成+人综合+亚洲专区| 国产爱豆传媒在线观看 | 日韩国内少妇激情av| 视频区欧美日本亚洲| 久久国产精品人妻蜜桃| 黄色女人牲交| 久久久久精品国产欧美久久久| 久久天堂一区二区三区四区| 搡老熟女国产l中国老女人| 欧美在线一区亚洲| 高潮久久久久久久久久久不卡| 他把我摸到了高潮在线观看| 欧美丝袜亚洲另类 | 国产午夜精品论理片| 亚洲熟妇熟女久久| 老司机福利观看| 国产精品久久久久久亚洲av鲁大| 国产成人aa在线观看| 亚洲av五月六月丁香网| 国产一区二区在线av高清观看| 最新美女视频免费是黄的| 又粗又爽又猛毛片免费看| 欧美在线一区亚洲| 一区二区三区国产精品乱码| 中文字幕人成人乱码亚洲影| 久久香蕉激情| 动漫黄色视频在线观看| 99国产精品一区二区蜜桃av| 99精品在免费线老司机午夜| 亚洲熟妇熟女久久| xxx96com| 中文字幕人成人乱码亚洲影| 热99re8久久精品国产| 亚洲欧美日韩东京热| 美女高潮喷水抽搐中文字幕| 国产v大片淫在线免费观看| 动漫黄色视频在线观看| 男女之事视频高清在线观看| 九色国产91popny在线| 欧美黑人精品巨大| 午夜激情福利司机影院| 亚洲成人免费电影在线观看| 老鸭窝网址在线观看| 日日摸夜夜添夜夜添小说| 国产高清激情床上av| 中文字幕精品亚洲无线码一区| 999久久久精品免费观看国产| 日本免费a在线| 成人一区二区视频在线观看| 日韩欧美三级三区| 免费看日本二区| 每晚都被弄得嗷嗷叫到高潮| av视频在线观看入口| 欧美 亚洲 国产 日韩一| 日本免费a在线| 熟妇人妻久久中文字幕3abv| 在线观看午夜福利视频| 久久香蕉精品热| 亚洲人成电影免费在线| 亚洲av成人精品一区久久| 日韩精品免费视频一区二区三区| 久久伊人香网站| 九色成人免费人妻av| 久久久久久久久中文| 国内毛片毛片毛片毛片毛片| 亚洲av成人一区二区三| 久久精品91无色码中文字幕| 久久精品91蜜桃| 热99re8久久精品国产| 国产亚洲精品第一综合不卡| 久久伊人香网站| 91麻豆精品激情在线观看国产| 成人精品一区二区免费| 欧美性猛交╳xxx乱大交人| 777久久人妻少妇嫩草av网站| 精品久久久久久久末码| 变态另类成人亚洲欧美熟女| 亚洲成人久久爱视频| 成人国产一区最新在线观看| √禁漫天堂资源中文www| 男女下面进入的视频免费午夜| 高潮久久久久久久久久久不卡| 女生性感内裤真人,穿戴方法视频| 不卡av一区二区三区| 欧美zozozo另类| 身体一侧抽搐| 99久久国产精品久久久| 91字幕亚洲| 国产成人啪精品午夜网站| av欧美777| 国产探花在线观看一区二区| 亚洲av日韩精品久久久久久密| 又黄又爽又免费观看的视频| 麻豆av在线久日| 波多野结衣高清无吗| 免费在线观看亚洲国产| 免费av毛片视频| 亚洲中文日韩欧美视频| 哪里可以看免费的av片| 国产激情偷乱视频一区二区| 色综合站精品国产| 老司机在亚洲福利影院| 午夜免费激情av| 国产av在哪里看| 久久久国产成人精品二区| 久久久久久久久久黄片| 窝窝影院91人妻| 亚洲av熟女| a级毛片a级免费在线| 欧美乱色亚洲激情| 啪啪无遮挡十八禁网站| 午夜福利高清视频| 99久久精品国产亚洲精品| 国产精品一及| 丁香六月欧美| 色播亚洲综合网| 国产精品亚洲一级av第二区| 免费在线观看成人毛片| 97人妻精品一区二区三区麻豆| 99热这里只有是精品50| 国产麻豆成人av免费视频| 午夜a级毛片| 亚洲人成网站高清观看| 免费观看人在逋| 三级男女做爰猛烈吃奶摸视频| 99热6这里只有精品| 精品久久久久久,| 久久久久国产一级毛片高清牌| 日本 av在线| 国产精品野战在线观看| 亚洲七黄色美女视频| 亚洲一区二区三区色噜噜| 免费在线观看成人毛片| 99热只有精品国产| 国产1区2区3区精品| 最近视频中文字幕2019在线8| 国产精品亚洲av一区麻豆| 日本在线视频免费播放| 国产亚洲精品第一综合不卡| 男人的好看免费观看在线视频 | 亚洲乱码一区二区免费版| 免费看日本二区| 欧美av亚洲av综合av国产av| 最近在线观看免费完整版| 国产精品久久久久久亚洲av鲁大| 在线观看日韩欧美| 亚洲精品在线美女| 午夜精品在线福利| 免费在线观看影片大全网站| av欧美777| 亚洲欧美精品综合久久99| 午夜福利在线在线| 国产精品亚洲av一区麻豆| 2021天堂中文幕一二区在线观| 欧美日韩福利视频一区二区| 欧美又色又爽又黄视频| 无人区码免费观看不卡| 在线观看www视频免费| 亚洲精品中文字幕一二三四区| 亚洲av五月六月丁香网| 免费搜索国产男女视频| 夜夜看夜夜爽夜夜摸| 国产麻豆成人av免费视频| 国产一区在线观看成人免费| 久久草成人影院| 草草在线视频免费看| 国产激情欧美一区二区| 日本免费a在线| 国产精品久久久久久久电影 | 又紧又爽又黄一区二区| 两性夫妻黄色片| 人人妻人人澡欧美一区二区| 亚洲第一欧美日韩一区二区三区| 搡老岳熟女国产| 欧美成人免费av一区二区三区| 国产精品,欧美在线| 天堂av国产一区二区熟女人妻 | 又爽又黄无遮挡网站| 日韩精品免费视频一区二区三区| 午夜久久久久精精品| 美女高潮喷水抽搐中文字幕| 久久 成人 亚洲|