三维网格骨架提取

在面网格、体网格骨架提取的工作中，网格本身的拓扑连接关系是我们的关注点，而非网格各体素的坐标，体素就是体网格中每一个小方块，与黑白像素一样，0代表体素为空，1代表有体素，因此六面体网格是二进制体素模型的一种编码方式。在体素中也可以分别用实（solid）、虚（empty）来表示。骨架提取的工作在于，对实集合中大部分的体素逐个转换为虚集合内的体素，直到实集只剩下一个逐体素展开的链条，在数字拓扑关系的规则下，每一次对虚实集合的交界处进行处理。

1.数字拓扑

欧拉示性数是网格的内蕴度量，保持处理前后欧拉示性数的不变，从而间接保持网格的拓扑关系。欧拉示性数有两种变现方式，对一个体网格来说，欧拉示性数等于体网格内点的数量-边的数量+面的数量-体素的数量，对体网格求次导是这个网格的表面，等于网格表面点的数量-边的数量+面的数量。表面的欧拉示性数等于体网格欧拉示性数的两倍。体网格的欧拉示性数也和网格整体的连通分量数、把手数、窟窿数相关。比如球体的是1（1-0+0），环的是0（1-1+0）
关于体素有两种邻接表示，基于体素6个面相邻的6邻居，基于体素8个顶点相邻的26邻居（立方体333的中心）。实集中一般为26邻，虚集一般为6邻。如果两个集合用了一样的邻居定义会导致结果的模糊性，比如平面情况下的Jordan曲线理论：一条闭合的曲线能把平面分成内部区域和外部区域。
欧拉示性数的计算一般是通过遍历体素的顶点，一个顶点最多会被8个体素共有，这8个体素称为顶点的八分圆（Octant），八分圆之间相互重叠，边、面、体素都有可能在多个八分圆中重复出现。在一个立方体（2*2）中，八分圆的组合有22种（考虑对称性，C（0，4）+C（1，4）+C（2，4）=（1+4+6）*2），利用每个点的八分圆计算欧拉示性数，如下式：

2.骨架提取

骨架提取是一个连续的实体素删除过程。边界实体素的6邻居中必然有一个虚体素，每一个边界实体素也被称为简单点，删除它不会改变网格的拓扑关系，也就是网格的连通分量数、把手数（handle，h）、窟窿数（cavity，c）不会变。一个可以删除的简单点，在删除前后，必然满足下列式子：

我们也可以标记一些体素为不可删除的，在拓扑优化的结果中，我们可以把这些体素定义为满足第一边界条件和第二边界条件。当没有体素满足上面的式子或者是被标记为不可删除，剩余的实体素就是提取的骨架。

参考自《Topologically robust CAD model generation for structural optimisation》