DeformPointSet
VTKExamples/Cxx/Meshes/DeformPointSet
Description¶
Deform a point set with a polyhedra. Supply a .vtp file and see it deformed.
Code¶
DeformPointSet.cxx
#include "vtkSmartPointer.h" #include "vtkActor.h" #include "vtkPolyData.h" #include "vtkPolyDataMapper.h" #include "vtkRenderWindow.h" #include "vtkRenderWindowInteractor.h" #include "vtkRenderer.h" #include "vtkRenderer.h" #include "vtkSphereSource.h" #include "vtkElevationFilter.h" #include "vtkProperty.h" #include "vtkCellArray.h" #include "vtkPoints.h" #include "vtkDeformPointSet.h" #include "vtkCamera.h" #include "vtkXMLPolyDataReader.h" int main( int argc, char *argv[] ) { vtkSmartPointer<vtkPolyData> input = vtkSmartPointer<vtkPolyData>::New(); double bounds[6]; if(argc == 1) { // Create a sphere to warp vtkSmartPointer<vtkSphereSource> sphere = vtkSmartPointer<vtkSphereSource>::New(); sphere->SetThetaResolution(51); sphere->SetPhiResolution(17); sphere->Update(); sphere->GetOutput()->GetBounds(bounds); // Generate some scalars on the polydata vtkSmartPointer<vtkElevationFilter> ele = vtkSmartPointer<vtkElevationFilter>::New(); ele->SetInputConnection(sphere->GetOutputPort()); ele->SetLowPoint(0,0,-0.5); ele->SetHighPoint(0,0,0.5); ele->SetLowPoint((bounds[1] + bounds[0]) / 2.0, (bounds[3] + bounds[2]) / 2.0, -bounds[5]); ele->SetHighPoint((bounds[1] + bounds[0]) / 2.0, (bounds[3] + bounds[2]) / 2.0, bounds[5]); ele->Update(); input->ShallowCopy(ele->GetOutput()); } else { std::string inputFilename = argv[1]; vtkSmartPointer<vtkXMLPolyDataReader> reader = vtkSmartPointer<vtkXMLPolyDataReader>::New(); reader->SetFileName(inputFilename.c_str()); reader->Update(); input->ShallowCopy(reader->GetOutput()); input->GetBounds(bounds); } // Now create a control mesh, in this case a octagon that encloses // the point set vtkSmartPointer<vtkPoints> pts = vtkSmartPointer<vtkPoints>::New(); pts->SetNumberOfPoints(6); pts->SetPoint(0, bounds[0] - .1 * (bounds[1] - bounds[0]), (bounds[3] + bounds[2]) / 2.0, (bounds[5] + bounds[4]) / 2.0); pts->SetPoint(1, bounds[1] + .1 * (bounds[1] - bounds[0]), (bounds[3] + bounds[2]) / 2.0, (bounds[5] + bounds[4]) / 2.0); pts->SetPoint(2, (bounds[1] + bounds[0]) / 2.0, bounds[2] - .1 * (bounds[3] - bounds[2]), (bounds[5] + bounds[4]) / 2.0); pts->SetPoint(3, (bounds[1] + bounds[0]) / 2.0, bounds[3] + .1 * (bounds[3] - bounds[2]), (bounds[5] + bounds[4]) / 2.0); pts->SetPoint(4, (bounds[1] + bounds[0]) / 2.0, (bounds[3] + bounds[2]) / 2.0, bounds[4] - .1 * (bounds[5] - bounds[4])); pts->SetPoint(5, (bounds[1] + bounds[0]) / 2.0, (bounds[3] + bounds[2]) / 2.0, bounds[5] + .1 * (bounds[5] - bounds[4])); vtkSmartPointer<vtkCellArray> tris = vtkSmartPointer<vtkCellArray>::New(); tris->InsertNextCell(3); tris->InsertCellPoint(2); tris->InsertCellPoint(0); tris->InsertCellPoint(4); tris->InsertNextCell(3); tris->InsertCellPoint(1); tris->InsertCellPoint(2); tris->InsertCellPoint(4); tris->InsertNextCell(3); tris->InsertCellPoint(3); tris->InsertCellPoint(1); tris->InsertCellPoint(4); tris->InsertNextCell(3); tris->InsertCellPoint(0); tris->InsertCellPoint(3); tris->InsertCellPoint(4); tris->InsertNextCell(3); tris->InsertCellPoint(0); tris->InsertCellPoint(2); tris->InsertCellPoint(5); tris->InsertNextCell(3); tris->InsertCellPoint(2); tris->InsertCellPoint(1); tris->InsertCellPoint(5); tris->InsertNextCell(3); tris->InsertCellPoint(1); tris->InsertCellPoint(3); tris->InsertCellPoint(5); tris->InsertNextCell(3); tris->InsertCellPoint(3); tris->InsertCellPoint(0); tris->InsertCellPoint(5); vtkSmartPointer<vtkPolyData> pd = vtkSmartPointer<vtkPolyData>::New(); pd->SetPoints(pts); pd->SetPolys(tris); // Display the control mesh vtkSmartPointer<vtkPolyDataMapper> meshMapper = vtkSmartPointer<vtkPolyDataMapper>::New(); meshMapper->SetInputData(pd); vtkSmartPointer<vtkActor> meshActor = vtkSmartPointer<vtkActor>::New(); meshActor->SetMapper(meshMapper); meshActor->GetProperty()->SetRepresentationToWireframe(); meshActor->GetProperty()->SetColor(0,0,0); // Do the intitial weight generation vtkSmartPointer<vtkDeformPointSet> deform = vtkSmartPointer<vtkDeformPointSet>::New(); deform->SetInputData(input); deform->SetControlMeshData(pd); deform->Update(); // this creates the initial weights // Now move one point and deform double controlPoint[3]; pts->GetPoint(5, controlPoint); pts->SetPoint(5, controlPoint[0], controlPoint[1], bounds[5] + .8 * (bounds[5] - bounds[4])); pts->Modified(); // Display the warped polydata vtkSmartPointer<vtkPolyDataMapper> polyMapper = vtkSmartPointer<vtkPolyDataMapper>::New(); polyMapper->SetInputConnection(deform->GetOutputPort()); vtkSmartPointer<vtkActor> polyActor = vtkSmartPointer<vtkActor>::New(); polyActor->SetMapper(polyMapper); vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New(); vtkSmartPointer<vtkRenderWindow> renWin = vtkSmartPointer<vtkRenderWindow>::New(); renWin->AddRenderer(renderer); vtkSmartPointer<vtkRenderWindowInteractor> iren = vtkSmartPointer<vtkRenderWindowInteractor>::New(); iren->SetRenderWindow(renWin); renderer->AddActor(polyActor); renderer->AddActor(meshActor); renderer->GetActiveCamera()->SetPosition(1,1,1); renderer->ResetCamera(); renderer->SetBackground(.2, .3, .4); renWin->SetSize(300,300); renWin->Render(); iren->Start(); return EXIT_SUCCESS; }
CMakeLists.txt¶
cmake_minimum_required(VERSION 2.8) PROJECT(DeformPointSet) find_package(VTK REQUIRED) include(${VTK_USE_FILE}) add_executable(DeformPointSet MACOSX_BUNDLE DeformPointSet.cxx) target_link_libraries(DeformPointSet ${VTK_LIBRARIES})
Download and Build DeformPointSet¶
Danger
The generation of tar files has not been ported to the new VTKExamples. Some tarballs may be missing or out-of-date.
Click here to download DeformPointSet and its CMakeLists.txt file. Once the tarball DeformPointSet.tar has been downloaded and extracted,
cd DeformPointSet/build
If VTK is installed:
cmake ..
If VTK is not installed but compiled on your system, you will need to specify the path to your VTK build:
cmake -DVTK_DIR:PATH=/home/me/vtk_build ..
Build the project:
make
and run it:
./DeformPointSet
WINDOWS USERS PLEASE NOTE: Be sure to add the VTK bin directory to your path. This will resolve the VTK dll's at run time.