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