低潮线是指海水在海岸上退落的最低线界，在海洋划界、海岸工程建设和滩涂资源管理等应用中具有重要意义。根据《联合国海洋法公约》，沿岸低潮线是领海基线及基点确定的重要一环，其位置的准确与否对海洋划界具有直接影响。此外，低潮线的位置及长度也是海岸带调查中最基础的数据之一，其对滩涂面积计算分析等具有重要意义。

低潮线是指海水在海岸上退落的最低线界，在海洋划界、海岸工程建设和滩涂资源管理等应用中具有重要意义。根据《联合国海洋法公约》，沿岸低潮线是领海基线及基点确定的重要一环，其位置的准确与否对海洋划界具有直接影响。此外，低潮线的位置及长度也是海岸带调查中最基础的数据之一，其对滩涂面积计算分析等具有重要意义。

如何提取较高精度的沿岸低潮线，一直是海洋划界、平缓滩涂监测等领域面临的难题。低潮线提取主要有3种方法：

①实地测量方法；

②基于影像水边线的推算方法；

③基于潮滩DEM 和低潮面模型的简单求交方法。

实地测量方法是利用测量仪器在低潮线附近每隔一定距离采集特征点，标记在数字地形图上并连接成低潮线。该方法必须进行大量野外工作，效率低，而且，由于低潮线并无明显标记，实际中也很难测到其准确位置；基于影像水边线的推算方法，首先基于影像数据提取出瞬时水位线，然后对水位线进行潮位校正得到低潮线，但当前的推算做法基本上都假定潮滩地形是规则的斜坡，这明显不适用于地形起伏较大的区域；简单求交方法根据潮汐模型计算当地低潮面（深度基准面）的高程，以此高程面与潮滩DEM 横切，采用等值线追踪等技术得到低潮线。该方法所需人工干预少，在小范围区域内提取低潮线效果相对较好，目前得到了较为广泛的应用。然而，该方法在提取低潮线时将低潮面简化为平面处理，而事实上，对于范围较大的沿岸区域，特别是潮汐变化差异较大时，沿岸低潮面是一个曲面而不是简单的平面，由此会引起低潮线提取的较大偏差。尽管实践中也可以采用分区提取低潮线再进行拼接的方法，但这种方法很难得到最合理的分区，若分区粗略，会使区域内的低潮面呈不连续的阶梯状分布，不足以充分反映区域内的潮汐变化差异，可能导致低潮线提取精度不高及不连续现象；若分区较细，虽然能保证低潮线精度，但分区及低潮线拼接等处理复杂。

针对以上问题，提出一种基于曲面求交的低潮线提取方法。首先采用单向分区策略构建潮滩数字高程模型，然后将低潮面视为曲面进行处理，构建沿岸网格低潮面模型，表达曲面形态的沿岸低潮面，最后，对所构的两个网格模型进行曲面求交，提取沿岸低潮线。

潮滩数字高程模型的构建需要利用多源潮滩数据，而沿岸网格低潮面模型的构建需要利用已有的沿岸低潮面数据，为完成曲面求交处理，需要将不同源数据统一到同一基准，本文数据平面基准采用CGCS2000坐标系统，垂直基准采用1985国家高程基准。

潮滩数据来源复杂多样，需要进行多源潮滩数据的集成与融合，将数据统一到CGCS2000坐标系统和1985国家高程基准，潮滩DEM 结构选择三角网，建模的基本方法采用逐点插入法。

从大范围潮滩数据呈现狭长带状分布的特点出发，为优化潮滩模型的构建和方便曲面求交及拼接，提高整体建模效率，在构建潮滩DEM 时采用基于分割线的单向分区策略。

（1）建立分割线缓冲区。以分割线为中心线构建缓冲区，缓冲距h可以设为源数据间距的若干倍；将位于缓冲区的点存入分割线类中，并确定分割线缓冲区的x 值范围。通过构建分割线缓冲区，可以提高共公边界点的确定速度。

（2）构建相邻分区的公共边界，即计算分割线上新加点的坐标。根据源数据密度确定合适的间距，确定分割线上新加点的数目及坐标。新加点的高程值可以通过内插法计算获得。内插计算时，在分割线缓冲区内搜索离散点，采用反距离加权等插值方法进行计算，这样可以保证新加点的精度、密度与原始数据基本相符。

（3）分区构建TIN－DEM 及各分区TINDEM的合并。基于分割线的单向分区策略，充分利用潮滩数据狭长带状分布的特点，通过构建相邻分区的公共边界使相邻分区TIN－DEM 可以直接拼接成整体TIN－DEM，避免了传统方法在分块拼接时需要搜索子网边界公共切线及三角网拓扑关系重构。在进行曲面求交提取沿岸低潮线时，通过该方法构建分区DEM，分别进行曲面求交，然后由分区低潮线拼接组成整体低潮线。由于相邻分区之间具有公共边，因此，分区低潮线可以直接相连组成整体低潮线。

针对常用低潮面模型在较大范围区域内呈阶梯状断续分布的不足，为实现曲面求交提取低潮线，需要构建连续、曲面形态的沿岸低潮面模型。

虽然结合卫星测高数据和潮汐动力学方程，可以构建曲面形态的低潮面模型，但沿岸海域潮波传播复杂，该方法构建的低潮面模型在沿岸海域精度可能不高，且可能存在不少的无效网格点。若单纯地利用外推法估计无效网格点的低潮面值，在某些区域可能会产生较大误差。为此，提出的TCARI方法，顾及边界条件约束，综合利用沿岸验潮站和网格潮汐模型数据信息（低潮面源数据利用调和方法或潮位传递模型计算获得），采用拉普拉斯插值构建精度更高的曲面形态的沿岸网格低潮面模型。

构建的潮滩数字高程模型和沿岸低潮面模型分别采用三角网和规则格网，由于规则格网可以转化成三角网，两模型的曲面求交的核心问题是两三角网格曲面的求交问题。

关于三角网格曲面求交，已有学者进行了研究。在已有算法的基础上，结合潮滩数字高程模型和沿岸低潮面模型的特点，采用以下策略对两三角网格曲面求交进行优化。具体步骤如下：

（1）建立三角网格曲面的层次包围盒，优化曲面的碰撞检测。层次包围盒的构造过程为：

① 构造初始包围盒B0=B1∩B2，B1、B2分别是两个三角网格曲面的轴向包围盒。

② 确定两个曲面中与初始包围盒B0相交的三角形集合SB1、SB2。

③ 根据三角网格曲面的三角形数量，确定包围盒剖分层次（顾及算法效率，剖分层次根据经验选择5层），对初始包围盒B0进行剖分，得到子包围盒Bi。

（2）利用拓扑关系追踪交线段，完成曲面交线的快速重构。

三角网格曲面间的交线是由三角形对的交线段所组成的交线链或交线环。由于沿着曲面上交线的行进方向，相邻交线段所在的三角形也是相邻的。根据这一特性，利用拓扑关系追踪三角形的邻接三角形，得到三角网格曲面的交线。依次对两三角网格曲面的交线段进行追踪，得到最终的曲面交线。

选择南海某一较大范围沿岸区域的相关数据进行试验验证与分析。潮滩高程数据由部分实测数据和相关海图数据集成融合后经过适当内插得到。

低潮面数据采用区域内的验潮站低潮面数据和基于卫星测高数据构建的网格低潮面数据，边界约束数据采用试验区域内的岸线数据，构建的沿岸低潮面模型网格分辨率为1′。将潮滩数据和低潮面数据统一到同一基准（CGCS2000平面基准，1985国家高程基准）。

低潮线并无明显标记，实际中难以测到其准确位置，只能采用对比试验和理论分析的方式对其精度进行评估。为方便叙述，本文提出的基于曲面求交的低潮线提取方法简称曲面求交方法。在试验区域内分别采用曲面求交方法和简单求交方法提取低潮线，并对提取结果进行对比及分析，以验证曲面求交方法提取低潮线的有效性。具体的试验设计思路如下：

（1）利用潮滩DEM 构建方法构建潮滩DEM；利用TCARI方法构建沿岸低潮面模型；对潮滩DEM 和沿岸低潮面模型进行曲面求交提取出低潮线。

（2）建立一类对比区域，即小范围对比区域。该类区域是以各验潮站为中心，范围为±2km 的矩形区域，在区域内可以将低潮面认为是一平面。在对比区域内构建局部潮滩DEM，采用简单求交方法提取局部低潮线。

（3）在一类对比区域内，对潮滩DEM 作一系列剖面线（间隔为50m），求得剖面线与曲面求交方法提取低潮线的交点和剖面线与简单求交方法提取低潮线的交点，计算两交点（即检查点对）之间的距离偏差；计算多个检查点对，得到统计结果。

（4）建立二类对比区域，即较大范围对比区域；该类区域是以各验潮站为中心，范围为±5km的矩形区域，模拟分区提取低潮线方法的分区区域。在对比区域内构建局部潮滩DEM，采用简单求交方法提取局部低潮线。

（5）在二类对比区域内，对潮滩DEM 作一系列剖面线（间隔为50m），求得剖面线与曲面求交方法提取低潮线的交点和剖面线与简单求交方法提取低潮线的交点，计算两交点（即检查点对）之间的距离偏差；计算多个检查点对，得到统计结果。

分别采用曲面求交方法和简单求交方法提取6块对比区域的低潮线并沿X 方向作剖面线，计算两种低潮线与每个剖面线的离散交点对（检查点对）的距离偏差Δx，得出统计结果。通过结果对比和理论分析，可以进一步说明两种方法提取低潮线的精度差异。表中，平均低潮面差值指检查点对低潮面差值的平均值。

在一类对比区域内，两种方法提取的低潮线偏差平均值和中误差都较小。这表明，在小范围区域内，曲面求交方法与简单求交方法提取的低潮线精度相当。其原因是，在小范围区域内，潮汐变化较小，低潮面可以简化为一个平面，低潮面与潮滩DEM 曲面求交转换为平面与DEM 求交，即前文提到的简单求交。但如果区域内潮汐变化差异较大，则需要更精细的潮汐模型。在二类对比区域内，两种方法提取的低潮线偏差的平均值和中误差都较大。这表明，在该类区域内，两种方法提取的低潮线偏差较大。其原因是，该类区域相当于分区提取方法的子分区，在分区范围内，潮汐变化差异仍可能较大（从表中也可以看出，此时低潮面差值较大），简单求交方法将低潮面简化为平面处理，可能不足以充分反映潮汐变化的这种差异，而曲面求交方法通过构建曲面形态的沿岸低潮面模型，并直接对潮滩DEM 和低潮面模型进行曲面求交，可以更有效地提取低潮线。

可以看出，两种方法提取的低潮线在验潮站附近非常吻合，而在距离验潮站较远的区域偏差较大（如右侧的局部放大图）。这表明，在提取沿岸低潮线时，尤其是大范围区域内，需要构建更为精细的低潮面模型，以充分顾及区域内潮汐变化的差异。

由上述试验和分析可得，当提取低潮线的区域较大时，简单求交方法与曲面求交方法提取的低潮线偏差较大。由于曲面求交方法通过构建曲面形态的沿岸低潮面模型，直接采用曲面求交的方法提取低潮线，相比简单求交方法，更能充分顾及区域内的潮汐变化差异，因此，从理论上分析，曲面求交方法提取的低潮线精度更高。

低潮线位置偏差较大的检验点多分布于地形坡度较小的区域，而坡度较大区域的低潮线位置偏差较小。这表明，低潮线提取精度与地形坡度有密切关系，地形坡度较小的区域，提取的低潮线精度较差。其主要原因是，低潮面模型在垂直方向的误差，通过一定方式转化成水平位置误差，坡度越小，偏差的放大效应越明显。因此，对于地形平坦区域的低潮线提取，需要格外重视沿岸低潮面模型的精度。

针对当前常用的低潮线提取方法不能适应较大区域范围的现状，提出一种基于曲面求交的低潮线提取方法。通过分析、计算及试验比对，得出如下结论：

（1）大范围的沿岸区域内，潮汐变化差异可能较大，相比基于简单求交的低潮线提取方法，本文所提方法通过构建曲面形态的沿岸低潮面模型，直接采用曲面求交的方法提取低潮线，能够更充分地顾及大区域内潮汐变化的差异，可以提取精度更高的沿岸低潮线。

（2）低潮线提取精度与地形坡度有密切关系，相比坡度较大的区域，坡度较小的区域提取的低潮线精度较差。对于地形平坦区域的低潮线提取，需要格外重视沿岸低潮面模型的精度。