第五章 光学模型及其算法实现

1、121第五章 光学模型及其算法实现第五章 光学模型及其算法实现一、复习要求1简单光反射模型2增量式光反射模型3局部光反射模型4光源模型5简单光透射模型6光线跟踪显示技术二、内容提要1简单光反射模型（1）基本光学原理照度定律普通物理学中的照度定律(Lambert余弦定律)：I=KfIL cos,q Î 0, 式中，反射常数Kf与物体表面性质有关，也描述物体的颜色。注意：按一般规定，入射光L是从表面上一点指向光源的矢量。 材质分配事实上，一个物体的颜色就是它所反射出的光的颜色，取决于光源的颜色和该物体对光的反射性。例如，将一个阳光下的红球放到只有黄灯照明的室内，它就变成黑的了，因为那里没

2、有任何红色光可被反射，而所有黄色光都被吸收了。在使用光照时，原有的“绘图颜色”概念已不能适用，而采用“材质”一词。定义：材质(material)被定义为一个物体对环境光、漫反射光、镜面反射光的反射性。它们分别以一个对应的RGB值表示，称为材质的Ambient, Diffuse, Specular分量（即光学定律中的反射系数Ka, Kd, Ks）。材质还可以包括另一种辐射性，用于描述自身发光的物体，例如汽车尾灯或夜光表。通常，灯具的表面也被看成是一个自发光体。 折射和透射Snell正弦定律（或称折射定律）属于几何光学原理，用于确定两个物体间的入射角与折射角的关系：式中，q1表示在物体1表面处的入

3、射角，而q2表示在物体2内部的折射角；h1和h2分别是这两个物体的折射率。如图5-1所示，除了从同侧光源射过来的反射光外，观察者还将会看到从另一侧光源穿过物体后射出的透射光：图5-1 光的折射和透射（2） 简单光反射模型(Phong模型)的导出图5-2 光的反射简单光反射模型只模拟物体表面对光的反射作用，并不考虑物体表面的透射和散射作用。在简单光反射模型中一个点光源照射到物体表面一点，再反射出来的光，可分为三部分：环境光（泛光）、漫反射光镜面反射光。在Phong模型中，物体表面的光照效果是环境反射光、光源的漫反射光和镜面反射光的合成效果。· 环境反射光(Ambient Light)：

4、从背景物体散射出来的光。可用下式计算：I1=KaIa 其中Ka为环境光系数· 光源的漫反射光(Diffuse Light)：从一个粗糙的、无光泽的表面反射出，它与入射角有关（即符合Lamber定律）I2=KdILcos 其中Kd为粗糙度· 镜面反射光(Specular Reflection)：从一个光滑的、明亮的表面反射出（或称单反射光）：I3=KsILcosn 其中，Ks为单向反射率；n控制反射光束的锥度，其值越大，光束越窄，表示该表面越明亮（如磨光金属面 大理石）（较小可表示苹果、木质等） 当= 0（即视线位于反射方向上）时，将呈现一个高光点(hotspot)（如白光下

5、的铜球上的亮白点）根据Phong模型的定义，在多边形所有的点处按下式计算光强： I=I1+I2+I3 = KaIa + KdILcos+ Kscosn注意：在具体计算中，反射系数Kd、Ks和Ka也可按RGB分解成三个分量：(Kdr, Kdg, Kdb)、(Ksr, Ksg, Ksb)和(Kar, Kag, Kab)，各分量范围均为0.0 1.0。2增量式光反射模型（1） 双线性光强插值法(Gourand Shading)在图像空间中，Gouraud光照模型（双线性光强插值法）计算各像素的光强如下：图5-3 求法矢的平均值(1) 只在多边形的顶点处，按下式计算光强：I = KaIa + KdIL

6、cos(只考虑环境光和漫反射光) 而顶点的法矢等于各邻面法矢的平均值：(2) 对于多边形边上和内部的各点(像素)，用顶点明暗度的线性插值计算出：图5-4 双线性亮度插值法令 DIAB = ，DIAD = 故 IS = IA + DIAB·Dy，IT = IA + DIAD·Dy对于多边形内部的各点（像素），用两边端点光强的线性插值法计算出：令 DIST = 故 IP = IS + DIST·Dx另一种计算多边形边上和内部的各点光强的改进方法是增量法，它利用扫描线的相关性来简化计算如下：对于边上的点：IS, y+1 = IS, y + DIAB，IT, y+1 =

7、IT, y + DIAD对于内部的点：Ix+1, P = Ix, P + DIST，Gouraud算法优缺点：算法简单，计算量小。但不能再现高光。适用粗糙表面（2） 双线性法向插值法(Phong Shading) 在图像空间中，Phong光照模型（双线性法向插值法）计算各像素的度如下：顶点的法矢等于各邻面法矢的平均值，即对于多边形上和内部的各点(像素)，用顶点法矢的线性插值。计算公式如下：图5-5 双线性法矢插值法令DNS = ，DNT = 故NS, y+1 = NS, y + DNS，NT, y+1 = NT, y + DNT令DNP = 故 DNx+1, P = DNx, P + DNP注

8、意：这里，法矢Np就是实际曲面上该点的切平面法矢的近似值Phong算法优缺点：真实感强，但计算量较大，适用光滑表面目前大多数渲染软件采用Gouraud和Phong光照模型作为Quick(快速)和Full(完全) Rendering，能满足一般真实感显示要求。Gouraud和Phong光照模型两者比较如下：Gouraud ShadingPhong Shading双线性光强插值法双线性法向插值法适用于漫反射光（粗糙表面）适用于镜面反射光（特别是有高光的金属表面）马赫带效应较重马赫带效应较轻生成多面体真实感图象效果差高光域准确速度较快（算法简单、计算量小）速度较慢（计算量较大）然而，这类简单光照模型

9、共同存在下列的问题：¨ 由于光强函数的一阶导数不连续，线性插值法均有马赫带效应（光照效果在数值上不连续）。¨ 透明效果需要另作处理才行，而且不能产生阴影。注：实用Phong光照模型的计算公式要考虑随距离的线性或平方反比衰减（3） 加速算法 在简单的光线跟踪算法中，每条射线都要和所有物体求交，然后再对所得的全部交点进行排序，才能确定可见点。因此对于每条射线，计算复杂度都是O(n)以上。当物体个数较少时，计算量还可以接受。但是当环境较复杂，物体个数较多时(例如成干个原子组成的分子模型或成万个多边形构成的建筑物模型)，简单的光线跟踪算法的处理速度就无法接受了。然而，对于每条射线，

10、它实际上只与少数几个物体有交，与绝大多数物体根本不相交。也就是说，对于复杂环境，大多数求交计算都是无效运算。我们可以通过一些途经来减少这种无效运算。这就是加速算法所要考虑的主要问题。一般说来，加速算法是通过预先把物体按空间位置适当地组织起来，以便在绘制时缩小搜索范围，减少排序时间来提高效率的。常用的组织方法有包围体树、自适应八叉树和三维立方体阵列。3局部光反射模型局部光反射是由光源直接照射到射入点，经漫反射和镜面反射到视点而形成。4光源模型在计算机图形学中，光源模型的三个基本因素是几何形状、光强分布和光谱分布。¨ 有向光源(directed lights)：方向保持不变，即为平行光（

11、发射出半圆柱面光）；光源强度不作衰减。如太阳光、远光源。¨ 点光源(point lights)：方向是从光源位置到面片中心，发射出圆球面光；光强分布是各向同性的，要考虑随距离的线性或平方反比衰减(falloff)；可以产生阴影。¨ 聚光灯(spotlights)：发射出圆锥形光束；要考虑随距离的线性或平方反比衰减；可以产生软阴影区域。从几何形状上讲，一般电灯（钨丝白炽灯）理想化成点光源，而且光强分布是各向同性的；一般日光灯（长管灯）可理想化成线光源；经过特殊设计的吸顶灯可理想化成面光源。5简单光透射模型（1） 透明效果的模拟方法 透射计算公式最简单的透明算法不考虑折射的影响

12、，即假定光线穿过透明体表面时不改变方向，此时则有：（如图12.11C所示）ITR = (1 - KT)Ireflect + KTItransmit， 其中KT Î 0, 1KT称为物体的透明度(Transparency)，可取值范围从0（不透明）至1.0（全透明）。光线照射到透明体表面的两个交点A和B处的光强分别是反射光强Ireflect和透射光强Itransmit。图5-6 光的透射显然，在全透明场合下（即KT = 1），只有光的透射效果：ITR = KTItransmit这也称为规则透射。注意：1. 在颜色四元组RGBA中，ALPHA = 1 - KT，称作不透明度。 2. 对不

13、透明物体，将会在地面上投射出阴影(Shadow)。 漫透射、理想透明和半透明一个点光源照射在透明体表面上，将可能产生如下三种光透射作用：¨ 对于表面粗糙的透明物体上，将形成漫透射。透射出来的光将朝着各方向均匀地散射，此时透射光强It与入射光强IL及入射角的余弦cosq成正比，即：It = KtILcosq式中，Kt称为透射系数，可取0 1.0。¨ 进入理想的透明物体的透射光只沿折射方向传播，而在其他方向都看不见透射光。¨ 进入一般较透明（或半透明）物体的透射光也往其他方向传播一部分透射光，并且某方向透射光的强度随着偏离透射方向（即下述b角的增大）而明显地衰减。It

14、 = KtILcosntb式中，IL是光源的强度或入射光强；Kt称为透射系数，可取0 1.0。幂值nt用于控制透射光束的锥度，其范围1 2000；其值越大，光束越窄。b为视线V与透射方向T之间的夹角。（2） Wittded光透射模型Whitted光透射模型是1980年Whitted提出的，并运用在光线跟踪算法的实现中。它考虑了环境光、漫反射光和镜面反射光（即Phong模型）以及从环境反射或透视到视线的光。具体计算如下：Ip = KaIa + + KsIs + KtIt = Ia·Ka + Kd + Ks + KsIs + KtIt式中，第一部分包含了Phong模型中的所有参数；第二部

15、分模拟了场景光投射在光滑表面上产生的理想镜面反射作用（即a = 0），Ks为镜面反射率；第三部分模拟了场景光通过透明体表面产生的规则透射作用（即不考虑折射时的全透明），Kt为透射系数。Is和It分别表示场景中其他物体反射和透射过来的光在该物体表面上的镜面映像和透射映像，如图所示。图5-7 透明体表面上的镜面映像和透射映像（3） Hall光透射模型Hall整体光照模型要比Whitted模型更为完善，还考虑沿着反射方向和折射（透射）方向的散射现象，以模拟“透射高光”。在Hall整体光照模型中，假定光源和视点都在无穷远处，使光线和视线方向简化为常向量。计算公式如下： Ip = I局部 + I全局(1

16、)I局部 = KaIa + = Ia·Ka + Kd+ Ks+ Kt即包括环境光、漫反射光、镜面反射光和透射光。(2)I全局 = KRIR + KTIT即由其他物体反射或透射在透明体表面上产生的沿着视线方向的光。6光线跟踪显示技术（1） 基本光线跟踪算法Ray-Tracing算法的基本思想如下：在物体空间中，假定视点取在Z轴上，XOY平面取作为投影屏幕(其与显示器屏幕的像素布局一一对应)，通过跟踪多个光源对各像素的贡献而计算出它们的色彩明暗度。在光线遇到某一物体（不是背景）时，应继续跟踪光线，而不管光线的强度或深度如何。图5-8 光线跟踪法光线跟踪算法包括下列的三个步骤：步骤一、 从

17、视点穿过某像素P发出一条射线R，它逆着光线方向跟踪所有光源在这个可见点（即离视点最近的交点）上产生的色彩明暗度：IA = I局部 + I全局其中I局部和I全局的具体计算按步骤二。步骤二、 首先计算出可见点处的曲面法矢（并存储它，以供备用），然后查找表面数据表（其含有表面的颜色属性、粗糙度、反射率、透明度等）进行计算：（Hall整体光照模型）(1)I局部 = KaIa + IL(Kdcosq + Kscosnsa + Ktcosntb)即包括环境光、漫反射光、镜面反射光和透射光。(2)I全局 = KRIR + KTIT即由其他物体反射或折射到视线的反射或透视方向的光。步骤三、判断可见点是否处在阴

18、影(shadow)中：从该点向光源引射线。若射线与某个不透明的物体相交，则该点在阴影中（此时只取环境光），而返回到步骤一。光线跟踪法的特点：a. 同时处理消隐和渲染，并有透明效果和阴影生成。b. Ray-Tracing本质上是一种递归算法（层次等于物体个数），效率极高算法加速方法：事实上每条射线只与少数几个物体相交。若预先将物体按空间位置适当地组织起来（如采用包围体），则可缩小搜索范围，避免不必要的求交运算。树：仅当射线与根结点有交时，才进行它与子结点求交优缺点：光线跟踪法是目前最常用的高级渲染技术之一，但没有考虑物体表面之间的漫反射（如辐射度计算方法）。在光线跟踪算法中，跟踪光线的终止条件有

19、以下四个：1) 被跟踪的光线与光源相交；2) 被跟踪的光线遇到背景；3) 被跟踪的光线对某交点处的光强贡献趋于零；4) 光线跟踪的深度已经很深了。（2） 光线与物体求交与坐标系平面求交光线与坐标系平面的交点可以由下面向量方程求出：s+td=xiI+yiJ 由于s=xsI+ysJ+zsK，d=xdI+ydJ+zdK，可推出：下面两种情况，光线与平面不相交：zd=0时，光线平行于xy平面zs=0时，光线起始点在xy平面上zd0，且zs0时如果t>0，光线与平面相交如果t<0，光线的反向延长线与平面相交与任意平面求交右图中,n=xnI+ynJ+znK，为给定平面的法线向量，P0(x0,y

20、0,z0)为平面上一点，则平面向量方程为：xn(x-x0)+yn(y-y0)+zn(z-z0)=0令r(t)=s+td与该平面相交，又s=xsI+ysJ+zsK，d=xdI+ydJ+zdK，则xn(xs+txd-x0)+yn(ys+tyd -y0)+zn(zs+tzd -z0)=0可推出下式：另，则当=0时，光线平行于该平面，即不相交。如果t=0，光线起始点在平面上如果t>0，光线与平面相交如果t<0，光线的反向延长线与平面相交与球面求交右图为一半径为R，中心点为(xc,yc,zc) 的球体，(x,y,z)是球面上的任意一点。设向量c=xcI+ycJ+zcK，p=xI+yJ+zK则

21、球体的向量方程为：|p-c|=R令r(t)=s+td与该平面相交，则|s+td-c|2=R2展开得：|d|2t2+2|s-c|d|t+|s-c|2=R2令A=|d|2，B=|s-c|d|，C=|s-c|2-R2，则At2+2Bt+C=0（3） 光线跟踪中的简单阴影定义5-1：阴影(shadow)是指景物中那些没有被光源直接照射到的暗区。图5-10 阴影生成阴影生成从原理上来说非常简单：只有那些从视点可见的、而从光源不可见的面（多边形）才会落在阴影内。换句话说，若将视点放在光源处，那末它看不到所该光源产生的阴影。因此，阴影生成的一般性算法如下：步骤一、将光源位置作为视点，利用消隐算法求出它的不可

22、见面（包括背面和被其他面遮挡的面）和可见面。这样，将所有面分为两大类。图5-11 半阴影区步骤二、在实际的视点处，对上述面进行消隐，并选用一种光照模型双重可见面的明暗度，而对上述不可见面的可见部分采用环境光，形成阴影。步骤三、将不可见面的可见部分投射到地面上，即表示影子(shadow)。注意：仅当存在某一投影面（如地面、墙面）时，物体才有可能在该平面上投射出其影子。定义5-2：阴影可分为全阴影(whole-shadow)和半阴影(half-shadow)两种。全阴影就是物体表面上那些没有被所有光源直接照射的部分，而半阴影是指物体表面上那些只被一部分光源直接照射的部分。定义5-3：全阴影加上在它

23、周围的半阴影而组成一个软阴影(soft-shadow)区域，使影子边缘有较缓慢的过渡。显然，单个点光源只能形成全阴影，而半阴影则是由多个点光源或线光源、面光源产生的，如图5-11所示。由于半阴影的计算量较大，许多场合下（除聚光灯外）我们只考虑由点光源产生的全阴影。因为光源在物体表面上产生的阴影区域均为它们的不可见面。若取光源作为视点，那么任何隐藏面消除算法都能用于阴影计算。在实际应用中，还需根据阴影计算的特点考虑如何减少其时间和空间耗费。三、例题分析例1设n=I+J+2K为通过点P0=(1,1,0)的平面的法线矢量。试判断s=-2I+J+2K，d=I-K的光线是否与平面相交。解答：因为 故而可

24、以引入： >0，光线与平面相交所以交点是（xs+xd，ys+yd，zs+zd）=（-1，1，1）四、习 题1写出简单光反射模型近似公式，并说明其适用范围及能产生的光照效果。2写出线光源的光强公式及其积分算法。 3试描述Witted光透射反射模型和Hall光透射模型。4简单叙述光线跟踪算法。5试描述光线与几种常见物体面的求交算法。五、解题指导与习题解答1答案：简单光反射模型只模拟物体表面对光的反射作用，并不考虑物体表面的透射和散射作用。在简单光反射模型中一个点光源照射到物体表面一点，再反射出来的光，可分为三部分：环境光（泛光）、漫反射光镜面反射光。根据Phong模型的定义，在多边形所有的点处按下式计算光强： I=I1+I2+I3 = KaIa + KdILcos+ Kscosn在Phong模型中，物体表面的光照效果是环境反射光、光源的漫反射光和镜面反射光的合成效果。2答案:线光源的光强公式如下：其中 表示漫反射系数， 表示漫反射率， 是来自光源的光亮度， 是 朝辐射方向的投影角。 是被照射表面的法向， 是单位光线向量。在一般情况下，上述公式只能用数值算法做近似的计算。其实质是用多个点光源来近似线光源。3答案:Whitted光照模型是一种全局光照模型，它模拟了由环境引起的镜面反射分量和规则透射分量。在这个模型中，假定被照明物体的表面非常