CN102254347A - 电磁环境体数据三维可视化方法 - Google Patents

Abstract

本发明属于电磁环境可视化技术领域，公开了一种电磁环境体数据三维可视化方法。本发明采用电波传播物理模型确定电磁环境体数据，把电磁环境体数据剖分成四面体结构，按照四面体顶点和顶点索引的方式组织成二维纹理，输入四面体顶点函数值纹理和顶点索引纹理到图形处理器，根据设定的等值面阈值判断生成等值面顶点，将同一阈值的等值面顶点连接成等值面，构造电磁环境多等值面，根据等值面顶点法向量和视点关系突出电磁环境模型边缘显示，并将多等值面半透明融合嵌入到三维环境中。采用本发明可以使三维电磁环境表现得形象直观，提供更多的电磁环境内部信息表现，提高了电磁环境表现的正确性、完备性以及边界轮廓的清晰度。

Description

电磁环境体数据三维可视化方法

技术领域

本发明涉及电磁环境可视化技术领域，尤其涉及一种电磁环境体数据三维可视化方法。

背景技术

随着电磁技术在社会各个领域中的广泛应用，电磁环境呈现出越来越复杂的态势，由于电磁环境的不可见，人们难以直观分析、判断电磁环境。目前使用的用频设备类型广泛，电磁设备数量多，占用频谱范围宽，频率需求量大，电磁环境极为复杂。但是电磁环境看不见、摸不着，只能通过各种侦测设备去分析空间中的电磁信息，而且仅能反映局部区域的电磁分布特征，难以满足用户对整个电磁态势的把握，制约了对电磁态势做出快速准确的判断。

电磁环境是指存在于给定场所的电磁现象的总和，它是由各种电台、雷达等电磁设备辐射电磁波到三维空间中形成的。电磁波在空间传播，由于受大气、地形等多方面因素的影响，发射的能量不可避免地会随着传播距离的增加而逐渐衰减，通过电波传播衰减值和发射的能量值就能得到空间某点处的电磁波能量值。电波传播模型通常用于预测发射机到接收机之间的电波传播情况，估计各种环境因素影响下无线通信性能。HuffordG.A给出了Longley-Rice电波传播模型(参见Hufford G.A.The ITS Irregular TerrainModel，version 1.2.2，The Algorithm[R].http://flattop.its.bldrdoc.gov/itm/itm_alg.pdf，National Telecommunications and Information Administration Institute forTelecommunication Sciences，USA.p6-12)。电波传播模型作为设计各种无线电***的基础，可以通过各种电波传播模型来构建电磁环境。虽然电波传播模型可以估计电磁环境，但是得到的只是大量描述电磁环境的数据值即电磁环境体数据，人们难以从这些纷繁复杂的电磁环境体数据中直观察觉到电磁环境在三维空间中的分布情况，因此把电磁环境体数据以直观可视的方式反映给用户，并且能最大程度地完整地通过三维形式描述这些电磁环境体数据，进而让用户准确详尽地把握电磁环境成为目前电磁环境可视化的迫切需求。

传统对电磁环境的可视化主要是在获取电磁环境数据后，通过二维等值线的方式勾勒某高度上的电磁环境分布情况，或者通过二维数据链路的方式表现电磁环境的辐射连通情况。二维等值线通过提取电磁环境等值的二维线条并填充线条覆盖区域的方式来表现电磁环境在某高度上的分布情况，等值二维线条及其覆盖区域表示了电磁环境某个能量值的分布范围。二维数据链路的方式则是通过二维线条表示电磁设备间的连通情况，通过线条颜色、线条粗细等表示连接通信能力。通过这些二维可视化方式，仅能表示电磁环境的部分信息，不能完整描述电磁环境在三维空间中的分布情况，丢失了电磁环境三维分布信息。

等值面技术是空间数据场可视化领域中广泛使用的一种技术，通过提取体数据特定大小属性的轮廓信息，构建起三维网格模型，能清晰展示体数据特定属性值的轮廓范围，而且由于目前图形硬件支持网格模型绘制加速，绘制效率高，因此该方法是一种高效可视化体数据的技术。通过等值面提取技术，能构建起三维环境中电磁环境的三维网格模型，并实时绘制到三维环境中，对于电磁环境三维特征轮廓的直观表达有着重要的意义。从目前最新的研究报告看(参见：陈鹏.虚拟战场环境中雷达作用范围表现技术研究[D].湖南长沙：国防科学技术大学，2007.周桥，陈景伟等.电磁环境三维可视化技术[J].计算机工程，2008.)，对雷达、电台等典型装备的电磁环境覆盖范围的可视化已经采用了等值面的可视化技术，通过提取电磁环境某个特征值的等值面，能对该特征值的电磁环境在三维环境中的分布情况进行表现，但是由于它采用的传统等值面技术只通过单个轮廓信息表达电磁环境信息，这样只能将原始数据的部分属性映射为网格轮廓模型，不能反映整个数据全貌及细节，这样会造成大量空间数据场中的内部细节信息丢失，不能完整描述电磁环境信息，不利于电磁环境这样的空间数据场的信息可视化。

发明内容

本发明要解决的技术问题是：实现电磁环境的三维可视化表现，能实时与电磁环境交互，动态掌握电磁环境的变化，提高电磁环境表现的正确性、精确性以及完备性。

为解决上述技术问题，本发明提出了一种电磁环境体数据三维可视化方法，该方法包括如下步骤：

步骤一、确定电磁环境体数据

根据已知电磁设备参数、地形高程数据和气象环境数据，通过电波传播物理模型确定电磁环境体数据。

步骤二、电磁环境体数据六面体结构剖分

把电磁环境体数据六面体结构剖分成四面体结构，遍历每一个四面体并获得该四面体各个顶点的顶点位置和电磁环境体数据值，并按照四面体顶点和顶点索引的方式组织成二维纹理。

步骤三、判断生成等值面顶点

输入四面体顶点函数值纹理和顶点索引纹理到图形处理器，根据设定的等值面阈值大小判断生成等值面顶点。

步骤四、构造电磁环境多等值面

遍历所有四面体，将同一阈值的等值面顶点连接成等值面，构造电磁环境多等值面。步骤五、突出电磁环境多等值面边界轮廓

根据多等值面顶点法向量和视点的关系设定多等值面的颜色透明度，以多等值面颜色透明度来突出电磁环境多等值面边界轮廓。

步骤六、获得电磁环境多等值面半透明结果

采用快速分层伪排序融合绘制算法获得电磁环境多等值面半透明结果，并将电磁环境多等值面半透明结果融合嵌入到三维环境中。

进一步，本发明所述步骤一中所述电磁环境体数据确定的具体方法是：

(1)用户初始设定纬度θ＝θ₀、经度

高度r＝r₀，计算终止的纬度设为θ₁，经度为

高度为r₁，并设定纬度、经度和高度步长为Δθ、

Δr；

(2)输入电磁设备参数，地形高程数据和气象环境数据，采用电波传播物理模型，得到电磁环境体数据值

(3)经度递增

重复步骤(2)，直到经度

(4)纬度θ递增Δθ，重复步骤(2)，直到纬度θ＞θ₁；

(5)纬度r递增Δr，重复步骤(2)，直到高度r＞r₁。

本发明所述步骤二中所述按照四面体顶点和索引的方式组织成二维纹理的具体方法是：

(1)将电磁环境体数据剖分成四面体结构；

(2)按顺序将电磁环境体数据所有的顶点坐标存储下来，每个顶点仅存储一次，采用RGBA32F格式的二维纹理，称为顶点函数值纹理，每个顶点对应一个纹理纹素，纹素的R、G、B、A分量对应顶点的三维坐标和电磁环境体数据值，并将该纹素的顺序号作为顶点的索引；

(3)存储每个四面体顶点的顶点索引到另一个二维纹理，称为顶点索引纹理，其格式为RGBA32UI，每个四面体对应一个纹理纹素，纹理纹素有R、G、B、A分量，每个纹理纹素分量对应四面体的一个顶点索引。

本发明所述步骤三中所述根据设定的等值面阈值大小生成等值面顶点的方法是：

1、对于电磁环境体数据中其中一个四面体，判断该四面体的一条边上是否与给定阈值δ的等值面有交点，有交点即等值面的顶点。具体判断方法为：

(1)根据四面体其中一条边的端点，在顶点索引纹理中得到顶点索引；

(2)根据顶点索引在四面体顶点函数值纹理中获取该边两个端点对应的纹素，该纹素的R、G、B、A分量对应顶点的三维坐标和电磁环境体数据值，从而可以获得两个端点的顶点坐标v_a、v_b和电磁环境体数据值P_a、P_b；

(3)比较多等值面的阈值δ与端点的电磁环境体数据值，如果δ＞P_a且δ＞P_b，或者δ＜P_a且δ＜P_b则该边上不存在该阈值大小δ的等值面的顶点；否则存在等值面的顶点，插值生成等值面顶点为：

v＝(1-t)v_a+tv_b

其中， $t=\frac{\mathrm{\δ}-P_a}{P_b-P_a}$

(4)重复步骤(1)、(2)、(3)判断生成该四面体所有边上阈值δ的等值面的顶点；

2、遍历所有四面体完成步骤1，生成该阈值的等值面的所有顶点。

3、遍历所有等值面阈值，完成步骤1、2，生成所有阈值的等值面顶点。

本发明所述步骤四中所述电磁环境多等值面构造方法如下：

(1)对于每个四面体，根据其六条边的排列顺序，顺序连接与边相交的同一阈值的等值面顶点，即可得到某一阈值的电磁环境等值面的一个面片；

(2)遍历多等值面的所有阈值完成步骤(1)，即可得到该四面体的电磁环境多等值面面片；

(3)对所有四面体完成步骤(1)、(2)，即可得到由多个等值面面片构成的电磁环境多等值面。

本发明所述步骤五中所述多等值面颜色透明度设置的具体方法是：

(1)遍历所有多等值面的顶点，对每个顶点计算该顶点法矢量与视点方向的夹角余弦：

cosα＝n^Tv

其中n为等值面顶点法矢量的归一化结果，v为视点与该顶点的视线方向归一化结果；

(2)设置该顶点的透明度值为：

A＝1-k|cosα|

其中k为边界轮廓控制正数因子，取值越大透明区域越大，边界轮廓区域越突出(3)对多层等值面模型的所有顶点完成步骤(1)、(2)。

本发明所述步骤六中所述快速分层伪排序融合绘制算法如下：

(1)依据电磁环境多等值面提取阈值的大小，对每个等值面按照从内向外排序；

(2)绘制所有等值面的背面面片，实现等值面的背面面片按照从后往前的顺序融合；

(3)绘制所有等值面的正面面片，实现等值面的正面面片按照从后往前的顺序融合；

(4)融合多等值面正面面片和背面面片绘制结果，得到半透明结果；

(5)根据电磁环境多等值面的顶点法矢量，采用phong着色模型计算光照情况，为电磁环境多等值面的所有顶点着色，最终将带有光照效果的电磁环境多等值面显示在三维环境中。

综上所述，本发明有如下优点

(1)由于本方法采用步骤四方法构成了电磁环境多等值面，所以能够完整的表现电磁环境体数据的内部细节信息；

(2)采用步骤五的颜色透明度设置方法，能够突出电磁环境多等值面的边界轮廓信息；

(3)采用本发明的方法，能从多个角度查看电磁环境，能够调整电磁设备参数动态查看电磁环境，提高电磁环境表现的正确性、精确性以及完备性。

附图说明

图1为电磁环境体数据三维可视化方法总体流程图。

图2为电磁环境体数据中一个六面体剖分成六个四面体的示意图。

图3为电磁环境等值面与图2中一个四面体可能的交点个数示意图，

(a)与四面体无交点示意图；

(b)与四面体三个交点示意图；

(c)与四面体四个交点示意图。

图4为电磁环境多等值面其中一个顶点生成流程图。

图5为电磁环境等值面与一个四面体六条边的16种相交情况示意图。

图6为电磁环境多等值面边界轮廓增强示意图。

图7为电磁环境多等值面拓扑结构示意图。

图8为多等值面半透明融合绘制流程图。

图9为三维环境中展示的受环境影响后的电磁环境多等值面示意图，

(a)表示从侧面观察电磁环境多等值面，

(b)表示从顶部观察电磁环境多等值面。

具体实施方式

下面结合附图所描述的实施方式对本方法作进一步详细说明。

图1为本发明总体流程图，本发明提出的电磁环境体数据三维可视化方法，包括如下步骤：

步骤一、确定电磁环境体数据

根据电磁设备参数、地形高程数据和气象环境数据，通过电波传播物理模型确定电磁环境体数据。

电磁环境体数据是由电磁设备在三维空间中辐射电磁波形成的电磁能量强度场，它的大小主要由电磁波传播衰减决定的，电磁波传播衰减是指电磁波在传播路径上损失的能量，传播路径是电磁波从发射机到接收机之间的直线距离。

电磁环境体数据具体确定方法是：

1、由Longley-Rice电波传播模型，确定电磁设备辐射电磁波到纬度θ、经度

高度r点处的损失值为

其中，d_ls是平坦球面视距距离，d_x是散射区域计算起始距离。A_ed是一个考虑了大气、刃形衰减和地形遮挡等因素的电波传播损失值，m_d是表示在绕射区域内，衰减值与距离成线性关系。A_el是A_ed在d_ls处的值，K₁、K₂是在视距内插值控制系数。A_es是A_ed在d_x处的值，m_s是在散射区域衰减值线性插值因子。

Longley-Rice电波传播模型参见Hufford G.A.The ITS Irregular Terrain Model，version 1.2.2，The Algorithm[R].http://flattop.its.bldrdoc.gov/itm/itm_alg.pdf，NationalTelecommunications and Information Administration Institute for TelecommunicationSciences，USA.p6-12.

2、计算电磁设备辐射电磁波到纬度θ、经度

高度r点处的电磁环境体数据值为

其中P_t是电磁设备发射的功率强度值，G是天线增益，f是电磁波的频率，R是电磁波的传播路径长度，L_s是***损耗，取3-5dB。

通过其他电波传播物理模型都可求得电磁环境体数据值。

3、确定三维空间的电磁环境体数据

电磁环境体数据包含的范围是从纬度θ＝θ₀、经度高度r＝r₀，到纬度为θ₁，经度为

高度为r₁，并且纬度、经度和高度步长为Δθ、

Δr。电磁环境体数据中的每个数据点值由公式(2)确定，具体方法是：

(1)用户初始纬度θ＝θ₀、经度高度r＝r₀，计算终止的纬度为θ₁，经度为

高度为r₁，并设定纬度、经度和高度步长为Δθ、

Δr；

(2)输入电磁设备参数、地形高程数据和气象环境数据，根据公式(2)计算电磁环境体数据值

(3)经度

递增

重复步骤(2)，直到经度

(4)纬度θ递增Δθ，重复步骤(2)，直到纬度θ＞θ₁；

(5)高度r递增Δr，重复步骤(2)，直到高度r＞r₁。

步骤二、电磁环境体数据六面体结构剖分

图2(以立方体为例)为一个六面体剖分成六个四面体的示意图，电磁环境体数据中空间相邻的8个数据点v₀、v₁、v₂、v₃、v₄、v₅、v₆、v₇构成的六面体，其剖分成的6个四面体为：

T0＝(v0，v1，v4，v 7)，

T1＝(v1，v4，v5，v7)，

T2＝(v0，v1，v3，v7)，

T3＝(v1，v5，v6，v7)，

T4＝(v1，v2，v6，v7)，

T5＝(v1，v2，v3，v7)；

把电磁环境体数据六面体结构剖分成6个四面体结构，遍历每一个四面体并获得该四面体各个顶点的顶点位置和电磁环境体数据值，并按照四面体顶点和顶点索引的方式组织成二维纹理。

将电磁环境体数据中空间相邻的v₀、v₁、v₂、v₃、v₄、v₅、v₆、v₇8个数据点相连所构成的六面体剖分成6个四面体。由于相邻六面体的剖分方式完全相同，相邻六面体的共面上不会出现交叉剖分，这样能保证相邻六面体共面的剖分一致性，因此提取的等值面不会出现裂缝现象，具体方法如下：

1、将电磁环境体数据剖分成四面体结构，其具体方法为：

(1)用户初始纬度θ＝θ₀、经度

高度r＝r₀，计算终止的纬度为θ₁，经度为高度为r₁，并设定纬度、经度和高度步长Δθ、

Δr；

(2)对于电磁环境体数据中空间相邻的8个数据点为

其剖分成的6个四面体为：

T0＝(v0，v1，v4，v7)，

T1＝(v1，v4，v5，v7)，

T2＝(v0，v1，v3，v7)，

T3＝(v1，v5，v6，v7)，

T4＝(v1，v2，v6，v7)，

T5＝(v1，v2，v3，v7)；

(3)经度

递增

重复步骤(2)，直到经度

(4)纬度θ递增Δθ，重复步骤(2)，直到纬度θ＞θ₁-Δθ；

(5)纬度r递增Δr，重复步骤(2)，直到高度r＞r₁-Δr。

2、把剖分成的电磁环境体数据四面体结构组织成二维纹理，分别存储电磁环境体数据四面体的顶点和四面体连接关系两部分，四面体连接关系即四面体的四个顶点索引。具体方法为：

(1)按顺序将电磁环境体数据所有的顶点坐标存储下来，每个顶点仅存储一次，采用RGBA32F格式的二维纹理，称为顶点函数值纹理，每个顶点对应一个纹理纹素，纹素的R、G、B、A分量对应顶点的三维坐标和电磁环境体数据值，并将该纹素的顺序号作为顶点的索引；

(2)存储每个四面体顶点的顶点索引到另一个二维纹理，称为顶点索引纹理，其格式为RGBA32UI，每个四面体对应一个纹理纹素，纹理纹素有R、G、B、A分量，每个纹理纹素分量对应四面体的一个顶点索引。

步骤三、判断生成等值面顶点

输入四面体顶点函数值纹理和顶点索引纹理到图形处理器，根据设定的等值面阈值大小判断生成等值面顶点。

针对每个四面体的每条边，通过顶点索引纹理存储的顶点索引可以在四面体顶点函数值纹理中查找到该边两个端点的顶点坐标和电磁环境体数据值，然后通过等值面阈值与端点的电磁环境体数据值比较，确定该边是否存在等值面与边的交点，如果存在等值面与边的交点即等值面上的顶点。图3为等值面与四面体可能的交点个数示意图，分别为无交点(图3a)、三个交点(图3b)和四个交点(图3c)三种情况，其中“+”表示该顶点处的电磁环境体数据值大于等值面阈值，“-”号表示该顶点处的电磁环境体数据值不大于等值面阈值，斜划线表示由交点构成的等值面面片。多等值面其中一个相交点的生成流程图如图4所示。

等值面顶点生成具体方法是：

1、对于电磁环境体数据中其中一个四面体，判断该四面体的一条边上是否与给定阈值δ的等值面有交点，有交点即等值面的顶点。具体判断方法为：

(1)根据四面体其中一条边的端点，在顶点索引纹理中得到顶点索引；

(2)根据顶点索引在四面体顶点函数值纹理中获取该边两个端点对应的纹素，该纹素的R、G、B、A分量对应顶点的三维坐标和电磁环境体数据值，从而可以获得两个端点的顶点坐标v_a、v_b和电磁环境体数据值P_a、P_b；

(3)比较多等值面的阈值δ与端点的电磁环境体数据值，如果δ＞P_a且δ＞P_b，或者δ＜P_a且δ＜P_b则该边上不存在该阈值大小δ的等值面的顶点；否则存在等值面的顶点，插值生成等值面顶点为：

v＝(1-t)v_a+tv_b               (3)

其中， $t=\frac{\mathrm{\δ}-P_a}{P_b-P_a}$

(4)重复步骤(1)、(2)、(3)判断生成该四面体所有边上阈值δ的等值面的顶点；

2、遍历所有四面体完成步骤1，生成该阈值的等值面的所有顶点。

3、遍历所有等值面阈值，完成步骤1、2，生成所有阈值的等值面顶点。

步骤四、构造电磁环境多等值面

遍历所有四面体，将同一阈值的等值面顶点连接成等值面，构造电磁环境多等值面。

参见图5为等值面与一个四面体六条边的16种相交情况，其中包括两种没有相交点的情况(见图5中四面体0000和四面体1111)。根据这个相交情况示意图，可以快速确定相交点连接构成等值面面片。

构造电磁环境多等值面的具体方法如下：

1、对于每个四面体，根据其六条边的排列顺序，顺序连接与边相交的同一阈值的等值面顶点，即可得到某一阈值的电磁环境等值面的一个面片；

2、遍历多等值面的所有阈值完成步骤1，即可得到该四面体的电磁环境多等值面面片；

3、对所有四面体完成步骤1、2，即可得到由多个等值面面片构成的电磁环境多等值面。

步骤五、突出电磁环境多等值面边界轮廓

根据多等值面顶点法向量和视点的关系设定多等值面的颜色透明度，以多等值面颜色透明度来突出电磁环境多等值面边界轮廓；

为了增强多等值面边界轮廓信息，提高融合绘制结果的信息表现力，采用了逐像素着色的方法，提高电磁环境边缘轮廓显示强度，具体方法是：

1、计算等值面的所有顶点法向量：

(1)针对当前要计算法线的模型顶点，计算与该顶点邻接的多边形的法线矢量：

(a)任取等值面面片多边形不在一条直线上的三个顶点v₁、v₂和v₃，该多边形的法线矢量为这两条直线的叉乘(v₁-v₂)×(v₂-v₃)；

(b)继续对该顶点邻接的其它多边形继续步骤(a)，直到完成所有邻接多边形的法线矢量计算。

(2)把邻接多边形的法线矢量求平均，计算结果就是该顶点的法线矢量。

(3)继续对其它顶点完成步骤(1)、(2)，直到多等值面的所有顶点都完成法线矢量计算。

2、分别对多等值面的每个顶点计算其法矢量与视点方向的夹角，通过夹角大小设置等值面颜色透明度。图6为多等值面边界轮廓增强示意图，图中表示了等值面的面片片段，其中n为等值面顶点法矢量的归一化结果，v为视点与该顶点的视线方向归一化结果，通过多等值面的顶点法矢量与视点方向的夹角α来设置等值面颜色透明度。

具体方法为：

(1)遍历所有多等值面的顶点，对每个顶点计算夹角余弦为：

cosα＝n^Tv                   (4)

其中n为等值面顶点法矢量的归一化结果，v为视点与该顶点的视线方向归一化结果。

(2)设置该顶点的透明度值为：

A＝1-k|cosα|                (5)

其中k为边界轮廓控制正数因子，取值越大透明区域越大，边界轮廓区域越突出。

(3)对多等值面的所有顶点完成步骤(1)、(2)。

步骤六、获得电磁环境多等值面半透明结果

采用快速分层伪排序融合绘制算法获得电磁环境多等值面半透明结果，并将电磁环境多等值面半透明结果融合嵌入到三维环境中。

由于多等值面相互嵌套，如果采用不透明的绘制方式，则内层的等值面会被外层的等值面遮挡，只能展示最外层的电磁环境等值面信息，不能完全显示出多等值面信息，因此为了在三维环境中清晰地展示电磁环境多等值面，需要采用半透明融合绘制方法，让内层等值面透过外层等值面显现出来，提出了快速分层伪排序融合绘制算法。图7为多等值面拓扑结构，分别用不同颜色表示了四层不同阈值的等值面。示意图左部分是多等值面三维绘制效果图，右部分是对应多等值面的纵剖面二维效果图，四层等值面从外向里分别用G、Y、O、R表示，相应的等值面阈值大小是δ_G、δ_Y、δ_O、δ_R，则等值面阈值是按照δ_G、δ_Y、δ_O、δ_R有序排列。因此，根据多等值面的阈值大小可以确定多等值面的相对关系，能方便对G、Y、O、R等值面进行排序。图8为多等值面快速伪排序半透明融合绘制流程图，整个绘制流程分成预处理和绘制两个阶段。在预处理阶段，依据电磁环境多等值面的相对排列关系，对每层等值面按照其提取阈值的大小从内向外排序；在绘制阶段，采用图形处理器的模型面片裁剪功能，分别绘制电磁环境多等值面的正面和背面面片，实现电磁环境多等值面的近似由后及前的排序绘制。快速分层伪排序融合绘制算法具体为：

1、依据电磁环境多等值面提取阈值的大小，对每个等值面按照从内向外排序；

2、绘制所有等值面的背面面片，实现等值面的背面面片按照从后往前的顺序融合，具体方法为：

(1)开启图形处理器的正面裁剪功能；

(2)采用步骤五的颜色透明度设置方法，绘制最外层等值面，再绘制次外层的等值面，以此类推，直到绘制完所有等值面。

3、绘制所有等值面的正面面片，实现等值面的正面面片按照从后往前的顺序融合，具体方法为：

(1)开启图形处理器的背面裁剪功能；

(2)采用步骤五的颜色透明度设置方法，绘制内层的等值面，再绘制次内层的等值面，同样以此类推，直到绘制完所有等值面。

4、融合多等值面正面面片和背面面片绘制结果，得到最终半透明结果。

5、根据电磁环境多等值面的顶点法矢量，采用phong着色模型计算光照情况，为电磁环境多等值面的所有顶点着色，最终将带有光照效果的电磁环境多等值面显示在三维环境中。最终三维环境中展示的受环境影响后的电磁环境多等值面示意图，如图9所示，其中图9(a)表示从侧面观察电磁环境多等值面，图9(b)表示从顶部观察电磁环境多等值面。

Claims (7)

1.一种电磁环境体数据三维可视化方法，包括如下步骤：

步骤一、确定电磁环境体数据

根据已知电磁设备参数、地形高程数据和气象环境数据，通过电波传播物理模型确定电磁环境体数据；

步骤二、电磁环境体数据六面体结构剖分

把电磁环境体数据六面体结构剖分成四面体结构，遍历每一个四面体并获得该四面体各个顶点的顶点位置和电磁环境体数据值，并按照四面体顶点和顶点索引的方式组织成二维纹理；

步骤三、判断生成等值面顶点

输入四面体顶点函数值纹理和顶点索引纹理到图形处理器，根据设定的等值面阈值大小判断生成等值面顶点；

步骤四、构造电磁环境多等值面

遍历所有四面体，将同一阈值的等值面顶点连接成等值面，构造电磁环境多等值面；

步骤五、突出电磁环境多等值面边界轮廓

根据多等值面顶点法向量和视点的关系设定多等值面的颜色透明度，以多等值面颜色透明度来突出电磁环境多等值面边界轮廓；

步骤六、获得电磁环境多等值面半透明结果

采用快速分层伪排序融合绘制算法获得电磁环境多等值面半透明结果，并将电磁环境多等值面半透明结果融合嵌入到三维环境中。

2.根据权利要求1所述的电磁环境体数据三维可视化方法，其特征在于：所述步骤一中所述电磁环境体数据确定的具体方法是：

(1)用户初始设定纬度θ＝θ₀、经度

高度r＝r₀，计算终止的纬度设为θ₁，经度为

高度为r₁，并设定纬度、经度和高度步长为Δθ、

Δr；

(2)输入电磁设备参数、地形高程数据和气象环境数据，采用电波传播物理模型，得到电磁环境体数据值

(3)经度

递增

重复步骤(2)，直到经度

(4)纬度θ递增Δθ，重复步骤(2)，直到纬度θ＞θ₁；

(5)纬度r递增Δr，重复步骤(2)，直到高度r＞r₁。

3.根据权利要求1所述的电磁环境体数据三维可视化方法，其特征在于：所述步骤二中所述按照四面体顶点和顶点索引的方式组织成二维纹理的具体方法是：

(1)将电磁环境体数据剖分成四面体结构；

(2)按确定电磁环境体数据的顺序将电磁环境体数据所有的顶点坐标存储下来，每个顶点仅存储一次，采用RGBA32F格式的二维纹理，称为顶点函数值纹理，每个顶点对应一个纹理纹素，纹素的R、G、B、A分量对应顶点的三维坐标和电磁环境体数据值，并将该纹素的顺序号作为顶点的索引；

(3)存储每个四面体顶点的顶点索引到另一个二维纹理，称为顶点索引纹理，其格式为RGBA32UI，每个四面体对应一个纹理纹素，纹理纹素有R、G、B、A分量，每个纹理纹素分量对应四面体的一个顶点索引。

4.根据权利要求1所述的电磁环境体数据三维可视化方法，其特征在于：所述步骤三中所述根据设定的等值面阈值大小生成等值面顶点的方法是：

1、对于电磁环境体数据中其中一个四面体，判断该四面体的一条边上是否与给定阈值δ的等值面有交点，有交点即等值面的顶点。具体判断方法为：

(1)根据四面体其中一条边的端点，在顶点索引纹理中得到顶点索引；

(2)根据顶点索引在四面体顶点函数值纹理中获取该边两个端点对应的纹素，该纹素的R、G、B、A分量对应顶点的三维坐标和电磁环境体数据值，从而可以获得两个端点的顶点坐标v_a、v_b和电磁环境体数据值P_a、P_b；

(3)比较多等值面的阈值δ与端点的电磁环境体数据值，如果δ＞P_a且δ＞P_b，或者δ＜P_a且δ＜P_b则该边上不存在该阈值大小δ的等值面的顶点；否则存在等值面的顶点，插值生成等值面顶点为：

v＝(1-t)v_a+tv_b               (3)

其中， $t=\frac{\mathrm{\δ}-P_a}{P_b-P_a}$

(4)重复步骤(1)、(2)、(3)判断生成该四面体所有边上阈值δ的等值面的顶点；

2、遍历所有四面体完成步骤1，生成该阈值的等值面的所有顶点。

3、遍历所有等值面阈值，完成步骤1、2，生成所有阈值的等值面顶点。

5.根据权利要求1所述的电磁环境体数据三维可视化方法，其特征在于：所述步骤四中所述电磁环境多等值面构造方法如下：

(1)对于每个四面体，根据其六条边的排列顺序，顺序连接与边相交的同一阈值的等值面顶点，即可得到某一阈值的电磁环境等值面的一个面片；

(2)遍历多等值面的所有阈值完成步骤(1)，即可得到该四面体的电磁环境多等值面面片；

(3)对所有四面体完成步骤(1)、(2)，即可得到由多个等值面面片构成的电磁环境多等值面。

6.根据权利要求1所述的电磁环境体数据三维可视化方法，其特征在于：所述步骤五中所述多等值面颜色透明度设置的具体方法是：

(1)遍历所有多等值面的顶点，对每个顶点计算该顶点法矢量与视点方向的夹角余弦：

cosα＝n^Tv                        (4)

其中n为等值面顶点法矢量的归一化结果，v为视点与该顶点的视线方向归一化结果；

(2)设置该顶点的透明度值为：

A＝1-k|cosα|                     (5)

其中k为边界轮廓控制正数因子，取值越大透明区域越大，边界轮廓区域越突出

(3)对多层等值面模型的所有顶点完成步骤(1)、(2)。

7.根据权利要求1所述的电磁环境体数据三维可视化方法，其特征在于：所述步骤六中所述快速分层伪排序融合绘制算法如下：

(1)依据电磁环境多等值面提取阈值的大小，对每个等值面按照从内向外排序；

(2)绘制所有等值面的背面面片，实现等值面的背面面片按照从后往前的顺序融合；

(3)绘制所有等值面的正面面片，实现等值面的正面面片按照从后往前的顺序融合；

(4)融合多等值面正面面片和背面面片绘制结果，得到半透明结果；

(5)根据电磁环境多等值面的顶点法矢量，采用phong着色模型计算光照情况，为电磁环境多等值面的所有顶点着色，最终将带有光照效果的电磁环境多等值面显示在三维环境中。

Images