第01章 15 vtkDataArray的特点与应用

在 VTK（Visualization Toolkit）中，vtkDataArray 是一个用于存储数组的基类，可以用于存储不同类型的数值数据。虽然 vtkDataArray 本身是一个抽象类，但 VTK 提供了多种具体的实现类，如 vtkFloatArray、vtkDoubleArray、vtkIntArray 等，用于存储不同类型的数值数据。

vtkDataArray 的特点

• 多类型支持：每个 vtkDataArray 可以存储相同类型的数值数据（例如，所有元素都是浮点数），但不同下标的 vtkDataArray 可以是不同类型的单元数据。
• 灵活性：vtkDataArray 可以用于存储点数据、单元数据、场数据等多种属性数据。
• 高效性：VTK 提供了高效的内存管理和访问方法，适用于大规模数据的处理和可视化。

不同下标的 vtkDataArray 可以是不同类型的单元数据

在一个 vtkDataSet 中，不同类型的属性数据可以使用不同类型的 vtkDataArray 来表示。例如，同一 vtkDataSet 中可以包含标量数据、矢量数据和张量数据，这些数据可以分别存储在不同类型的 vtkDataArray 中。

示例应用

以下是一个具体的示例，说明如何在一个 vtkUnstructuredGrid 中使用不同类型的 vtkDataArray 来存储不同类型的属性数据。

#include <vtkNew.h>
#include <vtkUnstructuredGrid.h>
#include <vtkPoints.h>
#include <vtkCellArray.h>
#include <vtkFloatArray.h>
#include <vtkDoubleArray.h>
#include <vtkIntArray.h>

int main() {
    // 创建点和单元
    vtkNew<vtkPoints> points;
    points->InsertNextPoint(0.0, 0.0, 0.0);
    points->InsertNextPoint(1.0, 0.0, 0.0);
    points->InsertNextPoint(0.0, 1.0, 0.0);
    points->InsertNextPoint(1.0, 1.0, 0.0);

    vtkNew<vtkCellArray> cells;
    cells->InsertNextCell(4);
    cells->InsertCellPoint(0);
    cells->InsertCellPoint(1);
    cells->InsertCellPoint(3);
    cells->InsertCellPoint(2);

    // 创建 unstructured grid
    vtkNew<vtkUnstructuredGrid> unstructuredGrid;
    unstructuredGrid->SetPoints(points);
    unstructuredGrid->SetCells(VTK_QUAD, cells);

    // 添加标量数据（浮点型）
    vtkNew<vtkFloatArray> scalarData;
    scalarData->SetName("Temperature");
    scalarData->InsertNextValue(20.0);
    scalarData->InsertNextValue(25.0);
    scalarData->InsertNextValue(22.0);
    scalarData->InsertNextValue(23.0);
    unstructuredGrid->GetPointData()->SetScalars(scalarData);

    // 添加矢量数据（双精度型）
    vtkNew<vtkDoubleArray> vectorData;
    vectorData->SetName("Velocity");
    vectorData->SetNumberOfComponents(3); // 3D vector
    vectorData->InsertNextTuple3(1.0, 0.5, 0.0);
    vectorData->InsertNextTuple3(0.5, 1.0, 0.0);
    vectorData->InsertNextTuple3(0.0, 1.0, 0.5);
    vectorData->InsertNextTuple3(0.0, 0.5, 1.0);
    unstructuredGrid->GetPointData()->SetVectors(vectorData);

    // 添加张量数据（整型）
    vtkNew<vtkIntArray> tensorData;
    tensorData->SetName("Stress");
    tensorData->SetNumberOfComponents(6); // Symmetric tensor in 3D
    tensorData->InsertNextTuple6(1, 0, 0, 0, 0, 0);
    tensorData->InsertNextTuple6(0, 1, 0, 0, 0, 0);
    tensorData->InsertNextTuple6(0, 0, 1, 0, 0, 0);
    tensorData->InsertNextTuple6(0, 0, 0, 1, 0, 0);
    unstructuredGrid->GetPointData()->SetTensors(tensorData);

    // 输出数据
    unstructuredGrid->PrintSelf(std::cout, 0);

    return 0;
}

解释

1. 创建点和单元：首先创建点（vtkPoints）和单元（vtkCellArray），并定义一个四边形单元。
2. 创建 vtkUnstructuredGrid：使用创建的点和单元构建 vtkUnstructuredGrid。
3. 添加标量数据：使用 vtkFloatArray 存储标量数据（温度），并将其添加到点数据中。
4. 添加矢量数据：使用 vtkDoubleArray 存储矢量数据（速度），并将其添加到点数据中。
5. 添加张量数据：使用 vtkIntArray 存储张量数据（应力），并将其添加到点数据中。
6. 输出数据：最后，打印 vtkUnstructuredGrid 的信息，以验证数据的添加。

应用场合

• 多物理场分析：在一个 vtkUnstructuredGrid 中，可以同时存储温度（标量）、速度（矢量）和应力（张量）等不同类型的物理场数据。
• 多属性可视化：不同的属性数据可以分别用不同类型的 vtkDataArray 存储，便于进行多属性的可视化和分析。
• 复杂数据集处理：例如在流体力学模拟中，可以同时存储速度场、压力场和应力张量等不同类型的属性数据。

通过这种方式，VTK 提供了灵活的数据存储和管理机制，能够满足不同领域的需求。