micromine基本原理与方法
4.1 矿床三维模型构建方法
运用计算机技术建立矿床三维模型的研究工作从六十年代为解决浸染状矿床建模问题而采用三维块段模型以来,至今已经历了近四十年的发展。建模方法也由早期简单的方块模型,发展到如今的实体模型。下面就三维矿床模型建模方法分别进行简要的介绍。
4.1.1 线框模型
矿体的地质形态复杂多变,很难用规则的几何体来描述。它需要一种灵活、简便、快速的方法来建立矿体的不规则几何模型。目前,比较知名的采矿CAD 系统均是采用表面模型来描述矿体的几何模型。这种表面模型通常是由一系列的三角面围成的表面。如 MICL的MICROMINE 的线框模型、Maptek 的Vulcan的模型等均是表面模型。在不同的系统中表面模型的名称不同,但实质都是表面模型。由于这种表面模型在未渲染前看似由线框构成,因而在采矿CAD系统中多称为线框模型。不过,这种表面模型可以进行体积估算、表面渲染、切制剖面、快速三维显示等操作,比计算机图形学中的表面模型有所扩展。能满足矿山设计、生产中地质制图的基本要求,也是建立矿体三维实体模型的基础。
线框模型的构建主要是采用了TIN技术(不规则三角网模型)中的Voronoi图与Delaunay三角形算法。TIN是一种表示数字高程模型的方法,它既减少规则网格方法带来的数据冗余,同时在计算效率(如坡度)方面又优于纯粹基于等高线的方法。
TIN模型根据区域有限个点集将区域划分为相连的三角面网络,区域中的任
意点落在三角面的顶点、边上或三角形内。如果点不在顶点上,该点的高程值通常通过线性插值的方法得到(在边上用边的两个顶点的高程,在三角形内则用三个顶点的高程)。所以TIN是一个三维空间的分段线性模型,在整个区域内连续但不可微。
TIN的数据存储方式不仅要存储每个点的高程,还要存储其平面坐标、节点连接的拓扑关系,三角形及邻接三角形等关系。TIN模型在概念上类似于多边形网络的矢量拓扑结构,只是TIN模型不需要定义“岛”和“洞”的拓扑关系。
有许多种表达TIN拓扑结构的存储方式,一个简单的记录方式是:对于每一个三角形、边和节点都对应一个记录,三角形的记录包括三个指向它三个边的记录的指针;边的记录有四个指针字段,包括两个指向相邻三角形记录的指针和它的两个顶点的记录的指针;也可以直接对每个三角形记录其顶点和相邻三角形,见图4-1 三角网存储方式。每个节点包括三个坐标值的字段,分别存储X,X,Z坐标。这种拓扑网络结构的特点是对于给定一个三角形查询其三个顶点高程和相邻三角形所用的时间是定长的,在沿直线计算地形剖面线时具有较高的效率。当然可以在此结构的基础上增加其它变化,以提高某些特殊运算的效率,例如在顶点的记录里增加指向其关联的边的指针。
图4-1 三角网存储方式
不规则三角网数字高程由连续的三角面组成,三角面的形状和大小取决于不规则分布的测点,或节点的位置和密度。对于不规则分布的高程点,可以形式化的描述为平面的一个无序的点集P,点集中每个点p,对应于它的高程值。将该点集转成TIN,然后采用Delaunay三角剖分方法。 4.1.1.1 Delaunay三角网和Voronoi图
在介绍Delaunay三角网之前,首先需要知道Voronoi图,Voronoi图又叫泰森多边形或Dirichlet图,它是由一组由连接两邻点直线的垂直平分线组成的连续多边形组成。N个在平面上有区别的点,按照最邻近原则划分平面;每个点与它的最近邻区域相关联。Delaunay三角形是由与相邻Voronoi多边形共享一条边的相关点连接而成的三角形。Delaunay三角形的外接圆圆心是与三角形相关的Voronoi多边形的一个顶点。Voronoi三角形是Delaunay图的偶图,如图 4-2所示。
对于给定的初始点集P,有多种三角网剖分方式,其中Delaunay三角网具有以下特征:
(1)Delaunay三角网是唯一的;