VTK实现面模型切割

有任意vtkPolyData数据，要通过vtkCutter切割，形成切面。使用简单，但细节处理复杂：

vtkPolyData中的几何数据，可能为封闭面，也可能为不封闭面或线，要实现自适应切割，需要进行特殊处理。

通常的步骤如下：

1.vtkCutter切割目标vtkPolyData

           vtk.vtkCutter cutEdges = new vtk.vtkCutter();
            cutEdges.SetInput(srcObject);
            cutEdges.SetCutFunction(clipPlane);

2.得到vtkPolyData对象，是一段一段的线。

    vtkPolyData cutPanel = cutEdges.GetOutput();

3.使用vtkStripper将vtkPolyData进行条带化，形成多义线Polyline.

     vtk.vtkStripper cutStrips = new vtk.vtkStripper();
     cutStrips.SetInput(cutPanel);

4.对Polyline进行封装，形成Ploygon。

    vtk.vtkPolyData cutPoly = new vtk.vtkPolyData();
     cutPoly.DeepCopy(cutStrips.GetOutput());

     cutPoly.SetPolys(cutStrips.GetLines());

5.对Polygon进行三角化，实现三角面的显示。

   vtkTriangleFilter cutTriangles = new vtk.vtkTriangleFilter();
     //cutTriangles.PassLinesOff();
    //cutTriangles.PassVertsOn();
    cutTriangles.SetInput(cutPoly);
    return cutTriangles.GetOutput();

 

对第4步，如果为非封闭体，就需要判断cutStrips.GetLines()中单元的属性，如果单元的首节点与尾节点的序号一致，才认为是封闭的多变形，否则就是非封闭的多义线，不能直接加入到Polys中，仍然要加载在Lines单元中。

代码如下：

       vtkCellArray polyAry = new vtkCellArray();
                    vtkCellArray lineAry = new vtkCellArray();
                    for (int i = 0; i < iCells; i++)
                    {
                        vtkCell cell = cutStrips.GetOutput().GetCell(i);
                        if (cell.GetNumberOfPoints() > 0)
                        {
                            int iId = cell.GetPointId(0);
                            int iId1 = cell.GetPointId(cell.GetNumberOfPoints() - 1);
                            if (iId != iId1)
                            {
                                lineAry.InsertNextCell(cell);
                            }
                            else
                            {
                                polyAry.InsertNextCell(cell);
                            }
                        }
                    }

                    cutPoly.SetLines(lineAry);
                    cutPoly.SetPolys(polyAry);

 

为什么要第5步： 第5步的三角化必不可少，因为第4步得到的多变形可能是凹多边形，渲染时将会出现错误，必须处理。