最美情侣中文字幕电影,在线麻豆精品传媒,在线网站高清黄,久久黄色视频

歡迎光臨散文網(wǎng) 會員登陸 & 注冊

輻射度量學(xué)與蒙特卡諾積分

2023-03-11 17:26 作者:EnemyIncoming  | 我要投稿

背景:

Phong和Blinn-Phong的經(jīng)驗?zāi)P瓦^于經(jīng)驗,沒有嚴(yán)謹(jǐn)?shù)奈锢矶x,所以計算機圖形學(xué)急需要一個能夠符合物理基本規(guī)律的光照知識,這就是輻射度量學(xué)提出的背景(基于物理光學(xué))。

概念:

1. 輻射能Q,單位是焦耳,是能量,表示穿過一個曲面(類似于電磁學(xué)的高斯面)的光能。

2. 輻射通量%5Cphi

%0A%5Cphi%3D%5Cfrac%7BdQ%7D%7Bdt%7D%20

3. 立體角%5Comega(針對球面坐標(biāo)系),類比于平面角(針對極坐標(biāo)系):

%5Ctheta%3D%5Cfrac%7Bl%7D%7Br%7D%20

可以得到立體角的公式:

%20%5Comega%3D%5Cfrac%7BA%7D%7Br%5E2%7D%20

如圖所示,設(shè)球面上的單位面積為dA,那么可以算出:

%20%5Cbegin%7Balign%7D%20dA%26%3D(r%5C%20d%5Ctheta)(r%5Csin%5Ctheta%5C%20d%5Cphi)%3Dr%5E2%5Csin%5Ctheta%20%5C%20d%5Ctheta%20d%5Cphi%5C%5C%20d%5Comega%26%3D%5Cfrac%7BdA%7D%7Br%5E2%7D%3D%5Csin%5Ctheta%5C%20d%5Ctheta%20d%5Cphi%20%5Cend%7Balign%7D%20

4. 輻射強度E%3D%5Cphi%2F4%5Cpi(Radiant Intensity),單位為cd(坎德拉),在光學(xué)中叫做光強。他表示單位立體角的輻射通量:

%20I%3D%5Cfrac%7Bd%5Cphi%7D%7Bd%5Comega%7D%20

5. 輻照度L%3D%5Cfrac%7Bd%5Cphi%20d%5Cphi%7D%7Bd%5Comega%20dA%20%5Ccos%5Ctheta%7D%3D%5Cfrac%7BEd%5Cphi%7D%7Bd%5Comega%20%5Ccos%5Ctheta%7D%3D%5Cfrac%7BdE%7D%7Bd%5Comega%20%5Ccos%5Ctheta%7D%20(Irradiance),表示單位面積的輻射通量:

E%3D%5Cfrac%7Bd%5Cphi%7D%7BdA%7D%20

可以看出,這個量類似于我們在電磁學(xué)學(xué)到的磁通量定義公式:

%5CPhi%3D%5Cint_%5COmega%20%5Cvec%7BE_e%7D%5Ccdot%20d%5Cvec%7BA%7D%20

所以,E可以表示一個球面(或者叫高斯面)上單位面積的接收到光強。那么EI有什么區(qū)別呢?E表示的是一個球面上單位面積的接收到光強,而I表示的是一束光的光強而已。

這個理論可以解釋為什么在Phong模型中,光強按照平方階次衰減。

如圖,離光源最近的輻照度(半徑為1),而距離光源為r的輻照度E'%3D%5Cphi%2F4%5Cpi%20r%5E2%3DE%2Fr%5E2,所以兩者的比值是1%2Fr%5E2,這就證明了衰減是平方衰減。(可以看出來這里和電磁學(xué)非常相似,道理都是相通的)

6. 輻射率L%0A(Radiance),又叫光亮度。這是最常用(最重要)的一個物理量,也是構(gòu)建渲染方程的主要物理量。他每表示單位立體角,每單位垂直面積的輻射通量。很像是Intensity和irradiance的結(jié)合。它同時指定了光的方向與照射到的表面所接受到的亮度。(就是這個單位面積在某個方向上吸收的光能)

%20L(%7B%5Crm%20p%7D%2C%5Comega)%3D%5Cfrac%7Bd%5E2%5Cphi(%7B%5Crm%20p%7D%2C%20%5Comega)%7D%7Bd%5Comega%20dA%20%5Ccos%5Ctheta%7D%20

其中p是入射點。

那么L%0AE的關(guān)系是啥呢?作如下變形:

L%3D%5Cfrac%7Bd%5Cphi%20d%5Cphi%7D%7Bd%5Comega%20dA%20%5Ccos%5Ctheta%7D%3D%5Cfrac%7BEd%5Cphi%7D%7Bd%5Comega%20%5Ccos%5Ctheta%7D%3D%5Cfrac%7BdE%7D%7Bd%5Comega%20%5Ccos%5Ctheta%7D

這樣進一步得到:

%20%5Cbegin%7Balign%7D%20dE(%7B%5Crm%20p%7D%2C%5Comega)%26%3DL_i(%7B%5Crm%20p%7D%2C%5Comega)%5Ccos%5Ctheta%20d%5Comega%20%5C%5C%20E(%7B%5Crm%20p%7D)%26%3D%5Cint_%7BH%5E2%7D%20L_i(%7B%5Crm%20p%7D%2C%5Comega)%5Ccos%5Ctheta%20d%5Comega%20%5Cend%7Balign%7D%20

其中H%5E2是上半球。所以這個式子很好解釋了Radiance是單位面積在某個方向上接受到的光強,而Irradiance是整個Radiance的求和(即所有方向上接收到的光強的總和)。

BRDF函數(shù)與渲染方程:

考慮一個位置為x的單位面積dA,接收到某個方向的光強dE(%5Comega_i),經(jīng)過材質(zhì)的吸收能量后,反射出的光強為dL_r(x%2C%5Comega_r),為了衡量反射出的光強與入射的光強的關(guān)系,類比于反射率的定義,我們定義一個函數(shù)來表示這個關(guān)系:

%20f_r(%5Comega_i%20%5Crightarrow%20%5Comega_r)%3D%5Cfrac%7BdL_r(x%2C%5Comega_r)%7D%7BdE(%5Comega_i)%7D%3D%5Cfrac%7BdL_r(x%2C%5Comega_r)%7D%7BL_i(x%2C%5Comega_i)%5Ccos%5Ctheta_i%20d%5Comega_i%7D%20

這個函數(shù)f,就是BRDF函數(shù)。因為這個函數(shù)的物理意義就是反射率,所以又叫做反射率函數(shù)。

因此,反射出的光強可以解出來是:

%20L_r(x%2C%5Comega_r)%3D%5Cint_%7BH%5E2%7Df_r(%5Comega_i%20%5Crightarrow%20%5Comega_r)%5C%20L_i(x%2C%5Comega_i)%5C%20%5Ccos%5Ctheta_i%5C%20d%5Comega_i

這就是反射方程。

考慮到物體還能自身發(fā)光,所以反射出的光強還包括物體自身的光強:

%20L_r(x%2C%5Comega_r)%3DL_e(x%2C%5Comega_r)%2B%5Cint_%7BH%5E2%7Df_r(%5Comega_i%20%5Crightarrow%20%5Comega_r)%5C%20L_i(x%2C%5Comega_i)%5C%20%5Ccos%5Ctheta_i%5C%20d%5Comega_i

這就是渲染方程。

渲染方程的物理意義

渲染方程可以這樣寫:

L_o(p%2C%5Comega_o)%3DL_e(p%2C%5Comega_o)%2B%5Cint_%7B%5COmega%5E%2B%7Df_r(p%2C%5Comega_i%2C%5Comega_o)%5C%20L_i(p%2C%20%5Comega_i)%5C%20(%5Cvec%7Bn%7D%5Ccdot%20%5Cvec%7B%5Comega_i%7D)%5C%20d%5Comega_i%20

大部份物理意義上面已經(jīng)說過了,因為這里的積分其實本質(zhì)是一種卷積,所以可以用算子來簡化,這里講解一下算子理解,為了簡化方程,我們可以這樣寫一個方程大致表示渲染方程:

L%3DE%2BKL

這是一種遞歸的寫法,$L$為反射光強,$E$為光源光強,而$K$是反射算子(一種矩陣)。由線性代數(shù)的知識有:

%20%5Cbegin%7Balign%7D%20L%26%3DE%2BKL%5C%5C%20L-KL%26%3DE%5C%5C%20IL-KL%26%3DE%5C%5C%20(I-K)L%26%3DE%5C%5C%20L%26%3D(I-K)%5E%7B-1%7DE%20%5Cend%7Balign%7D%20

由泰勒展開,有:

%20(I-K)%5E%7B-1%7D%3DI%2BK%2BK%5E2%2BK%5E3%2B...%20

所以反射光強:

L%3DE%2BKE%2BK%5E2E%2BK%5E3E%2B...%3D%5Csum_%7Bi%3D0%7DK%5EiE%20

所以可以理解為,反射光強等于光源光強,加上光源經(jīng)一次反射后的光強,加上光源經(jīng)兩次反射后的光強,以此類推。這就是真實的全局光照了。

蒙特卡諾積分

現(xiàn)在有了渲染方程,但如果直接數(shù)值方法計算定積分,效率太低了。所以引入概率學(xué)的蒙特卡諾方法,類似于幾何概型。

蒙特卡諾方法會將:

%20%5Cint_a%5Ebf(x)dx%20

以概率的角度,求出他的近似值。

考慮一個函數(shù)圖像:

選定一個隨機變量X_i%5Csim%20P(x),遵循某種分布P(x)。從這個分布隨機挑選出一個值x_i%5Cin%20X_i%5Cin%5Ba%2Cb%5D,因為這個集合關(guān)系,所以就會導(dǎo)致P(x_i)的范圍為:

1%2FP(x_i)%5Cin%5Ba%2C%20b%5D%20

所以1%2FP(x_i)決定了矩形的寬度,可以用來計算矩形面積。

然后取f(x_i)作為矩形的高度,然后計算矩形面積:

%20S%3Df(x_i)%2FP(x_i)%20

重復(fù)$N$次,最后取算術(shù)平均數(shù),就是這個積分的近似值了(是不是很像幾何概型):

F_N%3D%5Cfrac%7B1%7D%7BN%7D%5Csum_%7Bi%3D1%7D%5EN%5Cfrac%7Bf(x_i)%7D%7BP(x_i)%7D%20

于是:

%20%5Cint_b%5Eaf(x)dx%5Capprox%20%5Cfrac%7B1%7D%7BN%7D%5Csum_%7Bi%3D1%7D%5EN%5Cfrac%7Bf(x_i)%7D%7BP(x_i)%7D%20

如果我們采用均勻分布,那么X_i%5Csim%20U(a%2Cb),就有P(x_i)%5Cequiv1%2F(b-a),這個時候每次的獲取的矩形面積的寬度都是(b-a),于是每次隨機計算出的矩形面積為:

%20S%3Df(x_i)(b-a)%20

于是:

%20F_%7BN%7D%3D%5Cfrac%7Bb-a%7D%7BN%7D%5Csum_%7Bi%3D1%7D%5ENf(x_i)%20

證明蒙特卡諾積分的正確性

因為蒙特卡諾是概率模型,我們計算數(shù)學(xué)期望,如果值就是積分本身,那么這個模型就是正確的:

%0AE%5Cleft%5B%5Cfrac%7Bf(x_i)%7D%7BP(x_i)%7D%5Cright%5D%3D%5Cint%20%5Cfrac%7Bf(x_i)%7D%7BP(x_i)%7DP(x_i)dx%3D%5Cint%20f(x_i)dx%20

證明完畢


輻射度量學(xué)與蒙特卡諾積分的評論 (共 條)

分享到微博請遵守國家法律
安化县| 汤阴县| 敦煌市| 鄱阳县| 祁门县| 江都市| 克山县| 勃利县| 和平区| 灵武市| 古丈县| 津市市| 姜堰市| 青川县| 鹤壁市| 建平县| 罗田县| 定州市| 英山县| 霍州市| 德化县| 阿勒泰市| 衡山县| 准格尔旗| 金沙县| 台安县| 昌都县| 图们市| 旬邑县| 和龙市| 宁安市| 曲沃县| 丰城市| 安岳县| 大厂| 林西县| 崇仁县| 射洪县| 景宁| 绍兴市| 蚌埠市|