本发明提供了一种识别和计算地下病害信息的方法，属于计算机图形图像处理领域。其特点在于采用不规则三角网构建地下病害三维形体，和用空间几何法估算病害体积大小。本发明通过探地雷达获取地下地质数据，对地质数据中的病害信息进行交互式拾取，再通过不规则三角网自动构建地下病害三维形体，最后利用空间几何法求取三维形体的体积。本发明具有良好的适应性、准确性和通用性。
1.一种基于不规则三角网和空间几何法的地下病害识别方法，该方法建立在地质雷达探测理论基础上，其特征在于，通过不规则三角网构建地下病害三维形体，利用空间几何法计算地下病害三维形体的体积，在识别和计算地下病害信息的过程中，依次含有以下步骤：步 骤 ( A 1 ) ：通 过 雷 达 数 据 构 造 相互 平 行的 地 质 数 据 切 片 图 像其中p为 的编号，u为地质数据切片图像的个数， 为地质数据切片图像上的地质数据点， 为 的三维坐标，为 的地质属性值，m和n为 在 上的行数和列数；
步骤(A2)：在 上构造拾取线 完成对地质病害进行交互式拾取，其中， 为 上第r个地质数据点，b＝1,2,…,t为 对应的地质病害体的编号，t为地质病害体的个数，c＝1,2，…,Nb为 在第b个地质病害体所对应的Nb个拾取线中的编号，Nb为第b个地质病害体所对应的拾取线的个数， 为 中地质数据点的个数；
步骤(A3)：构造 与 间的拾取线对 且p≠q，c＝1,2，…,Nb-1，其中q为地质数据切片图像 的编号；
步骤(A4)：基于不规则三角网构造病害体全部拾取线对间的三维形体，累计得整个病害体；
步骤(A5)：基于空间几何法求病害体各拾取线对间不规则多面体的体积，累计求得病害体的总体积；
采用基于不规则三角网构造拾取线对间不规则多面体，步骤如下：步骤(B1)：从 中选取地质数据点个数较少的一条拾取线作为基准线 且另一条拾取线作为非基准线 且
为 上的第r2个地质数据点，令s＝1，k＝1，其中，s为步骤(B2)基准边 的编号，k为步骤(B2)拾取三角形 的编号；
步骤(B2)：将以 与 为端点的线段作为基准边 根据Delaunay准则从非基准线上选取点 与基准边 构成三角形，如果k>1且jk
步骤(B3)：如果s＝1，j1≠1且 是闭合曲线,则以 为新的起始点，按由jk到jk+1的方向对 中的点重新编号，即 令s＝s+1；
步骤(B4)：如果s＝1，j1≠1，则将点 分别与点 相连，k＝k+j1-1；
步骤(B5)：如果k>1且jk-jk-1>1，则将点 分别与点 相连，k＝k+jk-jk-1-1；
步 骤 ( B 6 ) ：如 果 且 则 将 点 分 别 点 相 连 ，步骤(B7)：如果 且 则将点 分别与点 相连，步骤(B8)：如果 转步骤(B2)；否则完成拾取线对 间的不规则多面体的构造；
基于空间几何法求拾取线对间不规则多面体的体积，步骤如下：步骤(C1)：拾取线对 对应的拾取三角形集为 拾取三角形的三个顶点分别 和
步骤(C2)：找到 上相距最远的两个点记为 分别以 为圆心，以间的距离为半径在平行于地质数据切片的平面上作圆，两圆的交点记为构建四边形 使得 上的点尽可能位于 的内部，同理构造四边形使得 上的点尽可能位于 的内部；
步骤(C3)：用三角形 与 对 进行分割，用三角形与 对 进行分割；
步骤(C4)：计算 其中 表示向量 和 的混合积的绝对值，其中o为坐标原点；
步骤(C6)：计算拾取线对 对应的空间多面体的体积近似值为
一种基于不规则三角网和空间几何法的地下病害识别方法
技术领域
[0001] 本发明属于计算机图形图像处理领域，针对探地雷达数据的特点，利用不规则三角网和空间几何法精确识别和计算地下病害信息。
背景技术
[0002] 本发明对地下病害信息的识别基于探地雷达探测数据。探地雷达发射机向地下介质发射脉冲电磁波，遇到不同介质的地下目标和界面时，部分电磁波会被反射回地面并被接收天线接收。当探地雷达发射机在探测区域沿测线方向移动时，接收天线会按固定频率接收反射电磁波数据，并形成探测区域测线数据集。为提高识别精度，测线数据经过插值处理后可以生成探测区域三维空间离散数据集，通过对该数据集进行切片处理可形成相互平行的地质数据切片图像，在二维地质数据切片图像上拾取的病害信息通过不规则三角网方法可以构造出三维病害体，再结合空间几何法可以估算该病害体的体积信息。
[0003] 1908年，G.Voronoi首先在数学上限定了离散数据点的有效作用范围，并定义了二维平面上的Voronoi图，1934年，Delaunay由Voronoi图演化出更易分析和应用的Delaunay三角网。常见的Delaunay三角网的剖分准则包括空外接圆准则和最大最小角准则，同一种剖分方法可以保证构建的三角网的唯一性，并且数学上已经证明外接圆准则和最大最小角准则是等价的。
发明内容
[0004] 本发明的目的是通过对探地雷达数据按层拾取平面病害信息，利用不规则三角网构建地下病害三维形体，进而采用空间几何法估算病害体积大小。本方法具有良好的适应性、准确性和通用性。
[0005] 本发明所采用的基于不规则三角网和空间几何法的地下病害识别方法步骤如下：
[0006] 步骤 (A1) ：通过雷达数 据构造相互平行的 地质数据切片图 像(p＝1,2,…,u)，其中p为 的编号，u为地质数据切片图像的个数， 为地质数据切片图像上的地质数据点， 为 的三维坐
标， 为 的地质属性值，m和n为 在 上的行数和列数；
[0007] 步骤(A2)：在 上通过构造拾取线 (p＝1,2,…,u)对地质病害进行交互式拾取，其中， 为 上第r个地质数据点，b＝1,2,…,t为 对应的地质病害体的编号，t为地质病害体的个数，c＝1,2，…,Nb为 在第b个地质病害体所对应的Nb个拾取线中的编号，Nb为第b个地质病害体所对应的拾取线的个数， 为 中地质数据点的个数；
[0008] 步骤(A3)：构造 与 间的拾取线对 p，q∈{1，2，…，u}且p≠q，c＝1，2，…，Nb-1，其中q为地质数据切片图像 的编号；
[0009] 步骤(A4)：基于不规则三角网构造病害体全部拾取线对间的不规则多面体，累计得整个病害体；
[0010] 步骤(A5)：基于空间几何法求病害体各拾取线对间不规则多面体的体积，累计求得病害体的总体积。
[0011] 本发明所采用的基于不规则三角网构造拾取线对间不规则多面体的步骤如下：
[0012] 步骤(B1)：从 中选取地质数据点个数较少的一条拾取线作为基准线 且另一条拾取线作为非基准线 且其中 为 上的第r1个地质数据
点， 为 上的第r2个地质数据点，令s＝1，k＝1，其中，s为步骤(B2)基准边 的编号，k为步骤(B2)拾取三角形 的编号；
[0013] 步骤(B2)：将以 和 为端点的线段作为基准边 根据Delaunay准则从非基准线上选取点 与基准边 构成三角形，如果k＞1且jk＜jk-1则 令构造拾取三角形 令k＝k+1；
[0014] 步骤(B3)：如果s＝1，j1≠1且 是闭合曲线,则以 为新的起始点，按由jk到jk+1的方向对 中的点重新编号，即 令s＝s+1；
[0015] 步骤(B4)：如果s＝1，j1≠1，则将点 分别与点 e∈{1,2,…,j1-1}相连，k＝k+j1-1；
[0016] 步骤(B5)：如果k＞1且jk-jk-1＞1，则将点 分别与点 相连，e∈{jk-1+1,jk-1+2,…,jk-1}，k＝k+jk-jk-1-1；
[0017] 步 骤 ( B 6) ：如 果 且 则 将 点 分 别 点 相 连 ，
[0018] 步骤 (B7) ：如果 且 则将点 分别与点 相连 ，
[0019] 步骤(B8)：如果 转步骤(B2)；否则完成拾取线对 间的不规则多面体的构造。
[0020] 本发明采用的基于空间几何法求拾取线对间不规则多面体体积的步骤如下：
[0021] 步骤(C1)：拾取线对 对应的拾取三角形集为 拾取三角形 的三个顶点分别 和
[0022] 步骤(C2)：找到 上相距最远的两个点记为 分别以 为圆心，以间的距离为半径在平行于地质数据切片的平面上作圆，两圆的交点记为
构建四边形 使得 上的点尽可能位于 的内部，同理构造四边形
使得 上的点尽可能位于 的内部；
[0023] 步骤(C3)：用三角形 与 对 进行分割，用三角形 与 对 进行分割；
[0024] 步骤(C4)：计算 其中 表示向量 和 的混合积的绝对值(o为坐标原点)；
[0025] 步骤(C5)：计算 和
[0026] 步骤(C6)：计算拾取线对 对应的空间多面体的体积近似值为
[0027] 本发明具有以下优点：
[0028] 1、本发明能够识别任意形状的三维病害形体，适应性强。
[0029] 2、本发明能够自动计算三维病害形体整体或部分体积，鲁棒性和准确性强。
[0030] 3、本发明适合各种探地雷达数据的处理，通用性强。
附图说明
[0031] 图1本发明不规则三角网构造拾取线对间不规则多面体流程图
[0032] 图2本发明空间几何法求拾取线对间不规则多面体的体积流程图
[0033] 图3本发明简明流程图
具体实施方式
[0034] 本发明对探地雷达获取的地质数据中的病害信息进行交互式拾取，通过不规则三角网构造三维病害体，利用空间几何法求不规则多面体的体积，进而得到地下地质灾害的信息。
[0035] 不规则三角网构造拾各取线对间不规则多面体流程如下：
[0036] (1)如图1所示，在拾取线对 选包含地质数据点个数较少的拾取线作为基准线另一条拾取线作为非基准线
[0037] (2)把 中相邻的两个点间的连线作为拾取边 按照Delaunay准则从 中选取一点 与 两端点相连构造拾取三角形 最终得到 间的三维病害体。
[0038] 空间几何法求拾取线对间不规则多面体的体积流程如下：
[0039] (1)如图2所示，求拾取线 和 的对应的四边形 与 并用三角形分别对两外界矩形进行分割，得到三角形
[0040] (2)利用三角形 和拾取线对 对应的拾取三角形集计算 对应的空间几何体的体积。
