TW201401225A - 異質參與媒體之一點所接收光量之估計方法及所構成裝置 - Google Patents

異質參與媒體之一點所接收光量之估計方法及所構成裝置 Download PDF

Info

Publication number
TW201401225A
TW201401225A TW102120018A TW102120018A TW201401225A TW 201401225 A TW201401225 A TW 201401225A TW 102120018 A TW102120018 A TW 102120018A TW 102120018 A TW102120018 A TW 102120018A TW 201401225 A TW201401225 A TW 201401225A
Authority
TW
Taiwan
Prior art keywords
light
point
points
amount
projection coefficient
Prior art date
Application number
TW102120018A
Other languages
English (en)
Inventor
Pascal Gautron
Cyril Delalandre
Jean-Eudes Marvie
Original Assignee
Thomson Licensing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Thomson Licensing filed Critical Thomson Licensing
Publication of TW201401225A publication Critical patent/TW201401225A/zh

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/50Lighting effects
    • G06T15/506Illumination models

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Image Generation (AREA)

Abstract

本發明係關於異質參與媒體(10)M點(34)所接收光量之估計方法,光係由光環境所發射。為使計算時間最佳,此方法包括步驟為:對諸點第一集合之各點(401至40n),估計第一數值,代表該考慮點和沿複數特別發光方向(421至427)界限媒體(10)的第一表面(40)間之光衰減,並以球面函數之正交基礎,投射光強度還原之第一數值,估計第一投射係數;估計第二數值,代表點(34)和沿方向(421至427)的第二表面(41)間之光衰減,第二表面包括在點(34)鄰區的某些點(402,403,414,415);使用第一投射係數和光強度還原之第二數值,估計點所接收之光量。

Description

異質參與媒體之一點所接收光量之估計方法及所構成裝置
本發明係關於合成影像組合領域,尤指異質參與媒體內之光散射(擴散)模擬領域。本發明亦在現場組合之特效範圍內。
按照先前技術,為模擬參與媒體內之光擴散,諸如霧、煙、塵、雲,有不同的方法。參與媒體相當於風媒顆粒組成之媒體,與光互動,以修飾其途徑,尤其是強度。
參與媒體可分成二部份,即同質媒體,諸如水,和異質媒體,諸如煙或雲。以同質參與媒體而言,可按分析方式,計算光源所發射光的衰減。事實上,由於其同質性能,此等媒體具有參數,諸如在媒體任何點之光吸收係數或光擴散係數常數。相反地,在異質參與媒體中,光吸收和散射(擴散)性質逐點變化。模擬在此等異質參與媒體中光擴散所需之計算,成本高昂,因此不可能以分析方式計算,並使異質參與媒體所散射(擴散)之光量有現場感。此外,未受擴散的媒體(意即媒體的擴散為各向異性),被媒體擴散的光量也會因光的擴散方向,意指人員觀看媒體的方向而異。則對人員觀察媒體的各方向,必須重複計算估計所擴散光量,以獲得媒體之現實性描繪。
為產生異質參與媒體之現場顯示,有些方法進行預先計算代表異質參與媒體之某些參數。雖然此等方法例如在後期製作中,完美適合攝影棚使用,並提供良好品質的顯示,但此等方法不適於異質參與媒體的現場互應概念和組合之脈絡。此一方法例如載於微軟公司專利申請案WO 2009/003143,於2008年12月31日公告。WO 2009/003143申請案發明之目的,為現場軟體,供描繪異質媒體,並說明使用徑向基本函數之解法。惟此解法不能視為現場描繪解法,因為有些預處理操作必須在參與媒體之 離線應用,才能計算展現媒體之投射係數,用於圖像合成現場計算。
由於出現互應式模擬遊戲和應用,特別是三維度(3D),感到亟需有現場模擬方法,提供異質參與媒體之現時性顯示。
本發明之目的,在於克服先前技術上此等缺點至少其一。
更具體而言,本發明之目的在使異質參與媒體內組成光擴散之現實性現場顯示所需計算時間和/或所需計算功率達最佳。
本發明係關於一種異質參與媒體(10)之M點所接收光量之估計方法,光係由包括複數光源的光環境所發射。方法包括步驟為:選擇媒體之某些點,形成媒體諸點之第一集合,第一集合由M點構成;為第一集合內之各點,估計光強度折減之複數第一數值,代表所考慮點和第一表面間沿複數特別發光方向之光衰減,第一表面界限媒體;估計光強度折減之複數第二數值,代表M點和第二表面間沿複數特別發光方向之光衰減,第二表面包括屬於第一集合的M點鄰區某些點;利用投射球面函數正交基礎中光強度折減之第一數值,為第一集合之各點,估計第一投射係數,代表在所考慮點之光強度折減;使用所估計第一投射係數,和光強度折減之複數第二數值,估計M點所接收之光量。
按照特定特徵,本方法包括步驟為,使用對第二表面諸點所估計第一投射係數,估計屬於第二表面的第一集合內各點所接收光量。
本方法宜包括步驟為,估計第二表面和複數特別發光方向間之交叉點,當交叉點不屬於諸點之第一集合時,利用屬於諸點第一集合的第二表面至少二點所接收光量之內插法,估計交叉點所接收光量。
按照特定特徵,本方法包括步驟為,於該球面函數之正交基礎,投射光強度折減之第二數值,估計第二投射係數,代表在M點之光強度折減,M點所接收之光量,係使用所估計第一投射係數和所估計第二投射係數,加以估計。
光強度還原之第一和第二數值估計,宜經由特別發光方向之抽樣進行。
按照特定特徵,本方法包括步驟為,估計該球面函數正交基礎之第三投射係數,此第三投射係數代表光環境諸點集合之入射光度。
按照另一特徵,本方法包括步驟為,估計球面函數正交基礎之第四投射係數,此第四投射係數係代表媒體諸點第二集合之相位函數。
本方法宜包括步驟為,從M點接收之光量,估計M點擴散之光量。
按照特定特徵,M點擴散光量之估計,係使用光線行進演算法進行。
本發明又關於為異質參與媒體(10)之M點所接收光量估計所構成裝置,光係由包括複數光源的光環境所發射,此裝置包括:選擇機構,從媒體選擇若干點,形成媒體諸點之第一組合,第一組合包括該M點;估計機構,為第一集合內之各點,估計光強度折減之複數第一數值,代表所考慮點和沿複數特別發光方向的第一表面間之光衰減,第一表面界限媒體;估計機構,估計光強度還原之複數第二數值,代表M點和沿複數特別發光方向的第二表面間之光衰減,第二表面包括屬於第一集合的M點鄰區若干點;估計機構,利用以球面函數正交基礎,投射光強度還原之第一數值,為第一集合之各點,估計第一投射係數,代表在所考慮點之光強度折減;估計機構,使用所估計第一投射係數,和光強度還原之複數第二數值,估計M點所接收之光量。
按照特定特徵,本裝置包括估計機構,以球面函數正交基礎投射光強度折減之第二數值,估計第二投射係數,代表M點光強度折減,M點所接收光量,係使用所估計第一投射係數和所估計第二投射係數估計。
此裝置宜包括估計機構,以球面函數正交基礎,估計第三投射係數,第三投射係數代表若環境諸點即合之入射光度。
按照特別特徵,本裝置包括估計機構,以球面函數正交基礎,估計第四投射係數,第四投射係數係代表媒體諸點第二即合之相位函數。
本發明又關係電腦程式製品,包括程式碼指令,當程式在電腦執行時,可執行方法步驟,估計異質參與媒體M點所接收光量。
本發明亦涉及電腦可讀式儲存機構,儲存電腦可執行指令集合,以實施異質參與媒體M點所接收光量之估計方法。
2,3‧‧‧光環境
4‧‧‧界限容積
5‧‧‧裝置
10‧‧‧參與媒體
11‧‧‧光源
12‧‧‧觀看者
13‧‧‧M點
23,24,25‧‧‧光源
31,32,33‧‧‧光源
34‧‧‧M點
40‧‧‧第一表面
41‧‧‧第二表面
A21,B22‧‧‧光源照射點
51‧‧‧微處理器
52‧‧‧圖形卡
53‧‧‧顯示器裝置
54‧‧‧輸入/輸出裝置
55‧‧‧位址和資料匯流排
56‧‧‧ROM非無常性記憶器
57‧‧‧RAM隨機存取記憶器
58‧‧‧電源供應器
60‧‧‧啟動步驟
61‧‧‧選擇參與媒體第一集合之步驟
62‧‧‧估計光強度折減第一數值之步驟
63‧‧‧估計光強度折減第二數值之步驟
64‧‧‧估計基礎函數第一投射係數之步驟
65‧‧‧估計媒體在M點接收光量之步驟
520‧‧‧圖形處理單位
521‧‧‧圖形隨機存取記憶器
530‧‧‧暫存器
570,571‧‧‧代表媒體之參數
5210‧‧‧代表媒體之參數
5211‧‧‧第一投射係數
5212‧‧‧折減數值
5213‧‧‧第二投射係數
5214‧‧‧第三投射係數
5215‧‧‧第四投射係數
5216‧‧‧交叉點
5217‧‧‧擴散光量之數值
ωin 110‧‧‧入射方向
ωout 120‧‧‧擴散方向
σs‧‧‧媒體擴散係數
σa‧‧‧媒體吸收係數
C‧‧‧空間
401至40n‧‧‧M以外之點
δs‧‧‧尺寸
ω‧‧‧方向
P,Mmax‧‧‧媒體與方向在觀看者方向之交叉點
K‧‧‧媒體和入射光線間之交叉點
ω1A 211,ω2A 212,ω3A213,ω1B 221,ω2B 222,ω3B 223‧‧‧照射方向
401,402,403,410,414,415‧‧‧M鄰區各點
421,422,423,424,425,426,427‧‧‧入射方向ω
第1圖簡略繪示按照本發明特別具體例擴散光之異質參與媒體;第2A和2B圖繪示按照本發明特別具體例包括若干光源之光環境;第3圖簡略繪示按照本發明特別具體例利用第2A和2B圖光環境照明的第1圖媒體所接收和擴散光量之估計方法;第4圖簡略繪示按照本發明特別具體例在第1圖媒體諸點集合的各點投射係數之估計方法;第5圖繪示按照本發明特別實施具體例實施對接收光量估計方法所用之裝置;第6圖繪示按照本發明特別具體例在第5圖裝置內所實施於第1圖媒體一點所接收光量之估計方法。
本發明由以下參照附圖說明,即可更為明白其他特點和優點。
第1圖表示異質參與媒體10,例如雲、霧或煙。參與媒體係由許多懸浮顆粒組成的媒體,會吸收、發射和/或擴散光。最簡單的形式是,參與媒體只例如吸收光,從例如太陽等光源11接收之光。意即通過媒體10的光會衰減,衰減程度視媒體密度而定。媒體係異質,意即媒體之物理特性,諸如組成之顆粒密度,在媒體中各點不同。由於參與媒體是與入射光互應的小粒組成,意即從光源11按照方向ωin110接收的光,不但被吸收,也會被擴散。在各向同位擴散的參與媒體中,光在四面八方均勻擴散。於各向異性擴散的參與媒體,諸如第1圖所示雲10,光擴散視光的入射方向ωin110和擴散方向ωout120間之角度而定。媒體10之M點13在擴散方向ωout120擴散的光量,利用下式計算:Q(M,ω out )=D(M).σ s .p(M,ω out in ).L ri (M,ω in ) 方程式1
媒體M點13按方向ωout120擴散到達位於空間C點的觀看者12眼睛之光量,即M點所擴散被媒體10經途徑M-P所衰減之光量(P點係位於媒體10與方向ωout在觀看者12方向之交叉點),則為: 其中:σs為媒體擴散係數;σa為媒體吸收係數;σtsa為媒體消光係數;D(M)為媒體在指定點之密度,因媒體10為異質性,密度各點不同;p(M,ωoutin)為相位函數,說明來自入射方向ωin的光如何在M點的擴散方向ωout擴散:Lri(M,ωin)為M點來自入射方向ωin110的折減光密度,代表媒體10因在K-M段的光軌道衰減後,到達M點的入射光量,K係媒體10和入射光線ωin110間之交叉點,其值為:
代表沿P15至M13途徑吸收和擴散,所致折減光度衰減。
由方程式2可計算M點所擴散而達到位於方向ωout的觀看者12眼睛之光量。為計算朝方向ωout的觀看者所接收光量,必須加上位在方向ωout的媒體全部諸點,亦即位在P-Mmax段上諸點之全部貢獻,其中P和Mmax係媒體10和方向ωout120間之二交叉點。由於簡單擴散,從方向ωout 120到達P15之此項總擴散光度則為:
於此假設,涵蓋C-P途徑之光不衰減,因媒體10外側的光不衰減之故。
把位在以ωout方向的半徑上P和Mmax間全部諸點之貢獻積分,即得此總擴散光度。此等積分方程式一般無法以分析方式解答,更不用說現場估計擴散之光量。積分是使用稱為光線行進之方法,以數位方式評估。在此方法中,積分界域分立成尺寸δM之許多間距,得下列方程式:
異質參與媒體10宜為三維度要件,第1圖中以二維度表示,為了簡明起見。
按照變化例,媒體10係利用複數光源,例如1000、100,000或1,000,000光源照明。
第2A和2B圖繪示光環境2,包括若干光源23,24,25。第2A和2B圖上相同的要件使用相同的參照符號。第2A圖更特別繪示由三個光源23,24,25照射的二點A21和B22。A21點是由第一光源23沿方向ω1A211、第二光源24沿方向ω2A 212和第三光源25沿方向ω3A 213照射。B22點是由第一光源23沿方向ω1B 221、第二光源24沿方向ω2B 222和第三光源25沿方向ω3B 223照射。如此複雜光環境因擁有若干光源,造成的問題是,估計媒體入射光的計算很耗費,因為光源與媒體一點間之光方向,在媒體上各點均不同。誠然,第一光源23所發射光採取方向,對A和B不同,第二光源24所發射光採取方向,對A和B不同,而第三光源25所發射光採取方向,對A和B不同。為解決此問題,估計來自若干遠程光源的光,是按照本發明特別具體例之環境映射方向,諸如第2B圖所示為之。並不考慮一方面是A和B點,另方面是和光源23,24,25間之光的正確方向(諸如第2A圖所示),稱為環境映射之方法是把環境2的所有光源23,24,25,就A和B點言視為位在光學無限遠。如此可把光源23,24或25所發射光採取之方向,視為一致,與媒體所考慮之點A和B無關。以此方式可忽略A和B點分開距離引起之視差效應。連接A點至第一光源23之方向ω1A 211,可視為與連接B點至第一光源23之方向ω1B 221一致。同此,連接A點至第二光源24之方向ω2A 212,可視為與連接B點至第二光源24之方向ω2B 222一致,而連接A點至第三光源25之方向ω3A 213,可視為與連接B點至第三光源25之方向ω3B223一致。
按照一變化例,光環境包括二或三個光源以上,例如1,000、100,000或1,000,000光源。
第3圖表示媒體10所擴散光量之估計方法,按照本發明特別具體例,光所由來之光環境3,包括若干光源31,32,33。正如參照第1圖 所述,媒體10在M點13擴散之光,係由媒體10從光源11(或光環境3)所接收光衰減,和媒體10所接收此衰減光量之組合。起先,關於第3圖,估計方程式1中代表媒體10從光環境3所接收光衰減之項。為此,抽樣包圍M點34之球面Ω,對來自球面Ω中心M之各方向ω,使用如下方程式,相等於方程式3,估計沿包括M和媒體外部間途徑之光衰減: 其中R(M,ω)是按照方向ω在M點13的光強度衰減,並代表衰減後到達M點之入射光量。
D(s)為媒體密度;σt為媒體消光係數,相當於媒體σs擴散係數和媒體σa吸收係數之合計(σtsa);K35係遵循離開M點的方向ω,在媒體10和媒體10外部間之交叉點。
方程式6提供在指定方向ω上一點之光衰減。為估計按照方向ω在M點之光衰減,位在入射方向ω的積分界域,分立成一系列尺寸δs之間距,密度每一間距不同,因媒體10為異質之故。應用光線行進法,可在按照方向ω的M點獲得光衰減數值。此數值記錄在關聯GPU(圖形處理單位)的記憶器之列表內。對於在包括從M點出發的N方向ω集合中抽樣以M為中心之球面Ω各方向ω,重複此項估計在M點34光衰減之運算,其中N係任何自然正整數。代表按照特別方向ω在M點光衰減之此等數值,即儲存於與GPU關聯之記憶器列表內。
函數空間內之各函數,可寫成基本函數之線性組合,基本函數是函數空間之基本要素。使用球面函數之正交基礎,可利用下式在M點展示光衰減函數: 其中R(M)為在M點之光衰減函數,Crj(M)為基本函數Bj(M)(在全部Nc係數之上)第j個投射係數,其中Crj(M)係利用球面Ω上積分加以界定,意即:Cr j (M)=ʃΩ R(M,ω)B j (ω) 方程式8
如此計算的基本函數投射係數集合Nc,即儲存於GPU之記憶器列表內。此等係數代表若干光源的光環境所發射光之衰減,由於來自媒體10特別指定M點的任何方向之光單純擴散所致。
對媒體10諸M點之集合,宜重複上述運算。對諸點集合之各點,因此算出代表按照四面八方的光衰減之投射係數,並儲存在列表內,稱為衰減記錄。計算投射係數的M點數量愈多,因媒體10簡單擴散所致光衰減之代表性愈正確。
同理,媒體10之相位函數,可用球面函數正交基礎表示。媒體10的相位函數,可視為在媒體10任何點均相同。把相位函數p(ωoutin)投射入球面函數基本B的投射係數Cpj(稱為第三投射係數)集合Nc內,可得下式:Cp j (ω out )=ʃΩ p(ω out in )B(ω in ) in 方程式9
此項投射是為從球面Ω上M點出發方向ω之集合Nc進行。代表媒體相位函數之投射係數,儲存於與GPU關聯的記憶器列表內。按照變化例,相位函數是各點M不同,投射係數係對代表媒體10的諸點集合算出。
按類似方式,描述代表媒體10入射光的環境映射3之函數,係使用球面函數正交基礎表示。把代表環境映射3的函數投射於球面函數基本B的投射係數Clj(稱為第二投射係數)集合Nc內,可得下式:Cl j (ω in )=ʃΩ L i (ω in )B(ω in ) in 方程式10其中Li(ωin)為按照入射方向ωin的入射光函數。
光衰減函數、相位函數和入射光函數,宜以單一之同樣球面函數正交基礎表示。
如前就第1圖所述,方程式1載明在包括單一光源11的光環境之擴散方向ωout 120,媒體10的M點13所擴散光量。以此方程式開始,適於包括複數光源之光環境3,並藉入射光函數說明,可得如下簡單擴散方程式,說明在按照方向ωout的M點,由媒體10傳輸之光: Q(M,ω out )=ʃΩ σ s (M).p(M,ω out in ).L i (ω in )R(M,ω in ) in 方程式11
函數R、Li和p分別投射入球面基本函數投射係數集合、第三係數集合和第四係數集合內。投射係數分別標明為{Crj}j=0..n,{Clj}j=0..N和{Cpj}j=0..N。說明已為媒體10諸點集合且並非為媒體M諸點全部估計的光衰減函數之投射係數,未經方程式8估計的諸點函數R之投射係數,係利用內插法計算。欲得位在觀看者觀察方向的M點衰減函數之第一投射係數,進行為位在M鄰區的複數點(例如4,8或16點)所估計函數R第一投射函數之內插法。
使用代表媒體內光衰減之第一投射係數。來自光環境之光(藉方程式8,就第2圖所述,並利用空間內插法估計),代表入射光之第三投射係數,和代表媒體10相位函數之第四投射係數,即可藉分析方式估計觀看者36接收時M點的全局光衰減,所需計算功率資源遠比積分形式的方程式分析分解所需者少。為計算所投射函數的三重積數,使用技術專家公知之張量積數,且使用方程式2,8和11,即可估計媒體M點13所發射和觀看者12按方向ωout看所接收之光量,因此得: 其中代表張量積數,係按照下式計算:Γ jki Ω B j (ω)B k (ω)B i (ω)) 方程式13
方程式12代表M點所發射和觀看者所接收之光量。欲得位在C點的觀看者朝方向ωout120看所接收之總光量,只要加上Mi點集合從P進到Mmax所發射單元光量之和即可,得如下方程式:
欲得媒體10所擴散和觀看者12所接收總光量,對離開使用者和越過媒體10的四面八方,重複上述估計。觀看者按照各觀察方向接收之光量合計,提供觀看者12從媒體10所接收光量。
第4圖表示本發明非限制性之具體例中M點按照若干入射 方向ω 421,422,423,424,425,426,427光衰減之估計方法。容積4亦稱為界限箱,宜在參與媒體10周圍界定,即參與媒體10形成之容積,以宜為長方形或四方平行六面體形狀,或任何其他形狀之箱包圍。界限容積4可快速估計實際舞台上參與媒體10所佔有容積。使用界限容積可簡化並加速涉及參與媒體10之所有幾何形計算,界限容積之幾何形較參與媒體簡單,後者會很複雜。產生界限容積以包含參與媒體之全部幾何形,在界限容積4內表面和參與媒體外表面間存在之空間,至少在適當處。界限容積分成複數分格(例如40×40×40或50×50×50或100×100×100分格),各分格宜利用點401,402,403...,414,415,40n界定。界限容積4宜以其外表面40(以下稱第一表面)界定,包圍界限容積和參與媒體10。為決定M點在各方向ω之光衰減,要估計各方向ω與第一表面40之交叉。使用方向ω 427為例,說明在M點光衰減之決定過程。按照此例,方向ω 427與第一表面間之交叉為P1點401。利用技術專家公知之任何幾何學方法,計算與第一表面40之交叉點P1。就第4圖說明,M點與P1點401分開的距離,分立成一系列尺寸δs的間距。從方程式6和應用光線行進法,得M點按照方向ω 427之光衰減數值(和為光衰減第一數值),並例如記錄在與GPU關聯的記憶器列表內。對所有方向421至427和界定界限容積4方格之所有點401至40n,重複此運算,M點相當於諸點401至40n之一。此等第一數值可決定媒體10內由於簡單光擴散引起之光衰減,因而決定M點從光環境2接收之光量。
按照變化例,第一表面相當於參與媒體10外表面,可防止分格落到參與媒體外面,因而更佳估計媒體10所接收光量。
要更佳估計M點所接收光量,宜決定媒體10與M不同的諸點401至40n,對M點光明之貢獻。界定第二表面41,包括屬於媒體10和M點鄰區之諸點,鄰區各點接收之光量,已照M點同樣方式決定,即利用決定代表光環境沿複數方向ω所發射光衰減之數值。形成第二表面的M點鄰區諸點,例如相當於M點的最近點,即包括M點的分格在內之諸分格的諸點(尤指402,403,414,415)。第二表面界定M點周圍之第一界限容積。對於估計M點直接從光環境接收光之第一衰減數值所用各方向ω 421至427,決定第二表面41和方向ω 421至427間之交叉點。以使用方向ω 427為例,方向ω 427和第二表面41間之交叉,相當於P2點410。第一交叉點 P1 401和第二交叉點P2 410,係以涵蓋ω 427決定,一旦決定此二交叉點,宜即涵蓋方向。當交叉點P2 410不相對應於形成界限容積4的方格諸點之一,亦即已計算過所接收光量的媒體諸點,此點P2 410所接收光量,即利用屬於第二表面41的媒體10諸點414,415所接收光量之內插法決定,例如按照分開第二表面41諸點的各交叉點P2 410之距離,為此決定所接收光量,並儲存於記憶器內。決由方程式6,用於決定M點(跨越距離[MP1])的第一光衰減數值之方向ω 427分立,是用來決定M點代表在途徑[MP2]衰減的第二光衰減數值。易言之,第一和第二光衰減數值是單通決定,包含M和P1 401間通過P2 410方向ω 427一次。因此,為各方向ω 421至427決定第二光衰減數值,並儲存於與GPU關聯之記憶器列表內。
M點接收的總光量,係由代表第一光衰減數值和第二光衰減數值之第一投射係數計算。為此,把M點34利用簡單擴散所接收光量(即由媒體在第一表面和M點34間所包括途徑上接收並衰減之光量),和多次擴散之光量合計。多次擴散相當於鄰區諸點貢獻之合計,並由媒體10在此等諸點和M點34間所包括途徑上所衰減第二表面各點之簡單擴散相加而得。由如下方程式得總光量: 其中Pi相當於M點鄰區之第i點,即第二表面之第i點。
又為改進M點所接收光量之估計,宜按上段所述同樣方式,決定屬於界限第二表面之一或以上第三表面的媒體10諸點貢獻,為了決定把M點與諸方向ω和包括界限M點的媒體諸點之表面間各交叉分開的途徑關聯之所有光衰減數值,只能實施一次不同方向ω 421至427之分立。
第二光衰減數值之集合(代表第二交叉點和M點間之光衰減),形成在M點的第二光衰減函數,M點的第一光衰減函數係由第一光衰減數值集合形成。M點之第二光衰減函數係由方程式7,以表示第一光衰減函數所用之同樣球面函數正交基礎表示,Nc第二投射係數之第二集合,是由方程式8決定,第二投射係數係利用以球面函數正交基礎投射光強度折減之第二數值,代表M點的光強度折減。以第二投射係數代表M點的第 二光衰減函數,可減少儲存第二光強度折減數值所需記憶器空間,第二投射係數之Nc數量,比M點的光強度折減的第二數值之數量,少很多。
M點所接收總光量Qt,是由第一投射係數和第二投射係數估計。媒體在M點的簡單擴散,是使用技術專家公知的張量積數,由代表媒體內第一光衰減函數之第一投射係數,代表光環境的第三投射係數,和代表媒體相位函數之第四投射係數估計。M點接收的總光量Qt,是使用簡單擴散,和代表相對於鄰區的第二光衰減之第二投射係數(相當於多數擴散)估計。
M點所接收再朝方向ωout發射給觀看者12之總光量,是由如下方程式,由簡單擴散數值和多數擴散數值相加而得: 其中Q(M,ω out )相當於媒體10在M點34的簡單擴散,相當於第二係數,其Pi為屬於第二表面的鄰區一點。
第5圖簡略表示裝置5之硬體具體例,適於估計異質參與媒體10擴散之光量,並創作一或以上影像之顯示訊號。裝置5例如相當於個人電腦PC、膝上型或遊戲器控制台。
裝置5包括下列元件,利用位址和資料匯流排55連接在一起,也可以輸送時計訊號:微處理器51(或CPU);圖形卡52,包括:若干圖形處理單位520(或GPU);GRAM型隨機存取記憶器521(圖形隨機存取記憶器);ROM型(唯讀記憶器)之非無常性記憶器56;RAM隨機存取記憶器57;一或以上I/O(輸入/輸出)裝置54,諸如鍵盤、滑鼠、網路攝影機;以及電源供應器58。
裝置5亦包括顯示幕型之顯示器裝置53,直接連接到圖形卡52特別顯示圖形卡內所計算和組成的電腦產生圖形之描繪,例如在現 場。使用專用匯流排把顯示器裝置53連接到圖形卡52,其優點是有遠較高的資料傳輸率,並以此方式減少圖形卡所組成影像顯示之潛候時間。按照一變化例,顯示器裝置在裝置5之外部,並利用傳輸顯示訊號之電纜連接至裝置5。裝置5例如圖形卡52,包括傳輸媒體或連接器(第5圖上未示),適於傳輸顯示器訊號至外部顯示機構,例如LCD,或電漿顯示幕、投影機。
須知記憶器52,56,57說明所用「暫存器」,在所述各記憶器內指定為低容量記憶區(有些二進資料)和大容量記憶區(可儲存全部程式,或是代表所計算或顯示資料之全部或部份資料)。
開機時,微處理器51載錄並運轉RAM57內所含程式之指令。
隨機存取記憶器57特別包括:在暫存器530內的微處理器51運算程式,負責裝置之交換;代表異質參與媒體10之參數571(例如密度、光吸收係數、光擴散係數等參數)。
實施本發明專屬的方法步驟之下述演算法,儲存於與實施此等步驟的裝置5相關聯圖形卡52之GRAM記憶器57內。開機時,代表媒體之參數570一旦載錄於RAM 57內,圖形卡52之圖形處理單位520即把此等參數載錄於GRAM 521,並例如使用HLSL(高階漸影器語言)、GLSL(開放GL漸影語言)等「漸影器」型微程式之形式,執行此等演算法之指令。
GRAM隨機處理記憶器521特別包括:在暫存器5210內代表媒體10之參數;第一投射係數5211,代表媒體10各點之光強度第一折減數值;對媒體10各點之光強度第一和第二折減數值5212;第二投射係數5213,代表媒體10各點之光強度第二折減數值;第三投射係數5214,代表媒體10各點之入射光;第四投射係數5215,代表媒體10各點之相位函數;諸方向ω與第一和第二表面間之第一和第二交叉點5216;代表媒體10沿一或以上觀察方向所擴散光量之數值5217。
按照變化例,若GRAM 521可用之記憶器儲存空間不夠, 可由CPU 51指定RAM 57之一部份,供儲存係數5211,5213,5214,5215和數值5212,5216,5217。惟此變化例,在包括代表GPU所含微程式組成的媒體10之影像組合中,實現較長的潛候時間,因為必須利用傳輸容量一般比在圖形卡內把資料在GPU和GRAM間來回傳送所用為低之匯流排55,把資料從圖形卡傳輸到隨機存取記憶器57。
按照另一變化例,電源供應器58在裝置5之外部。
第6圖繪示在裝置5內實施異質參與媒體一點所接收光量之估計方法,按照本發明較佳非限制性之第一實施具體例。
於啟動步驟60,把裝置5之各種參數更新。具體言之,可以任何方式啟動表示異質參與媒體10之參數。
其次,於步驟61中,選擇參與媒體10諸點之第一集合。此諸點之第一集合例如相當於形成界限參與媒體10的第一容積之方格的中心。按照描繪包括媒體10的場景所需品質,決定諸點數量。諸點數量愈多,描繪品質愈佳。然而,諸點數量愈大,估計被諸點所接收光量所需計算愈多。選擇諸點數量是例如求得描繪品質和獲得描繪所需計算數量間之最佳妥協。
然後在第二步驟62,對所選擇諸點集合之各點,估計光強度之第一折減數值。於此,於媒體10的指定點M,就複數之指定入射方向,計算代表光強度折減之複數第一數值。代表光強度折減之各第一數值,相當於複數入射方向當中之指定入射方向,和媒體10沿入射方向之光途徑,即在M點34和界限媒體10的第一表面間包括之表面。對於從光環境進行到結合M點34的複數入射方向之各入射方向,使用技術專家公知之任何抽樣方向,宜用所謂光線行進演算法,計算代表M點34光強度折減之數值。據以計算代表M點光強度折減的第一數值之複數入射方向ω(或ωin),形成球面Ω,中點為M。選擇計算代表M點光強度折減的第一數值之方向ω數量,其方式在於裝置5使用者所需要,即計算此等數值所需計算功率和M點光強度折減估計準確性間之最佳妥協。選擇方向ω數量,意味抽樣以M為中心之球面Ω,宜對代表媒體10的諸點第一集合(例如50,100或1000點)之各點,重複步驟61。
接著,在步驟63估計代表沿光發射方向ω在M點34和第 二表面間光衰減之光強度折減第二數值。第二表面包括M點34鄰區諸點之一部份,屬於在步驟62估計光強度折減第一數值選用之諸點第一集合。按光強度折減之第一數值同樣方式,藉決定方向ω和第二表面間之交叉,決定光強度折減之第二數值。方向ω之分立是步驟62和63通用,惟M點34分別與第一表面和第二表面交叉之距離不同,以分別決定第一數值和第二數值。
再於第四步驟64,估計球面函數正交基礎中基礎函數之第一投射係數,此等第一投射係數代表在異質參加媒體10任何點之光強度折減。第一投射係數宜由球面函數正交基礎中光強度折減之第一數值投射的結果。由球面函數正交基礎中第一數值之投射,可知代表球面函數正交基礎中之此等第一數值。投射涉及從函數之第一基礎(第一數值起先沿入射方向表示,即在到達計算第一數值之M點以前,按照參與媒體沿入射方向所涵蓋距離),通到函數之第二基礎,即球面函數正交基礎。投射係數是球面函數正交基礎中第一數值表示(即表現)的結果。因此,在諸點第一集合之各點,得到代表光強度折減之第一投射係數。換句話說,第一投射係數之集合,相當於媒體之一點,在代表媒體10的諸點第一集合中有多少點,即計算多少集合之第一投射係數。為獲得已計算此等第一係數的媒體諸點之第一投射係數,使用所計算第一投射係數之內插法。
最後,在步驟65,使用先前估計的第一投射係數,和光強度折減之第二數值,估計媒體10在M點34接收之光量,代表M點34之鄰區諸點,對M點光亮之貢獻。
決定所接收光量,宜對諸點第一集合之各點實施。
按照變化例,利用投射球面函數基礎中光強度折減第二數值,為M點34估計代表光強度折減之第二投射係數。此等第二投射係數宜用來決定M點34所接收光量,代替迄今用來代表光衰減第一函數的光強度第二折減數值。
按照變化例,為光環境3諸點集合,估計代表入射光之第三投射係數,集合之諸點係代表光環境。於此,使用已用來代表第一光衰減函數的球面函數正交基礎,表示描述代表媒體10入射光的環境映射3的函數(稱為入射光函數)。因此,把代表環境映射3的函數,投射於球面函數 基礎之第三投射係數集合內。此等第三係數宜儲存於與裝置5的圖形卡之一或以上GPU關聯的記憶器列表內。第三投射係數在媒體任何點宜一致。按照變化例,第三投射係數各點不同,諸點之一集合與屬於光環境之另一集合也不同。
按照另一變化例,為媒體10諸點之第二集合(集合之諸點係代表光環境),估計代表媒體10相位函數之第四投射係數。於此,使用已用來代表第一光衰減函數和入射光函數之球面函數正交基礎,代表描述媒體10相位函數之函數。因此,相位函數投射於球面函數基本之第四投射係數集合。此等第四投射係數宜儲存於與裝置5的圖形卡之一或以上GPU關聯的記憶器列表內。第四投射係數在媒體任何點宜一致。按照變化例,第四投射係數各點不同,諸點之一集合與屬於媒體10之另一集合亦不同。
當然,本發明不限於前述具體例。
尤其是本發明不限於異質參與媒體所擴散光量之估計方法,亦可延伸至實施此方法之任何裝置,尤指包括至少一個GPU之所有裝置。就第1至4圖所述方程式對估計第一、第二和第三投射係數,光強度折減,所擴散光量之實施,也不限於漸影器型微程式之實施,亦可延伸到任何程式型之實施。例如可在CPU型微處理器內執行之程式。
估計投射係數所用基礎函數,宜為球諧型或球面子波型函數。
先前估計之第一和第二投射係數,宜記錄並儲存於由與GPU關聯的記憶器內所記錄列表組成之一或以上資料結構內。此等記錄亦稱為衰減記錄。衰減記錄列表宜包括媒體10諸點之第一和第二投射係數,無論此等第一和第二係數係由方程式8計算,或由經方程式8計算的第一係數內插而得。有利用媒體10諸點之第一投射係數集合,或按照變化例,為代表媒體10的諸點集合。第一和第二投射係數之儲存,優點是加速估計由媒體10擴散和觀看者所接收光量,代表入射光強度折減之第一投射係數,是任何時間和在方程式12和13中使用時立刻可得。同理,第三投射係數和第四投射係數記錄於專用列表內,可加速計算,尤其是當第三和第四係數各點不同(分別來自光環境3和媒體10)。
本發明使用不在於現場使用,亦可延伸到其他用途,例如在 記錄工作室內得後製作處理操作,例如供描繪電腦產生之圖像。本發明在後製作之實施,優點是尤其在賦予卓越之視覺描繪寫實,又能減少必要之計算時間。
本發明亦關係到媒體10之一點(或媒體10全體)在一或以上指定方向,所擴散光量之估計方法。媒體10按照發射方向120擴散之光量,係用媒體之一點或多點,例如M點34,所接收光量估計。於此,線段相當於發射方向120與媒體10之交叉,意即[PMmax]段在空間上分立成代表此段之諸多點或元件。為此段之各點(或各元件),應用方程式15或方程式16。按照變化例,實施光線行進法,以估計考慮的段之一點和在發射方向120位於媒體10周邊的P點15間之光強度折減。
本發明亦關係到視訊影像在二維度或三維度之組成方法,為此計算異質參與媒體所擴散光量,而所得代表光之資訊,即用來顯示影像之圖元,各圖元相當於按照觀察方向ωout之觀察方向。為利用影像的各圖元顯示所計算光數值,重新計算以適應觀看者之不同視點。
本發明可例如用於視訊遊戲用途,不論是利用PC或膝上型電腦內可執行的程式,或製作和顯示現場影像之專用遊戲控制台。就第5圖所述裝置5宜具有互應機構,諸如鍵盤和/或搖桿,亦可用其他模態以輸入命令,例如語音辨識。
4‧‧‧界限容積
10‧‧‧參與媒體
34‧‧‧M點
40‧‧‧第一表面
401至40n‧‧‧M以外之點
41‧‧‧第二表面
401,402,403,410,414,415‧‧‧M鄰區各點
421,422,423,424,425,426,427‧‧‧入射方向ω

Claims (15)

  1. 一種異質參與媒體(10)之M點(34)所接收光量之估計方法,光係由包括複數光源(31,32,33)的光環境(3)所發射,其特徵為,方法包括步驟為:選擇(61)該媒體之某些點,形成該媒體(10)的諸點(401至40n)之第一集合,該第一集合包括該M點(34);為該第一集合內之各點(401至40n),估計(62)光強度折減之複數第一數值(5212),代表該所考慮點和第一表面(40)間沿複數特別發光方向(421至427)之光衰減,第一表面(40)界限媒體(10);估計(63)光強度折減之複數第二數值(5212),代表該M點(34)和至少一個第二表面(41)間沿複數特別發光方向(421至427)之光衰減,至少一個第二表面包括屬於第一集合的該M點(34)鄰區某些點(402,403,414,415);利用投射球面函數正交基礎中光強度折減之該第一數值,為該第一集合之各點(401至40n),估計(64)第一投射係數(5211),代表在該點之光強度折減;使用該所估計第一投射係數,和光強度折減之複數第二數值,估計(65)該M點(34)所接收之光量者。
  2. 如申請專利範圍第1項之方法,其中該M點(34)所接收光量,相當於該光環境(3)所發射,並沿M點(34)和第一表面(40)間複數特別發光方向(421至427)衰減之光量,與至少一個第二表面(41)諸點集合再發射,並沿M點(34),和第二表面(41)間複數發光特別方向(421至427)衰減之光量合計者。
  3. 如申請專利範圍第1至2項之一項方法,其中包括步驟為,使用為第二表面(41)諸點所估計之第一投射係數,估計屬於第二表面(41)的第一集合中各點(402,403,414,415)所接收光量者。
  4. 如申請專利範圍第1至3項之一項方法,其中包括步驟為,估計第二表面(41)和該複數特別發光方向(421至427)間之交叉點(410),當該交叉點(410)不屬於諸點之第一集合,該交叉點所接收光量,係由屬於諸點第一集合的第二表面(41)至少二點(414,415)所接收光量之內插法估計者。
  5. 如申請專利範圍第1至4項之一項方法,其中包括步驟為,利用投射在該球面函數正交基礎中光強度折減之該第二數值,估計代表在M點(34)光強度折減之第二投射係數,該M點(34)所接收光量係使用該所估計第一 投射係數和所估計第二投射係數估計者。
  6. 如申請專利範圍第1至5項之一項方法,其中光強度折減的該第一和第二數值之估計,係利用該特別發光方向(421至427)之光線行進所為者。
  7. 如申請專利範圍第1至6項之一項方法,其中包括步驟為,估計該球面函數正交基礎中之第三投射係數,該第三投射係數代表該光環境(3)諸點集合之入射光者。
  8. 如申請專利範圍第1至7項之一項方法,其中包括步驟為,估計該球面函數正交基礎中之第四投射係數,該第四投射係數代表該媒體諸點第二集合之相位函數者。
  9. 如申請專利範圍第1至8項之一項方法,其中包括步驟為,從M點(34)接收之光量,估計M點(34)擴散之光量者。
  10. 一種從異質參與媒體(10)估計M點(34)接收光量所構成裝置(5),光係由包括複數光源(31,32,33)之光環境(3)所發射,其特徵為,裝置包括至少一處理器,其構成係為:選擇該媒體(10)之某些點,形成該媒體(10)的諸點(401至40n)之第一集合,該第一集合包括該M點(34);為該第一集合內之各點(401至40n),估計光強度折減之複數第一數值(5212),代表該所考慮點和第一表面(40)間沿複數特別發光方向(421至427)之光衰減,第一表面(40)界限媒體(10);估計光強度折減之複數第二數值(5212),代表該M點(34)和至少一個第二表面(41)間沿複數特別發光方向(421至427)之光衰減,至少一個第二表面包括屬於第一集合的該M點(34)鄰區某些點(402,403,414,415);利用投射球面函數正交基礎中光強度折減之該第一數值,為該第一集合之各點,估計第一投射係數(5211),代表在該點之光強度折減;使用該所估計第一投射係數,和光強度折減之複數第二數值,估計該M點(34)所接收之光量者。
  11. 如申請專利範圍第10項之裝置,其中至少一處理器又構成為,利用投射該球面函數正交基礎中光強度折減之該第二數值,估計第二投射係數,代表在M點(34)之光強度折減,該M點(34)所接收光量係使用該所估計第一投射係數,和所估計第二投射係數估計者。
  12. 如申請專利範圍第10至11項之一項裝置,其中至少一處理器又構成為,估計該球面函數正交基礎中之第三投射係數,該第三投射係數代表該光環境(3)諸點集合之入射光者。
  13. 如申請專利範圍第10至12項之一項裝置,其中至少一處理器又構成為,估計該球面函數正交基礎中之第四投射係數,該第四投射係數代表該媒體(10)諸點第二集合之相位函數者。
  14. 一種電腦程式製品,其特徵為,包括碼指令,當該程式在電腦內執行時,可供執行申請專利範圍第1至9項之一項方法者。
  15. 一種電腦可讀式儲存機構,儲存該電腦可執行之指令集合,以實施申請專利範圍第1至9項之任一項方法者。
TW102120018A 2012-06-22 2013-06-06 異質參與媒體之一點所接收光量之估計方法及所構成裝置 TW201401225A (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
FR1255940 2012-06-22

Publications (1)

Publication Number Publication Date
TW201401225A true TW201401225A (zh) 2014-01-01

Family

ID=46963861

Family Applications (1)

Application Number Title Priority Date Filing Date
TW102120018A TW201401225A (zh) 2012-06-22 2013-06-06 異質參與媒體之一點所接收光量之估計方法及所構成裝置

Country Status (6)

Country Link
US (1) US9235663B2 (zh)
EP (1) EP2677504A3 (zh)
JP (1) JP2014006901A (zh)
KR (1) KR20140000170A (zh)
CN (1) CN103514624A (zh)
TW (1) TW201401225A (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2948800A1 (fr) * 2009-11-16 2011-02-04 Thomson Licensing Procede d'estimation de diffusion de la lumiere
US20130100135A1 (en) * 2010-07-01 2013-04-25 Thomson Licensing Method of estimating diffusion of light
US9767601B2 (en) * 2014-08-28 2017-09-19 Disney Enterprises, Inc. Residual ratio tracking for estimating attenuation in heterogeneous volumes
EP3057067B1 (en) * 2015-02-16 2017-08-23 Thomson Licensing Device and method for estimating a glossy part of radiation
US10395422B2 (en) 2017-10-13 2019-08-27 Activision Publishing, Inc. Participating media baking
US11004253B2 (en) * 2019-02-21 2021-05-11 Electronic Arts Inc. Systems and methods for texture-space ray tracing of transparent and translucent objects

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3996673A (en) * 1975-05-29 1976-12-14 Mcdonnell Douglas Corporation Image generating means
US4737921A (en) * 1985-06-03 1988-04-12 Dynamic Digital Displays, Inc. Three dimensional medical image display system
US5649078A (en) * 1995-03-22 1997-07-15 International Business Machines Corporation Efficient two-pass rasterization scheme utilizing visibility information
US6891538B1 (en) * 1999-07-01 2005-05-10 Microsoft Corporation Dual mode device and method for generating vector cross products or dot products
JP3975892B2 (ja) * 2002-05-02 2007-09-12 富士ゼロックス株式会社 位置計測システム
US7924306B2 (en) 2006-09-15 2011-04-12 Hewlett-Packard Development Company, L.P. Videoconferencing with enhanced illusion of physical presence in a common space
US7990377B2 (en) * 2007-06-26 2011-08-02 Microsoft Corporation Real-time rendering of light-scattering media
US8009168B2 (en) * 2007-06-26 2011-08-30 Microsoft Corporation Real-time rendering of light-scattering media
US7940269B2 (en) * 2007-06-29 2011-05-10 Microsoft Corporation Real-time rendering of light-scattering media
FR2948800A1 (fr) 2009-11-16 2011-02-04 Thomson Licensing Procede d'estimation de diffusion de la lumiere
FR2964776A1 (fr) * 2010-09-14 2012-03-16 Thomson Licensing Procede d’estimation de diffusion de la lumiere dans un milieu homogene

Also Published As

Publication number Publication date
JP2014006901A (ja) 2014-01-16
KR20140000170A (ko) 2014-01-02
EP2677504A2 (en) 2013-12-25
CN103514624A (zh) 2014-01-15
US20130346041A1 (en) 2013-12-26
US9235663B2 (en) 2016-01-12
EP2677504A3 (en) 2014-04-16

Similar Documents

Publication Publication Date Title
TW201401225A (zh) 異質參與媒體之一點所接收光量之估計方法及所構成裝置
CN111340928A (zh) 一种结合光线跟踪的Web端实时混合渲染方法、装置及计算机设备
JP5873672B2 (ja) 仮想環境のポイントにおいて受光された光の量の推定方法
JP5873683B2 (ja) 仮想環境でのオクルージョンを推定する方法
US20210366176A1 (en) Shadow culling
JP2015114775A (ja) 画像処理装置および画像処理方法
JP5916764B2 (ja) 仮想環境における隠蔽の推定方法
US20180040155A1 (en) Cache friendly jittered hemispherical sampling
KR20120086739A (ko) 광 산란을 추정하는 방법
EP2831846B1 (en) Method for representing a participating media in a scene and corresponding device
US20120232830A1 (en) Method for estimating light scattering
US8842275B2 (en) Method for estimating light scattering
US20130100135A1 (en) Method of estimating diffusion of light
EP2428935B1 (en) Method for estimating the scattering of light in a homogeneous medium
US20150006113A1 (en) Method and device for estimating light scattering
Al-Oraiqat et al. Specialized Computer systems for environment visualization
Johansson et al. From BIM to VR
Perez Soler Volume Rendering Simulation in Real-Time
Petkov Lattice-based Immersive Visualization
Liu Efficient Algorithms for Realistic Lens Effect Simulation
Chelliah Parallel Methods for Projection on Strongly Curved Surfaces
Brownlee Parallel ray tracing in scientific visualization