网格切割原理及流程

最新推荐文章于 2024-04-18 07:51:57 发布

写程序员的bug

最新推荐文章于 2024-04-18 07:51:57 发布

阅读量1.1k

点赞数

分类专栏：网格运算算法文章标签：算法

————来自高贵的CSDN写Bug的程序员

本文链接：https://blog.csdn.net/wy1334526911/article/details/131619499

版权

网格运算同时被 2 个专栏收录

2 篇文章 0 订阅

订阅专栏

算法

2 篇文章 0 订阅

订阅专栏

文章介绍了网格切割的概念，特别是流形网格的定义，包括流行顶点和流行边。它详细阐述了如何通过检测和处理非流形点和边将非流形网格拆分为多个流形块，并使用CGAL库进行相交检测和约束剖分。流程包括合并网格、检测相交、局部remesh、提取流形块、分类Patch以及根据过滤器输出。文章强调了浮点运算精度、空间三角形相交等问题，并提到了在SharpCAD中的应用。

摘要由CSDN通过智能技术生成

流程步骤

假设输入的都是流形网格,自相交理论不会影响正确性，但会影响切割性能。

输入：源网格 A 、切割网格 B 及过滤器。

输出：得到满足过滤器的最终Patch。

步骤 1. 合并 A 和 B

步骤 2. 检测并计算相交对象及记录所在三角形(可并行)

步骤 3. 根据步骤2的结果重新进行局部remesh(可并行)

步骤 4. 将remesh后的所有网格按照非流行点及非流行边拆分成多个patch；

步骤 5. 判断每个patch与输入网格的关系(可并行)；

步骤 6. 根据输入的过滤器输出patch

步骤 1 Merge

合并源网格A 和切割网格B，这一步仅是单纯的合并网格，无其他复杂操作，需要记录源网格A的三角形个数。

步骤 2 Resolve Self-Intersections

这一步其实分为两步

1、检测并记录相交对象

可用CGAL::box_intersection_d及CGAL::box_self_intersection_d检测。这里需要注意一点，如果输入网格存在自相交需要检测所有的三角形，即A中三角形不仅需要与B中所有三角形检测，还需与A中除自身外的其余三角形检测，非常耗时。如果存在自相交但不这样检测可能会影响切割结果。如果不存在自相交则只需要检测网格A的三角形与网格B的三角形。

2、计算并记录相交对象及其三角形。

三角形相交也可用CGAL::intersection来计算，后续文章会介绍CGAL三角形相交的算法。

三空间角形相交的可能：

不在同一平面：点、线段

在同一平面：点、、线段、三角形、四边形、五边形、六边形，以下未展示相交未为点及线段的情况。

需要注意共平面三角形，所有共平面的三角形在后续局部约束剖分时需要一起剖分。

步骤 3 Remesh Intersections

根据步骤2的结果重新进行局部remesh。将三角形边界与相交对象作为约束剖分的约束做三角剖分=>局部remesh,得到重新三角化的局部三角网。可用CGAL 约束剖分处理。相关示例：

#include <CGAL/Exact_predicates_exact_constructions_kernel.h>
#include <CGAL/Constrained_Delaunay_triangulation_2.h>
#include <CGAL/Constrained_triangulation_plus_2.h>

typedef CGAL::Exact_predicates_exact_constructions_kernel K;
typedef CGAL::Constrained_Delaunay_triangulation_2<K> CDT;
typedef CGAL::Constrained_triangulation_face_base_2<K> Face;
typedef CGAL::Triangulation_vertex_base_2<K> Vertex;
typedef CGAL::Triangulation_data_structure_2<Vertex, Face> TDS;
typedef CGAL::Constrained_triangulation_plus_2<CDT, TDS> CTP;

int main() {
  CTP ctp;
  // 添加约束边
  ctp.insert_constraint(K::Point_2(0,0), K::Point_2(1,1));
  ctp.insert_constraint(K::Point_2(1,1), K::Point_2(2,0));
  ctp.insert_constraint(K::Point_2(2,0), K::Point_2(0,0));
  // 插入点
  ctp.insert(K::Point_2(1,0.5));
  // 输出所有三角形
  for (CTP::Finite_faces_iterator fit = ctp.finite_faces_begin(); fit != ctp.finite_faces_end(); ++fit) {
    std::cout << "Triangle: " << ctp.triangle(fit) << std::endl;
  }
  return 0;
}

后续文章会介绍约束剖分算法。

步骤 4 Extract manifold patch

将步骤3中得到的局部三角网替换对应原三角形，按照非流行边拆分为多个Patch。提取流行块时需要记录Patch的临边情况以及非流形边周围的三角形情况，此外哪些流行块包含哪些三角形也是需要记录的。

步骤 5 Classify Patch

判断每个Patch与输入网格的关系。这一步是所有布尔运算的难点。我们采用相交三角形法向来判断Patch与原始网格的关系(可能会有一些特殊情况需要处理)，即Patch的产生是因为三角形相交或非流形边导致的。下图两个三角形相交产生的公共边即为非流形边，这条非流形边将Tri2分为上下两部分，Tri1完全切割Tri2，但Tri1不能被Tri2完全分割，所以会有Undefine这个过滤器。判断关系时我们不需要对每个Patch判断，只需要判断非流形边邻接，当然会有一些Patch不是由非流形边提取的，比如输入存在多个连通区域，其中一些连通区域的所有三角形均不与其他区域三角形相交，像这种情况，如果另一个输入网格是闭合的则可以判断当前Patch与其关系，如果不是闭合的则归为Undefine.