[发明专利]一种基于GPGPU的天然气泄漏扩散的动态可视化方法

权利要求书

1.一种基于GPGPU的天然气泄漏扩散的动态可视化方法，其特征在于：包括以下过程，

首先，初始化绘制环境，使用OpenGL创建两个PBO和一个FBO，创建三维纹理对象和二维纹理对象用于存储GPGPU的计算结果，定义扩散坐标系和体素坐标系，初始化视点、光源在WGS84坐标系中的位置；

设当前的扩散时刻为t，先基于CPU进行以下泄露扩散范围的顶点阵列生成处理，

步骤a1，由CPU在扩散坐标系中，计算扩散时刻t时天然气分别在x轴、y轴、z轴的扩散距离Rx、Ry、Rz和扩散范围在x轴的最小值Minx；

步骤a2，由CPU在扩散坐标系中，将扩散范围内的采样距离与三维场景的LOD级别相关联，在当前三维场景的LOD级别下以deltaX、deltaY、deltaZ的距离分别在x轴、y轴、z轴方向采样，生成t时刻泄露扩散范围的顶点阵列，并重新调整Rx、Ry、Rz，使得x轴、y轴、z轴方向的采样点个数均为2的整次幂；将扩散坐标系中泄漏扩散范围的体包围盒转换到WGS84坐标系，继而转换到世界坐标系中，得到世界坐标系中泄漏扩散范围的体包围盒；由两个平行平面与世界坐标系中泄漏扩散范围的体包围盒构建扩大之后的新的体包围盒，其中一个平面经过((Rx+Minx)/2，0，Rz)在世界坐标系中的对应点，另一个平面经过世界坐标系中泄漏扩散范围体包围盒底面的四个顶点；

步骤a3，由CPU将t时刻泄露扩散范围的顶点阵列载入一个空闲的PBO，令t＝t+deltaT，deltaT为时间间隔，针对下一个当前的扩散时刻返回步骤a1，执行步骤a1-a3，将新的扩散时刻t泄露扩散范围的顶点阵列载入另一个PBO；

当CPU将t时刻泄露扩散范围的顶点阵列载入PBO后，针对下一扩散时刻执行步骤a1-a3时，由GPU将顶点阵列从PBO传输到与顶点阵列大小相同的三维纹理对象，并完成以下绘制处理，

绘制两个长方体，将浓度值体数据和顶点光照透射率体数据通过纹理映射的方式在片段着色器中实现GPGPU计算，把计算结果渲染并输出到FBO，则浓度值体数据和顶点光照透射率体数据分别作为三维纹理被写入纹理缓存；绘制一个矩形，将屏幕点颜色值数据通过纹理映射的方式，结合对浓度值体数据三维纹理和顶点光照透射率体数据三维纹理的采样，在片段着色器中实现GPGPU计算，计算结果写入二维纹理并直接渲染到屏幕，实现了扩散浓度数据生成和可视化的一体化。

2.如权利要求1所述一种基于GPGPU的天然气泄漏扩散的动态可视化方法，其特征在于：

步骤1中，对于扩散坐标系和体素坐标系的定义如下，

扩散坐标系定义为以泄漏源p₀在地面的投影点为原点，以下风向方向为x轴，以地面上垂直于x轴方向为y轴，以垂直于地面竖直向上方向为z轴；在扩散坐标系中，扩散范围由长方体包围盒(Minx，-Ry，0)、(Rx，-Ry，0)、(Rx，Ry，0)、(Minx，Ry，0)、(Minx，-Ry，Rz)、(Rx，-Ry，Rz)、(Rx，Ry，Rz)、(Minx，Ry，Rz)确定，Rx、Ry、Rz分别为x轴、y轴、z轴的最大扩散距离，Minx为扩散范围在x轴的最小值；

体素坐标系定义为以扩散坐标系中(Minx，-Ry，0)为原点，x轴、y轴、z轴方向均与扩散坐标系相同；在体素坐标系中，扩散范围由长方体包围盒(0，0，0)、(M-1，0，0)、(M-1，N-1，0)、(0，N-1，0)、(0，0，L-1)、(M-1，0，L-1)、(M-1，N-1，L-1)、(0，N-1，L-1)确定，M-1、N-1、L-1分别为x轴、y轴、z轴的最大扩散距离，M、N、L分别表示x轴、y轴、z轴的采样点个数。