VTK 三角剖分 Delaunay2D

Description

随机生成一系列网格点，进行三角剖分。同时获取边界点集，进行边界的缝合。

Code

Delaunay2D.cxx

#include <vtkSmartPointer.h>

#include <vtkPoints.h>
#include <vtkPolyData.h>
#include <vtkTriangle.h>
#include <vtkCellArray.h>

#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkProperty.h>
#include <vtkNamedColors.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>

#include <vtkDelaunay2D.h>
#include <vtkFeatureEdges.h>
#include <vtkStripper.h>
#include <vtkCleanPolyData.h>
#include <vtkAppendPolyData.h>
#include <vtkVertexGlyphFilter.h>


int main(int, char *[])
{
  //生成第一个三角网
  vtkSmartPointer<vtkPoints> ps1 =vtkSmartPointer<vtkPoints>::New();
 
  unsigned int GridSize = 10;

  for(unsigned int x = 0; x < GridSize; x++)
  {
    for(unsigned int y = 0; y < GridSize; y++)
    {
		ps1->InsertNextPoint(x+ vtkMath::Random(.1, .2), y+ vtkMath::Random(.1, .2), vtkMath::Random(-.25, .25));
	}
  }

  vtkSmartPointer<vtkPolyData> polydata1 =vtkSmartPointer<vtkPolyData>::New();
  polydata1->SetPoints(ps1);

  vtkSmartPointer<vtkDelaunay2D> delaunay1 =vtkSmartPointer<vtkDelaunay2D>::New();
  delaunay1->SetInputData(polydata1);
  delaunay1->Update();

  //获取三角网边界点集
  vtkSmartPointer<vtkFeatureEdges> featureEdges1 =vtkSmartPointer<vtkFeatureEdges>::New();
  featureEdges1->SetInputConnection(delaunay1->GetOutputPort());
  featureEdges1->BoundaryEdgesOn();
  featureEdges1->FeatureEdgesOff();
  featureEdges1->ManifoldEdgesOff();
  featureEdges1->NonManifoldEdgesOff();
  featureEdges1->Update();
  
  //使点集有序
  vtkSmartPointer<vtkStripper> stripper1 =vtkSmartPointer<vtkStripper>::New();
  stripper1->SetInputConnection(featureEdges1->GetOutputPort());

  vtkSmartPointer<vtkCleanPolyData> cleanPolyData1 =
	  vtkSmartPointer<vtkCleanPolyData>::New();
  cleanPolyData1->SetInputConnection(stripper1->GetOutputPort());
  cleanPolyData1->Update();
  
  //拷贝到边界点集中
  vtkSmartPointer<vtkPoints> boundary = vtkSmartPointer<vtkPoints>::New();
  vtkPoints *pt = cleanPolyData1->GetOutput()->GetPoints();
  int num = pt->GetNumberOfPoints();
  for (int i = 0; i < num; ++i) {
	  double *p = pt->GetPoint(i);
	  boundary->InsertNextPoint(p);
  }
 //生成第二个三角网
  vtkSmartPointer<vtkPoints> pts2 = vtkSmartPointer<vtkPoints>::New();
  for (unsigned int x = 0; x < GridSize; x++)
  {
	  for (unsigned int y = 0; y < GridSize; y++)
	  {
		  pts2->InsertNextPoint(x+ vtkMath::Random(.1, .2), y+ vtkMath::Random(.1, .2), vtkMath::Random(.5, .75));
	  }
  }
  vtkSmartPointer<vtkPolyData> polydata2 = vtkSmartPointer<vtkPolyData>::New();
  polydata2->SetPoints(pts2);
  vtkSmartPointer<vtkDelaunay2D> delaunay2 = vtkSmartPointer<vtkDelaunay2D>::New();
  delaunay2->SetInputData(polydata2);
  delaunay2->Update();


  vtkSmartPointer<vtkFeatureEdges> featureEdges2 =vtkSmartPointer<vtkFeatureEdges>::New();
  featureEdges2->SetInputConnection(delaunay2->GetOutputPort());
  featureEdges2->BoundaryEdgesOn();
  featureEdges2->FeatureEdgesOff();
  featureEdges2->ManifoldEdgesOff();
  featureEdges2->NonManifoldEdgesOff();
  featureEdges2->Update();

  vtkSmartPointer<vtkStripper> stripper2 =vtkSmartPointer<vtkStripper>::New();
  stripper2->SetInputConnection(featureEdges2->GetOutputPort());

  vtkSmartPointer<vtkCleanPolyData> cleanPolyData2 =vtkSmartPointer<vtkCleanPolyData>::New();
  cleanPolyData2->SetInputConnection(stripper2->GetOutputPort());
  cleanPolyData2->Update();
  //将第二个边界点集拷贝到边界点集中
  pt = cleanPolyData2->GetOutput()->GetPoints();
  num = pt->GetNumberOfPoints();
  for (int i = 0; i < num; ++i) {
	  double *p = pt->GetPoint(i);
	  boundary->InsertNextPoint(p);
  }

  //边界点集生成新的拓扑单元
  vtkCellArray* cells = vtkCellArray::New();
  vtkTriangle* vtkTriangle;

  for (int i = 0; i < num; ++i) {
	  
	  vtkTriangle = vtkTriangle::New();
	  vtkTriangle->GetPointIds()->SetId(0,i);
	  vtkTriangle->GetPointIds()->SetId(1, ((i+1)%num));
	  vtkTriangle->GetPointIds()->SetId(2, (num+i));

	  cells->InsertNextCell(vtkTriangle);

	  vtkTriangle = vtkTriangle::New();
	  vtkTriangle->GetPointIds()->SetId(0, ((i + 1) % num));
	  vtkTriangle->GetPointIds()->SetId(1, ((i + 1) % num + num));
	  vtkTriangle->GetPointIds()->SetId(2, (num + i));

	  cells->InsertNextCell(vtkTriangle);
  }
  vtkPolyData* trianglePolyData = vtkPolyData::New();

  trianglePolyData->SetPoints(boundary);
  trianglePolyData->SetPolys(cells);

  // combine two poly data
  // 上下三角网和侧面三角网组合
  vtkSmartPointer<vtkAppendPolyData> appendFilter =vtkSmartPointer<vtkAppendPolyData>::New();
  appendFilter->AddInputData(delaunay1->GetOutput());
  appendFilter->AddInputData(delaunay2->GetOutput());
  appendFilter->AddInputData(trianglePolyData);
  appendFilter->Update();

  // Visualize
  vtkSmartPointer<vtkNamedColors> colors =vtkSmartPointer<vtkNamedColors>::New();

  vtkSmartPointer<vtkPolyDataMapper> meshMapper =vtkSmartPointer<vtkPolyDataMapper>::New();
  meshMapper->SetInputConnection(appendFilter->GetOutputPort());
  /*meshMapper->SetInputData(appendFilter->GetOutput());
  meshMapper->Update();*/

  vtkSmartPointer<vtkActor> meshActor =vtkSmartPointer<vtkActor>::New();
  meshActor->SetMapper(meshMapper);
  meshActor->GetProperty()->SetColor(colors->GetColor3d("Banana").GetData());
  meshActor->GetProperty()->EdgeVisibilityOn();

  

  vtkSmartPointer<vtkVertexGlyphFilter> glyphFilter =vtkSmartPointer<vtkVertexGlyphFilter>::New();
  glyphFilter->SetInputData(appendFilter->GetOutput());

   
  vtkSmartPointer<vtkPolyDataMapper> pointMapper =vtkSmartPointer<vtkPolyDataMapper>::New();
  pointMapper->SetInputConnection(glyphFilter->GetOutputPort());

  vtkSmartPointer<vtkActor> pointActor =vtkSmartPointer<vtkActor>::New();
  pointActor->GetProperty()->SetColor(colors->GetColor3d("Tomato").GetData());
  pointActor->GetProperty()->SetPointSize(5);
  pointActor->SetMapper(pointMapper);

  

  vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();
  vtkSmartPointer<vtkRenderWindow> renderWindow =vtkSmartPointer<vtkRenderWindow>::New();
  renderWindow->AddRenderer(renderer);
  vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =vtkSmartPointer<vtkRenderWindowInteractor>::New();
  renderWindowInteractor->SetRenderWindow(renderWindow);

  renderer->AddActor(meshActor);
  renderer->AddActor(pointActor);

  renderer->SetBackground(colors->GetColor3d("Mint").GetData());

  renderWindow->Render();
  renderWindowInteractor->Start();

  return EXIT_SUCCESS;
}

CMakeLists.txt

cmake_minimum_required(VERSION 3.3 FATAL_ERROR)

project(Delaunay2D)

find_package(VTK COMPONENTS 
  vtkCommonColor
  vtkCommonCore
  vtkCommonDataModel
  vtkFiltersCore
  vtkFiltersGeneral
  vtkInteractionStyle
  vtkRenderingCore
  vtkRenderingFreeType
  vtkRenderingOpenGL2 QUIET)
if (NOT VTK_FOUND)
  message("Skipping Delaunay2D: ${VTK_NOT_FOUND_MESSAGE}")
  return ()
endif()
message (STATUS "VTK_VERSION: ${VTK_VERSION}")
if (VTK_VERSION VERSION_LESS "8.90.0")
  # old sysem
  include(${VTK_USE_FILE})
  add_executable(Delaunay2D MACOSX_BUNDLE Delaunay2D.cxx )
  target_link_libraries(Delaunay2D PRIVATE ${VTK_LIBRARIES})
else ()
  # include all components
  add_executable(Delaunay2D MACOSX_BUNDLE Delaunay2D.cxx )
  target_link_libraries(Delaunay2D PRIVATE ${VTK_LIBRARIES})
  # vtk_module_autoinit is needed
  vtk_module_autoinit(
    TARGETS Delaunay2D
    MODULES ${VTK_LIBRARIES}
    )
endif () 

File Tree

Build Delaunay2D

cd Delaunay2D/build
cmake .. -G "Visual Studio 14 2015 Win64"