CN102306401B - 一种包含模糊反射效果的3d虚拟场景的左右眼立体画面绘制方法 - Google Patents

Abstract

本发明公开一种包含模糊反射效果的3D虚拟场景的左右眼立体画面绘制方法，属于真实感三维虚拟场景绘制技术领域。传统的三维虚拟场景的左右眼立体画面绘制方法采用分别独立绘制左右眼虚拟相机的图像画面的方式实现立体画面的绘制。本发明方法在双核CPU上分别同时执行左右眼虚拟相机的图像画面的绘制程序，左眼虚拟相机的图像画面从第一行像素开始绘制，右眼虚拟相机的图像画面从最后一行像素开始绘制，保证了左右眼虚拟相机的图像画面绘制过程中的入射辐射亮度计算结果的可重用性；通过基于梯度的插值方法，在绘制过程中重用了入射辐射亮度计算结果。利用本发明方法绘制三维虚拟场景立体画面比传统方法快20％～35％。

Description

一种包含模糊反射效果的3D虚拟场景的左右眼立体画面绘制方法

技术领域

本发明属于真实感3D虚拟场景绘制技术领域，涉及一种包含模糊反射效果的3D虚拟场景的左右眼立体画面绘制方法。

背景技术

三维动画电影作为近年来新兴的计算机艺术，发展势头非常迅猛，已经在许多行业得到了广泛的应用。为了产生更强的视觉冲击力，目前很多三维动画电影都具有立体视觉效果。立体视觉效果的基本原理是，观众在观看动画电影时，左右眼能分别互不干扰地看到左右眼虚拟相机拍摄的图像画面，再通过人脑的合成，最终产生画面的入屏或者出屏效果。

为了提高三维动画电影画面的真实感，往往要求在绘制的图像中加入全局光照效果。实现全局光照效果的关键是，正确地模拟3D虚拟场景中的不同物体表面对光的反射过程。在三维图形绘制中，常见的光照反射类型有：镜面反射(Specular Reflection)、漫反射(DiffuseReflection)和模糊反射(Glossy Reflection)。镜面反射和漫反射实际上可以看成是模糊反射的特例。通常使用双向反射分布函数(BRDF)来建模物体表面对光的反射特性。2004年Pascal Gautron等人提出用半球谐和基函数来表示半球域函数(Hemispherical Functions)(见2004年欧洲图形学学会出版的Eurographics Symposium on Rendering会议论文集中的论文《A Novel Hemispherical Basis for Accurate and Efficient Rendering》，作者Pascal Gautron，Jaroslav Krivanek，Sumanta Pattanaik，Kadi Bouatouch)。由于3D虚拟场景中的物体表面的BRDF和入射辐射亮度等物理量实际上都属于半球域函数，因此可以用半球谐和基函数来表示。

传统的3D虚拟场景立体画面绘制方法采用分别独立绘制左右眼虚拟相机的图像画面的方式来实现立体画面的绘制。模糊反射是真实感3D虚拟场景中的一种常见光照反射类型。然而，由于模糊反射的出射光照值与观察角密切相关，而3D虚拟场景中的同一点相对于左右眼虚拟相机的观察角通常并不相等。因此，3D虚拟场景中的物体表面上的同一点在左右眼虚拟相机的图像画面的像素中所表现出的光照值通常并不完全相等，这使得左眼虚拟相机的图像画面的模糊反射光照计算结果难以直接在右眼虚拟相机的图像画面中复用。本发明利用半球谐和基函数来表示3D虚拟场景中的物体表面的BRDF和入射辐射亮度，通过入射辐射亮度缓存实现已有计算结果的重用，以便提高包含模糊反射效果的3D虚拟场景的左右眼立体画面绘制的速度。

发明内容

本发明的目的在于提供一种包含模糊反射效果的3D虚拟场景的左右眼立体画面绘制方法。本发明方法需要的设备包括：带有双核CPU和共享内存的计算机***。如图1所示，当左眼虚拟相机(101)和右眼虚拟相机(102)能同时拍摄到3D虚拟场景中的物体(105)的表面上的点p₁、点p₂和点p₃时，点p₁、点p₂和点p₃的光照值将分别保存在左眼虚拟相机(101)的像平面(103)和右眼虚拟相机(102)的像平面(104)上的像素中。首先将3D虚拟场景中的漫反射和模糊反射表面的双向反射分布函数(Bi-directional ReflectanceDistribution Function，缩写为BRDF)展开成低阶半球谐和基函数(Hemispherical HarmonicsBasis Function，缩写为HSHBF)的累加和的形式，把各阶HSHBF的系数保存在共享内存中并与该表面的BRDF相关联。然后在双核CPU上同时绘制左眼虚拟相机(101)的图像画面和右眼虚拟相机(102)的图像画面，将漫反射和模糊反射表面点的入射辐射亮度展开成低阶HSHBF的累加和的形式，并保存在入射辐射亮度缓存中，通过基于梯度的入射辐射亮度插值计算，实现入射辐射亮度计算结果的重用，从而提高3D虚拟场景的左右眼立体画面的绘制速度。

本发明方法提供一种入射辐射亮度记录项数据结构，其包括位置分量、局部坐标系分量、HSHBF展开系数分量、HSHBF展开系数平移梯度分量等成员变量。

本发明方法将入射辐射亮度记录项保存在位于共享内存中的入射辐射亮度缓存中；为了便于根据空间位置检索入射辐射亮度缓存中的入射辐射亮度记录项和向入射辐射亮度缓存中添加新的入射辐射亮度记录项，入射辐射亮度缓存使用三维空间均匀网格数据结构；根据给定的三维空间位置，可以直接检索出对应网格中的所有入射辐射亮度记录项，以及向网格中添加新的入射辐射亮度记录项。

本发明方法的第一部分计算3D虚拟场景中的所有漫反射和模糊反射表面的BRDF的低阶HSHBF展开系数，具体步骤如下：

步骤S101：设置最低HSHBF阶数n_l、最高HSHBF阶数n_h和HSHBF展开的误差阈值E_t；

步骤S102：将漫反射看作是模糊反射的特例，对于3D虚拟场景中的各漫反射和模糊反射表面的BRDF(A101)，标记BRDF(A101)能表示成低阶HSHBF的累加和的形式，令n_v＝n_l，做如下计算：

①按余弦分布对正半球空间进行空间角度采样，在每个角度采样

处计算(其中θ_o为极角分量，

为方位角分量)：

l＝0，1，...，n_v，m＝-l，...，l，其中

为BRDF(A101)在出射角为

入射角为

时的值，

为一个半球谐和基函数，为BRDF(A101)在角度采样上与

相关联的HSHBF展开系数；

②计算如果且n_v＜n_h，则n_v＝n_v+1，转①；如果

且n_v＝n_h，则标记BRDF(A101)不能表示成低阶HSHBF的累加和的形式；如果

则在共享内存中保存所有HSHBF展开系数并将其与BRDF(A101)及角度采样

相关联。

本发明方法的第二部分在带有双核CPU和共享内存的计算机***上并行地绘制左右眼虚拟相机的图像画面，通过对已经计算出的入射辐射亮度值进行缓存并利用基于梯度的插值方法，实现入射辐射亮度计算结果的复用，以提高左右眼虚拟相机的图像画面的绘制速度，具体步骤如下：

步骤S201：清空共享内存中的入射辐射亮度缓存；

步骤S202：在第一个CPU计算内核上运行左眼虚拟相机的图像画面的绘制程序，同时在第二个CPU计算内核上运行右眼虚拟相机的图像画面的绘制程序；左眼虚拟相机的图像画面的绘制程序执行步骤S203，右眼虚拟相机的图像画面的绘制程序执行步骤S204；

步骤S203：令I＝1，J＝1，执行如下计算：

①从左眼虚拟相机光心位置向左眼虚拟相机的像平面上的第I行、J列像素发射主光线(B01)，测试主光线(B01)与3D虚拟场景的表面是否相交；如果不相交，则将主光线(B01)对应的像素的颜色设置为背景光照值，转⑦，否则计算最近的交点位置p_i，利用蒙特卡洛方法对光源进行采样，计算交点p_i处的直接光照值；

②如果交点p_i所在表面的反射类型为镜面反射，则递归跟踪镜面反射二次光线，计算镜面反射间接光照值；如果交点p_i所在表面的反射类型为漫反射，则转③；如果交点p_i所在表面的反射类型为模糊反射，则在共享内存中查询交点p_i所在表面的BRDF是否能表示成低阶HSHBF的累加和的形式，若能，则转③，否则使用蒙特卡洛方法对交点p_i的局部坐标系(U_i，V_i，N_i)的正N_i半球空间进行重要性采样(N_i与交点p_i的法线方向同向)，直接计算模糊反射间接光照值，转⑥；

③对入射辐射亮度缓存进行访存加锁操作；从入射辐射亮度缓存中查找满足如下条件的入射辐射亮度记录项(B02)：

条件A：其中p为入射辐射亮度记录项(B02)中的位置分量成员变量，n为入射辐射亮度记录项(B02)的位置分量成员变量所表示的表面点处的法线向量，n_i为交点p_i的单位法线向量，R_i为交点p_i到所有可见对象的谐和函数平均距离，a为预先给定的精度阈值；对入射辐射亮度缓存进行访存解锁操作；

④将查找到的所有满足条件A的入射辐射亮度记录项(B02)放入一个集合S中，如果S非空，则转⑤，否则利用蒙特卡洛方法对交点p_i的局部坐标系(U_i，V_i，N_i)的正N_i半球空间进行空间角度采样，计算各角度采样方向上的入射辐射亮度值

k表示角度采样编号，根据它们计算入射辐射亮度的HSHBF展开系数

及其平移梯度

其中l＝0，1，...，n_v，m＝-l，...，l，N为角度采样数；创建一个新的入射辐射亮度记录项，其位置分量成员变量赋值为p_i，局部坐标系分量成员变量赋值为(U_i，V_i，N_i)，HSHBF展开系数分量成员变量赋值为

HSHBF展开系数平移梯度分量成员变量赋值为

对入射辐射亮度缓存进行访存加锁操作；将该入射辐射亮度记录项添加到入射辐射亮度缓存中的与位置p_i相对应的网格单元中，同时将其加入集合S中；对入射辐射亮度缓存进行访存解锁操作；

⑤根据集合S中的各元素的局部坐标系分量成员变量，旋转各元素的局部坐标系以便与交点p_i的局部坐标系对齐，利用基于梯度的插值方法，根据集合S中的入射辐射亮度记录项插值计算交点p_i的入射辐射亮度的HSHBF展开系数计算主光线(B01)在交点p_i的局部坐标系(U_i，V_i，N_i)的正N_i半球空间中所对应的角度

将

的值作为交点p_i的模糊反射间接光照值，其中

为交点p_i所在表面的BRDF的HSHBF展开系数；

⑥将交点p_i处的直接光照值和反射间接光照值相加作为交点p_i的最终光照值计算结果；

⑦J＝J+1，如果J＞N_pix，则J＝1，I＝I+1，其中N_pix表示左眼虚拟相机的图像画面的像素列数；

⑧如果I＜＝M_pix，其中M_pix表示左眼虚拟相机的图像画面的像素行数，则转①，否则转步骤S205；

步骤S204：令I′＝M_pix，J′＝N_pix，其中M_pix表示右眼虚拟相机的图像画面的像素行数，N_pix表示右眼虚拟相机的图像画面的像素列数，执行如下计算：

①从右眼虚拟相机光心位置向右眼虚拟相机的像平面上的第I′行、J′列像素发射主光线(B03)，测试主光线(B03)与3D虚拟场景的表面是否相交；如果不相交，则将主光线(B03)对应的像素的颜色设置为背景光照值，转⑦，否则计算最近的交点位置p′_i，利用蒙特卡洛方法对光源进行采样，计算交点p′_i处的直接光照值；

②如果交点p′_i所在表面的反射类型为镜面反射，则递归跟踪镜面反射二次光线，计算镜面反射间接光照值；如果交点p′_i所在表面的反射类型为漫反射，则转③；如果交点p′_i所在表面的反射类型为模糊反射，则在共享内存中查询交点p′_i所在表面的BRDF是否能表示成低阶HSHBF的累加和的形式，若能，则转③，否则使用蒙特卡洛方法对交点p′_i的局部坐标系(U′_i，V′_i，N′_i)的正N′_i半球空间进行重要性采样(N′_i与交点p′_i的法线方向同向)，直接计算模糊反射间接光照值，转⑥；

③对入射辐射亮度缓存进行访存加锁操作；从入射辐射亮度缓存中查找满足如下条件的入射辐射亮度记录项(B04)：

条件A′：

p′为入射辐射亮度记录项(B04)中的位置分量成员变量，n′为入射辐射亮度记录项(B04)的位置分量成员变量所表示的表面点处的法线向量，n′_i为交点p′_i的单位法线向量，R_i为交点p′_i到所有可见对象的谐和函数平均距离，a为预先给定的精度阈值；对入射辐射亮度缓存进行访存解锁操作；

④将查找到的所有满足条件A′的入射辐射亮度记录项(B04)放入一个集合S′中，如果S′非空，则转⑤，否则利用蒙特卡洛方法对交点p′_i的局部坐标系(U′_i，V′_i，N′_i)的正N′_i半球空间进行空间角度采样，计算各角度采样方向上的入射辐射亮度值

k表示角度采样编号，根据它们计算入射辐射亮度的HSHBF展开系数

及其平移梯度

其中l＝0，1，...，n_v，m＝-l，...，l，N为角度采样数；创建一个新的入射辐射亮度记录项，其位置分量成员变量赋值为p′_i，局部坐标系分量成员变量赋值为(U′_i，V′_i，N′_i)，HSHBF展开系数分量成员变量赋值为

HSHBF展开系数平移梯度分量成员变量赋值为

对入射辐射亮度缓存进行访存加锁操作；将该入射辐射亮度记录项添加到入射辐射亮度缓存中的与位置p′_i相对应的网格单元中，同时将其加入集合，中；对入射辐射亮度缓存进行访存解锁操作；

⑤根据集合S′中的各元素的局部坐标系分量成员变量，旋转各元素的局部坐标系以便与交点p′_i的局部坐标系对齐，利用基于梯度的插值方法，根据集合S′中的入射辐射亮度记录项插值计算交点p′_i的入射辐射亮度的HSHBF展开系数

计算主光线(B03)在交点p′_i的局部坐标系(U′_i，V′_i，N′_i)的正N′_i半球空间中所对应的角度

将的值作为交点p′_i的模糊反射间接光照值，其中

为交点p′_i所在表面的BRDF的HSHBF展开系数；

⑥将交点p′_i处的直接光照值和反射间接光照值相加作为交点p′_i的最终光照值计算结果；

⑦J′＝J′-1，如果J′＜1，则J′＝N_pix，I′＝I′-1；

⑧如果I′≥1，则转①，否则转步骤S205；

步骤S205：左右眼虚拟相机的图像画面绘制完毕。

有益效果

本发明提供了一种包含模糊反射效果的3D虚拟场景的左右眼立体画面绘制方法。由于使用双核CPU同时绘制左右眼虚拟相机的图像画面，并且采用左眼虚拟相机的图像画面从第一行像素开始绘制，右眼虚拟相机的图像画面从最后一行像素开始绘制的绘制顺序，进一步保证了左右眼虚拟相机的图像画面绘制过程中的入射辐射亮度计算结果的可重用性；通过基于梯度的插值方法，在绘制过程中重用了入射辐射亮度计算结果，因此减少了绘制3D虚拟场景左右眼立体画面的总时间。利用本发明方法绘制3D虚拟场景立体画面比传统方法快20％～35％。

附图说明

图1是3D虚拟场景中的物体表面上的点在左右眼虚拟相机的像平面上的成像示意图。

具体实施方式

为了使本发明的特征和优点更加清楚明白，下面结合具体实施例对本发明作进一步的描述。

在本实施例中，选择Intel(R)Xeon^TM双核CPU，通过加锁的方法实现两个CPU计算内核对共享内存的互斥访问。入射辐射亮度缓存设置在共享内存中，左右眼虚拟相机的图像画面绘制程序都将计算出的入射辐射亮度结果存放到入射辐射亮度缓存中，并从入射辐射亮度缓存查找临近点的入射辐射亮度进行插值计算。因此左眼虚拟相机的图像画面绘制程序的入射辐射亮度计算结果，可以直接被右眼虚拟相机的图像画面绘制程序使用；同理，右眼虚拟相机的图像画面绘制程序的入射辐射亮度计算结果，也可以直接被左眼虚拟相机的图像画面绘制程序使用。

本发明的目的在于提供一种包含模糊反射效果的3D虚拟场景的左右眼立体画面绘制方法。本发明方法需要的设备包括：带有双核CPU和共享内存的计算机***。如图1所示，当左眼虚拟相机(101)和右眼虚拟相机(102)能同时拍摄到3D虚拟场景中的物体(105)的表面上的点p₁、点p₂和点p₃时，点p₁、点p₂和点p₃的光照值将分别保存在左眼虚拟相机(101)的像平面(103)和右眼虚拟相机(102)的像平面(104)上的像素中。首先将3D虚拟场景中的漫反射和模糊反射表面的双向反射分布函数(Bi-directional ReflectanceDistribution Function，缩写为BRDF)展开成低阶半球谐和基函数(Hemispherical HarmonicsBasis Function，缩写为HSHBF)的累加和的形式，把各阶HSHBF的系数保存在共享内存中并与该表面的BRDF相关联。然后在双核CPU上同时绘制左眼虚拟相机(101)的图像画面和右眼虚拟相机(102)的图像画面，将漫反射和模糊反射表面点的入射辐射亮度展开成低阶HSHBF的累加和的形式，并保存在入射辐射亮度缓存中，通过基于梯度的入射辐射亮度插值计算，实现入射辐射亮度计算结果的重用，从而提高3D虚拟场景的左右眼立体画面的绘制速度。

本发明方法提供一种入射辐射亮度记录项数据结构，其包括位置分量、局部坐标系分量、HSHBF展开系数分量、HSHBF展开系数平移梯度分量等成员变量。

本发明方法将入射辐射亮度记录项保存在位于共享内存中的入射辐射亮度缓存中；为了便于根据空间位置检索入射辐射亮度缓存中的入射辐射亮度记录项和向入射辐射亮度缓存中添加新的入射辐射亮度记录项，入射辐射亮度缓存使用三维空间均匀网格数据结构；根据给定的三维空间位置，可以直接检索出对应网格中的所有入射辐射亮度记录项，以及向网格中添加新的入射辐射亮度记录项。

本发明方法的第一部分计算3D虚拟场景中的所有漫反射和模糊反射表面的BRDF的低阶HSHBF展开系数，具体步骤如下：

步骤S101：设置最低HSHBF阶数n_l、最高HSHBF阶数n_h和HSHBF展开的误差阈值E_t；

步骤S102：将漫反射看作是模糊反射的特例，对于3D虚拟场景中的各漫反射和模糊反射表面的BRDF(A101)，标记BRDF(A101)能表示成低阶HSHBF的累加和的形式，令n_v＝n_l，做如下计算：

①按余弦分布对正半球空间进行空间角度采样，在每个角度采样

处计算(其中θ_o为极角分量，

为方位角分量)：

l＝0，1，...，n_v，m＝-l，...，l，其中

为BRDF(A101)在出射角为

入射角为

时的值，

为一个半球谐和基函数，

为BRDF(A101)在角度采样上与

相关联的HSHBF展开系数；

②计算

如果

且n_v＜n_h，则n_v＝n_v+1，转①；如果

且n_v＝n_h，则标记BRDF(A101)不能表示成低阶HSHBF的累加和的形式；如果

则在共享内存中保存所有HSHBF展开系数并将其与BRDF(A101)及角度采样

相关联。

本发明方法的第二部分在带有双核CPU和共享内存的计算机***上并行地绘制左右眼虚拟相机的图像画面，通过对已经计算出的入射辐射亮度值进行缓存并利用基于梯度的插值方法，实现入射辐射亮度计算结果的复用，以提高左右眼虚拟相机的图像画面的绘制速度，具体步骤如下：

步骤S201：清空共享内存中的入射辐射亮度缓存；

步骤S202：在第一个CPU计算内核上运行左眼虚拟相机的图像画面的绘制程序，同时在第二个CPU计算内核上运行右眼虚拟相机的图像画面的绘制程序；左眼虚拟相机的图像画面的绘制程序执行步骤S203，右眼虚拟相机的图像画面的绘制程序执行步骤S204；

步骤S203：令I＝1，J＝1，执行如下计算：

①从左眼虚拟相机光心位置向左眼虚拟相机的像平面上的第I行、J列像素发射主光线(B01)，测试主光线(B01)与3D虚拟场景的表面是否相交；如果不相交，则将主光线(B01)对应的像素的颜色设置为背景光照值，转⑦，否则计算最近的交点位置p_i，利用蒙特卡洛方法对光源进行采样，计算交点p_i处的直接光照值；

②如果交点p_i所在表面的反射类型为镜面反射，则递归跟踪镜面反射二次光线，计算镜面反射间接光照值；如果交点p_i所在表面的反射类型为漫反射，则转③；如果交点p_i所在表面的反射类型为模糊反射，则在共享内存中查询交点p_i所在表面的BRDF是否能表示成低阶HSHBF的累加和的形式，若能，则转③，否则使用蒙特卡洛方法对交点p_i的局部坐标系(U_i，V_i，N_i)的正N_i半球空间进行重要性采样(N_i与交点p_i的法线方向同向)，直接计算模糊反射间接光照值，转⑥；

③对入射辐射亮度缓存进行访存加锁操作；从入射辐射亮度缓存中查找满足如下条件的入射辐射亮度记录项(B02)：

条件A：

其中p为入射辐射亮度记录项(B02)中的位置分量成员变量，n为入射辐射亮度记录项(B02)的位置分量成员变量所表示的表面点处的法线向量，n_i为交点p_i的单位法线向量，R_i为交点p_i到所有可见对象的谐和函数平均距离，a为预先给定的精度阈值；对入射辐射亮度缓存进行访存解锁操作；

④将查找到的所有满足条件A的入射辐射亮度记录项(B02)放入一个集合S中，如果S非空，则转⑤，否则利用蒙特卡洛方法对交点p_i的局部坐标系(U_i，V_i，N_i)的正N_i半球空间进行空间角度采样，计算各角度采样方向上的入射辐射亮度值

k表示角度采样编号，根据它们计算入射辐射亮度的HSHBF展开系数

及其平移梯度

其中l＝0，1，...，n_v，m＝-l，...，l，N为角度采样数；创建一个新的入射辐射亮度记录项，其位置分量成员变量赋值为p_i，局部坐标系分量成员变量赋值为(U_i，V_i，N_i)，HSHBF展开系数分量成员变量赋值为

HSHBF展开系数平移梯度分量成员变量赋值为

对入射辐射亮度缓存进行访存加锁操作；将该入射辐射亮度记录项添加到入射辐射亮度缓存中的与位置p_i相对应的网格单元中，同时将其加入集合S中；对入射辐射亮度缓存进行访存解锁操作；

⑤根据集合S中的各元素的局部坐标系分量成员变量，旋转各元素的局部坐标系以便与交点p_i的局部坐标系对齐，利用基于梯度的插值方法，根据集合S中的入射辐射亮度记录项插值计算交点p_i的入射辐射亮度的HSHBF展开系数计算主光线(B01)在交点p_i的局部坐标系(U_i，V_i，N_i)的正N_i半球空间中所对应的角度

将

的值作为交点p_i的模糊反射间接光照值，其中

为交点p_i所在表面的BRDF的HSHBF展开系数；

⑥将交点p_i处的直接光照值和反射间接光照值相加作为交点p_i的最终光照值计算结果；

⑦J＝J+1，如果J＞N_pix，则J＝1，I＝I+1，其中N_pix表示左眼虚拟相机的图像画面的像素列数；

⑧如果I＜＝M_pix，其中M_pix表示左眼虚拟相机的图像画面的像素行数，则转①，否则转步骤S205；

步骤S204：令I′＝M_pix，J′＝N_pix，其中M_pix表示右眼虚拟相机的图像画面的像素行数，N_pix表示右眼虚拟相机的图像画面的像素列数，执行如下计算：

①从右眼虚拟相机光心位置向右眼虚拟相机的像平面上的第I′行、J′列像素发射主光线(B03)，测试主光线(B03)与3D虚拟场景的表面是否相交；如果不相交，则将主光线(B03)对应的像素的颜色设置为背景光照值，转⑦，否则计算最近的交点位置p′_i，利用蒙特卡洛方法对光源进行采样，计算交点p′_i处的直接光照值；

②如果交点p′_i所在表面的反射类型为镜面反射，则递归跟踪镜面反射二次光线，计算镜面反射间接光照值；如果交点p′_i所在表面的反射类型为漫反射，则转③；如果交点p′_i所在表面的反射类型为模糊反射，则在共享内存中查询交点p′_i所在表面的BRDF是否能表示成低阶HSHBF的累加和的形式，若能，则转③，否则使用蒙特卡洛方法对交点p′_i的局部坐标系(U′_i，V′_i，N′_i)的正N′_i半球空间进行重要性采样(N′_i与交点p′_i的法线方向同向)，直接计算模糊反射间接光照值，转⑥；

③对入射辐射亮度缓存进行访存加锁操作；从入射辐射亮度缓存中查找满足如下条件的入射辐射亮度记录项(B04)：

条件A′：

p′为入射辐射亮度记录项(B04)中的位置分量成员变量，n′为入射辐射亮度记录项(B04)的位置分量成员变量所表示的表面点处的法线向量，n′_i为交点p′_i的单位法线向量，R_i为交点p′_i到所有可见对象的谐和函数平均距离，a为预先给定的精度阈值；对入射辐射亮度缓存进行访存解锁操作；

④将查找到的所有满足条件A′的入射辐射亮度记录项(B04)放入一个集合S′中，如果S′非空，则转⑤，否则利用蒙特卡洛方法对交点p′_i的局部坐标系(U′_i，V′_i，N′_i)的正N′_i半球空间进行空间角度采样，计算各角度采样方向上的入射辐射亮度值

k表示角度采样编号，根据它们计算入射辐射亮度的HSHBF展开系数

及其平移梯度

其中l＝0，1，...，n_v，m＝-l，...，l，N为角度采样数；创建一个新的入射辐射亮度记录项，其位置分量成员变量赋值为p′_i，局部坐标系分量成员变量赋值为(U ′_i，V′_i，N′_i)，HSHBF展开系数分量成员变量赋值为

HSHBF展开系数平移梯度分量成员变量赋值为

对入射辐射亮度缓存进行访存加锁操作；将该入射辐射亮度记录项添加到入射辐射亮度缓存中的与位置p′_i相对应的网格单元中，同时将其加入集合S′中；对入射辐射亮度缓存进行访存解锁操作；

⑤根据集合S′中的各元素的局部坐标系分量成员变量，旋转各元素的局部坐标系以便与交点p′_i的局部坐标系对齐，利用基于梯度的插值方法，根据集合S′中的入射辐射亮度记录项插值计算交点p′_i的入射辐射亮度的HSHBF展开系数

计算主光线(B03)在交点p′_i的局部坐标系(U′_i，V′_i，N′_i)的正N′_i半球空间中所对应的角度

将

的值作为交点p′_i的模糊反射间接光照值，其中

为交点p′_i所在表面的BRDF的HSHBF展开系数；

⑥将交点p′_i处的直接光照值和反射间接光照值相加作为交点p′_i的最终光照值计算结果；

⑦J′＝J′-1，如果J′＜1，则J′＝N_pix，I′＝I′-1；

⑧如果I′≥1，则转①，否则转步骤S205；

步骤S205：左右眼虚拟相机的图像画面绘制完毕。

Claims (1)

1.一种包含模糊反射效果的3D虚拟场景的左右眼立体画面绘制方法，其特征在于，所需的设备、数据结构以及实现步骤如下：

本方法的目的在于提供一种包含模糊反射效果的3D虚拟场景的左右眼立体画面绘制方法；本方法需要的设备包括：带有双核CPU和共享内存的计算机***；当左眼虚拟相机(101)和右眼虚拟相机(102)能同时拍摄到3D虚拟场景中的物体(105)的表面上的点p₁、点p₂和点p₃时，点p₁、点p₂和点p₃的光照值将分别保存在左眼虚拟相机(101)的像平面(103)和右眼虚拟相机(102)的像平面(104)上的像素中；首先将3D虚拟场景中的漫反射和模糊反射表面的双向反射分布函数展开成低阶半球谐和基函数的累加和的形式，双向反射分布函数展的英文为Bi-directional Reflectance Distribution Function，缩写为BRDF，半球谐和基函数的英文为Hemispherical Harmonics Basis Function，缩写为HSHBF，把各阶HSHBF的系数保存在共享内存中并与该表面的BRDF相关联；然后在双核CPU上同时绘制左眼虚拟相机(101)的图像画面和右眼虚拟相机(102)的图像画面，将漫反射和模糊反射表面点的入射辐射亮度展开成低阶HSHBF的累加和的形式，并保存在入射辐射亮度缓存中，通过基于梯度的入射辐射亮度插值计算，实现入射辐射亮度计算结果的重用，从而提高3D虚拟场景的左右眼立体画面的绘制速度；

本方法提供一种入射辐射亮度记录项数据结构，其包括位置分量、局部坐标系分量、HSHBF展开系数分量、HSHBF展开系数平移梯度分量；

本方法将入射辐射亮度记录项保存在位于共享内存中的入射辐射亮度缓存中；为了便于根据空间位置检索入射辐射亮度缓存中的入射辐射亮度记录项和向入射辐射亮度缓存中添加新的入射辐射亮度记录项，入射辐射亮度缓存使用三维空间均匀网格数据结构；根据给定的三维空间位置，可以直接检索出对应网格中的所有入射辐射亮度记录项，以及向网格中添加新的入射辐射亮度记录项；

本方法的第一部分计算3D虚拟场景中的所有漫反射和模糊反射表面的BRDF的低阶HSHBF展开系数，具体步骤如下：

步骤S101：设置最低HSHBF阶数n_l、最高HSHBF阶数n_h和HSHBF展开的误差阈值E_t；

步骤S102：将漫反射看作是模糊反射的特例，对于3D虚拟场景中的各漫反射和模糊反射表面的BRDF(A101)，标记BRDF(A101)能表示成低阶HSHBF的累加和的形式，令n_v＝n_l，做如下计算：

①按余弦分布对正半球空间进行空间角度采样，在每个角度采样

处计算，其中θ_o为极角分量，

为方位角分量：

l＝0，1，...，n_v，m＝-l，...，l，其中

为BRDF(A101)在出射角为入射角为

时的值，

为一个半球谐和基函数，

为BRDF(A101)在角度采样

上与

相关联的HSHBF展开系数；

②计算

如果

且n_v＜n_h，则n_v＝n_v+1，转①；如果

且n_v＝n_h，则标记BRDF(A101)不能表示成低阶HSHBF的累加和的形式；如果

则在共享内存中保存所有HSHBF展开系数

并将其与BRDF(A101)及角度采样

相关联；

本方法的第二部分在带有双核CPU和共享内存的计算机***上并行地绘制左右眼虚拟相机的图像画面，通过对已经计算出的入射辐射亮度值进行缓存并利用基于梯度的插值方法，实现入射辐射亮度计算结果的复用，以提高左右眼虚拟相机的图像画面的绘制速度，具体步骤如下：

步骤S201：清空共享内存中的入射辐射亮度缓存；

步骤S202：在第一个CPU计算内核上进行左眼虚拟相机的图像画面的绘制，同时在第二个CPU计算内核上进行右眼虚拟相机的图像画面的绘制；步骤S203为绘制左眼虚拟相机的图像画面，步骤S204为绘制右眼虚拟相机的图像画面；

步骤S203：令I＝1，J＝1，执行如下计算：

①从左眼虚拟相机光心位置向左眼虚拟相机的像平面上的第I行、J列像素发射主光线(B01)，测试主光线(B01)与3D虚拟场景的表面是否相交；如果不相交，则将主光线(B01)对应的像素的颜色设置为背景光照值，转⑦，否则计算最近的交点位置p_i，利用蒙特卡洛方法对光源进行采样，计算交点p_i处的直接光照值；

②如果交点p_i所在表面的反射类型为镜面反射，则递归跟踪镜面反射二次光线，计算镜面反射间接光照值；如果交点p_i所在表面的反射类型为漫反射，则转③；如果交点p_i所在表面的反射类型为模糊反射，则在共享内存中查询交点p_i所在表面的BRDF是否能表示成低阶HSHBF的累加和的形式，若能，则转③，否则使用蒙特卡洛方法对交点p_i的局部坐标系(U_i，V_i，N_i)的正N_i半球空间进行重要性采样，N_i与交点p_i的法线方向同向，直接计算模糊反射间接光照值，转⑥；

③对入射辐射亮度缓存进行访存加锁操作；从入射辐射亮度缓存中查找满足如下条件的入射辐射亮度记录项(B02)：

条件A：

其中p为入射辐射亮度记录项(B02)中的位置分量，n为入射辐射亮度记录项(B02)的位置分量所表示的表面点处的法线向量，n_i为交点p_i的单位法线向量，R_i为交点p_i到所有可见对象的谐和函数平均距离，a为预先给定的精度阈值；对入射辐射亮度缓存进行访存解锁操作；

④将查找到的所有满足条件A的入射辐射亮度记录项(B02)放入一个集合S中，如果S非空，则转⑤，否则利用蒙特卡洛方法对交点p_i的局部坐标系(U_i，V_i，N_i)的正N_i半球空间进行空间角度采样，计算各角度采样方向上的入射辐射亮度值

k表示角度采样编号，根据它们计算入射辐射亮度的HSHBF展开系数

及其平移梯度

其中l＝0，1，...，n_v，m＝-l，...，l，N为角度采样数；创建一个新的入射辐射亮度记录项，其位置分量赋值为p_i，局部坐标系分量赋值为(U_i，V_i，N_i)，HSHBF展开系数分量赋值为

HSHBF展开系数平移梯度分量赋值为对入射辐射亮度缓存进行访存加锁操作；将该入射辐射亮度记录项添加到入射辐射亮度缓存中的与位置p_i相对应的网格单元中，同时将其加入集合S中；对入射辐射亮度缓存进行访存解锁操作；

⑤根据集合S中的各元素的局部坐标系分量，旋转各元素的局部坐标系以便与交点p_i的局部坐标系对齐，利用基于梯度的插值方法，根据集合S中的入射辐射亮度记录项插值计算交点p_i的入射辐射亮度的HSHBF展开系数

计算主光线(B01)在交点p_i的局部坐标系(U_i，V_i，N_i)的正N_i半球空间中所对应的角度

将的值作为交点p_i的模糊反射间接光照值，其中

为交点p_i所在表面的BRDF的HSHBF展开系数；

⑥将交点p_i处的直接光照值和反射间接光照值相加作为交点p_i的最终光照值计算结果；

⑦J＝J+1，如果J＞N_pix，则J＝1，I＝I+1，其中N_pix表示左眼虚拟相机的图像画面的像素列数；

⑧如果I＜＝M_pix，其中M_pix表示左眼虚拟相机的图像画面的像素行数，则转①，否则转步骤S205；

步骤S204：令I′＝M_pix，J′＝N_pix，其中M_pix表示右眼虚拟相机的图像画面的像素行数，N_pix表示右眼虚拟相机的图像画面的像素列数，执行如下计算：

①从右眼虚拟相机光心位置向右眼虚拟相机的像平面上的第I′行、J′列像素发射主光线(B03)，测试主光线(B03)与3D虚拟场景的表面是否相交；如果不相交，则将主光线(B03)对应的像素的颜色设置为背景光照值，转⑦，否则计算最近的交点位置p′_i，利用蒙特卡洛方法对光源进行采样，计算交点p′_i处的直接光照值；

②如果交点p′_i所在表面的反射类型为镜面反射，则递归跟踪镜面反射二次光线，计算镜面反射间接光照值；如果交点p′_i所在表面的反射类型为漫反射，则转③；如果交点p′_i所在表面的反射类型为模糊反射，则在共享内存中查询交点p′_i所在表面的BRDF是否能表示成低阶HSHBF的累加和的形式，若能，则转③，否则使用蒙特卡洛方法对交点p′_i的局部坐标系(U′_i，V′_i，N′_i)的正N′_i半球空间进行重要性采样，N′_i与交点p′_i的法线方向同向，直接计算模糊反射间接光照值，转⑥；

③对入射辐射亮度缓存进行访存加锁操作；从入射辐射亮度缓存中查找满足如下条件的入射辐射亮度记录项(B04)：

条件A′：p′为入射辐射亮度记录项(B04)中的位置分量，n′为入射辐射亮度记录项(B04)的位置分量所表示的表面点处的法线向量，n′_i为交点p′_i的单位法线向量，R_i为交点p′_i到所有可见对象的谐和函数平均距离，a为预先给定的精度阈值；对入射辐射亮度缓存进行访存解锁操作；

④将查找到的所有满足条件A′的入射辐射亮度记录项(B04)放入一个集合S′中，如果S′非空，则转⑤，否则利用蒙特卡洛方法对交点p′_i的局部坐标系(U′_i，V′_i，N′_i)的正N′_i半球空间进行空间角度采样，计算各角度采样方向上的入射辐射亮度值k表示角度采样编号，根据它们计算入射辐射亮度的HSHBF展开系数及其平移梯度

其中l＝0，1，...，n_v，m＝-l，...，l，N为角度采样数；创建一个新的入射辐射亮度记录项，其位置分量赋值为p′_i，局部坐标系分量赋值为(U′_i，V′_i，N′_i)，HSHBF展开系数分量赋值为

HSHBF展开系数平移梯度分量赋值为对入射辐射亮度缓存进行访存加锁操作；将该入射辐射亮度记录项添加到入射辐射亮度缓存中的与位置p′_i相对应的网格单元中，同时将其加入集合S′中；对入射辐射亮度缓存进行访存解锁操作；

⑤根据集合S′中的各元素的局部坐标系分量，旋转各元素的局部坐标系以便与交点p′_i的局部坐标系对齐，利用基于梯度的插值方法，根据集合S′中的入射辐射亮度记录项插值计算交点p′_i的入射辐射亮度的HSHBF展开系数

计算主光线(B03)在交点p′_i的局部坐标系(U′_i，V′_i，N′_i)的正N′_i半球空间中所对应的角度

将的值作为交点p′_i的模糊反射间接光照值，其中为交点p′_i所在表面的BRDF的HSHBF展开系数；

⑥将交点p′_i处的直接光照值和反射间接光照值相加作为交点p′_i的最终光照值计算结果；

⑦J′＝J′-1，如果J′＜1，则J′＝N_pix，I′＝I′-1；

⑧如果I′≥1，则转①，否则转步骤S205；

步骤S205：左右眼虚拟相机的图像画面绘制完毕。

Images