VTK 实例49:计算三角网络模型面积、体积

发布时间 2023-08-16 08:34:47作者: 一杯清酒邀明月

 1 #include "vtkAutoInit.h" 
 2 VTK_MODULE_INIT(vtkRenderingOpenGL2);
 3 VTK_MODULE_INIT(vtkInteractionStyle);
 4 
 5 #include <vtkSmartPointer.h>
 6 #include <vtkPointData.h>
 7 #include <vtkPolyDataMapper.h>
 8 #include <vtkActor.h>
 9 #include <vtkScalarBarActor.h>
10 #include <vtkRenderWindow.h>
11 #include <vtkRenderer.h>
12 #include <vtkRenderWindowInteractor.h>
13 #include <vtkPolyData.h>
14 #include <vtkCubeSource.h>
15 #include <vtkProperty.h>
16 #include <vtkTriangleFilter.h>
17 #include <vtkMassProperties.h>
18 
19 int main(int argc, char *argv[])
20 {
21     vtkSmartPointer<vtkCubeSource> cubeSource = 
22         vtkSmartPointer<vtkCubeSource>::New();
23     cubeSource->Update();
24 
25     vtkSmartPointer<vtkTriangleFilter> triFilter = 
26         vtkSmartPointer<vtkTriangleFilter>::New();
27     triFilter->SetInputData(cubeSource->GetOutput());
28     triFilter->Update();
29 
30     vtkSmartPointer<vtkMassProperties> massProp = 
31         vtkSmartPointer<vtkMassProperties>::New();
32     massProp->SetInputData(triFilter->GetOutput());
33     float vol = massProp->GetVolume();
34     float area= massProp->GetSurfaceArea();
35     float maxArea = massProp->GetMaxCellArea();
36     float minArea = massProp->GetMinCellArea();
37 
38     std::cout<<"Volume      :"<<vol<<std::endl;
39     std::cout<<"Surface Area:"<<area<<std::endl;
40     std::cout<<"Max Area    :"<<maxArea<<std::endl;
41     std::cout<<"Min Area    :"<<minArea<<std::endl;
42 
43     vtkSmartPointer<vtkPolyDataMapper> mapper = 
44         vtkSmartPointer<vtkPolyDataMapper>::New();
45     //mapper->SetInput(cubeSource->GetOutput());
46     mapper->SetInputData(triFilter->GetOutput());
47 
48     vtkSmartPointer<vtkActor> actor = 
49         vtkSmartPointer<vtkActor>::New();
50     actor->SetMapper(mapper);
51     actor->GetProperty()->SetColor(0,1,0);
52     actor->GetProperty()->SetEdgeColor(1,0,0);
53     actor->GetProperty()->SetEdgeVisibility(1);
54 
55     vtkSmartPointer<vtkRenderer> renderer = 
56         vtkSmartPointer<vtkRenderer>::New();
57     renderer->AddActor(actor);
58     renderer->SetBackground(1.0,1.0,1.0);
59 
60     vtkSmartPointer<vtkRenderWindow> renderWindow = 
61         vtkSmartPointer<vtkRenderWindow>::New();
62     renderWindow->AddRenderer(renderer);
63     renderWindow->SetSize( 640, 480 );
64     renderWindow->Render();
65     renderWindow->SetWindowName("PolyDataMassProperty");
66 
67     vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor = 
68         vtkSmartPointer<vtkRenderWindowInteractor>::New();
69     renderWindowInteractor->SetRenderWindow(renderWindow);
70 
71     renderWindow->Render();
72     renderWindowInteractor->Start();
73 
74     return EXIT_SUCCESS;
75 }