VTK入门：vtkUnstructuredGrid——非结构化网格的“自由组合乐高盒”

原创已于 2025-10-24 16:25:31 修改 · 718 阅读

8 ·

CC 4.0 BY-SA版权

文章标签：

#算法 #VTK #1024程序员节

于 2025-10-24 16:24:00 首次发布

VTK 同时被 2 个专栏收录

47 篇文章

订阅专栏

网格

10 篇文章

订阅专栏

VTK入门：vtkUnstructuredGrid——非结构化网格的“自由组合乐高盒”

如果你之前接触过VTK的vtkImageData（规则体素）或vtkPolyData（离散面），会发现它们的网格都是“整齐排列”的——要么是像魔方一样的规则体素，要么是像拼图一样的离散多边形。但现实世界中，很多几何体是“不规则”的：比如带凹槽的机械零件、起伏的地质模型、复杂的人体器官（如心脏血管），这些用规则网格根本装不下。这时候，vtkUnstructuredGrid就该登场了——它像一个“自由组合的乐高盒”，里面的“积木”（VTK里叫“细胞”）可以是三角形、四面体、六面体等任意形状，还能随意组合，专门用来装非结构化的复杂网格。
在这里插入图片描述

一、先搞懂：vtkUnstructuredGrid到底是什么？

用一句话总结：vtkUnstructuredGrid是VTK中专门处理“非结构化网格”的数据集，支持任意组合的细胞类型（0D点、1D线、2D面、3D体），没有固定的排列规则，完全根据几何体的形状灵活定义。

它和之前学的“结构化数据集”（如vtkImageData）的区别，就像“定制积木”和“标准魔方”的区别：

特性	vtkUnstructuredGrid（非结构化）	vtkImageData（结构化）
网格排列	无固定规则，灵活适配形状	3D规则排列（体素）
细胞类型	支持所有VTK细胞（三角形、四面体等）	仅六面体（体素）
适用场景	复杂不规则几何体	规则体数据（如CT/MRI）
拓扑关系管理	需手动建立（如点-细胞关联）	自动计算（索引即可定位）

简单说：只要你的几何体“不规整”，比如机械零件的孔、地质模型的断层，就用vtkUnstructuredGrid。

二、vtkUnstructuredGrid的“核心零件”：3个关键组成

非结构化网格之所以能“自由组合”，全靠三个核心部分协同工作——就像乐高的“积木形状”“拼接说明”“特殊连接件”，缺一不可。

组成部分	通俗解释	作用举例
Connectivity（连接数组）	记录“每个细胞由哪些点组成”的“拼接说明”	三角形细胞由点0、1、2组成，就记在这里
Types（细胞类型数组）	标记“每个细胞是什么形状”的“积木标签”	标记细胞0是三角形（VTK_TRIANGLE）、细胞1是四面体（VTK_TETRA）
Faces（面数组）	专门为“多面体”准备的“特殊连接件”，记录多面体的面信息	五面体由5个面组成，每个面的点连接关系存在这里

此外，还有两个辅助部分新手也需要了解：

CellLocations（细胞位置数组）：快速定位每个细胞在Connectivity中的起始位置，比如细胞0从索引0开始，细胞1从索引3开始，提高查询效率；
Links（拓扑链接）：记录“每个点被哪些细胞使用”，比如点0被细胞0和细胞2使用，用于快速拓扑查询（如找相邻细胞）。

三、新手必学：vtkUnstructuredGrid的4个常用功能

入门阶段，不用掌握所有功能，先学会这4个最常用的，就能应对大部分基础场景。

1. 内存预分配（避免频繁扩容）

非结构化网格的细胞数量通常很大（比如几十万个），如果每次插入细胞都临时扩容，会严重影响效率。vtkUnstructuredGrid提供两种内存分配方法：

AllocateEstimate(numCells, maxCellSize)：估算分配，适合知道“大概有多少细胞”和“单个细胞最多用多少点”的场景，比如预计1000个细胞，每个最多用4个点（四面体）；
AllocateExact(numCells, connectivitySize)：精确分配，适合知道“确切细胞数”和“总连接点数量”的场景，比如1000个细胞，总共有3000个连接点。

示例（估算分配）：

// 实例化非结构化网格
vtkSmartPointer<vtkUnstructuredGrid> ug = vtkSmartPointer<vtkUnstructuredGrid>::New();
// 估算分配：1000个细胞，每个最多4个点（如四面体）
ug->AllocateEstimate(1000, 4);

2. 插入细胞（组装“乐高积木”）

插入细胞是核心操作，需要指定“细胞形状”和“组成该细胞的点索引”。VTK支持的细胞类型很多，新手先记住常用的几种：

细胞类型	VTK宏定义	点数量	用途举例
三角形	VTK_TRIANGLE	3	2D不规则面
四面体	VTK_TETRA	4	3D不规则体（如地质模型）
六面体	VTK_HEXAHEDRON	8	3D规则块（如机械零件局部）
线段	VTK_LINE	2	1D曲线（如血管中心线）

示例（插入一个三角形和一个四面体）：

// 1. 设置点（先定义3D空间中的点坐标）
vtkSmartPointer<vtkPoints> points = vtkSmartPointer<vtkPoints>::New();
// 点0：(0,0,0)，点1：(1,0,0)，点2：(0,1,0)，点3：(0,0,1)
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(0.0, 0.0, 1.0);
ug->SetPoints(points); // 把点传给非结构化网格

// 2. 插入三角形细胞（由点0、1、2组成）
vtkIdType triPoints[3] = {0, 1, 2};
ug->InsertNextCell(VTK_TRIANGLE, 3, triPoints);

// 3. 插入四面体细胞（由点0、1、2、3组成）
vtkIdType tetraPoints[4] = {0, 1, 2, 3};
ug->InsertNextCell(VTK_TETRA, 4, tetraPoints);

3. 拓扑查询（找“积木的邻居”）

有时候需要知道“某个点被哪些细胞使用”或“某个细胞的邻居是谁”，这就需要先建立拓扑链接（BuildLinks()），再进行查询：

BuildLinks()：建立“点→细胞”的关联关系，必须在查询前调用；
GetPointCells(ptId, cellIds)：查询“点ptId被哪些细胞使用”，结果存在cellIds中。

示例（查询点0被哪些细胞使用）：

// 1. 建立拓扑链接（必须先调用）
ug->BuildLinks();

// 2. 查询点0关联的细胞
vtkSmartPointer<vtkIdList> cellIds = vtkSmartPointer<vtkIdList>::New();
ug->GetPointCells(0, cellIds);

// 3. 输出结果（应该是细胞0和细胞1）
std::cout << "点0被以下细胞使用：";
for (int i = 0; i < cellIds->GetNumberOfIds(); i++) {
    std::cout << cellIds->GetId(i) << " ";
}
// 输出：点0被以下细胞使用：0 1

4. 多面体支持（处理复杂形状）

对于更复杂的形状（如五面体、八面体），普通细胞类型不够用，vtkUnstructuredGrid支持“多面体（VTK_POLYHEDRON）”，需要额外指定“面信息”——比如五面体由5个面组成，每个面由哪些点组成。

示例（简单多面体插入，了解即可）：

// 多面体需要指定点、面的连接关系，步骤稍复杂，入门阶段可先跳过，后续深入
vtkIdType faceStream[] = {5, // 多面体有5个面
                          4, 0,1,2,3, // 第1个面：4个点（四边形）
                          4, 0,3,4,5, // 第2个面：4个点
                          4, 0,5,6,1, // 第3个面：4个点
                          4, 1,6,7,2, // 第4个面：4个点
                          4, 2,7,4,3};// 第5个面：4个点
ug->InsertNextCell(VTK_POLYHEDRON, 8, // 多面体共8个点
                  vtkIdList::New(), // 点列表（此处简化）
                  5, faceStream); // 5个面，面流信息

四、实战：创建你的第一个非结构化网格

下面用一段完整代码，创建一个包含“三角形（2D）”和“四面体（3D）”的非结构化网格，步骤清晰，新手可直接复制运行。

#include <vtkSmartPointer.h>
#include <vtkUnstructuredGrid.h>
#include <vtkPoints.h>
#include <vtkIdList.h>
#include <vtkCellTypes.h>

int main() {
    // 1. 实例化非结构化网格
    vtkSmartPointer<vtkUnstructuredGrid> ug = vtkSmartPointer<vtkUnstructuredGrid>::New();

    // 2. 定义点坐标（4个点：构成三角形和四面体）
    vtkSmartPointer<vtkPoints> points = vtkSmartPointer<vtkPoints>::New();
    points->InsertNextPoint(0.0, 0.0, 0.0); // 点0
    points->InsertNextPoint(1.0, 0.0, 0.0); // 点1
    points->InsertNextPoint(0.0, 1.0, 0.0); // 点2
    points->InsertNextPoint(0.0, 0.0, 1.0); // 点3
    ug->SetPoints(points); // 把点添加到网格

    // 3. 预分配内存（2个细胞，每个最多4个点）
    ug->AllocateEstimate(2, 4);

    // 4. 插入第一个细胞：三角形（VTK_TRIANGLE）
    vtkIdType triPoints[3] = {0, 1, 2}; // 三角形由点0、1、2组成
    ug->InsertNextCell(VTK_TRIANGLE, 3, triPoints);

    // 5. 插入第二个细胞：四面体（VTK_TETRA）
    vtkIdType tetraPoints[4] = {0, 1, 2, 3}; // 四面体由点0、1、2、3组成
    ug->InsertNextCell(VTK_TETRA, 4, tetraPoints);

    // 6. 建立拓扑链接，查询点0关联的细胞
    ug->BuildLinks();
    vtkSmartPointer<vtkIdList> cellIds = vtkSmartPointer<vtkIdList>::New();
    ug->GetPointCells(0, cellIds);

    // 7. 输出网格信息
    std::cout << "非结构化网格信息：" << std::endl;
    std::cout << "总点数：" << ug->GetNumberOfPoints() << std::endl; // 输出4
    std::cout << "总细胞数：" << ug->GetNumberOfCells() << std::endl; // 输出2
    std::cout << "点0关联的细胞ID：";
    for (int i = 0; i < cellIds->GetNumberOfIds(); i++) {
        std::cout << cellIds->GetId(i) << " "; // 输出0 1
    }
    std::cout << std::endl;

    // 8. 输出每个细胞的类型
    vtkSmartPointer<vtkCellTypes> cellTypes = vtkSmartPointer<vtkCellTypes>::New();
    ug->GetCellTypes(cellTypes);
    std::cout << "细胞类型：" << std::endl;
    for (int i = 0; i < ug->GetNumberOfCells(); i++) {
        std::cout << "细胞" << i << "：" << cellTypes->GetTypeName(ug->GetCellType(i)) << std::endl;
        // 输出：细胞0：TRIANGLE，细胞1：TETRA
    }

    return 0;
}

运行后，你会看到网格的点数、细胞数、细胞类型，以及点0关联的细胞——这就是一个最简单的非结构化网格！

五、新手避坑指南：3个常见错误

忘记调用BuildLinks()，导致拓扑查询失败
如果你用GetPointCells()或GetCellNeighbors()查询拓扑关系，却没先调用ug->BuildLinks()，会返回空结果或崩溃。因为拓扑链接不会自动建立，需要手动触发。
解决：查询前必须调用ug->BuildLinks()。
细胞类型与点数量不匹配
比如三角形（VTK_TRIANGLE）需要3个点，你却传了2个点；四面体（VTK_TETRA）需要4个点，传了3个点，会导致细胞创建失败。
解决：牢记常用细胞的点数量（三角形3、四面体4、六面体8），传参前核对。
没有预分配内存，效率低下
如果直接插入大量细胞（如10万个），不调用AllocateEstimate()或AllocateExact()，VTK会频繁扩容数组，速度很慢。
解决：根据预计的细胞数量和最大点数量，提前预分配内存。

六、适用场景：什么时候用vtkUnstructuredGrid？

记住3个核心场景，新手就能判断是否该选它：

复杂不规则几何体：如带凹槽的机械零件、有孔洞的建筑模型、分支复杂的血管；
有限元分析（FEA）结果：有限元软件（如ANSYS）输出的网格通常是四面体/六面体的非结构化网格，需要用它存储；
地质/气象模型：如起伏的地形、不均匀的大气网格，无法用规则体素表示。

不适用场景：

规则体数据（如CT/MRI）：用vtkImageData更高效；
简单离散面（如立方体表面）：用vtkPolyData更轻便。

七、总结

vtkUnstructuredGrid是VTK处理“非结构化复杂网格”的“万能工具”，它的核心优势是“灵活”——支持任意细胞类型组合，能适配各种不规则形状。但灵活性也意味着需要手动管理更多细节（如拓扑链接、内存分配），新手入门时先掌握“点→细胞→拓扑查询”的基础流程，再逐步深入多面体、大数据量优化等高级功能。

建议新手先跑通文中的实战代码，再尝试修改参数（如增加一个六面体细胞），慢慢熟悉它的特性。随着你处理的几何体越来越复杂，会发现vtkUnstructuredGrid是VTK中最强大的网格处理工具之一！