VTK：对输入的三维模型在某个方向等间距提取模型的切面轮廓线

测试的三维模型为bunny，即斯坦福兔子

代码示例：

#include <vtkOBJReader.h>#include <vtkSmartPointer.h>#include <vtkCutter.h>#include <vtkPolyDataMapper.h>#include <vtkPlane.h>#include <vtkProperty.h>#include <vtkActor.h>#include <vtkRenderer.h>#include <vtkRenderWindow.h>#include <vtkRenderWindowInteractor.h>#include <vtkSphereSource.h>#include <vtkPolyDataWriter.h>#include <vtkPolyData.h>#include <vtkStripper.h>#include <vtkContourFilter.h>#include <vtkPoints.h>#include <iostream>using namespace std;double Distance(double *a,double *b){double o=(a[0]-b[0])*(a[0]-b[0]);double p=(a[1]-b[1])*(a[1]-b[1]);double q=(a[2]-b[2])*(a[2]-b[2]);return o+p+q;}int main(){//读取obj文件vtkSmartPointer<vtkOBJReader> reader =vtkSmartPointer<vtkOBJReader>::New();  reader->SetFileName("C:\\Users\\HuangWang\\Desktop\\bunny-chuli.obj");  reader->Update();vtkSmartPointer<vtkPolyData> inputPolyData=vtkSmartPointer<vtkPolyData>::New();inputPolyData = reader->GetOutput();vtkSmartPointer<vtkPolyDataMapper> inputMapper =vtkSmartPointer<vtkPolyDataMapper>::New();inputMapper->SetInput(inputPolyData);//创建切割平面vtkSmartPointer<vtkPlane> plane =vtkSmartPointer<vtkPlane>::New();plane->SetOrigin(inputPolyData->GetCenter());//设置切割平面起点plane->SetNormal(1,0,0);//设置切割方向为X方向//得到输入的obj模型的最小坐标double minBound[3];minBound[0] = inputPolyData->GetBounds()[0];minBound[1] = inputPolyData->GetBounds()[2];minBound[2] = inputPolyData->GetBounds()[4];//得到输入的obj模型的最小大坐标double maxBound[3];maxBound[0] = inputPolyData->GetBounds()[1];maxBound[1] = inputPolyData->GetBounds()[3];maxBound[2] = inputPolyData->GetBounds()[5];//得到输入的obj模型的中心坐标double center[3];center[0] = inputPolyData->GetCenter()[0];center[1] = inputPolyData->GetCenter()[1];center[2] = inputPolyData->GetCenter()[2];double distanceMin=sqrt(Distance(minBound,center));double distanceMax=sqrt(Distance(maxBound,center));//创建模型切割器vtkSmartPointer<vtkCutter> cutter =vtkSmartPointer<vtkCutter>::New();cutter->SetCutFunction(plane);//设置切割平面cutter->SetInput(inputPolyData);//设置模型cutter->GenerateValues(30, -distanceMin, distanceMax);//在模型的最大最小范围内等间距创建30个切面，得到轮廓线//将切线结果输出为vtk文件格式vtkSmartPointer<vtkPolyData> ResultPoly=cutter->GetOutput();//输出为polydatavtkSmartPointer<vtkPolyDataWriter> vtkWriter = vtkSmartPointer<vtkPolyDataWriter>::New();vtkWriter->SetInput(cutter->GetOutput());vtkWriter->SetFileName("C:\\Users\\HuangWang\\Desktop\\result.vtk");vtkWriter->Write();vtkSmartPointer<vtkPolyDataMapper> cutterMapper =vtkSmartPointer<vtkPolyDataMapper>::New();cutterMapper->SetInputConnection( cutter->GetOutputPort());cutterMapper->ScalarVisibilityOff();vtkSmartPointer<vtkActor> planeActor =vtkSmartPointer<vtkActor>::New();planeActor->GetProperty()->SetColor(1,0,0);planeActor->GetProperty()->SetLineWidth(5);planeActor->SetMapper(cutterMapper);vtkSmartPointer<vtkActor> inputActor =vtkSmartPointer<vtkActor>::New();inputActor->GetProperty()->SetColor(0,1,0);inputActor->SetMapper(inputMapper);vtkSmartPointer<vtkRenderer> renderer =vtkSmartPointer<vtkRenderer>::New();renderer->AddActor(planeActor); renderer->AddActor(inputActor);vtkSmartPointer<vtkRenderWindow> renderWindow =vtkSmartPointer<vtkRenderWindow>::New();renderWindow->AddRenderer(renderer);renderWindow->SetSize(600, 600);vtkSmartPointer<vtkRenderWindowInteractor> interactor =vtkSmartPointer<vtkRenderWindowInteractor>::New();interactor->SetRenderWindow(renderWindow);renderer->SetBackground(0,0,0);renderWindow->Render();interactor->Start();return 0;}

结果示例：

最后输出的结果VTK文件包含了每一条轮廓线的点数据