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步得到的多变形可能是凹多边形，渲染时将会出现错误，必须处理。