基于Skyline的地质模型属性识别算法的研究与实现

2015-12-11史少维姚鹏君曹宏涛林巨超

测绘通报 2015年2期

史少维，姚鹏君，曹宏涛，林巨超

(1.北京东方新星石化工程股份有限公司(BNEC)，北京100070;2.北京建筑大学测绘与城市空间信息学院，北京100070)

一、前言

在地下地质相关的工程项目建设过程中，为便于工程管理，需要建立三维模型系统，实现地质模型数据浏览显示的功能是必不可少的，进而要求对特定区域的地质岩性情况进行实时的检测查询，以对工程地质情况进行分析与决策。

目前，随着三维激光扫描技术的发展，已存在利用三维影像及点云数据应用到地质建模及空间识别的领域中;在软件平台方面，QuantyView平台是比较成熟的三维地质建模系统，其利用地质钻孔数据、平面和断面剖图等多源数据进行建模，以及对地质模型进行分析管理等［1］。

本文主要针对地质模型的岩石属性检测的问题，基于空间几何拓扑关系，在Skyline的三维平台中，利用其强大的空间分析算法接口进行二次开发，实现了地质岩石模型空间识别的功能开发。

二、空间几何拓扑关系

空间几何拓扑关系反映了空间目标的逻辑结构，对空间目标查询、分析和空间目标重建具有重要意义。空间几何数据有点、线、面、体4类，它们之间有相离、相等、相接、相交、包含于、包含、交叠、覆盖、被覆盖、进入、穿越和被穿越12种基本空间拓扑关系。其中穿越与被穿越关系在城市GIS、矿山GIS中进行空间曲面分析时运用较为常见，具有重要的意义。

面与面的空间相交穿越关系识别规则为:判断两个面对象是否相交，可以先判断它们的最小包围盒是否相交，如果不相交，则可以直接判断这两个面对象不相交;如果最小包围盒相交，可以通过依次判断面对象与构成另一个面对象的每条线段是否相交得到。其空间面间相交关系如图1所示。

图1 面与面之间相交穿越检测示意图

本文基于空间几何体之间拓扑关系，提出了利用空间面与面相交穿越关系运算来识别目标岩石区域面，进而显示属性信息的方法。

三、Skyline地质模型识别方案

在地下地质工程系统中实现对地下巷道岩石属性的查询模块，当鼠标单击巷道时，能精确显示岩石的属性。通过3ds Max或Skyline builder建立的巷道模型结合运用CorelDRAW等图形软件处理巷道的地质纹理素描图，生成具有逼真效果的地质模型。不过在Skyline中不识别具有地质纹理识别相应模块，而是将三维模型作为一个整体存储显示，因此不能查询不同区域模型表面的岩石属性。

1.总体思路方案

首先，在Skyline中地质模型是作为一个整体的，因此不能查询不同区域的岩石属性。由此设计将不同的地质岩石在模型的基础上进行碎片化，形成岩石曲面体数据图层。

其次，由于实际地质情况的复杂性及区域不规则的特点，鼠标点触时准确选择岩石曲面体变得困难，往往视觉直观选择了目标岩石区域，但是在三维空间中，会出现假选或选空的情况。为了应对这种情况，选择设计空间查询体进行空间几何面的相交来侦测目标岩石曲面体。

最后，基于空间面相交穿越分析，利用Skyline的接口实现岩石属性识别功能的开发。程序上体现为:获取岩性区域目标点，在目标点生成查询体，即以目标点为中心生成一个包含3个相互垂直正方形面的立方面体，再利用查询体与查询目标体曲面层数据进行空间分析计算，判断是否有曲面体被侦测到，返回结果进而显示曲面体的岩石属性值。

根据总体思路方案，基于Skyline的模型识别算法流程设计如图2所示。

图2 模型识别算法流程

2.空间查询体设计

正立方体面(PolygonGeometry)定义为:xyz轴上设置一定的距离间隔生成3个正方形面，每一个正方形面定义为查询面。每个查询面都有以下特征:以鼠标点坐标为中心，分别垂直于x、y、z轴;距离中心间隔相等。间隔设置要结合FeatureLayer中曲面多边形的大小来设置，过大则侦测范围大，则返回的结果数量多且不准确，过小则侦测范围小会出现没有结果返回的情况，一般采用经验法进行测试。测试的结果表明，间隔设置为最小的曲面体的一半左右即可，效果快且准确。空间查询体设计如图3所示。

3.岩石数据图层制造

岩石数据图层(FeatureLayer)(如图4所示)将不同的地质岩石三维模型在3ds Max中切成空间三维几何体集，每一个三维几何曲面代表着一个不同的地质岩石区域，转换数据格式导入到ArcSence中或在Skyline中进行属性编辑。这一步主要是按照对地质模型的勘察报告逐一进行岩性描述编辑，在三维空间多边面的岩性属性字段中进行赋值。在算法实现时，对其岩石属性字段进行结果显示即可。具体步骤如下:

1)根据巷道纹理贴图在3ds Max里进行碎片剪切的工作，将岩石模型的属性图层抽离出来生成空间多边形几何体，添加定义岩石属性字段，导出SHP等格式文件。

图3 查询体结构设计示意图

图4 岩石数据图层示意图

2)导入到ArcSence中结合巷道勘察素描图数据对每一个碎片的岩石属性字段进行描述编辑保存。在进行属性编辑的过程时，将对应勘察岩石报告整理属性到Excel表中，打开SHP图层的属性表进行批量导入。

3)保存成Skyline的格式数据作为一个图层加载到Skyline中。将岩性FL图层以SHP格式准确地添加到FLY数据的信息树中，在具体的操作过程中要注意FeatureLayer高程、精度及字段选择的设置，保证岩石属性图层与原有的地质模型相吻合。这样保证对不同纹理选择目标时准确对应Feature Layer的岩石曲面，进而通过算法空间分析显示相应岩石属性。

选中FeatureLayers中的一个目标面进行岩性属性描述字段进行编辑。图5中显示的立方体表示选中的目标面Polygon，在勘察素描图中的相应区域为:花岗岩C类Q值为1.67。