D:\install\CGAL-5.3\examples\Linear_cell_complex\linear_cell_complex_3_with_colored_vertices.cpp

    image.png
    image.png
    image.png
    image.png

    1. #include <CGAL/Linear_cell_complex_for_combinatorial_map.h>
    2. #include <iostream>
    3. #include <algorithm>
    5. #ifdef _DEBUG
    6. #include <CGAL/draw_linear_cell_complex.h>
    7. #endif // _DEBUG
    10. struct Average_functor
    11. {
    12.     template<class CellAttribute>
    13.     void operator()(CellAttribute& ca1, CellAttribute& ca2)
    14.     {
    15.         ca1.info() = (ca1.info() + ca2.info()) / 2;
    16.     }
    17. };
    19. struct Myitem
    20. {
    21.     template<class Refs>
    22.     struct Dart_wrapper
    23.     {
    24.         typedef CGAL::Cell_attribute_with_point< Refs, int, CGAL::Tag_true,
    25.             Average_functor >
    26.             Vertex_attribute;
    28.         typedef std::tuple<Vertex_attribute> Attributes;
    29.     };
    30. };
    32. typedef CGAL::Linear_cell_complex_traits
    33. <3, CGAL::Exact_predicates_inexact_constructions_kernel> Traits;
    34. typedef CGAL::Linear_cell_complex_for_combinatorial_map<3, 3, Traits, Myitem> LCC_3;
    35. typedef LCC_3::Dart_handle                               Dart_handle;
    36. typedef LCC_3::Point                                     Point;
    37. typedef LCC_3::FT                                        FT;
    39. //创建正方体
    40. //长方体中心点：basepoint
    41. //lg：长方体长度的一半
    42. Dart_handle make_iso_cuboid(LCC_3& lcc, const Point& basepoint, FT lg)
    43. {
    44.     return lcc.make_hexahedron(basepoint,
    45.         Traits::Construct_translated_point()
    46.         (basepoint, Traits::Vector(lg, 0, 0)),
    47.         Traits::Construct_translated_point()
    48.         (basepoint, Traits::Vector(lg, lg, 0)),
    49.         Traits::Construct_translated_point()
    50.         (basepoint, Traits::Vector(0, lg, 0)),
    51.         Traits::Construct_translated_point()
    52.         (basepoint, Traits::Vector(0, lg, lg)),
    53.         Traits::Construct_translated_point()
    54.         (basepoint, Traits::Vector(0, 0, lg)),
    55.         Traits::Construct_translated_point()
    56.         (basepoint, Traits::Vector(lg, 0, lg)),
    57.         Traits::Construct_translated_point()
    58.         (basepoint, Traits::Vector(lg, lg, lg)));
    59. }
    61. int main()
    62. {
    63.     LCC_3 lcc;
    65.     std::cout << "==========创建两个正方体==========" << std::endl;
    67.     // Create two iso_cuboids.
    68.     Dart_handle d1 = make_iso_cuboid(lcc, Point(-2, 0, 0), 1);
    69.     Dart_handle d2 = make_iso_cuboid(lcc, Point(0, 0, 0), 1);
    71. #ifdef _DEBUG
    72.     //绘制
    73.     CGAL::draw(lcc);
    74. #endif // _DEBUG
    76.     std::cout << "==========为两个正方体设置颜色==========" << std::endl;
    78.     //设置第一个正方体所有顶点颜色为1
    79.     // Set the "color" of all vertices of the first cube to 1.
    80.     for (LCC_3::One_dart_per_incident_cell_range<0, 3>::iterator
    81.         it = lcc.one_dart_per_incident_cell<0, 3>(d1).begin(),
    82.         itend = lcc.one_dart_per_incident_cell<0, 3>(d1).end(); it != itend; ++it)
    83.     {
    84.         lcc.info<0>(it) = 1;
    85.     }
    87.     //设置第二个正方体所有顶点颜色为2
    88.     // Set the "color" of all vertices of the second cube to 19.
    89.     for (LCC_3::One_dart_per_incident_cell_range<0, 3>::iterator it =
    90.         lcc.one_dart_per_incident_cell<0, 3>(d2).begin(),
    91.         itend = lcc.one_dart_per_incident_cell<0, 3>(d2).end(); it != itend; ++it)
    92.     {
    93.         lcc.info<0>(it) = 19;
    94.     }
    96. #ifdef _DEBUG
    97.     //绘制
    98.     CGAL::draw(lcc);
    99. #endif // _DEBUG
    101.     std::cout << "==========缝合==========" << std::endl;
    103.     // 3-Sew the two cubes along one facet.
    104.     lcc.sew<3>(lcc.beta(d1, 1, 1, 2), lcc.beta(d2, 2));
    106. #ifdef _DEBUG
    107.     //绘制
    108.     CGAL::draw(lcc);
    109. #endif // _DEBUG
    111.     std::cout << "==========对两个Cube的公共面做Barycentric三角剖分==========" << std::endl;
    113.     //对两个Cube的公共面做Barycentric三角剖分
    114.     // Barycentric triangulation of the facet between the two cubes.
    115.     Dart_handle d3 = lcc.insert_barycenter_in_cell<2>(lcc.beta(d2, 2));
    117.     //设置新顶点的颜色为5。
    118.     // Set the color of the new vertex to 5.
    119.     lcc.info<0>(d3) = 5;
    121.     // Display all the vertices of the map.
    122.     for (LCC_3::Vertex_attribute_range::iterator
    123.         it = lcc.vertex_attributes().begin(),
    124.         itend = lcc.vertex_attributes().end();
    125.         it != itend; ++it)
    126.     {
    127.         std::cout << "point: " << lcc.point_of_vertex_attribute(it) << ", " << "color: "
    128.             << lcc.info_of_attribute<0>(it) << std::endl;
    129.     }
    131. #ifdef _DEBUG
    132.     //绘制
    133.     CGAL::draw(lcc);
    134. #endif // _DEBUG
    136.     return EXIT_SUCCESS;
    137. }