如何实现Delaunay三角剖分外推至预设边界并计算新增顶点高度？

这类场景在测绘、地形建模、CAD几何处理领域是非常常见的常规需求，一线做相关工作的从业者基本都碰到过，对应你提到的两项核心工作，实际经过项目验证的通用做法、避坑要点如下：

三角剖分延伸至预设矩形边界的实操流程 #

• 全程严格保留原凸包内部的三角连接关系不做改动，避免打乱已经有可靠高度值的内部网格结构，只在原三角网凸包和外围矩形边界之间的环形空隙里补网
• 先整理所有参与补网的边界候选点：包括矩形的四个角点、原凸包每条外边界边的延长线与矩形边的所有交点、以及矩形边上按均匀步长补的辅助点（避免补出来的三角形过于狭长影响后续计算精度）
• 从原凸包的每一条外边界边出发，沿着向外的法向方向逐段连接相邻的边界候选点，形成连续的三角带，逐次往矩形边界方向推进，保证补出来的三角形没有重叠、没有空洞，最终整个网格刚好完整覆盖到矩形边界

新增边界顶点的高度值计算方案 #

实际工程里会根据精度、效率要求选不同的外推方法，对应的避坑逻辑很明确：

• 首选低误差的工程方案：自然邻点外推。这是从经典自然邻点插值演化出来的凸包外专用算法，对每个新增的边界点，只取距离最近的、原凸包上的3-5个已知高度顶点，基于这些点组成的局部三角单元的高度梯度做外推，不会出现全局外推带来的边界高度翘边问题。
• 效率优先的简化方案：调高张力系数的最小曲率样条外推。用原凸包内所有已知高度点拟合低次全局样条曲面，直接读取新增边界点位置的曲面高度即可，适合快速出图的场景，注意一定要把样条的张力值调高，不然外推区域很容易出现无意义的高度震荡。
• 必须避开的常见坑：不要直接用反距离加权（IDW）、或者单条边的线性斜率外推。前者在距离原三角网稍远的边界位置会因为权重分配混乱出现高度失真，后者碰到倾斜边界的时候很容易出现相邻边界点的高度不连续，让整个外推网格出现高度断层。

常见同类应用场景 #

这类问题在实际工作里出现频率很高：比如野外采集的地形散点形成的Delaunay三角网，凸包范围刚好比要求的标准图幅矩形边界小一圈，必须把高程延伸到边界才能完成图幅裁切；还有室内点云建模时，扫描得到的三角网范围比房间墙体的预设矩形轮廓小，需要把网格延伸到墙线位置才能和墙面模型无缝拼接。

内容的提问来源于stack exchange，提问作者gamcad