天气雷达三维回波实时渲染算法研究.docxVIP

? ? 天气雷达三维回波实时渲染算法研究 ? ? 王 兴， 钱代丽， 朱 彬，， 安俊琳， 吕晶晶 （南京信息工程大学 a.大气科学与环境气象国家级实验教学示范中心，b.大气科学与气象信息国家级虚拟仿真实验教学中心，南京 210044） 0 引 言 曲面重建是对离散的空间数据进行三维渲染的主要技术手段，其基本思想是按照一组给定的阈值从三维点云数据中提取相应的等值面，利用曲面渲染算法生成一组由不同颜色相互叠加的等值面［5-6］。曲面重建的算法很多，主要可分为两类：一是基于断层轮廓线的表面重建，即将三维数据分割成若干个二维切面，分别计算各切面的等值线，再将这些等值线连接成等值曲面；二是基于体素的等值面重建，即对三维点云构成的体数据场中的各个体素加以分析，提取等值面信息，再将这些等值面连接成等值曲面［7-8］。其中，移动立方体法（Marching Cubes，MC）［9］是沿用至今且颇具影响力的体素等值面重建算法之一。 曲面重建技术在医疗、建筑等工程领域的研究最为深入，近年来，也有一些研究将该技术应用于天气雷达、卫星和数值模式等资料，实现气象数据的三维可视化交互［10-11］。由于雷达探测数据的空间分布极不规则，距离雷达近端与远端的相邻离散点之间，间隔距离相差超过100倍，在重构均匀三维格点场雷达回波时，传统空间插值算法并不适用。雷达远端的插值点不仅难以还原真实的探测值，还极大地增加了点云的密度，造成面片的计算规模陡增，网络结构复杂化，大幅增加了三维渲染的时长，影响三维可视化的实时性和可交互性。 为解决上述技术难题，针对多普勒天气雷达探测方式及其数据特征，在移动立方体算法基础上，提出一系列改进，并对改进前、后的渲染效果做出比对分析。 1 基本原理与技术路线 1.1 移动立方体（MC）算法 MC算法最早是Lorensen等［9］提出，用于解决三维物体的面绘制问题。其本质是将三维坐标系下的物体分解为体素集合，每个体素包含8个空间上紧相邻的顶点，构成一个正立方体，将该正立方体视为基本单位。曲面重建时，先设定一组大小不同的阈值，通过计算每个体素内的梯度变化，并与各个阈值分别进行比较，找出所有被阈值穿过的正立方体，再利用插值算法计算出这些表面（即等值面），再通过绘制一系列阈值的等值面，实现三维矢量图像的绘制。 受制于20世纪80年代有限的计算机图形处理能力，MC算法巧妙地利用三角面片近似表示待求取的等值面，极大地提升了计算效率。由于等值面与体素边界面的交线是双曲线，在特定情况下对相同的等值点可以采取不同的连接方式，这种近似表示的算法无法保证三角片所构成的等值面的拓扑一致性，造成很多三角面片衔接处出现孔隙和不连续现象。为克服这一问题，很多学者对MC算法提出了改进，有的是利用插值和双曲渐近线将体素进一步拆分［12-13］，有的是将MC的正立方体进一步简化，采用四面体剖分［14］等。 1.2 MC应用于雷达回波渲染的技术缺陷 气象雷达探测得到的数据，其空间分布是一个不规则的三维点云结构，随着探测距离的增加，点云愈发稀疏，距离雷达近端与远端的相邻离散点之间，间隔距离相差可达100倍以上。在重构均匀三维各点场雷达回波时，诸如双线性插值、最近邻插值、三次内插等传统空间插值算法并不适用。 宴会中，马老当场宣布，这次活动的成功举办，在座的各位功不可没。公司将论功行赏，给予在座的各位同仁以物质奖励，对公司的高管层将另行奖赏。奖励田卓小姐、高潮先生赴巴厘岛六日游，冯可儿小姐本来也在去巴厘岛之列，但考虑到公司不可一日无主，就留下来代理田卓行使管理权，但公司会额外补偿冯可儿小姐不能成行的损失。 尽管MC利用三角面片拟合曲面的算法设计非常巧妙，在曲面渲染时也表现出不错的性能，但对于面片的生成，却需要消耗“过多”的内存和计算时长。这主要是因为MC算法需要遍历三维点云数据中的所有体素，而与等值面相交的体素数量仅占所有体素的1/10甚至更低。雷达探测的空间范围大（数百公里），数据密度高，一次计算需要处理数百万个体素，而近9成的时间被用于计算这些无效体素。雷达远端的插值点不仅难以还原真实的探测值，还极大地增加了点云的密度，使得体素产生的面片面积过小，造成面片的计算规模陡增，网络结构复杂化，大幅增加了三维渲染的时长，影响了三维回波渲染的实时性和可交互性。 1.3 雷达三维回波实时渲染算法流程 为克服上述不足，本算法以传统MC算法为基础，针对雷达基数据的空间结构特点，着重在三角面片顶点位置的计算、三角剖分构型唯一性的判定以及包含等值面体素的加速遍历，这3方面提出修改和优化。雷达三维回波实时渲染的总体算法流程如图1所示。 图1 总体算法流程图 2 雷达三维回波实时渲染算法 2.1 数据准备 国内常用的天气雷达分为S波段、C波段和X波段3类，每个波段各有多种型号，每个型号的雷