三维虚拟地球中实时视域分析方法研究.doc

三维虚拟地球中实时视域分析方法研究 1引言 基于网络环境的三维虚拟地球软件推动了地理信息大众化发展。但在专业GIS丿应用方 而，由于缺乏三维空间分析功能，三维虚拟地球软件往往还不具备实用性。视域分析作为 重要的空间分析方法，广泛应用在旅游屮的风景评价，房地产屮视线遮挡判断，以及通讯 屮的信号覆盖，或军事上的火力覆盖等多方面［1］。在传统三维GIS软件屮，视域分析基于 通视分析进行，即通过判断给定对彖与一纟R目标点之间的通视性，从而获得指定对象的视 域可见范围。主要通视分析算法有以下4种：JAXUS算法，DWACS(Dynamic Tactical Simulation)算法，ModSAF算法和Bresenham算法［2］。以上方法部是基于平面基准，未考 虑地球曲率的影响。而基于三维虚拟地球的空间视域分析，地球曲率对分析结果影响很 大。例如，美军“宙斯盾”系统的核心SPY-1相控阵雷达，最大探测距离超过300 T?米， 由于地球曲率的影响，对低空掠海10米飞行目标的探测距离却只有约31千米［3］ o 对此，木文主要探讨了三维虚拟地球屮基于球面基准下的视域分析方法，对JANUS算 法进行改进，对全球地形数据采用金字塔分层分块数据组织基础上，顾及地球曲率条件， 实现大范FBI视域分析，同时基于NSCGlobe平台将分析结果进行了可视化。 2三维虚拟地球中实时视域分析方法 2.1视域分析中地形数据组织 在三维虚拟地球屮，构建三维场景的地形数据具有海量特性。因此，在进行视域分析 时，无法将大范围地形数据一次性载入内存进行计算分析。同时，对于视域分析，根据应 用特点灵活采用不同尺度的地形，既不影响分析精度，还可以提高分析处理的效率。例 如，对于地面火炮系统，采用较高分辨率的地形有利于精确打击。而对于那些与细节模型 关系不大的应用屮，高分辨率模型反而会影响分析处理效率。如对于离地10000米以上作 战的飞机，分析其火力覆盖的视域时就可以采用较低分辨率的地形，以便火力覆盖的实时 计算。 因此，根据视域分析的多样性，应该灵活采用多尺度地形进行分析计算。这是由于多 尺度地形避免了信息兀余，并且易于支持任何细节层次的地形表达［3］。基于上述考虑，综 合已有全球空间数据组织方法，在三维虚拟地球屮进行视域分析，木文以金字塔结构分层 分块数据组织方法对三维虚拟地球屮海址存储的数字高程模型DEM (Digital Elevatio n Model）进行组织，其核心是而向全球范即I，采用Ottoson提出的用椭球四叉树结构建立全 球等经纬度的空间数据索引［4］,分别建立东西半球两个层次四叉树结构，每个层次四叉树 结构屮四个节点记录节点的空间范|韦|,同时存储指向该范用的地形数据集指针，如图1。 -iwroMI8Q -iwr oM I8Q ?9（户 0° 90 180° 基于上述数据纟R织方法，在进行全球范帀视域分析时，根据视域分析特点，采用介适 尺度的地形数据层。确定地形数据层示，根据视域覆盖，计算需要的地形块索引。然后根 据索引调度数据，逐块进行视域分析，建立多尺度分块的地形数据的视域分析计算分析。 2. 2全球基准下的视域分析方法 在传统基于平面基准的视域分析方法屮，无需考虑地球曲率的影响。在全球基准下进 行大范I韦］视域分析，受到地球曲率影响明显。若将上述方法运用在全球大范围条件下，将 导致球面下不可见的目标点判断成可见点,如图2、图3。 图2球面下A、B不能通视 B图3展成平面示 B 图3展成平面示A、B可以通视 针对上述问题，木文对JANUS算法进行改进，在三维球面空间屮进行视域分析。主要 流程如下： ?在三维虚拟地球屮实时确定观察点经纬度B, L和高程Z,并转换为以地心为原点 的直角坐标系统中V(xO, yO, zO); ?根据观察点经纬度坐标(lat, Ion)和分析半径R(分析半径),计算需要分析的DEM 块索引编码，计算方法如公式仃)： DEM块半径:2?(DEM) = 72(分析半径)/N DE碱层号 LevelJD = [los(；?(DEM)/2?(Lev4_Zero))]+l DE碱编号 BlockJD = RoyvjCol (1) Row = (9 0. (Hlat)! ti 1 esize , Col = (180.(Hlon)! til esize tilesize = 2*7?(Le\rel__Zero) * 0.5^ z 地形块索引编码由Level_ID+ BlockJD两部分纟fl成。式屮，N为经验值，依显不屏幕 分辨率而定，取值范围为3—6, R(Level_Zero)为0层(即最上层)地形块半径,tilesize为 第Level.ID层地形块经纬度跨度。 ?根据地形块索引编码，将地形块逐块载入内存进行视域分析；