德劳内三角网在多波束测深中的应用.doc

德劳内三角剖分算法在多波束测深中的应用 Gert Brouns1; Alain De Wulf 2; and Denis Constales3 摘要：德劳内三角剖分对于地形面的曲面表示而言，是一种极具研究价值的模型。本文在简要的理论说明之后，分析列举了六种建立德劳内三角网的基本算法，并研究了德劳内三角网可以用于多波束测深数据处理的一些特性。除此之外，文章还着重介绍了两种理论：“分而治之”的构造算法和“增量法”，前者的计算速度之快使其成为多波束数据初始三角化时的首选算法，而与“增量法”相比，它的运行效率也证明了这一点。“分而治之”的构造算法中的融合步骤在代替已布设三角网区域中的由新数据生成的既有三角形方面也发挥着巨大的作用。另一方面，虽然“增量法”的构造效率不高，但应用此方法，使得在既有三角网中内插顶点和编辑数据变得更加容易。 DOI：10.1061/（ASCE）733-9453~（2003）129:2~（79） CE数据库主题词：三角剖分；算法；流量过程线；数学模型；地形 引言 获取海底地形的精确模型是疏浚工作中需要重点考虑的一个问题。一方面，目前的海道测量工具，特别是多波束测深仪，对海底地形的采样量非常大。举个典型的例子，以Altas Fansweep 20（目前分辨率较高的一套系统）为例，进行多波束测深时的更新速率为每秒1-12次，每个测幅中的波束脚印个数为200-800个，如此，每小时就可以采样50万至3400万个点。另一方面，为生成精确的地形模型，如此庞大的数据量将按照用户的特定需要来进行处理。 以下几点就在处理疏浚用多波束数据时显得尤为重要： ? 体积计算的最终结果越精确越好且要保证在限差之内； ? 原始测点经常被要求包含在求解体积的数字地面模型之中； ? 项目的目的是生成实时的或近期的数字地面模型来对在过去已经实现的拖曳进行实时的问责。 在数字地面模型中保存原始观测点或反映原始的观测值是经常会遇到的一种情况，常用的解决方法就是构建不规则三角网（TIN）。另一种常用的数字地面模型构造方法就是用规则格网覆盖采样点，但是这种方法有两个明显的缺点。第一，由于测量结果不是在规则格网中获得的，所以要在规则格网中为每个采样点分配一个独立的位置是极为困难的，而成果取决于测量船自身的测量系统（是否为等间隔量测）和姿态(横偏，纵偏和艏偏)。其次，格网值不能反映真是的测量值，因为格网化或者要在测量区间大于格网范围处设值内插，或者要在格网值大于测量值处删除原有信息并重采样。 还有其他的模型，要通过观测获得数据，并以模型通过已知点，例如MQ-B方法（多次曲面插值法，Cressie 1993），各种形式的克里金插值法（Cressie 1993），或者用基本函数拟合采样点以便描绘样条曲线。但是通常情况下，在应用这些函数的时候，计算出的数值被分配个规则格网中的每个单元，然后这些格网单元就将独立的代表地表面而不是进一步的参照于衍生格网的数学模型。这种对于数据集的操作简便快捷却失去了实施实际观测过程的优势。 德劳内三角剖分不仅在处理庞杂数据集方面十分高效，而且是为数不多的支持增量法内插的方法之一。 德劳内三角剖分：理论概述 图1：（a）点集；（b）自由边；（c）三角表面 一般的，将一个点集V之中的点三角化会生成一个图形（见图1），这个图形将点集V之中的点用相互独立的自由边连接起来。所谓的相互独立，是指这些边没有相交，而且以点集V之中的点作为顶点。这样一来，整个平面就被剖分成许多互不相交的自由三角形，我们称之为三角面片。面片，边，点都被凸多边形的边界所围起来，多边形和它的外围区域将整个平面完全填充。 德劳内三角剖分是一种特殊的三角化方法，其准则就是，完成三角化之后，任意一个三角形的外接圆内都不能存在其他点。 这样的一个几何特性，通常会被成为内接圆检验法，很容易用数学知识加以解释说明，也正因如此，在进行德劳内三角剖分时也是很常用的方法。假设有这样一个数据集{a，b，c，d}（见图2），其中a，b，c，d是既不共圆也不共线的四个点，将要进行三角剖分的多边形是一个凸四边形，这时有两种方法实现这个数据集的三角剖分，做四边形的一条对角线：方案（1）做对角线ac；方案（2）做对角线bd。这时，我们就可以采用内接圆检验法来确定究竟是方案（1）还是方案（2）才是正确的进行了德劳内三角剖分。 图2：（a）非德劳内三角形；（b）德劳内三角形；（c）四个共圆点连接成边 还有一个可供参考的进行德劳内三角剖分的方法：Sibson（1978）推导出了在一组点中构建三角网的极值法则，Lawson（1972）用这个法则证明在局部等角的三角化方法是唯一的，这个结论同样适用于对