VTK多边形拉伸开发：孔洞添加方法与曲线异常问题排查

问题根因 #

• 孔洞实现逻辑错误：VTK原生VTK_POLYGON单元不支持直接附加孔洞定义，官方规范中带孔洞的多连通面需要先按绕序规则组织轮廓，经三角化生成拓扑合法的面网格后才能做后续操作。现有代码仅创建了单外轮廓多边形单元，未处理孔洞环逻辑，也未做三角化，拉伸过滤器无法识别非法拓扑结构，会生成自交扭曲面。
• 多边形异常原因：代码未校验轮廓绕序、未做环的拓扑关联，直接按传入点序连接单元，点序不符合VTK多边形的绕序要求时就会出现面连接错误。

修正方案 #

带洞多边形的VTK标准组织规则 #

1. 所有点（外轮廓点+所有孔洞轮廓点）统一存入vtkPoints
2. 外轮廓点按逆时针方向排序，每个孔洞的内轮廓点按顺时针方向排序（与外轮廓绕序相反）
3. 不直接创建VTK_POLYGON单元，而是将每个闭合环（外轮廓+每个孔洞）分别创建为闭合的VTK_POLY_LINE单元存入vtkCellArray

合法面生成 #

使用vtkContourTriangulator处理组织好的轮廓线，自动生成带孔洞的三角化曲面，这是VTK官方推荐的带洞面构建方式，不存在直接给多边形单元添加孔洞的公开接口。

修正后可运行代码 #

import vtk

def extrude_polygon(outer_vertices, holes, extrude_height):
    """
    带孔洞多边形拉伸
    :param outer_vertices: 外轮廓顶点列表，格式[(x,y), ...]，需为逆时针绕序
    :param holes: 孔洞列表，每个元素为单个孔洞的顶点列表[(x,y), ...]，需为顺时针绕序
    :param extrude_height: Z轴方向拉伸高度
    """
    points = vtk.vtkPoints()
    points.SetDataTypeToDouble()
    lines = vtk.vtkCellArray()

    # 写入外轮廓点与闭合多段线
    for x, y in outer_vertices:
        points.InsertNextPoint(x, y, 0)
    outer_line = vtk.vtkIdList()
    for i in range(len(outer_vertices)):
        outer_line.InsertNextId(i)
    outer_line.InsertNextId(0)  # 闭合外轮廓
    lines.InsertNextCell(outer_line)

    # 写入所有孔洞的点与闭合多段线
    current_point_id = len(outer_vertices)
    for hole in holes:
        start_id = current_point_id
        for x, y in hole:
            points.InsertNextPoint(x, y, 0)
        hole_line = vtk.vtkIdList()
        for i in range(len(hole)):
            hole_line.InsertNextId(start_id + i)
        hole_line.InsertNextId(start_id)  # 闭合孔洞轮廓
        lines.InsertNextCell(hole_line)
        current_point_id += len(hole)

    # 构建轮廓数据集
    contour_pd = vtk.vtkPolyData()
    contour_pd.SetPoints(points)
    contour_pd.SetLines(lines)

    # 三角化生成拓扑合法的带洞面
    triangulator = vtk.vtkContourTriangulator()
    triangulator.SetInputData(contour_pd)
    triangulator.Update()

    # 线性拉伸
    extr = vtk.vtkLinearExtrusionFilter()
    extr.SetInputConnection(triangulator.GetOutputPort())
    extr.SetVector(0, 0, extrude_height)
    extr.SetExtrusionTypeToNormalExtrusion()
    extr.CappingOn()  # 自动封闭拉伸体上下底面

    # 计算法向
    pn = vtk.vtkPolyDataNormals()
    pn.SetInputConnection(extr.GetOutputPort())
    pn.AutoOrientNormalsOn()
    pn.ConsistencyOn()
    pn.SplittingOff()
    pn.Update()

    return pn

使用注意事项 #

• 必须保证轮廓绕序正确：外轮廓逆时针、内轮廓顺时针。如果原始点序不符合要求，可以通过计算多边形有向面积判断绕序，不符合时反转点列表即可，也可使用vtkContourLoopExtraction自动校正绕序。
• 禁止跳过三角化步骤直接拉伸带洞轮廓：vtkLinearExtrusionFilter仅支持拉伸无孔洞的单连通简单多边形，多连通带洞面必须先三角化保证拓扑合法。
• 拉伸时开启CappingOn()选项，否则生成的拉伸体上下底面为镂空状态，无法得到封闭实体。

效果参考 #

• 错误输出：
  
• 预期输出：
  

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