ITK学习笔记：1.Data presentation-1.3Mesh网格（三）：拓扑和k复形

最新推荐文章于 2022-09-25 11:18:37 发布
最新推荐文章于 2022-09-25 11:18:37 发布
阅读量364 收藏
点赞数
分类专栏： ITK指北
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44671418/article/details/107032972
版权

Mesh类支持形式化拓扑的表示。特别地，k复形的概念可以在网格中正确地表示。
k复形是一种拓扑结构，其中对于每一个维数N的单元格，其边界面为维数N−1的单元格也属于该结构。

下面演示如何使用网格实例化K-Complex结构。

  • 该实例结构由一个四面体、四个三角形面、六个线边和四个顶点组成
#include "itkMesh.h" 
#include "itkVertexCell.h" 
#include "itkLineCell.h" 
#include "itkTriangleCell.h" 
#include "itkTetrahedronCell.h"

//定义PixelType并实例化网格类型。
//注意，在这种情况下，空间的维度是3。
typedef float PixelType; 
typedef itk::Mesh< PixelType, 3 > MeshType;

//现在可以使用从网格中获取的特征实例化单元格类型
typedef MeshType::CellType CellType; 
typedef itk::VertexCell< CellType > VertexType; 
typedef itk::LineCell< CellType > LineType; 
typedef itk::TriangleCell< CellType > TriangleType; 
typedef itk::TetrahedronCell< CellType > TetrahedronType;

//网格被创建，与顶点相关联的点被插入。
//正四面体节点的几何坐标可以通过从正立方体中取出每一个节点来得到。
MeshType::Pointer mesh = MeshType::New();
MeshType::PointType point0;
MeshType::PointType point1;
MeshType::PointType point2; 
MeshType::PointType point3;
point0[0] = -1; point0[1] = -1; point0[2] = -1;
point1[0] = 1; point1[1] = 1; point1[2] = -1;
point2[0] = 1; point2[1] = -1; point2[2] = 1;
point3[0] = -1; point3[1] = 1; point3[2] = 1;
mesh->SetPoint( 0, point0 ); mesh->SetPoint( 1, point1 ); mesh->SetPoint( 2, point2 ); mesh->SetPoint( 3, point3 );
/*
注意，网格中的点和itk::VertexCell概念之间有一个重要的区别。
VertexCell是维度为零的单元格。它与点的主要区别在于，cell可以感知到与其他cell的邻域关系。
而点不知道cell的存在。事实上，从纯粹的拓扑角度来看，网格中点的坐标是完全无关的。它们也可能完全不存在于网格结构中。
VertexCells对于表示k -复形上有邻域关系的完整集合是必要的。
*/

//现在我们继续创建单元格，将它们与点关联起来，并将它们插入网格中。
//从四面体开始，我们编写以下代码
CellType::CellAutoPointer cellpointer;
cellpointer.TakeOwnership( new TetrahedronType );
cellpointer->SetPointId( 0, 0 ); 
cellpointer->SetPointId( 1, 1 );
cellpointer->SetPointId( 2, 2 ); 
cellpointer->SetPointId( 3, 3 );
 mesh->SetCell( 0, cellpointer );
//现在创建了四个三角形面并与网格相关联。第一个三角形连接点0,1,2。
cellpointer.TakeOwnership( new TriangleType ); 
cellpointer->SetPointId( 0, 0 );
cellpointer->SetPointId( 1, 1 ); 
cellpointer->SetPointId( 2, 2 ); 
mesh->SetCell( 1, cellpointer );
//第二个三角形连接点0,2,3
cellpointer.TakeOwnership( new TriangleType ); 
cellpointer->SetPointId( 0, 0 ); 
cellpointer->SetPointId( 1, 2 ); 
cellpointer->SetPointId( 2, 3 ); 
mesh->SetCell( 2, cellpointer );
//第三个
cellpointer.TakeOwnership( new TriangleType ); 
cellpointer->SetPointId( 0, 0 ); 
cellpointer->SetPointId( 1, 3 ); 
cellpointer->SetPointId( 2, 1 ); 
mesh->SetCell( 3, cellpointer );
//第四个三角形连接点3,2,1 
cellpointer.TakeOwnership( new TriangleType ); 
cellpointer->SetPointId( 0, 3 ); 
cellpointer->SetPointId( 1, 2 ); 
cellpointer->SetPointId( 2, 1 ); 
mesh->SetCell( 4, cellpointer );
/*注意每次CellAutoPointer是如何重用的。
提醒:当itk::AutoPointer作为SetCell()方法的参数传递时，它失去了对单元格的所有权。
AutoPointer通过使用TakeOwnership()方法连接到一个新的单元格。*/

//现在继续构建K-Complex，在四面体的边缘上创建六条线。
cellpointer.TakeOwnership( new LineType ); 
cellpointer->SetPointId( 0, 0 ); 
cellpointer->SetPointId( 1, 1 ); 
mesh->SetCell( 5, cellpointer );

cellpointer.TakeOwnership( new LineType );
 cellpointer->SetPointId( 0, 1 ); 
 cellpointer->SetPointId( 1, 2 ); 
 mesh->SetCell( 6, cellpointer );
 
cellpointer.TakeOwnership( new LineType ); 
cellpointer->SetPointId( 0, 2 ); 
cellpointer->SetPointId( 1, 0 ); 
mesh->SetCell( 7, cellpointer );

cellpointer.TakeOwnership( new LineType ); 
cellpointer->SetPointId( 0, 1 ); 
cellpointer->SetPointId( 1, 3 ); 
mesh->SetCell( 8, cellpointer );

cellpointer.TakeOwnership( new LineType ); 
cellpointer->SetPointId( 0, 3 ); 
cellpointer->SetPointId( 1, 2 ); 
mesh->SetCell( 9, cellpointer );

cellpointer.TakeOwnership( new LineType );
cellpointer->SetPointId( 0, 3 ); 
cellpointer->SetPointId( 1, 0 ); 
mesh->SetCell( 10, cellpointer );

//最后，创建并插入由itk::VertexCell表示的零维单元格
cellpointer.TakeOwnership( new VertexType );
cellpointer->SetPointId( 0, 0 );
mesh->SetCell( 11, cellpointer );

cellpointer.TakeOwnership( new VertexType );
cellpointer->SetPointId( 0, 1 ); 
mesh->SetCell( 12, cellpointer );

cellpointer.TakeOwnership( new VertexType ); 
cellpointer->SetPointId( 0, 2 ); 
mesh->SetCell( 13, cellpointer );

cellpointer.TakeOwnership( new VertexType );
cellpointer->SetPointId( 0, 3 ); 
mesh->SetCell( 14, cellpointer );

此时网格包含4个点和从0到14枚举的15个单元格。

  • 可以使用PointContainer迭代器访问这些点
typedef MeshType::PointsContainer::ConstIterator PointIterator; 
PointIterator pointIterator = mesh->GetPoints()->Begin(); 
PointIterator pointEnd = mesh->GetPoints()->End();
while( pointIterator != pointEnd ) 
{ 
std::cout << pointIterator.Value() << std::endl; 
++pointIterator;
 }
  • 可以使用CellsContainer迭代器访问单元格
typedef MeshType::CellsContainer::ConstIterator CellIterator;
CellIterator cellIterator = mesh->GetCells()->Begin(); 
CellIterator cellEnd = mesh->GetCells()->End();
while( cellIterator != cellEnd )
 { 
 CellType * cell = cellIterator.Value();
 std::cout << cell->GetNumberOfPoints() << std::endl; 
 ++cellIterator; 
 }

请注意，单元格存储为指向泛型单元格类型的指针，泛型单元格类型是所有特定单元格类的基类。

  • 这意味着在这个级别上我们只能访问CellType中定义的虚拟方法。

可以使用CellType特征中定义的迭代器访问与单元格关联的点标识符。
下面的代码演示了PointIdIterators的使用。PointIdsBegin()方法将迭代器返回到单元格中的第一个点标识符。PointIdsEnd()方法将迭代器返回到单元格中的过期结束点标识符：

typedef CellType::PointIdIterator PointIdIterator;
PointIdIterator pointIditer = cell->PointIdsBegin();
 PointIdIterator pointIdend = cell->PointIdsEnd();
while( pointIditer != pointIdend )
 {
  std::cout << *pointIditer << std::endl; 
  ++pointIditer;
   }

到这里为止，k复形的拓扑还没有完全定义，因为我们只引入了单元格。ITK允许用户显式地定义单元格之间的睦邻关系。很明显，对点标识符的巧妙探索可以让用户找出邻居关系。

  • 例如，共享相同的两个点标识符的两个三角形单元格可能是相邻单元格。这种邻域关系的隐式发现的一些缺点是，它需要计算时间，而且一些应用程序可能不接受同样的假设。

在网格中，邻域关系用边界特征的概念表示。

每个cell都有一个单元标识符的内部列表，指向被认为是它的邻居的其他cell。边界特征按维数进行分类。

  • 例如，一条线有两个维度为零的边界特征，对应于它的两个顶点。一个四面体将具有维度为0、1和2的边界特征，对应于它的四个顶点、六个边和四个三角形面。由用户指定单元格之间的连接。

我们以当前的例子为例，将与单元标识符0相关联的四面体单元格分配给它四个顶点作为零维度的边界。这是通过调用网格类上的SetBoundaryAssignment(()方法来完成的：

MeshType::CellIdentifier cellId = 0; // the tetrahedron
int dimension = 0; // vertices
MeshType::CellFeatureIdentifier featureId = 0;
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 11 );
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 12 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 13 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 14 );

featureId只是一个与特定单元中相同维度的边界单元序列相关联的数字。

  • 例如，一个四面体的零维特征是它的四个顶点。然后，此单元格的零维特性id的范围将从0到3。四面体的一维特征是它的六个边，因此它的一维特征id将从0到5。四面体的二维特征是它的四个三角形面。二维特征id的范围将从0到3。

下表总结：
在这里插入图片描述现在让我们指定这个四面体的一维边界特征。

  • 这些是具有标识符{5、6、7、8、9、10}的行单元格。

  • 注意，特性标识符被重新初始化为零，因为计数独立于每个维度。

cellId = 0; // still the tetrahedron 
dimension = 1; // one-dimensional features = edges 
featureId = 0; // reinitialize the count
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 5 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 6 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 7 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 8 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 9 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 10 );

最后确定四面体的二维边界特征。

  • 这是四个三角形单元格，标识符为{1,2,3,4}。featureId被重置为零，因为特性id独立于每个维度。
cellId = 0; // still the tetrahedron 
dimension = 2; // two-dimensional features = triangles 
featureId = 0; // reinitialize the count
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 1 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 2 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 3 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 4 );

此时，我们可以查询四面体单元格的边界特征信息。例如，可以通过getnumberofborderyfeatures()方法得到每个维度的边界特征个数：

cellId = 0; // still the tetrahedron
MeshType::CellFeatureCount n0; // number of zero-dimensional features 
MeshType::CellFeatureCount n1; // number of one-dimensional features 
MeshType::CellFeatureCount n2; // number of two-dimensional features
n0 = mesh->GetNumberOfCellBoundaryFeatures( 0, cellId ); 
n1 = mesh->GetNumberOfCellBoundaryFeatures( 1, cellId ); 
n2 = mesh->GetNumberOfCellBoundaryFeatures( 2, cellId );

边界分配可以通过 GetBoundaryAssigment()方法重新获得。
例如，可以通过以下代码获得四面体的零维特征：

dimension = 0; 
for(unsigned int b0=0; b0 < n0; b0++) 
{ 
MeshType::CellIdentifier id; 
bool found = mesh->GetBoundaryAssignment( dimension, cellId, b0, &id ); 
if( found ) std::cout << id << std::endl; }

//下面的代码演示了如何为一个三角形面设置边缘边界。
cellId = 2; // one of the triangles 
dimension = 1; // boundary edges 
featureId = 0; // start the count of features
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 7 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 9 ); 
mesh->SetBoundaryAssignment( dimension, cellId, featureId++, 10 );
人工智睿
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
专栏目录
《itk实用demo》-itkMesh
就这样吧
09-12 331
itk中的网格Mesh #include &quot;itkMesh.h&quot; MeshType ::Pointer PointsToMesh(ImageType::Pointer sliceimage) { MeshType::Pointer mesh = MeshType::New(); int size[2]; size[0]=sliceimage-&amp;gt;...
ITK学习笔记：1.Data presentation-1.3Mesh网格（二）：自定义mesh
今天写代码了吗
06-29 343
网格的一般概念集成了许多不同的元素。原则上可以为每个这样的元素使用独立的类型。泛型编程中用于指定概念中涉及的许多不同类型的机制称为traits。它们基本上是与当前类交互的所有类型的列表。 网格通过三个参数模板化。到目前为止，我们只讨论了两个，即像素类型和维度。第三个参数是一个提供网格所需特征集的类。当省略第三个参数时，使用默认的类。这个默认类是itk::DefaultStaticMeshTraits。如果您想定制网格使用的类型，接下来的方法是修改默认特征，并将它们作为网格类实例化的第三个参数提供。有两种方法
ITK学习笔记：1.Data presentation-1.3Mesh网格（七）：访问单元
今天写代码了吗
07-01 188
为了方便访问特定的单元类型，itk::Mesh中内置了一个机制。该机制基于Visitor Pattern。 该模式旨在促进遍历共享公共基类的异构对象列表的过程。 //使用CellVisitor机制的第一个要求是包含CellInterfaceVisitor头文件 #include "itkCellInterfaceVisitor.h" typedef float PixelType; typedef itk::Mesh< PixelType, 3 > MeshType; typedef M
ITK学习笔记：1.Data presentation-1.3Mesh网格（五）：简化网格的创建
今天写代码了吗
06-30 266
Mesh类是非常通用和灵活的。 可以通过itk::AutomaticTopologyMeshSource类，根据我们所添加的单元自动生成一个显性复合体。 它包括了所有的边界信息，因此产生的网格可以很容易地遍历。 它合并了所有共享的边、顶点和面，因此几何特征不会重复出现。 本文展示如何使用AutomaticTopologyMeshSource实例化表示K-Complex的网格。我们将首先从上文中生成相同的四面体，之后我们将添加一个空心的四面体来说明网格源的一些额外特性。 首先包含需要的头文件： #incl
ITK 实例 （四）
weidadeyaqhhya的博客
07-19 331
ITK 官方网址：https://itk.org/ITK/project/project.html （四）ITK 中的网格 itk::Mesh 类用来表示空间中的形状。它派生自 itk::PoitnSet 类，因此它集成了与点相关的所有功能以及与点相关的所有像素数据的访问。网格类是 N 维的，它具有很大的灵活性。在实际应用中，网格类也可以被看作是一个点集，其中添加了许多不同维度和形状的单元格...
itk软件指南：介绍与开发准则-翻译版.pdf
06-14
TK（ Insight Segmentation and Registration Toolkit）是美国国家卫生院下属的国立医学图书馆开发的一款医学图像处理软件包，是一个开源的、跨平台的影像分析扩展软件工具。
itk-vtk-viewer:使用itk.js和vtk.js的2D 3D Web图像，网格和点集查看器
03-09
ITK / VTK Viewer是用于医学和科学图像，网格和点集可视化的开源软件系统。 报告问题 如果您想讨论错误或可能的改进： 如果您有补丁程序，请阅读文档。 在打开一个条目。 要求 通常，ITK / VTK Viewer尝试尽可能地...
PyPI 官网下载 | itk_filtering-5.1.1.post1-cp38-cp38-win_amd64.whl
02-07
《PyPI官网下载的itk_filtering-5.1.1.post1-cp38-cp38-win_amd64.whl：Python图像处理库的深入解析》 在Python编程环境中，PyPI（Python Package Index）是全球最大的Python软件包仓库，开发者可以在这里找到并...
simpleitk-feedstock：simpleitk的conda-smithy存储库
02-25
关于simpleitk 主页：http :: / 软件包许可证：Apache-2.0 原料许可证： 摘要：用于图像配准和分段的Insight Toolkit的简化界面 当前构建状态 蔚蓝 变体 状态 linux_64_python3.6 .____ cpython linux_64_...
Python库 | itk_core-5.2.0.post1-cp38-cp38-win_amd64.whl
03-21
《Python库itk_core详解与应用》 在Python编程领域，库是不可或缺的资源，它们提供了丰富的功能，帮助开发者高效地完成各种任务。本文将详细探讨一个名为`itk_core`的Python库，该库在版本`5.2.0.post1`中进行了...
网格剖切功能强大的视图插值
03-05
现有的视图插值方法（如DIBR）需要精确的视差图，这极大地限制了它们的应用。 为此，我们提出了一种新颖的基于网格的视图插值算法，该算法能够合成视觉匹配的虚拟视图，并通过立体匹配算法估计出粗糙的视差图。 我们采用edgeaware.mesh切割方法来显式处理遮挡并保留尖锐的深度不连续性。 在Middlebridge数据集和3D-HEVC测试序列上进行的实验表明，该方法在视觉质量和PSNR方面优于DIBR和基于最新网格的视图插值算法
Java求geometry的面积最小外接矩形
qq_40985985的博客
09-25 2744
Java求geometry的面积最小外接矩形 geom.getEnvelope() 得到外接矩形，不一定是面积最小；可以对多边形的每一条边求外接矩形，然后比较得到面积最小外接矩形。这篇博客将分为3步进行求解； 1. 获取Envelope的外接矩形，获取面积； 2. 获取凸包，旋转获取每一条边对应的外接矩形，得到凸包面积最小外接矩形； 3. 将1获取的矩形与2获取的矩形面积进行比较，得到最终的面积最小外接矩形；
ITK学习笔记：1.Data presentation-1.3Mesh网格（一）基操：创建mesh，嵌入cell，管理数据
今天写代码了吗
06-25 438
1.创建网格 网格类用于表示空间中的形状。 它源于itk::PointSet类，因此继承了与点相关的所有功能以及对与点关联的像素数据的访问。 Mesh类也是n维的，这使得它的使用具有很大的灵活性。 实际上，Mesh类可以看作是一个点集，其中添加了许多不同维度和形状的单元格(也称为元素)。网格中的单元格使用它们的点标识符根据现有的点定义 与PointSet相同，ITK中有静态和动态两种基本的网格样式。 第一个用于使用点的数量可预知,并不需要对其修改操作。 第二个支持插入和删除点。 区分这两种样式的原因是
itk中创建点集itk::PointSet
yuxing55555的博客
05-11 434
#include "itkPoint.h" #include "itkPointSet.h" int main(int, char* []) { using PixelType = float; constexpr unsigned int Dimension = 3; using PointSetType = itk::PointSet&lt; PixelType, Dimens...
ITK学习笔记：1.Data presentation-1.3Mesh网格（六）：遍历单元
今天写代码了吗
07-01 203
单元格以指针的方式存储在itk::Mesh中，指针指向通用单元格itk::CellInterface。 这意味着只能调用在这个基单元格类上定义的虚方法。 为了使用特定于每个单元格类型的方法，必须向下转换为指向单元格的实际类型指针。这可以通过利用允许识别单元格的实际类型的GetType()方法安全地完成。 下例，假设建立了一个由一个四面体及其所有边界面定义的网格。也就是四个三角形，六个边和四个顶点。 可以使用CellsContainer迭代器访问单元格。迭代器Value()对应于CellType基类的原始
ITK学习笔记（一）
qq_34760715的博客
06-19 4245
ITK 中表示数据的类有：itk::Image，itk::Mesh 和 itk::PointSet。 itk::Image 是遵循范型编程思想的类，其类型是由算法行为的类演化而来。ITK 支持任何像素类型和空间维的图像。 一 创建图像 #include “itkImage.h” //包含头文件 //设置好模板参数，并给模板起别名，（两种写法都可） //typedef itk::Image...
蓝牙mesh-拓扑结构
qq_37932504的博客
12-27 4257
多个节点可以形成mesh网络。mesh网络的示意图如下所示。 上图显示了三个中继节点:Q、R、S。支持Friend feature的三个节点:N、O、P，而N没有任何friendships;因此，只有O和P是Friend节点。有5个Low Power nodes:I、J、K、L和M。节点I、J和K的friend是P；L和M的friend是O。NodesT仅...
什么是Mesh网络结构？
热门推荐
zpznba的博客
10-27 1万+
本文是对什么是Mesh技术？网状网络（Mesh）和传统网络结构相比有什么优势？这篇文章的简化和总结，原文说的比较啰嗦，想看的点链接自己看。 1.先来看传统局域网的问题： 局域网的基本结构包括星型结构、环型结构、总线结构、分布式结构、树型结构等。最常见最广泛使用的星型结构是以中央节点作为核心，其他节点都连接至中央节点上，这种结构的成本较高、可靠性较低，但是其延迟小、结构简单便于管理。总线型结...
ITK学习笔记：1.Data presentation-1.2PointSet（一）点集：创建、获取与访问
今天写代码了吗
06-24 627
1.创建 PointSet是一个基本类，旨在以n维空间中一组点的形式表示几何。 它是itk::Mesh的基类，提供了操作点集所需的方法。点可以有与之相关的值。这些值的类型是由itk::PointSet类的模板参数定义的。 ITK中有两种基本的切入点集交互样式。这些样式分别称为静态样式（static ）和动态样式（ dynamic）。点集的数量是固定并且不希望改变时用第一种，要增改删则用第二种。 #include "itkPointSet.h" //我们必须决定什么类型的值与这些点相关联。 //为了使术语与
itk软件指南:介绍与开发准则-翻译版下载
最新发布
07-29
itk软件指南是一份全面介绍ITK（Insight ...总之，ITK软件指南的翻译版下载为使用ITK软件的开发人员提供了一个全面介绍和学习ITK软件的工具，帮助他们更好地理解和应用ITK软件，并在科学研究和应用中发挥更大的作用。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值