Pymeshfix/Pyvista网格修复异常：牙桥模型被拆分问题求助

问题描述 #

我尝试修复一个存在孔洞与相交三角面问题的牙桥模型网格，使用tin.fill_small_boundaries()填充孔洞，检测到相交三角面后调用tin.clean(max_iters=10, inner_loops=3)处理，但修复后牙桥模型被拆分，最终仅保留单个牙体。

原始代码 #

# -*- coding: utf-8 -*-
import pymeshfix


def plotMesh(tin):
    
    # 返回顶点和面数组
    vertices, faces = tin.return_arrays()
    
    # 从顶点和面数组创建MeshFix对象
    meshfix = pymeshfix.MeshFix(vertices, faces)

    # 绘制输入网格
    meshfix.plot()


# 创建TMesh对象
tin = pymeshfix.PyTMesh()

# 加载stl文件
tin.load_file("bridge08.stl")

# 可视化输入网格
plotMesh(tin)


###############################################################################
# 清理流程

# 尝试合并邻近组件
# tin.join_closest_components()

# 填充孔洞
tin.fill_small_boundaries()
print('当前边界数量: {:d}'.format(tin.boundaries()))

# 检测相交三角面
faces = [tin.select_intersecting_triangles()]

# 如果检测到相交三角面
if faces:
    # 清理自相交
    tin.clean(max_iters=10, inner_loops=3)
###############################################################################


# 可视化输出网格
plotMesh(tin)

解决方案 #

问题根源在于操作顺序不合理和clean参数过于激进，导致有效网格被误删，可通过以下步骤修正：

1. 调整操作顺序：先清理相交面，再填充孔洞
   原流程先填充孔洞，新增的填充面会与原有相交面形成更复杂的冲突，触发clean时大量删除有效网格。改为先处理相交三角面，再填充孔洞，能减少冲突概率。

2. 降低clean操作的激进程度
   inner_loops=3和max_iters=10会让清理迭代次数过多，导致过度删除。建议降低参数，并启用小相交区域的温和处理：

   tin.clean(max_iters=5, inner_loops=1, resolve_small_intersections=True)
   

3. 启用组件合并
   取消注释tin.join_closest_components()，在清理后合并可能被拆分的组件，避免模型碎片化。

4. 分步验证每一步结果
   在每一步操作后添加可视化，定位哪一步导致模型异常，方便调试。

修改后的完整代码 #

# -*- coding: utf-8 -*-
import pymeshfix


def plotMesh(tin):
    vertices, faces = tin.return_arrays()
    meshfix = pymeshfix.MeshFix(vertices, faces)
    meshfix.plot()


# 加载模型
tin = pymeshfix.PyTMesh()
tin.load_file("bridge08.stl")
plotMesh(tin)

# 调整后的清理流程
# 1. 先检测并清理相交三角面
intersecting_faces = tin.select_intersecting_triangles()
if intersecting_faces:
    tin.clean(max_iters=5, inner_loops=1, resolve_small_intersections=True)
    plotMesh(tin)  # 验证清理后的模型状态

# 2. 填充孔洞
tin.fill_small_boundaries()
print('剩余边界数量: {:d}'.format(tin.boundaries()))
plotMesh(tin)  # 验证填充后的模型状态

# 3. 合并拆分的组件
tin.join_closest_components()
plotMesh(tin)  # 查看最终修复结果

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