【VTK手册037】vtkFeatureEdges：多边形网格特征提取利器
0. 概述
在医学图像处理与三维建模中，识别网格的边界、孔洞、非流形结构以及尖锐棱角是拓扑分析与可视化增强的关键环节。vtkFeatureEdges是 VTK 中专门用于从多边形数据（vtkPolyData）中提取特定类型拓扑边的核心滤镜。它能够根据邻域拓扑关系和几何夹角，精确分离出边界边、特征边、流形边及非流形边。
1. 快速示例
以下示例展示了如何提取一个网格的所有边界边（Boundary Edges）和特征边（Feature Edges）：
#include<vtkFeatureEdges.h>#include<vtkSmartPointer.h>#include<vtkPolyData.h>#include<vtkPolyDataMapper.h>// 1. 初始化滤镜autofeatureEdges=vtkSmartPointer<vtkFeatureEdges>::New();featureEdges->SetInputData(inputPolyData);// 2. 配置提取参数featureEdges->BoundaryEdgesOn();// 提取边界边（如孔洞边缘）featureEdges->FeatureEdgesOn();// 提取尖锐棱角featureEdges->SetFeatureAngle(30.0);// 设置特征角阈值为 30 度featureEdges->NonManifoldEdgesOff();// 关闭非流形边提取featureEdges->ManifoldEdgesOff();// 关闭普通流形边提取// 3. 配置辅助功能featureEdges->ColoringOn();// 按边类型给输出赋予标量值featureEdges->Update();// 4. 获取结果vtkPolyData*result=featureEdges->GetOutput();2. 基本原理与核心公式
vtkFeatureEdges的核心逻辑基于边与单元（Cell）的邻接关系以及二面角（Dihedral Angle）的计算。
2.1 边类型定义
该滤镜通过统计每条边被多少个多边形（Polygon）共享来分类：
- 边界边 (Boundary Edges)：仅被 1 个多边形共享的边，或属于线单元（Line cell）的边。
- 流形边 (Manifold Edges)：恰好被 2 个多边形共享的边。
- 非流形边 (Non-Manifold Edges)：被 3 个或更多个多边形共享的边。
- 特征边 (Feature Edges)：属于流形边，但其相邻两个多边形的法向量夹角（二面角）超过了指定的阈值。
2.2 特征角计算公式
对于共享同一条边的两个三角形T1T_1T1和T2T_2T2，设其单位法向量分别为n⃗1\vec{n}_1n1和n⃗2\vec{n}_2n2。它们之间的二面角θ\thetaθ定义为：
θ=arccos(n⃗1⋅n⃗2)\theta = \arccos(\vec{n}_1 \cdot \vec{n}_2)θ=arccos(n1⋅n2)
当θ>FeatureAngle\theta > \text{FeatureAngle}θ>FeatureAngle时，该边被判定为特征边。
3. 源码逻辑深度分析
vtkFeatureEdges的内部实现主要在RequestData中完成：
- 拓扑构建：算法首先调用
BuildLinks()建立点到单元的索引。 - 边遍历：遍历输入网格的每一个多边形单元。对于每一条边，查找其共用单元的数量NNN。
- 分类逻辑：
- 若N=1N=1N=1，标记为
Boundary。 - 若N>2N>2N>2，标记为
NonManifold。 - 若N=2N=2N=2，计算相邻面法线夹角。若夹角>FeatureAngle> \text{FeatureAngle}>FeatureAngle，标记为
Feature；否则标记为Manifold。
- 若N=1N=1N=1，标记为
- 坐标合并：利用
vtkIncrementalPointLocator（默认为vtkMergePoints）来合并提取出的边段中的重复点，保证拓扑连续性。 - 着色机制：如果开启
Coloring，算法会为提取出的每一条边分配CellData标量值，不同类型的边对应不同的数值。
4. 关键接口列表
以下是开发中常用的核心 API：
4.1 提取开关控制
| 函数接口 | 说明 |
|---|---|
void ExtractAllEdgeTypesOn() | 开启所有类型的边提取。 |
void ExtractAllEdgeTypesOff() | 关闭所有类型的边提取。 |
void SetBoundaryEdges(bool) | 设置是否提取边界边（由 1 个多边形使用的边）。 |
void SetFeatureEdges(bool) | 设置是否提取特征边（二面角 > FeatureAngle）。 |
void SetNonManifoldEdges(bool) | 设置是否提取非流形边（由 3 个以上多边形使用的边）。 |
void SetManifoldEdges(bool) | 设置是否提取流形边（由 2 个多边形使用的边）。 |
4.2 参数配置
| 函数接口 | 说明 |
|---|---|
void SetFeatureAngle(double) | 设置特征角阈值，取值范围[0,180][0, 180][0,180]，默认 30°。 |
void SetPassLines(bool) | 是否直接透传输入数据中的线单元（Lines）。 |
void SetColoring(bool) | 是否根据边类型对输出结果进行标量着色。 |
void SetRemoveGhostInterfaces(bool) | 是否移除 Ghost 单元界面上的边（并行计算常用）。 |
void SetOutputPointsPrecision(int) | 设置输出点数据的精度（Single/Double）。 |
4.3 基础设施
| 函数接口 | 说明 |
|---|---|
void SetLocator(vtkIncrementalPointLocator*) | 设置用于合并点的定位器。 |
void CreateDefaultLocator() | 创建默认的vtkMergePoints定位器。 |
5. 开发建议与注意事项
- 数据预处理：
vtkFeatureEdges对输入网格的拓扑质量较为敏感。如果输入网格存在微小的几何重叠但未焊接，会导致错误的边界识别。建议先执行vtkCleanPolyData。 - 可视化提示：开启
ColoringOn()后，若要观察到不同颜色，需设置 Mapper 的标量模式：mapper->SetScalarModeToUseCellData()。 - 性能考量：对于超大规模网格，
FeatureEdges的计算量较大，特别是计算二面角的过程。如非必要，可关闭ManifoldEdges以减少结果集大小。
