VTK_Learning_三角网格体积、表面积、测地距离、包围盒

2023-10-13 08:30
文章标签 体积 距离 包围 learning 表面积 三角 网格 vtk 测地

本文主要是介绍VTK_Learning_三角网格体积、表面积、测地距离、包围盒，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

1.基本图形操作意义

图形处理，比如图形平滑、多分辨率分析、特征提取等都离不开一些基本的图形操作。掌握这些基本的图形操作有助于理解和深入学习图形处理和分析方法。
VTK中提供了多种图形的基本操作，其中最简单的是点的欧氏距离计算，可以使用vtkMath进行计算，也可以直接计算向量的模。一些图元类提供了许多可以方便使用的静态函数，如
vtkLine提供了点与线间的距离计算；
vtkTriangle提供了面积、外接圆、法向量的计算，点与三角形位置关系判断等；
vtkTetra中实现了四面体体积，重心计算等。
有了这些函数，可以实现很多其他功能，如计算一个三角网格模型的表面积，只需要遍历每个三角形单元并计算其面积即可。
另外，还有一个办法是vtkMassProperties。这个类可以实现三角网格的表面积和体积计算，但是要求网格必须是封闭的三角形网格数据。网格的封闭性计算在后面会有更加详细的讨论。对于非三角形网格，需要先将网格转换为三角形网格。vtkTriangleFilter可以实现多边形网格数据向三角形网格数据转换。
 

2.三角网格模型面积、体积计算 

利用vtkMassProperties计算三角网格模型面积、体积代码如下：

#include <vtkAutoInit.h>
VTK_MODULE_INIT(vtkRenderingOpenGL);
VTK_MODULE_INIT(vtkInteractionStyle);#include <vtkSmartPointer.h>
#include <vtkCubeSource.h>
#include <vtkTriangleFilter.h> //其他网格类型转换成三角网格类型
#include <vtkMassProperties.h> //计算三角网格的基本属性 面积。体积等
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkProperty.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>int main()
{vtkSmartPointer<vtkCubeSource> cubeSource =vtkSmartPointer<vtkCubeSource>::New(); //vtkPolyData类型数据cubeSource->Update();vtkSmartPointer<vtkTriangleFilter> triFilter =vtkSmartPointer<vtkTriangleFilter>::New();triFilter->SetInputData(cubeSource->GetOutput());triFilter->Update();vtkSmartPointer<vtkMassProperties> massProp =vtkSmartPointer<vtkMassProperties>::New();massProp->SetInputData(triFilter->GetOutput());float Volume = massProp->GetVolume();float SurfaceArea = massProp->GetSurfaceArea();float maxArea = massProp->GetMaxCellArea();float minArea = massProp->GetMinCellArea();std::cout << "the Volume   : " << Volume << std::endl;std::cout << "Surface Area : " << SurfaceArea << std::endl;std::cout << "MaxAreaofCell: " << maxArea << std::endl;std::cout << "MinAreaofCell: " << minArea << std::endl;///vtkSmartPointer<vtkPolyDataMapper> mapper =vtkSmartPointer<vtkPolyDataMapper>::New();mapper->SetInputData(triFilter->GetOutput());vtkSmartPointer<vtkActor> actor =vtkSmartPointer<vtkActor>::New();actor->SetMapper(mapper);actor->GetProperty()->SetColor(0, 0, 1);actor->GetProperty()->SetEdgeColor(1, 0, 0);actor->GetProperty()->SetEdgeVisibility(1);vtkSmartPointer<vtkRenderer> render =vtkSmartPointer<vtkRenderer>::New();render->AddActor(actor);render->SetBackground(0, 0, 0);vtkSmartPointer<vtkRenderWindow> rw =vtkSmartPointer<vtkRenderWindow>::New();rw->AddRenderer(render);rw->SetSize(480, 420);rw->SetWindowName("Calculating Area and Volume of Triangle grid");vtkSmartPointer<vtkRenderWindowInteractor> rwi =vtkSmartPointer<vtkRenderWindowInteractor>::New();rwi->SetRenderWindow(rw);rwi->Initialize();rwi->Start();return 0;
}

其输出结果为：

 

3.三维网格测地距离 

对于三维网格模型来讲，测地距离也是一种重要的距离度量。与欧氏距离不同，一个三维模型上的亮点测地距离是指沿着模型表面两者之间的最短距离。测地距离通常采用Dijkstra算法类近似求解。VTK中的vtkDijkstraGraphGeodesicPath类就可以实现测地距离的求解。该类的使用如下例所示：
 

#include <vtkAutoInit.h>
VTK_MODULE_INIT(vtkRenderingOpenGL);
VTK_MODULE_INIT(vtkInteractionStyle); #include <vtkSmartPointer.h>
#include <vtkSphereSource.h>
#include <vtkDijkstraGraphGeodesicPath.h>
#include <vtkProperty.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>int main()
{vtkSmartPointer<vtkSphereSource> sphereSource =vtkSmartPointer<vtkSphereSource>::New();sphereSource->Update();vtkSmartPointer<vtkDijkstraGraphGeodesicPath> dijstra =vtkSmartPointer<vtkDijkstraGraphGeodesicPath>::New();dijstra->SetInputData(sphereSource->GetOutput());dijstra->SetStartVertex(0);dijstra->SetEndVertex(10);dijstra->Update();///vtkSmartPointer<vtkPolyDataMapper> mapper =vtkSmartPointer<vtkPolyDataMapper>::New();mapper->SetInputData(sphereSource->GetOutput());vtkSmartPointer<vtkPolyDataMapper> pathMapper =vtkSmartPointer<vtkPolyDataMapper>::New();pathMapper->SetInputData(dijstra->GetOutput());vtkSmartPointer<vtkActor> actor =vtkSmartPointer<vtkActor>::New();actor->SetMapper(mapper);vtkSmartPointer<vtkActor> pathActor =vtkSmartPointer<vtkActor>::New();pathActor->SetMapper(pathMapper);pathActor->GetProperty()->SetColor(1, 0, 0);pathActor->GetProperty()->SetLineWidth(5);vtkSmartPointer<vtkRenderer> renderer =vtkSmartPointer<vtkRenderer>::New();renderer->AddActor(actor);renderer->AddActor(pathActor);renderer->SetBackground(0, 0, 0);vtkSmartPointer<vtkRenderWindow> rw =vtkSmartPointer<vtkRenderWindow>::New();rw->AddRenderer(renderer);rw->SetSize(640, 480);rw->SetWindowName("Calculating Geodesic Path");vtkSmartPointer<vtkRenderWindowInteractor> rwi =vtkSmartPointer<vtkRenderWindowInteractor>::New();rwi->SetRenderWindow(rw);rwi->Initialize();rwi->Start();return 0;
}

 

4.三维图像的包围盒 

 包围盒是指能够包围模型的最小立方体，常常用于模型的碰撞测量中。vtkPolyData中定义了函数GetBounds()来获取包围盒的参数。即三个坐标轴方向上的最大、最小值。仅仅获取这些参数并不直观，有时候还需要显示包围盒。vtkOutlineFilter则提供了一个方便的方法来生成包围盒，其输入为一个vtkPolyData莫形数据，输出同样为一个vtkPolyData类型数据，因此非常容易进行可视化显示。
实例的参考代码如下所示：
 

#include <vtkAutoInit.h>
VTK_MODULE_INIT(vtkRenderingOpenGL);
VTK_MODULE_INIT(vtkInteractionStyle);#include <vtkSmartPointer.h>
#include <vtkSphereSource.h>
#include <vtkProperty.h>
#include <vtkOutlineFilter.h>
#include <vtkPolyData.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>int main()
{vtkSmartPointer<vtkSphereSource> sphereSource =vtkSmartPointer<vtkSphereSource>::New();sphereSource->SetCenter(0.0, 0.0, 0.0);sphereSource->SetRadius(5.0);sphereSource->Update();vtkPolyData* sphere = sphereSource->GetOutput();vtkSmartPointer<vtkOutlineFilter> outline =vtkSmartPointer<vtkOutlineFilter>::New();outline->SetInputData(sphere);outline->Update();//算法执行完毕，必须更新！！！vtkSmartPointer<vtkPolyDataMapper> mapper =vtkSmartPointer<vtkPolyDataMapper>::New();mapper->SetInputData(sphere);vtkSmartPointer<vtkPolyDataMapper> outlineMapper =vtkSmartPointer<vtkPolyDataMapper>::New();outlineMapper->SetInputData(outline->GetOutput());vtkSmartPointer<vtkActor> actor =vtkSmartPointer<vtkActor>::New();actor->SetMapper(mapper);vtkSmartPointer<vtkActor> outlineActor =vtkSmartPointer<vtkActor>::New();outlineActor->SetMapper(outlineMapper);outlineActor->GetProperty()->SetColor(0, 1, 0);outlineActor->GetProperty()->SetLineWidth(3);vtkSmartPointer<vtkRenderer> renderer =vtkSmartPointer<vtkRenderer>::New();renderer->AddActor(actor);renderer->AddActor(outlineActor);renderer->SetBackground(0, 0, 0);vtkSmartPointer<vtkRenderWindow> rw =vtkSmartPointer<vtkRenderWindow>::New();rw->AddRenderer(renderer);rw->SetSize(640, 480);;rw->SetWindowName("PolyData Bounding Box");rw->Render();vtkSmartPointer<vtkRenderWindowInteractor> rwi =vtkSmartPointer<vtkRenderWindowInteractor>::New();rwi->SetRenderWindow(rw);rwi->Start();return 0;
}

其输出结果如下图所示：

这篇关于VTK_Learning_三角网格体积、表面积、测地距离、包围盒的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


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

相关文章

简单的Q-learning|小明的一维世界(3)

简单的Q-learning|小明的一维世界(3)

简单的Q-learning|小明的一维世界(1) 简单的Q-learning|小明的一维世界(2) 一维的加速度世界 这个世界，小明只能控制自己的加速度，并且只能对加速度进行如下三种操作：增加1、减少1、或者不变。所以行动空间为： { u 1 = − 1 , u 2 = 0 , u 3 = 1 } \{u_1=-1, u_2=0, u_3=1\} {u1​=−1,u2​=0,u3​=1}
阅读更多...
简单的Q-learning|小明的一维世界(2)

简单的Q-learning|小明的一维世界(2)

上篇介绍了小明的一维世界模型 、Q-learning的状态空间、行动空间、奖励函数、Q-table、Q table更新公式、以及从Q值导出策略的公式等。最后给出最简单的一维位置世界的Q-learning例子，从给出其状态空间、行动空间、以及稠密与稀疏两种奖励函数的设置方式。下面将继续深入，GO！ 一维的速度世界 这个世界，小明只能控制自己的速度，并且只能对速度进行如下三种操作：增加1、减
阅读更多...
线性代数|机器学习-P35距离矩阵和普鲁克问题

线性代数|机器学习-P35距离矩阵和普鲁克问题

文章目录 1. 距离矩阵2. 正交普鲁克问题3. 实例说明 1. 距离矩阵 假设有三个点 x 1 , x 2 , x 3 x_1,x_2,x_3 x1​,x2​,x3​,三个点距离如下： ∣ ∣ x 1 − x 2 ∣ ∣ 2 = 1 , ∣ ∣ x 2 − x 3 ∣ ∣ 2 = 1 , ∣ ∣ x 1 − x 3 ∣ ∣ 2 = 6 \begin{equation} ||x
阅读更多...
初步了解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
阅读更多...
005:VTK世界坐标系中的相机和物体

005:VTK世界坐标系中的相机和物体

VTK医学图像处理---世界坐标系中的相机和物体 左侧是成像结果                                                    右侧是世界坐标系中的相机与被观察物体 目录 VTK医学图像处理---世界坐标系中的相机和物体 简介 1 在三维空间中添加坐标系 2 世界坐标系中的相机 3 世界坐标系中vtkImageData的参数 总结：
阅读更多...
模拟退火求n个点到某点距离和最短

模拟退火求n个点到某点距离和最短

/*找出一个点使得这个店到n个点的最长距离最短，即求最小覆盖圆的半径用一个点往各个方向扩展，如果结果更优，则继续以当前步长扩展，否则缩小步长*/#include<stdio.h>#include<math.h>#include<string.h>const double pi = acos(-1.0);struct point {double x,y;}p[1010];int
阅读更多...
黑神话：悟空》增加草地绘制距离MOD使游戏场景看起来更加广阔与自然，增强了游戏的沉浸式体验

黑神话：悟空》增加草地绘制距离MOD使游戏场景看起来更加广阔与自然，增强了游戏的沉浸式体验

《黑神话：悟空》增加草地绘制距离MOD为玩家提供了一种全新的视觉体验，通过扩展游戏中草地的绘制距离，增加了场景的深度和真实感。该MOD通过增加草地的绘制距离，使游戏场景看起来更加广阔与自然，增强了游戏的沉浸式体验。 增加草地绘制距离MOD安装 1、在%userprofile%AppDataLocalb1SavedConfigWindows目录下找到Engine.ini文件。 2、使用记事本编辑
阅读更多...
SimD：基于相似度距离的小目标检测标签分配

SimD：基于相似度距离的小目标检测标签分配

摘要 https://arxiv.org/pdf/2407.02394 由于物体尺寸有限且信息不足，小物体检测正成为计算机视觉领域最具挑战性的任务之一。标签分配策略是影响物体检测精度的关键因素。尽管已经存在一些针对小物体的有效标签分配策略，但大多数策略都集中在降低对边界框的敏感性以增加正样本数量上，并且需要设置一些固定的超参数。然而，更多的正样本并不一定会带来更好的检测结果，事实上，过多的正样本
阅读更多...
”CSS 网格“二维布局系统（补充）——WEB开发系列32

”CSS 网格“二维布局系统（补充）——WEB开发系列32

CSS 网格布局是一种二维布局系统，用于网页设计。通过使用网格，你可以将内容以行和列的形式进行排列。此外，网格布局还能够简便地实现一些复杂的布局结构。 一、什么是网格布局？ CSS网格布局是一种二维布局系统，它允许我们创建复杂的网页布局，既可以处理行也可以处理列。与传统的布局方法不同，网格布局将网页分成多个可控的区域，这些区域可以任意排列、对齐和调整大小。网格布局使得创建灵活且响应
阅读更多...
Matlab）实现HSV非等间隔量化--相似判断：欧式距离--输出图片-

Matlab）实现HSV非等间隔量化--相似判断：欧式距离--输出图片-

%************************************************************************** %                                 图像检索——提取颜色特征 %HSV空间颜色直方图(将RGB空间转化为HS
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2