VTK 数据处理：特征边提取

2024-05-26 16:04
文章标签 提取 特征 数据处理 vtk

本文主要是介绍VTK 数据处理：特征边提取，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

VTK 数据处理：特征边提取

  • VTK 数据处理：特征边提取
    • 原理
    • 实例 1：边界边提取
    • 实例 2：模型特征边提取
    • 实例 3：利用 vtkFeatureEdges 提取的边界补洞
    • 实例 4：利用 vtkFillHolesFilter 补洞

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_OBJECTpublic: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_OBJECTpublic: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 数据处理：特征边提取的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/1004930。

相关文章

Python xmltodict实现简化XML数据处理

Python xmltodict实现简化XML数据处理

《Pythonxmltodict实现简化XML数据处理》Python社区为提供了xmltodict库,它专为简化XML与Python数据结构的转换而设计,本文主要来为大家介绍一下如何使用xmltod... 目录一、引言二、XMLtodict介绍设计理念适用场景三、功能参数与属性1、parse函数2、unpa
阅读更多...
Java后端接口中提取请求头中的Cookie和Token的方法

Java后端接口中提取请求头中的Cookie和Token的方法

《Java后端接口中提取请求头中的Cookie和Token的方法》在现代Web开发中,HTTP请求头（Header）是客户端与服务器之间传递信息的重要方式之一,本文将详细介绍如何在Java后端（以Sp... 目录引言1. 背景1.1 什么是 HTTP 请求头？1.2 为什么需要提取请求头？2. 使用 Spr
阅读更多...
Python数据处理之导入导出Excel数据方式

Python数据处理之导入导出Excel数据方式

《Python数据处理之导入导出Excel数据方式》Python是Excel数据处理的绝佳工具,通过Pandas和Openpyxl等库可以实现数据的导入、导出和自动化处理,从基础的数据读取和清洗到复杂... 目录python导入导出Excel数据开启数据之旅：为什么Python是Excel数据处理的最佳拍档
阅读更多...
使用Java解析JSON数据并提取特定字段的实现步骤(以提取mailNo为例)

使用Java解析JSON数据并提取特定字段的实现步骤(以提取mailNo为例)

《使用Java解析JSON数据并提取特定字段的实现步骤(以提取mailNo为例)》在现代软件开发中,处理JSON数据是一项非常常见的任务,无论是从API接口获取数据,还是将数据存储为JSON格式,解析... 目录1. 背景介绍1.1 jsON简介1.2 实际案例2. 准备工作2.1 环境搭建2.1.1 添加
阅读更多...
python解析HTML并提取span标签中的文本

python解析HTML并提取span标签中的文本

《python解析HTML并提取span标签中的文本》在网页开发和数据抓取过程中,我们经常需要从HTML页面中提取信息,尤其是span元素中的文本,span标签是一个行内元素,通常用于包装一小段文本或... 目录一、安装相关依赖二、html 页面结构三、使用 BeautifulSoup javascript
阅读更多...
OmniGlue论文详解（特征匹配）

OmniGlue论文详解（特征匹配）

OmniGlue论文详解（特征匹配） 摘要1. 引言2. 相关工作2.1. 广义局部特征匹配2.2. 稀疏可学习匹配2.3. 半稠密可学习匹配2.4. 与其他图像表示匹配 3. OmniGlue3.1. 模型概述3.2. OmniGlue 细节3.2.1. 特征提取3.2.2. 利用DINOv2构建图形。3.2.3. 信息传播与新的指导3.2.4. 匹配层和损失函数3.2.5. 与Super
阅读更多...
Python：豆瓣电影商业数据分析-爬取全数据【附带爬虫豆瓣，数据处理过程，数据分析，可视化，以及完整PPT报告】

Python：豆瓣电影商业数据分析-爬取全数据【附带爬虫豆瓣，数据处理过程，数据分析，可视化，以及完整PPT报告】

**爬取豆瓣电影信息，分析近年电影行业的发展情况** 本文是完整的数据分析展现，代码有完整版，包含豆瓣电影爬取的具体方式【附带爬虫豆瓣，数据处理过程，数据分析，可视化，以及完整PPT报告】   最近MBA在学习《商业数据分析》，大实训作业给了数据要进行数据分析，所以先拿豆瓣电影练练手，网络上爬取豆瓣电影TOP250较多，但对于豆瓣电影全数据的爬取教程很少，所以我自己做一版。 目
阅读更多...
ROS - C++实现RosBag包回放/提取

ROS - C++实现RosBag包回放/提取

文章目录 1. 回放原理2. 回放/提取 多个话题3. 回放/提取数据包，并实时发布 1. 回放原理 #include <ros/ros.h>#include <rosbag/bag.h>#include <std_msgs/String.h>int main(int argc, char** argv){// 初始化ROS节点ros::init(argc, argv,
阅读更多...
《计算机视觉工程师养成计划》 ·数字图像处理·数字图像处理特征·概述~

《计算机视觉工程师养成计划》 ·数字图像处理·数字图像处理特征·概述~

1 定义         从哲学角度看：特征是从事物当中抽象出来用于区别其他类别事物的属性集合，图像特征则是从图像中抽取出来用于区别其他类别图像的属性集合。         从获取方式看：图像特征是通过对图像进行测量或借助算法计算得到的一组表达特性集合的向量。 2 认识         有些特征是视觉直观感受到的自然特征，例如亮度、边缘轮廓、纹理、色彩等。         有些特征需要通
阅读更多...
初步了解VTK装配体

初步了解VTK装配体

VTK还不太了解，根据资料， vtk.vtkAssembly 是 VTK库中的一个重要类，允许通过将多个vtkActor对象组合在一起来创建复杂的3D模型。 import vtkimport mathfrom vtk.util.colors import *filenames = ["cylinder.stl","sphere.stl","torus.stl"]dt = 1.0renW
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2