cgal读取obj格式的三维线段

CGAL（Computational Geometry Algorithms Library）是一个计算几何算法库，支持多种类型的几何计算，包括三维几何、二维几何、三角剖分等。要读取OBJ格式的三维线段，可以使用CGAL中的Polyhedron_3类来实现。

以下是一个示例程序，演示如何使用CGAL读取OBJ格式的三维线段，并输出每个线段的起点和终点坐标：

#include <iostream>
#include <fstream>
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Polyhedron_3.h>
#include <CGAL/Polyhedron_traits_with_normals_3.h>
#include <CGAL/IO/Polyhedron_iostream.h>

typedef CGAL::Exact_predicates_inexact_constructions_kernel Kernel;
typedef CGAL::Polyhedron_3<Kernel> Polyhedron;

int main()
{
    // 读取OBJ文件
    std::ifstream input("lines.obj");
    if (!input)
    {
        std::cerr << "Failed to open file!" << std::endl;
        return -1;
    }

    // 创建Polyhedron对象
    Polyhedron poly;

    // 读取OBJ文件内容到Polyhedron对象
    input >> poly;

    // 遍历Polyhedron对象的所有线段
    for (Polyhedron::Edge_iterator eit = poly.edges_begin(); eit != poly.edges_end(); ++eit)
    {
        // 输出线段的起点和终点坐标
        std::cout << "Start point: (" << eit->vertex()->point().x() << ", "
                  << eit->vertex()->point().y() << ", " << eit->vertex()->point().z() << ")" << std::endl;
        std::cout << "End point: (" << eit->opposite()->vertex()->point().x() << ", "
                  << eit->opposite()->vertex()->point().y() << ", " << eit->opposite()->vertex()->point().z() << ")" << std::endl;
    }

    return 0;
}

在上述示例程序中，我们首先打开一个OBJ格式的文件，然后创建Polyhedron对象，并将OBJ文件中的内容读取到该对象中。然后，我们遍历Polyhedron对象的所有边，输出每个边的起点和终点坐标。

请注意，上述代码需要引入CGAL的头文件和库文件，并且需要使用C++11或更高版本的编译器进行编译。