基于三维GIS的场数据模型研究与实践

　　地理信息系统 (geographic information sys-tem , GIS) 的灵魂是空间分析，而空间分析的基础是空间数据模型。空间数据模型是人们对现实世界地理空间实体、现象以及它们之间相互关系的认识和理解，是现实世界在计算机中的抽象与表达。
　　GIS 的概念数据模型分为 3 类：对象模型 (ob-ject model)、网络模型 ( network model) 和场模型 (field model)。比如基于平面图的点、线、面数据模型和基于连续铺盖的栅格数据模型，分别为对象模型和场模型，描述道路网络、地下管线等需要网络模型；定义了一种基于对象的模型和一种基于域的模型，分别对应对象模型和场模型；OGC (open GIS consortium) 的简单对象模型规范定义了点、线、面和表面模型，其中表面模型可以用来表达实体，也可以表达地形等 2. 5 维的连续表面。
　　但是随着三维 GIS 技术的发展，已有的对象模型和场模型在地学表达上并不充足，这些空间数据模型都缺乏对三维空间中连续、非匀质的场数据的表达能力。
　　随着三维 GIS 技术的发展以及三维 GIS 在建筑信息模型 (building information modeling, BIM)+GIS、数字化孪生、透明城市、智慧城市、城市设计等领域的广泛应用，GIS 对现实世界的表达全面扩展到空 / 天、地表和地下，并从室外走进室内，空间数据模型的二三维一体化是三维 GIS 发展的必经之路，早期基于平面图的空间数据模型已经不能满足对三维属性场的表示，比如电磁场：三维空间中 5G 信号的分布；空气的属性场：污染、温度、湿度场的表达；地质的属性场更多更复杂，如密度、孔隙度、杨氏模量、泊松比等。
　　对于三维场数据的表达，在理论研究方面，基于 TEN (tetrahedron network) 的数据模型，对其数据结构进行了简要描述；基于八叉树和 TEN 的混合数据结构，用于平衡数据量和数据表达精度之间的矛盾；基于体元的三维场数据表达方法；构建了泛在知识数据表示模型 (mubiquitous knowledgeable data representation model，UKR) 用来表达复杂的三维空间事物，包括对场数据的表达。
　　在应用研究方面，指出了三维场空间的有效表达对三维 GIS 在地学应用中的重要性，比如对地质、矿山、地下水文等真三维现象的表达与模拟；研究了三维场数据模型在地质行业的具体应用；提出了面向广义三菱柱的三维地质模型空间剖切方法，不规则四面体是其最小剖分单元；总结了地质建模软件中各类方法，包括 TIN (triangulated irregular network)、四面体等，但仍局限于地质建模领域；基于 ZTMAPGIS 软件实现了地下岩土管理系统，对地层的模拟采用多面体模型，仍未涉及三维场数据模型，无法地表达连续变化元的地质属性场。
　　综上所述，三维场数据模型在地质领域得到了一定的研究和应用，而在其他领域的应用还未见大范围展开。在通用 GIS 基础软件领域，除 2018-08SuperMap 发布了支持三维场数据模型的产品外，鲜见其他国内外主流 GIS 基础软件支持三维场数据模型的文献和报道。
　　长期以来，GIS 基础软件不支持场三维数据模型，严重制约了此类数据的应用推广，如空中 5G 信号强度场、空气污染场、水体污染场等领域的应用研究。
　　本文在 ZTMAPGIS 平台中，研究实现了不规则四面体网格 (tetrahedralized irregular mesh，TIM) 和体元栅格 (voxel grid) 两种三维场数据模型，并构建了相关的数据管理、可视化和空间分析方法，解决了三维属性场数据在 GIS 基础软件中的存储、管理、可视化和分析应用问题。
　　1 场模型从二维到三维的升维表达
　　空间数据概念模型中，网络模型描述对象之间的联通关系，包括二维和三维网络数据模型，后者可应用于三维的道路与管线拓扑关系表达与分析。对象模型采用对象的方式描述离散的空间实体，包括二维点 / 线 / 面对象、三维点 / 线 / 面对象及三维体数据模型，后者用封闭的面表达有明确边界、内部匀质的三维对象，广泛应用于 BIM 和地质体等领域。场模型则用于描述和管理连续变化的空间属性，在地形表达中常用的 TIN 和 Grid 属于场模型，但这两个模型仅能表达连续变化的表面场，常被称为 2. 5 维，本文研究的 TIM 和体元栅格则实现了对 TIN 和 Grid 的升维，用于表达和管理连续变化的三维空间属性。
　　GIS 基础软件此前实现了二维和三维网络模型、二维和三维对象模型，以及 2. 5 维持场模型，但尚不支持三维场模型。基于本文研究，ZTMAPGIS 率先实现对三维空间模型的支持，从而全覆盖支持 3 空间数据概念模型，为空间 / 天、地表、地下全要素的一体化表达、模拟与分析计算提供了技术基础，如下图所示。
　　臻图信息--GIS模型
　　场模型表达的是空间中连续的、非匀质的数据，Goodchild 将场模型的空间最基本的信息定义为元组:
　　臻图信息--GIS模型
　　式中，每个元组表示在空间位置 (x，y ) 的 n 个变量的数值，连续 x、y 表达的元组 T 的集合即可表达场数据。实际上，该元组局限于二维空间，对每个位置对应的各空间变量取唯一的数值。常见的 Grid 和 TIN 就是将地表属性，比如高程值，作为空间变量的场数据模型，严格地说，Grid 和 TIN 应属于 2. 5 维模型或表面模型。
　　以此类推，三维场数据模型表达的是三维空间中的连续的、非匀质的属性数据，三维场数据模型的基本信息可以定义为元组:
　　臻图信息--GIS模型
　　式中，元组 T 表示在三维空间位置（x，y，z）的 n 个空间变量的数值；ωi 是该位置的属性，ωi 可以是三维空间中的电磁信号、温度、湿度、污染分布等，或者是地质体的孔隙度、渗透率、含水饱和度等属性值。
　　2. 5 维空间中，栅格数据模型是对场的规则剖分表达，TIN 是对场的不规则部分表达。对应到三维空间，体元栅格为三维场的规则剖分表达，是对 Grid 的升维；不规则四面体是对三维场的不规则部分表达，是对 TIN 的升维，本文参考 TIN 将其命名为 TIM（图 2）。
　　臻图信息--GIS模型
　　2 实现难点与关键技术
　　2. 1 数据结构设计
　　在理论研究方面，国内早有学者提出采用体元 (voxel) 表达三维场数据。本文进一步明确了体元可表现为不规则和规则的多种形态，如 TIM 由不规则四面体体元构成，而体元栅格则由规则体元构成。
　　1）TIM 的最小单元是不规则四面体，如图 3（a）所示。在空间结构上，拓扑相连的不规则四面体形成的网格就是对三维空间的划分，如图 3（b）所示；在属性结构上，属性值可以依附于顶点、边或四面体本身。
　　臻图信息--GIS模型
　　TIM 可以通过带属性的三维离散点基于 3D-Delaunay 方法来构建，该方法具有以下特点:①保留原始的离散点作为不规则四面体的顶点，不会产生新的顶点；②离散点中任意五点不会在同一个球体上。因此，3D-Delaunay 剖分方法保留了原始离散点的属性值从而保证了精度，且组成四面体的三角形近似等边或等角，四面体体元的组合更逼近真实的目标实体，在插入新的顶点或者修改已有顶点时，对周围点的影响最小。
　　2）体元栅格。体元栅格的最小单元是规则的体对象，比如立方体或正六棱柱，如图 4（a）所示。通过体元对象的规则排列构成阵列，形成对三维空间的一个划分；属性值保留在体元本身，如图 4（b）所示。
　　臻图信息--GIS模型
　　元栅格是规则体元的阵列，在数据结构上比 TIM 简单。基于离散点可以插值成体元栅格，但离散点本身很可能不会被保留。在表达精度上，体元栅格弱于不规则四面体网格，对于需要表达特定特征或分布的应用，不规则四面体有一定优势。
　　2.2 可视化表达
　　与 TIN 和栅格只能表达三维表面不同，TIM 体元栅格可以表达三维空间中的任意 (x，y，z) 位置的属性分布，可以把三维场数据作为一个实体，采用不同的方式对数据内部进行直观地表达，包括:①剖切显示，对 TIM 和体元栅格，采用体绘制技术，设置不同裁剪面进行剖切显示。②分层设色，针对体元代表的属性值进行分级分类，对三维场数据的基本单元赋予不同的颜色值。具体地，TIM 采用三角网表达，颜色可以附着在顶点上，在片元阶段插值、着色。体元栅格则采用体绘制的方法，将不同属性值对应的颜色绘制到三维纹理上，实现在三维场数据在场景中的可视化。③过滤显示，渲染时根据属性值进行的过滤显示，过滤掉不满足条件的对象，实现体元栅格的过滤显示，如图 5 所示。
　　臻图信息--GIS模型
　　2.3 基于属性值的分析与计算
　　三维场数据的表达本质上是通过规则或不规则的体元来表达三维实体空间的属性场分布。因此，基于体元属性值的分析计算是实际应用的重要内容，以下主要探讨三维场数据的插值计算和针对体元栅格的统计分析和代数运算。
　　1）TIM 的构建及分析计算。 TIM 由不规则四面体构成，属性值可以依附于顶点、边或者体元本身。实际应用中，通常是通过带属性的三维观测点数据构建 TIM，属性值一般附着于不规则四面体的顶点。针对不同的应用，本文采用 3D-Delaunay 方法构建。
　　对于四面体内部任意点的属性值则需要通过插值方法计算。TIM 可以通过降维运算，获取三维空间任意剖面的属性值。对于精度要求不高的应用，可以将 TIM 数据转成体元栅格。
　　2）体元栅格的分析计算。 体元栅格由规则体元构成，属性依附于体元本身。基于体元栅格数据可以进行统计分析、代数运算、提取剖面、提取等值线等计算。
　　基于体元栅格的统计分析和代数运算，类似于二维的栅格数据，可以按基本构成单元对属性值进行各类分析计算；对于相同地理范围且空间分辨率相同的不同场数据，比如相同属性特征、不同时间点或者同一区域不同属性特征的体元栅格，可以进行叠加运算，在实际应用中可以得到特征值随时间变化的趋势或特征值之间的相关性等指标。
　　基于体元栅格数据，可以给一个剖面，提取出该剖面上的属性特征值 (图 6)，从而得到体元栅格内部属性值的分布情况。
　　臻图信息--GIS模型
　　基于三维场数据的属性特征值，可以根据体元与点空间位置关系，提取三维点、线、面及模型对象的各特征点的属性值，如图 7（a）所示；可以提取属性场中的等值线，如图 7（b）所示。
　　臻图信息--GIS模型
　　3 应用案例
　　本文基于以上关键技术，在三维基础软件中实现了 TIM 和体元栅格数据模型的存储、可视化及分析计算。以日照时长分析、5G 通讯信号强度分析、地震场的表达分析和地层属性场的表达分析为例，验证了 TIM 和体元栅格在实际应用中的可行性和有效性。
　　1）日照时长分析。 日照时长分析是指特定时间范围内，指定区域中每个位置能够被太阳光照射的总时长。通常采用规则分布的离散采样点来模拟连续区域，按 X、Y、Z 方向指定采样点的步长，分析结果为每个采样点能够获取的总日照时长。体元栅格能够很好地表达分析结果，如图 8（a）所示。其中每个立方体的属性值代表该区域的日照时长，可以对体元栅格按属性值过滤显示，查看区域内部日照时长的分布情况，如图 8（b）所示。此外，也可以对体元栅格进行剖切查看，提取任意位置的日照时长属性值。
　　臻图信息--GIS模型
　　日照分析与 BIM 数据结合可实现住宅楼日照分析应用，如图 9 所示。先基于 BIM 数据的窗户构件提取每个窗户的采样点位置，再从表达日照时长的三维场数据中提取采样点位置的属性值，即代表每个窗户的日照时长，可以进行专题表达直观地展示采光情况；通过 BIM 数据上窗户的归属关系，能够得到每户的总日照时长，为户型定价和选择提供决策依据。
　　臻图信息--GIS模型
　　2）5G 信号强度分析。 近年来，第五代移动通信系统 5G 已经成为通信业和学术界探讨的热点。全国各省市均在积极推进 5G 通信网络建设。其中，网络规划 (简称网规) 贯穿整个网络建设的全部过程，在通信网络建设中，往往会涉及许多与地理环境相关的最优方案选择问题，如选址、管线的铺设等。利用三维 GIS 技术可以实现各种渲染、分析功能，进行科学的选址和传输线路规划。网络优化 (简称网优) 是在网络安全稳定运行的基础上，通过三维 GIS 技术对网络的性能数据、测试数据等进行分析和诊断，识别出影响网络资源利用的潜在因素和问题，有利于运营商及时采取相应的技术措施，对网络进行系统地调整和控制，从而实现网络资源合理分配。
　　如图 10 所示，体元栅格将 5G 信号强度可视化，将体元栅格贴附到建筑物表面，表达区域内通信天线覆盖与信号强度，运营商可以根据可视化结果调整选址，优化网络服务，通过调整信号强弱得到不同的效果图，然后根据可视化效果选择最优的网络效果。未来，相信随着 5G 的不断发展，相关应用会得到进一步的深化。
　　臻图信息--GIS模型
　　3）地震场的表达与分析。 地震研究中的地震层析成像技术，是一种用地震数据来反演地下结构的物质属性并逐层剖析绘制其图像的技术，可以通过地震波速度模型映射地下结构构造，那么如何将地震波传播速度转换为直观的三维地质结构模型。
　　如图 11 所示，地震台网中心制作了高精度地震速度采样数据，可以将该采样数据构造成三维矩阵点，并导入 ZTMAPGIS 平台，然后插值成体元栅格数据，发布到三维 WebGL 客户端，进行三可视化表达。
　　臻图信息--GIS模型
　　在可视化层面，可以对体元栅格进行剖切，实时查看地下不同深度和纵剖面上地震波传播速度情况，如图 12（a）所示。还可以根据不同数值进行过滤和筛选，动态显示不同波速值区间的 3D 结构，基于此可以发现同一地下深度中的异常值，为地球内部结构研究、地震预报等提供参考，如图 12（b）所示。
　　臻图信息--GIS模型
　　4）地层属性场的表达。 钻孔数据是地质研究中常见的探测数据，每个钻孔点可以表达为带属性的三维点，属性可以是实际探测的土壤湿度、渗水率等。以地层渗水率为例，根据离散的钻孔数据构建 TIM，如图 13 所示，每个不规则四面体的顶点用四元组 (x,y,z,w) 表示，w 代表该点的渗水率。在实际渲染中 w 值可以对应成顶点的颜色，从而表达属性值在 TIM 中的变化。
　　在可视化层面，可以对 TIM 进行剖切，实时查看任意剖面的属性分布；在分析计算层面，可以进行插值计算，得到任意位置的渗水率。
　　臻图信息--GIS模型
　　4 结束语
　　关于三维场数据模型的理论研究较为成熟，但通用 GIS 基础软件中对三维场数据的表达和处理能力不足限制了相关理论在实际中的应用。本文在 ZTMap GIS 中实现了 TIM 和体元栅格两种三维属性场，填补了 GIS 基础软件对三维场支持的空白，有效解决了通用 GIS 平台对空中电磁信号场 (如 5G 信号)、空气与水体污染场、地震场、地下地质属性场等连续、非匀质的三维空间属性场数据进行表达、分析和应用的难题。
　　结合三维 GIS 在不同行业应用的需要，进一步探索 TIM 和体元栅格在不同领域的深层次应用，完善三维空间属性场数据的多样性表达和分布式空间分析计算能力是今后研究的重点。