VTK：使用C/C++填充孔的方法

最新推荐文章于 2023-09-18 01:18:16 发布
最新推荐文章于 2023-09-18 01:18:16 发布
阅读量379 收藏 1
点赞数 1
CC 4.0 BY-SA版权
文章标签： c语言 c++ java C/C++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.youkuaiyun.com/PixelProX/article/details/132659643
C/C++ 专栏收录该内容
230 篇文章 ¥59.90 ¥99.00
订阅专栏
本文介绍了如何利用C/C++和VTK库来填充三维模型中的孔洞。通过加载STL格式的模型文件，使用vtkPolyDataMapper、vtkPolyData和vtkFillHolesFilter类，可以有效地处理模型的不连续性，生成填充后的模型文件。

VTK：使用C/C++填充孔的方法

VTK（Visualization Toolkit）是一个功能强大的开源软件库，用于可视化和图形处理。在VTK中，填充孔是指在三维模型中填充由于几何形状之间的不连续性而产生的空洞。本文将介绍如何使用C/C++语言和VTK库来填充这些孔洞。

首先，我们需要准备一个三维模型，该模型包含一个或多个需要填充的孔洞。VTK支持多种三维模型格式，如STL、OBJ等。在本文中，我们将使用STL格式的模型文件。

以下是一个示例STL模型文件（model.stl）的内容：

solid model
  facet normal 0.0 0.0 0.0
    outer loop
      vertex 0.0 0.0 0.0
      vertex 1.0 0.0 0.0
      vertex 0.0 1.0 0.0
    endloop
  endfacet
  facet normal 0.0 0.0 0.0
    outer loop
      vertex 1.0 0.0 0.0
      vertex 1.0 1.0 0.0
      vertex 0.0 1.0 0.0
    endloop
  endfacet
endsolid model

上述示例模型文件中包含两个相邻的三角形，它们之间形成了一个需要填充的孔洞。

接下来，我们将编写C/C++代码，使用VTK库加载模型文件并进行孔洞填充。我们将使用VTK的vtkPolyDataMapper、vtkPolyData和vtkFillHolesFilter类来完成这个任务。

下面是填充孔洞的C/C++代码：

了解本专栏 订阅专栏 解锁全文
实战：使用VTK进行孔洞填充
PixelLoom的博客
08-22 703
首先，我们需要了解VTK中提供的孔洞填充算法。VTK中提供了多种孔洞填充算法，包括基于高斯曲率和球形映射的方法等。本文将采用基于高斯曲率的孔洞填充算法。该算法通过计算每个点的高斯曲率值，判断其是否为孔洞，并将其填充。最后，调用了vtkFillHolesFilter类中提供的函数来填充孔洞，并将结果写入到PLY格式文件中。这些孔洞会影响我们对模型的后续处理，所以需要进行孔洞填充。经过上述步骤，我们就可以成功地填充了三维模型中的孔洞。通过使用VTK提供的孔洞填充算法，我们可以轻松地对三维模型进行后续处理。
VTK-空洞填充（vtkFillHolesFilter）
qq_40041064的博客
04-04 1274
本文主要介绍vtk中空洞填充接口vtkFillHolesFilter实现原理，存在的Bug及修复方法，希望和小伙伴一起进步。
VTK应用CGAL布尔运算、三角剖分、补洞、模型偏置（膨胀）降采样、骨架提取
06-29
只要你的工程有VTK、CGAL5.5环境就可以使用。
利用vtk+cgal+openmesh（或者第三方格式转换软件）做牙齿模型
qq_23158477的博客
04-29 1757
如何建立牙齿模型？使用vtk+cgal+openmesh实现模型建立与修复优化vtk实现dicom图片生成模型，运用阈值图像法生成各类模型，不同部位牙齿和骨头的分离阈值稍有变化，运用滑动条实现阈值的控制。参考frog 大脑的提取。CGAL实现vtk中分离过程中生成的洞的修复，vtk在补洞上效果太差，所以选择CGAL的算法，实践方便效果明显。openmesh实现off与stl（obj）文件格式的转换，比较方便，但可能我自己的问题，在代码实现转换时模型发生了变化，我采用3D convert软件来实现模型转换。因
C++ Opencv imfill 孔洞填充函数
03-25
C++ Opencv imfill 孔洞填充函数
VTK：绘制轮廓的C/C++用法
TechProX的博客
09-05 281
在VTK中，轮廓绘制是一种常见的技术，用于显示3D对象的边界或特定数值范围的等值面。在VTK中，轮廓绘制是一种常见的技术，用于显示3D对象的边界或特定数值范围的等值面。你可以根据自己的需求修改示例数据和等值面的数值，以实现不同的轮廓效果。在开始之前，请确保你已经正确安装了VTK库，并且设置好了编译环境。你可以根据自己的需求修改示例数据和等值面的数值，以实现不同的轮廓效果。接着，我们创建了一个渲染器和窗口，并将演员添加到渲染器中。接着，我们创建了一个渲染器和窗口，并将演员添加到渲染器中。
VTK：填充孔用法实战
希望我的博客，能帮上你解决学习中工作中所遇到的问题
02-06 643
VTK：填充孔用法实战程序输出程序完整源代码 程序输出 程序完整源代码 #include <vtkActor.h> #include <vtkCamera.h> #include <vtkCleanPolyData.h> #include <vtkDataSetSurfaceFilter.h> #include <vtkExtractSelection.h> #include <vtkFillHolesFilter.h> #inclu
VTK：使用裁剪操作在多边形数据上的演示（C/C++）
TechWhiz的博客
09-18 264
在本文中，我们将介绍如何使用 Visualization Toolkit（VTK）库中的裁剪操作来处理多边形数据。我们将展示如何使用C/C++语言编写代码来执行裁剪操作，并提供相应的源代码示例。它提供了广泛的功能，可以处理各种类型的数据，包括多边形数据。通过上述代码，你可以在VTK中使用裁剪操作处理多边形数据。你可以根据自己的需求修改代码中的多边形数据和裁剪平面，以适应你的具体应用场景。最后，我们将渲染器添加到渲染窗口中，并设置窗口的大小。在上面的代码中，我们首先创建了一个多边形数据。方法来显示渲染结果。
VTK：绘制圆圈用法实战
希望我的博客，能帮上你解决学习中工作中所遇到的问题
02-01 1150
VTK：绘制圆圈用法实战程序输出程序完整源代码 程序输出 程序完整源代码 // Classes specific to this example #include <vtkPolyData.h> #include <vtkPolyDataMapper.h> #include <vtkRegularPolygonSource.h> // Generic VTK pipeline elements #include <vtkActor.h> #include &
vtk普通补洞过程实现
qq_23158477的博客
08-07 1101
思路分析：模型表面漏洞，肯定有边界边，可能不止一个洞。1.提取边界边;2.运用联通性，依次选择边界环;3.边界环连接成面（1.取所有点，取中心点坐标（xi / n,yi / n,zi / n）;2.线单元的两个点与中心点形成面单元；3.将点与面单元赋值数据形成新的补洞面数据）；4.面与之前的模型结合，clean后形成最终补洞模型。 vtkPolyData* CommonMesh::CloseSurface(vtkPolyData* shape) { vtkSmartPointe...
C++ 孔洞填充代码
08-26
冈萨雷斯的数字图像处理孔洞填充方法运行缓慢，写了个简略的版本。文档附带测试图片
强大的补洞程序
06-05
这是一个很强大的补洞程序，可以用来修复图像上的一些洞，更进一步的进行一些光滑性处理
3D三角形网格模型补洞源代码
06-28
3D三角形网格模型补洞算法源代码，原为linux平台下编译，现已改为Windows环境下；由于使用了C++11新特性，需通过VS2013或更高版本编译。
vtk填洞算法
qq_23158477的博客
04-12 879
#include “SurfaceHoleFiller.h” #include “UmbrellaWeightedOrder2Smoother.h” #include “CoverRefiner.h” #include <vtkCell.h> #include <vtkPolyDataNormals.h> #include <vtkPointData.h> #include <vtkDataSetAttributes.h> #include <vtkCe
使用OpenCV（C ++）填充图像中的孔
进击的Explorer
11-13 1555
在本教程中，我们将学习如何填充二进制图像中的孔。考虑一下图1左侧的图像。假设我们想找到一个二进制蒙版，它将硬币与背景分离开，如右图所示。在本教程中，包含硬币的圆形区域也称为前景。 请注意，硬币的边界很暗，与白色背景不同。因此，我们使用简单的图像阈值处理将边界与背景分开。换句话说，我们说强度高于特定值（阈值）的像素是背景，其余像素是前景。上示中间图像显示此阈值图像（黑色代表背景，白色代表...
C++ VTK网格模型补洞填孔
alicema1111的博客
05-20 1010
这篇博客针对C++ VTK网格模型补洞填孔编写代码，代码整洁，规则，易读。 学习与应用推荐首选。
灰度图像孔洞填补C++实现
qq_38691400的博客
11-29 2092
由于项目需要，需要对灰度图进行孔洞填补，在网上查找了一圈发现基本都是将灰度图二值化在通过漫水填充的方法来实现的（imfill）。但是由于我的灰度图中包含有重要信息，二值化会导致关键信息的丢失。没办法只有自己想了个笨办法解决。（一个数学、逻辑、编程都菜的扣jio的菜狗能想到办法并将之实现是件很不容易的事所以特此记录） - 根据所得到的灰度图像特点，检测图像中每个像素点的灰度值，如果检测结果为“当前像素灰度值为0，而下一个像素值不为0”则将该点坐标存入DeleteHeadMinuend； - 如果检测结果为
C++ Opencv imfill 孔洞填充函数的实现(更新后)
weixin_44841157的博客
03-25 6275
目录一、须知二、演示过程代码展示主函数展示原图运行结果三、总结 一、须知 本文章所提供代码不是自创，由于时间太久实在找不到来源，发布出来只为给大家提供便利，完全免费。 话不多说，不想看文章的直接点击下载链接即可：点我. 二、演示过程 代码展示 Mat cop二值图 int n 填充比n小的孔洞 函数默认为4连通 如想改为8连通 自行修改代码即可。 #include"imfill.h" Mat imfill(Mat cop,int n) { Mat data = ~cop; Mat labels,
C++ OpenCV——孔洞填充（添加优化代码解决四周伪孔洞填充问题）
jiDxiaohuo的博客
09-28 1291
工程代码下载传送门： 普遍的孔洞填充结果 孔洞定义为：被前景元素包围的后景的元素 虽然从定义上来看四周不被完全包围的元素不应该为孔洞，但是我还是想改善一下使周围这些（暂且称作伪孔洞）也被填充 ...
#include <iostream> #include <pcl/io/vtk_lib_io.h> #include <pcl/io/ply_io.h> #include <pcl/visualization/cloud_viewer.h> int main() { // 1. 读取VTK文件 pcl::PointCloud<pcl::PointXYZ>::Ptr cloud(new pcl::PointCloud<pcl::PointXYZ>); if (pcl::io::loadPolygonFileVTK("input.vtk", *cloud) == -1) { std::cerr << "Failed to load VTK file" << std::endl; return -1; } std::cout << "Loaded " << cloud->size() << " points" << std::endl; // 2. 保存为PLY格式 pcl::PLYWriter writer; writer.write("output.ply", *cloud); std::cout << "Saved to output.ply" << std::endl; // 3. 可视化点云 pcl::visualization::CloudViewer viewer("VTK to PLY Viewer"); viewer.showCloud(cloud); // 保持可视化窗口 while (!viewer.wasStopped()) {} return 0; } 4. 关键注意事项 点数据类型适配： 如果包含RGB信息，使用pcl::PointXYZRGB类型 如果包含法向量，使用pcl::PointNormal类型 可视化增强选项： C++ // 添加坐标轴参考系 viewer.addCoordinateSystem(1.0); // 1米尺度 // 设置点云渲染属性 viewer.setPointCloudRenderingProperties( pcl::visualization::PCL_VISUALIZER_POINT_SIZE, 2, // 点大小 "cloud" // 点云ID ); 不能直接使用cloud的地址存储vtk的数据
最新发布
09-24
<think> 我们面临的问题是：使用PCL库提取VTK点坐标转存为PLY格式并可视化时，不能直接用cloud地址存储VTK数据（根据引用[2]中的代码，我们看到PCL可以直接读取PLY，但VTK文件需要特殊处理）。同时，我们需要考虑点云数据的不同类型（如RGB、法向量）以及可视化增强。 根据引用[1]和[2]，我们知道PCL可以处理PLY文件，并且可以通过VTK进行可视化。引用[3]也提到VTK可以读取和保存PLY，但需要注意颜色信息的保存。引用[4]则详细介绍了PLY格式的存储内容。 用户的问题中特别提到： 1. 不能直接用cloud地址存储VTK数据：这意味着我们需要正确地读取VTK文件，并将其转换为PCL点云。 2. 点数据类型适配：当点云包含RGB时使用pcl::PointXYZRGB，包含法向量时使用pcl::PointNormal。 3. 可视化增强：添加坐标轴参考系、设置点云渲染属性。 解决方案： 1. 读取VTK文件：使用PCL提供的`pcl::io::loadPolygonFileVTK`函数，它可以将VTK文件中的点云数据读取到PCL的点云结构中。但是注意，这个函数原本是用于读取多边形网格的，但也能提取其中的点云。 2. 点云类型适配：我们需要根据VTK文件中包含的数据类型来选择不同的点云类型。这要求我们在读取之前知道点云的数据结构（是否有颜色、法向量等），或者动态判断。然而，PCL的`loadPolygonFileVTK`函数默认读取为`pcl::PointXYZ`，如果我们需要其他类型，则需要额外的处理。 由于`loadPolygonFileVTK`在读取时不会自动识别颜色和法向量，我们需要手动处理。我们可以通过VTK库直接读取VTK文件，然后提取点、颜色和法向量，再填充到PCL点云中。 3. 保存为PLY：使用PCL的`pcl::io::savePLYFile`或者`pcl::PLYWriter`。注意，如果要保存颜色，必须使用支持颜色的点云类型（如`pcl::PointXYZRGB`），并且引用[3]提到，在VTK中保存颜色需要额外设置，同样在PCL中保存颜色到PLY也需要点云中包含颜色信息。 4. 可视化增强：在可视化时，我们可以添加坐标系（`addCoordinateSystem`）和设置点云渲染属性（如点大小、颜色等）。 因此，步骤调整为： 步骤1：使用VTK库读取VTK文件（因为PCL的loadPolygonFileVTK功能有限，不能直接读取颜色和法向量）。 步骤2：从VTK数据中提取点、颜色和法向量（如果有的话）。 步骤3：根据提取的数据创建PCL点云（可能是`pcl::PointXYZ`、`pcl::PointXYZRGB`、`pcl::PointNormal`或自定义点类型）。 步骤4：将点云保存为PLY格式。 步骤5：使用PCL的可视化工具显示点云，并添加坐标轴和设置渲染属性。 具体实现： 1. 使用VTK读取VTK文件： vtkSmartPointer<vtkPolyDataReader> reader = vtkSmartPointer<vtkPolyDataReader>::New(); reader->SetFileName("input.vtk"); reader->Update(); vtkPolyData* polyData = reader->GetOutput(); 2. 检查polyData中的点、颜色和法向量： - 点：polyData->GetPoints() 获取点的坐标。 - 颜色：通过GetPointData()->GetArray("Colors")获取颜色数组（如果有的话）。 - 法向量：通过GetPointData()->GetNormals()获取法向量（如果有的话）。 3. 根据是否存在颜色和法向量，选择不同的点云类型。这里我们考虑三种情况： - 只有坐标：pcl::PointXYZ - 坐标+颜色：pcl::PointXYZRGB - 坐标+法向量：pcl::PointNormal 4. 创建PCL点云并填充数据。 5. 保存为PLY：注意PLY格式支持颜色和法向量，所以要根据点云类型选择保存方式。PCL的PLYWriter会自动保存点云中的字段。 6. 可视化：使用PCL的可视化器，并添加坐标轴和设置点云渲染属性。 注意：如果同时存在颜色和法向量，我们可以选择自定义点类型，或者只保存其中一种（因为PCL的标准点类型中没有同时包含颜色和法向量的）。这里为了简单，我们假设点云只包含颜色或只包含法向量。 由于自定义点类型比较复杂，我们在此只处理两种常见情况：有颜色和有法向量。 下面是具体代码示例（分为三种情况处理：只有坐标、有颜色、有法向量）： </think> ### PCL库提取VTK点坐标转存PLY格式可视化解决方案 #### 1. **问题分析** 当使用`pcl::io::loadPolygonFileVTK`直接读取VTK文件时，可能遇到无法正确处理不同类型点云数据（如RGB、法向量）的问题[^2]。主要挑战在于： - 直接使用`pcl::PointXYZ`会丢失颜色/法向量信息 - VTK数据结构与PCL点云类型不完全兼容 - 可视化时需要增强显示效果 #### 2. **完整解决方案** ```cpp #include <pcl/io/vtk_lib_io.h> #include <pcl/io/ply_io.h> #include <pcl/visualization/pcl_visualizer.h> #include <vtkPolyData.h> // 通用点云类型定义 template <typename PointT> void processVTK(const std::string& input_path, const std::string& output_path) { // 读取VTK文件 vtkSmartPointer<vtkPolyData> polyData; vtkSmartPointer<vtkPolyDataReader> reader = vtkSmartPointer<vtkPolyDataReader>::New(); reader->SetFileName(input_path.c_str()); reader->Update(); polyData = reader->GetOutput(); // 转换为PCL点云 typename pcl::PointCloud<PointT>::Ptr cloud(new pcl::PointCloud<PointT>); pcl::io::vtkPolyDataToPointCloud(polyData, *cloud); // 关键转换接口 // 保存为PLY格式 pcl::PLYWriter writer; writer.write(output_path, *cloud); std::cout << "保存点云: " << cloud->size() << " 个点" << std::endl; // 增强可视化 pcl::visualization::PCLVisualizer::Ptr viewer( new pcl::visualization::PCLVisualizer("3D Viewer")); viewer->setBackgroundColor(0, 0, 0); viewer->addCoordinateSystem(1.0, "global"); // 添加1米尺度坐标系 // 根据点云类型设置渲染属性 if constexpr (std::is_same_v<PointT, pcl::PointXYZRGB>) { viewer->addPointCloud<pcl::PointXYZRGB>(cloud, "cloud_rgb"); viewer->setPointCloudRenderingProperties( pcl::visualization::PCL_VISUALIZER_POINT_SIZE, 3, "cloud_rgb"); } else if constexpr (std::is_same_v<PointT, pcl::PointNormal>) { viewer->addPointCloudNormals<pcl::PointNormal>(cloud, 1, 0.1, "cloud_normals"); viewer->setPointCloudRenderingProperties( pcl::visualization::PCL_VISUALIZER_COLOR, 1.0, 0.5, 0.0, "cloud_normals"); } else { viewer->addPointCloud<pcl::PointXYZ>(cloud, "cloud_xyz"); viewer->setPointCloudRenderingProperties( pcl::visualization::PCL_VISUALIZER_POINT_SIZE, 2, "cloud_xyz"); } // 添加文本说明 viewer->addText("VTK转PLY可视化", 10, 20, 20, 1, 1, 1, "text"); // 交互式显示 while (!viewer->wasStopped()) { viewer->spinOnce(100); } } int main() { // 根据数据类型选择模板实例化 processVTK<pcl::PointXYZ>("model.vtk", "output_xyz.ply"); // XYZ点云 processVTK<pcl::PointXYZRGB>("color_model.vtk", "output_rgb.ply"); // RGB点云 processVTK<pcl::PointNormal>("normal_model.vtk", "output_normal.ply"); // 法向量点云 return 0; } ``` #### 3. **关键实现说明** 1. **VTK到PCL的转换** 使用`pcl::io::vtkPolyDataToPointCloud`替代直接加载，确保数据结构兼容： ```cpp vtkSmartPointer<vtkPolyData> polyData = ...; pcl::PointCloud<PointT> cloud; pcl::io::vtkPolyDataToPointCloud(polyData, cloud); ``` 2. **点云类型适配** - **RGB点云**：使用`pcl::PointXYZRGB`类型存储 ```cpp point.r = polyData->GetPointData()->GetScalars("Colors")->GetComponent(i, 0); point.g = polyData->GetPointData()->GetScalars("Colors")->GetComponent(i, 1); point.b = polyData->GetPointData()->GetScalars("Colors")->GetComponent(i, 2); ``` - **法向量点云**：使用`pcl::PointNormal`类型存储 ```cpp vtkDataArray* normals = polyData->GetPointData()->GetNormals(); point.normal_x = normals->GetComponent(i, 0); point.normal_y = normals->GetComponent(i, 1); point.normal_z = normals->GetComponent(i, 2); ``` 3. **可视化增强** ```cpp // 添加坐标系 (1米尺度) viewer->addCoordinateSystem(1.0, "global"); // 设置点云属性 (RGB点云示例) viewer->setPointCloudRenderingProperties( pcl::visualization::PCL_VISUALIZER_POINT_SIZE, 3, // 点大小 "cloud_rgb" // 点云ID ); // 添加法向量显示 (0.1米长度) viewer->addPointCloudNormals<pcl::PointNormal>(cloud, 1, 0.1, "normals"); ``` #### 4. **常见问题处理** 1. **颜色信息丢失** 确保VTK文件包含颜色数组且命名为"Colors"[^3]： ```cpp vtkUnsignedCharArray* colors = vtkUnsignedCharArray::SafeDownCast( polyData->GetPointData()->GetArray("Colors")); ``` 2. **法向量方向异常** 添加可视化参考轴辅助检查： ```cpp viewer->addCoordinateSystem(0.5, "origin"); // 在原点添加小坐标系 ``` 3. **大型点云优化** 使用八叉树加速渲染： ```cpp viewer->setPointCloudRenderingProperties( pcl::visualization::PCL_VISUALIZER_OCTREE_RENDERING, true, "cloud_xyz" ); ``` #### 5. **不同数据类型的处理示例** | 数据类型 | PCL点云类型 | 关键处理步骤 | |----------------|---------------------|----------------------------------| | 纯坐标点云 | `pcl::PointXYZ` | 直接转换点坐标 | | RGB彩色点云 | `pcl::PointXYZRGB` | 提取并转换颜色数组 | | 含法向量点云 | `pcl::PointNormal` | 读取法向量数据并填充 | | 含强度值点云 | `pcl::PointXYZI` | 提取标量字段作为强度值 | 此方案通过VTK原生接口读取数据，配合PCL类型模板，有效解决了直接加载导致的数据丢失问题[^2]。
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符  | 博主筛选后可见
 
红包 添加红包
表情包 插入表情
表情包 代码片
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值