Using Lipschitz subdomain sizing with a mesh generated from a 3D labeled image

[oarcelus]

I want to use Lipschitz subdomain sizing on some particular labels for my 3D image to make the surface mesh finer with respect to the 'bulk' mesh. But I am not managing to compile my example correctly. I am finding errors related to AABB_Tree.

I am not an expert in CGAL and I kind of mixed both examples from documentation together.

  #define CGAL_MESH_3_VERBOSE 1
  
  #include "mesh_from_labeled_3d_array_lipschitz.hpp"
  
  #include <cassert>
  #include <iostream>
  #include <string.h>
  
  #include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
  #include <CGAL/Image_3.h>
  
  #include <CGAL/Mesh_triangulation_3.h>
  #include <CGAL/Mesh_complex_3_in_triangulation_3.h>
  #include <CGAL/Mesh_criteria_3.h>
  
  #include <CGAL/Implicit_mesh_domain_3.h>
  #include <CGAL/Mesh_domain_with_polyline_features_3.h>
  #include <CGAL/make_mesh_3.h>
  
  #include <CGAL/Mesh_3/experimental/Lipschitz_sizing_polyhedron.h>
  
  
  typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
  typedef CGAL::Labeled_mesh_domain_3<K> Mesh_domain;
  
  // Triangulation
  typedef CGAL::Mesh_triangulation_3<Mesh_domain>::type Tr;
  typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
  
  // Mesh Criteria
  typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
  typedef Mesh_criteria::Facet_criteria Facet_criteria;
  typedef Mesh_criteria::Cell_criteria Cell_criteria;
  
  typedef CGAL::Mesh_3::Lipschitz_sizing<K, Mesh_domain, Mesh_domain::AABB_tree> Lip_sizing;
  typedef CGAL::Mesh_constant_domain_field_3<Mesh_domain::R,
                                             Mesh_domain::Index> Sizing_field_cell;
  
  void
  generate_from_inr_lipschitz(
      const std::string & inr_filename,
      const std::string & outfile,
      const bool lloyd,
      const bool odt,
      const bool perturb,
      const bool exude,
      const double max_edge_size_at_feature_edges,
      const double min_facet_angle,
      const double max_radius_surface_delaunay_ball,
      const double max_facet_distance,
      const double max_circumradius_edge_ratio,
      const double max_cell_circumradius,
      const double exude_time_limit,
      const double exude_sliver_bound,
      const bool verbose,
      const int seed
      )
  {
    CGAL::get_default_random() = CGAL::Random(seed);
  
    CGAL::Image_3 image;
    const bool success = image.read(inr_filename.c_str());
    if (!success) {
      throw "Could not read image file";
    }
    Mesh_domain cgal_domain = Mesh_domain::create_labeled_image_mesh_domain(image);
  
    Lip_sizing lip_sizing(cgal_domain);
  
    lip_sizing.add_parameters_for_subdomain(2, 0.3, 0.02, 0.5);
  
    Mesh_criteria criteria(
        CGAL::parameters::edge_size=max_edge_size_at_feature_edges,
        CGAL::parameters::facet_angle=min_facet_angle,
        CGAL::parameters::facet_size=max_radius_surface_delaunay_ball,
        CGAL::parameters::facet_distance=max_facet_distance,
        CGAL::parameters::cell_radius_edge_ratio=max_circumradius_edge_ratio,
        CGAL::parameters::cell_size=lip_sizing
        );
  
    // Mesh generation
    if (!verbose) {
      // suppress output
      std::cerr.setstate(std::ios_base::failbit);
    }
    C3t3 c3t3 = CGAL::make_mesh_3<C3t3>(
        cgal_domain,
        criteria,
        lloyd ? CGAL::parameters::lloyd(CGAL::parameters::time_limit(0)) : CGAL::parameters::no_lloyd(),
        odt ? CGAL::parameters::odt(CGAL::parameters::time_limit(0)) : CGAL::parameters::no_odt(),
        perturb ? CGAL::parameters::perturb() : CGAL::parameters::no_perturb(),
        exude ?
          CGAL::parameters::exude(
            CGAL::parameters::time_limit = exude_time_limit,
            CGAL::parameters::sliver_bound = exude_sliver_bound
          ) :
          CGAL::parameters::no_exude()
        );
    if (!verbose) {
      std::cerr.clear();
    }
  
    // Output
    std::ofstream medit_file(outfile);
    c3t3.output_to_medit(medit_file);
    medit_file.close();
    return;
  }
  
  int main(int argc, char* argv[])
  {
      // std::string mode = "image";
      // const char *modes = &mode[0];
      const std::string inr_filename = "NMC75_LMO25_p0.51_Thick44um_Loading1.00.inr";
      const std::string outfile = "NMC75_LMO25_p0.51_Thick44um_Loading1.00.mesh";
      const double default_max_cell_circumradius = 0.0;
      const std::vector<double> max_cell_circumradiuss{ 0.0, 0.0 };
      const std::vector<int> cell_labels{ 1, 2 };
      const bool lloyd { true };
      const bool odt { true };
      const bool perturb { true };
      const bool exude { false };
      const double relative_error_bound_weights = 1e-6;
      const double max_edge_size_at_feature_edges = 0.0;
      const double min_facet_angle = 25.0;
      const double max_radius_surface_delaunay_ball = 1.0;
      const double max_facet_distance = 0.0;
      const double max_circumradius_edge_ratio = 0.0;
      const double max_cell_circumradius = 2.0;
      const double exude_time_limit = 0.0;
      const double exude_sliver_bound = 0.0;
      const bool verbose { true };
      const int seed = 1;
  
      generate_from_inr_lipschitz(inr_filename,
                          outfile,
                          lloyd,
                          odt,
                          perturb,
                          exude,
                          max_edge_size_at_feature_edges,
                          min_facet_angle,
                          max_radius_surface_delaunay_ball,
                          max_facet_distance,
                          max_circumradius_edge_ratio,
                          max_cell_circumradius,
                          exude_time_limit,
                          exude_sliver_bound,
                          verbose,
                          seed);
  
      return 0;
  }

The error message reads as follows

  C:\Users\oarcelus\ElectrodeBuilder\meshing_cgal\mesh_from_labeled_3d_array_lipschitz.cpp(35,69): error C2039:
  'AABB_tree': is not a member of 'CGAL::Labeled_mesh_domain_3<K,int,std::pair<int,int>>' [C:\Users\oarcelus\E 
  lectrodeBuilder\meshing_cgal\build\mesh_from_labeled_3d_array_lipschitz.vcxproj]
  
  C:\Users\oarcelus\ElectrodeBuilder\meshing_cgal\mesh_from_labeled_3d_array_lipschitz.cpp(35,78): error C2065: 
  'AABB_tree': undeclared identifier [C:\Users\oarcelus\ElectrodeBuilder\meshing_cgal\build\mesh_from_labeled_ 
  3d_array_lipschitz.vcxproj]
  C:\Users\oarcelus\ElectrodeBuilder\meshing_cgal\mesh_from_labeled_3d_array_lipschitz.cpp(35,23): error C2923: 
  'CGAL::Mesh_3::Lipschitz_sizing': 'AABB_tree' is not a valid template type argument for parameter 'AABBTreeT 
  emplate' [C:\Users\oarcelus\ElectrodeBuilder\meshing_cgal\build\mesh_from_labeled_3d_array_lipschitz.vcxproj] 
  
  C:\Users\oarcelus\ElectrodeBuilder\meshing_cgal\mesh_from_labeled_3d_array_lipschitz.cpp(68,36): error C2514:
  'CGAL::Mesh_3::Lipschitz_sizing': class template cannot be constructed [C:\Users\oarcelus\ElectrodeBuilder\m 
  eshing_cgal\build\mesh_from_labeled_3d_array_lipschitz.vcxproj]
  
  C:\Users\oarcelus\ElectrodeBuilder\meshing_cgal\mesh_from_labeled_3d_array_lipschitz.cpp(70,60): error C2662: 
  'void CGAL::Mesh_3::Lipschitz_sizing<Kernel,MeshDomain,AABBTreeTemplate>::add_parameters_for_subdomain(const 
  MeshDomain::Subdomain_index &,const Kernel::FT &,const Kernel::FT &,const Kernel::FT &)': cannot convert 'th 
  is' pointer from 'Lip_sizing' to 'CGAL::Mesh_3::Lipschitz_sizing<Kernel,MeshDomain,AABBTreeTemplate> &' [C:\U 
  sers\oarcelus\ElectrodeBuilder\meshing_cgal\build\mesh_from_labeled_3d_array_lipschitz.vcxproj]
  
  C:\Users\oarcelus\ElectrodeBuilder\meshing_cgal\mesh_from_labeled_3d_array_lipschitz.cpp(78,25): error C2678:
  binary '=': no operator found which takes a left-hand operand of type 'const CGAL::parameters::Boost_paramet 
  er_compatibility_wrapper<CGAL::internal_np::cell_size_param_t,false,false>' (or there is no acceptable conver
  sion) [C:\Users\oarcelus\ElectrodeBuilder\meshing_cgal\build\mesh_from_labeled_3d_array_lipschitz.vcxproj] 

Any help or tip would be much appreciated.