OpenMesh入门2（译自OpenMesh6.3 Documents）
使用迭代器和环绕器
本例主要说明以下问题：
1.如何使用迭代器
2.如何使用环绕器
这个例子是第一版mesh光滑器。这里，我们将介绍迭代器和环绕器，这两个概念能提两个线性枚举功能：一是线性枚举mesh上的所有节点；二是枚举环绕某个节点的所有节点（例如，环绕与某个节点相邻的所有一环节点）。更多的细节描述，可参考Mesh IteratorsandCirculators。
首先，我们要定义我们自己的mesh类型。这次，我们使用三角面替代上节提到的多边形面。
#include
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
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
#include
// -------------------- OpenMesh
#include
#include
typedef OpenMesh::TriMesh_ArrayKernelT<>MyMesh;
int main(int argc, char **argv)
{
MyMesh mesh;
// check command line options
if (argc != 4)
{
std::cerr << "Usage: " <
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
std::vector
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;
}