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"