我从未使用过 CGAL，几乎没有 C/C++ 经验。但以下 谷歌但是我设法编译了示例“Alpha_shapes_3” (\CGAL-4.1-beta1\examples\Alpha_shapes_3) 在 Windows 7 64 位机器上使用 Visual Studio 2010。

现在，如果我们检查程序“ex_alpha_shapes_3”的源代码，我们 请注意，名为“bunny_1000”的数据文件是红色的，其中 3d 点 集群驻留。 现在我的问题是如何更改源代码以便在 alpha 之后 计算给定点的形状，alpha 形状的表面网格是 保存/写入外部文件。它可以只是多边形列表和 它们各自的 3D 顶点。我想这些多边形将定义 alpha 形状的表面网格。如果我能做到，我可以看到输出 我熟悉的外部工具中的 alpha 形状生成程序。

我知道这很简单，但我无法用我的 对 CGAL 的了解有限。

我知道你们有代码，但我再次粘贴它以完成。

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Delaunay_triangulation_3.h>
#include <CGAL/Alpha_shape_3.h>

#include <fstream>
#include <list>
#include <cassert>

typedef CGAL::Exact_predicates_inexact_constructions_kernel Gt;

typedef CGAL::Alpha_shape_vertex_base_3<Gt>          Vb;
typedef CGAL::Alpha_shape_cell_base_3<Gt>            Fb;
typedef CGAL::Triangulation_data_structure_3<Vb,Fb>  Tds;
typedef CGAL::Delaunay_triangulation_3<Gt,Tds>       Triangulation_3;
typedef CGAL::Alpha_shape_3<Triangulation_3>         Alpha_shape_3;

typedef Gt::Point_3                                  Point;
typedef Alpha_shape_3::Alpha_iterator               Alpha_iterator;

int main()
{
  std::list<Point> lp;

  //read input
  std::ifstream is("./data/bunny_1000");
  int n;
  is >> n;
  std::cout << "Reading " << n << " points " << std::endl;
  Point p;
  for( ; n>0 ; n--)    {
    is >> p;
    lp.push_back(p);
  }

  // compute alpha shape
  Alpha_shape_3 as(lp.begin(),lp.end());
  std::cout << "Alpha shape computed in REGULARIZED mode by default"
            << std::endl;

  // find optimal alpha value
  Alpha_iterator opt = as.find_optimal_alpha(1);
  std::cout << "Optimal alpha value to get one connected component is "
            <<  *opt    << std::endl;
  as.set_alpha(*opt);
  assert(as.number_of_solid_components() == 1);
  return 0;
} 

在网上搜索了很多之后我发现我们可能需要使用类似

std::list<Facet> facets;
alpha_shape.get_alpha_shape_facets
(
  std::back_inserter(facets),Alpha_shape::REGULAR
);

但我仍然完全不知道如何在上面的代码中使用它!

最佳答案

记录在案here , 刻面是一对 (Cell_handle c, int i) 定义为 c 中与索引 i 的顶点相对的刻面。 在 this page , 你有一个单元格的顶点索引的描述。

在下面的代码示例中，我添加了一个小输出，用于在 cout 上打印一个 OFF 文件。通过复制顶点。做一些干净的事情，你可以使用 std::map<Alpha_shape_3::Vertex_handle,int>为每个顶点关联一个唯一索引或向顶点添加信息，如 those examples .

/// collect all regular facets
std::vector<Alpha_shape_3::Facet> facets;
as.get_alpha_shape_facets(std::back_inserter(facets), Alpha_shape_3::REGULAR);

std::stringstream pts;
std::stringstream ind;

std::size_t nbf=facets.size();
for (std::size_t i=0;i<nbf;++i)
{ 
  //To have a consistent orientation of the facet, always consider an exterior cell
  if ( as.classify( facets[i].first )!=Alpha_shape_3::EXTERIOR )
    facets[i]=as.mirror_facet( facets[i] );
  CGAL_assertion(  as.classify( facets[i].first )==Alpha_shape_3::EXTERIOR  );

  int indices[3]={
    (facets[i].second+1)%4,
    (facets[i].second+2)%4,
    (facets[i].second+3)%4,
  };

  /// according to the encoding of vertex indices, this is needed to get
  /// a consistent orienation
  if ( facets[i].second%2==0 ) std::swap(indices[0], indices[1]);


  pts << 
  facets[i].first->vertex(indices[0])->point() << "\n" <<
  facets[i].first->vertex(indices[1])->point() << "\n" <<
  facets[i].first->vertex(indices[2])->point() << "\n"; 
  ind << "3 " << 3*i << " " << 3*i+1 << " " << 3*i+2 << "\n";
}

std::cout << "OFF "<< 3*nbf << " " << nbf << " 0\n";
std::cout << pts.str();
std::cout << ind.str();

关于c++ - 保存 CGAL alpha 形状表面网格，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/15905833/