使用CGAL isotropic_remeshing进行约束求精

sudo apt-get install libcgal-dev

cmake_minimum_required(VERSION 3.1)
project(IsotropicRemeshing)

find_package(CGAL REQUIRED)

add_executable(isotropic_remeshing main.cpp)
target_link_libraries(isotropic_remeshing CGAL::CGAL)

#include <CGAL/Simple_cartesian.h>
#include <CGAL/Surface_mesh.h>
#include <CGAL/Polygon_mesh_processing/remesh.h>
#include <CGAL/Polygon_mesh_processing/IO/polygon_mesh_io.h>

namespace PMP = CGAL::Polygon_mesh_processing;

typedef CGAL::Simple_cartesian<double> Kernel;
typedef Kernel::Point_3 Point;
typedef CGAL::Surface_mesh<Point> Mesh;
typedef boost::graph_traits<Mesh>::vertex_descriptor vertex_descriptor;
typedef boost::graph_traits<Mesh>::halfedge_descriptor halfedge_descriptor;
typedef boost::graph_traits<Mesh>::edge_descriptor edge_descriptor;

int main(int argc, char* argv[]) {
    if (argc != 2) {
        std::cerr << "Usage: " << argv[0] << " input.off" << std::endl;
        return EXIT_FAILURE;
    }

    Mesh mesh;
    if (!PMP::IO::read_polygon_mesh(argv[1], mesh) || !CGAL::is_triangle_mesh(mesh)) {
        std::cerr << "Invalid input mesh." << std::endl;
        return EXIT_FAILURE;
    }

    // 标记要保留的边（例如边界边）
    Mesh::Property_map<edge_descriptor, bool> is_constrained_map =
        mesh.add_property_map<edge_descriptor, bool>("e:is_constrained", false).first;

    for (edge_descriptor e : edges(mesh)) {
        if (CGAL::is_border(e, mesh)) {
            is_constrained_map[e] = true;
        }
    }

    // 进行各向同性网格重构
    double target_edge_length = 0.1;
    unsigned int nb_iter = 3;
    PMP::isotropic_remeshing(
        faces(mesh),
        target_edge_length,
        mesh,
        PMP::parameters::number_of_iterations(nb_iter)
                        .edge_is_constrained_map(is_constrained_map)
    );

    // 输出重构后的网格
    PMP::IO::write_polygon_mesh("output.off", mesh, CGAL::parameters::stream_precision(17));

    return EXIT_SUCCESS;
}

mkdir build
cd build
cmake ..
make

./isotropic_remeshing ../input.off