OpenMesh入门2（译自OpenMesh6.3 Documents）

使用迭代器和环绕器

本例主要说明以下问题：

1.如何使用迭代器

2.如何使用环绕器

这个例子是第一版mesh光滑器。这里，我们将介绍迭代器和环绕器，这两个概念能提两个线性枚举功能：一是线性枚举mesh上的所有节点；二是枚举环绕某个节点的所有节点（例如，环绕与某个节点相邻的所有一环节点）。更多的细节描述，可参考Mesh IteratorsandCirculators。

首先，我们要定义我们自己的mesh类型。这次，我们使用三角面替代上节提到的多边形面。

#include <OpenMesh/Core/Mesh/TriMesh_ArrayKernelT.hh>

typedef OpenMesh::TriMesh_ArrayKernelT<>MyMesh;

我们这次从文件中读取一个mesh，使用OpenMesh::IO::read_mesh函数。
if ( ! OpenMesh::IO::read_mesh(mesh,argv[2]) )

构造用于光滑处理的迭代器需要两个步骤：

1.对每个节点：计算其一环临近节点的重心；

2.对每个节点：将其移动到刚刚计算出的重心位置。

通过节点迭代器，可以轻松实现。mesh通过vertices_begin()和vertices_end()函数提供迭代器的开始和结束位置。

MyMesh::VertexIter v_it, v_end(mesh.vertices_end());

for (v_it=mesh.vertices_begin(); v_it!=v_end; ++v_it)

为了计算重心，我们要通过VertexVertexIter提供的迭代器环绕当前节点的一环临近节点。

MyMesh::VertexVertexItervv_it;

for (vv_it=mesh.vv_iter( *v_it ); vv_it.is_valid(); ++vv_it)

现在，我们可以计算存储在数组cogs中的所有节点的重心，代码如下：

std::vector<MyMesh::Point> cogs;

for (v_it=mesh.vertices_begin(); v_it!=v_end; ++v_it)

{

cog[0] = cog[1] = cog[2] = valence = 0.0;

for (vv_it=mesh.vv_iter( *v_it ); vv_it.is_valid(); ++vv_it)

{

cog += mesh.point( *vv_it );

++valence;

}

cogs.push_back(cog / valence);

}

重心计算完成后，剩下的就是将每个节点移动到其对应的重心位置，完整的代码如下：

#include <iostream>

#include <vector>

// -------------------- OpenMesh

#include <OpenMesh/Core/IO/MeshIO.hh>

#include <OpenMesh/Core/Mesh/TriMesh_ArrayKernelT.hh>

typedef OpenMesh::TriMesh_ArrayKernelT<>MyMesh;

int main(int argc, char **argv)

{

MyMesh mesh;

// check command line options

if (argc != 4)

{

std::cerr << "Usage: " <<argv[0] << " #iterations infile outfile\n";

return 1;

}

// read mesh from stdin

if ( ! OpenMesh::IO::read_mesh(mesh,argv[2]) )

{

std::cerr << "Error: Cannot read mesh from " << argv[2] << std::endl;

return 1;

}

// this vector stores the computed centers of gravity

std::vector<MyMesh::Point> cogs;

std::vector<MyMesh::Point>::iterator cog_it;

cogs.reserve(mesh.n_vertices());

// smoothing mesh argv[1] times

MyMesh::VertexIter v_it, v_end(mesh.vertices_end());

MyMesh::VertexVertexItervv_it;

MyMesh::Pointcog;

MyMesh::Scalarvalence;

unsignedint i,N(atoi(argv[1]));

for (i=0; i < N; ++i)

{

cogs.clear();

for (v_it=mesh.vertices_begin(); v_it!=v_end; ++v_it)

{

cog[0] = cog[1] = cog[2] = valence = 0.0;

for (vv_it=mesh.vv_iter( *v_it ); vv_it.is_valid(); ++vv_it)

{

cog += mesh.point( *vv_it );

++valence;

}

cogs.push_back(cog / valence);

}

for (v_it=mesh.vertices_begin(), cog_it=cogs.begin();

v_it!=v_end; ++v_it, ++cog_it)

if ( !mesh.is_boundary( *v_it ) )

mesh.set_point( *v_it, *cog_it );

}

// write mesh to stdout

if ( ! OpenMesh::IO::write_mesh(mesh,argv[3]) )

{

std::cerr << "Error: cannot write mesh to " << argv[3] << std::endl;

return 1;

}

return 0;

}