VTK 数据处理：特征边提取

最新推荐文章于 2025-03-18 20:26:45 发布
最新推荐文章于 2025-03-18 20:26:45 发布
阅读量1.2k 收藏 23
点赞数 19
分类专栏： VTK 文章标签： C++ VTK
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.youkuaiyun.com/ProgramNovice/article/details/139032583
版权

VTK 数据处理：特征边提取

VTK 数据处理：特征边提取

原理

VTK 的特征边提取只针对 PolyData，属于拓扑操作。

可以提取出 4 种类型的边：

  1. 边界边：只被 1 个网格使用的边；
  2. 流形边：被 2 个网格使用的边；
  3. 非流形边：被 2 个以上网格使用的边，一般不会在 PolyData 中出现；
  4. 特征边：属于流形边的一种，当一条流形边关联的两个面片的夹角大于特征角的话，它就是特征边。

实例 1：边界边提取

本实例创建了一个 vtkDiskSource 对象 disk，通过 vtkFeatureEdges 类提取 disk 的边界边，并设置不提取流形边、非流形边和特征边，并显示边界边的颜色以作区分。

VTKBoundaryEdge.h：

#pragma once

#include <QtWidgets/QMainWindow>
#include "ui_VTKFeatureEdge.h"

#include <QVTKOpenGLNativeWidget.h>

class VTKBoundaryEdge : public QMainWindow
{
	Q_OBJECT

public:
	VTKBoundaryEdge(QWidget* parent = nullptr);
	~VTKBoundaryEdge();

private:
	Ui::VTKFeatureEdgeClass ui;

	QVTKOpenGLNativeWidget* _pVTKWidget = nullptr;
};

VTKBoundaryEdge.cpp：

#include "VTKBoundaryEdge.h"

#include <vtkDiskSource.h>
#include <vtkPointData.h>
#include <vtkCellData.h>
#include <vtkPolyData.h>
#include <vtkFeatureEdges.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>

VTKBoundaryEdge::VTKBoundaryEdge(QWidget* parent) : QMainWindow(parent)
{
	ui.setupUi(this);

	_pVTKWidget = new QVTKOpenGLNativeWidget();
	this->setCentralWidget(_pVTKWidget);

	vtkNew<vtkRenderer> renderer;
	this->_pVTKWidget->renderWindow()->AddRenderer(renderer);
	this->_pVTKWidget->renderWindow()->Render();

	vtkNew<vtkDiskSource> disk;
	disk->Update();

	vtkNew<vtkFeatureEdges> featureEdges;
	featureEdges->SetInputConnection(disk->GetOutputPort());
	featureEdges->BoundaryEdgesOn(); // 提取边界边设置打开
	featureEdges->ManifoldEdgesOff();
	featureEdges->NonManifoldEdgesOff();
	featureEdges->FeatureEdgesOff();
	featureEdges->ColoringOn(); // 显示颜色

	vtkNew<vtkPolyDataMapper> diskMapper;
	diskMapper->SetInputConnection(disk->GetOutputPort());
	vtkNew<vtkPolyDataMapper> edgeMapper;
	edgeMapper->SetInputConnection(featureEdges->GetOutputPort());

	vtkNew<vtkActor> diskActor;
	diskActor->SetMapper(diskMapper);
	vtkNew<vtkActor> edgeActor;
	edgeActor->SetMapper(edgeMapper);

	renderer->AddActor(diskActor);
	renderer->AddActor(edgeActor);
}

VTKBoundaryEdge::~VTKBoundaryEdge()
{
}

运行结果：

在这里插入图片描述

实例 2：模型特征边提取

前面提到特征边的识别取决于两个面片的夹角，如下图所示，指针指着的角就是夹角，VTK 默认将该角大于 30 度的都视为特征边。可以看出，连续性不是很好的边都被视为特征边。

在这里插入图片描述

本实例使用 vtkBYUReader 读取了一个牛的模型，并用 vtkFeatureEdges 类提取了模型的特征边，并设置不提取流形边、非流形边和边界边，并显示特征边的颜色以作区分。

VTKFeatureEdge.h：

#pragma once

#include <QtWidgets/QMainWindow>
#include "ui_VTKFeatureEdge.h"

#include <QVTKOpenGLNativeWidget.h>

class VTKFeatureEdge : public QMainWindow
{
	Q_OBJECT

public:
	VTKFeatureEdge(QWidget* parent = nullptr);
	~VTKFeatureEdge();

private:
	Ui::VTKFeatureEdgeClass ui;

	QVTKOpenGLNativeWidget* _pVTKWidget = nullptr;
};

VTKFeatureEdge.cpp：

#include "VTKFeatureEdge.h"

#include <vtkBYUReader.h>
#include <vtkPointData.h>
#include <vtkCellData.h>
#include <vtkPolyData.h>
#include <vtkFeatureEdges.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>

VTKFeatureEdge::VTKFeatureEdge(QWidget* parent)
	: QMainWindow(parent)
{
	ui.setupUi(this);

	_pVTKWidget = new QVTKOpenGLNativeWidget();
	this->setCentralWidget(_pVTKWidget);

	vtkNew<vtkRenderer> renderer;
	this->_pVTKWidget->renderWindow()->AddRenderer(renderer);
	this->_pVTKWidget->renderWindow()->Render();

	vtkNew<vtkBYUReader> reader;
	reader->SetFileName("cow.g");
	reader->Update();

	vtkNew<vtkFeatureEdges> featureEdges;
	featureEdges->SetInputConnection(reader->GetOutputPort());
	featureEdges->BoundaryEdgesOff();
	featureEdges->ManifoldEdgesOn(); // 提取特征边设置打开
	featureEdges->NonManifoldEdgesOff();
	featureEdges->FeatureEdgesOff();
	featureEdges->ColoringOn(); // 显示颜色

	vtkNew<vtkPolyDataMapper> cowMapper;
	cowMapper->SetInputConnection(reader->GetOutputPort());
	vtkNew<vtkPolyDataMapper> edgeMapper;
	edgeMapper->SetInputConnection(featureEdges->GetOutputPort());

	vtkNew<vtkActor> cowActor;
	cowActor->SetMapper(cowMapper);
	vtkNew<vtkActor> edgeActor;
	edgeActor->SetMapper(edgeMapper);

	renderer->AddActor(cowActor);
	renderer->AddActor(edgeActor);
}

VTKFeatureEdge::~VTKFeatureEdge()
{}

实例 3：利用 vtkFeatureEdges 提取的边界补洞

建立一个 vtkPlane 类对象作为截面， 设置好 plane 的法向量。利用 vtkClipPolyData 类将模型作为输入，将 plane 设为 clipper 的截取函数。再用 vtkFeatureEdges 类提取截取后模型 clipper 的边界边。将边界边作为 vtkStripper 类对象 boundaryStrips 的输入，将 featureEdges 生成的三角片连接成三角带，再转换成多边形数据 boundaryPolyData。最后把 clipper、featureEdges、boundaryPolyData 全部显示出来。

#include "VTKFeatureEdge.h"

#include <vtkBYUReader.h>
#include <vtkPlane.h>
#include <vtkPointData.h>
#include <vtkCellData.h>
#include <vtkPolyData.h>
#include <vtkClipPolyData.h>
#include <vtkFeatureEdges.h>
#include <vtkStripper.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>

VTKFeatureEdge::VTKFeatureEdge(QWidget* parent)
	: QMainWindow(parent)
{
	ui.setupUi(this);

	_pVTKWidget = new QVTKOpenGLNativeWidget();
	this->setCentralWidget(_pVTKWidget);

	vtkNew<vtkRenderer> renderer;
	this->_pVTKWidget->renderWindow()->AddRenderer(renderer);
	this->_pVTKWidget->renderWindow()->Render();

	vtkNew<vtkBYUReader> reader;
	reader->SetFileName("cow.g");
	reader->Update();

	// 截面
	vtkNew<vtkPlane> plane;
	plane->SetOrigin(reader->GetOutput()->GetCenter());
	plane->SetNormal(1.0, -1.0, -1.0); // 设置截面的法向量

	// 截取模型
	vtkNew<vtkClipPolyData> clipper;
	clipper->SetInputData(reader->GetOutput());
	clipper->SetClipFunction(plane); // 将 plane 设为 clipper 的截取函数
	clipper->SetValue(0.0);

	vtkNew<vtkFeatureEdges> featureEdges;
	featureEdges->SetInputConnection(clipper->GetOutputPort());
	featureEdges->BoundaryEdgesOn();
	featureEdges->ManifoldEdgesOff(); // 提取特征边设置打开
	featureEdges->NonManifoldEdgesOff();
	featureEdges->FeatureEdgesOff();
	featureEdges->ColoringOn(); // 显示颜色

	// 建立三角带对象
	vtkNew<vtkStripper> boundaryStrips;
	boundaryStrips->SetInputConnection(featureEdges->GetOutputPort()); // 将 featureEdges 生成的三角片连接成三角带
	boundaryStrips->Update();

	vtkNew<vtkPolyData> boundaryPolyData;
	boundaryPolyData->SetPoints(boundaryStrips->GetOutput()->GetPoints());
	boundaryPolyData->SetPolys(boundaryStrips->GetOutput()->GetLines());

	vtkNew<vtkPolyDataMapper> clipperMapper;
	clipperMapper->SetInputConnection(clipper->GetOutputPort());
	vtkNew<vtkPolyDataMapper> edgeMapper;
	edgeMapper->SetInputConnection(featureEdges->GetOutputPort());
	vtkNew<vtkPolyDataMapper> boundaryMapper;
	boundaryMapper->SetInputData(boundaryPolyData);

	vtkNew<vtkActor> clipperActor;
	clipperActor->SetMapper(clipperMapper);
	vtkNew<vtkActor> edgeActor;
	edgeActor->SetMapper(edgeMapper);
	vtkNew<vtkActor> boundaryActor;
	boundaryActor->SetMapper(boundaryMapper);

	renderer->AddActor(clipperActor);
	renderer->AddActor(edgeActor);
	renderer->AddActor(boundaryActor);
}

VTKFeatureEdge::~VTKFeatureEdge()
{}

运行结果：

在这里插入图片描述

实例 4：利用 vtkFillHolesFilter 补洞

vtkFillHolesFilter 是 VTK 中用于自动填充三角网格模型中洞的过滤器。vtkFillHolesFilter 支持多种数据格式，包括 vtkPolyData、vtkUnstructuredGrid、vtkStructuredGrid、vtkImageData 等。它还支持多种洞填充算法，包括 Delaunay 三角剖分、路径填充、平面填充等。使用 vtkFillHolesFilter 进行填洞时，可以保持原始模型的拓扑结构，避免生成不良的三角形。

vtkFillHolesFilter 的内部执行过程是首先检测出网格中的所有边界边，然后找出这些边界边中的每一个闭合回路，最后将这些闭合回路进行三角化（即生成三角网格）以实现填补的目的。这个类也是非常简单的，只需要设置需要填补的网格数据即可。

需要注意的是，有些边界的闭合回路是不需要三角化的，例如一个平面网格，若填补其四周的边界边，则会与原网格产生覆盖。vtkFillHolesFilters 中的 SetHoleSize() 函数可用于控制需要修补的漏洞面积的最大值，大于该值的漏洞则不需要填补处理。

#include "VTKFeatureEdge.h"

#include <vtkBYUReader.h>
#include <vtkPlane.h>
#include <vtkPointData.h>
#include <vtkCellData.h>
#include <vtkPolyData.h>
#include <vtkClipPolyData.h>
#include <vtkFeatureEdges.h>
#include <vtkStripper.h>
#include <vtkFillHolesFilter.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>

VTKFeatureEdge::VTKFeatureEdge(QWidget* parent)
	: QMainWindow(parent)
{
	ui.setupUi(this);

	_pVTKWidget = new QVTKOpenGLNativeWidget();
	this->setCentralWidget(_pVTKWidget);

	vtkNew<vtkRenderer> renderer;
	this->_pVTKWidget->renderWindow()->AddRenderer(renderer);
	this->_pVTKWidget->renderWindow()->Render();

	vtkNew<vtkBYUReader> reader;
	reader->SetFileName("cow.g");
	reader->Update();

	// 截面
	vtkNew<vtkPlane> plane;
	plane->SetOrigin(reader->GetOutput()->GetCenter());
	plane->SetNormal(1.0, -1.0, -1.0); // 设置截面的法向量

	// 截取模型
	vtkNew<vtkClipPolyData> clipper;
	clipper->SetInputData(reader->GetOutput());
	clipper->SetClipFunction(plane); // 将 plane 设为 clipper 的截取函数
	clipper->SetValue(0.0);

	/*
	vtkNew<vtkFeatureEdges> featureEdges;
	featureEdges->SetInputConnection(clipper->GetOutputPort());
	featureEdges->BoundaryEdgesOn();
	featureEdges->ManifoldEdgesOff(); // 提取特征边设置打开
	featureEdges->NonManifoldEdgesOff();
	featureEdges->FeatureEdgesOff();
	featureEdges->ColoringOn(); // 显示颜色

	// 建立三角带对象
	vtkNew<vtkStripper> boundaryStrips;
	boundaryStrips->SetInputConnection(featureEdges->GetOutputPort()); // 将 featureEdges 生成的三角片连接成三角带
	boundaryStrips->Update();

	vtkNew<vtkPolyData> boundaryPolyData;
	boundaryPolyData->SetPoints(boundaryStrips->GetOutput()->GetPoints());
	boundaryPolyData->SetPolys(boundaryStrips->GetOutput()->GetLines());

	vtkNew<vtkPolyDataMapper> clipperMapper;
	clipperMapper->SetInputConnection(clipper->GetOutputPort());
	vtkNew<vtkPolyDataMapper> edgeMapper;
	edgeMapper->SetInputConnection(featureEdges->GetOutputPort());
	vtkNew<vtkPolyDataMapper> boundaryMapper;
	boundaryMapper->SetInputData(boundaryPolyData);

	vtkNew<vtkActor> clipperActor;
	clipperActor->SetMapper(clipperMapper);
	vtkNew<vtkActor> edgeActor;
	edgeActor->SetMapper(edgeMapper);
	vtkNew<vtkActor> boundaryActor;
	boundaryActor->SetMapper(boundaryMapper);

	renderer->AddActor(clipperActor);
	renderer->AddActor(edgeActor);
	renderer->AddActor(boundaryActor);
	*/

	vtkNew<vtkFillHolesFilter> fillHolesFilter;
	fillHolesFilter->SetInputConnection(clipper->GetOutputPort());
	fillHolesFilter->SetHoleSize(1000.0); // 大于该值的漏洞则不需要填补处理

	vtkNew<vtkPolyDataMapper> fillHolesMapper;
	fillHolesMapper->SetInputConnection(fillHolesFilter->GetOutputPort());

	vtkNew<vtkActor> fillHolesActor;
	fillHolesActor->SetMapper(fillHolesMapper);

	renderer->AddActor(fillHolesActor);
}

VTKFeatureEdge::~VTKFeatureEdge()
{}

运行结果：

在这里插入图片描述

VTK：边缘提取的使用方法及编程实现
CyberBladeX的博客
09-13 369
在VTK（Visualization Toolkit）中，我们可以使用各种滤波器和算法来提取图像或模型的边缘。本文将介绍VTK中边缘提取的使用方法，并提供相应的源代码示例。在进行边缘提取之前，我们需要加载要处理的图像或模型数据。通过以上步骤，我们可以在VTK中实现边缘提取的功能。根据需要选择合适的滤波器和参数，可以实现不同类型的边缘提取效果。最后，我们需要创建一个渲染器，并将边缘结果添加到渲染器中进行显示。要使用VTK进行边缘提取，首先需要导入必要的库和模块。，并将边缘提取结果添加到渲染器中的。
VTK笔记——提取边界(vtkFeatureEdges)
Andy的专栏
04-04 7224
vtkFeatureEdges，VTK的一个类，主要用于边缘提取。 提取的边缘包括： 1)边界(boundary)，只被一个多边形使用的边，或者直线单元； 2)非流行(non-manifold)，被三个以上的多边形共用的边； 3)特征边(feature edges)，被两个三角形使用且二面角大于特征角的边； 4)流行边(manifold edges)，只被两个多边形使用的边。
五、vtkFeatureEdges边过滤器
最新发布
weixin_51795597的博客
03-18 462
1、vtkFeatureEdges是VTK中的一个过滤器类，用于从输入的多边形数据（vtkPolyData）中提取特征边（Feature Edge）。特征边通常是指几何模型中的边界边、非流行边、锐边或轮廓边等。2、vtkFeatureEdges是一个非常有用的工具，特别适合用于几何处理、模型分析和可视化任务中。通过提取特征边，可以更好地理解和处理复制地几何模型。
vtk Edges 特征边 提取 网格封闭性检测
q610098308的专栏
10-20 1439
vtkFeatureEdges，VTK的一个类，用于特征边界提取 提取的边界包括： 1)边界(boundary)，只被一个多边形使用的边，或者直线单元； 2)非流行(non-manifold)，被三个以上的多边形共用的边； 3)特征边(feature edges)，被两个三角形使用且二面角大于特征角的边； 4)流行边(manifold edges)，只被两个多边形使用的边。
vtkFeatureEdges 类简介
Jane_yuhui的专栏
08-31 1885
vtkFeatureEdges 【类说明】 提取多边形数据（polygonal data）的边界，非流行、锋利的边缘。 该类是针对输入数据为多边形几何数据的滤波器，提取特定类型的边缘数据。 这些边缘可以是： 一个多边形的边界 boundary 或者一个 line cell;非流形数据 - 三个或者更多的多边形数据；特征边（两个三角形的边缘， 三角形的二面角大于特征角）流形边（
VTK笔记——提取边界(vtkFeatureEdges)和检测缺面
weixin_45793544的博客
12-06 160
看到一篇文章：https://blog.youkuaiyun.com/minmindianzi/article/details/88951927，做下记录vtkFeatureEdges，VTK的一个类，用于边界提取。提取的边界包括：1)边界(boundary)，只被一个多边形使用的边，或者直线单元；2)非流行(non-manifold)，被三个以上的多边形共用的边；3)特征边(feature edges)，被两个三角形使用且二面角大于特征角的边；
vtkFeatureEdges提取多边形数据“指定”的轮廓
Jane_yuhui的专栏
01-18 2641
该类的主要目的就是为了从多边形数据中提取边界、非流形边、锐利边缘/轮廓。 该类是一个过滤器，为输入的多边形数据提取特殊类型/指定的边缘，边界类型可以是 ① 分界线（一个多边形的）或者是一条边； ② 非流形的 （被三个或者更多的多边形使用）； ③ 特征边（被两个三角形使用的边界，且二面角大于特征角）； ④ 流形边（只是被两个多边形使用的边） 提取的边界可能是任意连接的。
VTK寻找网格边界
m0_67357141的博客
07-25 1222
3.特征边(featureedges)被两个三角形使用，并且两个三角形之间的夹角大于特征角度（可自行设置此角度）；1.边界(boundary)只被一个多边形使用的边，或者直线单元；2.非流行(non-manifold)被三个以上的多边形共用的边；4.流行边(manifoldedges)只被两个多边形使用的边。.........
vtk-examples:python vtk示例集合
05-18
- **数据处理**：通过VTK滤波器对数据进行操作，如体积分割、表面提取等。 - **渲染与显示**：如何设置光源、颜色和透明度，以及创建窗口和相机视图来展示3D模型。 - **交互控制**：实现用户与3D模型的交互，如鼠标...
VTK图像处理：子采样与vtkExtractVOI滤波器
系统设计部分，VTK遵循模块化和层次化的架构，包含多个层次的组件，如数据模型、数据处理、渲染等，每个组件都有特定的职责，协同工作以实现高效的可视化。通过理解这些组件的工作原理和交互方式，开发者可以更好地...
VTK点云数据读取：点云增强与特征提取，专业技能提升
本文全面系统地介绍了VTK（Visualization Toolkit）在点云数据处理中的应用，涵盖了从数据读取、预处理到增强技术、特征提取及高级处理技术的各个方面。通过对点云增强、特征提取、三维重建和压缩存储等技术的详细...
基于VTK的医学三维图像模型构建与切割
07-12
：虚拟外科手术是医学领域的一个重要研究方向。基于面绘制和体绘制两种 可视化算法实现医学图像的三维模型重建，针对不同重建模型分别实现虚拟任意切割：面切 割和体切割，并讨论了虚拟切割的交互操作实现方法。用VTK(Visualization Toolkit)if)7步实 现了模拟手术系统中的虚拟切割脊柱体功能。
深入理解VTK数据结构：3D渲染性能优化技巧
![深入理解VTK数据结构：3D渲染性能优化技巧]... VTK数据结构基础 VTK（Visualization Toolkit）是一种开源的软件系统，用于3D计算机图形学、图像处理和可视化领
VTK网格处理：从入门到高级建模的全面指南
本文全面介绍了VTK（Visualization Toolkit）在数据表示、网格处理、可视化、分析以及建模实践中的应用。首先，介绍了VTK的基本概念、安装配置方法，进而深入探讨了网格处理的基础，包括数据的表示方法、基本操作和...
VTK轮廓提取
OneSea的专栏
06-06 4276
参考资料 [1]VTK教程之九 可视化基础算法-二维轮廓线（等值线）提取
VTK中对模型进行封闭性检测
黄岛主的博客
09-16 2566
由于受原始数据、重建方法的限制，得到的网格模型并不是封闭的。有时为了显示或者处理某些要求，需要网格必须是封闭的。比如一个球面网格就是封闭性网格。 如果一条边只被一个多边形包含，那么这条边就是边界边。是否存在边界边是检测一个网格模型是否封闭的重要特征。vtkFeatureEdges是一个非常重要的类，该类能够提取多边形网格模型中四种类型的边： 边界边：即只被一个多边形或者一条边包含的边。 非...
提取物体轮廓实例：通过VTK库实现
DevProZ的博客
08-23 586
提取物体轮廓实例：通过VTK库实现本文将介绍如何通过使用VTK库，来实现物体轮廓的提取。在本例中，我们将使用VTK的Marching Cubes算法来提取物体表面的网格，在此基础上计算出物体的轮廓边缘。步骤一：导入VTK库并读入数据首先，我们需要导入VTK库，并读入待处理的数据。在本例中，我们将使用示例数据文件“cow.vtp”，该文件可以在VTK官方网站上下载得到。请确保已经正确安装VTK库并且Python环境已经配置好。import vtk。
VTK 获取物体外轮廓 ExtractOutsideSurface
q610098308的专栏
02-19 1254
有时我们需要从多个表面体中，提取 最后层哪个，简单的方法 可以使用这种：
VTK随笔十：VTK图形处理（封闭性检测、联通区域分析、多分辨率处理）
u012959478的博客
08-30 1509
VTK随笔十：VTK图形处理（封闭性检测、联通区域分析、多分辨率处理）
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

UestcXiye

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值