VTK 中提取等值面的等值线，使用C++代码实现

在VTK中提取等值面的等值线，可以使用vtkContourFilter类。vtkContourFilter是一种数据过滤器，可以从输入数据集中提取等值线，并生成等值面。以下是使用vtkContourFilter提取等值面的等值线的C++代码实现：

#include <vtkSmartPointer.h>
#include <vtkContourFilter.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkXMLImageDataReader.h>

int main(int argc, char* argv[])
{
    // 读取图像数据
    vtkSmartPointer<vtkXMLImageDataReader> reader =
        vtkSmartPointer<vtkXMLImageDataReader>::New();
    reader->SetFileName("volume.vti");
    reader->Update();

    // 提取等值面的等值线
    vtkSmartPointer<vtkContourFilter> contourFilter =
        vtkSmartPointer<vtkContourFilter>::New();
    contourFilter->SetInputConnection(reader->GetOutputPort());
    contourFilter->SetValue(0, 200); // 等值参数为200

    // 显示等值线
    vtkSmartPointer<vtkPolyDataMapper> mapper =
        vtkSmartPointer<vtkPolyDataMapper>::New();
    mapper->SetInputConnection(contourFilter->GetOutputPort());
    mapper->ScalarVisibilityOff();

    vtkSmartPointer<vtkActor> actor =
        vtkSmartPointer<vtkActor>::New();
    actor->SetMapper(mapper);

    // 创建渲染器和窗口
    vtkSmartPointer<vtkRenderer> renderer =
        vtkSmartPointer<vtkRenderer>::New();
    renderer->AddActor(actor);

    vtkSmartPointer<vtkRenderWindow> renderWindow =
        vtkSmartPointer<vtkRenderWindow>::New();
    renderWindow->AddRenderer(renderer);

    // 创建交互器并启动渲染
    vtkSmartPointer<vtkRenderWindowInteractor> interactor =
        vtkSmartPointer<vtkRenderWindowInteractor>::New();
    interactor->SetRenderWindow(renderWindow);
    interactor->Initialize();
    renderWindow->Render();
    interactor->Start();

    return 0;
}

在上面的代码中，我们首先使用vtkXMLImageDataReader读取了一个VTI格式的图像数据集。然后，我们使用vtkContourFilter提取了等值面的等值线，并将其输出连接到vtkPolyDataMapper的输入端。最后，我们创建了一个vtkActor对象，并将vtkPolyDataMapper设置为其Mapper属性。我们使用vtkRenderWindow和vtkRenderer显示了等值线，并使用vtkRenderWindowInteractor启动了渲染。