Delaunay三角网在土方量计算中的应用研究.pdf

中国测绘学会第十七届海洋测绘综合性学术研讨会论文集 97 Delaunay三角网在土方量计算中的应用 汪连贺，董江 (天津海事局海测大队，天津300220) 摘要：不规则三角网在各个行业和领域中的应用越来越广泛，简要介绍了建立不规则三角网的主要过程以及在此基础上 进行土方量计算的主要算法。 关键词：不规则三角网；-I-_方量计算；面积公式 L()P算法优化过程见图1。 1引言 3建立初始Delaunay三角网 在港口航道施工中经常要计算挖泥土方量，某 些施工单位只是简单的计算施工区域平均水深，然 3．1寻找离散数据点的凸包多边形 后用此值作为深度，再根据设计水深和施工区域的 为了加快寻找凸包，以最小的z、Y和最大的 面积计算土方量。无疑，这种计算方法与实际情况 z、Y为坐标的数据点的连线为界(以下称为A线)， 相差太大，不能准确地反映施工工作量。笔者根据 把所有离散数据点分为两部分，并标记每个数据点 以往的编程经验，在实现Delaunay三角网的基础上相对于A线左右区的标识，判别公式如下： 计算挖泥土方量。 F(z，Y)=(Y—Ax—B)×(3：2一z1)(1) 式中A：羔2二羔!R：羔!釜三2二羔2垒苎 2建立Delaunay三角网主要过程 ’、’”一 Z2一Z1’一 Z2一Zl。 如果F(z，Y)0，则点位于A线左区；如果F Voronoi多边形和Delaunay三角网的定义许多 书籍或者文章里都有详细地介绍，在这里就不再重 点位于A线右区。作者对原来的判别公式进行了 复，这里简要介绍一下逐点插入法的基本步骤： (1)定义一个包含所有数据点的初始多边形； 改进，在原来的公式基础上乘以(z2一z。)，则可以 (2)在初始多边形中建立初始三角网，然后迭判断点在A线的上侧还是下侧。 寻找凸包多边形的点可以从最小坐标数据点开 代以下步骤，直至所有数据点被处理； 始，从左向右寻找凸包下边界，这时只需在负区的数 (3)插入一个数据点P，在三角网中找出包含 P点的三角形，把P点与三角形的三个顶点相连， 据点集合中去找，到达最大X后从右向左寻找凸包 生成新的三角形； 上边界，这时只需在正区的数据点集合中去找。在 ．，”～＼、 建立凸包的过程中，应该保证正在寻找下一点与当 前点之间的斜率最小，这样可以保证选取的点都是 ／ 、 ，／一、、 边界上的点。由最小的z到达最大的z后，再由最 大的X到达最小的z，找到了所有的边界点后，顺序 相连所有数据点，即完成了凸包的建立，并建立了三 角网的初始边的集合。斜率计算公式如下： K=必 (2) ＼一～—／／／ 3：2一Z1 图1 LOP算法优化过程图