VTK 实例60：等值面提取（表面重建）

 1 #include <vtkAutoInit.h>
 2 VTK_MODULE_INIT(vtkRenderingOpenGL2);
 3 VTK_MODULE_INIT(vtkRenderingFreeType);
 4 VTK_MODULE_INIT(vtkInteractionStyle);
 5 
 6 #include <vtkSmartPointer.h>
 7 #include <vtkMarchingCubes.h>
 8 #include <vtkContourFilter.h>
 9 #include <vtkVoxelModeller.h>
10 #include <vtkSphereSource.h>
11 #include <vtkImageData.h>
12 #include <vtkMetaImageReader.h>
13 #include <vtkInteractorStyleImage.h>
14 #include <vtkActor.h>
15 #include <vtkPolyDataMapper.h>
16 #include <vtkRenderWindowInteractor.h>
17 #include <vtkRenderWindow.h>
18 #include <vtkRenderer.h>
19 #include <vtkImageActor.h>
20 #include <vtkProperty.h>
21 
22 //测试：../data/HeadMRVolume.mhd 200
23 int main(int argc, char *argv[])
24 {
25     vtkSmartPointer<vtkMetaImageReader> reader =
26         vtkSmartPointer<vtkMetaImageReader>::New();
27     reader->SetFileName("C:\\Users\\Administrator\\Desktop\\VTK2\\hellovtk\\vtk_图像处理学习\\第六章_图像处理\\data\\HeadMRVolume.mhd");
28     reader->Update();
29 
30     double isoValue = atof("200");//设置等值面值为200
31 
32     vtkSmartPointer<vtkMarchingCubes> surface = 
33         vtkSmartPointer<vtkMarchingCubes>::New();
34     surface->SetInputData(reader->GetOutput());
35     surface->ComputeNormalsOn();
36     surface->SetValue(0, isoValue);
37     //surface->GenerateValues(5, 150,200);
38 
39     /*vtkSmartPointer<vtkContourFilter> surface = 
40     vtkSmartPointer<vtkContourFilter>::New();
41     surface->SetInput(reader->GetOutput());
42     surface->ComputeNormalsOn();
43     surface->SetValue(0, isoValue);*/
44 
45     vtkSmartPointer<vtkPolyDataMapper> surfMapper = 
46         vtkSmartPointer<vtkPolyDataMapper>::New();
47     surfMapper->SetInputConnection(surface->GetOutputPort());
48     surfMapper->ScalarVisibilityOff();
49 
50     vtkSmartPointer<vtkActor> surfActor = 
51         vtkSmartPointer<vtkActor>::New();
52     surfActor->SetMapper(surfMapper);
53     surfActor->GetProperty()->SetColor(1.0, 0.0, 0.0);
54 
55     vtkSmartPointer<vtkRenderer> renderer = 
56         vtkSmartPointer<vtkRenderer>::New();
57     renderer->AddActor(surfActor);
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("PolyDataMarchingCubes");
66     renderWindow->Render();
67 
68     vtkSmartPointer<vtkRenderWindowInteractor> interactor = 
69         vtkSmartPointer<vtkRenderWindowInteractor>::New();
70     interactor->SetRenderWindow(renderWindow);
71 
72     renderWindow->Render();
73     interactor->Start();
74 
75     return EXIT_SUCCESS;
76 }