第06章 12 VTK同时绘制几何信息及体绘制

VTK体绘制同时绘制几何信息及体绘制的用途

VTK（Visualization Toolkit）是一个强大的开源库，广泛应用于科学可视化领域。它提供了多种绘制技术，其中包括体绘制（Volume Rendering）和几何绘制（Geometry Rendering）。体绘制用于直接渲染三维数据体，而几何绘制则用于渲染表面几何信息，如等值面（Isosurface）。

用途

1. 医学影像可视化：在医学领域，CT、MRI等扫描数据通常包含大量的体数据。通过体绘制，可以直接观察到器官、骨骼等的内部结构。同时，几何绘制可以用于显示解剖结构的表面信息，如皮肤、骨骼轮廓等。这种方式可以帮助医生更全面地分析和诊断病情。

2. 工程与科学数据可视化：在工程领域，体绘制可以用于显示材料内部的应力、温度、密度等分布情况。几何绘制则可以用于显示模型的外部轮廓或特定等值面的形状，帮助工程师理解复杂的三维数据。

3. 教育与科研：通过同时显示体绘制和几何信息，研究人员可以更直观地理解数据的内部和外部特征，这对于教育和科研具有重要意义。

C++示例：同时绘制膝盖部骨头体绘制及剖分表面皮肤（使用等值面）

以下是一个简单的C++示例，展示如何使用VTK同时进行体绘制和几何绘制。假设我们有一个膝盖的DICOM数据集，我们希望同时显示骨头的体绘制和皮肤的等值面。

#include <vtkSmartPointer.h>
#include <vtkDICOMImageReader.h>
#include <vtkVolume.h>
#include <vtkVolumeProperty.h>
#include <vtkColorTransferFunction.h>
#include <vtkPiecewiseFunction.h>
#include <vtkFixedPointVolumeRayCastMapper.h>
#include <vtkMarchingCubes.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkInteractorStyleTrackballCamera.h>

int main(int argc, char *argv[])
{
    // 读取DICOM数据
    vtkSmartPointer<vtkDICOMImageReader> reader = vtkSmartPointer<vtkDICOMImageReader>::New();
    reader->SetDirectoryName("path_to_dicom_data");
    reader->Update();

    // 体绘制Mapper
    vtkSmartPointer<vtkFixedPointVolumeRayCastMapper> volumeMapper = vtkSmartPointer<vtkFixedPointVolumeRayCastMapper>::New();
    volumeMapper->SetInputConnection(reader->GetOutputPort());

    // 体绘制属性
    vtkSmartPointer<vtkVolumeProperty> volumeProperty = vtkSmartPointer<vtkVolumeProperty>::New();
    vtkSmartPointer<vtkColorTransferFunction> colorTransferFunction = vtkSmartPointer<vtkColorTransferFunction>::New();
    colorTransferFunction->AddRGBPoint(0, 0.0, 0.0, 0.0);
    colorTransferFunction->AddRGBPoint(500, 1.0, 0.0, 0.0);
    colorTransferFunction->AddRGBPoint(1000, 1.0, 1.0, 0.0);
    colorTransferFunction->AddRGBPoint(1500, 0.0, 1.0, 0.0);
    colorTransferFunction->AddRGBPoint(2000, 0.0, 0.0, 1.0);

    vtkSmartPointer<vtkPiecewiseFunction> opacityTransferFunction = vtkSmartPointer<vtkPiecewiseFunction>::New();
    opacityTransferFunction->AddPoint(0, 0.0);
    opacityTransferFunction->AddPoint(500, 0.1);
    opacityTransferFunction->AddPoint(1000, 0.3);
    opacityTransferFunction->AddPoint(1500, 0.6);
    opacityTransferFunction->AddPoint(2000, 1.0);

    volumeProperty->SetColor(colorTransferFunction);
    volumeProperty->SetScalarOpacity(opacityTransferFunction);
    volumeProperty->ShadeOn();

    // 体绘制Actor
    vtkSmartPointer<vtkVolume> volume = vtkSmartPointer<vtkVolume>::New();
    volume->SetMapper(volumeMapper);
    volume->SetProperty(volumeProperty);

    // 等值面提取（皮肤）
    vtkSmartPointer<vtkMarchingCubes> skinExtractor = vtkSmartPointer<vtkMarchingCubes>::New();
    skinExtractor->SetInputConnection(reader->GetOutputPort());
    skinExtractor->SetValue(0, 100); // 调整等值面的阈值以适应皮肤

    // 几何绘制Mapper
    vtkSmartPointer<vtkPolyDataMapper> skinMapper = vtkSmartPointer<vtkPolyDataMapper>::New();
    skinMapper->SetInputConnection(skinExtractor->GetOutputPort());
    skinMapper->ScalarVisibilityOff();

    // 几何绘制Actor
    vtkSmartPointer<vtkActor> skinActor = vtkSmartPointer<vtkActor>::New();
    skinActor->SetMapper(skinMapper);

    // 创建Renderer
    vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();
    renderer->AddVolume(volume);
    renderer->AddActor(skinActor);
    renderer->SetBackground(0.1, 0.2, 0.3);

    // 创建RenderWindow
    vtkSmartPointer<vtkRenderWindow> renderWindow = vtkSmartPointer<vtkRenderWindow>::New();
    renderWindow->AddRenderer(renderer);

    // 创建RenderWindowInteractor
    vtkSmartPointer<vtkRenderWindowInteractor> interactor = vtkSmartPointer<vtkRenderWindowInteractor>::New();
    interactor->SetRenderWindow(renderWindow);

    // 设置交互风格
    vtkSmartPointer<vtkInteractorStyleTrackballCamera> style = vtkSmartPointer<vtkInteractorStyleTrackballCamera>::New();
    interactor->SetInteractorStyle(style);

    // 启动渲染和交互
    renderWindow->Render();
    interactor->Start();

    return 0;
}

代码说明

1. 数据读取：使用vtkDICOMImageReader读取DICOM数据。
2. 体绘制：使用vtkFixedPointVolumeRayCastMapper进行体绘制，并通过vtkVolumeProperty设置颜色和透明度。
3. 等值面提取：使用vtkMarchingCubes提取皮肤的等值面。
4. 几何绘制：使用vtkPolyDataMapper和vtkActor将等值面渲染为几何表面。
5. 渲染和交互：使用vtkRenderer、vtkRenderWindow和vtkRenderWindowInteractor进行渲染和交互。

总结

通过VTK，我们可以轻松地将体绘制和几何绘制结合在一起，创建出丰富的三维可视化效果。这个示例展示了如何在同一窗口中同时显示骨头的体绘制和皮肤的等值面，适用于医学影像等领域的应用。