本人最近需要研究三维模型的孔洞修补算法，故上网看了一些大神写的资料，现汇总如下：

重建骨骼比较典型的方法有：用径向基函数从不完整的扫描数据生成连续网格；基于CT等值面数据生成曲面；用傅里叶级数拟合CT图像提取边缘轮廓曲线；基于形状的利用数学形态学算子进行骨架提取的插值算法。

常用的孔洞修补算法可以分为体素方法和几何方法，其中体素方法采用有想距离函数差分或体素滤波，以及在此基础上使用多种偏微分方程进行优化的方法，消除重建三维结构中因采样问题产生的孔洞。基于体素方法存在两方面的问题，首先，它针对的是点云数据，不能直接应用于Mesh模型；其次它是一种物理方法，通常只能够修复重建物体表面的小孔洞，并容易带来几何形变。

对于模型的孔洞修补算法主要采用以下几种策略：

1、直接三角剖分：早期的孔洞修补算法主要考虑如何对孔洞区域进行三角化。三角剖分本质上是一个非线性优化问题，为了获得该问题的最优解，可以使用动态规划算法、遗传算法、模拟退火算法等对孔洞区域进行剖分，以获得最优的三角化结果。已有方法：孔洞多边形的顶点构造新的三角片；各向异性孔洞修补算法。

2、新增采样点：已有方法：分段处理思想，根据孔洞的复杂性将孔洞分为若干个简单的子孔洞；三角片自然增长。

3、采样点调整：已有方法：根据曲率来调整新增点的位置；基于移动最小二乘法的补点方法；基于体素扩散方法来调整新增点的位置；基于径向基函数的调整方法；神经网络的方法。

4、点云数据三角化。

Delaunay三角剖分具有三角剖分最小内角为最大的性质，能够进行任意多连通域有限网格的自动生成可最大限度的保证网中三角形满足近似等边性，避免了过于狭长和尖锐的三角形的出现，是公认的最优三角网。其主流算法有三种：分割-归并法，逐点插入法，三角网生长法。

一个理想的孔洞修补算法具有如下性质：

（1）、自动性。用户如果选择了孔洞，交互应该尽可能简单。

（2）、效率性。合理的运行时间，良好的交互速度。

（3）、准确性。用修补网格来填补孔洞，应该尽可能得和周围网格相融合，网格密度和形状应该和原孔洞边界相匹配。

（4）、鲁棒性。能够处理任意网格的任意孔洞。这是很难被满足的。

《三维颅骨模型的孔洞修补算法研究》中提出一种算法，主要分为五个步骤：

（1）孔洞检测。对输入的颅骨三角网格模型进行孔洞分析，并对检测出的孔洞进行标记。

 孔洞检测室通过自动搜索模型边界点来完成的。

1）模型边界点的提取：判断边界点的过程就是判断该点是否是边界边上的点的过程。

//在一般情况下，空间曲面上点的存储是随机的，所以得到的所有边界点是无序的，对曲面上的点进行边界点判断的时候，可得到与当前边界点的邻接边界点，基于此，即可对所有边界点排序，得到完整的边界曲线。

2）模型边界点的分类和排序

（2）孔洞区域离散点的生成和三角剖分。用户选择需要修补的孔洞，在其中插入离散点，然后对插入的离散点进行三角剖分。

1）孔洞区域离散点插入：对较小的孔洞，直接进行三角剖分；对较大的孔洞，先在孔洞多边形内插入离散点，然后再进行三角化。

2）离散点的三角剖分

（3）网格细分。对三角剖分后的网格进行网格细分。

基本思想：对于每一个三角形，连接各边的中点生成四个新的三角形。

（4）孔洞隐式曲面的建立。利用RBF建立孔洞区域的隐式曲面，调整插入点的位置。

（5）网格平滑。用伞状算子对隐式曲面进行平滑处理。

《面向颅骨损伤修补的缺损重建算法》

以包含缺损结构的三维重建Mesh模型为操作对象，根据交互输入的引导点，通过三角面片的边特征及目标点位置搜索并向外扩展出孔洞上边界环，得到目标问题区域边界，其中内环为孔洞边界，此区域内的所有三角面片为目标问题区域；通过双环采样，投影构建代理面，结合代理面对内环的三角化细化，三角面进行基于代理曲线的 Freeform 变形， 然后在微分域进行保持边缘细节特征的平滑操作， 得到缺失实体孔洞的上缝合面； 通过向颅骨厚度方向扩展， 得到孔洞下边界轮廓， 并搜索出上下边界环之间的三角网格区域； 通过上缝合面的构建方法， 得到缺失实体孔洞的下缝合面； 最后， 连接上下缝合面以及两个边界轮廓环之间的三角网格区域， 获取修复实体。

《A robust hole-filling algorithm for triangular mesh》

主要步骤：

（1）确定三角网格中的孔洞；

（2）对每一个网格模型中的孔洞

1）利用AFM技术生成初始修补网格；

2）基于柏松等式来改善补丁网格：利用谐波等式或者测地线插值来计算法向量；利用局部旋转旋转三角片；解柏松等式得到每一顶点的新坐标；更新坐标得到平滑的补丁网格。

《A piecewise hole-filling algorithm in reverse engineering》

它将一个复杂的孔洞分成几个简单的子洞，子洞连续地进行填补。

结构：

步骤：

（1）对于每一个孔洞确定一个投射曲面或者孔洞向量

（2）投射边界边到曲面上，计算交点

（3）根据交点将孔洞分成几个子孔洞

（4）重复：对每个子孔洞应用平面三角测量算法；将计算的交点参数应用于原始边；细分

通用网格模型的孔洞修补算法BHRA能够较好地处理对不规则区域复杂孔洞的修补，但是不能通用与各种复杂孔洞。

对于大区域复杂孔洞的修补，利用IRS（向内递归求解法）

（1）内部控制辅助点插入

（2）隐式曲面建立

（3）辅助插入点调整到隐式曲面

（4）逐层插入辅助点并重新计算隐式曲面

对于特征区域复杂孔洞的修补，采用TMA（特征模型匹配法）

（1）特征模型匹配

（2）网格融合

孔洞修补研究总结