计算机图形学(一)——辐照度学概述

辐照度(irradiance)定义为电磁辐射入射于曲面时每单位面积的功率,直观理解为能量在入射点处的量化值。在计算机图形学中,就是要求解场景中任意位置的辐照度大小。如果求得辐照度的大小,剩下的事情就是如何利用所求解得到的辐照度渲染出一张静帧图像,而不同的求解方法,对应着不同的渲染算法,每一个算法都有其独特的风格。本部分就是要理清计算机图形学中最基本的辐照度物理量之间的定义与关系。

 

1.1 立体角

 

立体角的概念其实是平面几何中角的概念的扩展。数学上定义平面角度是在单位圆上,两射线张角的弧长即为张出的角度,单位为radian;则类比平面几何中对角度的定义,我们可以去如下定义立体角的概念:在单位球体上,球心处一个锥形包含的球表面积定义为这个锥形的立体角(如图1),单位为stereo radian,简写为sr。计算机图形学中常用粗体的来表示一个立体角。需要注意的是,立体角的形状并不仅限于像图1中那样的锥形,还可以是像那样的不规则形状,我们仅用单位圆内该立体角内部包含的球表面积来度量它的角度,由于单位球面的表面积为4π,所以一个立体角的取值范围也在[0,4π]之间。有了上述的介绍,我们就可以继续定义立体角的微分了。即当一个立体角所包含的单位球表面积趋于0时,我们称其为立体角的微分,写为。为表述简便,后续有关“方向”的概念,如果没有特殊说明,默认都是在说指向这个方向的一个立体角。

 

1.2 辐照度中的基本量及其关系

辐照度学中,定义了如下表1的常见基本量,详细请参见相关的教材。

上述物理量中,较为常用的几个物理量为辐射通量Φe、辐照度Ee和辐射率Le。为表述简便,后续的内容都会省略角标,即简写为ΦEL。现定义它们之间满足的基本关系。

(1)E=dΦdA
式中之所以用单位面积dA上入射的功率Φ而不是能量Q来定义辐照度E是因为前者考虑了时间因素,把时间度量统一到单位时间。如果按照后者来定义的话,会造成物理量E与时间相关,这样做并不是一个明智的选择。

 


定义垂直入射光方向的辐照度
(2)E=dΦdA
根据图2,上式中的AA有如下的关系
(3)A=Acosθ
式(3)中,θ是入射方向与法线的夹角,注意这里的θ取值范围为0θπ2,使得cosθ总是非负数。为简便起见,我们规定入射方向是沿着入射点朝向入射源的,直观上这与“入射”的含义相悖,但是并不影响后续的分析。这样一来,所有有关“方向”的概念都是沿物体表面入射点朝向表面外,反而会使很多式子的形式得到简化。我们进而又能很简单地得到EE之间的关系
(4)E=Ecosθ

有了式(1),我们可以进一步定义辐射率
(5)L=dEdω
所以
(6)L=dΦdAdω=dΦdAdωcosθ
由(4)、(5)可得
(7)dE=Lcosθdω
之所以定义辐射率的概念,是因为我们经常要衡量单位球表面上任意一点的辐照度,但是点的面积是不存在的,所以要用辐照度微分与立体角微分之比的形式来定义辐射率。式(6)、(7)是很重要的公式,在之后对表面处的辐射率估计的推导中会经常用到。

 

1.3 符号约定

由于后续的推导过程中运用了大量的数学概念,我们在正式介绍这些内容之前,首先要约定符号的写法以及各种物理量的表达方式,否则很容易造成混淆。我们规定:

1) 立体角用粗体小写的ω来表示,其微分形式为dω,这在1.1节中已经陈述。
2) 以物体表面任意点x处,沿表面局部法线建立的半球面区域写为Ωx,其立体角恒为2π
3) 双向反射分布函数(下文会详细介绍)用fr来表示,且其第一个参数x为参考点的位置,第二个参数为入射方向,第三个参数为出射方向。
4) 表1中的辐照度(irradiance)和辐射度(radiosity)的符号应该予以统一,以符号Ei表示辐照度,以符号Eo表示辐射度。类似地,规定入射辐射率符号为Li,出射辐射率符号为Lo


1.4 双向反射分布函数(BRDF)

1.4.1 BRDF的定义

双向反射分布函数(Bidirectional Reflectance Distribution Function, BRDF),是用于描述物体表面性质的重要函数之一,它被定义为沿出射光方向ωo出射的辐射率Lo与沿入射光方向ωi入射的辐照度Ei的比值,单位为sr1。由上述表述,则在物体表面x点处的BRDF的定义式即为

(8)fr(x,ωi,ωo)=dLo(x,ωo)dEi(x,ωi)=dLo(x,ωo)dEi,(x,ωi)cosθi=dLo(x,ωo)Li(x,ωi)cosθidωi

 

如图3所示,为形式简便,将入射方向规定为由参考点指向源点,这样所有有关方向的物理量均统一为沿表面指向外界。n为物体局部法线,一束入射方向为ωi的光线击中物体表面x点处,其能量会在以x点处形成的半球面沿一定分布特征辐射出去。那么以ωi为入射方向,以ωo为我们要考察的出射方向,沿ωo方向出射的能量dLo与入射光线在x点处对物体表面的辐照度贡献dEi之比,即被定义为在x处,沿入射方向ωi和出射方向ωo的双向反射系数fr(x,ωi,ωo)

这里请注意式(8)中“d”的含义与数学上严格意义上的微分是有区别的,Lo(x,ωo)是指在物体表面x点处,半球面ωx区域内所有入射光的能量总和作为输入时,沿ωo方向出射的辐射率;而dLo(x,ωo)则特指以方向ωi处入射光的能量作为输入时,沿ωo方向出射的辐射率。dEi(x,ωi)特指方向ωi处的辐照度。dEi(x,ωi)Ei(x,ωi)的含义是相同的。

针对BRDF的定义式,有如下两个问题需要解释:

1) 为什么要用微分形式去定义?
2) 为什么不定义成dLodLi或者dEodEi,使得它变成一个纯粹的系数而不是一个具有量纲sr1这样奇怪的一个物理量?
接下来我们就来回答上面两个问题。之所以要用微分形式去定义,在于后面有关BRDF的所有计算几乎都涉及到积分,所以原始的BRDF定义式将BRDF定义为微分形式会更有利于后续对积分式的计算,而且形式美观简洁,但之后进行公式推导时,仍需要明确分子的Lo是特指ωi方向上对出射辐射率的贡献,以后为强调这个事实,也常用Lo(x,ωo)|ωi表示dLo(x,ωo);而针对BRDF的量纲问题,则可以这样解释:BRDF的定义是根据对实际现象的物理建模而得出的,这个量也最好应该能够通过实验仪器测得,那么定义式的形式,实际上就决定了对测量仪器的要求。直接测量物体表面某点沿某入射方向的入射光辐照率Li是有困难的,而想测量沿着某一方向的出射辐射率Lo则很简单,只需要将一个管状测量仪器对准入射点处,就能测得沿着管状测量仪出射的辐射率Lo。虽然测量入射光辐照率由困难,但是测量入射光的辐照度Ei则很简单,因为辐照度没有方向的概念,所以研究者仅需要制作一个半球状的光线接收器即可。

实际计算的时候,通常是以物体表面参考点x为原点,以物体表面x点处指向物体表面外侧的局部法线nz轴建立的求坐标系来进行积分的求解,所以式(8)中的BRDF定义式也可以写为fr(x,θi,φi,θo,φo),如果物体表面的反射特性与参考点x的位置无关(通常情况下),那么我们可以在fr中省略变量x,变为fr(θi,φi,θo,φo),图3中θiφiωi的球坐标系参数表示,θoφoωo的球坐标系参数表示。
很多实验室测得的BRDF数据都是一个巨大的四维表格(例如MERL BRDF数据库,每种材质表面测得的BRDF数据有34MB,具体详见http://www.merl.com/brdf/),用户在实际使用这张表格的时候,就是依据θiφiθoφo四个独立变量来查表得到对应的BRDF值,进而完成渲染的。


1.4.2 BRDF的性质

BRDF具有以下三条重要性质:
1) 非负性:对于物体表面上任意点x以及半球面内所有ωiωo,BRDF函数值均非负。这一性质很好理解,即便是绝对黑体,它的BRDF也是恒为0而不可能出现负值。
(9)fr(x,ωi,ωo)0
2) 互易性:由基本光路可逆原理,若交换入射方向和出射方向,BRDF函数值不变。此性质可将原始BRDF数据量降为原来的一半。
(10)fr(x,ωi,ωo)=fr(x,ωo,ωi)
3) 能量守恒性:沿半球面向外辐射的辐射率总和应小于等于入射光的辐射率总和。即对于ωi
(11)Ωxfr(x,ωi,ωo)cosθodωo1
BRDF的能量守恒性可以用来鉴别给定的BRDF函数是否合理,同时也可以给未归一化的BRDF进行归一化操作。下面我们就来简单证明一下上述的能量守恒性。

证明:假设有一束辐照度Ei(x,ωi)
的光线,沿ωi方向入射并击中物体表面x点处。设x点处朝物体表面外的半球面区域为Ωx,由能量守恒定律,在排除该物体表面为光源的前提下,沿半球面Ωx向外辐射的所有辐射度之和应小于等于入射光的辐照度。即满足
Ωxd[Eo(x,ωo)|ωi]Ei(x,ωi)
上式中的Eo(x,ωo)|ωi是特指沿方向ωi入射的光线所造成的辐射度贡献。
由式(7)可得
Ωx[Lo(x,ωo)|ωi]cosθodωoEi(x,ωi)
根据式(8),又有
Ωx[fr(x,ωi,ωo)Ei(x,ωi)]cosθodωoEi(x,ωi)

Ωxfr(x,ωi,ωo)cosθodωo1
对于ωi均成立。

 


 

posted on   烈日行者  阅读(7775)  评论(4编辑  收藏  举报

编辑推荐:
· 聊一聊 C#异步 任务延续的三种底层玩法
· 敏捷开发:如何高效开每日站会
· 为什么 .NET8线程池 容易引发线程饥饿
· golang自带的死锁检测并非银弹
· 如何做好软件架构师
阅读排行:
· 欧阳的2024年终总结,迷茫,重生与失业
· 聊一聊 C#异步 任务延续的三种底层玩法
· 上位机能不能替代PLC呢?
· 2024年终总结:5000 Star,10w 下载量,这是我交出的开源答卷
· .NET Core:架构、特性和优势详解

导航

< 2025年1月 >
29 30 31 1 2 3 4
5 6 7 8 9 10 11
12 13 14 15 16 17 18
19 20 21 22 23 24 25
26 27 28 29 30 31 1
2 3 4 5 6 7 8
点击右上角即可分享
微信分享提示