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

#include <vtkAutoInit.h>
VTK_MODULE_INIT(vtkRenderingOpenGL2);
VTK_MODULE_INIT(vtkRenderingFreeType);
VTK_MODULE_INIT(vtkInteractionStyle);

#include <vtkSmartPointer.h>
#include <vtkMarchingCubes.h>
#include <vtkContourFilter.h>
#include <vtkVoxelModeller.h>
#include <vtkSphereSource.h>
#include <vtkImageData.h>
#include <vtkMetaImageReader.h>
#include <vtkInteractorStyleImage.h>
#include <vtkActor.h>
#include <vtkPolyDataMapper.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkRenderWindow.h>
#include <vtkRenderer.h>
#include <vtkImageActor.h>
#include <vtkProperty.h>

//测试：../data/HeadMRVolume.mhd 200
int main(int argc, char *argv[])
{
    vtkSmartPointer<vtkMetaImageReader> reader =
        vtkSmartPointer<vtkMetaImageReader>::New();
    reader->SetFileName("C:\\Users\\Administrator\\Desktop\\VTK2\\hellovtk\\vtk_图像处理学习\\第六章_图像处理\\data\\HeadMRVolume.mhd");
    reader->Update();

    double isoValue = atof("200");//设置等值面值为200

    vtkSmartPointer<vtkMarchingCubes> surface = 
        vtkSmartPointer<vtkMarchingCubes>::New();
    surface->SetInputData(reader->GetOutput());
    surface->ComputeNormalsOn();
    surface->SetValue(0, isoValue);
    //surface->GenerateValues(5, 150,200);

    /*vtkSmartPointer<vtkContourFilter> surface = 
    vtkSmartPointer<vtkContourFilter>::New();
    surface->SetInput(reader->GetOutput());
    surface->ComputeNormalsOn();
    surface->SetValue(0, isoValue);*/

    vtkSmartPointer<vtkPolyDataMapper> surfMapper = 
        vtkSmartPointer<vtkPolyDataMapper>::New();
    surfMapper->SetInputConnection(surface->GetOutputPort());
    surfMapper->ScalarVisibilityOff();

    vtkSmartPointer<vtkActor> surfActor = 
        vtkSmartPointer<vtkActor>::New();
    surfActor->SetMapper(surfMapper);
    surfActor->GetProperty()->SetColor(1.0, 0.0, 0.0);

    vtkSmartPointer<vtkRenderer> renderer = 
        vtkSmartPointer<vtkRenderer>::New();
    renderer->AddActor(surfActor);
    renderer->SetBackground(1.0, 1.0, 1.0);

    vtkSmartPointer<vtkRenderWindow> renderWindow = 
        vtkSmartPointer<vtkRenderWindow>::New();
    renderWindow->AddRenderer(renderer);
    renderWindow->SetSize(640, 480);
    renderWindow->Render();
    renderWindow->SetWindowName("PolyDataMarchingCubes");
    renderWindow->Render();

    vtkSmartPointer<vtkRenderWindowInteractor> interactor = 
        vtkSmartPointer<vtkRenderWindowInteractor>::New();
    interactor->SetRenderWindow(renderWindow);

    renderWindow->Render();
    interactor->Start();

    return EXIT_SUCCESS;
}