在自定义Paraview过滤器中使用vtkPolyLine时出现分段错误

时间:2016-11-20 13:25:15

标签: c++ 3d segmentation-fault vtk paraview

我想使用vtkPolyLine显示多组3D点。 这些点在多维向量中存储为节点(自定义类):     vector<vector <Node> > criticalLines;节点有double posX; double posY; double posZ;来存储其位置。

对于以下部分,我尝试使用与此示例类似的vtkPolyLine: http://www.paraview.org/Wiki/VTK/Examples/Cxx/GeometricObjects/PolyLine

在向量填充节点后调用此函数:

void Algorithm::displayLines(vtkSmartPointer<vtkPoints> points,vtkSmartPointer<vtkCellArray> lines)
{
for(int i = 0; i<criticalLines.size(); i++)
{
    if(criticalLines[i].empty())
    {
        continue;
    }

    vtkSmartPointer<vtkPolyLine> polyLine =
        vtkSmartPointer<vtkPolyLine>::New()

    for(int j =0; j< criticalLines[i].size(); ++j)
    {

        vtkIdType idx=points->InsertNextPoint(criticalLines[i][j].posX,
                      criticalLines[i][j].posY,
                      criticalLines[i][j].posZ);
        //print posX,posY,posZ of current Node
        criticalLines[i][j].PrintSelf();


        //Seg. Fault occurs here 
        polyLine->GetPointIds()->SetId(j,idx);


    }
    lines->InsertNextCell(polyLine);

}

}

pointslines都在Algorithm.h文件中定义,并在构造函数中初始化,如下所示:

points = vtkSmartPointer<vtkPoints>::New();
lines = vtkSmartPointer<vtkCellArray>::New();

稍后又添加到vtkPolyData

vtkSmartPointer<vtkPolyData> opd=vtkSmartPointer<vtkPolyData>::New() ;
opd->SetPoints(algorithm.points);
opd->SetLines(algorithm.lines);

criticalLines[i][j].PrintSelf();的输出按预期显示值。

使用vtkSmartPointer<vtkTriangle> triangle = vtkSmartPointer<vtkTriangle>::New();代替vtkPolyLine时,一切正常。

这个create multiple polylines given a set of points using vtk某种相关问题的解决方案似乎并不是我想要的。

我不确定我的代码中缺少/错误。 如果您需要更多信息,请与我们联系。

非常感谢任何帮助!

1 个答案:

答案 0 :(得分:1)

您的vtkPolyLine需要为点ID分配一些空间,例如

polyLine->GetPointIds()->SetNumberOfIds(5);
您链接到的示例中的

。在您的情况下,您需要致电

polyLine->GetPointIds()->SetNumberOfIds(criticalLines[i].size());
创建polyLine后立即