GoogleEarth的出现给了人们全新的感受,数据动态装载、地形渐进描绘、多重细节层次(LOD: Levels of Detail)和虚拟现实表现等已成为虚拟地球的典型技术特征。构建数字地球有这么几个问题:

  • 全球角度空间数据模型的构造;
  • 海量地理空间数据(地形、影像等)的组织、管理和调度;
  • 空间数据的三维可视化。

    全球角度的空间数据模型

    投影理论解决了球面或椭球面不可展的难题(尽管不完美),从而建立了二维平面与实际球面三维空间的对应关系,球面网格理论则构建了全球空间数据高效访问的应用分析模型。

    (一)、球面网格

    球面网格的基本原理是将球面(参考椭球面)递归剖分为面积、形状近似相等且具有多分辨率层次结构的网格单元,采用单元对应的地址编码代替地理坐标在球面上进行操作。
    全球网格主要分为经纬线网格和多面体网格。
    经纬度网格符合人们的思维习惯,计算简单直观且与现有各类数据转换比较方便,容易进行纹理坐标的计算和坐标转换(地理坐标到三维直角坐标),在添加地理数据时,所做的工作也最少。缺点是网格单元形状不均匀,接近极点的时候,网格扭曲,数据浪费严重。如图所示
    数字地球实现_第1张图片
    多面体网格多限于学术研究,少有工程化的实现。其最大优点是每个网格单元扭曲很小,没有太大的数据浪费。这个特性有利于三维球面绘制。但缺点是在球面绘制时纹理映射比较复杂,另外在与地理坐标的转换上也比较麻烦。
    数字地球实现_第2张图片

    (二)、基于经纬度的球面网格

    基于经纬度的球面本身就隐含了一个矩形网格系统,为了改进在高纬度地区单元变形严重的缺陷,GeoFusion(2005)将全球划分为6个区域,其中4个处于南北纬45度之间,采用矩形网格,两个处于高纬度地区,采用三角形网格,避免了高纬度地区使用矩形网格带来的变形,但数学模型较复杂。
    目前主流的虚拟地球浏览器均采用基于经纬度的球面网格系统,其中GoogleEarth、NASA WorldWind采用统一的矩形网格,ArcGlobe则采用GeoFusion的全球网格技术。
    数字地球实现_第3张图片
    (world wind网格)
    数字地球实现_第4张图片
    (geofgeofushion网格)

    (三)、全球地形网格

    全球地形是虚拟地球的可视化基础。地形网格可以分为规则网格(RSG)和不规则三角网格(TIN)两类。
    1.规则网格
    2.不规则三角网格

    空间数据组织与管理

    虚拟地球的数据主要包括地形数据、二维矢量数据和模型数据。

    (一)、地形信息的管理

    地形数据又分为高程数据和影像数据,他们多采用分级组织、分块存储。
    1.多分辨率金字塔
    多分辨率金字塔就是用一系列网格表示统一地区的地形,相邻上下层次之间的网格的采样精度一般为1:4,因此在大范围地形显示时,采用较低精度的数据能大幅度缩短数据访问时间。
    如果在每层网格上采用“分块策略”,则数据访问时间可进一步缩短,其基本思想就是将各层网格都细分为大小相等的矩形“瓦片”(Tile),这种数据结构也称为“瓦片金字塔”。其优点是一次读/写操作就能访问一个瓦片,并且经过甄别可以只访问所需的瓦片;缺点是数据冗余更多。
    数字地球实现_第5张图片
    2.瓦片四叉树
    多分辨率的TIN不具备简单的跨越式空间结构,在执行诸如视景体裁减、碰撞检测等操作时显得更为复杂,并且对TIN的压缩算法也比较复杂,不利于数据的快速访问。
    采用“分层分块”策略的瓦片金字塔不但大大缩短了数据访问时间而且能够表示多分辨率数据,是海量地形数据组织的理想结构,通过与树结构的结合,具备了快速数据索引的功能,从而形成了大范围空间数据组织方式——瓦片四叉树。采用四叉树描述的瓦拓扑关系简单,其索引只涉及加、减、加倍和减半操作,可以通过效率较高的位运算实现。
    数字地球实现_第6张图片
    基于经纬的球面网格的瓦片的数据采用经纬度加高程描述,符合人们的习惯,也比较直观,而在显示时,则要转换为空间直角坐标系。
    目前流行的虚拟地球软件的数据组织方式,均采用瓦片四叉树结构。

    (二)、模型数据的管理

    在实现模型数据的管理时一般有两种方式,一是集成式存储,这种方式多用于存储自动建模或精细建模所生成的模型文件;另一种是以三维空间数据模型实现的模型数据(二维矢量数据是三维模型数据的一个子集),多以数据库的方式管理,其中三维空间数据模型由包含点、线、面、体这四类空间对象的空间实体构成,其索引可以采用三维R树来实现 。

    三维可视化

    在与地形、影像等栅格数据相关的各类应用中,大范围三维地形实时浏览对数据组织和呈现方式的要求较高,主要体现在两方面:海量数据处理和球面模型建立。

    (一)、可视化关键技术

    在计算机交互式图形处理中,实时动画往往要求每秒25~30帧的图形刷新频率,对数据调度机制和图形绘制策略都提出了很高要求。
    为了在保证地形外观平滑的前提下尽量减少需要渲染的数据,有多种图形生成加速方法被提出,数据裁减、多细节层次技术(Level of Detail ,LOD)、异步数据加载、数据缓存、请求预测已经被广泛用于各种大规模地形场景的快速绘制算法中。
    1.数据裁减
    数据裁减主要是视景体裁减(View Frustum Culling),视景体裁减的核心是计算视场的锥体裁减范围,即由视场角定义的上下左右四个面(类似照相机的视角)和由投影矩阵定义的远近剪切平面(类似照相机的景深),而后看三维实体(对于球面瓦片,可用其包容圆柱体)是否与视景体的相交或被视景体所包含。
    数字地球实现_第7张图片
    2.多细节层次技术
    LOD技术是在不影响画面视觉效果的前提条件下,通过逐次简化景物的表面细节来减少场景的几何复杂性,从而提高绘制算法的效率,其所遵循的原则就是能够模拟“越远越模糊,越近越清楚”这一视觉现象,通常利用视角高度(Eyealt)与瓦片尺寸的比来判定当前显示的瓦片是否需要分解为更高分辨率的子瓦片。
    数字地球实现_第8张图片
    3.异步数据加载
    异步数据加载就是通过不同的线程加载瓦片数据,线程之间相互独立而又彼此关联,其工作原理是:主线程通过LOD算法与数据裁减计算出所需瓦片,如果瓦片不在数据缓存中,则将其添加到请求队列,当数据加载完毕,通知主线程瓦片可以使用,刷新视图。为了在数据加载时不影响显示的流畅性,主线程会先采用其上一级直至顶级的数据来绘制。异步数据加载还需剔除请求队列中因场景变换已不再需要加载的瓦片,避免资源的浪费。
    4.数据缓存
    数据缓存有磁盘缓存和内存缓存,磁盘缓存用于将远程数据下载到本地硬盘,以提高数据加载速度,一般情况下磁盘缓存的大小是被限定的。
    内存缓存也叫高速缓存,是已经加载了的数据,可直接用于场景显示。高速缓存的大小设有上下限,用以控制缓存的瓦片数量不致太少也不要太多,缓存数据的置换算法可以用最近最少使用算法(Least Recent Use, LRU),LRU的主要出发点就是,如果某个瓦片被访问了,则它可能马上还会被访问。
    5.请求预测
    数据缓存机制可以提高系统的性能,当地形场景连续运动时,如果应用程序能够提前预测系统即将需要的数据,并能提前把这些数据读到缓存区,那么当应用程序真正需要这些数据的时候,便可以直接从缓存区读取,从而省去了从数据源读取这些数据的时间。在海量地形漫游时,数据调度和场景绘制需要同时进行,因此,可以设立一个独立的线程来进行数据的预测调度。

    (二)、地形的三维可视化

    1.地形网格的构造
    地形数据多采用地理坐标(如经纬度)加高程信息的地形数据,这需要转换为对计算机而言表达简单,处理效率较高的空间直角坐标。
    数字地球实现_第9张图片
    大地坐标系是统一的,能准确、唯一地描述地球上任一点的位置,并能动态适应历史的、现代的越来越精密的各种参考椭球体,是虚拟地球最为适当的坐标系。
    数字地球实现_第10张图片
    三维地形模型的几何构造要素多为三角形,即每个瓦片根据地形点阵数据,再构成三角形网格,这一过程是动态实时进行的,构造好的地形网格被存入地形高速缓存以备调用。
    数字地球实现_第11张图片
    (地形网格)
    数字地球实现_第12张图片
    (geofusion地形网格)
    2.纹理映射
    纹理映射能够有效增强图像的真实感,在地形可视化中,卫星影像是最能提高地形真实性纹理。目前,几乎所有支持纹理映射的图形处理硬件都在其硬件系统中设置了纹理内存来缓存纹理数据以提高纹理映射的效率。纹理内存的容量有限,图形处理硬件可能无法缓存全部纹理数据,每次纹理映射时都需从主存中传输部分纹理数据,因此也需要设置纹理高速缓存中以备调用。

    (三)、模型的三维可视化

    模型数据也可以用瓦片为单位建立高速缓存,来提高系统性能。
    在三维模型的展示上,点状对象(如树、路灯、管线点等)常采用三维模型符号来表示;线状对象(如管线)也是先建立三维线状符号(如圆柱体),获取线状对象的走向、长度、半径、高度方面的信息后,经过缩放、旋转、平移操作实现线状对象的绘制;面状对象绘制需建立模型并配以纹理贴图。

    (四)、三维注记

    注记是虚拟地球中一个不可缺少的要素,有了注记,才能保证对地理信息的正确认识和方便地进行查询分析。汉字的显示与单字节文字的显示不一样,在OpenGL中是以轮廓文本的方式输出的。

    测试程序设计与实现

    (一)、三维平台简介

    1.三维引擎
    主要有OpenGL和Driect3D。
    2.现有的虚拟地球产品
    基于OpenGL的虚拟地球产品主要有:NASA WorldWind Java版、OssimPlanet(以开源的Open Scene Graph为基础),ArcGlobe;基于Driect3D的虚拟地球产品主要有:NASA WorldWind C#版、Skyline。GoogleEarth则支持OpenGL、Driect3D两种三维引擎渲染方式。
    上述平台都支持二次开发,国内的一些虚拟地球产品多以开源的平台为基础进行二次研发。
    3.测试程序开发平台的选择
    作为一个测试程序,在三维平台上主要参考开源的NASA WorldWind Java版和基于OSG的OssimPlanet。
    NASAWorldWind Java版基于JOGL开发,是一款架构优秀的虚拟地球开源软件,其数据模型针对椭球球面网格而建,而在城市三维模型展示以及矢量数据处理上目前支持不足,且在渲染时需要由Java到C语言的一次内存拷贝,性能上会打折扣。
    数字地球实现_第13张图片
    (nasa world wind)
    OssimPlanet继承了OSG的场景数据管理结构且实现了与NASA WorldWind类似的全球网格模型,OSG基于C++语言,大量运用了C++模版技术,因而开发难度较大,在数据结构中主要运用树结构组织数据,用C++智能指针来实现数据的共享,用仿函数类(重载()操作符的C++类)实现事件映射。
    数字地球实现_第14张图片

    (二)、测试数据组织

    在全球高程和影像数据瓦片金字塔构建上,各虚拟地球产品基本一致,不同的是最上层(0层)的瓦片大小(Level Zero Tile Size),第0层可以是一个瓦片,经度跨360度,纬度跨180度,也可以是多个瓦片,NASA World Wind的第0层影像瓦片有10×5个,经纬度跨36度,其高程数据第0层瓦片为18×9个,经纬度跨20度。
    数字地球实现_第15张图片
    1.瓦片金字塔模型构建
    在构建地形金字塔时,首先把原始地形数据作为金字塔的底层(图5-4),如图第2 层,并对其进行分块,形成第2 层瓦片矩阵。在第2 层的基础上,按每2×2个像素合成为一个像素的方法生成第1 层,并对其进行分块,形成第1 层瓦片矩阵。如此下去,构成整个瓦片金字塔。
    数字地球实现_第16张图片
    对于全球的高程和瓦片数据,分多少级是合适的呢?三维引擎的透视原理是模拟相机在标准视角下拍摄物体的情形,相机的标准视角约为25度,此时的视角高度与瓦片尺寸的比值为2,根据这一原理,如果原始图像是1米精度,视角高度为500米的卫星影像,则赤道上每一瓦片的宽度应为250米,对于周长为40000公里的地球,如果第0层是1个瓦片的话,要建0至17总计为18层的瓦片金字塔。
    2.瓦片四叉树
    四叉树是一种每个非叶子节点最多只有四个分支的树型结构,也是一种层次数据结构,其特性是能够实现空间递归分解。下图是瓦片金字塔模型的四叉树结构示意图,其中矩形符号代表叶子节点,圆形符号代表非叶子节点。
    wpscs
    (瓦片四叉树)
    四叉树的层与瓦片的层对应,四叉树的节点编号是瓦片所在层的层号加上瓦片的网格坐标。因此,瓦片的拓扑关系包括同一层内邻接关系和上下层之间的父与子关系两个方面。邻接关系有边邻接和角邻接,上下层间有向上继承和向下分解。
    数字地球实现_第17张图片
    3.瓦片金字塔的编程实现
    在瓦片金字塔的编程实现上,可以有两种方式,一种是建立一个四叉树数据结构,保存瓦片上下层之间的拓扑关系,在一个节点中保存父节点和四个子节点的指针,另一种是用散列映射(hash_map)数据结构,其中瓦片的索引作为键,瓦片的智能指针作为值。
    四叉树数据结构在调入一个瓦片时必须先生成其父瓦片指针,在删除一个瓦片时必须先删除其所有子瓦片指针,不利于LRU算法的有效实施,而散列映射数据结构则要灵活得多。

    (三)、数据的调度与可视化

    地形网格的调度过程和纹理影像的调度过程是一样的,一先一后,假设地形数据与影像数据是一起加载,则数据调度与可视化流程表示为图所示。
    数字地球实现_第18张图片
    系统运行效果如图所示。
    数字地球实现_第19张图片