Dear all,
I try to extract the boundary nodes of an unstructured grids. Based on the example
https://vtk.org/Wiki/VTK/Examples/Cxx/Meshes/BoundaryEdges
I used the following code:
vtkSmartPointer grid = vtkSmartPointer::New();
vtkSmartPointer filter = vtkSmartPointer::New();
vtkSmartPointer edges = vtkSmartPointer::New();
/* load mesh into grid */
…
/* extract the boundary points */
filter->SetInputData(grid);
filter->Update();
edges->SetInputData(filter->GetOutput());
edges->BoundaryEdgesOn();
edges->FeatureEdgesOff();
edges->NonManifoldEdgesOff();
edges->ManifoldEdgesOff();
edges->Update();
vtkSmartPointer polyData = edges->GetOutput();
When I print out the number of the points in polydata (polydata->GetNumberOfPoints()) I obtain 0, i.e.
polydata is empty. When I comment out the lines
edges->NonManifoldEdgesOff();
edges->ManifoldEdgesOff()
then I get some subset of points in polydata extrated from the grid but when I visualize these points I see
that some boundary points are missing and there are some inner nodes in polydata. Thus
polydata does not contain the exact set of the boundary points of the grid.
Now I am not sure if this is the correct way to extract the boundary points from an unstructured grid
and I would like to ask if there is any. Thank you very much for your help!
Kind regards,
Zoltan